TY - JOUR

T1 - Implicit and explicit solvent models for the simulation of a single polymer chain in solution: Lattice Boltzmann versus Brownian dynamics

AU - Pham, Tri Thanh

AU - Schiller, Ulf D

AU - Jagadeeshan, Ravi Prakash

AU - Duenweg, Burkhard

PY - 2009

Y1 - 2009

N2 - We present a comparative study of two computer simulation methods to obtain static and dynamic properties of dilute polymer solutions. The first approach is a recently established hybrid algorithm based upon dissipative coupling between Molecular Dynamics and lattice Boltzmann (LB), while the second is standard Brownian Dynamics (BD) with fluctuating hydrodynamic interactions. Applying these methods to the same physical system (a single polymer chain in a good solvent in thermal equilibrium) allows us to draw a detailed and quantitative comparison in terms of both accuracy and efficiency. It is found that the static conformations of the LB model are distorted when the box length L is too small compared to the chain size. Furthermore, some dynamic properties of the LB model are subject to an L^ -1 finite size effect, while the BD model directly reproduces the asymptotic L \to \infty behavior. Apart from these finite size effects, it is also found that in order to obtain the correct dynamic properties for the LB simulations, it is crucial to properly thermalize all the kinetic modes. Only in this case, the results are in excellent agreement with each other, as expected. Moreover, Brownian Dynamics is found to be much more efficient than lattice Boltzmann as long as the degree of polymerization is not excessively large.

AB - We present a comparative study of two computer simulation methods to obtain static and dynamic properties of dilute polymer solutions. The first approach is a recently established hybrid algorithm based upon dissipative coupling between Molecular Dynamics and lattice Boltzmann (LB), while the second is standard Brownian Dynamics (BD) with fluctuating hydrodynamic interactions. Applying these methods to the same physical system (a single polymer chain in a good solvent in thermal equilibrium) allows us to draw a detailed and quantitative comparison in terms of both accuracy and efficiency. It is found that the static conformations of the LB model are distorted when the box length L is too small compared to the chain size. Furthermore, some dynamic properties of the LB model are subject to an L^ -1 finite size effect, while the BD model directly reproduces the asymptotic L \to \infty behavior. Apart from these finite size effects, it is also found that in order to obtain the correct dynamic properties for the LB simulations, it is crucial to properly thermalize all the kinetic modes. Only in this case, the results are in excellent agreement with each other, as expected. Moreover, Brownian Dynamics is found to be much more efficient than lattice Boltzmann as long as the degree of polymerization is not excessively large.

U2 - 10.1063/1.3251771

DO - 10.1063/1.3251771

M3 - Article

VL - 131

SP - 1

EP - 11

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 16

ER -