Multiscale coarse graining of liquid-state systems

A methodology is described to systematically derive coarse-grained (CG) force fields for molecular liquids from the underlying atomistic-scale forces. The coarse graining of an interparticle force field is accomplished by the application of a force-matching method to the trajectories and forces obtained from the atomistic trajectory and force data for the CG sites of the targeted system. The CG sites can be associated with the centers of mass of atomic groups because of the simplicity in the evaluation of forces acting on these sites from the atomistic data. The resulting system is called a multiscale coarse-grained (MS-CG) representation. The MS-CG method for liquids is applied here to water and methanol. For both liquids one-site and two-site CG representations without an explicit treatment of the long-ranged electrostatics have been derived. In addition, for water a two-site model having the explicit long-ranged electrostatics has been developed. To improve the thermodynamic properties (e.g., pressure and density) for the MS-CG models, the constraint for the instantaneous virial was included into the force-match procedure. The performance of the resulting models was evaluated against the underlying atomistic simulations and experiment. In contrast with existing approaches for coarse graining of liquid systems, the MS-CG approach is general, relies only on the interatomic interactions in the reference atomistic system.     

Systematic implicit solvent coarse graining of dimyristoylphosphatidylcholine lipids

We have used systematic structure‐based coarse graining to derive effective site–site potentials for a 10‐site coarse‐grained dimyristoylphosphatidylcholine (DMPC) lipid model and investigated their state point dependence. The potentials provide for the coarse‐grained model the same site–site radial distribution functions, bond and angle distributions as those computed in atomistic simulations carried out at four different lipid–water molar ratios. It was shown that there is a non‐negligible dependence of the effective potentials on the concentration at which they were generated, which is also manifested in the properties of the lipid bilayers simulated using these potentials. Thus, effective potentials computed at low lipid concentration favor to more condensed and ordered structure of the bilayer with lower average area per lipid, while potentials obtained at higher lipid concentrations provide more fluid‐like structure. The best agreement with the reference data and experiment was achieved using the set of potentials derived from atomistic simulations at 1:30 lipid:water molar ratio providing fully saturated hydration of DMPC lipids. Despite theoretical limitations of pairwise coarse‐grained potentials expressed in their state point dependence, all the resulting potentials provide a stable bilayer structure with correct partitioning of different lipid groups across the bilayer as well as acceptable values of the average lipid area, compressibility and orientational ordering. In addition to bilayer simulations, the model has proven its robustness in modeling of self‐aggregation of lipids from randomly dispersed solution to ordered bilayer structures, bicelles, and vesicles. © 2014 Wiley Periodicals, Inc.     

Semi-bottom-up coarse graining of water based on microscopic simulations

The generalized dissipative particle dynamics (DPD) equation derived from the generalized Langevin equation under Markovian approximations is used to simulate coarse-grained (CG) water cells. The mean force and the friction coefficients in the radial and transverse directions needed for DPD equation are obtained directly from the all atomistic molecular dynamics (AAMD) simulations. But the dissipative friction forces are overestimated in the Markovian approximation, which results in wrong dynamic properties for the CG water in the DPD simulations. To account for the non-Markovian dynamics, a rescaling factor is introduced to the friction coefficients. The value of the factor is estimated by matching the diffusivity of water. With this semi-bottom-up mapping method, the radial distribution function, the diffusion constant, and the viscosity of the coarse-grained water system computed with DPD simulations are all in good agreement with AAMD results. It bridges the microscopic level and mesoscopic level with consistent length and time scales.     

Effective thermostat induced by coarse graining of simple point charge water

We investigate how the transport properties of a united atom fluid with a dissipative particle dynamics thermostat depend on the functional form and magnitude of both the conservative and the stochastic interactions. We demonstrate how the thermostat strongly affects the hydrodynamics, especially diffusion, viscosity, and local escape times. As model system we use simple point charge (SPC) water, from which projected trajectories are used to determine the effective interactions in the united atom model. The simulation results support our argument that the thermostat should be viewed as an integral part of the coarse-grained dynamics rather than a tool for approaching thermal equilibrium. As our main result we show that the united atom model with the adjusted effective interactions approximately reproduces the diffusion constant and the viscosity of the underlying detailed SPC water model.     

