计算化学公社

标题: Tpssh泛函优化激发态出错。 [打印本页]

作者
Author: 小范范1989 时间: 2015-6-25 07:07
标题: Tpssh泛函优化激发态出错。
我用tpssh泛函，（需要10%的HF）  计算了分子的td，能正确结束，同时我用tpssh优化s1和t1，出现下面的错误，
Excitation energies and oscillator strengths:

Excited State 1:    Singlet-A"    3.2510 eV  381.38 nm  f=0.7139  <S**2>=0.000
   145 ->146       0.70501
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2046.78951508
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A'    3.5137 eV  352.86 nm  f=0.1131  <S**2>=0.000
   144 ->146       0.69519
   145 ->149       0.11226

Excited State 3:    Singlet-A"    3.7195 eV  333.34 nm  f=0.0092  <S**2>=0.000
   144 ->148       -0.12693
   145 ->147       0.68706
SavETr:  write IOETrn= 770 NScale= 10 NData=  16 NLR=1 LETran=    64.
The selected state is a singlet
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass= 1 NMax=    1.
CISAX will form    1 AO SS matrices at one time.
NMat=    1 NSing=    1.
GePol: Maximum number of non-zero 1st derivatives =    193
Invalid value of MaxDer in TPSSx
Invalid value of MaxDer in TPSSx
Error termination via Lnk1e in /home/software/g09/l914.exe at Tue Jun 16 12:34:30 2015.
Error termination via Lnk1e in /home/software/g09/l914.exe at Tue Jun 16 12:34:30 2015.
Job cpu time:  0 days  7 hours 58 minutes 27.0 seconds.
File lengths (MBytes):  RWF= 1180 Int=    0 D2E=    0 Chk=    27 Scr=    1
是泛函的问题，还是没有使用nosymm的问题，我的td文件中没有加nosymm，能正常结束。谢谢指点
同时，当时在群里问了一句sob老师，老师说用近似的03lyp计算，我用这个泛函在优化计算，还在计算着，临时没有出什么错误。
谢谢指点

作者
Author: sobereva 时间: 2015-6-25 08:27
和nosymm无关
和TPSS在程序里实现可能有关，也许是bug，所以让你换泛函再试

作者
Author: 小范范1989 时间: 2015-6-25 08:36

sobereva 发表于 2015-6-25 08:27
和nosymm无关
和TPSS在程序里实现可能有关，也许是bug，所以让你换泛函再试

谢谢sob老师，等03lyp计算完了，我在和老师汇报结果。

作者
Author: 小范范1989 时间: 2015-7-1 10:38

sobereva 发表于 2015-6-25 08:27
和nosymm无关
和TPSS在程序里实现可能有关，也许是bug，所以让你换泛函再试

sob老师你好，现在用o3lyp计算出来了，但是，我看o3lyp优化s1，计算的荧光波长结果为493nm。然而实验上是402nm。我感觉是不是要是换位tpssh，计算的荧光波长应该更接近于实验值，
所以我的问题是：
1：根据老师您的经验，文章中能说493nm个实验值符合的好吗？
2：有没有其他的方法来实现10%HF，我这边的机器用tpssh优化s1出错，我让别人帮我试了一下，他们那边的机器也是这个错误。可能真是个bug？
3：谢谢sob老师的指点。

作者
Author: sobereva 时间: 2015-7-1 11:07

小范范1989 发表于 2015-7-1 10:38
sob老师你好，现在用o3lyp计算出来了，但是，我看o3lyp优化s1，计算的荧光波长结果为493nm。然而实验上是 ...

计算结果应该更蓝一些才符合实验，所以HF成份需要更大。TPSSh比O3LYP的HF成份还小，肯定不合适。应该尝试B3LYP或PBE0

作者
Author: 小范范1989 时间: 2015-7-1 17:28

sobereva 发表于 2015-7-1 11:07
计算结果应该更蓝一些才符合实验，所以HF成份需要更大。TPSSh比O3LYP的HF成份还小，肯定不合适。应该尝试 ...

嗯，谢谢sob老师的指点，

作者
Author: 小范范1989 时间: 2015-9-2 20:15

sobereva 发表于 2015-6-25 08:27
和nosymm无关
和TPSS在程序里实现可能有关，也许是bug，所以让你换泛函再试

sob老师好，最近我又用pbe1kcis优化激发态出错： KCIS cannot do third derivatives.
是不是我的这个高斯也不支持pbe1kcis？

还一个问题就是，我最近看了一个文献，利用的是LC-wpbe来计算，默认的w是0.4？然后我看他是选择了不同的w来计算，我想咨询一下sob老师，这个w怎么调节？比如说我要用0.2的，
老师曾经给我的网站http://comp.chem.umn.edu/info/DFT.htm#mpw1k 我看过了，没有这个泛函。谢谢sob老师指点。

