Department of Chemical and Biomolecular Engineering, University of California Berkeley, Berkeley, California 94720, USA.
J Chem Phys. 2012 Jul 28;137(4):044117. doi: 10.1063/1.4738763.
This work illustrates that fluctuating hydrodynamics (FHD) simulations can be used to capture the thermodynamic and hydrodynamic responses of molecular fluids at the nanoscale, including those associated with energy and heat transfer. Using all-atom molecular dynamics (MD) trajectories as the reference data, the atomistic coordinates of each snapshot are mapped onto mass, momentum, and energy density fields on Eulerian grids to generate a corresponding field trajectory. The molecular length-scale associated with finite molecule size is explicitly imposed during this coarse-graining by requiring that the variances of density fields scale inversely with the grid volume. From the fluctuations of field variables, the response functions and transport coefficients encoded in the all-atom MD trajectory are computed. By using the extracted fluid properties in FHD simulations, we show that the fluctuations and relaxation of hydrodynamic fields quantitatively match with those observed in the reference all-atom MD trajectory, hence establishing compatibility between the atomistic and field representations. We also show that inclusion of energy transfer in the FHD equations can more accurately capture the thermodynamic and hydrodynamic responses of molecular fluids. The results indicate that the proposed MD-to-FHD mapping with explicit consideration of finite molecule size provides a robust framework for coarse-graining the solution phase of complex molecular systems.
这项工作表明,波动流体动力学 (FHD) 模拟可用于捕捉纳米尺度分子流体的热力学和流体动力学响应,包括与能量和热传递相关的响应。使用全原子分子动力学 (MD) 轨迹作为参考数据,将每个快照的原子坐标映射到欧拉网格上的质量、动量和能量密度场,以生成相应的场轨迹。在这个粗粒化过程中,通过要求密度场的方差与网格体积成反比,明确施加与有限分子大小相关的分子长度尺度。通过场变量的波动,从全原子 MD 轨迹中计算出编码的响应函数和输运系数。通过在 FHD 模拟中使用提取的流体性质,我们表明,流体场的波动和弛豫与在参考全原子 MD 轨迹中观察到的情况定量匹配,从而在原子和场表示之间建立兼容性。我们还表明,在 FHD 方程中包含能量传递可以更准确地捕捉分子流体的热力学和流体动力学响应。结果表明,提出的具有明确考虑有限分子大小的 MD 到 FHD 映射为粗粒化复杂分子系统的溶液相提供了一个稳健的框架。