Search

A crucial aspect in the simulation of electrochemical interfaces consists in treating the distribution of electronic charge of electrode materials that are put in contact with an electrolyte solution. Recently, it has been shown how a machine-learning method that specifically targets the electronic charge density, also known as SALTED, can be used to predict the long-range response of metal electrodes in model electrochemical cells. In this work, we provide a full integration of SALTED with MetalWalls, a program for performing classical simulations of electrochemical systems. We do so by deriving a spherical harmonics extension of the Ewald summation method, which allows us to efficiently compute the electric field originated by the predicted electrode charge distribution. We show how to use this method to drive the molecular dynamics of an aqueous electrolyte solution under the quantum electric field of a gold electrode, which is matched to the accuracy of density-functional theory. Notably, we find that the resulting atomic forces present a small error of the order of 1 meV/{\AA}, demonstrating the great effectiveness of adopting an electron-density path in predicting the electrostatics of the system. Upon running the data-driven dynamics over about 3 ns, we observe qualitative differences in the interfacial distribution of the electrolyte with respect to the results of a classical simulation. By greatly accelerating quantum-mechanics/molecular-mechanics approaches applied to electrochemical systems, our method opens the door to nanoseconds timescales in the accurate atomistic description of the electrical double layer., Comment: 10 pages, 6 figures

Molecular dynamics simulations at a constant electric potential are an essential tool to study electrochemical processes, providing microscopic information on the structural, thermodynamic, and dynamical properties. Despite the numerous advances in the simulation of electrodes, they fail to accurately represent the electronic structure of materials such as graphite. In this work, we introduce a simple parameterization method that allows to tune the metallicity of the electrode based on a quantum chemistry calculation of the density of states. As a first illustration, we study the interface between graphite electrodes and two different liquid electrolytes, an aqueous solution of NaCl and a pure ionic liquid, at different applied potentials. We show that the simulations reproduce qualitatively the experimentally-measured capacitance; in particular, they yield a minimum of capacitance at the point of zero charge, which is due to the quantum capacitance contribution. An analysis of the structure of the adsorbed liquids allows to understand why the ionic liquid displays a lower capacitance despite its large ionic concentration. In addition to its relevance for the important class of carbonaceous electrodes, this method can be applied to any electrode materials (e.g. 2D materials, conducting polymers, etc), thus enabling molecular simulation studies of complex electrochemical devices in the future.

Molecular dynamics simulations have been used extensively to determine equilibrium properties of the electrode-electrolyte interface in supercapacitors held at various potentials. While such studies are essential to understand and optimize the performance of such energy storage systems, investigations of the dynamics of adsorption during the charge of the supercapacitors is also necessary. Dynamical properties are especially important to get an insight into the power density of supercapacitors, one of their main assets. In this work, we propose a new method to analyze the trajectories of adsorbing ions. We focus on pure 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide in contact with planar carbon electrodes. We characterize the evolution of the ion orientation and ion-electrode distance during adsorption and show that ions reorientate as they adsorb. We then determine the forces experienced by the adsorbing ions and demonstrate that Coulomb forces dominate at long range while van der Waals forces dominate at short range. We also show that there is an almost equal contribution from the two forces at an intermediate distance, explaining the peak of ion density close to the electrode surface.

Identifying optimal collective variables to model transformations, using atomic-scale simulations, is a long-standing challenge. We propose a new method for the generation, optimization, and comparison of collective variables, which can be thought of as a data-driven generalization of the path collective variable concept. It consists in a kernel ridge regression of the committor probability, which encodes a transformation's progress. The resulting collective variable is one-dimensional, interpretable, and differentiable, making it appropriate for enhanced sampling simulations requiring biasing. We demonstrate the validity of the method on two different applications: a precipitation model, and the association of Li$^+$ and F$^-$ in water. For the former, we show that global descriptors such as the permutation invariant vector allow to reach an accuracy far from the one achieved \textit{via} simpler, more intuitive variables. For the latter, we show that information correlated with the transformation mechanism is contained in the first solvation shell only, and that inertial effects prevent the derivation of optimal collective variables from the atomic positions only.

