| 查看: 19780 | 回复: 56 | |||||||||||||||||
| 【奖励】 本帖被评价38次，作者sobereva增加金币 34.2 个 | |||||||||||||||||
| 本帖产生 1 个 QC强帖 ，点击这里进行查看 | |||||||||||||||||
[资源]
【sobereva个人文集】电子定域性的图形分析 已有2人参与
|
|||||||||||||||||
|
注：Multiwfn 1.4尚未正式发布，但可以从http://multiwfn.codeplex.com下载临时版本，文中涉及的功能皆已支持。 电子定域性的图形分析 文/Sobereva First release:2010-Jun-13 摘要：本文回顾了电子定域性分析的基本问题，介绍了拉普拉斯值函数、ELF函数、LOL函数的基本概念，利用一些实例介绍了这三种函数能够说明的问题以及如何对分子进行分析。同时穿插介绍了函数绘图的方法和经验，最后给出利用Multiwfn生成函数图像的简单例子，便于读者将图形化分析方法投入实际应用。希望本文能对推广分子图形化分析起到一定作用。 1. 电子定域性分析杂谈 波函数分析是计算量子化学方法中重要却常被忽视的一个领域，现在更多的人着眼于能量、几何结构、光谱等性质的计算。然而，为了获得精确的属性值而同时不得不“附带”得到的精确波函数却往往被忽视掉。波函数无异于一个黑箱，蕴含着分子一切信息，计算出来的分子性质是已经从中被提取的部分，然而还有无穷未被挖掘的信息，若不利用明显是巨大的浪费。波函数分析的目的是从波函数及衍生信息中提取易于被人所理解的概念，是有助于从更深层次认识分子本性的一大类方法，它既包含理论，又包含“技术”，也就是如何具体下手去解释计算结果。波函数分析的范畴并没有具体的定义，我认为它应包括拓扑分析（如AIM）、布居分析、键级分析、属性分解分析（如能量分解、电荷转移分解，轨道成分分解）、轨道定域化分析、图形化轨道相互作用分析（如前线轨道理论）、概念密度泛函理论等内容。然而绝大部分量化教材里缺乏单独的章节介绍这些方法，是不应该的，很多初学者遇到相应问题也毫无思路，例如总有人问怎么求每个原子轨道对某个轨道的贡献（多数人都说用系数的平方，以讹传讹）。 化学键、孤对电子、sigma/pi/delta键、轨道杂化、lewis结构、VSEPR这些经典理论诞生于大半个世纪前，属于“概念化学”范畴。这些模型化的表述从当今量子化学角度看似乎缺乏根据，是一些很粗俗的观点，但实际上它们并非是完全人为臆断的，它们从量子化学的角度是可以找到充分依据的，定域化分析正是力图通过波函数或其衍生信息还原这些信息。当然我们也不能固步自封，凭借量子化学的威力，在不失传统概念简洁、形象的优点下，可以得到更多的信息，比如精确的定量数据，还可以将应用范围扩展到难以通过经验分析的体系，比如团簇。目前已有不少轨道定域化方法力图实现这样的目的，其中很多是根据一定规则令原本离域的正则HF/DFT占据轨道作酉变换获得定域化轨道(LMO)，常用的比如Boys、Edmiston–Ruedenberg、Pipek–Mezey定域化方法，NBO分析尽管原理不同，但从目的上也可以归为这类。然而如何变换成LMO的方法不是唯一的，由此产生一些矛盾。比如双键，有的方法能得到sigma和pi键形状的LMO，然而有的方法却得到两个香蕉键。再比如苯环，lewis共振结构至少能写出两种，究竟LMO该表示哪种结构？可能分子取向、构型稍稍变化一点，根据程序的判断LMO就从描述一种共振式突然变成描述另一种共振式了，这样的不连续描述若用在某些场合会造成问题，而且从概念上也莫名其妙。 还有一类方法是使用实空间坐标为变量的函数来描述电子定域性特征，这不需要明确地计算定域化轨道，因而避免了上述问题。这些函数将复杂的3N维波函数所携带的关键信息约化、抽象成三维函数，从这些函数在分子附近空间中各点的数值可以分析分子各个位置定域性特征的强弱，使化学键、孤对电子、原子壳层这样电子定域性强的区域充分凸现出来。 PS: 回想基态电子密度包含了分子一切信息，利用各种形式的实空间函数，不仅可以做定域性分析，还有可能将构成分子各种属性的因素皆以图形化方式清晰、生动地展现出来，并由此总结出普遍规律用于分析未知分子的性质，这将有重要的意义。探索、发展这样的函数，可能是潜在的热点。然而现在多数人还是热衷于计算各种分子各种属性的数值，追求精确到小数点后的位数，这实属“暴力”的研究，理论研究应当比这更美。 另外，轨道定域化方法依赖于对基组的选取，假设用非定域的基函数，例如平面波、多项式函数，乃至离散的格点描述，由于无法直接变换成定域轨道，轨道定域化方法就不适用了，尽管也有人用正交投影的方法投影到原子中心的基组来解决这个问题，却不可避免地造成波函数质量的降低。基函数无非只是数学工具，依赖于它的轨道定域化方法从根本上缺乏物理意义。而使用函数的描述形式就避免了这一问题，只依赖于空间各点属性的值，而不管这是以何种数学形式表达的，无论对于适合定域基函数描述的分子体系还是适合离域基函数描述的周期性体系，这类方法都同样适用。 一般通过作图来研究这些函数数值的分布情况，表现切面上的情况主要通过填色图（也称赝色图）、等值线图、梯度线图、地形图(Relief map)。曾经作这些图操作十分复杂，缺乏简单、通用、免费、小巧、快速的软件，Multiwfn彻底解决了这一问题，使函数图形化分析变得十分容易，在1.4版中加入了更多的函数，包括下文的ELF和LOL函数，下文的曲线、切面图都是通过由Gaussian计算波函数（大部分为HF/6-31G*）并载入Multiwfn后直接获得的。函数的三维空间分布一般通过等值面图、颜色映射等值面图（一般将属性映射到电子密度等值面）研究，常规做法是先计算格点文件，它包含了离散化的空间上各点函数的数值，化学领域最常用的格点文件类型是高斯的cube文件，可以由Multiwfn等程序生成。然后载入可视化程序，GaussView、Chemcraft、VMD等都支持Cube文件，在Multiwfn中也可以在生成格点文件后直接观看等值面图。 下文主要谈谈这类函数中拉普拉斯值函数、ELF和LOL函数。 2. Laplacian函数 拉普拉斯值函数L(r)的值被定义为r处电子密度Hessian矩阵的迹，也就是拉普拉斯算符作用到电子密度的结果，亦即电子密度在三个正交方向曲率和。一般这三个正交方向就是指X,Y,Z方向，由于可以对Hessian做酉变换来旋转坐标系而保持迹不变，所以三个正交方向是任意的。在AIM理论的键临界点位置，可以将坐标系变换为一个平行于和两个垂直于键轴的方向而使其意义更明显，其各自的曲率在AIM分析中有重要意义，但本文不考虑拉普拉斯值的具体成分而只考虑总数值。拉普拉斯值函数在某处为正，说明此处电子密度以发散为主；若为负，说明以聚集为主。非极性共价键、孤对电子处由于电子呈聚集状态，因而它们都处于负值区域；极性键也有电子凝聚区域，相对于非极性键发生了偏移、变形。对于这样由价层电子聚集产生的效应称为VSCC(Valence Shell Charge Concentration)。对于闭壳层键（离子、氢键、VDW键）则不出现电子凝聚的成键区域。通过观察拉普拉斯值可以了解分子的特征，这比直接观察电子密度图清晰得多，电子密度从整体看都是以原子为中心向四周呈指数形衰减，成键、孤对电子引起的电子密度聚集仅能微弱影响这个总趋势，从电子密度图上很难分辨其效应。而求电子密度二阶导数就能将这样细部信息鲜明地展现出来，如同信息的放大镜（有兴趣可以做一下比如8*exp(-7*x)+exp(-10*(x-0.5)^2)的函数图和它的二阶导数图来对比一下）。  上图的右上是苯分子的拉普拉斯等值面图，这是由Multiwfn生成此属性的格点文件后直接显示出来的，蓝色为负值，绿色为正值。可见分子的整体轮廓被等值面勾勒了出来，成键区域处拉普拉斯值为负。原子由于具有壳层结构，在每个壳层之间都有拉普拉斯值数值为正的区域，C原子处显示的绿色球层就是这个区域。H是比较特殊的原子，仅有单个质子，其电荷与电子的束缚作用并不足以构成具有独立特征的区域，无法形成壳层结构，甚至可看成质子化的孤对电子，所以与成键区域连为一体而没有截面，对于ELF/LOL函数，氢原子也表现这样特征。 图的右侧中央是1-3丁二烯的拉普拉斯等值面图（这三幅图等值面并不一样，没有完全的可比性），众所周知中间的C-C键比两端的C-C键要弱，从图上也能看出这一点，两边的C-C呈椭圆形而中间的C-C键类似于柱形，这是因为两边的C-C键中的pi键更强所致。左图是此分子平面的等值线图，能更清楚地反映细节，箭头所示的是键之间的拉普拉斯值函数极大点所在区域（此处等值线为-1.19）。中间的C-C sigma键之间有两个成键极大点，两个成键极大点中间存在鞍点。两边的C-C键距离更近，成键极大点之间也更为接近，对于乙炔，C-C键更短，中间只有一个极大点。这与后文ELF/LOL函数的情况有所不同，ELF/LOL函数描述的sigma键总是只有一个极大点。 右下方的图是Sn2H4，它是反式弯曲而不像乙烯那样是平面构型，其原因曾有人指出是由于CH2的三重态能量低于单重态，s和p未成对电子可以用于形成sigma-pi键；而SnH2单重态能量低于三重态，呈sp2杂化，其孤对电子会与另一个SnH2的空着的p轨道配位，因而两个SnH2之间靠两个配位键链接，导致弯曲构型。从图上可看出Sn确实有sp2杂化的孤对电子，且冲着另一个Sn的p空轨。 用拉普拉斯值函数分析定域性问题的缺点是数值范围不确定，需要反复调整作图的数值范围，远没有ELF、LOL函数有确定的[0,1]区间那么方便。拉普拉斯值在原子内层区域数值波动极大，而且与价层区域不在一个数量级，做填色图的时候原子内层很不好看，做地形图时还必须将其数值在某处截断。虽然大部分能够用ELF、LOL分析的体系用拉普拉斯值函数也能分析出相似的结论，但是对具有定域性特征的区域分辨能力比ELF和LOL偏弱，有时不很明显，甚至几乎无法分辨，有文章指出拉普拉斯函数没法展现出价层主量子数很大的原子的壳层结构。拉普拉斯值函数难以凸显的特征都能用下面的ELF、LOL函数很好地表现。 3. ELF(Electron Localization Function): ρ_cond(r1,r2)称为条件电子对儿概率，代表已经知道有一个电子出现在r1处（称参考点），在r2处出现另一个相同自旋电子的概率，是电子对儿密度ρ(r1,r2)除以ρ(r1)的值。可以将之变换为ρ_cond(r,s)，称为球形平均电子对儿概率，代表令r为参考点，以其为中心在半径为s的球面上出现相同自旋电子的概率。将ρ_cond(r,s)对s做Taylor级数展开并截取第一个非零项就定义了D(r)函数，代表在r附近出现另一个自旋相同电子的概率，它在r处的值越小被认为r处电子定域性越强。D(r)的值一定大于等于0，为0时则表明完全定域。D(r)是三维函数，比六维函数ρ_cond(r1,r2)或四维函数ρ_cond(r,s)更易于考察。 D(r)的具体表达式为D(r)=1/2*∑[i ]N(i)*|▽ψ_i(r)|^2 - |▽ρ(r)|^2/(8*ρ(r))，这里N(i)为第i条轨道轨道占据数，▽为梯度算符，ψ_i是第i条轨道波函数（后HF波函数时为自然轨道），ρ(r)是总电子密度函数。根据Savin等人对其的解释，其中1/2*∑[i ]N(i)*|▽ψ_i(r)|^2是当前体系动能密度函数（也称G(r)或梯度动能密度），|▽ρ(r)|^2/(8*ρ(r))是Weizsacker泛函，表现了没有Pauli互斥时，如同玻色子体系的动能密度，因此D(r)本质上体现的是由于Pauli互斥造成的超额动能密度，也称Pauli动能密度，体系超额动能密度越低的区域定域程度越高。用D(r)描述电子定域情况有不妥之处，因为其数值越小代表越定域，这和观念相反，而且它的数值没有上限，观察时不便。所以Becke等人将之变换后定义了ELF函数，ELF(r)=1/(1+(D(r)/D_0(r))^2)，这样就把数值投影到了[0,1]区间，且数越大表明定域程度越高。这样将D(r)变换为ELF(r)的方式本身是任意的，并没有确定的物理含义。式中D_0(r)=3/10*(3*pi^2)^(2/3)*ρ(r)^(5/3)，是Thomas-Fermi均匀电子气动能泛函，表现了电子完全离域化的状态，所以当ELF(r)=1/2时，即D(r)=D_0(r)，说明此处电子没有定域性。 ELF函数的定义略有缺陷，当r趋近于无穷时，D(r)会渐进地接近0，所以在远离分子的位置ELF会等于1，成了完全定域的情况，这显然不对。但由于D_0(r)具有ρ(r)^(5/3)的行为，往往消失得更快，所以(D(r)/D_0(r))^2在远处一般为无穷大，故ELF在远处一般还是会等于0。然而有时D(r)比D_0(r)更早地接近0，此时ELF在分子外部区域还是会等于1。为了避免这种情况，可以在D(r)上加入极微小量，比如10^-5，这样可确保(D(r)/D_0(r))^2在电子密度已经很小的分子外部区域变得很大（因为D(r)至少有10^-5），故一定能令ELF(r)等于0，这只对分子内部区域ELF函数产生微不足道的影响。 从分子的ELF函数分布上可以清晰地看到原子壳层结构、孤对电子、化学键等特征。除Multiwfn外，CheckDen(http://atenea.montes.upm.es, J.Mol.Graph.Model. 28 (2009) 102-112)、DGrid、(http://www.cpfs.mpg.de/~kohout/dgrid.html)、TopMoD(http://www.lct.jussieu.fr/pagesperso/silvi/silvi_english.html, Comp.Chem. 23 (1999) 597-604)也可以进行ELF函数计算，TopMoD还支持拓扑分析。 下图为Kr (cc-pVTZ)的从0到6波尔半径的ELF函数曲线，可见主量子数为1,2,3,4的四个壳层清晰可见，若只是考察总电子密度的变化则只能看到呈指数衰减的曲线，获得不了任何细节信息。如果以曲线各个极小点对应的球面此对空间进行划分，并分别积分其中的电子密度，就能得到每个壳层内的布居数，数值近似等于相同主量子数原子轨道上电子数的加和值，差异来自轨道隧穿效应，而在量化中近似等于轨道指数处在相应壳层范围内的基函数的布居数加和。  下图中间的图是吡嗪(C4H4N2)分子平面ELF函数的地形图+投影图，投影图部分单独拿出来就是填色图。可见C-N、C-C化学键、N的孤对电子（与C-H区域形状较像）都清楚地显示了出来。  这里顺便说一下填色图的色彩刻度的设置问题。填色图就是根据每个点数值大小，用色彩刻度条上对应数值的颜色显示，最常用的色彩刻度从下限到上限是由蓝变绿再变红，类似彩虹。色彩刻度的上下限的设定是任意的，这会影响到图上格点的颜色，设置标准的关键是要清晰地反映数据的差异。此图将色彩刻度设为了0~1，这一方面是因为ELF函数数值在0~1范围内，但最主要是因为在0~1范围内填色图能比较清晰地反应各个定域化区域。左下角的图是将色彩刻度改为0~1.3的图，原本在红色区域（数值接近1）的位置都变成了黄色，成键区域与绿色区域不易区分，给人的感觉不够明显。右上角是将刻度设为0~0.8，那些超过0.8的数值通常就被自动设成了白色，这样0.8~1.0部分的数值变化就没法分辨了。所以上下限的设置应当比较接近数据集数值主体部分的上下限，这样就能充分利用各种颜色将不同数值区域的点分开。从经验上看，对于ELF函数的填色图适合用0~1的范围，但对于某些分子应进行适当调整，对于LOL函数用0~0.75比较好。 下图是HCN分子在Chemcraft下做的ELF函数等值面图，格点为80*80*150，由内往外的等值面分别为0.97、0.87、0.5。可以看出N的孤对电子以及H附近定域程度最高；随后是一个sigma键和两个pi键，在图中合在一起呈现环形等值面，尽管这看似与一般观念上sigma键和pi键的形状很不像，但实际上是等价的。继续降低等值面数值后pi键等值面范围扩大，成为血小板形状。  类似这样结合不同透明度、不同颜色、不同显示方式，可以同时在一张图里不受影响地表现多层等值面，表现丰富的信息。在Chemcraft里每绘制完一层点一下keep this surface，则已绘制的层会一直保留，不受之后操作的影响，调整设置后绘制新的层，反复如此就可以实现这样的效果。保存图像不宜直接截图，应当用file-save image保存，或者用View-draw export picture再截图，这样就会使用抗锯齿技术使边缘和线条平滑，尤其是用网格线表现等值面的时候效果会好很多。 对ELF函数也可以进行类似AIM的拓扑分析，是对AIM分析的推广，由此可以比图形化观察更细致地定性或定量对体系进行分析以得到更多信息。ELF拓扑分析也有ELF函数的极大点（也称点吸引子）、极小点、鞍点、零通量面这样的概念。零通量面划分了盆，不仅包括AIM原子盆那样空间中的一团的形状，还有环形、球层等形状。此外还多了一些额外概念，环吸引子就是比如前面在HCN中看到的C-N三键的圆环，球吸引子就是例如Kr径向ELF函数曲线中每个极大点所对应的一层球面。f-localization domain是指ELF函数值为f的等值面包含的空间区域，域内各点数值都大于f，若只包含一个吸引子称之为不可约的，否则称为可约的，对于可约的定域性化域，当f逐渐加大并超过某一个值后就会分裂为多个包含更少吸引子的定域化域。对于HCN的例子，若将等值面数值进一步降低，三个彼此不相接触的不可约域就会扩张并连通成一个可约域。随f的加大，由最初包含整个分子的可约域最终分裂为各个不可约域的过程可以绘制成树状图。 这些概念对分析问题很有用处。按Bader的分类化学键可分为共享电子相互作用和闭壳层相互作用两大类，前者包括共价、金属键；后者包括离子键、氢键、二氢键、范德华键。Silvi和Savin在Nature 371 (1994) 683指出，对于前一类，在两个原子（氢原子除外）的键径总会出现至少一个用于成键的点吸引子或环吸引子（不参与成键时的孤对电子、内层电子的吸引子不算成键吸引子）；而对后一类不出现这样的情况。这样的规则与前述拉普拉斯值函数的情况很类似，本文各个例子也都验证了这一规则，这给成键的分类提供了明确依据，极大点的数值也可以评估键的强度。对共价键的盆空间内的电子密度进行积分能够用于了解有多少电子用于形成此键。比如Theor.Chem.Acc.108 (2002) 150~156文中就以两金属原子间是否存在ELF函数极大点判断它们是否有直接的共价键作用，并且根据积分Sc-Sc键的盆内电子数推断出除了3d以外，3s、3p也参与了成键。还有人系统地研究了原子电负性大小对某一f值的f-定域化域体积的影响，发现与VSEPR的假设比较吻合，并且指出键的盆布居数与相连的两原子电负性有相关性，键的极性越大，盆布居数越小于2，详见The Quantum Theory of Atoms in Molecules第六章的讨论。 [ Last edited by yjcmwgk on 2010-6-16 at 21:23 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | 第一性原理相关文档 | 高斯计算总结 | 量化软件学习 |
| 量子化学+实验研究 | dragonlxd的专辑 | 量化软件-Multiwfn | 金属负载氮化碳的催化剂 |
| mulwiwfn | niiiiiiiiiiiiiiiiiiiiiiiiice | 帖子 | 杂七杂八 |
| 学海无涯 | 作图技巧 | 学习 |
» 猜你喜欢
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士
已经有1人回复
-
固态锂金属电池迁移数计算求助
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有175人回复
-
荷兰奈梅亨大学赵文博与Galimberti课题组招收理论计算CSC博士 2000 欧元/月+房补
已经有0人回复
-
26储能博士申请自荐
已经有22人回复
-
cif求助下载
已经有1人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有2人回复
-
上海交通大学--宁波东方理工大学电池方向博士招生
已经有10人回复
-
海南大学全国重点实验室杨金霖老师招收储能电池方向博士生（2026年3月报名，9月入学）
已经有1人回复
-
海南师范大学2026年博士研究生招收 (在职想提升学历人员可报考) 申请考核制+普通招考
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用Multiwfn做拓扑分析以及计算孤对电子角度
已经有18人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
- 求助一个word2003的图形填充问题（只对特定区域填充）
已经有6人回复
- 判断自由基位置
已经有9人回复
- 【sobereva个人文集】制作动画分析电子结构特征
已经有29人回复
- 【原创】Multiwfn 2.0 正式发布
已经有35人回复
- chemdraw 的复制粘贴问题！！图片不全，还有空白区域！！
已经有14人回复
|
4. LOL(Localized orbital locator) Schmider和Becke在J.Mol.Struct.(Theo) 527 (2000) 51-61中定义了LOL函数γ(r)，γ(r)=t(r)/(1+t(r))，t(r)=D_0(r)/(1/2*∑[i ]N_i(r)*|▽ψ_i(r)|^2)。可见t(r)就是均匀电子气动能密度与当前实际体系动能密度的比值，某处动能密度相对于均匀电子气越小，γ越大，定域性越强。所以LOL和ELF引入D_0作为参考函数以获得相对值，是为了令体系各处动能密度的差异更显著地表现出来，否则结果很难考察。将t(r)变换为γ(r)则是为了使数值在[0,1]区间内。 LOL函数的本质也可以与轨道定域化方法进行比较来说明，t(r)具有轨道酉变换不变性，即是说用MO轨道或是定域化后的轨道其结果是一致的，因此这里假设轨道已经是定域化的。成键轨道波函数极大点的梯度为0，因而动能密度1/2*∑[i ]N(i)*|▽ψ_i(r)|^2很小（由于定域轨道交叠很小，可把加和号去掉而只保留此成键轨道），导致t(r)很大，γ(r)接近于1，由于轨道波函数是平滑变化的，故轨道波函数极大点附近的数值也很大。而在轨道交叠处，轨道波函数梯度较大，致使t(r)分母较大，t(r)的值较小，所以γ(r)也很小。这就是为什么成键区域LOL的值较大，而在定域性强的区域相交处值较小。Burdett等人也在J.Phys.Chem.A 102 (1998) 6366-6372中也从类似的轨道波函数的角度对ELF函数的意义进行了分析。 由于LOL与ELF本质上都是对动能密度函数不同形式的反映，所表现的信息是相仿佛的。实际上由于存在这样的关系：h_bar^2/(4*m)*▽^2*ρ(r)=2*G(r)+V(r)，V(r)为势能密度（总为负值），可见拉普拉斯值越小动能密度G(r)越小，所以拉普拉斯值函数的结论与LOL/ELF函数也很相似。用LOL能给出与ELF定性一致的结果，此文对ELF与LOL的讨论也是可以互相套用的。总体来说LOL给出的图形比ELF的图形更为明晰，更容易分析，而且函数形式更简单，也没有像ELF在远距离时数值可能成为1的情况。所以拉普拉斯值函数、ELF、LOL这三者中最推荐使用LOL进行分析。 下图是乙烯的LOL函数图，格点为120*120*80，用Chemcraft绘制，填色图颜色范围为-0.05~0.7，等值面为0.5。同时取了两个切面显示填色图，由于乙烯有pi键，所以在垂直于分子平面方向比平行于分子平面方向红色区域更大。结合蓝色小点表现的等值面更清楚地显示出pi键的空间形状，由于LOL值为0.5时正是体系动能密度等于相同密度下均匀电子气时的值，所以常用0.5的等值面来划分成键区域。Becke提出可以将0.5等值面以内的体积作为计算键级的依据。  下图是ClO分子的LOL填色图，黑线代表LOL函数为0.5的等值线。上面是Cl下面是O。Cl的内部第一壳层（中心圆点）、第二壳层（黄圈）都明显地展现出来，并且它们并未参与成键而变形，表现了内层电子的惰性。分子中间成键区域，以及原子旁边的孤对电子也都清晰地显现了出来。  下图从左到右是F2分子、CF4分子与LiF分子的LOL图，它们分别是非极性共价分子、极性共价分子、近乎完全由离子键连接的分子。从它们的LOL图上可以看出，F2中F原子发生了sp3杂化，孤对电子和共价键的特征充分显现了出来；CF4中成键区域比F2中显著，说明键更强，这在一定程度上能解释C-H与F-F平均键能差异之大的原因（485KJ/Mol和153KJ/Mol）。由于C和F电负性不同，C-F键从LOL图上也显示出了极化效应。当两原子电负性差更大，也就是LiF的情况，电子几乎全被F吸走而形成满壳层的F-，F-附近LOL函数值分布基本接近球对称，只因为Li+的正电荷极化效果略有变形而已，二者之间主要靠静电力束缚而没有共价相互作用，没有出现LOL函数极大点而只有鞍点。有兴趣的读者可进一步计算并分析BeF2、BF3、NF3、OF2的LOL函数图，研究随极性的变化成键区域的变化。要注意如同乙烯的情况，有的分子成键区域没有键轴对称性，不能仅考察分子平面的情况。  上面的例子凭经验都能猜到原子间是如何成键的，然而对于新颖、复杂的情况，就必须借助ELF或LOL分析了，例如6个Li组成的团簇，凭经验难以分析，而通过LOL函数图（下图左侧）这个问题变得显而易见。从图中可见每个角上的Li与它相邻的两个边上的Li形成了很强的三中心键。而图中右上角的拉普拉斯值函数地形图并没有给予我们任何有用的关于成键的信息，各个Li之间拉普拉斯值几乎为零，所以用拉普拉斯值函数分析成键问题不如LOL。这实际上与Li之间电子密度太小有关，必须将拉普拉斯图的Z轴上下限设得非常小才能看出端倪，缺乏固定的数值范围令拉普拉斯值分析很不方便。图中右下角是这个体系的密度差图，红色和蓝色分别代表孤立原子形成体系后电子增加和减少的区域。由于此体系密度差数值非常小，为了表示清楚，色彩刻度的Z轴上下限已设得非常小，可见Li背后的电子在形成强三中心键的区域富集，而三个处在边上的Li之间电子没有发生富集。对这个体系进一步进行Mayer键级分析，发现处在边上的Li与相邻的处在角上的Li之间的键级为0.386898，而处在边上的Li与另一个处在边上的Li之间的键级为0.165081，从LOL图上可见，前一种成键方式跨越红色区域程度很大，而后一种情况跨越程度很小，这在一定程度上解释了Mayer键级的结果。这个问题还可以进一步与多中心键级分析结果相比较。  ELF/LOL不仅可以用于讨论分子的静态特征，还能用于研究化学反应过程，通过观察反应路径中成键区域和孤对电子的出现、消失、范围和形状的变化，可以得到很多重要的结论，限于笔者时间和精力不在此举例，十分建议阅读LOL函数原文中的几个例子，其中包括SN2反应、水分子解离、HCN异构化、1,3丁二烯环化反应。 计算上述函数所用波函数质量会对结果有所影响。不要使用极小基，这样的基组变分自由度太小，缺乏柔性，价层区域的实际情况难以正确描述。至少使用双分裂价+极化函数级别的基组，6-31G*是比较适合的，对于阴离子、弱相互作用、激发态体系需要加上弥散函数，这样才能较好描述距原子核较远处的信息。然而再大的基组没有必要，上述函数的空间分布不像能量等属性对基组那么敏感，6-31G*一般已经足够，试图加大基组对结果影响微乎甚微。对于大体系如果用6-31G*基组仍然困难，可以使用混合基组，只在感兴趣的区域用大基组。从计算方法看，使用后HF波函数显然最好，但考虑计算时间，通常使用DFT波函数就可以。是否使用含电子相关作用的方法计算波函数对分析上述函数一般不会有定性影响，HF波函数也是完全可以的。但如果要求分析的精度比较高，比如几种分子的上述函数差异很小，却需要从这微小的差别中分析它们的不同点，或者要求准确的定量数据，比如求函数极大点的位置和数值或者计算f定域化域的体积，则应当使用更高质量的波函数。 上述函数也可以用于配合物体系，这总要涉及到赝势。利用上述函数分析化学上感兴趣的区域并不需要全电子基组，使用赝势会虽然使得原子内部区域函数值为0，也因此不会显示壳层结构，但价层的行为是完全正确的，感兴趣的区域一般只与这部分相关，所以不必担心赝势会造成问题。在介绍拉普拉斯值函数部分的Sn2H4就对Sn用了赝势，其孤对电子、与H的成键都很好地表现了出来。但有些分子比如Li6Ca2[Mn2N6]的Mn的内层电子也涉及成键[Angew.Chem. 110 (2000) 1667]，或者分子处于电离态而使内层电子直接暴露，必须明确考虑亚价层的影响，所以建议使用包含亚价层轨道的赝势基组。如果必须要用全电子基组的话，也可以在几何优化的过程中使用赝势来节省时间。 5. 使用Multiwfn绘图简单实例 ELF/LOL函数的意义显著，分析方法容易掌握，适用性广泛，熟练运用可以使文章增光添彩，值得大力推广，也很值得作为黑箱方法纳入普通化学教材。当然ELF/LOL分析若要普及开来必须有方便的程序，Multiwfn正符合要求。下面是使用Multiwfn获得ClF3分子在分子平面内LOL函数填色图以及等值面图的具体步骤，操作一遍就会发现这十分容易。此例只是Multiwfn最基本的应用，其它功能请参阅Readme文件的功能介绍和实例部分。 Multiwfn需要利用.wfn或.fch文件中的波函数信息进行计算。首先用GaussView生成ClF3的输入文件，并在Route section部分写上out=wfn，在分子坐标末尾空一行写上.wfn文件输出的路径。也就是例如下面这样，用高斯执行之： ===================== # HF/6-31G* opt out=wfn Title Card Required 0 1 Cl 0.00000000 0.00000000 -0.32113636 F -1.57000000 0.00000000 -0.32113636 F 0.00000000 -0.00000000 1.24886364 F 1.57000000 0.00000000 -0.32113636 c:\ltwd\ClF3.wfn [空行] [空行] ===================== 这样就得到了优化后的ClF3的波函数文件。先做填色图，启动Multiwfn，按照以下顺序依次输入，//后面是注释。 c:\ltwd\ClF3.wfn //输入文件名 4 //选择绘制平面图 10 //选择LOL函数 1 //选择填色图 200,200 //两个方向格点数，越大越精细 3 //选择YZ平面 0 //设定YZ平面是X=0的YZ平面，也就是分子平面 立刻图像就显示出来了，如下图左图所示。从图中可见，在赤道上的F与Cl成键效应比轴向的F与Cl成键更强。如果在高斯输入文件的Route section部分写上iop(6/80=1)来进行Mayer键级分析，会看到轴向的F与Cl的键级为0.621760，而赤道的F与Cl的键级为0.848993，与LOL函数的结论十分吻合。 在图像上点击右键来关闭图像，之后输入0将此图像保存到当前目录下的以DISLIN为前缀的.png图像文件中。 做等值面图需要先用Multiwfn计算格点文件，启动Multiwfn然后依次输入： c:\ltwd\ClF3.wfn 5 //计算格点文件 10 //选择LOL函数 80,120,120 //XYZ方向格点数，格点越多等值面图越精细。因为X方向垂直于分子平面，没有丰富内容，所以此方向格点数可减少些。 n //用Multiwfn默认的格点设置，不手动输入格点设定。默认设置一般是很合适的。 现在格点开始计算。计算速度取决于总格点数和高斯函数数目，数目越多计算越慢。Multiwfn计算ELF/LOL函数格点文件是很快的，比CheckDen快10倍。计算完毕后格点文件以LOL.cub为文件名输出到了当前目录下，然后输入n，代表不预览等值面。 现在将LOL.cub用Chemcraft打开，设定Contour value为0.55，点Show isosurface，等值面即显示出来，如图中右边所示。从立体的等值面图能看到Cl的两个孤对电子，而F的三个孤对电子组成了环形等值面。  [ Last edited by sobereva on 2010-6-26 at 21:48 ] |
2楼2010-06-13 06:26:38
3楼2010-06-13 07:08:57
4楼2010-06-13 08:21:25
好贴要顶。 5楼2010-06-13 08:34:25
6楼2010-06-13 08:35:52
7楼2010-06-13 08:54:44
|
一些问题， 1. 在非临界点，电子密度的一阶导数并不为零，拉普拉斯量对应的是二阶导数的变化（再求迹），而文中的 “极性键则有一部分落在正值区域，一部分落在负值区域，或者说电子凝聚区域发生了偏移。” 涉及了两个部分，至少有一个在非临界点，那么在这些区域为什么可以直接考察拉普拉斯值跳过一阶导数？至少在文中没有找到一阶导数的考察 2. "h_bar^2/(4*m)*▽^2*ρ(r)=2*G(r)+V(r)，V(r)为势能密度（总为负值），可见拉普拉斯值越小动能密度G(r)越小，" 更小的拉普拉斯值也可能是更低的是势能密度产生的，为什么可见动能密度越小？ 3. “从它们的LOL图上可以看出，F2中F原子发生了sp3杂化，” 据我所知，杂化轨道首先是对价键理论定义的，本文是否修改了杂化的定义？ 我没有做VB的计算，在初步的NBO分析上，基于CASSCF(10,10)/cc-pVDZ. 平衡位置(1.41193 angstrom), NBO提供的结果是 (Occupancy) Bond orbital/ Coefficients/ Hybrids --------------------------------------------------------------------------------- 1. (1.89815) BD ( 1) F 1 - F 2 ( 50.00%) 0.7071* F 1 s( 4.26%)p22.46( 95.61%)d 0.03( 0.14%) -0.0002 -0.2050 -0.0233 0.0000 0.0000 0.0000 0.0000 0.9775 -0.0219 0.0000 0.0000 0.0000 0.0000 -0.0367 ( 50.00%) 0.7071* F 2 s( 4.26%)p22.46( 95.61%)d 0.03( 0.14%) -0.0002 -0.2050 -0.0233 0.0000 0.0000 0.0000 0.0000 -0.9775 0.0219 0.0000 0.0000 0.0000 0.0000 -0.0367 .... 0.0000 0.0000 0.0000 0.0000 4. (1.99559) LP ( 1) F 1 s( 95.79%)p 0.04( 4.21%)d 0.00( 0.00%) 0.0000 0.9787 -0.0022 0.0000 0.0000 0.0000 0.0000 0.2051 -0.0022 0.0000 0.0000 0.0000 0.0000 -0.0023 5. (1.99152) LP ( 2) F 1 s( 0.00%)p 1.00( 99.99%)d 0.00( 0.01%) 0.0000 0.0000 0.0000 1.0000 0.0050 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0079 0.0000 0.0000 0.0000 6. (1.98716) LP ( 3) F 1 s( 0.00%)p 1.00( 99.99%)d 0.00( 0.01%) 0.0000 0.0000 0.0000 0.0000 0.0000 0.9999 0.0084 0.0000 0.0000 0.0000 0.0000 -0.0078 0.0000 0.0000 7. (1.99559) LP ( 1) F 2 s( 95.79%)p 0.04( 4.21%)d 0.00( 0.00%) 0.0000 0.9787 -0.0022 0.0000 0.0000 0.0000 0.0000 -0.2051 0.0022 0.0000 0.0000 0.0000 0.0000 -0.0023 8. (1.99152) LP ( 2) F 2 s( 0.00%)p 1.00( 99.99%)d 0.00( 0.01%) 0.0000 0.0000 0.0000 1.0000 0.0050 0.0000 0.0000 0.0000 0.0000 0.0000 0.0079 0.0000 0.0000 0.0000 9. (1.98716) LP ( 3) F 2 s( 0.00%)p 1.00( 99.99%)d 0.00( 0.01%) 0.0000 0.0000 0.0000 0.0000 0.0000 0.9999 0.0084 0.0000 0.0000 0.0000 0.0000 0.0078 0.0000 0.0000 .... 并没有找到sp3杂化的迹象 4. F2是一个多组态的分子，UHF可以给出解离的结果， http://arxiv.org/PS_cache/arxiv/pdf/0712/0712.1935v1.pdf 本文的计算水平是什么？RHF/6-31G*在优化平衡位置的能量(electronic energy) -198.6778 已经高于两个F原子能量的加和(ROHF -99.3618 ,UHF -99.3650)，以这种方式计算键能是负的；而F2在10 angstrom的RHF/6-31G* 能量为-198.2960，以这种方式键能为1002kJ mol^{-1} 。可见Hartree-Fock对F2并不是一个良好的近似。 LOL函数在这个分子对方法的敏感度是？ 5. 一个一般性和两个具体的问题。据我所知，并没有对于化学键唯一严格的量子力学定义，存在很多种键级方案。Laplacian对应的是局域virial系数，ELF/LOL对应局域过剩动能。如果我们基本上认为，化学键是一种能够降低体系能量的作用，而就文中的信息，Laplacian/ELF/LOL并没有直接对应体系的总能量，如何面对这种的区别？ 文中称“CF4中成键区域比F2中显著，说明键更强“ 为什么？ 如果考察从氢原子和质子结合到H_2^+的过程，从无穷远到平衡位置，动能先减小再增加， J Comput Chem. 2007 Nov 15;28(14):2389. Figure 15 http://www3.interscience.wiley.com/journal/113493175/abstract 如何在Laplacian/ELF/LOL理解这一点？ 6. 听过一种说法，量子化学的三个境界是，不计算，就知道；计算了，知道了；计算了，还不知道。轨道相互作用分析一个优势是，通过图像不进行复杂的计算也可以定性预言一些反应， 从Laplacian/ELF/LOL如何预言，而不仅仅是在计算的结果之后分析一些分子的行为？ 谢谢 [ Last edited by fichte on 2010-6-13 at 18:00 ] |
8楼2010-06-13 09:34:51
|
赞美sobereva大师！ 实验化学上的种种迹象都表明分子的性质具有某种定域化的趋势，但是分子电子波函数却总是含有大量的离域成分，这确实是个矛盾，事实上有人证明(jacs，116(1994)12022):分子内一个小邻域的电子密度rho可以决定整个分子的rho。按照这个定理可以说：一个化学意义上的“官能团”，如C=O，CH3，在两个不同的分子，甚至同一分子不同构象中都不可能有同样甚至接近的电子密度。这样来看，使得C=O、CH3在不同的分子中表现出的几乎相同的、定域的化学性质的因素可能不是电子密度（波函数模的平方），而是其他的因素。波函数的拓扑性质？甚至是其他未知的性质？ 事实上现在欧洲一些研究计算方法的也把注意力集中到了定域化上。Weiner等提出的LMP2、LCCSD、LCCSD(T)的方法对“extend”型的分子非常有效。但是这些方法终归是数学技巧而已，未能从物理本质上解释清楚分子性质定域性的由来。甚至讽刺的是，定域化方法(LMP2等)对传统意义上的“定域分子”如烷烃类的计算结果不如苯这样的“离域分子”。看来想要解释清楚分子的性质，必须要从问题的物理本质上入手，提出新概念，而不能再像过去三十年陶醉于越来越大的基组和组态波函数上来凑变分了，呵呵。 像sob总结的ELF,LOL，AIM，以及Faratos等人用分叉、奇异点方法研究分子的应该算是这类方法的尝试，不过本质上都是基于波函数的。波函数究竟能不能彻底或者说完备的反应分子的自然属性呢？也许只有等21世纪新一代的天才来回答了 |
9楼2010-06-13 10:01:54
10楼2010-06-13 10:33:09
11楼2010-06-13 10:41:51
12楼2010-06-13 10:50:05
14楼2010-06-13 13:01:08
|
村雨纯夏，你好： 1 可能没有图使我这句话造成了误会，我是从空间范围角度说的，临界点是在拉普拉斯负值区域的。我已把前半句话删了。 2 我记得Bader是在06年有篇文章讨论了ELF与拉普拉斯值函数的关系，具体可以参考那篇。 3 我这里没把杂化的概念定义得那么精确，这个概念已经被泛用化了。我是从常规的VSEPR对AXE3型分子预测结构和杂化的关系说的，从LOL等值面图可看出这种说法是比较易于接受的描述。由于最近时间紧，对于对同一问题各种表述的差异，等有时间时再做分析讨论。 4 由于只是定性上讨论，本文用RHF波函数，对F2这样经典体系是没什么问题的，提高波函数质量不影响定性结论。 5 这个角度从上述函数不容易直接体现出来。这是个比较大、值得深究的话题，等有时间时我将试图把各种对成键本质的观点放到一起综合讨论，写个小文，理清其联系。 6 上述函数主要在于解释问题，但也可以通过类比的方法从分析过的已知体系的性质的成因去解释未知分子的性质及其成因。若与轨道相互作用分析方法相互结合，应当能够互补从而成为预测分子性质更有利的武器，这值得探索。 coolrainbow说得很好。常有感觉，用大基组、更多的组态波函数暴力破解分子信息是虽然是主流，确是偏路。这种做法在后人看来无异于用在疯狂地摇手摇计算器，我们由于诞生得早而吃了大亏，但智慧却是公平的，今人并不吃亏，我们应当认识这一点来引导正确的研究方向。计算化学和理论化学的称呼也应该独立起来，前者利用技术解决现实问题，而后者侧重于探索化学问题的根本原理，形成机制，需要更发散的思维，更多的灵感，应当多鼓励这类话题的讨论。 |
15楼2010-06-13 17:44:04
|
谢谢您的回复，我想我不应该多次编辑原回复，在这个帖子继续吧。 2. 我到实验室再查查文献去 4. 重新用Gaussian 03 E01计算了一下F2分子，RHF/6-31G*在优化平衡位置1.3349 angstrom的能量(electronic energy) -198.6778已经高于两个F原子能量的加和(ROHF -99.3618 ,UHF -99.3650)，以这种方式计算键能是负的；而F2在10 angstrom的RHF/6-31G* 能量为-198.2960，以这种方式键能为1002kJ mol^{-1} 。可见RHF对F2并不是一个良好的近似。参见J. Chem. Phys. 86, 887-907 (1987). 在VB的计算上F2也是一个比较特殊的体系，一些作者将它分类为离子键共价键之外的charge-shift bonding http://www.nature.com/nchem/journal/v1/n6/abs/nchem.327.html 关于”分析“和”计算“，有很多种处理方法。“分析”可以基于第一性原理计算的结果；在一些情况下，如SAPT所做到的，也可以把对化学的"理解”集成到计算过程中去，静电，诱导，色散及其交换修正直接通过微扰展开得到不需任何后续处理。Morokuma-Ziegler能量分析可以看作是一种对结果分解的方法，也可以看作把两个fragments的波函数作一个symmetric orthogonalization构造了整个体系的inital guess，作为一种改进的初始猜测，我的经验是在ADF确实减少一些SCF的迭代次数。当然这种把能量分解视为计算过程的角度更多的停留在概念阶段，实际计算并不比默认的initial guess节省时间 [ Last edited by fichte on 2010-6-14 at 00:38 ] |
16楼2010-06-13 18:17:31
17楼2010-06-13 19:48:42
19楼2010-06-14 00:12:21
20楼2010-06-14 00:34:07
|
以前的确忽视了F2的特殊性，一般把它看成典型共价键，看来说法得改改。刚发现有篇文章Chem. Eur. J. 2005, 11, 6358 – 6371正是通过ELF论证了CS键的存在，等有时间读完后再把文中相应讨论进行修改。 用RHF/UHF分别在6-31G*扫了一下F2势能面，蓝线是两个UHF的F原子能量和，RHF能量的情况如你所说。拿RHF研究解离问题众所周知不妥，换成UHF后曲线比较有意思，如预期地在偏离平衡距离一定范围后出现对称性破缺使能量降低，一直比自由F原子能量高，也如预期在很远处给出两个自由F原子的能量。在6.3埃处能量突然落到蓝线上成为自由原子状态，这个不连续点比较有趣。  |
21楼2010-06-14 00:57:48
|
我用G03 E01基本上成功地重现了您的结果。一般而言，通过guess=mix并不足以保证计算的结果是一个局域极小点，需要通过stable计算electronic Hessian进一步确认。从r=re=1.3449 - 6 angstrom, UHF + guess=mix得到的electronic Hessian的最低本征值仍然是负的，stable=opt 可以得到更低的能量。而在7 angstrom, UHF + guess=mix得到了本征值为正的electronic Hessian,表现为一个突然地下落，这实际就是stable=opt一直得到的，不成键的解离曲线。而且对于F2，它的Hartree-Fock的instability point是在r_e内部的，参见R. Bauernschmitt and R. Ahlrichs, JCP 104, 9047 (1996). RHF的缺陷不仅仅在于解离，即使在r_e，得到的结果也高于两个F原子的能量 以下是计算数据 UHF, guess=mix stable=opt 1.2 -198.649467 -198.649467 1.3 -198.6755614 -198.6755614 1.3449 -198.677757 -198.677833 1.4 -198.675042 -198.681276 1.5 -198.659881 -198.692337 1.6 -198.637198 -198.703302 1.8 -198.589541 -198.718155 2 -198.548154 -198.725142 2.5 -198.471665 -198.729953 3 -198.427393 -198.730367 4 -198.392033 -198.729957 5 -198.377384 -198.729910 6 -198.368729 -198.729912 7 -198.729913 -198.729913 8 -198.729913 -198.729913 10 -198.729914 -198.729914 输入文件 %chk=F2.chk # UHF/6-31G* guess=mix stable=opt title 0 1 F 0.0 0.0 0.0 F 0.0 0.0 2.0 输出文件部分 .... SCF Done: E(UHF) = -198.548154231 A.U. after 12 cycles Convg = 0.4046D-08 -V/T = 2.0043 S**2 = 0.7423 Annihilation of the first spin contaminant: S**2 before annihilation 0.7423, after 0.0881 ... *********************************************************************** Stability analysis using *********************************************************************** Eigenvectors of the stability matrix: Excited state symmetry could not be determined. Eigenvector 1: ?Spin -?Sym Eigenvalue=-0.4190574 5A -> 10A -0.11731 9A -> 10A -0.69251 5B -> 10B -0.11731 9B -> 10B 0.69251 Excited state symmetry could not be determined. Eigenvector 2: ?Spin -?Sym Eigenvalue= 0.0000042 8A -> 10A 0.69712 8B -> 10B -0.69712 Excited state symmetry could not be determined. Eigenvector 3: ?Spin -?Sym Eigenvalue= 0.0343533 6A -> 10A 0.69637 6B -> 10B 0.69637 The wavefunction has an internal instability. .... Restarting incremental Fock formation. SCF Done: E(UHF) = -198.725142035 A.U. after 22 cycles Convg = 0.3263D-15 -V/T = 2.0019 S**2 = 0.9618 ... *********************************************************************** Stability analysis using *********************************************************************** Eigenvectors of the stability matrix: Excited state symmetry could not be determined. Eigenvector 1: ?Spin -?Sym Eigenvalue= 0.0035922 6A -> 10A 0.14291 9A -> 10A 0.68381 6B -> 10B 0.14279 9B -> 10B 0.68363 The wavefunction is stable under the perturbations considered. The wavefunction is already stable. The wavefunction is already stable. .... [ Last edited by fichte on 2010-6-14 at 16:31 ] |
22楼2010-06-14 02:08:39
23楼2010-06-14 04:39:28
24楼2010-06-14 07:45:08
25楼2010-06-14 15:33:32
26楼2010-06-14 18:33:42
27楼2010-06-14 20:50:05
|
对于一个组态函数的波函数，对它的inactive orbital之间的unitary transformation并不改变体系的能量。这实际是许多orbital localization的方法。一个localized的波函数并不一定包含大量离域的成分。 看了一下这篇文献，作者计算了一些fragment的电子密度，我没有发现关于这个定理的证明，也许是看得不细致，您能具体说明一下在哪里推出的吗？ 谢谢 可以通过local correlation发现是哪些电子的相关对应哪些能量上的贡献， 比如 J. Chem. Phys. 110, 7210 (1999)； Figure 2中包括的相关模式来自每一个monomer，这在一定程度上说明了定域性的来源。 个人不认为“官能团”在量子化学上是一个难以理解的概念。 毕竟可以写下甲基、乙基等等的electronic Hamiltonian，这些electronic Hamiltonian在不同的分子中，虽然v_{Ne}项（N代表原子核，e代表电子）因为原子核位置的不同而有所不同，大体上都是很相近的。Hohenberg-Kohn定理的出发点来自v_{Ne}决定了整个体系的Hamiltonian，既然官能团在不同环境里的v_{Ne}都是接近的，虽然带一些局域电荷可能造成了一些麻烦，在定性的概念上很容易理解官能团在不同环境的行为是接近的。 在技术上有很多种方法实现这一点，VB的波函数是定域化的；可以选择一些fragment做energy decomposition/SAPT；可以做一个多级展开，这样相互作用的Hamiltonian就被分解为monomer的多级矩算符的乘积，对应电多级矩和多级极化率，虽然这个微扰展开并不是收敛的，但在一些情况下还是可以给出积极的结果的。或者我们选择原子的波函数作为initial guess, 像搭分子一样完成能量最小化的过程， J. Chem. Theory Comput., 2009, 5 (4), pp 1027–1032 Multilevel minimization的成功在一定程度上反映了分子中原子的加合代表了整体的性质。觉得移植到官能团上没有原则性的困难。 [ Last edited by fichte on 2010-6-15 at 16:37 ] |
28楼2010-06-14 23:00:32
29楼2010-08-13 18:32:05
30楼2010-08-13 22:26:45
★★★★★ 五星级,优秀推荐
|
大家都太厉害了，努力学习啊！ 本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2139884&fpage=2&view=&highlight=&page=3 |
31楼2010-08-18 09:07:23
32楼2010-09-07 18:57:53
33楼2010-09-28 20:09:58
34楼2010-11-03 10:14:01
35楼2010-11-05 10:12:03
36楼2011-04-13 21:47:56
|
Excited State 1: Triplet-A 3.4898 eV 355.27 nm f=0.0000 44A -> 46A 0.52398 45A -> 47A 0.39231 44B -> 46B -0.52398 45B -> 47B -0.39231 This state for optimization and/or second-order correction. Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Triplet-A 5.0150 eV 247.22 nm f=0.0000 44A -> 47A 0.30979 45A -> 46A 0.61551 44B -> 47B -0.30979 45B -> 46B -0.61551 Excited State 3: Triplet-A 5.2509 eV 236.12 nm f=0.0000 44A -> 46A -0.41155 45A -> 47A 0.55200 44B -> 46B 0.41155 45B -> 47B -0.55200 |
37楼2011-04-13 21:50:10
38楼2011-05-10 16:20:35
39楼2011-05-25 09:58:55
★★★★★ 五星级,优秀推荐
|
现在想使用multiwfn做图，遇到一些困难，向你请教： 我的目的是，比较两个分子中某一部分的差异，一个有杂原子与双键的p-pi共轭作用，另一个没有。这两个分子差异比较大，因此做电子密度差图好像不合适。我现在应用你编写的multiwfn做了lol，elf图，可以看出两分子在电子定域方面的区别，但是有其他电子（或者说其他轨道）的重叠，因此还不是很明显。 我的问题是：multiwfn中有指定某些轨道或者排除某些轨道的功能吧？（根据手册是有，或者我理解有误）但是我不知道这些轨道具体的信息，比如是哪个原子的，是内层占据轨道，还是外层的孤对电子轨道，所以无法准确选定。如何获得这些信息？我的分子共有41个原子，88个MO，比较复杂。谢谢！ |
40楼2011-07-26 14:31:44
41楼2011-08-05 12:04:17
42楼2011-08-06 17:40:21
43楼2011-08-20 10:32:42
45楼2011-11-28 09:45:44
46楼2011-12-13 20:03:57
liuying287
禁言 (正式写手)
|
本帖内容被屏蔽 |
47楼2012-03-11 16:42:58
49楼2012-09-29 21:20:01
50楼2012-11-20 22:49:07
简单回复
heyo_12313楼
2010-06-13 10:53
回复
qzhaosdu18楼
2010-06-13 19:52
回复
[ Last edited by qzhaosdu on 2010-6-13 at 19:54 ]
liazhi200844楼
2011-10-02 08:55
回复
五星好评 谢谢分享！
tider48楼
2012-09-21 03:29
回复
五星好评 赞
