Source author record

Mathieu Salanne

Mathieu Salanne appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2021arXiv

Multi-scale simulation of the adsorption of lithium ion on graphite surface: from Quantum Monte Carlo to Molecular Density Functional Theory

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.

preprint2019arXiv

A first-principles investigation of the structural and electrochemical properties of biredox ionic species in acetonitrile

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.

preprint2016arXiv

Structural Transitions at Ionic Liquids Interfaces

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.

preprint2015arXiv

Polarization effects in ionic solids and melts

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.

preprint2015arXiv

Simulations of room temperature ionic liquids: From polarizable to coarse-grained force fields

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.

preprint2015arXiv

Sparse Cyclic Excitations Explain the Low Ionic Conductivity of Stoichiometric Li$_{7}$La$_{3}$Zr$_{2}$O$_{12}$

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.

preprint2015arXiv

Tuning the Electronic Structure of Anatase Through Fluorination

A highly fluorinated anatase lattice has been recently reported, providing a new class of materials whose general chemical formula is $\rm Ti_{1- \mathit x}\square_{\mathit x}X_{4\mathit x}O_{2- 4\mathit x}$ (X$^-$ = F$^-$ or OH$^-$). To characterise the complex structural features of the material and the different F environments, we here apply a computational screening procedure. After deriving a polarisable force--field from DFT simulations, we screen in a step-wise fashion a large number of possible configurations differing in the positioning of the titanium vacancies ($\square$) and of the fluorine atoms. At each step only 10% of the configurations are retained. At the end of the screening procedure, a configuration is selected and simulated using DFT-based molecular dynamics. This allows us to analyse the atomic structure of the material, which is strongly disordered, leading to a strong decrease (by 0.8~eV) of the band gap compared to conventional anatase.

preprint2014arXiv

Ceria co-doping: Synergistic or average effect?

Ceria (CeO$_{2}$) co-doping has been suggested as a means to achieve ionic conductivities that are significantly higher than those in singly-doped systems. Rekindled interest in this topic over the last decade has given rise to claims of much improved performance. The present study makes use of computer simulations to investigate the bulk ionic conductivity of Rare Earth (RE) doped ceria, where RE = Sc, Gd, Sm, Nd and La. The results from the singly doped systems are compared to those from ceria co-doped with Nd/Sm and Sc/La. The pattern that emerges from the conductivity data is consistent with the dominance of local lattice strains from individual defects, rather than the synergistic co-doping effect reported recently and, as a result no enhancement in the conductivity of co-doped samples is observed.

preprint2014arXiv

Effect of dispersion interactions on the properties of LiF in condensed phases

Classical molecular dynamics simulations are performed on LiF in the framework of the polarizable ion model. The overlap-repulsion and polarization terms of the interaction potential are derived on a purely non empirical, first-principles basis. For the dispersion, three cases are considered: a first one in which the dispersion parameters are set to zero and two others in which they are included, with different parameterizations. Various thermodynamic, structural and dynamic properties are calculated for the solid and liquid phases. The melting temperature is also obtained by direct coexistence simulations of the liquid and solid phases. Dispersion interactions appear to have an important effect on the density of both phases and on the melting point, although the liquid properties are not affected when simulations are performed in the NVT ensemble at the experimental density.

preprint2014arXiv

The electric double layer has a life of its own

Using molecular dynamics simulations with recently developed importance sampling methods, we show that the differential capacitance of a model ionic liquid based double-layer capacitor exhibits an anomalous dependence on the applied electrical potential. Such behavior is qualitatively incompatible with standard mean-field theories of the electrical double layer, but is consistent with observations made in experiment. The anomalous response results from structural changes induced in the interfacial region of the ionic liquid as it develops a charge density to screen the charge induced on the electrode surface. These structural changes are strongly influenced by the out-of-plane layering of the electrolyte and are multifaceted, including an abrupt local ordering of the ions adsorbed in the plane of the electrode surface, reorientation of molecular ions, and the spontaneous exchange of ions between different layers of the electrolyte close to the electrode surface. The local ordering exhibits signatures of a first-order phase transition, which would indicate a singular charge-density transition in a macroscopic limit.

preprint2014arXiv

Thermal Conductivity of Simple Liquids:Origin of Temperature and Packing Fraction Dependences

