请教：Amber通过Gaussian调用xtb的问题

霜晨月 · 发表于 Post on 2019-12-23 23:43:06

用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的关键词相比，似乎没有什么异常。

请问这是什么原因？有没有老师试过这种调用方法？
多谢

sobereva · 发表于 Post on 2019-12-23 23:52:32

之前论坛里有专门的帖子
http://bbs.keinsci.com/thread-15304-1-1.html

霜晨月 · 发表于 Post on 2019-12-23 23:58:32

谢谢社长，之前看过那个帖子，只是orca和xtb编译都不熟悉，怕麻烦所以想用gaussian来调用

霜晨月 · 发表于 Post on 2019-12-23 23:59:32

这么说，通过gaussian来调用xtb怕是行不通了

sobereva · 发表于 Post on 2019-12-24 06:54:58

ORCA本来就不需要编译。
xtb的编译看
Grimme的xtb程序的编译方法
http://sobereva.com/521

霜晨月 · 发表于 Post on 2019-12-24 09:16:24

多谢社长，我试试看

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Amber] 请教：Amber通过Gaussian调用xtb的问题

浏览过的版块