24小时热门版块排行榜    

查看: 4389  |  回复: 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的回帖

zhurl123

新虫 (正式写手)

引用回帖:
8楼: Originally posted by Jasminer at 2012-04-18 18:47:57:
lz研究的混合物是在气相沉积的条件下得到的,那么情况可能会更复杂一些,因为此时得到的不见得就会是最稳定结构,据我所知一些CVD膜放置后有明显的相分离或晶化倾向。那么最保险的做法还是根据实验中相关表征结果 ...

老兄考虑的很周到,考虑到了结合实际的实验情况来建模,严格来讲,确实是在气相情况下,也很有可能得到的是非晶薄膜,如果用MD跑真的可能要复杂很多。

我把一楼中图中的计算细节贴出来,老兄再给点意见,我觉得他们应该是没考虑那么多的因素,还望老兄指点指点:

To identify the origins of the gap states found in valence
band, theoretical simulation with density functional theory
methods were used to compute the molecular orbitals
(MOs) and simulated density of states. The optimized geometric structures were obtained with the method of B3-
LYP and the basis set of 6-31G, and the single-point energy
was then applied to retrieve the detail information of
molecular energy using Gaussian03 software [17]. The
simulated valence-band spectra using calculated MOs
was extracted via GaussSum [18] with full-width-at-halfmaximum of 0.4 eV.

计算后的结果就是图一的结果,最上面和最下面的图是单个有机分子的,中间那个是两个一起计算后的结果。
10楼2012-04-19 08:51:23
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 13/650 2026-08-19 14:33 by 陈晨晨陈啊
[基金申请] filecode,4个jtjc了 +8 ziyangfang 2026-08-19 9/450 2026-08-19 13:08 by blacksandy
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
信息提示
请填处理意见