第1届高级班到本届高级班相隔约1年时间,在这期间讲授者卢老师煞费苦心、耗费了巨量时间精力将高级量子化学培训班内容打造得尽可能完美、扩充了多达的1/3内容!!!幻灯片已经从第1届的1500页左右增加到了2000页左右!!!本届增加的内容完整罗列在下面。本届培训还同时修复了第1届培训内容中存在的很多bug。
• “CASSCF、多参考方法的原理与应用”部分,扩展到了350多页,比第一届增加了近一倍的内容!!!现在这部分对CASSCF和多参考方法相关背景知识已经讲得极为系统和详细,精心编排的充分涵盖各种情况的计算例子已达到了30多个,从最基础的CASSCF(2,2)到高精度的MRCI/MRCC再到最新的NEVPT4(SD)和ICE-NEVPT2全覆盖、从单原子到有机分子再到过渡金属配合物各种体系全覆盖、从化学反应能/势垒计算到激发能计算到振动分析再到NMR计算等等各种问题全覆盖、MO/NO/NBO/LMO/UNO/有限温度MO/AVAS等各种初猜轨道的利用全覆盖...各个例子都特意强调对相应情况如何恰当定义活性空间和产生初猜轨道。目前这一节的含金量和内容丰富程度说是CASSCF和多参考计算的葵花宝典也不算太夸张,着实是这类方法计算从初识到精通的近乎完美的学习途径!本届给这一节新添加的内容:
- 方法学部分:CASPT2思想的具体介绍;DCD-CAS(2)方法的介绍;MR-CEPA(0)方法的简介;多参考四阶微扰方法NEVPT4(SD)的介绍;MR-EOM-CC的介绍;short-range DFT (srDFT)方法的介绍;大幅扩展极高精度的多参考耦合簇(MRCC)的介绍;大幅扩展了MC-PDFT的介绍;判断MCSCF的波函数的稳定性的方法;计算化学反应的活性空间选取的讨论;不同多参考方法对活性空间要求差异的讨论;OO-CCD(T)用于静态相关不强不弱体系的价值的讨论;扩展了ICE-CI的介绍
- 程序使用部分:大幅扩展了ORCA做MRCI的介绍;ORCA中的MRCC、NEVPT4(SD)、MR-EOM-CC的用法介绍;ORCA中用AVAS方法产生初猜轨道;对ORCA做CASSCF中涉及的CI系数优化的细节和不收敛的解决的说明;ORCA中RAS和ORMAS空间的设置以及用户自主定义CI空间中的组态方法的说明;actorbs、ActConstraints等重要选项的使用说明;ORCA做ICE-CI、ICE-SCF、ICE-NEVPT2的说明
- 实例部分(若未注明都是ORCA做的):CASSCF优化1,4-脱氢苯的双自由基态;标量X2C-NEVPT2对Cr2做势能面扫描得到结合能曲线;SC-NEVPT2计算吖嗪的反转S1-T1能隙;NEVPT2计算Diels-Alder反应势垒和反应能的例子和活性空间选取的详细讨论;QD-NEVPT2的计算LiF的最低两个1∑+态曲线考察避免交叉的例子和结果的深入讨论;MRCI计算CH2的最低单/三重态能量随键角的变化;MRCI+Q计算CH2的特定不可约表示的态的能量随键角的变化;MRCI+Q计算笑气NNO的异构化能并结合此例对计算异构化能所用的活性空间问题做深刻的讨论;NEVPT4(SD)计算乙酸-Cu双核配合物的单-三重态能量差;MRCC计算静态相关极强的H4体系;MRCI、FIC-MRCC、MR-EOM-CC计算Al原子的3p-4s激发能;FIC-CASPT2-IPEA做N2势能曲线扫描;借助AVAS做二苯甲烷的pi-pi*激发的NEVPT2计算;借助AVAS做顺铂的d-d跃迁的MRCI计算并体现double-shell的考虑;使用ICE-CI对LiH做Full CI的近似计算;ICE-SCF计算五并苯(用到CAS(22,22)活性空间);ICE-NEVPT2结合CAS(18,18)计算18轮烯;用Dalton对苯分子在CASSCF级别下计算磁化率和磁屏蔽值;用Dalton对C16碳环的三重态在CASSCF级别下计算NICS(1)ZZ芳香性指数(重复卢老师在《全面揭示16碳环(cyclo[16]carbon)非常奇特的激发态芳香性!》http://sobereva.com/741介绍的论文中的计算)
• 新加入了“势能面交叉与光化学”一节,超过130页:不同电子态的势能面在特定几何结构下可以出现相交或接近相交,称为势能面交叉,这与光化学机理、电子激发态的无辐退激过程、涉及多态的反应的产物分布等方面有直接联系,是做非绝热过程研究的关键。本节的目的是让学员充分具备研究这些问题的能力,给出了涵盖各种情况的大量例子使得学员能够很容易地举一反三,并充分领会实际计算中的各种重要细节和关键要点。本节有三部分内容:
(1) 势能面交叉的相关知识:这部分内容目的是使得学员对势能面交叉及各种相关背景知识有充分、深刻的了解,讲授内容包括:势能面交叉的概念;绝热态与透热态;势能面完全相交的条件与避免交叉;势能面交叉区域的势能面结构、圆锥交叉;MECP与MECI的概念;MECP的搜索算法;MECI的搜索算法;不同理论方法描述势能面交叉的能力;光化学的概念;大量的扩展阅读文章;势能面交叉与光化学的相关专著盘点
(2) MECP搜索例子:这部分详细示例和讲解怎么用具体程序以各种方式搜索各种情况的MECP,包含以下内容
• Gaussian用户搜索MECP的方法说明
• sobMECP的用法介绍(sobMECP是本培训讲授者卢天编写的被广泛使用的搜索MECP的程序)
• 用sobMECP在DFT下搜索FeO+的四-六重态MECP实例
• 用sobMECP在DFT下搜索Pt配合物的单-三重态MECP实例
• 用sobMECP在DFT下搜索O3的S0-T1的MECP实例
• 用sobMECP在DFT下搜索氮宾分子环化的MECP实例
• 用sobMECP在TDDFT/DFT下搜索乙酰苯的S1-T1态的MECP实例
• 用Gaussian在CASSCF下搜索CH2的T0-S1的MECP实例
• ORCA搜索MECP的功能的介绍
• ORCA在DFT下搜索[CH3O]+的单-三重态的MECP实例
• ORCA在DFT下搜索乙烯二酮解离路径的MECP实例
(3) MECI搜索例子:这部分详细示例和讲解怎么用具体程序以各种方式搜索各种情况的MECI,包含以下内容
• Gaussian用CASSCF搜索1,3-丁二烯的S1-S2的MECI实例
• Gaussian用CASSCF搜索苯-盆苯之间光异构化的MECI实例
• Gaussian用CASSCF搜索光诱导噻吩开环的MECI实例
• ORCA搜索MECI的功能的介绍
• ORCA用SF-TDDFT搜索亚乙基乙叉(ethylidene)的S0-S1间MECI实例
• ORCA用SF-TDDFT搜索亚乙基乙叉的S0-T1间的MECP实例
• ORCA用SF-TDDFT搜索偶氮苯的光异构化的MECI实例
• ORCA用SF-TDDFT搜索二苯乙烯的光环化的MECI实例
• ORCA用TDDFT搜索尿嘧啶的S1-S2间的MECI实例
• GAMESS-US搜索MECI的功能实例
• GAMESS-US用MRSF-TDDFT搜索乙烯结构扭曲的S0-S1间MECI实例
• 新加入了“X光光谱的计算”一节(正在制作中,下为暂定内容):量子化学研究者研究的最多的电子光谱的是价层电子激发相关的光谱,而内核电子的激发(对应X光吸收光谱,XAS)和内核电子的电离(对应X光电子谱,XPS)也同样十分重要,能够提供原子所处状态、局部化学环境方面的关键信息。本节讲授X光光谱的各种相关背景知识,详细盘点各种计算X光光谱的方法并进行比较(如TDDFT、STEOM-CCSD、IP-EOM-CCSD、CASCI/NEVPT2、DFT/ROCIS、结合MOM做ΔSCF、Koopmans定理结合QTP系列泛函等),并给出以不同方法计算各种体系的XAS和XPS的实际例子,使研究者能够举一反三、轻车熟路地做这类计算
• 新加入了“real-time TDDFT模拟”一节(40多页)。real-time TDDFT (RT-TDDFT)通过求解含时Kohn-Sham方程,可以模拟体系的能量和电子结构特征随时间的实时变化。常用于模拟随时间变化的电场对化学体系的能量和电子结构的影响、模拟价层和X光电子光谱、考察电子转移等超快过程。本节详细介绍RT-TDDFT的各种相关知识以及通过NWChem做RT-TDDFT模拟的方法,内容包括:RT-TDDFT的原理和用途介绍;支持RT-TDDFT的程序盘点;RT-TDDFT的相关学习资料;使用NWChem做RT-TDDFT的输入文件写法详解;做RT-TDDFT模拟电子光谱的流程和细节的详细介绍;实例1:施加kick场模拟水分子的全范围电子光谱并与LR-TDDFT的结果进行对比,顺带强调RT-TDDFT具有的模拟吸收光谱在强场下表现的非线型特征的能力、强调RT-TDDFT算大体系完整光谱相对于LR-TDDFT的效率优势;实例2:模拟高斯包络脉冲电场令水分子产生的共振激发;实例3:模拟外加的电子在TCNQ二聚体中的分子间瞬时转移过程
• 新加入了“对角Born-Oppenheimer校正(DBOC)的计算”一节(12页)。DBOC是计算极高精度计算热化学数据(达到亚化学精度)必须考虑的一项。此节详细介绍了DBOC的概念、重要性和计算的背景知识,并演示怎么用ORCA和CFOUR程序分别在DFT和耦合簇级别下计算DBOC
• 原先的“极高精度耦合簇计算”更名为“使用MRCC和eT做基态和激发态的高精度耦合簇计算”,并从原先的20多页扩增到了56页。相对于上一届的变化:
- MRCC部分:加入了ADC(2)、SOS-CIS(D)形式的TD-双杂化、SOS-ADC(2)形式的TD-双杂化的使用说明和实际计算例子。其中示例的RS-PBEP86/SOS-ADC(2)方法算任何类型电子激发都表现得很好,很适合想获得比TDDFT明显更高精度、更稳健结果且希望计算不过于昂贵时使用。
- eT部分:eT是近年来新开发的特别擅长做耦合簇计算的免费、高效的程序,介绍文章见J. Chem. Phys., 152, 184103 (2020)。这一节用20多页幻灯片专门介绍eT程序的功能、安装与使用,重点演示通过eT程序使用比TDDFT稳健精确得多的Cholesky分解加速的LR-CC2方法以及使用极高精度的EOM-CC3方法计算激发能和振子强度。还演示通过不对称Lanczos链方法结合EOM-CCSD精确高效计算UV-Vis全范围光谱的做法。
• “ORCA程序的使用”部分,新加入的内容有:轨道优化的微扰和耦合簇计算示例;输出UNO轨道;通过MOM方法做ΔSCF方式的激发态计算的原理和实例;ORCA结合Multiwfn预测物质的颜色;计算动态极化率;使用openCOSMO-RS 24a溶剂模型;使用ORCA做多尺度计算的介绍;ORCA做ONIOM(wB97X-3c:GFN2-xTB)搜索姜黄素的H转移过渡态的实例;B3LYP-3c的简介;CSF-ROHF与GS-ROCIS方法介绍;TDDFT做激发态优化和振动分析的例子;TDDFT计算激发态偶极矩的例子;ORCA的json文件的意义和产生方法、示例将AUTOCI模块产生的CCSD(T)级别的单电子密度矩阵写入json文件;背景电荷使用的例子。
同时针对近期发布的ORCA 6.1版的改进做了很多更新,如:令幻灯片里对手册内容的介绍与新版本对应;此版本已能够解析地计算Raman活性而重新修改了例子;加入了CCSD(T)解析计算极化率的例子;加入了新支持的超极化率计算的例子;加入了利用RIGIDBODYOPT在优化时保持所有分子为刚性的例子;加入了利用RigidFrags和RelaxFrags在笛卡尔坐标下优化时分别令片段处于刚性和柔性的例子。
• “
从头算分子动力学”部分:加入了metadynamics的介绍和ORCA做metadynamics的AIMD的实际例子。加入了“AIMD程序的选择建议”。加入了ORCA做18氮环分子解离过程的模拟例子(
http://sobereva.com/725博文报道的卢老师的研究里的体系)
• “低标度耦合簇方法”部分:对C8F8做DLPNO-CCSD(T)的CPS外推、做DLPNO-双杂化泛函计算的例子。加入了能够令MP2、耦合簇实现线性标度的cluster-in-molecule (CIM)方法的具体介绍,并给出在ORCA中做CIM-DLPNO-CCSD(T)计算的实际例子
• “相对论量子化学”部分:增加了CASSCF和MRCI计算时考虑旋轨耦合的例子,包括SOC-CASSCF计算旋轨耦合造成F原子基态的能级分裂,以及SOC-MRCI计算OH自由基基态的旋轨耦合分裂能。
• “spin-flip TDDFT方法”部分:加入了GAMESS-US做MRSF-TDDFT的具体例子,加入了做MRSF-TDDFT的OpenQP程序的介绍
• “磷光速率的计算”部分:加入了用Dirac程序做二分量相对论的TDDFT得到磷光速率的例子。加入了响应函数的背景知识的介绍
• “DFTB、GFN系列方法与xtb程序的使用”部分:加入了g-xTB方法的介绍。加入了ORCA 6.1开始原生支持的(sp)GFN1/2-xTB方法使用的说明
• “振动电子光谱的计算”部分:加入了ORCA做振动电子光谱计算的实际例子
• “分子振动问题的高级计算方法”部分:加入了做VPT2的独立的程序pyVPT2的简介。加入了ORCA通过GVPT2计算对NMR做振动校正的例子
• “双光子吸收的计算”部分:加入了三光子吸收的介绍和Dalton做这种计算的实例
• “交换-相关泛函的调控”部分:加入了利用Multiwfn做ELF调控的具体例子
• “密度拟合、RI与COSX技术”部分:加入介绍被很多程序用于加速库仑作用计算的Fast multipole moment (FMM)技术,以及ORCA 6.1开始引入的加速巨大体系库伦计算的Bubblepole (BUPO)技术