请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 1508|回复: 0

[Multiwfn资源与经验] 使用Multiwfn绘制态密度(DOS)图考察电子结构

[复制链接]

2万

帖子

25

威望

2万

eV
积分
45931

管理员

公社社长

发表于 2019-5-14 22:50:53 | 显示全部楼层 |阅读模式
使用Multiwfn绘制态密度(DOS)图考察电子结构

文/Sobereva@北京科音  2019-May-14

0 前言

基于主流量子化学程序产生的波函数信息,使用Multiwfn可以对孤立体系非常简单、快速、灵活地产生质量非常理想的态密度图,对于分析化学体系的电子结构特征很有帮助。本文将介绍各类态密度图(Total DOS、Partial DOS、Overlap population DOS、Local DOS)的原理,以及它们在Multiwfn程序中的绘制方法。DOS虽然是第一性原理计算范畴的重要概念,但本文对DOS的绘制只涉及量子化学计算范畴,Multiwfn也不支持基于第一性原理计算产生的数据绘制DOS。

Multiwfn从很早期的版本开始就支持DOS的绘制,后来又不断进行了改进和完善,如今已非常强大,显著强于其它任何能绘制DOS的程序。Multiwfn的DOS绘制功能已经被大量发表的文章所使用,比如J. Comput. Chem., 38, 1574 (2017)、RSC Adv., 5, 78192 (2015)、Nature Communication, 8, 14551 (2017)、Phys. Chem. Chem. Phys., 19, 23373 (2017)、J. Phys. Chem. C, 119, 8349 (2015)、J. Phys. Chem. A, 121, 4009 (2017)、Eur. J. Chem., 24, 17046 (2018)等等等等,这些文章也可以当成例子。

本文的内容非常深入细致,如果你很急着作出TDOS/PDOS/OPDOS图,可以先直接参考本文第3节的例子,之后若想全面了解Multiwfn的绘制DOS的功能,再看本文其它部分。本文内容对应2019-May-14更新的Multiwfn 3.6(dev)版以及之后发布的版本,绝对不要用更老的版本!Multiwfn可在其官网http://sobereva.com/multiwfn免费下载,相关基础知识看《Multiwfn入门tips》(http://sobereva.com/167)和《Multiwfn FAQ》(http://sobereva.com/452)。


1 态密度的物理意义与种类

态密度(Density-of-states,DOS)是一个以能量为变量的函数。DOS(E)代表在能量为E的位置,单位能量区间内态的数目。“态”具体指什么,看这个词具体用在什么地方。

一般说DOS这个词较多的场合是第一性原理领域,此时的态通常指的是解KS-DFT方程得到的单电子轨道。第一性原理研究的主要是周期性体系,有k点的概念,不同k点处轨道能级也不同。第一性原理计算可以得到物理意义严格的DOS,即通过对k空间进行积分得到。

量子化学主要研究的是孤立体系(分子、团簇等),当使用DOS考察体系电子结构时,态指的是分子轨道(可以由HF、半经验、KS-DFT等各种基于单电子近似的理论方法计算产生)。孤立体系计算的时候没有k点的概念,算出来的轨道能级都是离散分布的,因此DOS这个词对于孤立体系来说本身没什么意义,画出图来的话,就是下面这样,对应的公式也一起给出了(为了避免和后文说的其它类型的DOS混淆,这里用了Total DOS (TDOS)这个词)。
1.png
上图中虚线位置是HOMO的位置,公式里ε是轨道能量,δ是狄拉克δ函数(不懂的话Google)。从上图可见,孤立体系的TDOS图本身只不过相当于在分子轨道能级位置画一个竖线而已,这种图一点意思也没有,没法传递出比轨道能级分布本身更多的信息。

如果通过一些函数将上图中这些离散的竖线依次进行人为地展宽然后再累加,从图中将可以一目了然地看到不同能量区域内轨道分布的密集程度,如下图所示。注意这种展宽纯粹是人为的,仅仅是为了便于直观分析而做的,这样得到的DOS曲线并不是像第一性原理计算那样是以严格方式得到的。
2.png

