Do you want to publish a course? Click here

Advanced multicanonical Monte Carlo methods for efficient simulations of nucleation processes of polymers

166   0   0.0 ( 0 )
 Added by Stefan Schnabel
 Publication date 2012
  fields Physics
and research's language is English




Ask ChatGPT about the research

The investigation of freezing transitions of single polymers is computationally demanding, since surface effects dominate the nucleation process. In recent studies we have systematically shown that the freezing properties of flexible, elastic polymers depend on the precise chain length. Performing multicanonical Monte Carlo simulations, we faced several computational challenges in connection with liquid-solid and solid-solid transitions. For this reason, we developed novel methods and update strategies to overcome the arising problems. We introduce novel Monte Carlo moves and two extensions to the multicanonical method.



rate research

Read More

215 - Ajay Jasra , Kody Law , 2017
This article reviews the application of advanced Monte Carlo techniques in the context of Multilevel Monte Carlo (MLMC). MLMC is a strategy employed to compute expectations which can be biased in some sense, for instance, by using the discretization of a associated probability law. The MLMC approach works with a hierarchy of biased approximations which become progressively more accurate and more expensive. Using a telescoping representation of the most accurate approximation, the method is able to reduce the computational cost for a given level of error versus i.i.d. sampling from this latter approximation. All of these ideas originated for cases where exact sampling from couples in the hierarchy is possible. This article considers the case where such exact sampling is not currently possible. We consider Markov chain Monte Carlo and sequential Monte Carlo methods which have been introduced in the literature and we describe different strategies which facilitate the application of MLMC within these methods.
A coarse-grained simulation model eliminates microscopic degrees of freedom and represents a polymer by a simplified structure. A priori, two classes of coarse-grained models may be distinguished: those which are designed for a specific polymer and reflect the underlying atomistic details to some extent, and those which retain only the most basic features of a polymer chain (chain connectivity, short-range excluded-volume interactions, etc.). In this review we mainly focus on the second class of generic polymer models, while the first class of specific coarse-grained models is only touched upon briefly.
We present a novel Ensemble Monte Carlo Growth method to sample the equilibrium thermodynamic properties of random chains. The method is based on the multicanonical technique of computing the density of states in the energy space. Such a quantity is temperature independent, and therefore microcanonical and canonical thermodynamic quantities, including the free energy, entropy, and thermal averages, can be obtained by re-weighting with a Boltzmann factor. The algorithm we present combines two approaches: the first is the Monte Carlo ensemble growth method, where a population of samples in the state space is considered, as opposed to traditional sampling by long random walks, or iterative single-chain growth. The second is the flat-histogram Monte Carlo, similar to the popular Wang-Landau sampling, or to multicanonical chain-growth sampling. We discuss the performance and relative simplicity of the proposed algorithm, and we apply it to known test cases.
The Dynamic Monte Carlo (DMC) method is an established molecular simulation technique for the analysis of the dynamics in colloidal suspensions. An excellent alternative to Brownian Dynamics or Molecular Dynamics simulation, DMC is applicable to systems of spherical and/or anisotropic particles and to equilibrium or out-of-equilibrium processes. In this work, we present a theoretical and methodological framework to extend DMC to the study of heterogeneous systems, where the presence of an interface between coexisting phases introduces an additional element of complexity in determining the dynamic properties. In particular, we simulate a Lennard-Jones fluid at the liquid-vapor equilibrium and determine the diffusion coefficients in the bulk of each phase and across the interface. To test the validity of our DMC results, we also perform Brownian Dynamics simulations and unveil an excellent quantitative agreement between the two simulation techniques.
Understanding the rheology of colloidal suspensions is crucial in the formulation of a wide selection of industry-relevant products. To characterise the viscoelastic behaviour of these soft materials, one can analyse the microscopic dynamics of colloidal tracers diffusing through the host fluid and generating local deformations and stresses. This technique, referred to as microrheology, links the bulk rheology of fluids to the microscopic dynamics at the particle scale. If tracers are subjected to external forces, rather than freely diffusing, it is called active microrheology. Motivated by the impact of microrheology in providing information on local structure in complex systems such as colloidal glasses, active matter or biological systems, we have extended the dynamic Monte Carlo (DMC) technique to investigate active microrheology in colloids. The original DMC framework, able to accurately describe the Brownian dynamics of colloids at equilibrium, is here reconsidered and expanded to describe the effects of an external force pulling a tracer embedded in isotropic colloidal suspensions at different densities. To this end, we studied the dynamics of a spherical tracer dragged by a constant external force through a bath of spherical and rod-like particles of comparable size. We could extract valuable details on its effective friction coefficient, being constant at small and large values of the external force, but otherwise displaying a nonlinear behaviour that indicates the occurrence of a force-thinning regime. Our DMC simulation results are in excellent quantitative agreement with past Langevin dynamics simulations and theoretical works for the bath of spherical colloids. The bath of rod-like particles is studied in the isotropic phase, and displays an example where DMC is more convenient than Brownian or Langevin dynamics, in this case in dealing with particle rotation.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا