计算化学公社

标题: 求助激发态优化给出的激发能和用优化的结构单独做TD得到的激发能不同 [打印本页]

作者
Author: dreamyeye 时间: 2021-11-6 22:56
标题: 求助激发态优化给出的激发能和用优化的结构单独做TD得到的激发能不同
一个配合物体系，优化S2，td(nstates=5,root=2) opt 最后一步给出的激发能量和把优化后的结构拿出来单独做一步td td(nstates=5,50-50)，得到的同一个态的激发能不同，哪个更准一些呢？

作者
Author: wzkchem5 时间: 2021-11-6 23:08
是不是加了溶剂模型？

作者
Author: pwzhou 时间: 2021-11-6 23:09
气相的话，两者应该完全相同。液相的话，如果用的Gaussian的默认参数，两者必然不同。把你的详细参数和结果列出来看一下。

作者
Author: dreamyeye 时间: 2021-11-6 23:14

wzkchem5 发表于 2021-11-6 23:08
是不是加了溶剂模型？

是的。

作者
Author: dreamyeye 时间: 2021-11-6 23:16

pwzhou 发表于 2021-11-6 23:09
气相的话，两者应该完全相同。液相的话，如果用的Gaussian的默认参数，两者必然不同。把你的详细参数和结果 ...

前者#p opt=gdiis td=(nstates=5,root=2)/genecp scrf=(smd,solvent=dichloromethane) pbe1pbe
后者 #p td=(nstates=20,50-50)/genecp scrf=(smd,solvent=dichloromethane) pbe1pbe
两者基于同一个分子结构，后者用的是前者优化之后的结构
有没有可能和50-50有关？

作者
Author: wzkchem5 时间: 2021-11-7 04:26

dreamyeye 发表于 2021-11-6 16:16
前者#p opt=gdiis td=(nstates=5,root=2)/genecp scrf=(smd,solvent=dichloromethane) pbe1pbe
后者 #p ...

参见http://sobereva.com/314，2.5节

作者
Author: dreamyeye 时间: 2021-11-7 11:18

wzkchem5 发表于 2021-11-7 04:26
参见http://sobereva.com/314，2.5节

谢谢。

作者
Author: pwzhou 时间: 2021-11-7 21:53

dreamyeye 发表于 2021-11-6 23:16
前者#p opt=gdiis td=(nstates=5,root=2)/genecp scrf=(smd,solvent=dichloromethane) pbe1pbe
后者 #p ...

后面那一步在TD里面加上Eqsolv，就和前面的结果一致了，TD(Eqsolv,nstates=20,50-50)

作者
Author: dreamyeye 时间: 2021-11-8 07:41

pwzhou 发表于 2021-11-7 21:53
后面那一步在TD里面加上Eqsolv，就和前面的结果一致了，TD(Eqsolv,nstates=20,50-50)

谢谢。

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