计算化学公社

标题: 计算三重激发态能量时，T1数值为负值 [打印本页]

作者
Author: wza1477 时间: 2025-11-19 11:18
标题: 计算三重激发态能量时，T1数值为负值
各位老师好，本人在计算一个分子的三重激发态能量时，最终得到的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

作者
Author: wal 时间: 2025-11-19 11:42
基态波函数不稳定

作者
Author: wza1477 时间: 2025-11-19 13:55

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

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

作者
Author: wal 时间: 2025-11-19 14:46

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

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

作者
Author: sobereva 时间: 2025-11-19 23:54
不要写毫无意义的test关键词。direct也是多余的

参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）讲电子激发部分的幻灯片：
(, 下载次数 Times of downloads: 0)

作者
Author: wza1477 时间: 2025-11-20 08:23

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

好的，谢谢您！

作者
Author: wza1477 时间: 2025-11-20 08:24

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

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

好的，多谢sob老师指导

作者
Author: iamjjLin 时间: 2025-11-20 11:27
建议先仔细检查结构，特别是H原子，多了少了都会导致激发态为负。结构确定没问题后再尝试稳定波函数。

作者
Author: wza1477 时间: 2025-11-20 13:59

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

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

作者
Author: wza1477 时间: 2025-11-20 14:32
本帖最后由 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

作者
Author: wza1477 时间: 2025-11-20 14:56

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

作者
Author: sobereva 时间: 2025-11-20 21:07

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

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

作者
Author: wza1477 时间: 2025-11-20 22:10

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

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

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