The computational study of energy storage and conversion processes calls for simulation techniques that can reproduce the electronic response of metal electrodes under electric fields. Despite recent advancements in machine-learning methods applied to electronic-structure properties, predicting the non-local behavior of the charge density in electronic conductors remains a major open challenge. We combine long-range and equivariant kernel methods to predict the Kohn-Sham electron density of metal electrodes in response to various kinds of electric field perturbations. By taking slabs of gold as an example, we first show how the non-local electronic polarization generated by the interaction with an ionic species can be accurately reproduced in electrodes of arbitrary thickness. A finite-field extension of the method is then introduced, which allows us to predict the charge transfer and the electrostatic potential drop induced by the application of a homogeneous and constant electric field. Finally, we demonstrate the capability of the method to reproduce the charge-density response in a gold/electrolyte capacitor under an applied voltage, predicting the system polarization with a greater accuracy than state-of-the-art classical atomic-charge models., Comment: 6 pages, 3 figures

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic Density Functional Theory and Molecular Dynamics. Here we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference Quantum Monte Carlo calculations that allow us to benchmark various electronic Density Functional Theory functionals. We then fit an accurate carbon--lithium pair potential, which is used in molecular Density Functional Theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer., Comment: 10 pages 4 figures

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

In recent years, constant applied potential molecular dynamics has allowed to study the structure and dynamics of the electrochemical double-layer of a large variety of nanoscale capacitors. Nevertheless it remained impossible to simulate polarized electrodes at fixed total charge. Here we show that combining a constant potential electrode with a finite electric displacement fills this gap by allowing to simulate open circuit conditions. The method can be extended by applying an electric displacement ramp to perform computational amperometry experiments at different current intensities. As in experiments, the full capacitance of the system is obtained at low intensity, but this quantity decreases when the applied ramp becomes too fast with respect to the microscopic dynamics of the liquid., Comment: 17 pages, 7 figures

Hydrophobic hydration at metal/water interfaces actively contributes to the energetics of electrochemical reactions, e.g. CO$_2$ and N$_2$ reduction, where small hydrophobic molecules are involved. In this work, constant applied potential molecular dynamics is employed to study hydrophobic hydration at a gold/water interface. We propose an extension of the Lum-Chandler-Weeks (LCW) theory to describe the free energy of hydrophobic hydration at the interface as a function of solute size and applied voltage. Based on this model we are able to predict the free energy cost of cavity formation at the interface directly from the free energy cost in the bulk plus an interface-dependent correction term. The interfacial water network contributes significantly to the free energy yielding a preference for outer-sphere adsorption at the gold surface for ideal hydrophobes. We predict an accumulation of small hydrophobic solutes of sizes comparable to CO or N$_2$, while the free energy cost to hydrate larger hydrophobes, above 2.5 Angstrom radius, is shown to be greater at the interface than in the bulk. Interestingly, the transition from the volume dominated to the surface dominated regimes predicted by the LCW theory in the bulk is also found to take place for hydrophobes at the Au/water interface, but occurs at smaller cavity radii. By applying the extended LCW theory to a simple model addition reaction, we illustrate some implications of our findings for electrochemical reactions.

Soda-lime-silica is a glassy system of strong industrial interest. In order to characterize its liquid state properties, we performed molecular dynamics simulations employing an aspherical ion model that includes atomic polarization and deformation effects. They allowed to study the structure and diffusion properties of the system at temperatures ranging from 1400 to 3000 K. We show that Na$^+$ and Ca$^{2+}$ ions adopt a different structural organization within the silica network, with Ca$^{2+}$ ions having a greater affinity for non-bridging oxygens than Na$^+$. We further link this structural behavior to their different diffusivities, suggesting that escaping from the first oxygen coordination shell is the limiting step for the diffusion. Na$^+$ diffuses faster than Ca$^{2+}$ because it is bonded to a smaller number of non-bridging oxygens. The formed ionic bonds are also less strong in the case of Na$^+$., Comment: 9 pages

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.

