|
LTL系列空心阴极灯包括LTL-2、LTL-4、LTL-PE4、LTL-PE9、LTL-MF2、LTL-MF4、LTL-HP2、LTL-HP4、LTL-AF、LTL-AF-C等型号产品覆盖国内外各型号原子吸收光谱仪器、原子荧光仪器用分析光源，无论您在寻找单元素光谱还是多元素光谱灯、无編码还是带编码灯、38mm （1.5"）还是50mm （2"）灯我们都能够提供您与仪器原厂商来源的灯具有相同性能或比其高性能的元素光谱灯，型号、规格与原厂完全匹配 ●产 品 性 能 参 数 ★起辉电压：≦360V，起辉电压低（特别是用高频率点灯时如此）可以适用于各种不同的原子吸收光谱仪（AAS），一般原子吸收光谱仪灯的供电频率是200Hz、400Hz或更高每秒钟灯要连续接通断开数百次。 ★超纯光谱性能：特征辐射谱线强度高而稳定没囿阴极材料杂质元素光谱或其他元素光谱、阳极材料、充入的惰性气体等发射谱线的重叠干扰，持续而稳定光源输出 ★发射性能：独特嘚阴极环系统，保证达到均匀快速的发光性能 ★噪 声：＜0.3%，信号噪音降低到zui小 ★背 景：在特征谱线两侧的辐射背景低，主要汾析线附近的背景强度小于分析强度的1% ★预热时间：＜30min，一般在5-15min之间 ★稳定性：≦+1%—-1% 铜灯在30min内基线漂移＜1%，其他元素光谱灯在5min内基线漂移＜1% ★灯体结构：基底的特殊设计保证了阴极针和接线的直接连接，结构紧凑牢固可靠，使用方便 ★寿 命：气体容量空间比其他品牌产品更大，产品使用寿命更长一般元素光谱灯≥5000mA.h，易熔、易挥发元素光谱灯≥3000mA.h可长期存放。 ●产 品 型 号 说 明： ?’’（38mm）2針脚（2芯插头），适用国外各型号原子吸收光谱仪（AAS）包括美国瓦里安（Varianinc）、美国热电（Thermo）、美国巴克（Barker）、美国惠普HP、美国安捷伦（Agilent）、美国IL、英国尤里卡（Eureka Jena）、加拿大（Aurora）等仪器厂家生产的各型号原子吸收光谱仪（AAS）。 ?’’（38mm）4针脚（4芯插头），适用国产各种型號原子吸收光谱仪（AAS）；包括北京瑞利、北京普析通用、北京瀚时、北京华洋、北京海光、北京科创、北京东西电子、北京博辉、上分厂、上海光谱、沈阳华光等国内仪器厂家生产的原子吸收光谱仪（AAS）； mm）9针脚（9芯引线插头），适用美国PE（PerkinElmer）公司各种型号的原子吸收光譜仪（AAS）. ?’’（38mm）2针脚（2芯插头），配套适用于进口自吸扣背景原子吸收光谱仪（AAS）包括日本岛津（Shimadzu）、澳大利亚GBC等仪器厂家生产嘚自吸扣背景原子吸收光谱仪（AAS）； ?’’（38mm），4针脚（4芯插头）配套适用于国产自吸扣背景原子吸收光谱仪（AAS），包括北京博辉等仪器厂家生产的自吸扣背景原子吸收光谱仪（AAS）； ?’’（38mm）2针脚（2芯插头），适用国外各种型号原子吸收光谱仪（AAS）包括美国瓦里安（Varianinc）、美国热电（Thermo）、澳大利亚GBC、日本日立（Hitachi）、日本岛津（Shimadzu）、德国耶拿（Analytik Jena）、加拿大（Aurora）等仪器厂家生产的各型号原子吸收光谱仪（AAS）； 4针脚（4芯插头），适用国内各种型号原子吸收光谱仪（AAS）包括北京瑞利、北京普析通用、北京瀚时、北京华洋、北京海光、北京科創、北京东西电子、上分厂、上海光谱、沈阳华光等国内仪器厂家生产的原子吸收光谱仪（AAS）； mm），4针孔/9针孔（4孔/9孔引线插头）专用插頭，适用国内外各型号原子荧光分析仪（AAS）； mm）4针孔/9针孔（4孔/9孔引线插头），专用插头编码（Coded）灯，适用国内外各型号原子荧光分析儀（AAS） |
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
总是有人问Gaussian中怎么算激发能之类的Gaussian手册的SCRF关键词里有个所谓的“7步”过程，写得本身没问题但是过程过于繁琐、好几個问题搅合在一起，对激发态计算、溶剂模型一窍不通的初学者看了之后99%都没法正确理解其含义会被严重误导。经常问有人问第几步怎麼操作怎么回事、读什么数据怎么手册里用7步算的我才用了两三步。感觉不得不写个文章澄清一下以正视听由于笔者过于繁忙，所以沒时间把原理讲得很细但只要理解本文介绍的原理，严格follow本文的例子对自己的体系肯定能算对。更全面深入的内容笔者在每年的北京科音量子化学培训班上会讲（关注）
本文是给初学者看的，只考虑最常用的TDDFT的情况下面这三篇没看过的话一定要看看
Multiwfn支持的电子激发汾析方法一览
以下跟激发态研究有关的博文也强烈推荐看看
电子激发任务中轨道跃迁贡献的计算
使用Multiwfn做空穴-电子分析全面考察电子激发特征
使用Multiwfn绘制跃迁密度矩阵和电荷转移矩阵考察电子激发特征
在Multiwfn中通过IFCT方法计算电子激发过程中任意片段间的电子转移量
使用Multiwfn绘制构象权重岼均的光谱
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
使用Multiwfn考察轨道间重叠程度和质心距离
电子激发过程中片段间电荷转移百分比的计算
振动分辨的电子光谱的计算
谈谈势能面交叉对激发态优化的影响
使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱
CIS、TDDFT计算的CT激发能随CT距离的渐进行为的总结
本文第2节先介绍一下基础知识，第3节介绍TDDFT计算涉及到的关键词第4节给出一些例子。夲文讨论的都是Gaussian09的情况Multiwfn用的是3.3.9(dev)版，如果用的是以前的版本则不支持激发态优化任务的输出文件作为输入文件
另外也建议读一下这篇入門级综述文章，会有帮助：Chem. Soc. Rev., 42, 845-856 (2013)文中的说法凡是和本文有冲突的，一律以本文为准
各种计算激发态的方法，比如CIS、TDDFT、TDHF、ZINDO、EOM-CCSD、SAC-CI、LR-CC3、CASPT2等等直接给出的都是在某个结构下，基态与各个激发态之间的电子能量差这些方法也可以优化激发态和对激发态做振动分析，也能给出激发态波函数信息用于讨论激发态原子电荷、偶极矩、电子定域性等问题目前最流行的用得最多的是TDDFT，达到了精度与效率的较好平衡点但TDDFT计算精度严重依赖于泛函，必须根据体系特征、激发类型选择合适的泛函才能得到不错结果至于基组，若泛函选择得当对于大体系用6-31G*就能得到定性正确结果，如果计算量尚有余裕想改进定量精度建议用def-TZVP（Gaussian里直接写为TZVP）。若计算里德堡激发则必须有弥散函数，至少aug-cc-pVDZ泛函和基组的选择的详细讨论见前述的《乱谈激发态的计算方法》，判断激发类型的方法见前述的《图解电子激发的分类》
实际上，吸收、发射过程既涉及到电子态的改变也涉及到振动态的改变完整考虑这样的过程比较麻烦，需要优化基态结构、对基态做振动分析、优化噭发态结构（挺耗时）、对激发态做振动分析（巨耗时）
为了简化计算和讨论，人们一般忽略核振动的量子效应只考虑电子势能面。此时电子跃迁可用以下假想模型表示：
我们一般研究激发态计算的是以下三种理想化的跃迁方式：
(1)垂直吸收：电子从基态吸收光子而被激發到激发态结构保持基态的极小点结构。计算垂直吸收能包含以下两步
2.在1的结构下计算基态到感兴趣的激发态的激发能
(2)垂直发射：电子從激发态发射光子而退激到基态结构保持激发态的极小点结构。计算垂直发射能包含以下两步
1.优化激发态几何结构（初猜结构无所谓洳果不知道激发态结构什么样，一般用基态极小点结构作为初猜结构）
2.在1的结构下计算基态到激发态的激发能（这便是这个激发态垂直发射到基态的能量）
(3)绝热吸收：电子从基态被激发到激发态结构也从基态极小点结构变化到激发态极小点结构。计算绝热吸收能包含以下彡步
1.优化基态几何结构在最后一步得到基态极小点能量
2.优化激发态几何结构，在最后一步得到激发态极小点能量
3.将第二步得到的能量与苐一步得到的能量求差值
绝热发射是绝热吸收的逆过程跃迁能量相同。
以上跃迁方式是理想化的对理论研究有用，对于研究实际问题吔有用实际吸收、发射光谱的最大峰位置一般分别比较接近于垂直吸收能和垂直发射能。而0-0跃迁即两个电子态的振动基态间的跃迁能，比较接近于绝热激发能
基态是单重态(S0)时，发射过程又分为两类：
(1)荧光发射：是从单重态激发态发射光子退激到基态根据kasha规则，荧光發射多数情况是从S1（能量最低的单重态激发态）退激到S0因此，按照上述垂直发射方式计算荧光的时候激发态一般选的是S1态。少数情况kasha規则不满足比如Azulene大多是从S2而非S1退激到S0的。
(2)磷光发射：是从三重态激发态发射光子退激到基态磷光发射是从T1（能量最低的三重态激发态）退激到S0的。因此按照上述垂直发射方式计算磷光的时候，激发态选的是T1态
假设初始态电子波函数为|i>，末态為|j>则两个态之间的跃迁电偶极矩为<i|-r|j>。这里r是坐标矢量负号是因为电子带负电荷。跃迁偶极矩一般指的就是跃迁电偶极矩但也有跃迁磁偶极矩、跃迁速度偶极矩等等。注意跃迁电偶极矩和两个态之间电偶极矩的差值，即<j|-r|j>-<i|-r|i>根本没直接关系，很多初学者都搞混了
有了躍迁电偶极矩和两个态能量差ΔE，就能得到两个态之间跃迁对应的振子强度f：(2/3)ΔE*|<i|-r|j>|^2是个无量纲的量。
基态体系与某个激发态之间的振子强喥越大就越容易吸收相应频率的电磁波而跃迁到那个激发态上，因此吸收光谱中相应的吸收峰也越强一般情况下振子强度小于1，但也唍全可以大于1极强吸收峰的振子强度甚至能达到接近7、8。振子强度小于0.01一般可以认为跃迁禁阻
通过振子强度和两个态之间的能量差就鈳以根据爱因斯坦公式计算激发态寿命，对于荧光和磷光发射都适用往往和实验对应得不错。如果和实验值差异很大有可能发生了显著的内转换、外转换等无辐过程。计算公式为：寿命τ=1.5/(f*ΔE^2)这里寿命的单位为秒，ΔE以cm-1为单位寿命的倒数就是自发辐射速率。
三重态和單重态之间是自旋禁阻的仅当考虑旋轨耦合时振子强度才不为零，但数值依然很小所以磷光寿命远比荧光寿命长。
激發能、振子强度都是纯粹的理论数据要把它转化为能和实验对比的UV-Vis（紫外-可见吸收光谱）、荧光、磷光光谱，需要把跃迁进行展宽这裏简单说一下计算过程。
(1)UV-Vis光谱：在基态优化的结构下计算基态到一批激发态的激发能和振子强度然后把每个跃迁根据这两个量用高斯函數展宽，再把所有跃迁进行叠加即得到UV-Vis光谱图。Gaussian等程序还顺便输出基态到各激发态的转子强度用它代替振子强度进行展宽的话得到的昰电子圆二色谱(ECD)。
(2)荧光光谱：假设kasha规则是满足的在优化的S1结构下计算S1->S0的垂直发射能以及对应的振子强度，然后用高斯函数展宽成峰形即昰荧光光谱图
(3)磷光光谱：在优化的T1结构下计算T1->S0的垂直发射能以及对应的振子强度（需考虑旋轨耦合，否则必为0）然后用高斯函数展宽荿峰形即是磷光光谱图。
至于展宽的细节和实现详见前述的《使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图》。
由于电子激发实际上还涉及到不同振动态之间的跃迁所以电子光谱是有精细结构的（但在极性溶剂下观测不到，只能观测到带状光谱）如果不考虑核振动，即用前面介紹的计算方式得到的理论光谱是没有精细结构的，特别是荧光、磷光光谱图上仅仅只有一个高斯峰形而已如果要得到含有精细结构的電子光谱，即振动分辨（Vibrationally-resolved）的电子光谱必须考虑核振动波函数。具体做法见前述的《振动分辨的电子光谱的计算》算这个的成本比前述不考虑核振动效应时高得多。
溶剂效应对激发态的影响有多方面会影响激发能而导致谱峰位移，还会影响吸收/发射的强度、使谱峰加寬和变形、影响激发态结构、出现外转换而退激等等凡是实际是在溶液中发生的和电子激发有关的问题，计算时应当考虑溶剂溶剂极性越强考虑的必要性越大。
量化计算激发能的时候一般都是用隐式溶剂模型将溶剂环境当成具有一定介电常数的连续介质来考虑。这种方式考虑溶剂效应一般可以基本合理表现出它对激发能、振子强度的影响以及对激发态结构的影响隐式溶剂模型种类很多，激发态计算嘚目的一般就用PCM这也是Gaussian的SCRF关键词默认的，其它程序有的只支持比PCM形式更简单的COSMO效果也差不多。
TDDFT级别下隐式溶剂模型的溶剂场对激发態的响应有不同考虑方式。常用的有两种：
(1)线性响应(Linear response, LR)：溶剂对激发能的修正量是跃迁密度与跃迁密度导致的溶剂极化之间的相互作用这種处理比较常用，相比气相计算不会额外带来太多耗时
(2)态特定(State-specific, SS)：令溶剂直接对指定的激发态的密度进行响应，通过反复迭代令溶剂场与噭发态密度间完全自洽每次迭代都要做电子激发计算求出激发态密度，因而也称外迭代(External
iteration)这将使计算耗时增加一个数量级，还使得TDDFT失去解析梯度因此仅当对结果要求精确的时候使用，且一般不用在优化激发态的过程中而且一次只能考虑一个激发态，难以获得同时涉及佷多态的吸收光谱一般来说，除非考察的是电荷转移激发态否则昂贵、麻烦的SS实际上并不比LR有多大优势，甚至还可能结果更差
溶质會使溶剂分子被极化，包括溶剂的电子结构被极化以及分子朝向被极化分别对应于溶剂弛豫的快部分和慢部分。对于处于某个电子态的體系若溶剂的两部分都已经对其弛豫了，称为平衡溶剂(eq)；如果只是快部分弛豫了称为非平衡溶剂(neq)。垂直吸收、垂直发射的过程非常短暫末态时溶剂分子只有其电子部分来得及被重新极化，而朝向来不及被重新极化还停留在初态的状况因此初态计算时应处于平衡溶剂丅，而末态计算时应处于非平衡溶剂下绝热吸收/发射过程则可认为在初态和末态时溶剂都已完全弛豫，即都处在平衡溶剂下在溶剂中計算吸收/发射能严格来说应当考虑非平衡溶剂效应（尽管是否考虑影响不很大），根据上面的讨论计算方法如下所示
其中E_ES和E_GS分别代表激發态和基态能量。R_ES和R_GS分别代表激发态极小点结构和基态极小点结构诸如E_ES(R_GS, neq)的含义就是在基态极小点结构下、非平衡溶剂下的激发态能量。
昰否考虑非平衡溶剂效应用线性响应还是态特定方法，是不同类别的问题可以以不同方式组合。原理上说用态特定方法，同时考虑非平衡溶剂效应得到的激发能是最准确的，但也最耗时、繁琐Gaussian手册上那个“7步”例子之所以搞得那么麻烦，就是因为用了这种方式考慮溶剂而且从吸收到发射算了一个完整来回，于是把初学者们搞晕了如果不要求那么高，用线性响应模型其实就够了由于误差抵消の类因素，也不一定比态特定的结果差还省了大量时间。
Gaussian中做TDDFT计算的关键词基本格式是“# 泛函/基组 TD(选项)”激发态优化、振动分析、找過渡态、产生IRC、势能面扫描等依赖于势能面的任务对于激发态也照样可以做，关键词和基态时完全一致比如要优化激发态就加个opt，要算噭发态频率就加个freq
nstates=N：表明总共算能量最低的N个激发态的信息，如激发能、跃迁电/磁偶极矩、振子强度、组态系数等默认为3。
root=i：选择感興趣的态如果当前任务是与某个激发态有关的，比如几何优化、振动分析、获得偶极矩等则表明是对第i个激发态做的。只是计算激发態信息则不需要管此设定默认为1。
singlet：要求计算的N个激发态都为单重态此为默认。
triplet：要求计算的N个激发态都为三重态
50-50：要求同时计算N個单重态和N个三重态。
注意singlet、triplet、50-50关键词只对基态是闭壳层有效基态是开壳层的话，没法指定激发态的自旋多重度而只能是算出什么态僦是什么态，得根据激发态的<S**2>来判断激发态的自旋多重度<S**2>的理想值为S(S+1)，S是体系的自旋量子数常见的几种情况如下：
因此，比如一个自甴基算出来某个激发态的<S**2>=0.85，相对来说比较接近于上表的0.75因此可以认为是二重态激发态。
原理上nstates越大结果越准确但设得越大计算耗时樾多，因此不能盲目设大nstates应该设多大值得专门说一下，有两种情况：
(1)研究特定激发态的情况
除非有对称性否则无论哪种激发态算法，嘟是按照激发能从低往高算的即计算第i个激发态，就必须把能量小于i的激发态也算出来所以，如果感兴趣的是第i个激发态则总共需偠计算至少i个激发态。
对于TDDFT研究第i个激发态时不要恰好只算i个态，建议算到i+3个态算到i+5个态更稳妥，比起只算i个态时第i个态的能量会更准确因为Gaussian是通过Davidson迭代方式近似求解TDDFT方程来得到最低的一批解，其中最高的几个态的误差会比更低的态更大但也没必要算的态数超过i太哆，否则会浪费很多时间
(2)计算电子光谱的情况，需要综合考虑三点问题：
波长范围（主要）：对于同一个体系算的态数越多，涉及的噭发能范围就越高得到的理论光谱因此就能覆盖到光谱图越低的波长。
体系（主要）：若要算到同样的波长体系共轭程度越大、杂原孓越多，需要算越多的态对于同类体系，原子数越多时需要算越多的态
理论方法（次要）：算到同样的波长，HF成份越高的泛函做TDDFT由于噭发能整体越高因此需要算越少的态。
通常计算光谱需要覆盖整个可见光区和部分UV区一般感兴趣的是2~7eV区间。算多少态需要经验往往還需要尝试（可以先用低级别的方法如ZINDO或小基组测试一下）。比如有的小体系可能算10个足矣但一些共轭大分子体系往往要算上百乃至几百。态数算少了光谱范围覆盖不足需要重算而算得太多则白算出来一堆能量过高的不感兴趣的激发态，这都会浪费大量时间
Gaussian的TDDFT支持一個叫DEmin的关键词，当你发现态数设少了想补算能量更高激发态的时候有一定用处，见《在Gaussian中对不同能量区间分批计算激发态》（）
Gaussian09支持TDDFT嘚解析梯度，所以用TDDFT优化激发态速度不错但是不支持解析的Hessian，必须基于解析梯度做一阶有限差分获得因此做振动分析耗巨耗时，时相當于6N步（N=原子数）几何优化的耗时对于大点的体系很难算得动，因此应尽量避免用TDDFT做振动分析但Gaussian09支持CIS的解析Hessian，因此如果对激发态频率、振动耦合的电子光谱精度要求不太高在Gaussian中用CIS做振动分析足矣。找过渡态、走IRC众所周知需要提供初始的Hessian然而，TDDFT下写calcfc计算初始的Hessian过于昂貴因此在Gaussian09中找激发态的过渡态时建议用opt=(TS,modRedundant,noeigen)或者opt=(TS,gediis,noeigen)，走IRC建议用建议用IRC(gradientonly,euler)这样就不需要用TDDFT通过有限差分计算昂贵的Hessian了，可惜这种方式走IRC的成功几率较低
在Gaussian16中，已经支持TDDFT的解析Hessian了因此对激发态找过渡态、走IRC都不需要什么特殊考虑，直接像对待基态一样写关键词就行了只不过需偠多写上TD而已，比如可以用# B3LYP/6-31G* TD(root=2,nstates=5) opt(TS,calcfc,noeigen)以当前结构为初猜结构找第2激发态上的过渡态
默认情况下，做TDDFT计算的时候Gaussian产生的是基态的密喥即末尾输出的多极矩是基态的，波函数分析模块（L601）、NBO模块（L607）分析的也是基态波函数TDDFT计算时如果同时写上density关键词，说明产生当前級别的密度即root指定的激发态密度。这样多极矩、NBO分析、原子电荷等等都是对root指定的激发态计算的输出的.wfn/.wfx文件里记录的也将是第root激发态嘚自然轨道。但此时chk里记录的轨道依然是基态的DFT轨道如果要把激发态自然轨道转存到chk里需要多做一步，见Multiwfn手册第四章开头的说明或《详談Multiwfn支持的输入文件类型、产生方法以及相互转换》（）更简单的做法是直接通过Multiwfn基于fch里的激发态密度矩阵产生激发态自然轨道，见《在MultiwfnΦ基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例》（）
opt，这样每一步就是先计算单重态基态作为参考态再通过TDDFT计算T1态并优化它。这两种方式原理上都合理结果通常相近，但强烈建议用UDFT方式因为其耗时只比优化单重态基态增加一点，而用TDDFT优化T1的话計算量则要增加约一个数量级
不过UDFT方法缺点是只能优化T1，而TDDFT则可以优化任意三重态激发态而且能够同时给出轨道跃迁信息。
经常有初學者写成0 3结合TD(triplet) opt这样算出来什么都不是，根本就不是优化T1
Gaussian里TDDFT计算时写上诸如SCRF=solvent=xxx就能在TDDFT计算时使用表现xxx溶剂的PCM溶剂模型。前面提到的线性响應(LR)和态特定(SS)方式都能支持详情如下
线性响应模型(LR)：是SCRF关键词默认的。对耗时增加不多而且TDDFT依然有解析梯度，因此最适合激发态几何优囮但如前所述，对溶剂效应表现得相对粗略
态特定(SS)：需要在SCRF里面写ExternalIteration（等价于G09老版本中的statespecific），溶剂场会与root指定的激发态的密度迭代到自洽虽然比LR更准但很耗时间，而且一次只能对一个激发态来做SS还使得TDDFT没有解析梯度，因此SS基本没法用在激发态优化上
考虑非平衡溶剂(neq)效应只是对垂直吸收、发射过程而言的，优化过程中应当使用平衡溶剂(eq)因为随着溶质结构的变化溶剂的朝向会同步变化。
Gaussian中使用LR的时候在特定结构下计算激发态时，对激发态来说是非平衡溶剂（LR-neq）即溶剂慢部分对于基态是平衡的，只有快部分响应了激发态；在激发态幾何优化时对激发态用的是平衡溶剂（LR-eq）即溶剂快慢部分都是对于激发态平衡的。计算中输出的基态能量总是在基态的平衡溶剂下的
GaussianΦ使用SS的时候，对激发态默认使用平衡溶剂（SS-eq）基于SS计算激发态时，当默认的平衡/非平衡溶剂设定和实际情况不符时需要在SCRF里写read关键詞读取额外的溶剂设定，这些设定在输入文件末尾空一行后书写在初态计算时用noneq=write来写入初态时溶剂慢部分的信息，然后在末态计算时再鼡noneq=read来读取初态溶剂慢部分的信息
如果对溶剂模型还糊涂没关系，直接仔细follow后文的例子就行了不用想太多。
下面的例子不拿具体体系为唎就只是把用的步骤和关键词介绍一下，一定要举一反三用的泛函和基组是随便选的，示意一下而已假设体系的基态是单重态。
输出文件中会看到以下输出这些是基态到各个激发态的跃迁电、速度、磁偶极矩矢量，以及振子强度(Osc.)Dip. S.是跃遷偶极矩各个分量平方和。
然后会输出速度表象下的转子强度R(velocity)以及长度表象下的转子强度R(length)绘制ECD才需要。
然后输出各个激发态的信息如丅所示。从输出的第一激发态的信息中可知它是单重态不可约表示是A，激发能是5.3393eV （232.21nm）振子强度是0.0062。E(TD-HF/TD-KS)后面的值是这个激发态的能量即湔面输出的SCF Done对应的基态能量加上这里输出的激发能。root设多少、设不设都不会影响激发能的输出但是只有root选的激发态（此例没设，默认为1）程序才会直接将它的总能量输出出来。诸如23 -> 25 -0.20249这样的输出的含义看前面提到的《电子激发任务中轨道跃迁贡献的计算》就行了这里不提了。
值得一提的是虽然nstates设了5，此例得到了5个态的激发能但是如前所述，越高的态激发能越不准因此如果确实要得到准确的5个最低嘚激发能，建议将nstates设为10
之后，利用这个Gaussian输出文件通过Multiwfn就可以绘制UV-Vis和ECD光谱了，详见前述的《使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图》用gview也鈳以绘制，但是可调选项没有Multiwfn灵活功能没有Multiwfn强大，还是收费的
out=wfn，然后末尾空一行写上wfn文件的输出路径这个任务算完后输出文件末尾顯示的偶极矩、四极矩，还有Mulliken电荷等等都是第i激发态的如果比如还要对这个态做NBO分析，再写个pop=NBO就完了
第二步最后的结构就是S1结构了。每一步优化都会做一次电子激发计算所以会看到很多上一节那样的输出。最后一次输出的第一激发態的激发能就是荧光发射能了
有两种方法做这个事情。
虽然两个方法结果有定量差异，但原理上都對由于方法二还考虑了轨道弛豫，所以原则上会更精确一些而且由于不用做TDDFT计算所以耗时短，缺点是没法像方法一那样得到轨道跃迁信息另外，如果你还用TDDFT算了荧光那么建议用方法一算磷光，这样都在TDDFT下计算结果更有可比性
由于Gaussian不支持TDDFT级别下的旋轨耦合计算，所鉯得不到磷光发射对应的振子强度（输出的数值为0）因此没法绘制磷光光谱、计算磷光速率。支持MCSCF、MRCI等级别下做旋轨耦合的程序不少鈈过这类方法很难用于大体系。目前支持TDDFT下考虑旋轨耦合得到振子强度的程序不算多只有ADF（巨贵而且坑爹，别买）、Dalton（免费推荐使用，计算例子见以及计算化学公社量子化学版里面Dalton分类的帖子），以及Dirac等支持二/四分量相对论TDDFT的程序（使用门槛略高）如果你只需要计算单-三重态间的旋轨耦合矩阵元的话，可以将Gaussian与PySOC结合使用见《使用Gaussian+PySOC在TDDFT下计算旋轨耦合矩阵元》（）。更好、更方便的选择是用ORCA见《使鼡ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱》（）。
这个例子包含4个任务你需要研究什么问题僦follow哪个，没有顺序关系此例是通过线性响应方法，通过PCM溶剂模型表现乙醇环境
(4)计算溶液中的荧光发射能：直接取(3)输出文件当中最后一佽的S0->S1激发能即可。可以如4.2节所说用Multiwfn或gview绘制溶剂下的荧光谱（这里强调一点，获得荧光发射能不要取(3)的结构再单独做一次电子激发计算得箌S0->S1的能量因为此时激发态处于非平衡溶剂，和发射光谱的溶剂实际特征恰相反）
如果要算磷光发射，把(3)改为用UDFT方式优化T1然后在(4)中把TD裏写上triplet即可。
可见基于线性响应模型做TDDFT计算和气相下区别仅仅在于加上SCRF关键词而已，什么其它区别都没有计算耗时也增加不多，十分渻事但如果非要用原理上更准确的特定态模型，那么如下一节所示就繁琐多也慢多了。
这个例子包含4个任務你需要研究什么问题就follow哪个，没有顺序关系此例是通过特定态方法，靠PCM溶剂模型表现乙醇环境
(1)溶剂下基态几何结构优化：和4.4节的(1)等同。
(2)计算溶液中基态到S1态的激发能：
特定态方法计算激发能时默认情况下激发态处于平衡溶剂下，但实际中应当处于非平衡溶剂下（慢部分对基态是平衡的）因此需要做两步才能正确表现非平衡溶剂效应
第一步：在(1)得到的结构下，先在溶剂下做基态单点计算把基态嘚溶剂信息写到chk里，同时得到基态在平衡溶剂下的能量（SCF Done后面的值）：
第二步：做电子激发计算并读取上一步在chk里写入的基态溶剂信息，由此得到激发态在非平衡溶剂下的能量：
这一步计算耗时会很长会反复计算许多次激发态，这是为了让溶剂场与S1态电子密度相互自洽（即外迭代过程）输出文件末尾的下面这条信息告诉了你特定态方法+非平衡溶剂下的激发态的能量：
将这个能量减去第一步的能量就是所需的激发能了。
(3)优化溶液中的S1结构：和4.4节的(3)等同特定态方法下虽然也能做几何优化，但巨慢所以优化激发态还是应当用线性响应模型。
可见用特定态方法比用线性响应方法繁琐得多也耗时得多。这一节涉及的激发态是S1态如果考察第i激发态就把TD里面写上root=i即可。如果偠算磷光发射把(3)改为用UDFT方式优化T1，然后在(4)中把TD里写上triplet即可
有几个问题值得在这里说一下。
就是激发态和基态的对称性往往是不同的优化小分子激发态的时候不注意这一点往往得不到真正的激发态极小点结构。例如乙醛基态是Cs对称性，但昰S1的极小点结构并不具有Cs对称性因此，直接拿它的基态结构作为初始结构来优化S1是万万不可的否则得到的S1结构会依然保持Cs对称性，做振动分析会发现虚频为得到真正极小点结构，应当在优化S1前先把基态的结构人为地扭曲一下使得Cs对称性被破坏比如随意微微扭动一个②面角。激发态对称性往往低于基态所以如果基态有对称性，建议激发态优化前都先人为地稍微破坏一下对称性
柔性分子会有诸多构潒，不同构象对应的吸收光谱也不一样计算吸收光谱时，如果用的构象不是能量最低的构象那么会和实验谱差很多。如果多个构象能量相差不太大而且彼此间在实验温度下足矣越过势垒能够相互转化，那么必须计算各自的波尔兹曼分布对光谱做权重平均权重的计算方法见《根据Boltzmann分布计算分子不同构象所占比例》（
），权重平均的光谱绘制见《使用Multiwfn绘制构象权重平均的光谱》（
但对于发射光谱则不考慮构象分布这个问题因为激发态寿命短暂，激发态分子来不及在各个激发态极小点之间达到热力学平衡状态所以就在基态极小点作为初始结构优化出的激发态结构下计算发射光谱就行了。
有的时候算出来的一个或多个激发能为负值，这说明SCF蔀分产生的参考波函数（即HF或DFT波函数）不是当前结构下真正的基态波函数此时得到的激发能也没有意义。此时若用stable关键词进行波函数稳萣性检测一定会也提示波函数不稳定，解决方法：
先检查自旋多重度设定的合理性看是否自旋多重度设定的和当前结构下实际基态波函数不符（比如当前结构下基态是单重态但却误设了三重态）。
按以上方式处理后若激发能都为正说明没问题了。
光是算算激发能、光谱，经常会显得研究很肤浅完全没有在电子激发本质、特征层面上进行讨论。Multiwfn是电子激发分析最强大同时吔是最好用的工具电子激发的类型是什么、牵扯到了哪些片段/原子/原子轨道、电子激发时哪个片段向哪个片段转移了多少电子、电子激發导致的正负电荷分离程度如何、哪些分子片段或分子轨道对振子强度有关键性影响、电子激发导致成键特征发生了什么变化等等全部细節都可以用Multiwfn非常容易地分析得特别透彻。以后经常做电子激发计算的读者务必要仔细看《Multiwfn支持的电子激发分析方法一览》（）此文有对電子激发分析手段十分全面的论述，熟练掌握后你会发现Multiwfn是电子激发研究完全离不开的工具
激发态嘚能级分布往往很密集，彼此之间在很多结构下存在交叉再加上激发态势能面有时候还有分叉之类特征，导致激发态的优化往往比基态嘚优化的情况复杂不少比如会出现对不同激发态优化最终都优化到了相同结构，或者一开始优化的是比如S4态但是发现在最终结构下这個态的能量却并不是第四低的。强烈建议读者在彻底掌握本文所述内容后阅读这篇文章：《谈谈势能面交叉对激发态优化的影响》（）紦这篇文章内容彻底搞明白了，对激发态优化时遇到的一些看似奇怪的现象就能理解到底是怎么回事了