[论文解读] Simulating Van der Waals-interactions in water/hydrocarbon-based complex fluids

Igor Pasichnyk, Ralf Everaers+1 more|ArXiv.org|Nov 8, 2007

Material Dynamics and Properties参考文献 19被引用 5

一句话总结

本文表明，一种用于异质介电介质中静电作用的局部蒙特卡洛算法，其内在地捕捉了由Lifshitz理论预测的零频（Keesom）相互作用贡献，包括多体效应和盐的屏蔽效应。该方法准确再现了水/烃类复杂流体中的非加和性、长程相互作用，其结果与介电常数对比度高达78的三明治结构的解析结果一致。

ABSTRACT

In systems composed of water and hydrocarbons Van der Waals-interactions are dominated by the non-retarded, classical (Keesom) part of the Lifshitz-interaction; the interaction is screened by salt and extends over mesoscopic distances of the order of the size of the (micellar) constituents of complex fluids. We show that these interactions are included intrinsically in a recently introduced local Monte Carlo algorithm for simulating electrostatic interactions between charges in the presence of non-homogeneous dielectric media.

研究动机与目标

• 为解决标准隐式溶剂模型在模拟复杂流体中范德华相互作用时的局限性。
• 包含Lifshitz相互作用中占主导地位的非色散、经典（Keesom）部分，该部分在水/烃类体系中占主导。
• 考虑盐对范德华相互作用中经典（Keesom）部分的屏蔽效应，以及胶束或层状结构中的长程、介观尺度扩展的相互作用。
• 验证一种局部蒙特卡洛算法作为在介电常数空间变化的任意几何结构中模拟这些相互作用的物理上一致的方法。

提出的方法

• 采用一种新近发展的局部蒙特卡洛算法，通过局部电荷-场相互作用模拟非均匀介电介质中的静电作用。
• 使用热力学积分法，引入耦合参数λ，计算参考与目标介电构型之间的自由能差。
• 采用蠕虫算法（worm algorithm）采样涨落的电位移场D，从而高效采样受限的配分函数。
• 将该方法应用于介电常数为ε₁ = 2和ε₁ = 78的三明治结构系统，中间夹有低介电常数（ε₂ = 1）间隙。
• 在20个λ点上进行数值积分，每个λ点进行800次平衡步和800次生产步，共使用500组独立的自由能估计。
• 通过外推至h → ∞，利用Lifshitz相互作用的解析形式（公式13）校正自能，从而分离出相互作用能。

实验结果

研究问题

• RQ1 能否通过异质介电介质中静电作用的局部蒙特卡洛算法，隐式再现Lifshitz理论预测的多体、长程范德华相互作用？
• RQ2 该方法在高介电常数对比度体系（如水和烃类）中，对范德华力的非加和性、非成对性质的捕捉效果如何？
• RQ3 该算法在经典（Keesom）贡献中对盐屏蔽效应的考虑程度如何？
• RQ4 该方法能否可靠地应用于具有任意空间变化介电性质的复杂流体几何结构？

主要发现

• 该局部蒙特卡洛算法成功捕捉了Lifshitz相互作用的零频（Keesom）部分，包括所有多体效应，且无需显式引入多体项。
• 对于ε₁ = 2和ε₁ = 78的模拟结果与解析Lifshitz理论高度一致，尤其在Hamaker近似失效的非成对区域表现优异。
• 晶格离散化误差较小，表现为介电常数对比度为2和78时，a = 1.0与a = 0.5的计算结果一致。
• 该方法正确再现了高介电常数对比体系中相互作用的长程、介观尺度扩展特性。
• 经外推至h → ∞校正后的自由能差计算，成功分离出相互作用能，尽管存在因数值近似引起的微小偏差，但仍与解析曲线高度吻合。
• 该方法具有鲁棒性和可扩展性，适用于任意几何结构和空间变化的介电函数，可实现对胶束、层状结构等复杂流体结构的模拟。