The origin of both weak temperature dependence and packing fraction dependence of $T^{1/4}η^{3/2}$ in the thermal conductivity of the simple Lennard-Jones (LJ) liquid is explored. In order to discuss the relative contributions from attractive or repulsive part of the interaction potential separately, the thermal conductivity of a series of Weeks-Chandler-Anderson (WCA) fluids is calculated by molecular dynamics simulations. The results show that the repulsive part plays the main role in the heat conduction, while the attractive part has no direct effect on the thermal conductivity for a given packing fraction. By investigating WCA fluids with potentials of varying softness, we explain the difference observed between the LJ liquids such as argon and Coulombic liquids such as NaCl.

preprint2013arXiv

Charge fluctuations in nano-scale capacitors

The fluctuations of the charge on an electrode contain information on the microscopic correlations within the adjacent fluid and their effect on the electronic properties of the interface. We investigate these fluctuations using molecular dynamics simulations in a constant-potential ensemble with histogram reweighting techniques. This approach offers in particular an efficient, accurate and physically insightful route to the differential capacitance that is broadly applicable. We demonstrate these methods with three different capacitors: pure water between platinum electrodes, and a pure as well as a solvent-based organic electrolyte each between graphite electrodes. The total charge distributions with the pure solvent and solvent-based electrolytes are remarkably Gaussian, while in the pure ionic liquid the total charge distribution displays distinct non-Gaussian features, suggesting significant potential-driven changes in the organization of the interfacial fluid.

preprint2013arXiv

From molten salts to room temperature ionic liquids: Simulation studies on chloroaluminate systems

An interaction potential including chloride anion polarization effects, constructed from first-principles calculations, is used to examine the structure and transport properties of a series of chloroaluminate melts. A particular emphasis was given to the study of the equimolar mixture of aluminium chloride with 1-ethyl-3-methylimidazolium chloride, which forms a room temperature ionic liquid EMI-AlCl 4. The structure yielded by the classical simulations performed within the framework of the polarizable ion model is compared to the results obtained from entirely electronic structure-based simulations: An excellent agreement between the two flavors of molecular dynamics is observed. When changing the organic cation EMI+ by an inorganic cation with a smaller ionic radius (Li+, Na+, K+), the chloroaluminate speciation becomes more complex, with the formation of Al2Cl 7- in small amounts. The calculated transport properties (diffusion coefficients, electrical conductivity and viscosity) of EMI-AlCl4 are in good agreement with experimental data.

preprint2013arXiv

Influence of solvation on the structural and capacitive properties of electrical double layer capacitors

We use molecular dynamics simulations to explore the impact of a non-ionic solvent on the structural and capacitive properties of supercapacitors based on an ionic liquid electrolyte and carbon electrodes. The study is focused on two pure ionic liquids, namely 1-butyl-3-methylimidazolium hexafluorophosphate and 1-butyl-3-methylimidazolium tetrafluoroborate, and their 1.5 M solutions in acetonitrile. The electrolytes, represented by coarse-grained models, are enclosed between graphite electrodes. We employ a constant potential methodology which allows us to gain insight into the influence of solvation on the polarization of the electrodes as well as the structural and capacitive properties of the electrolytes at the interface. We show that the interfacial characteristics, different for two distinct pure ionic liquids, become very similar upon mixing with acetonitrile.

preprint2013arXiv

Structure and dynamics in yttrium-based molten rare earth alkali fluorides

The transport properties of molten LiF-YF3 mixtures have been studied by pulsed field gradient nuclear magnetic resonance spectroscopy, potentiometric experiments, and molecular dynamics simulations. The calculated diffusion coefficients and electric conductivities compare very well with the measurements accross a wide composition range. We then extract static (radial distribution functions, coordination numbers distributions) and dynamic (cage correlation functions) quantities from the simulations. Then, we discuss the interplay between the microscopic structure of the molten salts and their dynamic properties. It is often considered that variations in the diffusion coefficient of the anions are mainly driven by the evolution of its coordination with the metallic ion (Y3+ here). We compare this system with fluorozirconate melts and demonstrate that the coordination number is a poor indicator of the evolution of the diffusion coefficient. Instead, we propose to use the ionic bonds lifetime. We show that the weak Y-F ionic bonds in LiF-YF3 do not induce the expected tendency of the fluoride diffusion coefficient to converge toward the one of yttrium cation when the content in YF3 increases. Implications on the validity of the Nernst-Einstein relation for estimating the electrical conductivity are discussed.

