| 查看: 797 | 回复: 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优化后就出错了。 这个结构应该很简单吧，但是为什么会出现总电荷不是零的问题呢？ |
回复此楼» 猜你喜欢
ELISA试验中不容忽视的细节盘点
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有146人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
分子动力学模拟合作需求
已经有2人回复
5