We present an algorithm to reconstruct atomistic structures from their corresponding coarse-grained (CG) representations and its implementation into the freely available molecular dynamics (MD) program package GROMACS. The central part of the algorithm is a simulated annealing MD simulation in which the CG and atomistic structures are coupled via restraints. A number of examples demonstrate the...

We use the three-dimensional Mercedes-Benz model for water and Monte Carlo simulations to study the structure and thermodynamics of the hydrophobic interaction. Radial distribution functions are used to classify different cases of the interaction, namely, contact configurations, solvent separated configurations, and desolvation configurations. The temperature dependence of these cases is shown ...

The quasicontinuum method is a coarse-graining technique for reducing the complexity of atomistic simulations in a static and quasistatic setting. In this paper we give an a-priori error analysis for the quasicontinuum method in one dimension. We consider atomistic models with Lennard–Jones type long range interactions and a practical QC formulation. First, we prove the existence, the local uni...

Quantitative understanding of the kinetics of fibril formation and the molecular mechanism of transition from monomers to fibrils is needed to obtain insights into the growth of amyloid fibrils and more generally self-assembly multisubunit protein complexes. Significant advances using computations of protein aggregation in a number of systems have established generic and sequence-specific aspec...

Force-based atomistic-continuum hybrid methods are the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. For this reason, and due to their algorithmic simplicity, force-based coupling methods have become a popular class of atomistic-continuum hybrid models as well as other types of multiphysics models. However, the recently disco...

With computer simulations of >100,000 atoms, the mechanism for the hydrophobic collapse of an idealized hydrated chain was obtained by tiling space with (0.2 nm)(3) cubes and projecting the atomistic water molecule positions onto this grid. With the coarse-grained field thus defined, the string method in collective variables was used to compute a minimum free-energy pathway (MFEP) for the colla...