不同的展宽函数把离散的竖线展宽出的峰的形状不同。绘制DOS时展宽函数一般用Gaussian函数,在Multiwfn里也允许用Lorentzian函数,或者二者混合产生的Pseudo-Voigt函数。函数表达式这里就不列出了,读者可以在Multiwfn手册里的介绍DOS绘制功能的3.12节中看到。不管用的是哪个展宽函数,都会牵扯到半高全宽(FWHM)这个参数,即竖线展宽出来的峰在高度一半的位置峰的宽度。显然FWHM设得越大峰越宽,曲线看起来就越平滑、越连续,而FWHM越小则对峰的分辨率越高,绘制上图用的FWHM为0.05 a.u.。FWHM这个值没有什么物理意义,设多大是具有主观任意性的,一般的设定标准是让得到的图适合、便于说明自己想说明的问题。下图分别是FWHM改为0.1 a.u.和0.02 a.u.时候得到的图。
3.png

在《使用Multiwfn绘制光电子谱》(http://sobereva.com/478)一文中,得到光电子谱(PES)的方式在本质上和上述绘制TDOS其实如出一辙。

还可以定义Partial DOS (PDOS),可以叫做“局部DOS”或“分数DOS”,它体现的是特定的片段对TDOS贡献的曲线。若恰当定义片段,通过PDOS图就可以较好地把握不同能量区间的轨道的本质、主体构成。显然,如果定义的所有片段的并集等于整个体系的话,那么每个片段的PDOS曲线累加在一起正好就是TDOS。对于某个片段A,其PDOS定义如下。其中F是展宽函数,Θ是轨道成份。
4.png

计算轨道成份的方法很多,详见《谈谈轨道成份的计算方法》(http://sobereva.com/131)和《分子轨道成分的计算》(http://sioc-journal.cn/Jwk_hxxb/CN/abstract/abstract340458.shtml)里的介绍。在绘制DOS的功能中,Multiwfn支持Mulliken、SCPA、Hirshfeld和Becke方法计算轨道成份。

重叠布居DOS(Overlap population DOS,OPDOS)对于考察片段间的相互作用比较有用。笔者在《Multiwfn支持的分析化学键的方法一览》(http://sobereva.com/471)里专门介绍过什么叫Mulliken键级,这和重叠布居(Overlap population)是一码事,可以体现原子间是成键作用还是反键作用,没看过此文的话一定要仔细看一下。文中还提到,重叠布居可以精确分解为各个轨道的贡献,将这个思想和DOS图搞到一起,就是OPDOS的定义了。下式是片段A与B之间的OPDOS的表达式
5.png
式中的那个Χ项,也等价于假设i轨道是占据的,然后这个轨道对片段A和B之间的重叠布居的贡献量除以轨道占据数。在某个能量范围,若OPDOS数值明显为正,就说明这个能量范围里的轨道如果被电子占据了的话,会对两个片段间的结合起到积极贡献,换句话说,这些轨道对这两个片段而言是成键轨道;若数值为负,说明这些轨道如果被电子占据了的话,将对片段间的结合产生不利贡献,相当于起到反键轨道作用。

虽然OPDOS确实能展现片段间相互作用的情况,但是由于Mulliken键级本身对键的强度反映不好,所以也别太把OPDOS当回事。而且,有弥散函数的情况下,OPDOS结果会很离谱,可能数据完全没有意义。所以,若计算波函数时有弥散函数则不要讨论OPDOS,若必须讨论OPDOS且之前的计算用了弥散函数,应把弥散函数去掉算个单点重新得到波函数后再做分析。Multiwfn是分析成键问题的百宝箱,里面有大把方法可以很好地讨论化学键特征,所以千万别拘泥于拿OPDOS说事,强烈建议仔细看《Multiwfn支持的分析化学键的方法一览》(http://sobereva.com/471)。另外,之前笔者还见到有人试图拿两个片段各自的PDOS曲线讨论它们的成键,以为在某些能量区间它们同时较大就说明二者间有成键作用,这明显毫无意义,且不说它们是相位相同方式重叠还是相位相反重叠,首先这种情况都没法体现这两个片段间存在轨道重叠,因为PDOS图根本都体现不出片段的轨道分布在三维空间中的什么区域。

还有一种DOS叫Local DOS,这个留到本文第4节专门介绍。

我们来看一个典型的DOS例子,是Multiwfn的原文J. Comput. Chem., 33, 580 (2012)里的图,将二茂铁的TDOS、PDOS、OPDOS绘制在了一起。为了便于说清楚轨道特征,一部分轨道的等值面图是用Multiwfn的主功能0绘制出来后手动ps到恰当位置的。
6.png
图中左侧是对应TDOS、PDOS曲线的轴,右侧是对应OPDOS曲线的轴。DOS或OPDOS的数值的绝对大小是我们不需要关心的,刻度轴的标签都抹掉也无妨,需要关心的仅仅是不同能量区间之间曲线的相对高度。为了便于直观考察不同能量区间内轨道的本质,绘制时定义了三个片段:
• 片段1:Fe的所有轨道
• 片段2:碳的垂直于茂环的p轨道。由于茂环在XY平面上,因此它们也即碳的pz轨道,合在一起相当于茂环的pi轨道
• 片段3:茂环的碳的其它轨道,即s、px、py
由图可见,对于能量比较低区间的轨道,诸如-0.53 a.u.的那一批,从图上立马就可以判断出它们基本都是由碳的s、px、py轨道贡献的,结合轨道等值面图也看出确实如此。而对于能量为-0.25 a.u.附近的轨道,蓝色曲线和红色曲线都较高,说明这些轨道应当是由茂环的pi轨道与铁的轨道混合产生的,确实从轨道等值面上也可以证实是由铁的d轨道与茂环的pi轨道以相位相同叠加组合出的。在HOMO轨道位置,Fe片段的PDOS与TDOS非常接近,说明HOMO主要是由Fe的轨道构成的,茂环的参与较弱。我们再看绿色的OPDOS曲线。在比如-0.25 a.u.附近,绿色曲线明显为正,说明这些轨道的存在对于茂环与铁的结合起到积极作用,这也容易理解,毕竟是相位相同叠加,所以对于此二者而言起到的是显著的成键轨道作用。对于非占据轨道区域,可见OPDOS曲线都为显著的负值,这也就是说,如果它们被电子占了,倾向于让茂环与铁的结合散架。确实从LUMO的轨道图来看,是Fe的d轨道与茂环的pi轨道以相位相反方式组合出来的,对它俩而言无疑起到的是反键轨道作用。

从上面这个例子看出,DOS图很有助于清楚、直观地讨论化学体系的电子结构。虽然可能有人觉得,用Multiwfn方便的轨道成份计算功能(主功能8,操作见http://sobereva.com/131和手册4.8节的例子),批量计算片段对各个轨道的贡献值也能得到本质上同样的信息,但是这样的话你得去分析一大堆数据,而且写文章的时候得把一大堆数据堆上去,审稿人和读者看着也累,远没有DOS图一目了然、远不如通过DOS图讨论起来方便。


2 Multiwfn中的DOS绘制功能

下面对Multiwfn的绘制DOS的功能做一下系统的介绍。

2.1 基本使用流程

将输入文件载入Multiwfn,进入主功能10,然后会看到一个菜单,里面有很多选项可以调节作图设定。设好后,按0,DOS图马上就可以显示出来。关闭图像后,会看到一个后处理菜单,里面又有很多选项来调节绘制参数,修改后可以选1 Show graph again重新绘图查看当前的效果。

如果绘图前不设定片段,则绘制出来的图是TDOS。如果先用-1 Define fragments设定片段再绘图,则绘制出来的是TDOS+PDOS。如果片段1和片段2都定义了,则这俩片段间的OPDOS会一同绘制出来。

2.2 输入文件

.fch、.molden、.gms文件都可以作为输入文件绘制TDOS、PDOS、OPDOS以及后文提到的LDOS,怎么生成它们见《详谈Multiwfn支持的输入文件类型、产生方法以及相互转换》(http://sobereva.com/379)。也可以使用Gaussian单点任务的输出文件(必须带着pop=full)作为输入文件,但此时只能绘制TDOS。.wfn、.wfx格式不能作为输入文件,因为这俩格式不包含空轨道信息。

为灵活起见,对于只需要绘制TDOS的目的,Multiwfn也支持通过文本文件作为输入,自行把轨道能级等信息按照要求的格式写进去即可,格式见手册3.12.2。因此,原理上任意量子化学程序计算出的轨道能级都可以被Multiwfn展宽成TDOS曲线绘制出来。

2.3 DOS绘制界面里的选项

进入主功能10后看到的菜单里面的选项在这里说一下,用途特别显而易见的就不提了。一般来说默认选项下绘制出的图不会恰好满足需要,因此需要反复调节选项、反复作图直到满意为止。

-5 Customize energy levels, occupations, strengths and FWHMs for specific MOs:默认情况下,轨道能量、占据数都是从输入文件里读的,如果你想修改,可以用这个选项。默认情况下,每个轨道的强度值(strength)都为1,FWHM都是相同的,如果你想对不同的轨道进行独立的设定,也可以用这个选项。展宽出的峰的积分面积与轨道的强度值是呈正比的。

-4 Show all orbital information:显示所有轨道的能量、占据数、强度值、FWHM信息。

-3 Export energy levels, occupations, strengths and FWHMs to plain text file:把所有轨道的能量、占据数、强度值、FWHM信息导出到当前目录下的orginfo.txt中。这个文件的格式是满足作为Multiwfn绘制TDOS的输入文件的要求的,因此你可以自行修改这个文件里的数据后作为输入文件。

-1 Define fragments:设定片段。这会在本文2.5节专门说。

2 Set energy range:这用来设置绘制出的DOS图的横坐标,即能量下限、上限、标签间隔。

4 Set scale ratio for DOS curve:设定DOS曲线乘的系数。这个数值越小,曲线的数值显然就越小,而离散的竖线相对于曲线就显得越高。因此如果你觉得作出的图里面离散竖线的高度相对于曲线太高或太低,可以用这个选项调节。

6 Choose orbital spin:对于非限制性开壳层波函数的情况才会看到此选项,用于选择在绘图时考虑哪些轨道。可以只考虑alpha轨道,或只考虑beta轨道,或alpha和beta轨道同时考虑。

7 Set the method for calculating PDOS:设定绘制PDOS时用的轨道成份的计算方法。

8 Switch unit between a.u. and eV:切换横坐标的单位,默认是a.u.,对于能量而言就是指的Hartree,而eV也很常用。1 Hartree=27.2114 eV。

2.4 绘制PDOS用的轨道成份

绘制PDOS用的轨道成份的计算方法支持下面四种:
• Mulliken方法:此为默认。计算速度很快,但可能给出缺乏物理意义的结果(片段贡献超过100%或为负值),对于空轨道的可靠性低于占据轨道,而且绝对不能用于有弥散函数的情况!
• SCPA方法:速度和Mulliken一样快,相对于Mulliken方法唯一的好处是轨道成份不会算出来超过100%或负值,但其它问题仍在。
• Hirshfeld方法:基于Hirshfeld方式的实空间划分来计算轨道成份,比Mulliken和SCPA方法更为昂贵,优点是对于空轨道和占据轨道都非常可靠、稳健,而且不怕弥散函数。
• Becke方法:和Hirshfeld方法各方面都差不多,结果一般也都相近,只不过对原子空间的定义不同。

使用Mulliken和SCPA方法时,定义片段时最小的单位是基函数,也因此可以等效实现把特定的原子轨道、特定角动量的轨道、原子轨道壳层等定义为片段,非常灵活。而Hirshfeld和Becke方法只能计算原子在轨道中所占成分,因此定义片段的时候最小单位是原子,而且此时无法绘制OPDOS。

2.5 片段的定义方法

进入前述的片段定义界面后,你会看到一个片段列表,有10个片段可以定义。定义几个,PDOS就会绘制出几条曲线。片段1和片段2同时定义时也同会对二者绘制OPDOS。

要定义哪个片段,就在这个界面里输入哪个片段的序号即可,之后会开始片段的定义。如果比如2号片段已经定义了,输入-2就可以把这个片段清空,使之处于未定义状态。输入两个片段序号比如2,5可以交换两个片段(这个设置是为了便于把想考察OPDOS的两个片段切换到片段1和片段2)。如果你的片段定义比较复杂,以后还可能要用,可以输入字母e将片段定义导出为当前目录下的DOSfrag.txt,下次再进入这个界面时可以输入i来从当前目录下的DOSfrag.txt中导入片段定义,省得再重新设一遍了。

用Hirshfeld或Becke方法计算轨道成份时,设定片段的内容时就按提示输入原子序号即可,比如1,3-5,12,15就代表把1,3,4,5,12,13号原子定义为片段。

用Mulliken或SCPA方法计算轨道成份时,片段里记录的是基函数的序号。开始定义片段后,会看到一个命令很丰富的设定界面,值得细说一下。在这个界面里,你首先在屏幕上会看到一大堆命令的介绍以及用法示例。如果你想让这些帮助信息重新显示出来,可以随时输入help。下面对可以用的命令都解释一下:
q:保存片段并离开定义当前片段的界面,屏幕上会提示片段里都加入了哪些基函数
list:显示当前片段里包含的基函数序号
clean:清空片段。默认情况下,片段里是没有内容的,如果之前已经对这个片段进行过定义,再次定义时原先那些基函数还会在,因此想彻底重新定义之前应当输入clean清空一下
addall:将所有基函数加入这个片段,即这个片段等同于整个体系
inv:反转定义。原先在片段里的基函数将不再片段里,而不在片段里的将被加入片段里
a:按照原子(atom)定义片段。比如a 2,5-8,12代表把2,5,6,7,8,12原子上的基函数都加入片段
s:按照基函数壳层(shell)定义片段。比如s 2,5-8,12代表把2,5,6,7,8,12号基函数壳层里的基函数都加入片段
b:按照基函数(basis function)序号定义片段。比如b 2,5-8,12代表把序号为2,5,6,7,8,12的基函数都加入片段
l:按照角动量定义片段。后面跟着角动量名字。比如l s,d就代表把s和d角动量的基函数都加入片段
da、db:分别按照原子序号和基函数序号来从当前片段中删除基函数
cond:将满足条件的基函数加入片段,条件包括基函数所处的原子的序号范围、基函数所处的序号范围、基函数类型(如S、Y、XZ)或者基函数所处的壳层的类型(如S、P、D)。只有同时满足这三个条件的基函数才会被加入片段。如果想无视某个条件,设相应条件时输入a(代表all)即可。

a、s、b、l这些命令若反复用多次,每次都会把满足要求的基函数添加到片段,换句话说片段里包含的会是它们的并集。而cond这种方式加入的基函数是自己设的三个条件的交集。上面说到的所有这些往片段里加入基函数的命令都可以反复使用多次,效果是累加的(某些基函数加重复了也没关系,程序会自动去除重复的)。

如果你想查看基函数、壳层、原子的对应关系和特征,随时可以用all命令,会输出诸如下面这些,此例是6-31G*用于HF分子。如果某个基函数已经在片段里了,前面会显示星号。如果屏幕上显示不全,参考手册5.5节来加大窗口缓冲区。
  Basis:     1    Shell:    1    Center:    1(H )    Type: S
  Basis:     2    Shell:    2    Center:    1(H )    Type: S
  Basis:     3    Shell:    3    Center:    2(F )    Type: S
  Basis:     4    Shell:    4    Center:    2(F )    Type: S
  Basis:     5    Shell:    5    Center:    2(F )    Type: X
  Basis:     6    Shell:    5    Center:    2(F )    Type: Y
  Basis:     7    Shell:    5    Center:    2(F )    Type: Z
  Basis:     8    Shell:    6    Center:    2(F )    Type: S
  Basis:     9    Shell:    7    Center:    2(F )    Type: X
  Basis:    10    Shell:    7    Center:    2(F )    Type: Y
  Basis:    11    Shell:    7    Center:    2(F )    Type: Z
  Basis:    12    Shell:    8    Center:    2(F )    Type: XX
  Basis:    13    Shell:    8    Center:    2(F )    Type: YY
  Basis:    14    Shell:    8    Center:    2(F )    Type: ZZ
  Basis:    15    Shell:    8    Center:    2(F )    Type: XY
  Basis:    16    Shell:    8    Center:    2(F )    Type: XZ
  Basis:    17    Shell:    8    Center:    2(F )    Type: YZ
上面显示的X,Y,Z即PX,PY,PZ型基函数,XX,YY,ZZ,XY,XZ,YZ都是笛卡尔型D型基函数的标识,不懂的话看《谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识》(http://sobereva.com/51)。根据6-31G*定义,我们可以立刻判断出哪些基函数对应哪些原子轨道。比如对于F,由于每个内层原子轨道用一个基函数描述,每个价层原子轨道用两个基函数一起描述。因此你输入b 3,就等于把1s原子轨道加入片段了;如果输入b 4,8或者等价的s 4,6,就代表把2s原子轨道加进去了;如果输入b 5-7,9-11或者等价的s 5,7,就代表把2p原子轨道加进去了;如果输入b 12-17或s 8,就代表把给F加的d极化函数加进去了。

对于非Pople系列基组,往往基函数与原子轨道的对应关系没法根据基组名字简单判断,但笔者也发明了一种简单且普适的判断方法,参见《利用布居分析判断基函数与原子轨道的对应关系》(http://sobereva.com/418)。

2.6 后处理菜单中的选项

在绘制出图像并关闭图像后,后处理菜单中也有很多选项可以设定绘制参数,这里把其中一部分提一下,其它的选项含义不言自明。

-1 Set format of saved image file:用来设置通过选项2保存出来的图像的格式。对于除了二维LDOS(见后文)以外的DOS图,由于图中的主体都是线条,建议使用矢量格式,比如pdf、ps、wmf(可以嵌入Office的组件中)、svg(可以用网页浏览器打开),这样线条看着光滑且可以无损缩放。如果想修改默认的格式,可以改settings.ini里的graphformat。

1 Show graph again:后处理菜单修改过选项后,可以用这个选项重新绘制图像检查当前效果。

2 Save the graph to image file in current folder:将图像文件保存到当前目录,文件名前缀为DISLIN。如果不懂什么叫当前目录,看http://sobereva.com/237

3 Export curve and line data to plain text file in current folder:这个选项可以把DOS曲线数据和离散竖线数据导出到当前目录下的文本文件中,之后可以导入Origin之类程序中画DOS图,届时有更多可调选项可用。

选项5~10:用来设定是否显示TDOS、各个PDOS、OPDOS曲线和离散的竖线。

11 Set color for PDOS curves and lines:设定各个PDOS的曲线和离散竖线的颜色

12 Set position of legends:设定图例的X和Y位置。X数值越大越靠右,Y数值越大越靠下。通过调节这个可以避免图例把曲线挡住。用选项13也可以完全关闭图例。

14 Set scale factor of Y-axis for OPDOS:对OPDOS设置倍数因子。OPDOS对应图中右侧坐标轴,其范围是对左侧坐标轴(对应于TDOS和PDOS)的范围乘上这个倍数因子得到。因此通过设定这个值可以控制OPDOS曲线相对于TDOS/PDOS曲线的高低。

15 Toggle showing vertical dashed line to highlight HOMO level:默认情况绘制出的DOS图上有个垂直的虚线,用来标注HOMO轨道的位置以便分清楚哪些DOS区域是占据的和非占据的。用这个选项可以关闭这个虚线的显示。老有人以为这个虚线是费米能级的位置,这完全是误解。费米能级是固体物理里的概念,对于孤立体系,费米能级本身就是个ill-defined的概念,对于孤立体系绝对不要用费米能级这个词。

16 Set the texts in the legends:设置TDOS、各个PDOS、OPDOS显示在图例中的文字

选项17、18:设置DOS曲线和离散竖线的粗细。注意线条在屏幕上显示的粗细可能和保存出的图像文件里的粗细不同。

19 Toggle showing labels and ticks on Y-axis:用这个选项可以关闭纵轴上的标签和刻度。前面说过,DOS图中的绝对数值大小没什么实际意义,不想显示的话可以索性用这个关了。


3 TDOS/PDOS/OPDOS在Multiwfn中的绘制例子

下面给出几个有代表性的DOS图的绘制例子。这些例子显然不可能把所有情况都体现,读者务必将例子与上面介绍的原理和选项说明相结合,举一反三。手册4.10节还有更多的例子。下面例子用到的所有文件都可以在此下载:http://sobereva.com/attach/482/file.rar,计算都通过Gaussian 16实现。

3.1 绘制二茂铁的各种DOS图

这一节我们以二茂铁为例演示一下绘制TDOS、PDOS、OPDOS,复杂度由浅入深。

这个例子的Gaussian输入文件是本文文件包里的Ferrocene.gjf,可见关键词是BP86/genecp opt,对Fe用的SDD赝势基组,对茂环用的6-31G*。绘制DOS不需要用太高质量波函数,用这样比较便宜的基组就够了。计算产生的fch文件是文件包里的Ferrocene.fch。

我们先看怎么绘制最简单的TDOS。启动Multiwfn,载入Ferrocene.fch,依次输入
10  //DOS绘制模块
0   //绘制TDOS
马上就看到下图
f1.png

下面我们把两个茂环中的所有碳原子作为片段1来绘制它们的PDOS。为此我们需要先知道这些碳原子的序号,可以通过Multiwfn主功能0查看,但原子多的话一个一个记录数字比较费事。一个简单的做法是把Ferrocene.fch载入GaussView(笔者用的GaussView 6.0.16),用Builder界面里的选择工具把这些碳都选中成为黄色,然后进入Tools - Atom Selection界面,如下所示文本框里直接就显示了这些原子的序号,将它们拷出来备用。
f2.png

关闭之前显示的TDOS图像,在Multiwfn里依次输入
0  //从后处理菜单退回到DOS绘制界面
-1  //设置片段
1   //定义片段1
a 1-3,7-8,12-16  //将茂环上的碳加入片段
q  //保存片段1
0  //退出片段设定界面
0  //绘制TDOS+PDOS图

看到的图像如下,红色曲线就是茂环的PDOS。
f3.png

这一节最后,我们画个较复杂的,重现一下本文前面给出的二茂铁的TDOS+PDOS+OPDOS图。但由于计算级别和Multiwfn的版本与当时绘制那个图的情况有所不同,所以结果得到的图像不会完全相同。

退回到DOS绘制界面,然后依次输入
-1  //定义片段
-1  //撤销掉原先的片段1的定义(先进入片段1然后输入clean命令清空片段也可以实现相同目的)
1   //定义片段1
a 6  //将Fe作为片段1
q  //保存此片段
2   //定义片段2。我们要把所有碳的pz轨道加入
cond  //以条件方式定义向此片段添加基函数
1-3,7-8,12-16  //茂环上的碳的序号
a  //不对基函数序号范围做限制
Z  //必须是PZ型基函数,即垂直于茂环的P型基函数(Multiwfn主功能0里有笛卡尔坐标轴,一看就知道茂环正好是在XY平面上)
q  //保存此片段
3   //定义片段3。我们要把所有碳的s,px,py轨道加入,需要重复用三次cond命令来实现(大家可以把下面这一串内容直接复制到Multiwfn窗口中执行,省得一行一行敲进去)
cond
1-3,7-8,12-16
a
S
cond
1-3,7-8,12-16
a
X
cond
1-3,7-8,12-16
a
Y
q  //保存此片段
现在片段都已经定义好了,可以输入e把片段设定导出到当前目录下的DOSfrag.txt中,以便之后随时在这个界面里通过输入i来恢复之前好不容易设好的片段。然后输入
0  //退出片段定义界面
2  //设定横坐标能量范围
-0.6,0.1,0.1  //能量范围-0.6~0.1 a.u.,标签间隔0.1 a.u.
0  //绘制TDOS+DOS+OPDOS
当前看到的图像如下
f4.png

此图已经基本有模样了,但还需要做一些调整。图中显得FWHM数值有点大,显得本来挺大的gap看起来都有点连上了,而且纵坐标轴的范围可以再改改,让图像主体区域更加充满画面。图例的文字也需要设一下。关闭图像并输入0退回DOS绘制界面,然后依次输入
3