Highly parallelizable path sampling with minimal rejections using asynchronous replica exchange and infinite swaps.

Capturing rare yet pivotal events poses a significant challenge for molecular simulations. Path sampling provides a unique approach to tackle this issue without altering the potential energy landscape or dynamics, enabling recovery of both thermodynamic and kinetic information. However, despite its exponential acceleration compared to standard molecular dynamics, generating numerous trajectories can still require a long time. By harnessing our recent algorithmic innovations-particularly subtrajectory moves with high acceptance, coupled with asynchronous replica exchange featuring infinite swaps-we establish a highly parallelizable and rapidly converging path sampling protocol, compatible with diverse high-performance computing architectures. We demonstrate our approach on the liquid-vapor phase transition in superheated water, the unfolding of the chignolin protein, and water dissociation. The latter, performed at the ab initio level, achieves comparable statistical accuracy within days, in contrast to a previous study requiring over a year.

Effect of the Ion, Solvent, and Thermal Interaction Coefficients on Battery Voltage.

In order to increase the adoption of batteries for sustainable transport and energy storage, improved charging and discharging capabilities of lithium-ion batteries are necessary. To achieve this, accurate data that describe the internal state of the cells are essential. Several models have been derived, and transport coefficients have been reported for use in these models. We report for the first time a complete set of transport coefficients to model the concentration and temperature polarization in a lithium-ion battery ternary electrolyte, allowing us to test common assumptions. We include effects due to gradients in chemical potentials and in temperature. We find that the voltage contributions due to salt and solvent polarization are of the same order of magnitude as the ohmic loss and must be taken into account for more accurate modeling and understanding of battery performance. We report new Soret and Seebeck coefficients and find thermal polarization to be significant in cases relevant to battery research. The analysis is suitable for electrochemical systems, in general.

Ultrasound enhanced diffusion in hydrogels: An experimental and non-equilibrium molecular dynamics study.

Focused ultrasound has experimentally been found to enhance the diffusion of nanoparticles; our aim with this work is to study this effect closer using both experiments and non-equilibrium molecular dynamics. Measurements from single particle tracking of 40 nm polystyrene nanoparticles in an agarose hydrogel with and without focused ultrasound are presented and compared with a previous experimental study using 100 nm polystyrene nanoparticles. In both cases, we observed an increase in the mean square displacement during focused ultrasound treatment. We developed a coarse-grained non-equilibrium molecular dynamics model with an implicit solvent to investigate the increase in the mean square displacement and its frequency and amplitude dependencies. This model consists of polymer fibers and two sizes of nanoparticles, and the effect of the focused ultrasound was modeled as an external oscillating force field. A comparison between the simulation and experimental results shows similar mean square displacement trends, suggesting that the particle velocity is a significant contributor to the observed ultrasound-enhanced mean square displacement. The resulting diffusion coefficients from the model are compared to the diffusion equation for a two-time continuous time random walk. The model is found to have the same frequency dependency. At lower particle velocity amplitude values, the model has a quadratic relation with the particle velocity amplitude as described by the two-time continuous time random walk derived diffusion equation, but at higher amplitudes, the model deviates, and its diffusion coefficient reaches the non-hindered diffusion coefficient. This observation suggests that at higher ultrasound intensities in hydrogels, the non-hindered diffusion coefficient can be used.

Transport coefficients for ion and solvent coupling. The case of the lithium-ion battery electrolyte.

Transport properties are essential for the understanding and modeling of electrochemical cells, in particular complex systems like lithium-ion batteries. In this study, we demonstrate how a certain degree of freedom in the choice of variables allows us to efficiently determine a complete set of transport properties. We apply the entropy production invariance condition to different sets of electrolyte variables and obtain a general set of formulas. We demonstrate the application of these formulas to an electrolyte typical for lithium-ion batteries, 1M lithium hexafluoro-phosphate in a 1:1 wt. % mixture of ethylene and diethyl carbonates. While simplifications can be introduced, they provide inadequate predictions of conductivity and transport numbers, and we argue that a full matrix of Onsager coefficients is needed for adequate property predictions. Our findings highlight the importance of a complete set of transport coefficients for accurate modeling of complex electrochemical systems and the need for careful consideration of the choice of variables used to determine these properties.

The role of pressure and defects in the wurtzite to rock salt transition in cadmium selenide.

Using molecular dynamics and path sampling techniques we investigated the effect of pressure and defects in the wurtzite to rock salt transition in cadmium selenide (CdSe). In the pressure range 2-10 GPa, rate constants of transition are in the order of 10-23 to 105 s-1 for the transformation of a relatively small wurtzite crystal consisting of 1024 atoms with periodic boundary conditions. The transition paths predominantly evolve through an intermediate 5-coordinated structure, as reported before, though its typical lifetime within the transition paths is particularly long in the intermediate pressure range (4-6 GPa). The defects were created by removing Cd-Se pairs from an otherwise perfect crystal. The removals were either selected fully randomized or grouped in clusters (cavity creation). We find that the rate of transition due to the defects increases by several orders of magnitude even for a single pair removal. This is caused by a change in the transition mechanism that no longer proceeds via the intermediate 5-coordinated structure, when defects are present. Further, the cavity creation yields a lower rate than the fully randomized removal.

Local initiation conditions for water autoionization.

The pH of liquid water is determined by the infrequent process in which water molecules split into short-lived hydroxide and hydronium ions. This reaction is difficult to probe experimentally and challenging to simulate. One of the open questions is whether the local water structure around a slightly stretched OH bond is actually initiating the eventual breakage of this bond or whether this event is driven by a global ordering that involves many water molecules far away from the reaction center. Here, we investigated the self-ionization of water at room temperature by rare-event ab initio molecular dynamics and obtained autoionization rates and activation energies in good agreement with experiments. Based on the analysis of thousands of molecular trajectories, we identified a couple of local order parameters and show that if a bond stretch occurs when all these parameters are around their ideal range, the chance for the first dissociation step (double-proton jump) increases from [Formula: see text] to 0.4. Understanding these initiation triggers might ultimately allow the steering of chemical reactions.

PyRETIS 2: An improbability drive for rare events.

The algorithmic development in the field of path sampling has made tremendous progress in recent years. Although the original transition path sampling method was mostly used as a qualitative tool to sample reaction paths, the more recent family of interface-based path sampling methods has paved the way for more quantitative rate calculation studies. Of the exact methods, the replica exchange transition interface sampling (RETIS) method is the most efficient, but rather difficult to implement. This has been the main motivation to develop the open-source Python-based computer library PyRETIS that was released in 2017. PyRETIS is designed to be easily interfaced with any molecular dynamics (MD) package using either classical or ab initio MD. In this study, we report on the principles and the software enhancements that are now included in PyRETIS 2, as well as the recent developments on the user interface, improvements of the efficiency via the implementation of new shooting moves, easier initialization procedures, analysis methods, and supported interfaced software. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Non-equilibrium thermodynamics as a tool to compute temperature at the catalyst surface.

The surface temperature is computed for a heterogeneous catalytic reaction model, namely the oxidation of carbon monoxide on platinum. The surface temperature was found using non-equilibrium thermodynamic theory, a theory which provides the proper dependencies between heat and mass fluxes and the reaction rate. The theory predicts a possible coupling between the reaction rate and the thermal driving force and can help extend classical reaction kinetics. In the absence of direct measurements, we explore the coupling numerically. The results are able to capture experimental data reported in the literature, and give new insights into why Arrhenius plots may turn out to be non-linear.

PyRETIS: A well-done, medium-sized python library for rare events.

Transition path sampling techniques are becoming common approaches in the study of rare events at the molecular scale. More efficient methods, such as transition interface sampling (TIS) and replica exchange transition interface sampling (RETIS), allow the investigation of rare events, for example, chemical reactions and structural/morphological transitions, in a reasonable computational time. Here, we present PyRETIS, a Python library for performing TIS and RETIS simulations. PyRETIS directs molecular dynamics (MD) simulations in order to sample rare events with unbiased dynamics. PyRETIS is designed to be easily interfaced with any molecular simulation package and in the present release, it has been interfaced with GROMACS and CP2K, for classical and ab initio MD simulations, respectively. © 2017 Wiley Periodicals, Inc.

Rare event simulations reveal subtle key steps in aqueous silicate condensation.

A replica exchange transition interface sampling (RETIS) study combined with Born-Oppenheimer molecular dynamics (BOMD) is used to investigate the dynamics, thermodynamics and the mechanism of the early stages of the silicate condensation process. In this process, two silicate monomers, of which one is an anionic species, form a negatively charged five-coordinated silicate dimer. In a second stage, this dimer can fall apart again, forming the original monomers, or release a water molecule into the solution. We studied the association and dissociation reaction in the gas phase, and the dissociation and water removal step in the aqueous phase. The results on the aqueous phase dissociation suggest two possible mechanisms. The breakage of the bond between the intermediate oxygen and the five-coordinated silicon is sometimes accompanied by a proton transfer. After dissociation into silicate monomers, the anionic monomer is either the previously four-coordinated silicon or the previously five-coordinated silicon depending on whether the hydrogen transfer occurs or not. Our results show that the mechanism of proton transfer is highly predominant. Water removal simulations also show two possible mechanisms distinguished by the proton transfer reaction path. Proton transfer can occur either via a direct or via a water mediated reaction step. The calculations reveal that although both mechanisms contribute to the water removal process, the direct proton transfer is slightly favorable and occurs roughly in six out of ten occasions.

