Linking hydrophobicity and hydrodynamics by the hybrid fluctuating hydrodynamics and molecular dynamics methodologies

The development of a hybrid fluctuating hydrodynamics (FHD) and molecular dynamics (MD) simulation method that combines the molecular dynamics of moving particles with the fluctuating hydrodynamics of solvent fields on Eulerian grid cells is presented. This method allows resolution of solute-solvent interfaces and realization of excluded volumes of particles in the presence of hydrodynamic coupling. With these capabilities, we show that the ubiquitous forces mediated by the solvent, hydrophobicity and hydrodynamics, can be linked in a mesoscopic simulation. The strategies we devise to overcome the numerical issues of mixing variables in the Eulerian and Lagrangian coordinate systems, i.e., using a pair of auxiliary fluids to realize the excluded volumes of particles and assigning collocating gridding systems on solutes to interface with solvent fields, are also presented. Simulation results show that the hybrid FHD and MD method can reproduce the solvation free energies and scaling laws of particles dynamics for hydrophobes of different sizes. The collapse of two hydrophobic particles was also simulated to illustrate that the hybrid FHD and MD method has the potential to be generally applied to study nanoscale self-assembly and dynamics-structure-function relationships of biomolecules.