| 查看: 848 | 回复: 0 | ||
ra2ghgzh木虫 (正式写手)
学痴
|
[求助]
优化聚合物结构时出错Total charge is not ZERO
|
|
优化多聚噻吩的结构时出错，输入文件如下： %mem=600MW %chk=poly2.chk %nprocshared=16 # opt b3lyp/cc-pvtz Title Card Required 0 1 C 1.30606351 0.80252407 0.07216115 C 0.68735101 -0.43640823 -0.03787158 C -0.68979218 -0.43646160 -0.03914487 C -1.30860459 0.80232777 0.07059102 S -0.00153060 2.02649557 0.17958089 H 1.20500977 -1.38091472 -0.12202474 H -1.20800469 -1.38055460 -0.12333028 Tv 4.08921795 0.00000000 0.00000000 报错内容是： -------------------- ATTENTION-------------------- Total charge is not ZERO --- 0.000000001117 -------------------- ATTENTION-------------------- Possible problem detected with PBC diagonalization: FLenRc= 6.50D+00 SDif= 2.67D+07 STest= 1.73D+08 ThrSDi= 1.00D+01. Alpha NIter= 20 Fermi level: -50.2079518603 FoFJK: IHMeth= 1 ICntrl= 0 DoSepK=T KAlg= 1 I1Cent= 0 FoldK=F Possible problem detected with PBC diagonalization: FLenRc= 6.50D+00 SDif= 2.66D+07 STest= 1.73D+08 ThrSDi= 1.00D+01. Alpha NIter= 20 Fermi level: -69.0222808518 <P*S>= 41.999999994528 NE= 42.000000000000 Err=-5.47D-09 Lim= 4.20D-09 Possibly not enough cells included in PBC. Numerical problems in PBC. Error termination via Lnk1e in /share/apps/g09/l502.exe at Tue Feb 14 16:39:46 2017. Job cpu time: 0 days 16 hours 44 minutes 3.8 seconds. File lengths (MBytes): RWF= 480 Int= 0 D2E= 0 Chk= 28 Scr= 1 之前已经用pvdz成功优化出结构了（就是输入文件中的结构），但是改用pvtz优化后就出错了。 这个结构应该很简单吧，但是为什么会出现总电荷不是零的问题呢？ |
回复此楼» 猜你喜欢
分子动力学 Plumed教程求助
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有112人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有63人回复
20