Structural organization of sterol molecules in DPPC bilayers: a coarse-grained molecular dynamics investigation.

We investigate the structural organization of cholesterol (CHOL) analogues in 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers using coarse-grained molecular dynamics simulations and the MARTINI forcefield. Different sterol molecules are modelled by increasing (CHOLL) or decreasing (CHOLS) the diameter of the sterol beads employed in the MARTINI model of CHOL. At high sterol concentrations, (xsterol = 0.5), typical of liquid ordered phases, we find that the sterol arrangement and sterol-DPPC interactions strongly depend on the sterol size. Smaller sterols (CHOLS and CHOL) form linear clusters, while the larger sterols (CHOLL) arrange themselves into disc shaped clusters. By combining structural and dynamical properties we also investigate the So→ Ld transition for the CHOLL and CHOLS sterols. We show that small changes in the sterol size significantly affect the stability of the gel phase with the gel phase stabilized by the small sterols, but destabilized by large sterols. The general dependence of the phase behaviour of the membrane with sterol content is reminiscent of the one observed in naturally occurring membranes. The relevance of our results to understand current cholesterol-bilayer structural models is discussed.

Note: A new truncation correction for the configurational temperature extends its applicability to interaction potentials with a discontinuous force.

We present a simple truncation correction for the configurational temperature which, unlike previous corrections, works even at low truncation values for the shifted and truncated Lennard-Jones potential. The success of the new correction suggests that the expression for the configurational temperature is valid also for interaction potentials with a discontinuous force, given that the discontinuity is properly accounted for.

Michaelis-Menten kinetics under non-isothermal conditions.

We extend the celebrated Michaelis-Menten kinetics description of an enzymatic reaction taking into consideration the presence of a thermal driving force. A coupling of chemical and thermal driving forces is expected from the principle of non-equilibrium thermodynamics, and specifically we obtain an additional term to the classical Michaelis-Menten kinetic equation, which describes the coupling in terms of a single parameter. A companion equation for the heat flux is also derived, which actually can exist even in the absence of a temperature difference. Being thermodynamic in nature, this result is general and independent of the detailed mechanism of the coupling. Conditions for the experimental verification of the new equation are discussed.

Finite-size and truncation effects for microscopic expressions for the temperature at equilibrium and nonequilibrium.

Several expressions have been proposed for the temperature in molecular simulations, where some of them have configurational contributions. We investigate how their accuracy is influenced by the number of particles in the simulation and the discontinuity in the derivatives of the interaction potential introduced by truncation. For equilibrium molecular dynamics with fixed total volume and fixed average total energy per particle, all the evaluated expressions including that for the kinetic temperature give a dependence on the total number of particles in the simulation. However, in a partitioned simulation volume under the same conditions, the mean temperature of each bin is independent of the number of bins. This finding is important for consistently defining a local temperature for use in nonequilibrium simulations. We identify the configurational temperature expressions which agree most with the kinetic temperature and find that they give close to identical results in nonequilibrium molecular dynamics (NEMD) simulations with a temperature gradient, for high and low density bulk-systems (both for transient and steady-state conditions), and across vapor-liquid interfaces, both at equilibrium and during NEMD simulations. The work shows that the configurational temperature is equivalent to the kinetic temperature in steady-state molecular dynamics simulations if the discontinuity in the derivatives of the interaction potential is handled properly, by using a sufficiently long truncation-distance or tail-corrections.

Sorting particles with nanoscale thermophoretic devices: how efficient is it?

We investigate particle separation driven by thermal gradients across solid state nanopores using a combined molecular dynamics simulation, non-equilibrium thermodynamics theory and a kinetic model approach. The thermophoretic device, a thermal nanopump, exploits thermal gradients to sort particles of different mass, which accumulate preferentially in hot or cold reservoirs. We show that the large amount of energy dissipated by the thermal nanopump during the transport process leads in general to very low efficiencies, 0.01-0.15%. We find that the nanopump thermal conductivity and structure plays a crucial role in determining the efficiency and a route to enhance it. Doubling the pore radius, from 0.5-1 nm radius, leads to a large increase in the mass diffusion and to a 20 fold increase in the efficiency. Addition of nanoscale defects, without modification of the nanopore structure, leads to a large reduction of the nanopump thermal conductivity and to a large enhancement of the thermodynamic efficiency. We find that nanopumps with nanoscale defects are >3 times more efficient than those without defects. Finally, we identify the microscopic variables responsible for the enhancement of thermally induced transport across nanopores and discuss strategies to tune these variable in order to regulate transport efficiency.

