标题: 请教：Amber通过Gaussian调用xtb的问题 [打印本页] 作者Author: 霜晨月 时间: 2019-12-23 23:43 标题: 请教：Amber通过Gaussian调用xtb的问题 用Amber18调用Gaussian 09跑QM/MM MD，速度太慢，16核1.7ps/day。
我就想起用社长开发的Gaussian调用xtb的工具，从前使用过，运行没问题。所以我想用Amber调用g09，g09再调用xtb，方法是在“gau_job.tpl”文件中写上 # external='./xtb.sh' （此外什么都没写）。
运行之后秒报错退出，错误信息是：
Generating mol.tmp
Generating mol.xyz via genxyz
Running: xtb mol.xyz -P 4 --chrg -2 --uhf 0 --grad > xtbout
xtb running finished!
Extracting data from xtb outputs via extderi
Recovered energy= -94.2732038000 dipole= 0.000000000000 0.000000000000 0.000000000000
Energy= -94.2732038 NIter= 0.
Dipole moment= 0.000000 0.000000 0.000000 Unabled to recover density for method 0
Error termination via Lnk1e in /home/soft/g09/g09/l602.exe at Mon Dec 23 23:19:04 2019.
Job cpu time: 0 days 0 hours 0 minutes 0.7 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1
sander的.out文件末尾报错：
4. RESULTS
-------------------------------------------------------------------------------- SANDER BOMB in subroutine get_gau_forces (qm2_extern_gau_module)
Error with system call (executing Gaussian)
Will quit now.
检查了sander生成的gau_job.inp，关键词是：
# external='./xtb.sh' NoSymm Force Charge Prop=(Field,Read)
跟直接调用g09的关键词相比，似乎没有什么异常。