从GROMACS分子动力学模拟到电子吸收光谱：有机溶剂中小分子模拟教程

《Journal of Chemical Education》：From GROMACS Molecular Dynamics Simulations to Electronic Absorption Spectra: A Tutorial for Small Molecules in Organic Solvents

　　

经典分子动力学是一种用于研究非共价分子间相互作用和构象柔性的理论方法，能够模拟聚集/溶解和吸附等过程。分子动力学和量子化学正日益融入化学课程体系，然而学生在将理论概念转化为可靠的模拟方案时可能遇到困难，尤其是在生物分子系统之外的情境中。本文呈现一个教程，旨在引导学生完成非水溶剂中小有机分子的模拟制备、执行与分析，随后利用时间依赖密度泛函理论（time-dependent density functional theory, TD-DFT）计算紫外-可见（ultraviolet–visible, UV–vis）吸收光谱。本教程聚焦若干实际问题：力场参数的获取、模拟单元的构建、系综的选择以及关键模拟参数的解释。自定义溶剂模型通过与参考密度和径向分布函数（radial distribution function, RDF）的对照进行验证，展示了如何利用所选可观测量评估模拟的可靠性。研究人员通过展示分子柔性和溶剂相互作用展宽了靛蓝（indigo）在氯仿中的UV–vis吸收带，说明了分子动力学的实用性。本贡献旨在作为高年级本科生和研究生物理化学及计算化学课程的手把手教育资源，也适用于化学、物理学或材料科学领域初涉分子模拟的博士生。为促进传播，研究人员已在GitHub上发布了一个开箱即用的软件包，包含溶剂盒以及将新分子纳入GROMACS模拟所需的所有文件和数据：https://github.com/cippee/SolventSafari-for-GROMACS。

**论文解读：从GROMACS分子动力学模拟到电子吸收光谱——面向有机溶剂中小分子的教学教程**

**一、研究背景与目的**

原子尺度模拟正越来越多地融入物理化学、计算化学和材料科学的研究生课程，用以辅助解释核磁共振和紫外-可见（UV–vis）光谱等实验数据。分子动力学（MD）能够让学生探索构象柔性、溶剂化作用及集体分子行为，为研究小分子、大分子和材料提供了广阔的工具。尽管MD在生物分子系统（如蛋白质、脂质和碳水化合物）中已有成熟的力场、教程和工作流程，但当学生面对非生物分子范式（即非氨基酸基分子溶解于非水有机溶剂，这是材料科学和物理化学中常见场景）时，实用的教学指导相当稀缺。学生难以将理论概念转化为稳健的模拟方案，从零开始构建模拟输入文件（尤其是力场参数和拓扑文件）往往繁琐且令初学者望而生畏。为填补这一空白，研究人员在《Journal of Chemical Education》上发表了本教程，旨在提供一个通用的、可迁移的模拟协议，引导学生完成从MD模拟到UV–vis光谱计算的全流程。

**二、研究方法概述**

本教程采用Q-Force基于量子力学DFT计算（PBEPBE/6–31+G(d)水平）自动推导OPLS-aa力场参数，利用GROMACS执行MD模拟，结合时间依赖密度泛函理论（TD-DFT；B3LYP/6–31G(d,p)水平）计算UV–vis吸收光谱，并使用VMD进行轨迹可视化与分析。研究首先为15种常用有机溶剂（包括丙酮、乙腈、四氯乙烯、四氯化碳、氯仿、环己烷、氯苯、二硫化碳、二氯甲烷、二甲基亚砜、乙醇、甲醇、正己烷、苯、四氢呋喃和甲苯）构建了含9000个原子的立方模拟盒子，密度依据实验值确定，并依次进行能量最小化、NVT和NPT平衡及5 ns的NPT产出模拟。随后以靛蓝溶于氯仿为案例，提取20个包含第一溶剂化壳层的微观态，对每个微观态进行TD-DFT计算并等权重叠加获得整体光谱。所有输入输出文件已公开于GitHub仓库（SolventSafari-for-GROMACS）。

**三、研究结果**

**溶剂模型的验证：密度与径向分布函数。** 研究人员通过对比计算密度与实验数据，验证了力场参数的可靠性。5 ns产出模拟得到的密度平均值与实验值呈强相关，几乎所有溶剂的相对误差低于10%，部分溶剂（丙酮、四氯甲烷、四氯乙烯、苯）误差低于2%。回归线斜率接近1，截距仅26 kg/m³，决定系数R²接近1，表明所选力场参数和模拟条件能很好地预测不同分子液体的密度。此外，密度预测的绝对相对误差与溶剂沸点存在关联：低沸点溶剂的误差较大，这可能源于可压缩性数据的实验不确定性、液-气相变附近更强的密度涨落以及模拟该类条件的困难。溶剂模型的径向分布函数（g(r)）也与文献数据吻合，验证了溶剂微观结构的正确性。

**力场参数的质量检验。** 研究以乙醇为例，对MD轨迹进行短时密集采样（0.05–0.1 fs步长、数ps时长）并使用VMD计算红外光谱密度，与DFT计算的IR光谱对比，验证了Q-Force自定义OPLS-aa参数（尤其键合参数和力常数）的合理性。结果显示，Q-Force修正后的OPLS-AA参数（红色曲线）较原生OPLS-AA参数（蓝色曲线）与量子力学光谱（绿色曲线）的一致性更好，确认了该参数化流程的有效性。

**靛蓝在氯仿中UV–vis吸收光谱的模拟。** 研究人员将单个靛蓝分子溶于氯仿溶剂盒，经200 ps NVT和200 ps NPT平衡后进行2 ns NPT产出模拟（T = 300 K，P = 1 bar）。从MD轨迹中每100 ps提取一个微观态，共获得20个包含第一溶剂化壳层（约7 ?）的靛蓝微观态。对每个微观态进行B3LYP/6–31G(d,p)水平TD-DFT计算，等权重叠加各微观态光谱以获得整体吸收谱。结果表明，仅当所有微观态贡献叠加后（图10中图），模拟光谱与实验光谱（图10上图）在峰位和峰形上均高度吻合：560 nm附近的肩峰和主吸收峰的总体带宽均得到再现。相反，任一单个微观态的光谱都无法准确捕获靛蓝的吸收特征（图10下图），这凸显了MD采样在捕捉分子柔性和溶剂涨落对光谱展宽贡献中的关键作用。

**四、结论总结**

本教程带领学生完成了一个从分子模型构建到UV–vis光谱计算的全流程模拟工作，强调了方法论选择背后的物理依据，并提供了必要的软件命令。通过对多种有机溶剂的系统构建与验证，学生能够理解力场参数如何转化为密度和径向分布函数等可观测性质；靛蓝/氯仿案例进一步展示MD与TD-DFT联合使用解释UV–vis实验数据的可行性，凸显溶剂化效应在光谱成形中的关键角色。教程中所有文件、示例和溶剂盒均已公开，旨在降低分子模拟进入课堂的门槛。该协议可合理推广至新分子体系，希望这一资源能够促进计算化学和物理化学学生更广泛地采用MD模拟。