基于隐式溶剂化的分子高压模拟增强方法研究

《The Journal of Physical Chemistry A》：Enhancing Molecular High-Pressure Simulations by Implicit Solvation

　　

研究人员针对量子化学程序包 Q-Chem 中高压模拟方法 GOSTSHYP（Gaussians On Surface Tesserae Simulate HYdrostatic Pressure）与 X-HCFF（eXtended Hydrostatic Compression Force Field）此前无法稳定结合隐式溶剂模型 C-PCM（Conductor-like Polarizable Continuum Model）的问题，开发了独立的分子表面构建与离散化程序。该程序避免了原有共用表面算法导致的数据覆盖与计算不稳定，实现了 GOSTSHYP 或 X-HCFF 与 C-PCM 的稳定耦合，从而在单分子高压模拟中引入化学环境效应。研究人员选取可同时存在中性态与两性离子态的化合物（如甘氨酸、氨基磺酸、2-二甲基氨基苯甲酸）及原硅酸二聚反应体系，结合密度泛函理论（DFT）与 B3LYP-D3(BJ)/aug-cc-pVDZ 级别的计算，分析了压力下的能量差异、分子体积变化及拉曼光谱位移。结果表明，C-PCM 对两性离子结构的稳定性至关重要，尤其在高压条件下显著改变能量趋势与产物构型。与原硅酸二聚反应实验数据相比，耦合方法所得自由能变化与实验观测的加压导致解聚现象一致，验证了该方法在地球化学与极端条件物理化学研究中的适用性。

该研究发表于《The Journal of Physical Chemistry A》，旨在解决高压下单分子量子化学模拟缺乏可靠溶剂环境处理的问题。此前的高压模拟方法如 GOSTSHYP 与 X-HCFF 虽能施加静水压力，但无法与隐式溶剂模型稳定结合，限制了其在溶液相及高压化学环境中的应用。研究人员通过构建独立的分子表面离散化程序，实现了与 C-PCM 的无缝耦合，使高压模拟能够考虑溶剂极化效应，拓展了可研究的分子体系范围。

关键技术方法包括：在 Q-Chem 6.2 开发版中实现独立表面类，支持多种表面构建方式（如范德华表面、溶剂可及表面）；采用 Lebedev 网格离散化；结合 DFT/B3LYP-D3(BJ) 与不同基组（aug-cc-pVDZ、6–311++G(d,p)）进行结构与振动频率优化；使用 GOSTSHYP 与 X-HCFF 施加 0–1000 GPa 压力；通过 C-PCM 设置不同介电常数模拟溶剂环境；计算分子体积（PV 法）、拉曼位移及反应自由能变化。

研究结果分为四部分。第一部分代码验证表明，新表面程序在零压与高压条件下均能复现原有结果，并在与 C-PCM 联用时保持数值稳定，消除了旧版本中因临时文件覆盖导致的崩溃问题。第二部分两性离子结构研究表明，C-PCM 是获得并稳定甘氨酸、氨基磺酸及 2-二甲基氨基苯甲酸两性离子态的必要条件；在 GOSTSHYP 加压下，两性离子态因体积更小而能量更稳定，且在约 300 GPa 时甘氨酸转变为环状结构；拉曼光谱模拟结果与实验数据在主要振动模式的压力依赖趋势上吻合良好。第三部分氨基磺酸研究显示，C-PCM 在不同压力下均稳定两性离子态，并改变中性态与两性离子态的能量差趋势；无溶剂模型时，高压可导致中性态稳定或出现构象扭转。第四部分原硅酸二聚反应研究表明，仅用 GOSTSHYP 时高压促进二聚，而结合 C-PCM 后高压反而抑制二聚，产物形成水合配合物，与实验观测的解聚趋势一致。

讨论与结论部分指出，独立表面程序不仅解决了 GOSTSHYP 与 X-HCFF 和 C-PCM 的稳定耦合问题，还提高了高压模拟的物理准确性，尤其在处理电荷分离体系及溶剂化反应时优势明显。该方法为地球深部及外星极端条件下的生物分子和矿物化学反应模拟提供了可行工具，未来可通过加入显式溶剂层及激发态计算进一步提升精度。