计算三重激发态能量时，T1数值为负值

wza1477 · 发表于 Post on 2025-11-19 11:18:56

各位老师好，本人在计算一个分子的三重激发态能量时，最终得到的T1值为负值（-0.39 eV），我检查了输入文件，没有发现输入文件存在问题，请问这属于正常情况吗？谢谢各位老师解答。

输入文件：
%chk=Pt1_4_td_triplets.chk
%nprocshared=25
%mem=30000MB
#p BhandHLYP/gen pseudo=read td=(triplets, direct, nstates=10) scrf=(pcm, solvent=chloroform) test

td

0 1

C,xxxxxxxx,xxxxxxxx,xxxxxxxx
......
O,xxxxxxxx,xxxxxxxx,xxxxxxxx

C H N O S P 0
6-31G*
****
Pt 0
LANL2DZ
****

Pt 0
LANL2DZ

wza1477 · 发表于 Post on 2025-11-20 22:10:19

sobereva 发表于 2025-11-20 21:07
是
都说了，甭写test。下次再看见不再回复

好的，谢谢sob老师解答与提醒，我今后一定注意！

sobereva · 发表于 Post on 2025-11-20 21:07:39

wza1477 发表于 2025-11-20 14:56
sob老师您好，我按照您说的方法先加入stable=opt进行了单点计算，关键词如下：
%chk=Pt1_4_2_sp.chk
%n ...

是
都说了，甭写test。下次再看见不再回复

wza1477 · 发表于 Post on 2025-11-20 14:56:06

sobereva 发表于 2025-11-19 23:54
不要写毫无意义的test关键词。direct也是多余的

参考北京科音中级量子化学培训班（http://www.keinsci.c ...

sob老师您好，我按照您说的方法先加入stable=opt进行了单点计算，关键词如下：
%chk=Pt1_4_2_sp.chk
%nprocshared=25
%mem=30000MB
#p BhandHLYP/gen stable=opt pseudo=read scrf=(pcm, solvent=dichloromethane) test

sp

0 1

该结果文件中显示“The wavefunction is stable under the perturbations considered.”说明波函数应该已经稳定了。那再做激发态计算时加入guess=read读取波函数作为初猜时，输入文件的关键词是不是应该像下面这样呢？
%oldchk=Pt1_4_2_sp.chk
%chk=Pt1_4_3_td_triplets.chk
%nprocshared=25
%mem=30000MB
#p BhandHLYP/gen pseudo=read td=(triplets, nstates=10) guess=read scrf=(pcm, solvent=chloroform) test

td

0 1

wza1477 · 发表于 Post on 2025-11-20 14:32:37

本帖最后由 wza1477 于 2025-11-20 14:44 编辑

sobereva 发表于 2025-11-19 23:54
不要写毫无意义的test关键词。direct也是多余的

参考北京科音中级量子化学培训班（http://www.keinsci.c ...

sob老师您好，我先加入了stable=opt进行了单点计算，关键词行如下：
%chk=Pt1_4_2_sp.chk
%nprocshared=25
%mem=30000MB
#p BhandHLYP/gen stable=opt pseudo=read scrf=(pcm, solvent=dichloromethane) test

sp

0 1

后读取log文件发现“wavefunction is stable under the perturbations considered”，说明波函数应该已经稳定了。那么再做激发态计算时加入guess=read读取其作为初猜时，输入文件的关键词是不是像下面这样写呢？

%oldchk=Pt1_4_2_sp.chk
%chk=Pt1_4_3_td_triplets.chk
%nprocshared=25
%mem=30000MB
#p BhandHLYP/gen pseudo=read td=(triplets, nstates=10) guess=read scrf=(pcm, solvent=chloroform) test

td

0 1

wza1477 · 发表于 Post on 2025-11-20 13:59:48

iamjjLin 发表于 2025-11-20 11:27
建议先仔细检查结构，特别是H原子，多了少了都会导致激发态为负。结构确定没问题后再尝试稳úê�

感谢您的提醒！此前已经检查过结构，应该不存在问题。

iamjjLin · 发表于 Post on 2025-11-20 11:27:39

建议先仔细检查结构，特别是H原子，多了少了都会导致激发态为负。结构确定没问题后再尝试稳定波函数。

wza1477 · 发表于 Post on 2025-11-20 08:24:02

sobereva 发表于 2025-11-19 23:54
不要写毫无意义的test关键词。direct也是多余的

参考北京科音中级量子化学培训班（http://www.keinsci.c ...

好的，多谢sob老师指导

wza1477 · 发表于 Post on 2025-11-20 08:23:36

wal 发表于 2025-11-19 14:46
guess=read读取当前chk文件中的波函数，做stable=opt，然后再guess=read读取收敛波函数进行TDDFT计算

好的，谢谢您！

sobereva · 发表于 Post on 2025-11-19 23:54:10

不要写毫无意义的test关键词。direct也是多余的

参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）讲电子激发部分的幻灯片：

wal · 发表于 Post on 2025-11-19 14:46:10

wza1477 发表于 2025-11-19 13:55
感谢您的解答，请问有什么好的解决办法嘛？

guess=read读取当前chk文件中的波函数，做stable=opt，然后再guess=read读取收敛波函数进行TDDFT计算

wza1477 · 发表于 Post on 2025-11-19 13:55:55

wal 发表于 2025-11-19 11:42
基态波函数不稳定

感谢您的解答，请问有什么好的解决办法嘛？

wal · 发表于 Post on 2025-11-19 11:42:44

基态波函数不稳定

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

[新手求助] 计算三重激发态能量时，T1数值为负值

浏览过的版块