3. 光谱预测的物理量映射

量子化学计算给出的激发能、跃迁偶极矩等量，本身并不是实验直接测量的对象。实验上我们测到的是吸收峰位置、摩尔吸光系数 ε、荧光量子产率、圆二色（CD）谱的 Δε 曲线、以及积分吸收截面 σ(ν)。因此，光谱预测的核心任务不是"算出一个数"，而是建立一条从波函数输出到实验可观测量的映射链。这条链上任何一环的约定（单位、方向、简并度、展宽模型）出错，都会让最终谱图与实验产生系统性偏差。

一、映射链的起点：跃迁偶极矩

对初态 |i⟩ 与末态 |f⟩，电偶极跃迁矩定义为

μ_if = ⟨i| μ̂ |f⟩ = -e ⟨i| Σ r̂_k |f⟩

其中 μ̂ 是电偶极算符，求和遍及所有电子。在原子单位下，长度表象的跃迁偶极矩单位为 e·a₀（约 2.5417 D）。实际计算中常用三种等价表象：长度表象（r）、速度表象（∇）和加速度表象。对精确波函数三者严格相等，但对有限基组或近似泛函，三者数值会有差异。经验上，长度表象对基组更敏感但收敛后更可靠；速度表象对近简并态更稳健。若三者差异超过 10%，通常提示基组不足或存在强组态混合，需要检查活性空间。

跃迁偶极矩是矢量，其分量 μ_x、μ_y、μ_z 决定了跃迁的偏振方向。对取向固定的分子（如晶体中或取向样品），偏振方向直接决定哪些峰在特定偏振光下出现；对气相或溶液中的自由取向分子，则需对各向同性平均，得到标量形式的偶极强度。

二、从跃迁偶极矩到振子强度

振子强度 f 是无量纲量，表示跃迁的"强度权重"，其长度表象表达式为

f_if = (2/3) · ΔE_if · |μ_if|²

这里 ΔE_if 以 Hartree 为单位，μ_if 以 e·a₀ 为单位，系数 2/3 来自对三个空间方向的各向同性平均。若使用 SI 单位，则

f_if = (8π² m_e ν_if / (3 h e²)) · |μ_if|²

振子强度是连接计算与实验吸收谱的关键桥梁：在稀溶液近似下，积分吸收强度与 f 成正比。对高斯型吸收带，有

∫ ε(ν) dν ≈ (π e² N_A / (2 ln10 · m_e c²)) · f ≈ 1.13×10²⁰ · f  (单位 L·mol⁻¹·cm⁻²)

这意味着一个 f = 0.1 的跃迁，其积分摩尔吸光系数约为 1.1×10¹⁹ L·mol⁻¹·cm⁻²，对应峰值 ε 在 10⁴ 量级（假设半高宽约 3000 cm⁻¹）。这个数量级判断非常实用：若计算给出 f = 0.001 却预测出 ε ~ 10⁵ 的强峰，一定是单位或展宽设置出了问题。

三、旋转强度与手性光谱

对圆二色光谱，需要的是旋转强度 R（rotational strength），它同时涉及电偶极和磁偶极跃迁矩：

R_if = Im[ μ_if · m_if ]

其中 m_if = ⟨i| m̂ |f⟩ 是磁偶极跃迁矩，单位为 Bohr magneton（μ_B）。R 的单位是 e·a₀·μ_B，常用 10⁻⁴⁰ erg·esu·cm 或 cgs 单位表示。R 的符号决定 CD 峰的正负，绝对值决定峰强。与振子强度不同，R 对分子构型极其敏感——构型翻转或溶剂诱导的构象变化可能使 R 变号，这正是 CD 谱用于确定绝对构型的物理基础。

需要特别注意：磁偶极算符包含轨道角动量项，其原点选择会影响 m_if 的数值。对精确波函数，R 与原点无关；但近似计算中必须使用规范原点（通常取分子质心或电荷中心），否则 R 会出现虚假的大值。这是 CD 计算中最常见的陷阱之一。

四、从强度到截面：展宽与线型

实验测得的吸收截面 σ(ν) 是分立跃迁经展宽后的包络。对每个跃迁，常用高斯或洛伦兹线型：

σ(ν) = Σ_if (π e² / (2 m_e c ε₀)) · f_if · g(ν - ν_if)

其中 g 是归一化线型函数。高斯线型适合非均匀展宽（如气相多普勒展宽、固体中的非均匀应变），洛伦兹线型适合均匀展宽（如溶液中的碰撞展宽、寿命展宽）。实际溶液谱常用 Voigt 线型，即两者卷积。展宽参数 HWHM 通常取 0.1–0.3 eV（约 800–2400 cm⁻¹），需根据实验谱的峰宽反推，而不是随意设定。

发射截面则涉及 Einstein A 系数：

A_if = (16π³ ν_if³ / (3 ε₀ h c³)) · |μ_if|²

辐射寿命 τ = 1/Σ A_if。若已知荧光量子产率 Φ，则发射截面 σ_em(ν) 可由 A_if 和归一化发射线型得到。这里的关键是：吸收截面和发射截面通过振子强度和 A 系数共享同一套跃迁矩，因此吸收谱和发射谱的强度比可以反过来检验计算的内部一致性。

五、实操中的参数约定与检查清单

下面给出一份可直接对照使用的映射检查表，用于在完成激发态计算后逐项核对，避免单位、简并度和展宽约定不一致导致的谱图偏差。

六、一个完整的数值示例

假设某分子 S₀→S₁ 跃迁计算结果为：ΔE = 0.12 Hartree（约 3.27 eV，26360 cm⁻¹），μ_if = 1.8 a.u.。则

f = (2/3) × 0.12 × (1.8)² = 0.259

取高斯线型 HWHM = 0.15 eV，则峰值吸收截面可由 f 与线型归一化条件估算。若实验测得该带积分 ε 约为 3×10⁴ L·mol⁻¹·cm⁻²，对应 f ≈ 0.27，与计算值 0.259 吻合良好，说明映射链自洽。若实验 f 仅为 0.03，则需检查是否漏掉了溶剂导致的强度降低，或计算中是否误将振子强度按简并度平均而非求和。

对 CD 谱，若同一跃迁的 m_if = 0.4 μ_B 且与 μ_if 平行，则 R ≈ 1.8 × 0.4 = 0.72 a.u.，对应 Δε 峰值约在 10–100 L·mol⁻¹·cm⁻¹ 量级，具体数值取决于展宽。R 的符号直接给出 CD 峰的正负，是判断绝对构型的关键输出。

总结来看，光谱预测的物理量映射不是单一公式，而是一条从波函数到实验谱的完整链条：跃迁偶极矩提供方向和强度，振子强度提供吸收权重，旋转强度提供手性信息，截面和线型提供可观测谱形。每一步都有明确的单位约定和物理假设，只有逐项核对，才能让计算结果真正对应到实验可观测的量。