ADMP的bug还是我的输入有误

diaolanxinyu

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还是我的输入有误，麻烦大家看一下，谢谢！

sobereva

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

diaolanxinyu

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