计算化学公社

标题: ADMP的bug还是我的输入有误 [打印本页]

作者
Author: diaolanxinyu 时间: 2015-10-5 12:29
标题: ADMP的bug还是我的输入有误
ADMP计算时出错，麻烦大家看一下，谢谢！
输入文件为：
%mem=1GB
#p b3lyp/6-31g* admp=(StepSize=2000,maxpoint=1000,ReadMWVelocity) nosymm

Title Card Required

-1 2
C                -0.00153004 -0.03233370 -0.02607732
H                0.54723535 0.75821432 0.62498698
H                -0.65691642 -0.58663561 0.64645667
H                0.77087995 -0.55666386 -0.44363167
H                -0.64298093 0.77007701 -0.51731380

0
5.347549407213D+12 -2.018037913459D+13 -2.452356984382D+12
-3.583023315896D+13 -2.984296587042D+13 2.779314553014D+13
-4.928966973867D+12 1.433279503340D+13 8.431574743216D+13
8.455608905705D+13 -2.792759767792D+13 -4.091384745823D+13
-6.224928878988D+13 1.130727352774D+14 -6.273287559347D+13

出错信息为：
UHF open shell SCF:
Using DIIS extrapolation, IDIIS=  1040.
Two-electron integral symmetry not used.
Keep R1 and R2 ints in memory in canonical form, NReq=956475.
IVT=    20469 IEndB=    20469 NGot= 134217728 MDV= 134120808
LenX= 134120808 LenY= 134119838
Using Cholesky orthonormal basis.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on          energy=1.00D-06.
No special actions if energy rises.
Error in input density:  TNEA=************ NAE=    6 TNEB=************ NBE=    5.
Input densites are not normalized.
Error termination via Lnk1e in /opt/gaussian/g09d01/l502.exe at Mon Oct  5 11:21:39 2015.
Job cpu time:    0 days  0 hours  0 minutes 15.9 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    2 Scr=    1

当把关键字 ADMP换为 BOMD就可以计算了，不知是gaussian的bug还是我的输入有误，麻烦大家看一下，谢谢！

作者
Author: sobereva 时间: 2015-10-5 15:18
大抵是初始速度设得比较大，结构改变得大，导致按照ADMP的算法跑，到第二步结构的时候已经没法正常读取之前的波函数了。
加上FullSCF可解决，不过这样的话比起BOMD倒也没什么差异了。

作者
Author: diaolanxinyu 时间: 2015-10-6 22:35
谢谢sobereva老师，问题按照您的方法已经解决。
虽然每一步跟BOMD一样都会有迭代，但还是会比BOMD容易收敛的多，我想应该与局部二次曲面上校正有关。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)