中科科翼 专注模拟

一、激发态相关计算

（一）激发态的定义

- 在量子力学中，分子的状态分为基态和激发态。基态是分子在未吸收能量时所处的最稳定的电子状态。当分子吸收光子或其他形式的能量后，电子会从低能级轨道跃迁到高能级轨道，此时分子就处于激发态。例如，在荧光物质中，分子吸收紫外光后从基态跃迁到激发态，然后在一定时间内又回到基态并发射出可见光。

（二）激发态相关计算的目的

理解光物理过程：计算激发态能量、振子强度等性质，有助于解释光吸收、发射、荧光、磷光等光物理现象。例如，通过计算可以预测分子在什么波长下会吸收光，以及吸收光后的发射光谱特征。

研究光化学反应：许多化学反应是由光激发引发的，如光合作用、光催化反应等。激发态相关计算可以帮助我们了解这些光化学反应的起始步骤，即分子吸收光后的结构变化和能量变化，从而为研究光化学反应机理提供基础。

设计光功能材料：对于有机发光二极管（OLED）、太阳能电池等光功能材料，激发态计算可以指导材料的分子设计。例如，通过计算可以优化分子结构，使材料在特定波长下有更好的发光效率或者光吸收能力。

（三）激发态相关计算的方法

基于量子力学的方法

组态相互作用（CI）方法：该方法考虑了不同电子组态之间的相互作用。例如，在单激发组态相互作用（CIS）计算中，会考虑从基态组态到单激发组态的跃迁情况，从而可以计算激发态的能量和波函数。不过，CIS方法可能会高估激发态的能量，对于复杂体系计算量也较大。

含时密度泛函理论（TD - DFT）：这是一种广泛应用的计算激发态的方法。它以密度泛函理论为基础，通过考虑含时的外场（如光场）来计算激发态的性质。TD - DFT能够在相对较低的计算成本下，比较准确地计算激发态的能量、振子强度等参数，适用于较大的分子体系。

二、用计算模拟技术来计算化学反应机理

（一）计算模拟技术的基本原理

量子力学方法

从头算方法（ab initio）：这是一种从量子力学基本原理出发，不依赖实验参数来求解分子体系薛定谔方程的方法。它可以精确地计算分子的电子结构和能量，对于研究化学反应机理，能够提供非常准确的信息。例如，在研究小分子反应时，可以用从头算方法计算反应物、过渡态和产物的能量，从而确定反应的能垒和反应热。但是，对于复杂的大分子体系，从头算方法的计算量巨大，计算成本很高。

密度泛函理论（DFT）：它以电子密度为变量来描述分子体系的性质。通过选择合适的交换 - 关联泛函，可以在保证一定计算精度的同时，大大减少计算量。DFT在化学反应机理研究中应用广泛，例如在研究催化反应时，可以模拟反应物在催化剂表面的吸附、活化和反应过程。

分子力学方法：分子力学把分子看作是由原子通过化学键连接而成的力学系统。它通过经验力场来描述原子间的相互作用，包括键伸缩能、键角弯曲能、二面角扭转能和非键相互作用能等。分子力学方法计算速度快，但对电子结构的描述比较粗糙。它主要用于大分子体系（如蛋白质、聚合物）的构象分析和简单的动力学模拟，对于化学反应机理的精细电子结构变化描述能力有限。

（二）计算化学反应机理的步骤

1.构建反应物分子模型

根据实验已知的化学结构或者理论推测，构建反应物分子的三维结构模型。可以使用专业的分子建模软件来完成这一操作。例如，在研究有机化学反应时，需要准确地构建反应物有机分子的原子坐标、化学键长、键角等参数。对于复杂的分子，可能需要参考晶体结构数据或者其他实验测定的结构信息。

2.选择合适的计算方法和参数

根据研究体系的大小、性质和计算目的选择合适的计算方法（量子力学方法或者分子力学方法）。如果选择量子力学方法，还需要选择合适的基组（对于从头算方法和DFT）。基组是用于描述原子轨道的数学函数集合，不同的基组会影响计算结果的精度和计算量。例如，对于含有第一周期元素的小分子体系，常用的基组有6 - 31G(d)；对于过渡金属配合物，可能需要选择包含相对论效应校正的基组。同时，对于分子力学方法，需要选择合适的力场，如AMBER力场、CHARMM力场等，这些力场是根据不同类型的分子体系和研究目的而开发的。

3. 计算反应物分子的初始状态

计算反应物分子的各种性质，如能量、电荷分布、分子轨道等。这些信息可以帮助我们了解反应物的初始状态，为后续的反应过程分析提供基础。例如，通过分析反应物分子的最高占据轨道（HOMO）和最低未占据轨道（LUMO）的形状和能量，可以预测可能的反应位点和反应类型。在量子力学计算中，这些性质可以通过求解薛定谔方程或者基于电子密度的计算得到；在分子力学计算中，可以通过力场函数来计算分子的能量和几何结构相关的性质。

4. 寻找过渡态（TS）

过渡态是化学反应过程中的一个关键结构，它处于反应物和产物之间的能量最高点，代表着反应途径中的能垒。可以使用各种优化算法来寻找过渡态，如同步转变方法（QST）、二阶导数方法等。找到的过渡态需要通过振动分析来验证，真正的过渡态应该有且只有一个虚频，这个虚频对应的振动模式与反应坐标方向一致。例如，在研究酯化反应的反应机理时，通过寻找过渡态可以确定醇和酸在反应过程中形成酯键的中间状态，以及这个过程中化学键的断裂和形成情况。

5. 计算反应路径和能垒高度

通过内禀反应坐标（IRC）方法计算从过渡态到反应物和产物的反应路径，它可以追踪反应过程中分子结构和能量的连续变化。能垒高度是指过渡态与反应物之间的能量差，它是衡量化学反应难易程度的一个重要指标。能垒越高，反应越难进行。例如，在研究一个酶催化反应时，通过计算反应路径和能垒高度，可以了解酶是如何降低反应的能垒，从而使反应在温和的条件下快速进行的。

6.分析产物分子的性质

计算产物分子的性质并与反应物进行比较，以理解反应过程中发生的结构和性质变化。这包括能量变化、化学键变化、电荷分布变化等。例如，在一个氧化还原反应中，通过分析产物和反应物分子的电荷分布，可以确定电子的转移方向和程度，从而深入理解反应的本质。