24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1827  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

帅真帅真

新虫 (初入文坛)

[求助] [b]求助:做两个轨道的电子密度差图[/b]已有2人参与

看到了不少做差分密度图帖子,但都是激发态分子与基态分子整体的差。我想做两个轨道单独的电荷密度差图,不知道怎么做。求教各位大神~~~
比如说8→9轨道,这个图怎么做。先谢过各位
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

引用回帖:
10楼: Originally posted by 帅真帅真 at 2015-11-02 08:58:19
亲,能不能给一个后者的具体操作步骤就是激发态的那个,拜托~~~...

优化计算分子M1,多重度为1,作为开壳层算基态,加关键词output=wfn,得M1.wfn。
用M1的分子坐标,计算你要求的激发态(多重度为2)的单点能,加关键词nosymm,得分子M2,得M2.wfn。

将M1.wfn文件输入multiwfn,删去其它MO,只保留MO8,存一个MO8.wfn。
再将M2.wfn文件输入multiwfn,删去其它MO,只保留MO9,存一个MO9.wfn。
最后进入multiwfn,作MO9.wfn与MO8.wfn电子密度差。
这个电子密度差可以看到从M1到激发态M2,电子分布发生了什么变化。
更进一步可以在multiwfn中,定量这个电子变化的数值。

当然激发态并不是这么简单地只是某个轨道电子的跃迁,这里只是简单计算例子而已。只是对应MO之电子密度差的算例而已。

特别注意:
GS09 A版本,M1.wfn文件输入multiwfn后,只出现占有电子的MO。Alpha电子MO排完再排Beta电子MO。
而GS09 D版本,会将全部空轨道也列入,使轨道几倍增长,难以辨识。
注意M1之MO8与M2之MO9的对应与辨认。
12楼2015-11-02 15:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

帅真帅真

新虫 (初入文坛)

我是这样子的思路:用multiwfn算出8/9轨道的cub文件,然后放到高斯09里面去减。但是用multiwfn生成cub文件时会出现“input the orbital number,should within 1-8",然后9轨道的cub文件出不来,求助大神
2楼2015-10-30 11:21:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

保存wfn文件或molden文件,复制两份,一个只保留轨道8,一个只保留轨道9,然后在multiwfn里相减
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2015-10-30 13:10:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

multiwfn就可以,有相关的帖子可以看看。一搜就知道了
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
4楼2015-10-30 14:21:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见