| 查看: 849 | 回复: 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优化后就出错了。 这个结构应该很简单吧，但是为什么会出现总电荷不是零的问题呢？ |
回复此楼» 猜你喜欢
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
招收博士一名-宁波大学
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有61人回复
-
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生
已经有2人回复
-
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有1人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
20