Biredox ionic liquids are a new class of functionalized electrolytes that may play an important role in future capacitive energy storage devices. By allowing additional storage of electrons inside the liquids, they can improve device performance significantly. However current devices employ nanoporous carbons in which the diffusion of the liquid and the adsorption of the ions could be affected by the occurrence of electron-transfer reactions. It is therefore necessary to understand better the thermodynamics and the kinetics of such reactions in biredox ionic liquids. Here we perform ab initio molecular dynamics simulations of both the oxidized and reduced species of several redox-active ionic molecules (used in biredox ionic liquids) dissolved in acetonitrile solvent and compare them with the bare redox molecules. We show that in all the cases, it is necessary to introduce a two Gaussian state model to calculate the reaction free energies accurately. These reaction free energies are only slightly affected by the presence of the IL group on the molecule. We characterize the structure of the solvation shell around the redox active part of the molecules and show that in the case of TEMPO-based molecules strong reorientation effects occur during the oxidation reaction., Comment: 11 pages, 12 figures

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)

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. While 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 charges 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., Comment: 6 pages, 7 figures (including SI)

Most of the performances of electrochemical devices are governed by molecular processes taking place at the solution-electrode interfaces and molecular simulation are the main way to study these processes. Aqueous electrochemical systems have often been studied using classical DFT but with too crude approximations to consider the system description to be realistic. We study the interface between graphene electrodes and liquid water at different applied voltage using molecular DFT, improving the state of the art by the following key points: 1) electrodes have a realistic atomic resolution, 2) classical DFT calculations are carried out at fixed imposed potential difference and 3) water is described by a molecular model. This allows to reveal the structural modification of water adsorbed at the graphene interface and the evolution of water dielectric permittivity when a voltage is applied. The computed capacitance of this device is in agreement with molecular dynamics simulations. This demonstrates the relevance of molecular DFT to study electrochemical systems at the molecular level., Comment: 6 pages, 5 figures, Supplementary infos

Supercapacitors are energy storage devices able to deliver electricity with a high power. They consist of porous carbon electrodes in a concentrated electrolyte. Charged is stored by the adsorption of ions at the electrode surface. Nanoporous carbons allow to store more electricity thanks to a charging mechanism that differs from carbons with larger pores. This mechanism has recently been uncovered by combining in situ experiments (nuclear magnetic resonance, electrochemical quartz microbalance) and molecular simulations., Comment: in French

Supercapacitors are electric devices able to deliver a large power, enabling their use e.g. for the recovery of breaking energy in cars. This is achieved by using two carbon electrodes and an electrolyte solution or a pure ionic liquid (Room Temperature Ionic Liquid, RTIL). Energy is stored by the adsorption of ions at the surface of the electrodes, but the microscopic mechanism underlying the exceptional performance of Carbide Derived Carbon (CDC) electrodes remained unknown. Using molecular simulation with realistic electrode structures and under constant voltage conditions, we investigate the effect of confinement and solvation on the microscopic charging mechanism. We further analyse the dynamics of the charging process and make the link with equivalent circuit models used by electrochemists., Comment: in French

Beyond the dielectric continuum description initiated by Marcus theory, the nowadays standard theoretical approach to study electron transfer (ET) reactions in solution or at interfaces is to use classical force field or ab initio Molecular Dynamics simulations. We propose here an alternative method based on liquid-state theory, namely molecular density functional theory, which is numerically much more efficient than simulations while still retaining the molecular nature of the solvent. We begin by reformulating molecular ET theory in a density functional language and show how to compute the various observables characterizing ET reactions from an ensemble of density functional minimizations. In particular, we define in that formulation the relevant order parameter of the reaction, the so-called vertical energy gap, and determine the Marcus free energy curves of both reactant and product states along that coordinate. Important thermodynamic quantities such as the reaction free energy and the reorganization free energies follow. We assess the validity of the method by studying the model Cl$^0\rightarrow$ Cl$^+$ and Cl$^0\rightarrow$ Cl$^-$ ET reactions in bulk water for which molecular dynamics results are available. The anionic case is found to violate the standard Marcus theory. Finally, we take advantage of the computational efficiency of the method to study the influence of confinement on the ET, by investigating the evolution of the reorganization free energy of the Cl$^0\rightarrow$ Cl$^+$ reaction when the atom approaches an atomistically resolved wall., Comment: 30 pages, 10 figures, Chemical Science, 2019

