Coil-to-globule transition by dissipative particle dynamics simulation

The dynamics of a collapsing polymer under a temperature quench in dilute solution is investigated by dissipative particles dynamics. Hydrodynamic interactions and many-body interaction are preserved naturally by incorporating explicit solvent particles in this approach. Our simulation suggests a four-stage collapse pathway: localized clusters formation, cluster coarsening in situ, coarsening involving global backbone conformation change into a crumpled globule, and compaction of the globule. For all the quench depths and chain lengths used in our study, collapse proceeds without the chain getting trapped in a metastable "sausage" configuration, as reported in some earlier studies. We obtain the time scales for each of the first three stages, as well as its scaling with the quench depths ξ and chain lengths N. The total collapse time scales as τ(c) ～ ξ(-0.46 ± 0.04)N(0.98 ± 0.09), with the quench depth and degree of polymerization.     

Budding dynamics of individual domains in multicomponent membranes simulated by N-varied dissipative particle dynamics

We study the budding dynamics of individual domains in flat, multicomponent membranes using dissipative particle dynamics (DPD) simulations with varied bead number N, in which addition and deletion of beads based on their density at the membrane boundary is introduced. The budding process of a tubular bud, accompanied by a dynamical transition reflected in the energy and morphology evolutions, is investigated. The simulations show that budding duration is shortened with increasing line tension and depends on the domain size quadratically. At low line tension, increasing bending modulus accelerates budding at first, but suppresses the process as it increases further. In addition, the controlling role of the surface tension in the budding process is also explored. Finally, we use the N-varied DPD to simulate the experimentally observed multicomponent tubular vesicles, and the three bud growth modes are confirmed.     

Schmidt number effects in dissipative particle dynamics simulation of polymers

Simulation studies for dilute polymeric systems are presented using the dissipative particle dynamics method. By employing two different thermostats, the velocity-Verlet and Lowe's scheme, we show that the Schmidt number (S(c)) of the solvent strongly affects nonequilibrium polymeric quantities. The fractional extension of wormlike chains subjected to steady shear is obtained as a function of S(c). Poiseuille flow in microchannels for fixed polymer concentration and varying number of repeated units within a chain is simulated. The nonuniform concentration profiles and their dependence on S(c) are computed. We show the effect of the bounce-forward wall boundary condition on the depletion layer thickness. A power law fit of the velocity profile in stratified Poiseuille flow in a microchannel yields wall viscosities different from bulk values derived from uniform, steady plane Couette flow. The form of the velocity profiles indicates that the slip flow model is not useful for the conditions of these calculations.     

Flow-induced translocation of polymers through a fluidic channel: a dissipative particle dynamics simulation study

The dynamics of flow-induced translocation of polymers through a fluidic channel has been studied by dissipative particle dynamics (DPD) approach. Unlike implicit solvent models, the many-body energetic and hydrodynamic interactions are preserved naturally by incorporating explicit solvent particles in this approach. The no-slip wall boundary and the adaptive boundary conditions have been implemented in the modified DPD approach to model the hydrodynamic flow within a specific wall structure of fluidic channel and control the particles' density fluctuations. The results show that the average translocation time versus polymer chain length satisfies a power-law scaling of τ ～N(1.152). The conformational changes and translocation dynamics of polymers through the fluidic channel have also been investigated in our simulations, and two different translocation processes, i.e., the single-file and double-folded translocation events, have been observed in detail. These findings may be helpful in understanding the conformational and dynamic behaviors of such polymer and/or DNA molecules during the translocation processes.     

Mesoscopic simulation of self-assembly in surfactant oligomers by dissipative particle dynamics

Self-assembling properties of surfactant oligomers in an aqueous medium is simulated by dissipative particle dynamics (DPD). The critical micellar concentration (CMC) of dimeric (oligomerization = 2) and trimeric (oligomerization = 3) surfactant is much lower than their single-chain counterpart. All surfactants form spherical micelles at the concentration not far above their CMC. The transition from spherical to cylindrical micelles exhibits with increasing surfactant concentration. Lamellar micelles will appear with further increasing the surfactant concentration. For dimeric and trimeric surfactants, cylindrical micelles transform into extremely long “wormlike” or “threadlike” micelles before the transition to lamellar micelles. These results are in qualitative agreement with laboratory experiment. Average aggregation numbers (AN) of micelles increase with a power law of AN  c when the surfactant concentration c CMC. The self-diffusion coefficients will drop with a power law of D  c⁻ when wormlike micelles are formed.     

Finite-size effects in dissipative particle dynamics simulations

We have performed dissipative particle dynamics (DPD) simulations to evaluate the effect that finite size of transversal area has on stress anisotropy and interfacial tension. The simulations were carried out in one phase and two phases in parallelepiped cells. In one-phase simulations there is no finite-size effect on stress anisotropy when the simulation is performed using repulsive forces. However, an oscillatory function of stress anisotropy is found for attractive-repulsive interactions. In the case of liquid-liquid interfaces with repulsive interaction between molecules, there is only a small effect of surface area on interfacial tension when the simulations are performed using the Monte Carlo method at constant temperature and normal pressure. An important but artificial finite-size effect of interfacial area on surface tension is found in simulations in the canonical ensemble. Reliable results of interfacial tension from DPD simulations can be obtained using small systems, less than 2000 particles, when they interact exclusively with repulsive forces.     

