几天前咱们量化版新开了个子板“文章故事”,在版主的帮助下,我用两片JPC写了一个帖子(http://muchong.com/bbs/viewthread.php?tid=3780634),拿下第一滴血(500+金币和4个QC强帖指数)奖励标准可见看区长的帖子(http://muchong.com/bbs/viewthread.php?tid=3780342)真的很丰厚哦,可惜各位高手无人跟进,这就违背了版主让我写那个帖子的“示范带头作用”了呀。在版主要求下,我只好再写一个帖子。希望大家积极跟进呢,奖励很大的哟!
这个子板的本意就是“传授经验,植树造林,接受批评,奖掖后进”。我们很希望能够办下去。但是,如果大家不肯跟进,这个栏目很快就会被取消的。毕竟小木虫是个大站,它有自己的规则。一个刚刚试行的新版块,需要在一定时间内拥有一定的人气,否则就会撤版。如果要维持这个版块,我们很需要大家写原创帖。
好吧,我再写一个帖子来提升下人气。我们真的很需要大家踊跃参与的。
用很简单的DFT计算,为实验化学家答疑解惑
讲述人:密度泛函·小卒
北京科技大学 化学与生物工程学院
讲述论文: Phys Chem Chem Phys, 2011, DOI:10.1039/c1cp20056a
作者位次: 第一作者
计算方法: DFT、Fukui Index、概念密度泛函的分子软硬度理论
计算基组: 优化6-311+G(d,p)/SDD,性质6-311++G(2df,2pd)/SDD
所用软件: Gaussian03D01, NBO3
影响因子: 3.453
全文链接: http://pubs.rsc.org/en/Content/ArticleLanding/2011/CP/c1cp20056a
首先说一下我所在的课题组环境。我们课题组基本上是实验课题组,大老板并不太懂理论上得东西。之所以养着我们这凤毛麟角的几个专门做理论的人,是为了给实验提供支撑的。所以你们会发现我们的理论文章基本上总是围绕实验展开。如果想了解我们的课题组,可以去看看这些我们课题组的一些实验化学的论文(JACS,2009,17044; JACS,2008,11623; JACS,2003,11257; JACS,2010,16460; JACS,2005,15700; JACS,2009,30; 综述Acc.Chem.Res.2009,78。先列举这么多吧,)。不过,纯理论文章的档次就低多了,最好的一篇也仅仅是我大师兄做的(ChemEurJ,2007,9503),另加十几篇我们这两个师弟一个师妹,以及大师兄自己做的JPC\CPC\PCCP这个档次的文章而已。我们与实验结合的理论计算则大量充斥在我们组里很多人的实验文章中,太多了就不一一列举了。
好了,广告到此为止。下面说这篇文章的起源。
有一天大导师跟我说去让我解释一下这个问题。如图1,一个带负电的有机配体,也就是图中最右侧的双层配体,要和一个带正电的“金属离子-单层配体”复合物相结合,形成一个三层的化合物。看起来是没有什么问题的,但是实验上从未做成过(Jiang, et al: (1) Chem.Ber. 1996, 129, 933; (2) Inorg. Chem. 2004, 43, 4740; (3) Dalton Trans. 2010, 39, 1321)。反而是如图1最左侧的那个中性态的双层配体,要和一个带正电的“金属离子-单层配体”复合物相结合,形成一个三层的化合物,却是屡屡实现了的反应。这一点很让实验化学家疑惑。按理说,有正负电荷相互吸引这个好处,应该更容易反应才对啊。于是大导师让我去解释这个问题。
图1
很快,二导师跟我说去让我解释一下另一个问题。在晶体中,我们发现这种双层配合物的扭转角有很多种情况,从0度到45度都有可能,以及各种中间状态都有发现。(关于这一点,可以到PCCP,2011,13277去看它的参考文献第22-28号)那么,这两层之间互相旋转是任意的吗?它们这种旋转有没有规律可循?于是二导师让我去解释这个问题。
反正都是基于同一物质的计算,我为什么要算两次?我决定一次搞定,作一篇文章解释两个问题不就结了吗!
对于大导师提出的第一个问题,我可没兴趣去做反应历程的计算。那太扯了,反应过程无比复杂,我大师兄猜了半年都没猜出来,我才不去做。于是决定偷懒,用概念密度泛函的Fukui Index和分子软硬度理论来解释。计算完毕后发现,图1中带负电的那个有机双层配体,虽然可以把带正电的“金属离子-单层配体”复合物吸引到自己身边来,但是由于双层配体本身的电子结构趋于稳定,所以也仅仅就是把人家吸引到身边来,不反应!而如图1最左侧的那个中性态的双层配体,电子结构非常不稳定,在溶液中碰上那个“金属离子-单层配体”复合物的概率是比较低,但是只要碰上就非得把人家反应掉才好。所以在真实实验中,动辄煮上两个小时,大家也就都反应好了,产率一般不低于30%。
图2
对于二导师提出的第二个问题,我可没兴趣去做势能面的扫描。那太扯了,我的机器算一年也未必能算完,我才不去做。我直接扫了三条势能线,如图2,27个点的部分优化,就结束啦!发现这个旋转势垒仅仅30多kj/mol,这在室温下就可以热驱动旋转了啊。在晶体中随便几个取代基带来的空间位阻就能让它转过去。然后找了一下它们的前线轨道,解释了稳定性的来源,如图3,解释结束!
图3
于是文章就结束了。是不是很简单?事实上这里面我没有使用任何高深的理论或者高难度的计算技术,事实上我在这篇文章中使用的计算技术比版上各位高手的计算技术都更加简单粗糙。我总觉得计算技术虽然是我们必备的,但不是人家实验化学家看重的。能为实验化学家解释问题的理论计算,才是实验化学家所看重的。尽管简单,尽管粗糙,却是有用的。
前几天我写的那个《基于TDDFT设计有机分子,以期实现太阳能的全谱吸收》里面所讲述的两篇JPCA,里面确实涉及到一些计算技术的讨论。但是这一片PCCP,你们也看到了,几乎没有任何技术可言。当时coolrainbow和sobereva两位牛人都感觉这个文章无理论价值。我完全承认,这个文章确实没有任何理论价值,两位牛人的眼光是正确的!但是,就让我自我安慰一下吧,这个文章确实还是有一点儿指导实验的价值了。其实这也是我的无奈呢。毕竟我是在一个实验课题组里,大老板养着我们这凤毛麟角的几个专门做理论的人的唯一理由就是为了给实验提供支撑的。所以你们会发现我们的理论文章基本上总是围绕实验展开。如果说我其他论文还有“身在曹营心在汉”的感觉的话,那么这个PCCP文章纯粹是当了“理论化学的汉奸”,彻底为实验化学家提供支撑了。
最后,我再多说两句:几天前咱们量化版新开了个子板“文章故事”,在版主的帮助下,我用两片JPC写了一个帖子(http://muchong.com/bbs/viewthread.php?tid=3780634),拿下第一滴血(500+金币和4个QC强帖指数)可惜各位高手无人跟进,这就违背了版主让我写那个帖子的“示范带头作用”了呀。在版主要求下,我只好再写一个帖子。希望大家积极跟进呢,奖励很大的哟!
这个子板的本意就是“传授经验,植树造林,接受批评,奖掖后进”。我们很希望能够办下去。但是,如果大家不肯跟进,这个栏目很快就会被取消的。毕竟小木虫是个大站,它有自己的规则。一个刚刚试行的新版块,需要在一定时间内拥有一定的人气,否则就会撤版。如果要维持这个版块,我们很需要大家写原创帖。
好吧,就让我再次抛砖引玉,用这个新帖子来提升下人气。我们真的很需要大家踊跃参与的。
[ Last edited by yjcmwgk on 2011-11-6 at 20:37 ] |