Search

The relation between the frequency-dependent friction of a molecule in a liquid and the hydrodynamic properties of the liquid is fundamental for molecular dynamics. We investigate this connection for a water molecule moving in liquid water using all-atomistic molecular dynamics simulations and linear hydrodynamic theory. We analytically calculate the frequency-dependent friction of a sphere with finite surface slip moving in a viscoelastic compressible fluid by solving the linear transient Stokes equation, including frequency-dependent shear and volume viscosities, both determined from MD simulations of bulk liquid water. We also determine the frequency-dependent friction of a single water molecule moving in liquid water, as defined by the generalized Langevin equation from MD simulation trajectories. By fitting the effective sphere radius and the slip length, the frequency-dependent friction and velocity autocorrelation function from the transient Stokes equation and simulations quantitatively agree. This shows that the transient Stokes equation accurately describes the important features of the frequency-dependent friction of a single water molecule in liquid water and thus applies down to molecular length and time scales, provided accurate frequency-dependent viscosities are used. The frequency dependence of the shear viscosity of liquid water requires careful consideration of hydrodynamic finite-size effects to observe the asymptotic hydrodynamic power-law tail. In contrast, for a methane molecule moving in water, the frequency-dependent friction cannot be predicted based on a homogeneous model, which suggests, supported by the extraction of a frequency-dependent surface-slip profile, that a methane molecule is surrounded by a finite-thickness hydration layer with viscoelastic properties that are significantly different from bulk water., Comment: 25 pages, 13 figures

Finite-size effects are challenging in molecular dynamics simulations because they have significant effects on computed static and dynamic properties, in particular diffusion constants, friction coefficients and time- or frequency-dependent response functions. We investigate the influence of periodic boundary conditions on the velocity autocorrelation function and the frequency-dependent friction of a particle in a fluid and show that the long-time behavior (starting at the picosecond timescale) is significantly affected. We develop an analytical correction allowing to subtract the periodic boundary condition effects. By this we unmask the power-law long-time tails of the memory kernel and the velocity autocorrelation function in liquid water and a Lennard-Jones fluid from rather small simulation boxes., Comment: 6 pages, 4 figures, 9 pages of Appendices

Understanding the response of the surface of metallic solids to external electric field sources is crucial to characterize electrode-electrolyte interfaces. Continuum electrostatics offer a simple description of the induced charge density at the electrode surface. However, such a simple description does not take into account features related to the atomic structure of the solid and to the molecular nature of the solvent and of the dissolved ions. In order to illustrate such effects and assess the ability of continuum electrostatics to describe the induced charge distribution, we investigate the behaviour of a gold electrode interacting with sodium or chloride ions fixed at various positions, in vacuum or in water, using all-atom constant-potential classical molecular dynamics simulations. Our analysis highlights important similarities between the two approaches, especially in vacuum conditions and when the ion is sufficiently far from the surface, as well as some limitations of the continuum description, namely neglecting the charges induced by the adsorbed solvent molecules and the screening effect of the solvent when the ion is close to the surface. While the detailed features of the charge distribution are system-specific, we expect some of our generic conclusions on the induced charge density to hold for other ions, solvents and electrode surfaces. Beyond this particular case, the present study also illustrates the relevance of such molecular simulations to serve as a reference for the design of improved implicit solvent models of electrode-electrolyte interfaces.

The nuclear magnetic resonance (NMR) relaxation of quadrupolar nuclei is governed by the electric field gradient (EFG) fluctuations at their position. In classical molecular dynamics (MD), the electron cloud contribution to the EFG can be included via the Sternheimer approximation, in which the full EFG at the nucleus that can be computed using quantum DFT is considered to be proportional to that arising from the external, classical charge distribution. In this work, we systematically assess the quality of the Sternheimer approximation as well as the impact of the classical force field (FF) on the NMR relaxation rates of aqueous quadrupolar ions at infinite dilution. In particular, we compare the rates obtained using an ab initio parametrized polarizable FF, a recently developed empirical FF with scaled ionic charges and a simple empirical non-polarizable FF with formal ionic charges. Surprisingly, all three FFs considered yield good values for the rates of smaller and less polarizable solutes, provided that a model-specific Sternheimer parametrization is employed. Yet, the polarizable and scaled charge FFs yield better estimates for divalent and more polarizable species. We find that a linear relationship between the quantum and classical EFGs holds well in all of the cases considered, however, such an approximation often leads to quite large errors in the resulting EFG variance, which is directly proportional to the computed rate. We attempted to reduce the errors by including first order nonlinear corrections to the EFG, yet no clear improvement for the resulting variance has been found. The latter result indicates that more refined methods for determining the EFG at the ion position, in particular those that take into account the instantaneous atomic environment around an ion, might be necessary to systematically improve the NMR relaxation rate estimates in classical MD., Comment: 50 pages, 15 figures, 5 tables, supporting information

