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.

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.

Conditional and time-dependent radial distribution functions reveal the details of the water structure surrounding the hydronium during the proton mobility process. Using this methodology for classical multistate empirical valence bond (MS-EVB) and ab initio molecular dynamics trajectories, as well as quantal MS-EVB trajectories, we supply statistical proof that proton hops in liquid water occur by a transition from the H3O+[3H2O] Eigen-complex, via the H5O2+ Zundel-complex, to a H3O+[3H2O] centered on a neighboring water molecule. In the "resting period" before a transition, there is a distorted hydronium with one of its water ligands at a shorter distance and another at a longer distance than average. The identity of this "special partner" interchanges rapidly within the three first-shell water ligands. This is coupled to cleavage of an acceptor-type hydrogen bond. Just before the transition, a partner is selected by an additional translation of the H3O+ moiety in its direction, possibly enabled by loosening of donor-type hydrogen bonds on the opposite side. We monitor the transition in real time, showing how the average structure is converted to a distorted H5O2+ cation constituting the transitional complex for proton hopping between water molecules.