溶剂效应的自由能差定义及其对反应动力学与电子光谱的调控

说明：本文华算科技主要介绍溶剂效应如何由状态自由能差定义，混合溶剂的局域组成怎样偏离体相，以及溶剂怎样改变反应速率和电子吸收光谱。

溶剂效应怎样由自由能差定义？

溶剂效应指溶剂种类或组成改变后，某个结构、能量、速率或光谱量随之发生的差别。讨论这个差别时，温度、压力、浓度标度、溶质构型和电子态都要保持可比。同一物理量在两种溶剂条件下的差值，才是可计算、可与实验对应的溶剂效应。

单个物种从参照相转入溶液的自由能记为 ΔG_solv。对反应 R→P，溶剂 S 中的反应自由能为 ΔG_r^S=G_P^S−G_R^S。把溶剂从 S₁ 换成 S₂ 后，ΔΔG_r=ΔG_r^S2−ΔG_r^S1，也等于产物和反应物各自转移自由能之差。溶剂对产物的稳定量很大，若对反应物的稳定量同样大，反应平衡仍可能几乎不动。

速率对应另一组状态。活化自由能 ΔG^‡=G_TS−G_R，溶剂引起的势垒变化写成 ΔΔG^‡=ΔG_solv(TS)−ΔG_solv(R)。决定加速或减速的是过渡态相对反应物的溶剂化位移，产物是否稳定不能代替这项判断。

图1中的黑线对应水相参照，红线和绿线表示两种混合溶剂对各状态施加不同稳定量。红色情形把反应物—质子复合物压得更低，过渡态相对升高；绿色情形对过渡态的稳定更强，势垒随之下降。曲线强调的是三种状态各自的自由能移动，不能把“某溶剂更极性”直接换算成统一的加速方向。

带电物种跨溶剂比较时，标准态和单离子自由能约定会直接进入数值。1 atm 理想气体与 1 mol L⁻¹ 理想溶液之间的标准态换算在 298 K 约为 7.9 kJ mol⁻¹。同一电荷约定、同一力场或电子结构方法、同一标准态必须贯穿所有状态，否则数值差里会混入基准变化。

图2中 H₃O⁺转入 1,4-二氧六环、四氢呋喃和 γ-戊内酯的自由能为正，转入 N-甲基吡咯烷酮、丙酮和二甲基亚砜时可降为负；Cl⁻在六种转移中均付出正自由能。酸解离平衡取决于阳离子与共轭碱的合计变化。氢键受体能力和离子特异相互作用在这里比单个体相介电常数更接近实际能量来源。

混合溶剂为什么会形成局域偏聚？

二元溶剂的体相摩尔分数只描述整盒液体的平均组成。溶质附近若存在氢键、色散、静电吸引或疏溶剂作用，两种溶剂进入第一、第二溶剂壳层的概率会分开。某组分在局部区域的比例高于体相，称为优先溶剂化；比例降低，则是优先排斥。局域组成 x_i(r) 可以明显偏离体相组成 x_i^bulk，偏离方向由近程相互作用和分子尺寸共同决定。

分子动力学轨迹可用径向分布函数 g_ij(r)、配位数、局域/体相分配系数和 Kirkwood–Buff 积分描述这种偏聚。g(r) 的峰给出某组分在特定距离出现的相对概率，第一极小值内的积分给出平均配位数。单帧构型只显示一次涨落，稳定的局域偏聚要从足够长的轨迹和多个独立初态统计，才能与偶然密度涨落区分。

图3把纤维素表面分成非极性和极性区域。橙色等值面表示 THF 局域浓度达到体相的 3 倍，主要贴附在非极性表面；极性表面附近则形成富水区域。小角中子散射在 Q≈0.2 Å⁻¹处的结构信号与 MD 计算散射曲线对应，给出的局域团簇尺度约为 0.5 nm。混合溶剂在宏观上均一，材料表面附近仍可出现纳米尺度分相，局域可接触溶剂面积也随之改变。

溶剂在不同表面区域的竞争会改变链内接触、链间聚集和材料构象，自由能中既有溶质—溶剂作用，也有溶剂分子被置换回体相时的熵变化。把二元溶剂的介电常数按体积分数平均，只能近似长程极化，无法给出某类分子在特定位点的富集和停留时间。

图4中木质素在水里更紧缩，在 THF 和 THF/水中呈较展开的线团构象。混合溶剂曲线显示 THF 在木质素近邻的相对密度高于水，且对不同极性原子区域的富集程度不同。同一个“THF/水=1:1”标签下，局域溶剂环境随表面官能团改变；计算报告应列出体相配比、局域统计范围和采样温度。