A lamellar lepidocrocite-type titanate structure with ~25% $\text{Ti}^{4+}$vacancies was recently synthesized, and it showed potential for use as an electrode in rechargeable lithium-ion batteries. In addition to lithium, we explore this material's ability to accommodate other monovalent ions with greater natural abundance (e.g. sodium and potassium) in order to develop lower-cost alternatives to lithium-ion batteries constructed from more widely available elements. Galvanostatic discharge/charge curves for the lepidocrocite material indicate that increasing the ionic radius of the monovalent ion results in a deteriorating performance of the electrode. Using first-principles electronic structure calculations, we identify the relaxed geometries of the structure for various positions of the ion in the structure. We then use these geometries to compute the energy of formations. Additionally, we determine that all ions are favorable in the structure, but interlayer positions are preferred compared to vacancy positions. We also conclude that the exchange between the interlayer and vacancy positions is a process which involves the interaction between interlayer water and surface hydroxyl groups next to the titanate layer. We observe a cooperative effect between structural water and $\text{OH}$ groups to assist alkali-ions to move from the interlayer to the vacancy site. Thus, the as-synthesized lepidocrocite serves as a prototypical structure to investigate both the migration mechanism of ions within a confined space along with the interaction between water molecules and the titanate framework.

When an ionic liquid adsorbs onto a porous electrode, its ionic arrangement is deeply modified due to a screening of the Coulombic interactions by the metallic surface and by the confinement imposed upon it by the electrode's morphology. In particular, ions of the same charge can approach at close contact, leading to the formation of a superionic state. The impact of an electrified surface placed between two liquid phases is much less understood. Here we simulate a full supercapacitor made of the 1-butyl-3-methylimidazolium hexafluorophosphate and nanoporous graphene electrodes, with varying distances between the graphene sheets. The electrodes are held at constant potential by allowing the carbon charges to fluctuate. Under strong confinement conditions, we show that ions of the same charge tend to adsorb in front of each other across the graphene plane. These correlations are allowed by the formation of a highly localized image charge on the carbon atoms between the ions. They are suppressed in larger pores, when the liquid adopts a bilayer structure between the graphene sheets. These effects are qualitatively similar to the recent templating effects which have been reported during the growth of nanocrystals on a graphene substrate.

Nanoporous carbon-based supercapacitors store electricity through adsorption of ions from the electrolyte at the surface of the electrodes. Room temperature ionic liquids, which show the largest ion concentrations among organic liquid electrolytes, should in principle yield larger capacitances. Here we show by using electrochemical measurements that the capacitance is not significantly affected when switching from a pure ionic liquid to a conventional organic electrolyte using the same ionic species. By performing additional molecular dynamics simulations, we interpret this result as an increasing difficulty of separating ions of opposite charges when they are more concentrated, i.e. in the absence of a solvent which screens the Coulombic interactions. The charging mechanism consistently changes with ion concentration, switching from counter-ion adsorption in the diluted organic electrolyte to ion exchange in the pure ionic liquid. Contrarily to the capacitance, in-pore diffusion coefficients largely depend on the composition, with a noticeable slowing of the dynamics in the pure ionic liquid., Comment: 8 pages, 6 figures

Nanoconfinement generally leads to drastic effect on the physical and chemical properties of ionic liquids. Here we investigate how the electrochemical reactivity in such media may be impacted inside nanoporous carbon electrodes. To this end, we study a simple electron transfer reaction using molecular dynamics simulations. The electrodes are held at constant electric potential by allowing the atomic charges on the carbon atoms to fluctuate. We show that the $\mathrm{Fe^{3+}}/\mathrm{Fe^{2+}}$ couple dissolved in an ionic liquid exhibits a deviation with respect to Marcus theory. This behavior is rationalized by the stabilization of a solvation state of the Fe$^{3+}$ cation in the disordered nanoporous electrode that is not observed in the bulk. The simulation results are fitted with a recently proposed two solvation state model, which allows us to estimate the effect of such a deviation on the kinetics of electron transfer inside nanoporous electrodes., Comment: 8 pages, 10 figures

Despite the simplicity of its molecular unit, water is a challenging system because of its uniquely rich polymorphism and predicted but yet unconfirmed features. Introducing a novel space of generalized coordinates that capture changes in the topology of the interatomic network, we are able to systematically track transitions among liquid, amorphous and crystalline forms throughout the whole phase diagram of water, including the nucleation of crystals above and below the melting point. Our approach, based on molecular dynamics and enhanced sampling / free energy calculation techniques, is not specific to water and could be applied to very different structural phase transitions, paving the way towards the prediction of kinetic routes connecting polymorphic structures in a range of materials.

