24小时热门版块排行榜    

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

zhurl123

新虫 (正式写手)

[交流] 高斯计算两个有机分子 已有4人参与

请教下做量子化学的高手们

如何实现两个有机分子的相互作用计算,目的是算出calculated energy levels或者是density of state

比如像下面的图,图中最下面和最上面的曲线和短线图分别是单独的分子A和分子B的能级和DOS, 这个很简单就可以算出来,中间的是两个分子一起作用后的结果,这咋计算呀

请高手指教


[ Last edited by zhurl123 on 2012-4-16 at 16:01 ]
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
送鲜花一朵
小红豆: 金币+3, 谢谢回帖 2012-04-18 18:34:05
不是随便画出来放在一起,然后让他们自动优化就行了。

还是要从两个分子的结构特点搭配。特别是二者可能形成配键、氢键的可能部位,以此为根据来搭配。否则只是范德华作用。
进一步则要看A、B的HOMO、LUMO图形,来比对。看谁的HOMO图形会与另外的LUMO图形相搭配,这时要比较轨道对称性。

在这样的基础上再优化。
7楼2012-04-18 16:49:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

Jasminer

铁杆木虫 (著名写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
zhou2009: 金币+2 2012-04-18 16:26:52
这个关键是你要根据实际情况和说明问题的需求建立合适的混合物模型,算起来和单个分子没区别。
3楼2012-04-18 01:10:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhurl123

新虫 (正式写手)

引用回帖:
3楼: Originally posted by Jasminer at 2012-04-18 01:10:08:
这个关键是你要根据实际情况和说明问题的需求建立合适的混合物模型,算起来和单个分子没区别。

感谢楼上的老兄,建议很好。

对两种有机物分别而言,通过计算单个有机物的态密度,可以知道哪些地方有轨道,上图中最上面和最下面的单竖线和曲线是高斯计算得到的轨道和画出的相应的态密度,从图上可以看出,在费米能级到1.8eV没有能级(对单个的SubPc和NPB而言).

实际问题中,两种有机分子同时在高真空下蒸发沉积,然后测定光电子能谱,这时候需要通过计算看看那哪些地方有新的轨道出现(光电子能谱和轨道应该是有很好的吻合)。上图中中间的曲线和bar,就是一起计算SubPc和NPB后的轨道和态密度,计算后发现,相对于单个分子,在费米能级到1eV之间出现了ABC三个新轨道,同时实验中光电子能谱也有峰出现,而对单个的SubPc和NPB的光电子能谱,在0到1.8 之间是没有能级的。

请教如何建模比较合适,比如放置两个分子时的距离,取向等以什么标准来衡量,然后让他们计算

感谢帮助!!!!

[ Last edited by zhurl123 on 2012-4-18 at 09:26 ]
4楼2012-04-18 09:22:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

necturnal

铜虫 (初入文坛)

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
zhou2009: 金币+3 2012-04-18 16:29:14
引用回帖:
4楼: Originally posted by zhurl123 at 2012-04-18 09:22:27:
感谢楼上的老兄,建议很好。

对两种有机物分别而言,通过计算单个有机物的态密度,可以知道哪些地方有轨道,上图中最上面和最下面的单竖线和曲线是高斯计算得到的轨道和画出的相应的态密度,从图上可以看 ...

这个方面的计算我也没做过。我就随便说一下。
当分子二聚或者高聚的时候,两个分子的分子轨道相互作用使得HOMO和LUMO轨道分裂。于是光谱就出现变化了。
关于两个分子之间相互作用的取向,一般情况下都不止一种取向,如果分子结构比较简单刚性,可以根据分子间相互作用力来作二聚体的构想。如果想比较完整,就需要做Anealing--分子淬火来寻找二聚体的最低能量构想。
PS:你说的density of state,我不太懂到底是什么意思,我平时也有算density of state,不过我是算动力学RRKM的时候用,这个是统计方法算出来的,高斯频率分析可以提供参数,但是没法算density of state。如果你所说的density of state和我说的意思一样的话,你就得需要采用统计学方法,而且你这个情况很可能需要做一个温度下的density of state,这样就是做该温度下的boltzman温度分布的态的积分,感觉比较复杂。不过我也没做过类似的,不敢保证我说的是不是正确
5楼2012-04-18 14:15:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:33 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ph7bxbho8f 2026-07-09 3/150 2026-07-10 05:53 by j4frgf6zcw
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[基金申请] 祈祷青基必中 50+3 hadengry 2026-07-09 9/450 2026-07-09 17:37 by 文杰猪猪
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 关于如何从代码看上不上会 +11 且听虎啸 2026-07-09 14/700 2026-07-09 15:27 by libeny
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评?应用基础研究和基础研究,希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[基金申请] 今年E04面上 +15 布布和一二 2026-07-05 21/1050 2026-07-09 10:25 by 学员fOzRO9
[硕博家园] 大龄残疾硕士的一点执念 +14 3776737374 2026-07-06 22/1100 2026-07-09 06:07 by lovechizi
[基金申请] 生命口会评 +4 wuke100666 2026-07-07 4/200 2026-07-08 16:38 by 逍遥不够自在
[教师之家] 网络举报学术不端成为了男女之间情杀的手段了,真荒唐。。。 +3 zju2000 2026-07-04 3/150 2026-07-07 21:00 by hyx552200
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联,生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
[论文投稿] 刑法学论文投稿求助 +3 高源678 2026-07-04 5/250 2026-07-06 12:20 by 平天尺
[论文投稿] Journal of Environmental Chemical Engineering +3 雪花天涯 2026-07-05 3/150 2026-07-06 09:00 by bobvan
[基金申请] 最后一年,祈求好运 +3 pyw1992 2026-07-06 3/150 2026-07-06 08:53 by 苏_frp
信息提示
请填处理意见