24小时热门版块排行榜    

CyRhmU.jpeg
查看: 470  |  回复: 5
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 agou8888 的 13 个金币

agou8888

新虫 (小有名气)

[交流] 【求助】非d10过渡金属计算吸收谱

我本人不是做计算的,在这里向大家请教个问题,如果问题比较初级,请大家见谅

我一个含有Mn2+的配合物,最近用gaussian算了下它的吸收谱:
过程如下:
1:取晶体结构数据,首先进行基态结构优化
2:利用TDDFT  计算激发态的吸收光谱 命令如下(因是开壳层,所以采取默认的飞非限制开壳层计算):
       #p b3lyp/genecp TD(singlets) scfcyc=400 pop=full  

计算部分结果为:
Excitation energies and oscillator strengths:

Excited State   1:   ?Spin  -A      3.2240 eV  384.56 nm  f=0.0000
    175A ->181A        0.12072
    175A ->183A        0.10199
    178A ->180A       -0.12302
    178A ->184A        0.13707
    178A ->186A       -0.14122
    165B ->172B        0.14724
    165B ->174B        0.15046
    168B ->170B        0.79040
    168B ->171B       -0.36040
    168B ->173B        0.15314
    168B ->176B        0.10037
    168B ->182B        0.14175
    168B ->184B       -0.13344
This state for optimization and/or second-order correction.
Total Energy, E(RPA) =  -3789.74279082   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:   ?Spin  -A      3.2327 eV  383.53 nm  f=0.0003
    175A ->184A       -0.11131
    175A ->186A        0.11716
    178A ->181A       -0.31120
    178A ->183A       -0.26519
    165B ->170B       -0.21003
    165B ->171B        0.14634
    165B ->182B       -0.10512
    168B ->169B        0.43824
    168B ->172B       -0.42163
    168B ->174B       -0.46742
    168B ->175B        0.19210
    168B ->179B        0.20287
    168B ->183B        0.10797

Excited State   3:   ?Spin  -A      3.3171 eV  373.78 nm  f=0.0013
    176A ->181A       -0.19235
    176A ->183A       -0.15044
    177A ->180A       -0.16484
    177A ->184A        0.15798
    177A ->186A       -0.16716
    178A ->179A       -0.12872
    165B ->170B        0.10625
    166B ->172B        0.19362
    166B ->174B        0.20177
    166B ->175B       -0.10672
    166B ->179B       -0.11576
    167B ->170B        0.65511
    167B ->171B       -0.31386
    167B ->173B        0.12561
    167B ->182B        0.12734
    167B ->184B       -0.11933
    168B ->169B        0.31394

这里有2个问题:
1:我的计算过程有没有问题,到底gaussian03可不可以计算开壳层体系?
2:我只知道分子轨道分析中的A,B分别表示alpha和beta电子,我具体分析电子跃迁时是以alpha的为准还是beta或者两者都考虑。这点我很是不明白,还请专家尽可能提供与之相关的资料。

[ Last edited by agou8888 on 2010-3-24 at 22:53 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)

不要用QQ问我东西

仔细读读高斯说明书中关于TD的说明。
Gaussian.com就有说明书
Chemistry[]==[]Chem[]is[]try!!!
2楼2010-03-24 23:05:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agou8888

新虫 (小有名气)

引用回帖:
Originally posted by abbott at 2010-03-24 23:05:36:
仔细读读高斯说明书中关于TD的说明。
Gaussian.com就有说明书

说明书我是看过,TD针对闭壳体系,但是我最近查了些文献,对过度金属做TD的还真有不少。故到这里问个究竟。不知道专家您什么意思?
3楼2010-03-24 23:23:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

agou8888(金币+5):非常感谢,说实话你的解答有一些术语我都不懂,但是我会查阅文献理解的 2010-03-25 00:20
agou8888(金币+1):能否提供下你说的第一种方法有关的文献?非常感谢 2010-03-25 08:34
gaussian可以做开壳层tddft,同时忽略singlets、triplets、50-50选项。激发态的自旋多重度和基态相同,也就是说,得不到自旋反转态。

也可以加、减一个电子使之变成闭壳层,然后用sac-ci计算闭壳层参考态能量,并计算减、加一个电子后的离子态激发能。这个计算是CCSD级别的,你可以做个闭壳层参考态的CCSD计算,估计一下计算量。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2010-03-25 00:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shalene

至尊木虫 (著名写手)

agou8888(金币+1):非常感谢您的建议! 2010-03-25 22:07
gaussian不适合算开壳金属配合物,用adf吧
5楼2010-03-25 14:07:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agou8888

新虫 (小有名气)

此贴到此为止,谢谢大家!
6楼2010-03-25 22:09:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 agou8888 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见