[论文解读] Long-Range Order and Interactions of Macroscopic Objects in Polar Liquids

П. О. Федичев, L. I. Men’shikov|arXiv (Cornell University)|Jan 9, 2006

Advanced Chemical Sensor Technologies被引用 5

一句话总结

该论文提出了一种极性液体的连续矢量场模型，能够同时捕捉短程氢键作用与长程偶极-偶极相互作用，从而实现对水环境中宏观物体之间强作用力、长程作用力以及取向依赖性相互作用的描述。该模型解释了疏水界面处的自发极化现象、水合层中的铁电性有序结构，以及与氢键网络转变相关的拓扑激发。

ABSTRACT

We develop a phenomenological vector model of polar liquids capable to describe aqueous interactions of macroscopic bodies. It is shown that a strong, long-range and orientationally dependent interaction between macroscopic objects appears as a result of competition between short-range (hydrogen bonding) and the long-range dipole-dipole interactions of the solvent molecules. Spontaneous polarization of molecular dipoles next to a hydrophobic boundaries leads to formation of globally ordered network of hydrogen-bonded molecules with ferroelectric properties. The proposed vector model naturally describes topological excitations on the solute boundaries and can be used to explain the hydrogen bonds networks and order-disorder phase transitions in the hydration water layer.

研究动机与目标

• 开发一种现象学矢量模型，统一描述极性液体中的短程（氢键）与长程（偶极-偶极）相互作用。
• 解释水环境中宏观溶质之间强作用力、长程作用力及取向依赖性相互作用的起源。
• 描述疏水表面附近水合层中自发极化与铁电性有序结构的形成机制。
• 对溶质界面处氢键网络中的拓扑激发与相变进行建模。
• 为药物发现中的溶剂化自由能与pKa计算，提供一种计算效率更高的显式溶剂分子动力学替代方法。

提出的方法

• 基于连续矢量场 $\mathbf{s}(\mathbf{r}) = \langle \mathbf{d} \rangle / d_0 $ 的形式化方法，表示溶剂分子的局域平均极化。
• 将铁电体的朗道-金兹堡型理论适配用于建模同时具有标量密度与矢量极化自由度的极性液体的自由能。
• 通过幂律衰减 $ \langle s_\alpha s_\beta \rangle \sim r^{-3} $ 的方式引入长程偶极-偶极相互作用，该行为与高介电常数相关。
• 利用关联函数推导溶质之间的渐近相互作用势，揭示其取向依赖性。
• 线性化模型再现了关键分子动力学模拟结果：偶极子平行于疏水表面排列，以及溶质附近出现的涡旋状偶极结构。
• 将该模型应用于QUANTUM软件平台，实现pKa预测、质子化态分配以及药物设计中的结合亲和力计算。

实验结果

研究问题

• RQ1 竞争的短程氢键作用与长程偶极-偶极相互作用如何导致长程、取向依赖的溶剂化力？
• RQ2 水分子偶极子在疏水宏观表面附近自发极化的物理机制是什么？
• RQ3 氢键网络中的拓扑激发如何与水合层的相变相关联？
• RQ4 连续矢量模型能否在无需全原子分辨率的情况下，准确再现显式溶剂分子动力学模拟的关键特征？
• RQ5 介电屏蔽与长程关联在水环境中带电或中性溶质之间相互作用的介导中起何种作用？

主要发现

• 该模型预测极化关联函数呈现长程幂律衰减 $ \langle s_\alpha s_\beta \rangle \sim r^{-3} $，其与水的高静态介电常数直接相关。
• 在疏水边界附近，分子偶极子的自发极化导致水合层中形成全局有序、类似铁电性的网络结构。
• 由于短程氢键作用与长程偶极-偶极作用之间的竞争，宏观物体之间的相互作用变为长程、强作用且具有取向依赖性。
• 该模型自然地描述了溶质表面的拓扑激发，并提出其为分子动力学模拟中观察到的氢键网络渗透转变的物理起源。
• 该模型成功再现了分子动力学模拟中的关键现象，如偶极子平行于疏水表面排列，以及溶质周围出现的涡旋状偶极结构。
• 该矢量模型计算效率高，已在QUANTUM药物发现平台中实现，可准确预测溶剂化自由能、pKa值以及蛋白质-配体结合亲和力。