24小时热门版块排行榜    

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

faqianliu

银虫 (小有名气)

[交流] 【求助】关于adv mater 08年一篇论文中tddft计算的一些疑问,请高手指正 已有5人参与

关于adv mater 2008年一篇金属配合物td-dft计算的一些疑问

该论文下载地址如下:
http://www.brsbox.com/filebox/do ... 31e39a77b744bb53cbf

辅助材料信息下载地址如下:
http://www.brsbox.com/filebox/do ... 9674c91972a7f2c24d3

配合物1的cif文件下载地址:
http://www.brsbox.com/filebox/do ... 6ade2ec0968bcca8922

如不想麻烦下载,请参照后面的图形

此论文作者采用TD-DFT/PCM方式计算了三个3D金属簇合物的吸收光谱,并借此来解释其三阶非线性的差异,我只上载了配合物1的cif文件,我所关心的问题是:

1、配合物中cif的结构和计算模型并不相同,原始晶体信息如下图所示:


而论文中采用的模型为:


此处,模型和CIF文件明显不同,其中【WOS3Cu3】为+1价,也就是说,模型中的配合物电荷也不为0,计算无影响?

2、按照上图的模型,此处的自旋多重度应该为多少?从文献辅助材料提到的计算结构看,是按闭壳层计算的:
Results of Time-Dependent DFT Calculations: Cluster1-B3LYP-LanL2DZ TD-PCM (20states) in DMF (dielectric constant=36.7)
Excitation energies and oscillator strengths:
Excited State 1: Singlet-A 2.5930 eV 478.15 nm f=0.0004
101 ->106 -0.23472
102 ->106 -0.27424
103 ->104 0.57035
103 ->105 -0.13601
This state for optimization and/or second-order correction.
Total Energy, E(RPA) = -1616.83190592
Copying the excited state density for this state as the 1-particle RhoCI density.
Excited State 2: Singlet-A 2.6426 eV 469.18 nm f=0.0001
101 ->104 -0.16893
102 ->104 -0.37668
103 ->106 0.53940
Excited State 3: Singlet-A 2.7224 eV 455.42 nm f=0.0072
101 ->104 0.20670
102 ->104 0.47922
102 ->105 -0.10689
103 ->106 0.42014
Excited State 4: Singlet-A 2.7438 eV 451.87 nm f=0.0096
100 ->106 -0.21888
101 ->106 0.49112
102 ->106 0.26780
103 ->104 0.31098
Excited State 5: Singlet-A 2.7544 eV 450.14 nm f=0.0010
98 ->104 0.15400
98 ->105 0.62515
99 ->107 -0.15805
100 ->107 0.14758

辅助材料中配合物1给了20个态,我只列了前5个,显示均为单线态,可是我按照模型搭建的分子,按自旋多重度为1计算,根本无法计算,高斯出错:

The combination of multiplicity 1 and   207 electrons is impossible.


如果按开壳层计算的话,Excited State 1: Singlet-A 2.5930 eV 478.15 nm f=0.0004
计算结果中的singlet应该为?,也就是程序无法判断单线态或三线态,按作者提供的应该是闭壳层计算,且从文章的图形中也均以homo lumo来命名轨道,如下图:


也就是说,论文中采用的应该为闭壳层计算,作者是如何做到的?


以上2疑问,请指正

[ Last edited by yjcmwgk on 2010-3-25 at 09:55 ]
回复此楼

» 收录本帖的淘帖专辑推荐

小技巧

» 猜你喜欢

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

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

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

qzhaosdu

金虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流


[ Last edited by qzhaosdu on 2010-3-26 at 13:46 ]
8楼2010-03-26 13:08:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

abbott

金虫 (著名写手)

不要用QQ问我东西


小木虫(金币+0.5):给个红包,谢谢回帖交流
不懂 楼主写comment吧。
或者写信咨询作者?
Chemistry[]==[]Chem[]is[]try!!!
2楼2010-03-25 09:58:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+8):good 2010-04-01 17:53
作者在原文中已经明确说,他计算的是cluster。如果你看看这个cif文件的packing格式的话,您就会发现,每个分子都通过N=C键与3个分子相互连接。作者为了研究这种分子组建的cluster,所以取出了两个分子组建而成的Cluster。如下图:

当然这张图中显示了5个空的N=C键。这是通向其他晶胞的桥梁。
因此说,作者的计算模型不但是正确的,而且是准确的。如果作者仅仅考虑一个分子,才是滑天下之大稽。当然,如果作者有心,可以计算三个分子的cluster、四个分子cluster……等。下图显示了一种可能的cluster生长方式(沿c轴生长,绿色标记)。

当然如果沿着a轴生长或者沿着b轴生长的cluster,就不是这样的,您自己生长一下看看吧。不过,无论哪种生长方式,作者的二联cluster总是它们的一种基本结构,是这三种生长方式得到的大型cluster的共有的基本结构,因此我很佩服作者的cluster抽取能力。

第二个问题:如果两个分子已经桥联了,那么自然会形成闭壳层体系。

[ Last edited by yjcmwgk on 2010-3-25 at 11:26 ]
3楼2010-03-25 10:11:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

faqianliu

银虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yjcmwgk(金币+10):先向您道个歉,然后鼓励一下。周围没有量化的环境,自己研究,确实很艰苦呢 2010-03-25 11:34
erylingjet(金币+5):欢迎常来! 2010-04-01 17:53
引用回帖:
Originally posted by yjcmwgk at 2010-03-25 10:11:56:
楼主的问题真是莫名其妙啊。很明显楼主并没有经过仔细的思考,就匆忙上来发帖求助了。

作者在原文中已经明确说,他计算的是cluster。如果你看看这个cif文件的packing格式的话,您就会发现,每个分子都通过N=C键 ...

呵,首先感谢小卒耐心解答!

其实刚才也是心里有些郁闷,还是那句话,如果懂的话谁愿意在这发帖子提问,不懂的话提问又没人解答,只能发发牢骚了,其实有些问题对你们来说很简单,可对我们来说,如果有人稍微指点一下,就能想通很多事情

接触高斯也有一段时间了,软件是网上搜索的,资料也是从网上搜索的,中间有很多人热心人帮助,但是再帮助,一个TD-DFT计算荧光,我就纠结了很长很长时间,到现在结果还不太明白,为什么呢?因为从写gif就是自学,看结果自学,什么都是自学,如果我们有一个工作团队,师兄师弟都可以帮忙,那肯定是没问题,有问题问问,周围有懂的也行,问问,也没有。。。

我们现在看资料,有不懂的首先也是去搜索,去找答案,可是网络不是万能的,不可能我们需要的资料网上就能一搜就找到的,所以就可能浪费很多时间,我们提的问题在你们学的好的人来说非常简单,甚至非常可笑,为什么会问这种问题呢?

还有gaussian 03手册,看过啊,有事先翻翻,但是没基础的话能看懂多少?

所以我们还是在慢慢摸索着前行吧,只能一点一点的积累知识,然后再等自己某天突然豁然开朗,噢,原来如此。。。

说这么多,心里话,还是非常感谢楼上的指点!
4楼2010-03-25 11:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见