Active transport of the Ca(2+)-pump: introduction of the temperature difference as a driving force.

We analyse a kinetic cycle of the Ca(2+)-ATPase molecular pump using mesoscopic non-equilibrium thermodynamics. The pump is known to generate heat, and by analysing the operation on the mesoscopic level, we are able to introduce a temperature difference and the corresponding heat flux in the description. Integration over the internal coordinates then results in non-linear flux-force relations describing the operation of the pump on the macroscopic level. Specifically, we obtain an expression for the heat flux associated with the active transport and the coupling of heat effects to the transport of ions and the rate of the ATP-hydrolysis.

Computational modeling of the molecular basis for the calcium-dependence of the mannuronan C-5 epimerase AvAlgE6 from Azotobacter vinelandii.

The mannuronan C-5 epimerases catalyze epimerization of ß-d-mannuronic acid to α-l-guluronic acid in alginate polymers. The seven extracellular Azotobacter vinelandii epimerases (AvAlgE1-7) are calcium-dependent, and calcium is essential for the structural integrity of their carbohydrate binding R-modules. Ca2+ is also found in the crystal structures of the A-modules, where it is suggested to play a structural role. In this study, the structure of the catalytic A-module of the A. vinelandii mannuronan C-5 epimerase AvAlgE6 is used to investigate the role of this Ca2+. Molecular dynamics (MD) simulations with and without calcium reveal the possible importance of the bound Ca2+ in the hydrophobic packing of ß-sheets. In addition, a putative calcium binding site is found in the active site, indicating a potential direct role of this calcium in the catalysis. According to the literature, two of the residues coordinating calcium in this site are essential for the activity. MD simulations of the interaction with bound substrate indicate that the presence of a calcium ion in this binding site increases the binding strength. Further, explicit calculations of the substrate dissociation pathways with umbrella sampling simulations show and energetically higher dissociation barrier when calcium is present. The present study eludes to a putative catalytic role of calcium in the charge neutralizing first step of the enzymatic reaction. In addition to the importance for understanding these enzymes' molecular mechanisms, this could have implications for engineering strategies of the epimerases in industrial alginate processing.

On the thermodynamic efficiency of Ca²âº-ATPase molecular machines.

Experimental studies have shown that the activity of the reconstituted molecular pump Ca(2+)-ATPase strongly depends on the thickness of the supporting bilayer. It is thus expected that the bilayer structure will have an impact on the thermodynamic efficiency of this nanomachine. Here, we introduce a nonequilibrium-thermodynamics theoretical approach to estimate the thermodynamic efficiency of the Ca(2+)-ATPase from analysis of available experimental data about ATP hydrolysis and Ca(2+) transport. We find that the entropy production, i.e., the heat released to the surroundings under working conditions, is approximately constant for bilayers containing phospholipids with hydrocarbon chains of 18-22 carbon atoms. Our estimates for the heat released during the pump operation agree with results obtained from separate calorimetric experiments on the Ca(2+)-ATPase derived from sarcoplasmic reticulum. We show that the thermodynamic efficiency of the reconstituted Ca(2+)-ATPase reaches a maximum for bilayer thicknesses corresponding to maximum activity. Surprisingly, the estimated thermodynamic efficiency is very low, ∼12%. We discuss the significance of this result as representative of the efficiency of other nanomachines, and we address the influence of the experimental set-up on such a low efficiency. Overall, our approach provides a general route to estimate thermodynamic efficiencies and heat dissipation in experimental studies of nanomachines.

Kinetic and mesoscopic non-equilibrium description of the Ca(2+) pump: a comparison.

We analyse the operation of the Ca(2+)-ATPase ion pump using a kinetic cycle diagram. Using the methodology of Hill, we obtain the cycle fluxes, entropy production and efficiency of the pump. We compare these results with a mesoscopic non-equilibrium description of the pump and show that the kinetic and mesoscopic pictures are in accordance with each other. This gives further support to the mesoscopic theory, which is less restricted and also can include the heat flux as a variable. We also show how motors can be characterised in terms of unidirectional backward fluxes. We proceed to show how the mesoscopic approach can be used to identify fast and slow steps of the model in terms of activation energies, and how this can be used to simplify the kinetic diagram.