24小时热门版块排行榜

返回列表

让你变成回忆

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 4257.4
红花: 1
帖子: 77
在线: 58.8小时
虫号: 3297078
注册: 2014-06-28
性别: GG
专业: 理论和计算化学

[求助] 高斯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

回复此楼

» 猜你喜欢

华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员已经有11人回复
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员已经有10人回复
物理化学论文润色/翻译怎么收费? 已经有75人回复
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后已经有1人回复
200年来人们一直无法破解的难题：气液临界点状态方程已经有7人回复
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！已经有10人回复
求助立方晶体In(OH)3的CIF卡片已经有2人回复
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340 已经有2人回复
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料已经有0人回复
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建已经有0人回复

1楼2015-11-17 20:46:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主让你变成回忆的主题更新

返回列表