Statistical mechanics in rotationally inelastic scattering of molecules from surface within the dynamical Lie algebraic method

The dynamical Lie algebraic (DLA) method of Alhassid and Levine [Phys. Rev. A 18 (1978) 89] is applied to statistical mechanics in rotationally inelastic scattering of molecules from surfaces. Specifically, the method is generalized to include the motion of surface atoms, i.e., phonons. For given Hamiltonian and initial state, the set of constraints required to obtain the solution of the motion equations is determined by an algebraic procedure. It is furthermore found possible to derive the motion equations for the mean values of the constraints. Application of the method to the scattering of NO molecules from a Pt(1 1 1) surface is made. The mean values of the final energies of NO molecules scattered from the surface obtained using the DLA method are in good agreement with experimental results in qualitative trends. The DLA method thus appears to have a wide range of validity for describing the statistical mechanics of the gas-surface scattering.     

A coarse graining method for the identification of transition rates between molecular conformations

The coarse graining method to be advocated in this paper consists of two main steps. First, the propagation of an ensemble of molecular states is described as a Markov chain by a transition probability matrix in a finite state space. Second, we obtain metastable conformations by an aggregation of variables via Robust Perron Cluster Analysis (PCCA+). Up to now, it has been an open question as to how this coarse graining in space can be transformed to a coarse graining of the Markov chain while preserving the essential dynamic information. In this article, we construct a coarse matrix that is the correct propagator in the space of conformations. This coarse graining procedure carries over to rate matrices and allows to extract transition rates between molecular conformations. This approach is based on the fact that PCCA+ computes molecular conformations as linear combinations of the dominant eigenvectors of the transition matrix.     

Statistical mechanics of energy transfer in gas–surface scattering: A dynamical Lie algebraic approach

The dynamical Lie algebraic (DLA) method is used to describe statistical mechanics of energy transfer in rotationally inelastic molecule–surface scattering. Statistical average values of an observable for the scattering system are calculated in terms of density operator formalism in statistical mechanics. Employing a cubic expansion procedure of molecule–surface interaction potential leads to generation of a dynamical Lie algebra. Thus these statistical average values as a function of the group parameters can be obtained analytically in this formulation. The group parameters can be found from solving a set of coupled nonlinear differential equations. The DLA method, which has no need for determination of transition probabilities in advance as made routinely in the calculation, offers an efficient alternative to the method for computing the statistical average values. This method is much less computationally intensive because most of calculations can be analytically carried out. The average final rotational energies and their dependence on the main dynamic variables and the average interaction potential are presented for the rotationally inelastic scattering of NO molecules from a flat, static Ag(111) surface. Direct comparison is made between the predictions of this model calculation and experiment. The model reproduces well the degree of rotational excitation and correlation between the average final translational and the average rotational energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Topological coarse graining of polymer chains using wavelet-accelerated Monte Carlo. II. Self-avoiding chains

In the preceding paper [A. E. Ismail, G. C. Rutledge, and G. Stephanopoulos J. Chem. Phys. (in press)] we introduced wavelet-accelerated Monte Carlo (WAMC), a coarse-graining methodology based on the wavelet transform, as a method for sampling polymer chains. In the present paper, we extend our analysis to consider excluded-volume effects by studying self-avoiding chains. We provide evidence that the coarse-grained potentials developed using the WAMC method obey phenomenological scaling laws, and use simple physical arguments for freely jointed chains to motivate these laws. We show that coarse-grained self-avoiding random walks can reproduce results obtained from simulations of the original, more-detailed chains to a high degree of accuracy, in orders of magnitude less time.     

Particle-based multiscale coarse graining with density-dependent potentials: application to molecular crystals (hexahydro-1,3,5-trinitro-s-triazine)

