计算化学公社

标题: 请问用win64版Gaussian16进行casscf分子动力学模拟时卡住如何解决 [打印本页]

作者
Author: 柯柯爱喝热可可 时间: 2023-10-19 19:59
标题: 请问用win64版Gaussian16进行casscf分子动力学模拟时卡住如何解决
在模拟D2二聚体电离一个电子之后产生D3+的过程出现了这个问题，一开始搜了一下以为是默认内存不够，但是改变内存也还是每次都卡在第130步，然后去下了linux版的Gaussian计算还是卡在了这一步提示Segmentation fault (core dumped)，报错信息在最后，希望有老师能指点一下这个到底是什么原因，要怎么解决，十分感谢！
输入文件如下，
%chk=D:\WORK\H3\20231017\opt1jia10.chk
#p casscf(3,4)/aug-cc-pvtz bomd=(stepsize=1000,maxpoints=200,randomvelocity)

Title Card Required

1 2
H(Iso=2)          0.00000000 0.36747080 -1.88570370
H(Iso=2)          0.00000000 -0.36747080 -1.88570370
H(Iso=2)          0.00000000 0.00000000 1.51829070

H(Iso=2)          0.00000000 0.00000000 2.25311669

win64版的Gaussian16输出文件最后如下
Leave Link  202 at Wed Oct 18 15:29:50 2023, MaxMem= 104857600 cpu:             0.0 elap:             0.0
(Enter D:\GS\G16W\l301.exe)
Standard basis: Aug-CC-pVTZ (5D, 7F)
Ernie: Thresh=  0.10000D-02 Tol=  0.10000D-05 Strict=F.
92 basis functions, 108 primitive gaussians, 100 cartesian basis functions
   2 alpha electrons       1 beta electrons
   nuclear repulsion energy       0.0000000000 Hartrees.
IExCor= 0 DFT=F Ex=HF Corr=None ExCW=0 ScaHFX=  1.000000
ScaDFX=  1.000000  1.000000  1.000000  1.000000 ScalE2=  1.000000  1.000000
IRadAn=    5 IRanWt=    -1 IRanGd=          0 ICorTp=0 IEmpDi=  4
NAtoms= 4 NActive= 4 NUniq= 4 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F
Integral buffers will be 131072 words long.
Regular integral format.
Two-electron integral symmetry is turned off.
Leave Link  301 at Wed Oct 18 15:29:50 2023, MaxMem= 104857600 cpu:             0.0 elap:             0.0

linux版的log文件最后报错如下
DSYEV returned Info=          2 IAlg= 2 N=    3 NDim=    3 NE2=          2.
Diagonalization in DiagDN failed.
Error termination via Lnk1e in /home/yoyo/gaussian/g16/l118.exe at Thu Oct 19 18:12:22 2023.
Job cpu time:    0 days 15 hours 29 minutes 54.8 seconds.
Elapsed time:    0 days  2 hours 36 minutes 22.7 seconds.
File lengths (MBytes):  RWF=    15 Int=    0 D2E=    0 Chk=    2 Scr=    1

作者
Author: wzkchem5 时间: 2023-10-20 00:22
建议改用Molpro、ORCA等程序，做至少NEVPT2级别的AIMD，活性空间不要只选价层空间，多选一些轨道。这么小体系可能用FCI跑动力学都跑得动。高斯在多参考态计算方面远不如这些软件好，有解决这个问题的时间，早就把Molpro或ORCA学会甚至跑完了。
更好的做法是用FCI参考数据拟合势能面，然后做PIMD，考虑隧穿效应。三、四个轻原子的体系，不考虑核量子效应说不过去

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