Search

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 prominent role in this field since they provide a microscopic picture of the mechanisms involved. The current state-of-the-art consists of 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 differentiating 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. [ABSTRACT FROM AUTHOR]

Because magmatic liquids play a fundamental role in the evolution of the terrestrial planets, a precise knowledge of their physical properties is requisite to better understand the formation and the dynamics of planetary interiors. In using an improved force field for silicates (Dufils et al., 2017), we report the results of a series of molecular dynamics simulations (MD) aiming to evaluate the equation of state (EOS), the viscosity, the electrical conductivity, and the elemental self-diffusion coefficients of various planetary melts representative of the Earth, Mars, the Moon and Mercury. The agreement between MD calculations and experimental data (when they exist) is remarkable, a finding which suggests that the MD simulations can be used on trust to extend the existing laboratory data on planetary melts (as it is proposed here for some lunar, martian and mercurian basalts) or to predict the thermo-physical properties of more exotic compositions (magma oceans, lava planets). Moreover, the MD simulations show that the evolution of the viscosity and of the electrical conductivity with the pressure depends in a complex manner on temperature and composition and that is difficult to extract general trends from these behaviors. Another advantage of the MD simulations is that the transport coefficients (viscosity, conductivity, and diffusivities) being evaluated along the same numerical experiment, one is able to test the validity of the Eyring and the Stokes-Einstein equations relating viscosity and diffusivity, as also as the Nernst-Einstein equation expressing the conductivity as function of the ionic diffusivities. It is shown that the Stokes-Einstein equation is suitable to describe planetary melts of low viscosity (≤1 Pa·s), whereas the Eyring relation leads to a better estimation at high viscosities (>100 Pa·s). Concerning the electrical conductivity, the Nernst-Einstein relation fails to reproduce the conductivity values obtained by MD simulation, a result which can be explained by the crucial role played by ion-ion correlations in silicate melts. [ABSTRACT FROM AUTHOR]

A new atom–atom interaction potential is introduced for describing by classical molecular dynamics (MD) simulation the physical properties of natural silicate melts. The equation of state, the microscopic structure, the viscosity, the electrical conductivity, and the self-diffusion coefficients of ions in a mid-oceanic ridge basalt (MORB) melt are evaluated by MD over a large range of temperature and pressure (1673–3273 K and 0–60 GPa). A detailed comparison with experimental data shows that the model reproduces the thermodynamic, structural and transport properties of a MORB with an unprecedented accuracy. In particular, it is shown that the MORB melt crystallizes at lower mantle conditions into a perovskite phase whose the equation of state (EOS) is compatible with those proposed in the experimental literature. Moreover, in accordance with experimental findings, the simulation predicts not only that the MORB viscosity exhibits a (slight) minimum with the pressure, but also that the viscosity at high temperature remains very low (< 100 mPa.s for T > 2273 K) even at high pressure (up to 40 GPa). However the evolution of the electrical conductivity with temperature and pressure is not always the symmetrical of that of the viscosity. In fact, the relationship between viscosity and electrical conductivity shows a crossover at around 2073 K. [ABSTRACT FROM AUTHOR]

Despite its low abundance, water has a great influence on the geodynamics of the Earth's upper mantle. Indeed, water has the ability to modify the phase relations and to affect in a significant way the rheological properties of minerals and melts. However the mechanisms of water incorporation in silicate melts and the impact on the melt properties is still not fully understood. To improve our understanding of hydrous silicate melts, we have performed a series of molecular dynamics simulations to evaluate the H 2 O solubility, the liquid-vapour coexistence, the surface tension, the water speciation, the equation of state, the viscosity, the electrical conductivity, the diffusion of silicate elements and protonated species, as well as the melt structure of various magmatic liquids representative of the Earth's upper mantle (rhyolite, andesite, MORB, peridotite, and kimberlite). For that, we introduce a new force field for water, which is compatible with an accurate force field for silicates recently developed (Dufils et al., 2018). A comparison between MD calculations and experimental data (when they exist) shows that the MD simulations are reliable. Among all the results obtained in this study, the following points may be emphasized. (1) The solubility of water changes very little when the melt composition evolves from rhyolitic to andesitic and basaltic, but it is strongly enhanced in ultramafic melts. (2) When hydrous melt and aqueous fluid are coexisting with each other, the oxide content of the aqueous fluid increases rapidly with the pressure. (3) A consequence of point (2) is that water has a large influence on the surface tension, as the latter one drops by a factor of 2 ∼ 4 when the water pressure increases from 1 bar to a few kbar. (4) Concerning the water speciation, an important point is that the MD simulation probes the liquid phase, when most of the experimental studies are dealing with glasses. Thus at magmatic temperatures the concentration in hydroxyl groups and the one in molecular water are crossing for a water content of about 15 wt%, a value much higher than the one observed in glasses (∼ 3–4 wt%). (5) MD calculations show that the molar volume of the melt is a linear function of the water content, and so for all the chemical compositions investigated. Therefore the water partial molar volume (V H 2 O) is virtually independent of total water content and of water speciation. A by-product of this result is that an ideal mixing rule between water and the silicate component leads to an accurate estimate of the melt molar volume. (6) At fixed T and P, the melt viscosity decreases with water content, more depolymerized the melt the smaller the influence of water on the viscosity. However, at the high temperatures investigated in this study (T ≥ 1673 K), the decrease in viscosity induced by water does not exceed one or two orders of magnitude, as compared with many orders of magnitude near the glass transition temperature. (7) The diffusivity of ions increases exponentially with water content. As for the protonated species, it is found that, D O 2 - < D O H - < D H 2 O ≤ D H 3 O + , the lower the NBO/T ratio the smaller the ratio D O H - / D H 2 O. (8) A structural analysis shows that hydroxyl groups are more preferentially linked to metal cations than to structure makers. In contrast, H 2 O molecules (and H 3 O+ as well) are almost exclusively linked to metal cations. As for the melt polymerization, it decreases gradually with the water content in andesitic and basaltic melts, whereas it remains almost invariant in peridotitic melt. (9) O-H...O bonds (hydrogen bonding) taking place between the hydroxyl groups, the water molecules, and the oxygens of the silicate are characterized by O⋯O distances in the range 2.5 ∼ 3.2 A, and by O...H-O distances in the range 1.5 ∼ 2.2 A. But, because of the high temperature of investigation, these H-bonds are generally weak (weaker than in liquid water at ambient). [ABSTRACT FROM AUTHOR]

A better understanding of interfacial mechanisms is needed to improve the performances of electrochemical devices. Yet, simulating an electrode surface at fixed electrolyte composition remains a challenge. Here, we apply a finite electric field to a single electrode held at constant potential and in contact with an aqueous ionic solution, using classical molecular dynamics. The polarization yields two electrochemical interfaces on opposite sides of the same metal slab. The net charge on one electrode surface is the opposite of the net charge on the other, maintaining overall charge neutrality of the metal. The electrode surface charge fluctuations are compensated by the adsorption of ions from the electrolyte, forming a pair of electric double layers with aligned dipoles. This opens the way towards the efficient simulation of electrochemical interfaces using any flavor of molecular dynamics, from classical to first-principles-based methods. [ABSTRACT FROM AUTHOR]