A new approach is presented for obtaining coarse-grained (CG) force fields from fully atomistic molecular dynamics (MD) trajectories. The method is demonstrated by applying it to derive a CG model for the dimyristoylphosphatidylcholine (DMPC) lipid bilayer. The coarse-graining of the interparticle force field is accomplished by an application of a force-matching procedure to the force data obtained from an explicit atomistic MD simulation of the biomolecular system of interest. Hence, the method is termed a "multiscale" CG (MS-CG) approach in which explicit atomistic-level forces are propagated upward in scale to the coarse-grained level. The CG sites in the lipid bilayer application were associated with the centers-of-mass of atomic groups because of the simplicity in the evaluation of the forces acting on them from the atomistic data. The resulting CG lipid bilayer model is shown to accurately reproduce the structural properties of the phospholipid bilayer.

Coarse-grained ͑CG͒ models provide a computationally efficient method for rapidly investigating the long time-and length-scale processes that play a critical role in many important biological and soft matter processes. Recently, Izvekov and Voth introduced a new multiscale coarse-graining ͑MS-CG͒ method ͓J. Phys. Chem. B 109, 2469 ͑2005͒; J. Chem. Phys. 123, 134105 ͑2005͔͒ for determining the effective interactions between CG sites using information from simulations of atomically detailed models. The present work develops a formal statistical mechanical framework for the MS-CG method and demonstrates that the variational principle underlying the method may, in principle, be employed to determine the many-body potential of mean force ͑PMF͒ that governs the equilibrium distribution of positions of the CG sites for the MS-CG models. A CG model that employs such a PMF as a "potential energy function" will generate an equilibrium probability distribution of CG sites that is consistent with the atomically detailed model from which the PMF is derived. Consequently, the MS-CG method provides a formal multiscale bridge rigorously connecting the equilibrium ensembles generated with atomistic and CG models. The variational principle also suggests a class of practical algorithms for calculating approximations to this many-body PMF that are optimal. These algorithms use computer simulation data from the atomically detailed model. Finally, important generalizations of the MS-CG method are introduced for treating systems with rigid intramolecular constraints and for developing CG models whose equilibrium momentum distribution is consistent with that of an atomically detailed model.

In order to introduce flexibility into the simple point-charge (SPC) water model, the impact of the intramolecular degrees of freedom on liquid properties was systematically studied in this work as a function of many possible parameter sets. It was found that the diffusion constant is extremely sensitive to the equilibrium bond length and that this effect is mainly due to the strength of intermolecular hydrogen bonds. The static dielectric constant was found to be very sensitive to the equilibrium bond angle via the distribution of intermolecular angles in the liquid: A larger bond angle will increase the angle formed by two molecular dipoles, which is particularly significant for the first solvation shell. This result is in agreement with the work of Hochtl et al. [J. Chem. Phys. 109, 4927 (1998)]. A new flexible simple point-charge water model was derived by optimizing bulk diffusion and dielectric constants to the experimental values via the equilibrium bond length and angle. Due to the large sensitivities, the parametrization only slightly perturbs the molecular geometry of the base SPC model. Extensive comparisons of thermodynamic, structural, and kinetic properties indicate that the new model is much improved over the standard SPC model and its overall performance is comparable to or even better than the extended SPC model.

Ionic liquids have many promising industrial applications. The side-chain length of the cations has been known to significantly influence the physical and chemical properties of those liquids, especially liquid crystal formation. 1 Because there are numerous species of possible ionic liquids, it is of great interest to understand the general physical picture behind the effect of varying the length of side chains of cations.In this work, a multiscale coarse-graining (MS-CG) method 2 has been extended to explore the effect of various cation side-chain lengths in ionic liquids. This method allows for large systems to be simulated for long times, thus revealing features of the system that are difficult to see using conventional all-atom molecular dynamics (MD) simulations. Simulations with the MS-CG models show that, with sufficient side-chain length, neutral tail groups of cations aggregate to form spatially heterogeneous domains of the tails, while the charged headgroups and anions distribute as uniformly as possible due to the strong electrostatic interactions. The geometrical constraints for head and tail groups of cations result in a novel balanced liquid crystal-like structure at suitable temperatures. This physical picture can qualitatively explain the experimentally observed ionic liquid crystal formation, the transition from liquid crystal to isotropic ionic liquid, and the changes of structural, dynamic, and thermodynamic properties when varying the sidechain length.The MS-CG method 2 has been applied here to the EMIM + NO 3 -ionic liquid to develop coarse-grained (CG) models at T ) 400 and 700 K, respectively. Details of the CG models are given in the Supporting Information. As shown in Figure 1, the nitrate anion has been coarse-grained as Site D, while the aromatic ring of the cation as Site A, the methyl group as Site B, the methylene and methyl groups on the alkyl chain as Sites C and E, respectively. To study the effect of chain-length elongation, the alkyl chain is then extended with two more methylene groups (Sites C) to form BMIM + NO 3 -. The partial charges are assigned as the numbers in Figure 1. The bonded parameters for the sites on the alkyl chain are assigned with the parameters for the sp 3 carbon (CT) sites given in ref 3. Although this model is a coarse-grained one, the qualitative results should not depend on the details of the force field parameters. For convenience, the ionic liquid systems will be denoted by the number of carbons on the alkyl chain, e.g., EMIM + NO 3 -as C 2 and BMIM + NO 3 -as C 4 . The C 4 MS-CG systems with 64, 400, and 800 ion pairs have been simulated by using the DL_POLY program 4 at T ) 700 K. The aggregation of tail groups was determined by visual examination for all three sizes of simulation. One snapshot of the simulation for 400 ion pairs is shown in Figure 2. Headgroups and anions are seen to distribute relatively homogeneously, but tail groups aggregate together and form several spatially heterogeneous domains.The heterogeneity of the tail groups is char...

The dynamics and energetics of an excess proton in bulk phase water are examined computationally with a special emphasis on a quantum-dynamical treatment of the nuclear motion. The potential model used, the recently developed multistate empirical valence bond (MS-EVB) approach [U. W. Schmitt and G. A. Voth, J. Phys. Chem. B 102, 5547 (1998)], is also further refined and described in more detail. The MS-EVB model takes into account the interaction of an exchange charge distribution of the charge-transfer complex with the polar solvent, which qualitatively changes the nature of the solvated complex. Classical and quantum molecular dynamics simulations of the excess proton in bulk phase water reveal that quantization of the nuclear degrees of freedom results in an increased stabilization of the solvated H5O2+ (Zundel) cation relative to the H9O4+ (Eigen) cation, though the latter is still more stable, and that a species intermediate between the two also exists. The quantum proton transport rate, which is evaluated by the centroid molecular dynamics approach, is found to be on the order of two times faster compared to a purely classical treatment of the system and in good agreement with the experimental value. Calculation of the hydrogen-bonding lifetime beyond the first solvation shell of the excess proton reveals a similar quantum enhancement factor compared to the classical regime.

We propose and implement an alternative approach to the original Car–Parrinello method where the density matrix elements (instead of the molecular orbitals) are propagated together with the nuclear degrees of freedom. Our new approach has the advantage of leading to an O(N) computational scheme in the large system limit. Our implementation is based on atom-centered Gaussian orbitals, which are especially suited to deal effectively with general molecular systems. The methodology is illustrated by applications to the three-body dissociation of triazine and to the dynamics of a cluster of a chloride ion with 25 water molecules.

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.

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.