Constant-pressure simulations with dissipative particle dynamics

Dissipative particle dynamics (DPD) is a mesoscopic simulation method for studying hydrodynamic behavior of complex fluids. Ideally, a mesoscopic model should correctly represent the thermodynamic and hydrodynamic properties of a real system beyond certain length and time scales. Traditionally defined DPD quite successfully mimics hydrodynamics but is not flexible enough to accurately describe the thermodynamics of a real system. The so-called multibody DPD (MDPD) is a pragmatic extension of the classical DPD that allows one to prescribe the thermodynamic behavior of a system with only a small performance impact. In an earlier paper [S. Y. Trofimov, E. L. F. Nies, and M. A. J. Michels, J. Chem. Phys. 117, 9383 (2002)] we much improved the accuracy of the MDPD model for strongly nonideal systems, which are of most practical interest. The ability to correctly reproduce the equation of state of realistic systems in turn makes simulations at constant pressure sensible and useful. This situation of constant-pressure conditions is very common in experimental studies of (soft) condensed matter but has so far remained unexplored with the traditional DPD. Here, as a proof of concept, we integrate a modified version of the Andersen barostat into our improved MDPD model and make an evaluation of the performance of the new model on a set of single- and multicomponent systems. The modification of the barostat suppresses the "unphysical" volume oscillations after a sudden pressure change and simplifies the equilibration of the system.     

New parametrization method for dissipative particle dynamics

We introduce an improved method of parametrizing the Groot-Warren version of dissipative particle dynamics (DPD) by exploiting a correspondence between DPD and Scatchard-Hildebrand regular solution theory. The new parametrization scheme widens the realm of applicability of DPD by first removing the restriction of equal repulsive interactions between like beads, and second, by relating all conservative interactions between beads directly to cohesive energy densities. We establish the correspondence by deriving an expression for the Helmoltz free energy of mixing, obtaining a heat of mixing which is exactly the same form as that for a regular mixture (quadratic in the volume fraction) and an entropy of mixing which reduces to the ideal entropy of mixing for equal molar volumes. We equate the conservative interaction parameters in the DPD force law to the cohesive energy densities of the pure fluids, providing an alternative method of calculating the self-interaction parameters as well as a route to the cross interaction parameter. We validate the new parametrization by modeling the binary system SnI(4)SiCl(4), which displays liquid-liquid coexistence below an upper critical solution temperature around 140 degrees C. A series of DPD simulations were conducted at a set of temperatures ranging from 0 degrees C to above the experimental upper critical solution temperature using conservative parameters based on extrapolated experimental data. These simulations can be regarded as being equivalent to a quench from a high temperature to a lower one at constant volume. Our simulations recover the expected phase behavior ranging from solid-liquid coexistence to liquid-liquid coexistence and eventually leading to a homogeneous single phase system. The results yield a binodal curve in close agreement with the one predicted using regular solution theory, but, significantly, in closer agreement with actual solubility measurements.     

Many-body dissipative particle dynamics simulation of liquid/vapor and liquid/solid interactions

The combination of short-range repulsive and long-range attractive forces in many-body dissipative particle dynamics (MDPD) is examined at a vapor/liquid and liquid/solid interface. Based on the radial distribution of the virial pressure in a drop at equilibrium, a systematic study is carried out to characterize the sensitivity of the surface tension coefficient with respect to the inter-particle interaction parameters. For the first time, the approximately cubic dependence of the surface tension coefficient on the bulk density of the fluid is evidenced. In capillary flow, MDPD solutions are shown to satisfy the condition on the wavelength of an axial disturbance leading to the pinch-off of a cylindrical liquid thread; correctly, no pinch-off occurs below the cutoff wavelength. Moreover, in an example that illustrates the cascade of fluid dynamics behaviors from potential to inertial-viscous to stochastic flow, the dynamics of the jet radius is consistent with the power law predictions of asymptotic analysis. To model interaction with a solid wall, MDPD is augmented by a set of bell-shaped weight functions; hydrophilic and hydrophobic behaviors, including the occurrence of slip in the latter, are reproduced using a modification in the weight function that avoids particle clustering. The dynamics of droplets entering an inverted Y-shaped fracture junction is shown to be correctly captured in simulations parametrized by the Bond number, confirming the flexibility of MDPD in modeling interface-dominated flows.     

Synthesis and dissipative particle dynamics simulation of cross-linkable fluorinated diblock copolymers: self-assembly aggregation behavior in different solvents

