Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTowards Better Integrators for Dissipative Particle Dynamics Simulations
91
0
0.0
(
0
)
Authors
Gerhard Besold
Ask ChatGPT about the research
No Arabic abstract
Coarse-grained models that preserve hydrodynamics provide a natural approach to study collective properties of soft-matter systems. Here, we demonstrate that commonly used integration schemes in dissipative particle dynamics give rise to pronounced artifacts in physical quantities such as the compressibility and the diffusion coefficient. We assess the quality of these integration schemes, including variants based on a recently suggested self-consistent approach, and examine their relative performance. Implications of integrator-induced effects are discussed.
rate research
Read More
In this work we compare and characterize the behavior of Langevin and Dissipative Particle Dynamics (DPD) thermostats in a broad range of non-equilibrium simulations of polymeric systems. Polymer brushes in relative sliding motion, polymeric liquids in Poiseuille and Couette flows, and brush-melt interfaces are used as model systems to analyze the efficiency and limitations of different Langevin and DPD thermostat implementations. Widely used coarse-grained bead-spring models under good and poor solvent conditions are employed to assess the effects of the thermostats. We considered equilibrium, transient, and steady state examples for testing the ability of the thermostats to maintain constant temperature and to reproduce the underlying physical phenomena in non-equilibrium situations. The common practice of switching-off the Langevin thermostat in the flow direction is also critically revisited. The efficiency of different weight functions for the DPD thermostat is quantitatively analyzed as a function of the solvent quality and the non-equilibrium situation.
We show how to improve the molecular dynamics step of Hybrid Monte Carlo, both by tuning the integrator using Poisson brackets measurements and by the use of force gradient integrators. We present results for moderate lattice sizes.
We extend our previous study [J. Chem. Phys. 138, 204907 (2013)] to quantify the screening properties of four mesoscale smoothed charge models used in dissipative particle dynamics. Using a combination of the hypernetted chain integral equation closure and the random phase approximation, we identify regions where the models exhibit a real-valued screening length, and the extent to which this agrees with the Debye length in the physical system. We find that the second moment of the smoothed charge distribution is a good predictor of this behaviour. We are thus able to recommend a consistent set of parameters for the models.
In this paper we study the shape characteristics of star-like polymers in various solvent quality using a mesoscopic level of modeling. The dissipative particle dynamics simulations are performed for the homogeneous and four different heterogeneous star polymers with the same molecular weight. We analyse the gyration radius and asphericity at the bad, good and $theta$-solvent regimes. Detailed explanation based on interplay between enthalpic and entropic contributions to the free energy and analyses on of the asphericity of individual branches are provided to explain the increase of the apsphericity in $theta$-solvent regime.
The non-Newtonian behavior of a monodisperse concentrated dispersion of spherical particles was investigated using a direct numerical simulation method, that takes into account hydrodynamic interactions and thermal fluctuations accurately. Simulations were performed under steady shear flow with periodic boundary conditions in the three directions. The apparent shear viscosity of the dispersions was calculated at volume fractions ranging from 0.31 to 0.56. Shear-thinning behavior was clearly observed at high volume fractions. The low- and high-limiting viscosities were then estimated from the apparent viscosity by fitting these data into a semi-empirical formula. Furthermore, the short-time motions were examined for Brownian particles fluctuating in concentrated dispersions, for which the fluid inertia plays an important role. The mean square displacement was monitored in the vorticity direction at several different Peclet numbers and volume fractions so that the particle diffusion coefficient is determined from the long-time behavior of the mean square displacement. Finally, the relationship between the non-Newtonian viscosity of the dispersions and the structural relaxation of the dispersed Brownian particles is examined.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
