Enhanced Sampling of Configuration and Path Space in a Generalized Ensemble by Shooting Point Exchange.

The computer simulation of many molecular processes is complicated by long timescales caused by rare transitions between long-lived states. Here, we propose a new approach to simulate such rare events, which combines transition path sampling with enhanced exploration of configuration space. The method relies on exchange moves between configuration and trajectory space, carried out based on a generalized ensemble. This scheme substantially enhances the efficiency of the transition path sampling simulations, particularly for systems with multiple transition channels, and yields information on thermodynamics, kinetics and reaction coordinates of molecular processes without distorting their dynamics. The method is illustrated using the isomerization of proline in the KPTP tetrapeptide.

MetalWalls: Simulating electrochemical interfaces between polarizable electrolytes and metallic electrodes.

Electrochemistry is central to many applications, ranging from biology to energy science. Studies now involve a wide range of techniques, both experimental and theoretical. Modeling and simulations methods, such as density functional theory or molecular dynamics, provide key information on the structural and dynamic properties of the systems. Of particular importance are polarization effects of the electrode/electrolyte interface, which are difficult to simulate accurately. Here, we show how these electrostatic interactions are taken into account in the framework of the Ewald summation method. We discuss, in particular, the formal setup for calculations that enforce periodic boundary conditions in two directions, a geometry that more closely reflects the characteristics of typical electrolyte/electrode systems and presents some differences with respect to the more common case of periodic boundary conditions in three dimensions. These formal developments are implemented and tested in MetalWalls, a molecular dynamics software that captures the polarization of the electrolyte and allows the simulation of electrodes maintained at a constant potential. We also discuss the technical aspects involved in the calculation of two sets of coupled degrees of freedom, namely the induced dipoles and the electrode charges. We validate the implementation, first on simple systems, then on the well-known interface between graphite electrodes and a room-temperature ionic liquid. We finally illustrate the capabilities of MetalWalls by studying the adsorption of a complex functionalized electrolyte on a graphite electrode.

Erratum: Fluctuation relations for systems in a constant magnetic field [Phys. Rev. E 102, 030101(R) (2020)].

This corrects the article DOI: 10.1103/PhysRevE.102.030101.

Fluctuation Relations for Dissipative Systems in Constant External Magnetic Field: Theory and Molecular Dynamics Simulations.

We illustrate how, contrary to common belief, transient Fluctuation Relations (FRs) for systems in constant external magnetic field hold without the inversion of the field. Building on previous work providing generalized time-reversal symmetries for systems in parallel external magnetic and electric fields, we observe that the standard proof of these important nonequilibrium properties can be fully reinstated in the presence of net dissipation. This generalizes recent results for the FRs in orthogonal fields-an interesting but less commonly investigated geometry-and enables direct comparison with existing literature. We also present for the first time a numerical demonstration of the validity of the transient FRs with nonzero magnetic field via nonequilibrium molecular dynamics simulations of a realistic model of liquid NaCl.

Fluctuation relations for systems in a constant magnetic field.

The validity of the fluctuation relations (FRs) for systems in a constant magnetic field is investigated. Recently introduced time-reversal symmetries that hold in the presence of static electric and magnetic fields and of deterministic thermostats are used to prove the transient FRs without invoking, as commonly done, inversion of the magnetic field. Steady-state FRs are also derived, under the t-mixing condition. These results extend the predictive power of important statistical mechanics relations. We illustrate this via the nonlinear response for the cumulants of the dissipation, showing how the alternative FRs enable one to determine analytically null cumulants also for systems in a single magnetic field.

Adiabatic motion and statistical mechanics via mass-zero constrained dynamics.

In recent work [Coretti et al., J. Chem. Phys., 2018, 149, 191102], a new algorithm to solve numerically the dynamics of the shell model for polarization was presented. The approach, broadly applicable to systems involving adiabatically separated dynamical variables, employs constrained molecular dynamics to strictly enforce the condition that the force on the fast degrees of freedom, modeled as having zero mass, is null at each time step. The algorithm is symplectic and fully time reversible, and results in stable and efficient propagation. In this paper we complete the discussion of the mechanics of mass-zero constrained dynamics by showing how to adapt it to problems where the fast degrees of freedom must satisfy additional conditions. This extension includes, in particular, the important case of first principles molecular dynamics. We then consider the statistical mechanics of the mass-zero constrained dynamical system demonstrating that the marginal probability sampled by the dynamics in the physical phase space recovers the form of the Born-Oppenheimer probability density. The effectiveness of the approach and the favorable scaling of the algorithm with system size are illustrated in test calculations of solid Na via orbital-free density functional dynamics.

Charge fluctuations from molecular simulations in the constant-potential ensemble.

We revisit the statistical mechanics of charge fluctuations in capacitors. In constant-potential classical molecular simulations, the atomic charges of electrode atoms are treated as additional degrees of freedom which evolve in time so as to satisfy the constraint of fixed electrostatic potential for each configuration of the electrolyte. The present work clarifies the role of the overall electroneutrality constraint, as well as the link between the averages computed within the Born-Oppenheimer approximation and that of the full constant-potential ensemble. This allows us in particular to derive a complete fluctuation-dissipation relation for the differential capacitance, that includes a contribution from the charge fluctuations (around the charges satisfying the constant-potential and electroneutrality constraints) also present in the absence of an electrolyte. We provide a simple expression for this contribution from the elements of the inverse of the matrix defining the quadratic form of the fluctuating charges in the energy. We then illustrate numerically the validity of our results, and recover the expected continuum result for an empty capacitor with structureless electrodes at large inter-electrode distances. By considering a variety of liquids between graphite electrodes, we confirm that this contribution to the total differential capacitance is small compared to that induced by the thermal fluctuations of the electrolyte.

Communication: Constrained molecular dynamics for polarizable models.

A new algorithm to solve numerically the evolution of empirical shell models of polarizable systems is presented. It employs constrained molecular dynamics to satisfy exactly, at each time step, the crucial condition that the gradient of the potential with respect to the shell degrees of freedom is null. The algorithm is efficient, stable, and, contrary to the available alternatives, it is symplectic and time reversible. A proof-of-principle calculation on a polarizable model for NaCl is presented to illustrate its properties in comparison with the current method, which employs a conjugate-gradient procedure to enforce the null gradient condition. The proposed algorithm is applicable to other cases where a minimum condition on a function of an auxiliary set of driven dynamical variables must be satisfied.

Time-reversal symmetry for systems in a constant external magnetic field.

The time-reversal properties of charged systems in a constant external magnetic field are reconsidered in this paper. We show that the evolution equations of the system are invariant under a new symmetry operation that implies a new signature property for time-correlation functions under time reversal. We then show how these findings can be combined with a previously identified symmetry to determine, for example, null components of the correlation functions of velocities and currents and of the associated transport coefficients. These theoretical predictions are illustrated by molecular dynamics simulations of superionic AgI.