求助激发态优化给出的激发能和用优化的结构单独做TD得到的激发能不同

dreamyeye

一个配合物体系，优化S2，td(nstates=5,root=2) opt   最后一步给出的激发能量和把优化后的结构拿出来单独做一步td  td(nstates=5,50-50)，得到的同一个态的激发能不同，哪个更准一些呢？

wzkchem5

是不是加了溶剂模型？

pwzhou

气相的话，两者应该完全相同。液相的话，如果用的Gaussian的默认参数，两者必然不同。把你的详细参数和结果列出来看一下。

dreamyeye

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

是的。

dreamyeye

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有关？

wzkchem5

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节

dreamyeye

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

谢谢。

pwzhou

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)

dreamyeye

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

谢谢。