We describe the development of isotropic particle-based coarse-grain models for crystalline hexahydro-1,3,5-trinitro-s-triazine (RDX). The coarse graining employs the recently proposed multiscale coarse-graining (MS-CG) method, which is a particle-based force-matching approach for deriving free-energy effective interaction potentials. Though one-site and four-site coarse-grain (CG) models were parameterized from atomistic simulations of non-ordered (molten and ambient temperature amorphous) systems, the focus of the paper is a detailed study of the one-site model with a brief recourse to the four-site model. To improve the ability of the one-site model to be applied to crystalline phases at various pressures, it was found necessary to include explicit dependence on a particle density, and a new theory of local density-dependent MS-CG potentials is subsequently presented. The density-dependency is implemented through interpolation of MS-CG force fields derived at a preselected set of reference densities. The computationally economical procedure for obtaining the reference force fields starting from the interaction at ambient density is also described. The one-site MS-CG model adequately describes the atomistic lattice structure of α-RDX at ambient and high pressures, elastic and vibrational properties, pressure-volume curve up to P = 10 GPa, and the melting temperature. In the molten state, the model reproduces the correct pair structure at different pressures as well as higher order correlations. The potential of the MS-CG model is further evaluated in simulations of shocked crystalline RDX.     

Statistical mechanics of three-dimensional vesicles

We consider a lattice model of three-dimensional vesicles in which the boundary of the vesicle is a self-avoiding plaquette surface, homeomorphic to a sphere. Surfaces with fixed area can enclose a variety of different volumes and we associate a fugacity with the enclosed volume to mimic the effect of a pressure difference across the surface. Pairs of plaquettes which share a common edge can be in the same plane or normal to each other and we associate a fugacity with adjacent pairs of plaquettes at right angles to represent a surface stiffness term. We discuss the behaviour of the surfaces in the infinite surface area limit, as a function of these two fugacities.     

Topological coarse graining of polymer chains using wavelet-accelerated Monte Carlo. I. Freely jointed chains

We introduce a new, topologically-based method for coarse-graining polymer chains. Based on the wavelet transform, a multiresolution data analysis technique, this method assigns a cluster of particles to a coarse-grained bead located at the center of mass of the cluster, thereby reducing the complexity of the problem by dividing the simulation into several stages, each with a fraction of the number of beads as the overall chain. At each stage, we compute the distributions of coarse-grained internal coordinates as well as potential functions required for subsequent simulation stages. In this paper, we present the basic algorithm, and apply it to freely jointed chains; the companion paper describes its applications to self-avoiding chains.     

Statistical mechanics of worm-like polymers from a new generating function

We present a mathematical approach to the worm-like chain model of semiflexible polymers. Our method is built on a novel generating function from which all the properties of the model can be derived. Moreover, this approach satisfies the local inextensibility constraint exactly. In this paper, we focus on the lowest order contribution to the generating function and derive explicit analytical expressions for the characteristic function, polymer propagator, single chain structure factor, and mean square end-to-end distance. These analytical expressions are valid for polymers with any degree of stiffness and contour length. We find that our calculations are able to capture the fully flexible and infinitely stiff limits of the aforementioned quantities exactly while providing a smooth and approximate crossover behavior for intermediate values of the stiffness of the polymer backbone. In addition, our results are in very good quantitative agreement with the exact and approximate results of five other treatments of semiflexible polymers.     

Statistical mechanics of helical wormlike chain model

We investigate the statistical mechanics of polymers with bending and torsional elasticity described by the helical wormlike model. Noticing that the energy function is factorizable, we provide a numerical method to solve the model using a transfer matrix formulation. The tangent-tangent and binormal-binormal correlation functions have been calculated and displayed rich profiles which are sensitive to the combination of the temperature and the equilibrium torsion. Their behaviors indicate that there is no finite temperature Lifshitz point between the disordered and helical phases. The asymptotic behavior at low temperature has been investigated theoretically and the predictions fit the numerical results very well. Our analysis could be used to understand the statics of dsDNA and other chiral polymers.