Bernd Ensing

Energy Conservation in Adaptive Hybrid Atomistic/Coarse-Grain Molecular Dynamics.

Multiscale computer simulation algorithms are required to describe complex molecular systems with events occurring over a range of time and length scales. True multiscale simulations must solve the interface, or hand-shaking, problem of coupling together different levels of description in different spatial regions of the system. If the spatial regions of different resolution move over time, or if material is allowed to flow over the inter-region boundaries, a mechanism must be introduced into the multiscale algorithm to allow material to dynamically change its representation. While such a mechanism is highly desirable in many instances, it is fraught with technical difficulties. Here, we present a molecular dynamics simulation algorithm which is multiscale in both time and space. We supplement the potential and kinetic energy expressions with auxiliary terms in order to recover the total energy as a conserved quantity, even when the total number of degrees of freedom changes during the simulation. This is crucial for a proper assessment of the quality of adaptive hybrid algorithms, and in particular, it allows us to tune the hierarchy of RESPA levels to optimize the integration scheme.

Recent progress in adaptive multiscale molecular dynamics simulations of soft matter

Understanding mesoscopic phenomena in terms of the fundamental motions of atoms and electrons poses a severe challenge for molecular simulation. This challenge is being met by multiscale modeling techniques that aim to bridge between the microscopic and mesoscopic time and length scales. In such techniques different levels of theory are combined to describe a system at a number of scales or resolutions. Here we review recent advancements in adaptive hybrid simulations, in which the different levels are used in separate spatial domains and matter can diffuse from one region to another with an accompanying resolution change. We discuss what it means to simulate such a system, and how to enact the resolution changes. We show how to construct efficient adaptive hybrid quantum mechanics/molecular mechanics (QM/MM) and atomistic/coarse grain (AA/CG) molecular dynamics methods that use an intermediate healing region to smoothly couple the regions together. This coupling is formulated to use only the native forces inherent to each region. The total energy is conserved through the use of auxiliary bookkeeping terms. Error control, and the choice of time step and healing region width, is obtained by careful analysis of the energy flow between the different representations. We emphasize the CG → AA reverse mapping problem and show how this problem is resolved through the use of rigid atomistic fragments located within each CG particle whose orientation is preconditioned for a possible resolution change through a rotational dynamics scheme. These advancements are shown to enable the adaptive hybrid multiscale molecular dynamics simulation of macromolecular soft matter systems.

Toward a Practical Method for Adaptive QM/MM Simulations

We present an accurate adaptive multiscale molecular dynamics method that will enable the detailed study of large molecular systems that mimic experiment. The method treats the reactive regions at the quantum mechanical level and the inactive environment regions at lower levels of accuracy, while at the same time molecules are allowed to flow across the border between active and environment regions. Among many other things, this scheme affords accurate investigation of chemical reactions in solution. A scheme like this ideally fulfills the key criteria applicable to all molecular dynamics simulations: energy conservation and computational efficiency. Approaches that fulfill both criteria can, however, result in complicated potential energy surfaces, creating rapid energy changes when the border between regions is crossed. With the difference-based adaptive solvation potential, a simple approach is introduced that meets the above requirements and reduces fast fluctuations in the potential to a minimum. In cases where none of the current adaptive QM/MM potentials are able to properly describe the system under investigation, we use a continuous force scheme instead, which, while no longer energy conserving, still retains a related conserved quantity along the trajectory. We show that this scheme does not introduce a significant temperature drift on time scales feasible for QM/MM simulations.

DFT Study of the Active Intermediate in the Fenton Reaction

Density functional theory has been used to investigate the nature of the oxidizing agent in the Fenton reaction. Starting from the primary intermediate [FeII(H2O)5H2O2]2+, we show that the oxygen-oxygen bond breaking mechanism has a small activation energy and could therefore demonstrate the catalytic effect of the metal complex. The O-O bond cleavage of the coordinated H2O2, however, does not lead to a free hydroxyl radical. Instead, the leaving hydroxyl radical abstracts a hydrogen from an adjacent coordinated water leading to the formation of a second Fe-OH bond and of a water molecule. Along this reaction path the primary intermediate transforms into the [FeIV(H2O)4(OH)2]2+ complex and in a second step into a more stable high valent ferryl-oxo complex [FeIV(H2O)5O]2+. We show that the energy profile along the reaction path is strongly affected by the presence of an extra water molecule located near the iron complex. The alternative intermediate [FeII(H2O)4(OOH-)(H3O+)]2+ suggested in the literature has been also investigated, but it is found to be unstable against the primary intermediate. Our results support a picture in which an FeIV-oxo complex is the most likely candidate as the active intermediate in the Fenton reaction, as indeed first proposed by Bray and Gorin already in 1932.

Chemical Involvement of Solvent Water Molecules in Elementary Steps of the Fenton Oxidation Reaction We gratefully acknowledge the helpful discussions with Michiel Gribnau (Unilever-Vlaardingen) and we thank the Netherlands Organization for Scientific Research (NWO) for support through the PPM-CMS program and the NCF for providing computer time

