请问用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

wzkchem5

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