We present a hybrid computational method for simulating the dynamics of macromolecules in solution which couples a mesoscale solver for the fluctuating hydrodynamics (FH) equations with molecular dynamics to describe the macromolecule. The two models interact through a dissipative Stokesian term first introduced by Ahlrichs and Dünweg [J. Chem. Phys. 111, 8225 (1999)]. We show that our method correctly captures the static and dynamical properties of polymer chains as predicted by the Zimm model. In particular, we show that the static conformations are best described when the ratio , where is the Lennard-Jones length parameter and is the monomer bond length. We also find that the decay of the Rouse modes’ autocorrelation function is better described with an analytical correction suggested by Ahlrichs and Dünweg. Our FH solver permits us to treat the fluid equation of state and transport parameters as direct simulation parameters. The expected independence of the chain dynamics on various choices of fluid equation of state and bulk viscosity is recovered, while excellent agreement is found for the temperature and shear viscosity dependence of center of mass diffusion between simulation results and predictions of the Zimm model. We find that Zimm model approximations start to fail when the Schmidt number . Finally, we investigate the importance of fluid fluctuations and show that using the preaveraged approximation for the hydrodynamic tensor leads to around 3% error in the diffusion coefficient for a polymer chain when the fluid discretization size is greater than .
Skip Nav Destination
Article navigation
21 April 2007
Research Article|
April 18 2007
Hybrid method coupling fluctuating hydrodynamics and molecular dynamics for the simulation of macromolecules
G. Giupponi;
G. Giupponi
a)
Centre for Computational Science, Department of Chemistry,
University College London
, 20 Gordon Street, London WC1H 0AJ, United Kingdom
Search for other works by this author on:
G. De Fabritiis;
G. De Fabritiis
b)
Computational Biochemistry and Biophysics Laboratory (GRIB-IMIM/UPF)
, Barcelona Biomedical Research Park (PRBB), C/Doctor Aiguader 88, 08003 Barcelona, Spain
Search for other works by this author on:
Peter V. Coveney
Peter V. Coveney
c)
Centre for Computational Science, Department of Chemistry,
University College London
, 20 Gordon Street, London WC1H 0AJ, United Kingdom
Search for other works by this author on:
J. Chem. Phys. 126, 154903 (2007)
Article history
Received:
February 06 2007
Accepted:
March 02 2007
Citation
G. Giupponi, G. De Fabritiis, Peter V. Coveney; Hybrid method coupling fluctuating hydrodynamics and molecular dynamics for the simulation of macromolecules. J. Chem. Phys. 21 April 2007; 126 (15): 154903. https://doi.org/10.1063/1.2720385
Download citation file:
Sign in
Don't already have an account? Register
Sign In
You could not be signed in. Please check your credentials and make sure you have an active account and try again.
Sign in via your Institution
Sign in via your InstitutionPay-Per-View Access
$40.00
Citing articles via
Related Content
Lattice-Boltzmann simulations of the dynamics of polymer solutions in periodic and confined geometries
J. Chem. Phys. (February 2005)
Combining molecular dynamics with Lattice Boltzmann: A hybrid method for the simulation of (charged) colloidal systems
J. Chem. Phys. (May 2005)
Universal scaling and characterization of gelation in associative polymer solutions
Journal of Rheology (July 2021)
The raspberry model for hydrodynamic interactions revisited. I. Periodic arrays of spheres and dumbbells
J. Chem. Phys. (August 2015)
The Raspberry model for hydrodynamic interactions revisited. II. The effect of confinement
J. Chem. Phys. (August 2015)