作者
Author: sobereva 时间: 2015-9-3 00:34

小范范1989 发表于 2015-9-2 20:15
sob老师好，最近我又用pbe1kcis优化激发态出错： KCIS cannot do third derivatives.
是不是我的这个高 ...

你就别用相关泛函是KCIS的泛函了。

用IOp可以设，这是介绍：
IOp(3/107) 设定短程和长程HF交换项。0时使用标准HF交换项，MMMMMNNNNN时短程用NNNNN/10000的HF交换而长程用MMMMM/10000的HF交换。这两个数对应于输出文件中HFx后面的wShort和wLong。
IOp(3/108) 设定短程和长程DFT交换项。0时使用标准DFT交换项，MMMMMNNNNN时短程用NNNNN/10000的DFT交换而长程用MMMMM/10000的DFT交换。这两个数对应于输出文件中DFx后面的wShort和wLong。
（注：上面两个选项字面上的意思有点误导性。实际上，将3/107和3/108都设为MMMMM00000，就相当于用了ω参数为MMMMM/10000的范围分离泛函。对GGA加LC-时ω默认是0.47/Bohr，等价于0470000000。wB97X默认是ω=0.3，等价于0300000000，但注意wB97X短程是杂化泛函，其杂化比例有参数另设）

作者
Author: Tznameswu 时间: 2020-6-20 20:44
请问小范范老师解决了tpssh问题吗？求教，感谢

作者
Author: 小范范1989 时间: 2020-6-21 06:40

Tznameswu 发表于 2020-6-20 20:44
请问小范范老师解决了tpssh问题吗？求教，感谢

你说。按照之前说的，好像是不能优化激发态。

作者
Author: Tznameswu 时间: 2020-6-28 20:44

小范范1989 发表于 2020-6-21 06:40
你说。按照之前说的，好像是不能优化激发态。

我用的高斯09D，和您的版本相近，应改也是Tpssh无法计算激发态了。
看到o3lyp的HF成份和它相近，之前按照手册文献修改了b3lyp的参数，但是结果报错了，可能是我o3lyp的输入文件写错了？想问下小范范老师的o3lyp输入文件是如何写的呀？谢谢

作者
Author: 小范范1989 时间: 2020-6-29 07:37

Tznameswu 发表于 2020-6-28 20:44
我用的高斯09D，和您的版本相近，应改也是Tpssh无法计算激发态了。
看到o3lyp的HF成份和它相近，之前按 ...

应该直接就写O3LYP吧https://comp.chem.umn.edu/info/DFT.htm#mpw1k

作者
Author: 小范范1989 时间: 2020-6-29 07:45

Tznameswu 发表于 2020-6-28 20:44
我用的高斯09D，和您的版本相近，应改也是Tpssh无法计算激发态了。
看到o3lyp的HF成份和它相近，之前按 ...

O3LYP：直接用，高斯中有。但是用O3LYP来优化激发态出错，是g09D01的bug。所以需要换A01版本。

作者
Author: Tznameswu 时间: 2020-6-29 09:30

小范范1989 发表于 2020-6-29 07:45
O3LYP：直接用，高斯中有。但是用O3LYP来优化激发态出错，是g09D01的bug。所以需要换A01版本。

好的，谢谢老师

作者
Author: shuifangren 时间: 2021-3-5 10:01
老师，我想问一下tpssh在高斯中如何输入呢，没有找到相关的内容。。

作者
Author: snljty 时间: 2021-3-5 10:15

shuifangren 发表于 2021-3-5 10:01
老师，我想问一下tpssh在高斯中如何输入呢，没有找到相关的内容。。

https://gaussian.com/dft/
点Keywords: Hybrid Functionals
搜TPSSh
直接写TPSSh即可。
这个问题和本贴无直接关系，建议这类问题直接发到http://bbs.keinsci.com/thread-806-1-1.html

作者
Author: shuifangren 时间: 2021-3-5 10:28

snljty 发表于 2021-3-5 10:15
https://gaussian.com/dft/
点Keywords: Hybrid Functionals
搜TPSSh

好的，谢谢啦，我下次会注意~

作者
Author: GYJ2580 时间: 2025-3-10 09:36

shuifangren 发表于 2021-3-5 10:01
老师，我想问一下tpssh在高斯中如何输入呢，没有找到相关的内容。。

这个问题找到答案了吗

作者
Author: zjxitcc 时间: 2025-3-10 11:06

GYJ2580 发表于 2025-3-10 09:36
这个问题找到答案了吗

不用找答案，就写作TPSSh

作者
Author: GYJ2580 时间: 2025-3-11 15:26

zjxitcc 发表于 2025-3-10 11:06
不用找答案，就写作TPSSh

好的谢谢

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