Electrochemistry experiments have established that the capacitance of electrode-electrolyte interfaces is much larger for good metals such as gold and platinum than for carbon-based materials. Despite the development of elaborate electrode interaction potentials, to date molecular dynamics simulations were not able to capture this effect. Here we show that changing the width of the Gaussian charge distribution used to represent the atomic charges in gold is an effective way to tune its metallicity. Larger Gaussian widths lead to a capacitance of aqueous solutions (pure water and 1 molar NaCl) in good agreement with recent ab initio molecular dynamics results. For pure water, the increase in the capacitance is not accompanied with structural changes, while in the presence of salt the Na$^+$ cations tend to adsorb significantly on the surface. For a strongly metallic gold electrode, these ions can even form inner sphere complexes on hollow sites of the surface., Comment: 16 pages

We investigate the effect of the metallic character of solid substrates on solid-liquid interfacial thermodynamics using molecular simulations. Building on the recent development of a semi-classical Thomas-Fermi model to tune the metallicity in classical molecular dynamics simulations, we introduce a new thermodynamic integration framework to compute the evolution of the interfacial free energy as a function of the Thomas-Fermi screening length. We validate this approach against analytical results for empty capacitors and by comparing the predictions in the presence of an electrolyte with values determined from the contact angle of droplets on the surface. The general expression derived in this work highlights the role of the charge distribution within the metal. We further propose a simple model to interpret the evolution of the interfacial free energy with voltage and Thomas-Fermi length, which allows us to identify the charge correlations within the metal as the microscopic origin of the evolution of the interfacial free energy with the metallic character of the substrate. This new methodology opens the door to the molecular-scale study of the effect of the metallic character of the substrate on confinement-induced transitions in ionic systems, as reported in recent Atomic Force Microscopy and Surface Force Apparatus experiments., Comment: 12 pages, 9 figures including supporting information

The Gibbs-Thomson (GT) equation describes the shift of the crystallization temperature for a confined fluid with respect to the bulk as a function of pore size. While this century old relation is successfully used to analyze experiments, its derivations found in the literature often rely on nucleation theory arguments (i.e. kinetics instead of thermodynamics) or fail to state their assumptions, therefore leading to similar but different expressions. Here, we revisit the derivation of the GT equation to clarify the system definition, corresponding thermodynamic ensemble, and assumptions made along the way. We also discuss the role of the thermodynamic conditions in the external reservoir on the final result. We then turn to numerical simulations of a model system to compute independently the various terms entering in the GT equation, and compare the predictions of the latter with the melting temperatures determined under confinement by means of hyper-parallel tempering grand canonical Monte Carlo simulations. We highlight some difficulties related to the sampling of crystallization under confinement in simulations. Overall, despite its limitations, the GT equation may provide an interesting alternative route to predict the melting temperature in large pores, using molecular simulations to evaluate the relevant quantities entering in this equation. This approach could for example be used to investigate the nanoscale capillary freezing of ionic liquids recently observed experimentally between the tip of an Atomic Force Microscope and a substrate., Comment: 14 pages, 9 figures, contains appendix

Many key industrial processes, from electricity production, conversion and storage to electrocatalysis or electrochemistry in general, rely on physical mechanisms occurring at the interface between a metallic electrode and an electrolyte solution, summarized by the concept of electric double layer, with the accumulation/depletion of electrons on the metal side and of ions on the liquid side. While electrostatic interactions play an essential role on the structure, thermodynamics, dynamics and reactivity of electrode-electrolyte interfaces, these properties also crucially depend on the nature of the ions and solvent, as well as that of the metal itself. Such interfaces pose many challenges for modeling, because they are a place where Quantum Chemistry meets Statistical Physics. In the present review, we explore the recent advances on the description and understanding of electrode-electrolyte interfaces with classical molecular simulations, with a focus on planar interfaces and solvent-based liquids, from pure solvent to water-in-salt-electrolytes.

Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modelling their responses to the charge distribution of electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.

We revisit the statistical mechanics of charge fluctuations in capacitors. In constant-potential classical molecular simulations, the atomic charge 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 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.

Spurred by the increasing needs in electrochemical energy storage devices, the electrode/electrolyte interface has received a lot of interest in recent years. Molecular dynamics simulations play a proeminent role in this field since they provide a microscopic picture of the mechanisms involved. The current state-of-the-art consists in treating the electrode as a perfect conductor, precluding the possibility to analyze the effect of its metallicity on the interfacial properties. Here we show that the Thomas-Fermi model provides a very convenient framework to account for the screening of the electric field at the interface and differenciating good metals such as gold from imperfect conductors such as graphite. All the interfacial properties are modified by screening within the metal: the capacitance decreases significantly and both the structure and dynamics of the adsorbed electrolyte are affected. The proposed model opens the door for quantitative predictions of the capacitive properties of materials for energy storage., Comment: 10 pages, 10 figures (contains the SI)