| 查看: 1890 | 回复: 0 | ||
让你变成回忆木虫 (小有名气)
|
[求助]
高斯09 ccsd计算单点能出现能量振荡，
|
|
小弟在用ccsd方法计算一个体系（体系一般不算很大）的单点能，选择的基组较大，选的是aug-cc-pvdz，在计算完SCF后，出现了不收敛的问题，我将收敛标准改为了SCF(CONVER=6)这样能计算出需要的能量，这个值准确么？我看到有人说一般scf计算conver=6就可以了； 另外，在计算ccsd方法时，当iteration Nr 到默认最大的50时，还没收敛，我改成了ccsd（maxcyc=200）计算后 仍然计算到200时还是不收敛， 对于ccsd计算出现不收敛应该怎么解决呢？ 收敛标准改为多少合适？ 看了说明书，但是感觉不是很清楚。请各位大神帮助！ 这是我用grep命令查看log文件的能量值 DE(Corr)= -2.0055721 E(CORR)= -541.43415635 Delta=-5.48D-02 DE(Corr)= -2.0114485 E(CORR)= -541.44003280 Delta=-5.88D-03 DE(Corr)= -2.0165933 E(CORR)= -541.44517754 Delta=-5.14D-03 DE(Corr)= -2.0171844 E(CORR)= -541.44576866 Delta=-5.91D-04 DE(Corr)= -2.0175854 E(CORR)= -541.44616966 Delta=-4.01D-04 DE(Corr)= -2.0176548 E(CORR)= -541.44623906 Delta=-6.94D-05 DE(Corr)= -2.0177208 E(CORR)= -541.44630511 Delta=-6.60D-05 DE(Corr)= -2.0177449 E(CORR)= -541.44632916 Delta=-2.41D-05 DE(Corr)= -2.0177084 E(CORR)= -541.44629272 Delta= 3.64D-05 DE(Corr)= -2.0177151 E(CORR)= -541.44629940 Delta=-6.68D-06 DE(Corr)= -2.0177092 E(CORR)= -541.44629347 Delta= 5.93D-06 DE(Corr)= -2.0177074 E(CORR)= -541.44629169 Delta= 1.79D-06 DE(Corr)= -2.0177068 E(CORR)= -541.44629106 Delta= 6.29D-07 DE(Corr)= -2.0177077 E(CORR)= -541.44629199 Delta=-9.29D-07 DE(Corr)= -2.0177090 E(CORR)= -541.44629323 Delta=-1.24D-06 DE(Corr)= -2.0177087 E(CORR)= -541.44629301 Delta= 2.24D-07 DE(Corr)= -2.0177089 E(CORR)= -541.44629322 Delta=-2.08D-07 DE(Corr)= -2.0177092 E(CORR)= -541.44629346 Delta=-2.42D-07 DE(Corr)= -2.0177099 E(CORR)= -541.44629417 Delta=-7.12D-07 DE(Corr)= -2.0177102 E(CORR)= -541.44629447 Delta=-2.94D-07 DE(Corr)= -2.0177102 E(CORR)= -541.44629443 Delta= 3.11D-08 DE(Corr)= -2.0177102 E(CORR)= -541.44629450 Delta=-6.15D-08 DE(Corr)= -2.0177103 E(CORR)= -541.44629461 Delta=-1.17D-07 DE(Corr)= -2.0177107 E(CORR)= -541.44629495 Delta=-3.37D-07 DE(Corr)= -2.0177107 E(CORR)= -541.44629496 Delta=-7.01D-09 DE(Corr)= -2.0177106 E(CORR)= -541.44629485 Delta= 1.06D-07 DE(Corr)= -2.0177107 E(CORR)= -541.44629496 Delta=-1.13D-07 DE(Corr)= -2.0177107 E(CORR)= -541.44629494 Delta= 2.19D-08 DE(Corr)= -2.0177108 E(CORR)= -541.44629513 Delta=-1.88D-07 DE(Corr)= -2.0177106 E(CORR)= -541.44629492 Delta= 2.05D-07 DE(Corr)= -2.0177107 E(CORR)= -541.44629499 Delta=-6.81D-08 DE(Corr)= -2.0177108 E(CORR)= -541.44629506 Delta=-6.89D-08 DE(Corr)= -2.0177107 E(CORR)= -541.44629502 Delta= 4.01D-08 DE(Corr)= -2.0177108 E(CORR)= -541.44629508 Delta=-6.34D-08 DE(Corr)= -2.0177106 E(CORR)= -541.44629491 Delta= 1.78D-07 DE(Corr)= -2.0177109 E(CORR)= -541.44629516 Delta=-2.50D-07 DE(Corr)= -2.0177108 E(CORR)= -541.44629506 Delta= 9.95D-08 DE(Corr)= -2.0177107 E(CORR)= -541.44629500 Delta= 6.06D-08 DE(Corr)= -2.0177108 E(CORR)= -541.44629507 Delta=-7.08D-08 DE(Corr)= -2.0177107 E(CORR)= -541.44629497 Delta= 9.40D-08 DE(Corr)= -2.0177109 E(CORR)= -541.44629518 Delta=-2.04D-07 DE(Corr)= -2.0177106 E(CORR)= -541.44629488 Delta= 2.95D-07 DE(Corr)= -2.0177107 E(CORR)= -541.44629501 Delta=-1.28D-07 DE(Corr)= -2.0177108 E(CORR)= -541.44629508 Delta=-7.52D-08 DE(Corr)= -2.0177107 E(CORR)= -541.44629496 Delta= 1.21D-07 DE(Corr)= -2.0177108 E(CORR)= -541.44629512 Delta=-1.53D-07 DE(Corr)= -2.0177106 E(CORR)= -541.44629488 Delta= 2.40D-07 DE(Corr)= -2.0177109 E(CORR)= -541.44629521 Delta=-3.37D-07 DE(Corr)= -2.0177108 E(CORR)= -541.44629504 Delta= 1.77D-07 能量出现振荡，应该怎么处理呢。 这是log文件的最后几行： Iteration Nr. 50 ********************** DD1Dir will call FoFJK 3 times, MxPair= 444 NAB= 666 NAA= 0 NBB= 0 NumPrc= 16. DD1Dir: Pass 1 NRAB= 222 NRAA= 0 NRBB= 0 NR= 444. DD1Dir: Pass 2 NRAB= 222 NRAA= 0 NRBB= 0 NR= 444. DD1Dir: Pass 3 NRAB= 222 NRAA= 0 NRBB= 0 NR= 444. Norm of the A-vectors is 6.4985416D-03 conv= 1.00D-05. RLE energy= -2.0177107770 DE(Corr)= -2.0177108 E(CORR)= -541.44629504 Delta= 1.77D-07 NORM(A)= 0.13677643D+01 ************* *MAX. CYCLES* ************* Largest amplitude= 6.28D-02 Error termination via Lnk1e in /home/gaussian/g09/l913.exe at Fri Nov 13 13:21:14 2015. Job cpu time: 38 days 17 hours 38 minutes 18.8 seconds. File lengths (MBytes): RWF= 146091 Int= 0 D2E= 0 Chk= 4 Scr= 1 |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
求助火焰封管的时候管子炸了
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有275人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
七嗪类物质合成求助
已经有0人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
