24小时热门版块排行榜    

查看: 1892  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿河工大 081700 276求调剂 +4 地球绕着太阳转 2026-03-23 4/200 2026-03-26 14:27 by zzll406
[考研] 0856求调剂 +8 zhn03 2026-03-25 9/450 2026-03-26 13:42 by zzll406
[考研] 308求调剂 +5 墨墨漠 2026-03-25 5/250 2026-03-25 22:19 by 544594351
[考研] 材料与化工 322求调剂 +6 然11 2026-03-19 6/300 2026-03-25 18:37 by haxia
[考研] 各位老师您好:本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 300求调剂,材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
请填处理意见