Tpssh泛函优化激发态出错。

小范范1989

我用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计算，我用这个泛函在优化计算，还在计算着，临时没有出什么错误。
谢谢指点

sobereva

和nosymm无关
和TPSS在程序里实现可能有关，也许是bug，所以让你换泛函再试

小范范1989

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

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

小范范1989

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老师的指点。

sobereva

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

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

小范范1989

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

嗯，谢谢sob老师的指点，

小范范1989

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老师指点。

sobereva

小范范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短程是杂化泛函，其杂化比例有参数另设）

Tznameswu

请问小范范老师解决了tpssh问题吗？求教，感谢

小范范1989

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

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

Tznameswu

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

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

小范范1989

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

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

小范范1989

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

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

Tznameswu

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

好的，谢谢老师

shuifangren

老师，我想问一下tpssh在高斯中如何输入呢，没有找到相关的内容。。