Developing microstructures, such as low molecular aggregates, spherical micelles and multi-compartment micelles, is an expanding area of research in Materials Science. By applying an atom transfer radical polymerization (ATRP) process to cross-linkable fluorinated diblock copolymers and analyzing the data we are able to demonstrate the potential for developing films with different micro-structures for additional biological research. Applying the Dissipative Particle Dynamic (DPD) Method, Transmission Electron Microscopy (TEM) and Scanning Electron Microscopy (SEM) techniques to cross-linkable fluorinated diblock copolymers of (methyl methacrylate-co-hydroxyethyl methacrylate-co-butyl methacrylate)-b-2-(perfluoroalkyl)ethyl methacrylate (MMA-co-HEMA-co-BMA-b-FMA) we were able to analyze the structures and their relationships to the aggregation of various microstructure formations through the use of various solvents in the process. For the self-assembly of the cross-linkable diblock copolymer in solutions, the DPD simulation results are only in qualitative agreement with experimental data of aggregate morphologies and sizes. This suggests an improved approach to creating materials and methods necessary for studying microstructures in films used in other research areas. Our work examines whether using selective solvents can be easily extended to prepare aggregates with different morphologies, which is an effective shortcut to obtain films with different microstructures. DPD simulation can be considered as an adjunct to experiments and provides other valuable information for the experiment.     

Configurational temperature in membrane simulations using dissipative particle dynamics

The use of excessively long time steps in dissipative particle dynamics simulations may produce simulation artifacts due to the generation of configurations which are not representative of the desired canonical ensemble. The configurational temperature, among other quantities, may be used to assess the extent of the deviation from equilibrium. This paper presents results for simulations of models of water and lipid bilayer membranes to illustrate the nature of the problems.     

Constant-pressure and constant-surface tension simulations in dissipative particle dynamics

We present a method for constant-pressure and constant-surface tension simulations in dissipative particle dynamics using a Langevin piston approach. We demonstrate that the corresponding equations of motion lead to the relevant ensembles and propose an appropriate scheme of integration. After having identified a suitable set of parameters for the approach, we demonstrate the feasibility of the approach by applying it to two different systems, a simple isotropic fluid and an anisotropic fluid lipid-bilayer membrane in water. Results are presented for, respectively, isothermal bulk compressibility, tracer diffusion coefficient, lipid head-group area, and isothermal area compressibility. We find that our Langevin piston approach leads to improvements over other approaches in terms of faster equilibration and shorter correlation times of various system variables.     

Comparison between theoretical values and simulation results of viscosity for the dissipative particle dynamics method

In order to investigate the validity of the dissipative particle dynamics method, which is a mesoscopic simulation technique, we have derived an expression for viscosity from the equation of motion of dissipative particles. In the concrete, we have shown the Fokker-Planck equation in phase space, and macroscopic conservation equations such as the equation of continuity and the equation of momentum conservation. The basic equations of the single-particle and pair distribution functions have been derived using the Fokker-Planck equation. The solutions of these distribution functions have approximately been solved by the perturbation method under the assumption of molecular chaos. The expressions of the viscosity due to momentum and dissipative forces have been obtained using the approximate solutions of the distribution functions. Also, we have conducted nonequilibrium dynamics simulations to investigate the influence of the parameters, which have appeared in defining the equation of motion in the dissipative particle dynamics method. The theoretical values of the viscosity due to dissipative forces in the Hoogerbrugge-Koelman theory are in good agreement with the simulation results obtained by the nonequilibrium dynamics method, except in the range of small number densities. There are restriction conditions for taking appropriate values of the number density, number of particles, time interval, shear rate, etc., to obtain physically reasonable results by means of dissipative particle dynamics simulations.     

Improved dissipative particle dynamics simulations of lipid bilayers

The authors introduce a new parameterization for the dissipative particle dynamics simulations of lipid bilayers. In this parameterization, the conservative pairwise forces between beads of the same type in two different hydrophobic chains are chosen to be less repulsive than the water-water interaction, but the intrachain bead interactions are the same as the water-water interaction. For a certain range of parameters, the new bilayer can only be stretched up to 30% before it ruptures. Membrane tension, density profiles, and the in-plane lipid diffusion coefficient of the new bilayer are discussed in detail. They find two kinds of finite size effects that influence the membrane tension: lateral finite size effects, for which larger membranes rupture at a smaller stretch, and transverse finite size effects, for which tensionless bilayers are more compact in larger systems. These finite size effects become rather small when the simulation box is sufficiently large.     

Examination of membrane fusion by dissipative particle dynamics simulation and comparison with continuum elastic models

A dissipative particle dynamics model is applied to probe the lipidic membrane fusion. This model is verified by reproducing the lipid phase behavior. The classical stalk model has been visited and modified. The tilt deformation of the lipids and the noncircular shape of the stalk are supported. The stalk is shown to undergo asymmetric expansion to form the trans-monolayers contact (TMC). Unlike previous models, an energy barrier between the stalk and the TMC has been identified, implying that the TMC should be a metastable formation. This shows good agreement with the fusion experiments. Two typical elastic continuum models are compared with our result and possible modifications to the two elastic models are suggested. The effect of spontaneous curvature of lipid on selection of fusion pathway is also examined. It is observed that a bent stalk with pore or an inverted micellar intermediate will have more chance to occur than traditional stalk when the spontaneous curvature of the lipid becomes more negative.