24小时热门版块排行榜    

查看: 2206  |  回复: 14

wbn

银虫 (小有名气)

[求助] 多分子体系的体系的激发态计算--第一激发态能量特别低? 已有1人参与

大家好,我想计算两个分子间的电子转移,其中一个分子是闭壳层的电子受体,电子转移的过程是用光先将受体HOMO上的电子激发的LUMO上,然后供体上的电子填充到空出来的受体LUMO中。我先算了一下受体的第一激发态能量,计算的结果大概是在410 nm 左右,这与UV-vis测得的结果非常接近。但是将供体分子放入后,整个体系的LUMO和HOMO通过GaussView看起来和受体的LUMO和HOMO一模一样,所以判断体系的第一激发态就是受体的local excited state. 但是这个激发态能量非常低,只有0.4 eV,显得不太正常。第二激发态倒很像是电子转移态,我算了一下,供体上的电子密度减少了很多。然后我试着将两个分子间距离拉大了好几纳米,结果是第一激发态变成只有0.01eV了? 请问这是怎么回事呢?整个体系的构型我是从分子模拟的盒子中提取的。另外,我的供体有一个负电荷,我不太确定它是开壳层还是闭壳层,所以我multiplicity尝试了1和3,而且激发态我也算了singlets和triplets,结果虽然不一样,但是这个问题仍然存在。我是新手,如果哪位大侠指点迷津,不胜感激!

我的输入差不多是这样的:
# b3lyp/6-31g(d) td=(nstates=20,root=1,singlets) density=current

输出差不多是这样的:
Excitation energies and oscillator strengths:

Excited State   1:      Singlet-A      0.4050 eV 3061.25 nm  f=0.0394  <S**2>=0.000
      76 -> 77         0.74388
      76 <- 77        -0.23195
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -964.453011066
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      1.1760 eV 1054.26 nm  f=0.0018  <S**2>=0.000
      75 -> 77         0.70649

Excited State   3:      Singlet-A      2.1784 eV  569.15 nm  f=0.0127  <S**2>=0.000
      76 -> 78         0.70509
回复此楼

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

wbn

银虫 (小有名气)

wbn: 回帖置顶 2015-10-29 01:27:48
找到原因了,原来是我的输入文件是我用代码生成的,然后写到指定chk文件时的那一行时处理出了点问题,在等于号后面多了个空格,结果出来的chk文件就在给我瞎写...我去掉空格重新算了一遍,发现体系的HOMO是donor的HOMO,而LUMO是acceptor的LUMO,这样以上结论就没有问题了...之前说的没有电子转移可能是自旋不对,我把体系设成三重态以后发现第一激发态有0.8个电子转移了...

给高斯的输入文件自由度跪了,我把金币分分大家散了吧...
13楼2015-10-29 01:25:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

virtualzx

木虫 (著名写手)

电子数写对了吗,结构有没有画错

发自小木虫IOS客户端
2楼2015-10-28 06:18:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)

引用回帖:
2楼: Originally posted by virtualzx at 2015-10-28 06:18:31
电子数写对了吗,结构有没有画错

体系电荷是-1不会错。结构是我从MD的结果中提取的,由于两个分子正在撞击所以结构有一点扭曲,但是我把供体分子去掉以后再算第一激发态,和opt之前的410nm差别不大,应该不是这个问题
3楼2015-10-28 06:42:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

供体分子的结构呢?供体第一激发态应该是多少,在这个结构又是多少?

发自小木虫IOS客户端
4楼2015-10-28 06:45:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)

引用回帖:
4楼: Originally posted by virtualzx at 2015-10-28 06:45:37
供体分子的结构呢?供体第一激发态应该是多少,在这个结构又是多少?

我没算,等我算算
5楼2015-10-28 07:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

我是菜鸟,下面是我的一点讨论
1:但是将供体分子放入后,整个体系的LUMO和HOMO通过GaussView看起来和受体的LUMO和HOMO一模一样,所以判断体系的第一激发态就是受体的local excited state。这个地方,按说受体的L是整体体系的L,供体的H是整个体系的H,为什么你做出来会在一个分子上?是不是这个问题的原因。
2:但是这个激发态能量非常低,只有0.4 eV,显得不太正常。我怀疑你这个激发能是不是供体和受体的H之间的跃迁了?可能不对哈,我只是感觉。
3:咨询楼主一个问题:你是通过马修斯方程计算分子间的电荷转移速率?我知道这个方程,但是里面具体的参数,应该是怎么计算?谢谢楼主指点
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
6楼2015-10-28 07:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
wbn: 金币+20, ★★★很有帮助 2015-10-28 22:13:03
wbn: 金币+80, ★★★★★最佳答案 2015-10-29 01:27:03
“将供体分子放入后,整个体系的LUMO和HOMO通过GaussView看起来和受体的LUMO和HOMO一模一样,所以判断体系的第一激发态就是受体的local excited state. 但是这个激发态能量非常低,只有0.4 eV,显得不太正常。”

1、除了“形状”,LUMO、HOMO的能量是否和单一受体基本一致?
2、第一激发态是否对应于HOMO-LUMO跃迁?
7楼2015-10-28 11:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

“然后我试着将两个分子间距离拉大了好几纳米,结果是第一激发态变成只有0.01eV了?”

这个很诡异。。。
8楼2015-10-28 11:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)

引用回帖:
4楼: Originally posted by virtualzx at 2015-10-28 06:45:37
供体分子的结构呢?供体第一激发态应该是多少,在这个结构又是多少?

您好,谢谢帮助。我算得供体第一激发态是6.3个eV(190 nm),在这个结构是6.5 eV,看起来差不多,因为供体分子在可见光范围内确实没有吸收。我的供体分子是dicyanamide anion, 受体是9,10-dicyanoanthracene, 所以基本体系的LUMO往上几个轨道都是受体的。
9楼2015-10-28 21:37:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)

引用回帖:
6楼: Originally posted by 小范范1989 at 2015-10-28 07:28:50
我是菜鸟,下面是我的一点讨论
1:但是将供体分子放入后,整个体系的LUMO和HOMO通过GaussView看起来和受体的LUMO和HOMO一模一样,所以判断体系的第一激发态就是受体的local excited state。这个地方,按说受体的L是 ...

1. 照理来说应该是这样,我猜是不是因为两个分子没有摆放到恰好的位置,等我做一做别的构象看。也许这就是问题的关键。

2. 应该不太可能吧,不然我将两个分子拉开几纳米就不应该有作用了。

3. 是的。 Marcus方程里面总共有三个参数,reorganization energy 的公式就是Marcus理论的一部分,free  energy 可以用电化学的方法测出,见 http://goldbook.iupac.org/GT07388.html , 第三个electron coupling matrix element就是我现在想算的东西,用generalized mulliken-hush theory, 当然对于我这个体系得用一种改进的方法。 你具体是什么体系? 不同的体系方法千差万别,比如我的溶液双分子体系还得考虑diffusion-reaction的相互作用以及Zusman的solvent relaxation time.
10楼2015-10-28 21:52:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wbn 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 9/450 2026-08-04 07:46 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 06:07 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见