求助：Gaussian中BOMD动力学能量突变是正常的吗？

gezq

各位老师，我在用Gaussian中BOMD跑某一分子的动力学过程中，发生有一步发生能量突变（如下图所示）。随后，我又对该分子重复了5次测试，后面5次的计算过程中，均在该突变点(具体来说是第51步)由于SCF不收敛而停止了运行，只有第1次是收敛的并继续进行了运算且发生了能量突变（尝试了很多种SCF不收敛的解决办法，目前尚未解决），请问这一动力学过程是否存在什么异常？
【注】目标希望其自发解离，而解离过程发生在第57步，突变点在第51步；输出文件附件已上传

chands

先查结构是否有大变化，比如键断了，分子散架了之类。

gezq

chands 发表于 2025-6-2 16:56
先查结构是否有大变化，比如键断了，分子散架了之类。

没有散架，目标还真是希望其自发解离，而解离过程发生在第57步，突变点在第51步

sobereva

没输出文件，别人靠盲猜没法准确判断

chands

gezq 发表于 2025-6-2 17:01
没有散架，目标还真是希望其自发解离，而解离过程发生在第57步，突变点在第51步

那上传文件吧。

gezq

sobereva 发表于 2025-6-2 17:05
没输出文件，别人靠盲猜没法准确判断

好的，老师，输出文件已上传，1.log为第1次跑的结果，2.log为重复测试的结果。第一次计算因为一些原因而中途停止了计算

gezq

chands 发表于 2025-6-2 17:05
那上传文件吧。

好的，老师，文件已上传

sobereva

我就看了1.log，加粗的是第51次输出的S1激发能，也正是突越的那一帧，可见激发能近乎为0。此时TDDFT完全不适用。你需要用合适的方法做非绝热动力学。

47918  Excited State   1:  2.004-?Sym    0.1502 eV 8255.93 nm  f=0.0002  <S**2>=0.754
48925  Excited State   1:  2.003-?Sym    0.1105 eV 11223.95 nm  f=0.0002  <S**2>=0.753
50013  Excited State   1:  2.001-?Sym    0.0560 eV 22136.96 nm  f=0.0002  <S**2>=0.751
52671  Excited State   1:  1.998-?Sym    0.0236 eV 52484.76 nm  f=0.0004  <S**2>=0.748
53837  Excited State   1:  2.006-?Sym    0.2578 eV 4809.94 nm  f=0.0003  <S**2>=0.756
54793  Excited State   1:  2.004-?Sym    0.3970 eV 3123.32 nm  f=0.0002

gezq

sobereva 发表于 2025-6-2 18:21
我就看了1.log，加粗的是第51次输出的S1激发能，也正是突越的那一帧，可见激发能近乎为0。此时TDDFT完全不 ...

好的，老师。请问在哪个培训班课件的哪个部分里有讲非绝热动力学的部分，我去学一下，平时太忙，没有精力一点一点去听录音，麻烦老师直接帮我定位一下。

sobereva

gezq 发表于 2025-6-2 19:07
好的，老师。请问在哪个培训班课件的哪个部分里有讲非绝热动力学的部分，我去学一下，平时太忙，没有精力 ...

预计第2届高级量子化学培训班会讲非绝热动力学的模拟，会有专门的一节