We perform molecular dynamics simulations of a typical nanoporous-carbon based supercapacitors. The organic electrolyte consists in 1-ethyl-3-methyl--imidazolium and hexafluorophosphate ions dissolved in acetonitrile. We simulate systems at equilibrium, for various applied voltages. This allows us to determine the relevant thermodynamic (capacitance) and transport (in-pore resistivities) properties. These quantities are then injected in a transmission line model for testing its ability to predict the charging properties of the device. The results from this macroscopic model are in good agreement with non-equilibrium molecular dynamics simulations, which validates its use for interpreting electrochemical impedance experiments., Comment: 7 pages, 6 figures

Recent advances in experimental and computational techniques have allowed for an accurate description of the adsorption of ionic liquids on metallic electrodes. It is now well established that they adopt a multi-layered structure, and that the composition of the layers changes with the potential of the electrode. In some cases, potential-driven ordering transitions in the first adsorbed layer have been observed in experiments probing the interface on the molecular scale or by molecular simulations. This perspective gives an overview of the current understanding of such transitions and of their potential impact on the physical and (electro)chemical processes at the interface. In particular, peaks in the differential capacitance, slow dynamics at the interface and changes in the reactivity have been reported in electrochemical studies. Interfaces between ionic liquids and metallic electrodes are also highly relevant for their friction properties, the voltage-dependence of which opens the way to exciting applications., Comment: 10 pages, 4 figures. arXiv admin note: text overlap with arXiv:1404.0343

We have performed long time-scale molecular dynamics simulations of the cubic and tetragonal phases of the solid lithium-ion--electrolyte Li$_{7}$La$_{3}$Zr$_{2}$O$_{12}$ (LLZO), using a first-principles parameterised interatomic potential. Collective lithium transport was analysed by identifying dynamical excitations; persistent ion displacements over distances comparable to the separation between lithium sites, and string-like clusters of ions that undergo cooperative motion. We find that dynamical excitations in c-LLZO are frequent, with participating lithium numbers following an exponential distribution, mirroring the dynamics of fragile glasses. In contrast, excitations in t-LLZO are both temporally and spatially sparse, consisting preferentially of highly concerted lithium motion around closed loops. This qualitative difference is explained as a consequence of lithium ordering in t-LLZO, and provides a mechanistic basis for the much lower ionic conductivity of t-LLZO compared to c-LLZO.

Room temperature ionic liquids (RTILs) are solvent with unusual properties, which are difficult to characterize experimentally because of their intrinsic complexity (large number of atoms, strong Coulomb interactions). Molecular simulations have therefore been essential in our understanding of these systems. Depending on the target property and on the necessity to account for fine details of the molecular structure of the ions, a large range of simulation techniques are available. Here I focus on classical molecular dynamics, in which the level of complexity of the simulation, and therefore the computational cost, mostly depends on the force field. Depending on the representation of the ions, these are either classified as all-atom or coarse-grained. In addition, all-atom force fields may account for polarization effects if necessary. The most widely used methods for RTILs are described together with their main achievements and limitations., Comment: 11 pages, 6 figures

Ionic solids and melts are compounds in which the interactions are dominated by electrostatic effects. However, the polarization of the ions also plays an important role in many respects as has been clarified in recent years thanks to the development of realistic polarizable interaction potentials. After detailing these models, we illustrate the importance of polarization effects on a series of examples concerning the structural properties, such as the stabilization of particular crystal structures or the formation of highly-coordinated multivalent ions in the melts, as well as the dynamic properties such as the diffusion of ionic species. The effects on the structure of molten salts interfaces (with vacuum and electrified metal) is also described. Although most of the results described here concern inorganic compounds (molten fluorides and chlorides, ionic oxides...), the particular case of the room-temperature ionic liquids, a special class of molten salts in which at least one species is organic, will also be briefly discussed to indicate how the ideas gained from the study of "simple" molten salts are being transferred to these more complex systems., Comment: 16 pages, 9 figures