preprint2012arXiv

A transferable ab-initio based force field for aqueous ions

We present a new polarizable force field for aqueous ions (Li+, Na+, K+, Rb+, Cs+, Mg2+, Ca2+, Sr2+ and Cl-) derived from condensed phase ab-initio calculations. We use Maximally Localized Wannier Functions together with a generalized force and dipole-matching procedure to determine the whole set of parameters. Experimental data is then used only for validation purposes and a good agreement is obtained for structural, dynamic and thermodynamic properties. The same procedure applied to crystalline phases allows to parametrize the interaction between cations and the chloride anion. Finally, we illustrate the good transferability of the force field to other thermodynamic conditions by investigating concentrated solutions.

preprint2012arXiv

Diffusion coefficient and shear viscosity of rigid water models

We report the diffusion coefficient and viscosity of popular rigid water models: Two non polarizable ones (SPC/E with 3 sites, and TIP4P/2005 with 4 sites) and a polarizable one (Dang-Chang, 4 sites). We exploit the dependence of the diffusion coefficient on the system size [Yeh and Hummer, J. Phys. Chem. B 108, 15873 (2004)] to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang water overestimate the diffusion coefficient and underestimate the viscosity.

preprint2012arXiv

Including many-body effects in models for ionic liquids

Realistic modeling of ionic systems necessitates taking explicitly account of many-body effects. In molecular dynamics simulations, it is possible to introduce explicitly these effects through the use of additional degrees of freedom. Here we present two models: The first one only includes dipole polarization effect, while the second also accounts for quadrupole polarization as well as the effects of compression and deformation of an ion by its immediate coordination environment. All the parameters involved in these models are extracted from first-principles density functional theory calculations. This step is routinely done through an extended force-matching procedure, which has proven to be very succesfull for molten oxides and molten fluorides. Recent developments based on the use of localized orbitals can be used to complement the force-matching procedure by allowing for the direct calculations of several parameters such as the individual polarizabilities.

preprint2012arXiv

Thermal conductivity of MgO, MgSiO3 perovskite and post-perovskite in the Earth's deep mantle

We report lattice thermal conductivities of MgO and MgSiO3 in the perovskite and post-perovskite structures at conditions of the Earth's lower mantle, obtained from equilibrium molecular dynamics simulations. Using an advanced ionic interaction potential, the full conductivity tensor was calculated by means of the Green-Kubo method, and the conductivity of MgSiO3 post-perovskite was found to be significantly anisotropic. The thermal conductivities of all three phases were parameterized as a function of density and temperature. Assuming a Fe-free lower-mantle composition with mole fractions xMgSiO3 = 0.66 and xMgO = 0.34, the conductivity of the two-phase aggregate was calculated along a model geotherm. It was found to vary considerably with depth, rising from 9.5 W/(mK) at the top of the lower mantle to 20.5 W/(mK) at the top of the thermal boundary layer above the core-mantle boundary. Extrapolation of experimental data suggests that at deep-mantle conditions, the presence of a realistic amount of iron impurities lowers the thermal conductivity of the aggregate by about 50% (Manthilake et al., 2011). From this result and our thermal conductivity model, we estimate the heat flux across the core-mantle boundary to be 10.8 TW for a Fe-bearing MgO/MgSiO3 perovskite aggregate and 10.6 TW for a Fe-bearing MgO/MgSiO3 post-perovskite aggregate.

preprint2010arXiv

High-pressure behaviour of GeO2: a simulation study

In this work we study the high pressure behaviour of liquid and glassy GeO2 by means of molecular dynamics simulations. The interaction potential, which includes dipole polarization effects, was parameterized from first-principles calculations. Our simulations reproduce the most recent experimental data to a high degree of precision. The proportion of the various GeOn polyhedra is determined as a function of the pressure: a smooth transition from tetrahedral to octahedral network is observed. Finally, the study of high-pressure, liquid germania confirms that this material presents an anomalous behaviour of the diffusivity as observed in analog systems such as silica and water. The importance of penta-coordinated germanium ions for such behaviour is stressed.