溶剂为什么会改变活化自由能？

过渡态理论给出 k≈κ(k_BT/h)exp(−ΔG^‡/RT)。其中 κ 是透射系数，ΔG^‡包含构型采样、溶剂重排和熵。298 K 下势垒相差约 5.7 kJ mol⁻¹，速率常数便可相差一个数量级。氢键、离子配对、溶剂笼和介电极化分别稳定反应物与过渡态，净速率响应取决于这些贡献在两种状态之间是否抵消。

断键前若需要分子扭转、质子换位或溶剂交换，溶剂组成便会改变进入反应区的构型权重。断键坐标之外，分子扭转、质子位置、溶剂配位数或离子间距都可能构成慢变量。只沿一条静态内坐标优化，得到的是特定溶剂排布下的势能路径；显式溶剂增强采样得到的自由能面则把大量微观构型压缩到选定集体变量上。

图5以断裂 C—C 键距离 d_C–C和二面角 φ 为坐标。四个环境保留相似的反应盆地，势垒高度和谷底宽度却发生变化；反应物盆地中的扭转转换远快于断键，可在断键前充分采样。溶剂改变的不只是一条曲线的最高点，还包括到达过渡区域之前的构象权重。

把二维自由能面沿 φ 积分后得到边缘自由能曲线，可比较真空、二甲苯、水和单体环境下的 ΔG^‡。该体系中溶剂使断键自由能垒相对真空降低约 4 kcal mol⁻¹，不同液体的速率仍受活化熵和指前因子区分。势垒、温度斜率和透射系数需要分别检查。

静态 DFT 加连续介质适合识别长程极化造成的能量移动，具体氢键、溶剂交换和构象熵则要求离散溶剂与采样。若两条路径只相差数 kJ mol⁻¹ 的活化自由能，溶剂盒尺寸、力场、泛函、集体变量和轨迹收敛都可能改变排序。此时应报告误差区间，并核对过渡区域是否出现未纳入坐标的慢溶剂运动。

溶剂为什么会移动吸收峰和电子态？

吸收峰对应基态到激发态的垂直能隙 hν=E_exc−E_gs。溶剂若更强地稳定激发态，能隙缩小，吸收向长波移动；若基态获得更强稳定，能隙增大，峰位向短波移动。红移或蓝移由两个电子态的差分溶剂化决定，分子偶极矩增大只给出可能方向，还需确认具体态的电荷分布和轨道性质。

电子跃迁发生得很快，溶剂分子的核坐标来不及在吸收瞬间完成重排。计算垂直激发时，电子极化与取向极化处在不同时间尺度：前者可以响应，后者保持在基态分布；发射则从已弛豫的激发态溶剂环境出发。静态介电常数、光学介电响应和折射率对应不同响应过程，混用会把非平衡溶剂化误写成普通能量修正。

图7中 2-硫代胞嘧啶在乙酸乙酯、乙腈、二甲基亚砜、乙醇、甲醇和水中的低能吸收带位置与形状均不同。质子溶剂还能通过 N、S 位点附近的氢键改变几何和电子态混合。峰位没有沿一个体相极性数值作简单单调移动，特定氢键构型与连续介质极化需要分开辨认。

分别计算几何弛豫、介电响应、折射率和显式氢键，可定位溶致变色的能量来源。对图7所示体系，连续介质中的能级变化在 ε 从 1 增至约 10 时已大幅完成；多个常用极性溶剂的 ε 都高于这个范围，继续增加 ε 对计算跃迁能的影响可能很小。分子级氢键却仍会让不同溶剂产生明显差别。

图8a沿真空到水相优化构型插值，图8b单独改变 ε，图8c采用各溶剂的几何、ε 与折射率，图8d再加入显式溶剂分子。只有最后一组才重现较强的溶剂特异能级分裂，说明该体系的主要差别源于微观氢键。图谱解释要追踪相同电子态的组态与振子强度；若态排序发生交换，只按“第一个激发态”编号比较会接错能级。

溶剂效应的计算结果应写明所比较的量、标准态、温度、溶剂组成、局域构型采样和电子态身份。反应自由能、活化自由能、局域溶剂比例与垂直激发能属于不同观察窗口；这些量可能一起变化，也可能由不同分子作用控制。结论只在这些条件覆盖的状态空间内成立。