The spontaneous formation of the contested ferryl ion is evident in ab initio molecular dynamics calculations on the FeII /H2 O2 system in aqueous solution (Fenton reagents). Not only is the ferryl ion preferred over the hydroxyl radical as the active oxidative species, but the solvent water molecules play a crucial role in the mechanism. The picture shows the unit cell containing the iron complex surrounded by solvent water molecules 1.8 ps after the start of the simulation, when the ferryl ion is being formed (blue: iron, red: oxygen, white: hydrogen).

Hydrolysis of cisplatin—a first-principles metadynamics study

Cisplatin, or cis-[Pt(NH(3))(2)Cl(2)], was the first member of a new revolutionary class of anticancer drugs that is still used today for the treatment of a wide variety of cancers. The mode of action of cisplatin starts inside the cell with the hydrolysis of Pt-Cl bonds to form a Pt-aqua complex. The solvent environment plays an essential role in many biochemical processes in general, and is expected to have a particular strong effect on the activation (hydrolysis) of cisplatin and cisplatin derivatives. To investigate these solvent effects, we have studied the explicit solvent structures during cisplatin hydrolysis by means of Car-Parrinello molecular dynamics simulations. Since hydrolysis is an activated process, and thus a rare event on the simulation timescale, we have applied the metadynamics sampling technique to map out the free energy landscape from which the reaction mechanism and activation free energy are obtained. Our simulations show that hydrogen bonding between solvent water molecules and metal complexes in the hydrolyzed product systems is stronger than that in the reactant cisplatin system. In addition, the free energy profiles from our metadynamics simulations for the cisplatin hydrolysis shows that the second hydrolysis of cisplatin is thermodynamically favourable, which is in good agreement with experimental results and previous static density functional theory calculations. The reactant channels for both hydrolysis steps are rather wide and flat, indicative of a continuous spectrum of allowed mechanisms with no strong preference for either concerted dissociative or concerted associative pathways. Three or five coordinated metastable intermediates do not exist in aqueous solution.

Nonlinear reaction coordinate analysis in the reweighted path ensemble

We present a flexible nonlinear reaction coordinate analysis method for the transition path ensemble based on the likelihood maximization approach developed by Peters and Trout [J. Chem. Phys. 125, 054108 (2006)]. By parametrizing the reaction coordinate by a string of images in a collective variable space, we can optimize the likelihood that the string correctly models the committor data obtained from a path sampling simulation. The collective variable space with the maximum likelihood is considered to contain the best description of the reaction. The use of the reweighted path ensemble [J. Rogal et al., J. Chem. Phys. 133, 174109 (2010)] allows a complete reaction coordinate description from the initial to the final state. We illustrate the method on a z-shaped two-dimensional potential. While developed for use with path sampling, this analysis method can also be applied to regular molecular dynamics trajectories.

Complex Reaction Environments and Competing Reaction Mechanisms in Zeolite Catalysis: Insights from Advanced Molecular Dynamics

The methanol-to-olefin process is a showcase example of complex zeolite-catalyzed chemistry. At real operating conditions, many factors affect the reactivity, such as framework flexibility, adsorption of various guest molecules, and competitive reaction pathways. In this study, the strength of first principle molecular dynamics techniques to capture this complexity is shown by means of two case studies. Firstly, the adsorption behavior of methanol and water in H-SAPO-34 at 350 °C is investigated. Hereby an important degree of framework flexibility and proton mobility was observed. Secondly, the methylation of benzene by methanol through a competitive direct and stepwise pathway in the AFI topology was studied. Both case studies clearly show that a first-principle molecular dynamics approach enables unprecedented insights into zeolite-catalyzed reactions at the nanometer scale to be obtained.

On the Polarity of Buckminsterfullerene with a Water Molecule Inside

Since the recent achievement of Kurotobi and Murata to capture a water molecule in a C(60) fullerene (Science 2011, 333, 613), there has been a debate about the properties of this H(2)O@C(60) complex. In particular, the polarity of the complex, which is thought to be underlying the easy separation of H(2)O@C(60) from the empty fullerene by HPLC, was calculated and found to be almost equal to that of an isolated water molecule. Here we present our detailed analysis of the charge distribution of the water-encapsulated C(60) complex, which shows that the polarity of the complex is, with 0.5 ± 0.1 D, indeed substantial, but significantly smaller than that of H(2)O. This may have important implications for the aim to design water-soluble and biocompatible fullerenes.

The reweighted path ensemble

We introduce a reweighting scheme for the path ensembles in the transition interface sampling framework. The reweighting allows for the analysis of free energy landscapes and committor projections in any collective variable space. We illustrate the reweighting scheme on a two dimensional potential with a nonlinear reaction coordinate and on a more realistic simulation of the Trp-cage folding process. We suggest that the reweighted path ensemble can be used to optimize possible nonlinear reaction coordinates.