计算化学公社

标题: 关于绘制激发态的静电势着色的分子范德华表面图 [打印本页]

作者
Author: 量化小菜鸡 时间: 2023-9-15 01:29
标题: 关于绘制激发态的静电势着色的分子范德华表面图
各位同学老师好。最近有一些需要，想去做关于激发态（比如S1，S2）的静电势图。通过社长这个博文：http://sobereva.com/403得知，就算用了density的关键词，fchk文件里GTF信息记录的也只是基态的，所以假如直接用来静电势，也得不到激发态的。所以我在想，假如用这个博文的办法，生成含有激发态信息的new.mwfn文件，那么再用这个文件去算静电势是不是就得到激发态的了。
不知道理解的对不对，希望大家指正。

作者
Author: sobereva 时间: 2023-9-15 01:50
TDDFT计算直接写out=wfn产生的wfn文件记录的就是激发态波函数的自然轨道（如果是G09老版本，还需要同时加density关键词），直接用这个文件作为Multiwfn的各种静电势分析的输入文件即可考察激发态的静电势。

作者
Author: 量化小菜鸡 时间: 2023-9-15 02:26

sobereva 发表于 2023-9-15 01:50
TDDFT计算直接写out=wfn产生的wfn文件记录的就是激发态波函数的自然轨道（如果是G09老版本，还需要同时加de ...

谢谢社长，但假如我已经有了激发态的fchk文件，是不是这样子转一圈去转换成mwfn结果也是没问题的。

作者
Author: zjxitcc 时间: 2023-9-15 10:22

量化小菜鸡发表于 2023-9-15 02:26
谢谢社长，但假如我已经有了激发态的fchk文件，是不是这样子转一圈去转换成mwfn结果也是没问题的。

不如先讲讲你这个“激发态的fchk文件”怎么产生的，是不是有误

作者
Author: 量化小菜鸡 时间: 2023-9-15 14:42

zjxitcc 发表于 2023-9-15 10:22
不如先讲讲你这个“激发态的fchk文件”怎么产生的，是不是有误

就是正常写TD density=current得到fchk, 然后再用那个博文的方法，得到激发态的new.mwfn文件，但这个文件不也是根据之前得到的fchk文件里来的，我也不能说那个fchk文件里的信息完全是基态的吧。

作者
Author: sobereva 时间: 2023-9-16 01:18

量化小菜鸡发表于 2023-9-15 14:42
就是正常写TD density=current得到fchk, 然后再用那个博文的方法，得到激发态的new.mwfn文件，但这个文件 ...

Multiwfn载入mwfn文件，进主功能6，选4 Print detail information of an orbital，如果看到轨道占据数不都是整数，说明存的是自然轨道，可以用于分析激发态波函数

作者
Author: 量化小菜鸡 时间: 2023-9-16 01:54

sobereva 发表于 2023-9-16 01:18
Multiwfn载入mwfn文件，进主功能6，选4 Print detail information of an orbital，如果看到轨道占据数不 ...

好的，谢谢社长

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