Search

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, modeling 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. [ABSTRACT FROM AUTHOR]

A multiscale modeling approach is developed to compute the phase diagram of the RbF–CsF binary system. The mixing enthalpies of the (Rb,Cs)F solid and liquid solutions are evaluated using density functional theory and classical molecular dynamics calculations, respectively. For the solid solution, 18 different configurations are studied with density functional theory and the surrounded atom model is applied in order to compute the configurational partition function. We also measure the solidus and liquidus equilibria using differential scanning calorimetry. Finally the RbF–CsF phase diagram is constructed using the calculated excess free enthalpies of the solid and liquid solutions and a very good agreement with our experimental data is found. [ABSTRACT FROM AUTHOR]

Room-temperature ionic liquids, RTILs, based on alkylimidazolium cations build up an interesting family of products and are considered as a promising solvent class for chemistry and electrochemistry. Their features are chemical and thermal stability, nonvolatility, high ionic conductivity, and a wide electrochemical potential window and a new RTIL is synthesized, the acidic 1-ethyl-3-methylimidazolium fluoride, which shows unique properties.

We probe the possible inclusion of salt (NaCl) in the ice VII lattice over the pressure range from 2 to 4 gigapascal. We combine data from neutron diffraction experiments under pressure and from computational structure searches based on density functional theory. We observe that the high density amorphous precursor (NaCl·10.2D2O) crystallises during annealing at high pressure in the vicinity of the phase boundary between pure ices VII and VIII. The structure formed is very similar to that of pure ice VII. Our simulations indicate that substituting water molecules in the ice VII lattice with Na+ and Cl− ions would lead to a significant expansion of the lattice parameter. Since this expansion was not observed in our experiments, the ice crystallised is likely to be pure D2O or contains only a small fraction of the ions from the salt solution. [ABSTRACT FROM AUTHOR]

The structure of amorphous NaCl solutions produced by fast quenching is studied as a function of pressure, up to 4 GPa, by combined neutron diffraction experiments and classical molecular dynamics simulations. Similarly to LiCl solutions the system amorphizes at ambient pressure in a dense phase structurally similar to the e-HDA phase in pure water. The measurement of the static structure factor as a function of pressure allowed us to validate a new polarizable force field developed by Tazi et al., 2012, never tested under non-ambient conditions. We infer from simulations that the hydration shells of Na+ cations form well defined octahedra composed of both H2O molecules and Cl− anions at low pressure. These octahedra are gradually broken by the seventh neighbour moving into the shell of first neighbours yielding an irregular geometry. In contrast to LiCl solutions and pure water, the system does not show a polyamorphic transition under pressure. This confirms that the existence of polyamorphism relies on the tetrahedral structure of water molecules, which is broken here. [ABSTRACT FROM AUTHOR]

Supercapacitors are promising devices for energy storage. Being able to measure and predict their performances is a key step in order to optimize them. In the present study, we propose an original methodology to calculate the capacitance of a single nanoporous carbon electrode in contact with an ionic liquid, using molecular dynamics simulations. The results are compared to experimental electrochemical measurements conducted on the same systems at high temperature (close to 100°C). The two approaches are in qualitative agreement and show that, in the case of a butyl-methylimidazolium hexafluorophosphate electrolyte combined with a carbide-derived carbon with an average pore size of 0.9 nm, the positive electrode capacitance is fairly larger than the negative one. [ABSTRACT FROM AUTHOR]

Molten fluorides are known to show favourable thermophysical properties which make them good candidate coolants for nuclear fission reactors. Here we investigate the special case of mixtures of lithium fluoride and thorium fluoride, which act both as coolant and as fuel in the molten salt fast reactor concept. By usingab initioparameterised polarisable force fields, we show that it is possible to calculate the whole set of properties (density, thermal expansion, heat capacity, viscosity and thermal conductivity) which are necessary for assessing the heat transfer performance of the melt over the whole range of compositions and temperatures. We then deduce from our calculations several figures of merit which are important in helping the optimisation of the design of molten salt fast reactors. [ABSTRACT FROM AUTHOR]

The construction of a reliable potential for GeO2 from first principles is described. The obtained potential, which includes dipole polarization effects, is able to reproduce all the studied properties (structural, dynamical and vibrational) to a high degree of precision with a single set of parameters. In particular, the infrared spectrum was obtained using the expression proposed for the dielectric function of polarizable ionic solutions reported by Weis et al. [J. Chem. Phys. 91, 5544 (1989)]. The agreement with the experimental spectrum is very good, with three main bands that are associated with tetrahedral modes of the GeO2 network. Finally, we give a comparison with a simpler pair-additive potential. [ABSTRACT FROM AUTHOR]

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 or 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. [ABSTRACT FROM AUTHOR]

We perform molecular dynamics simulations of a typical nanoporous-carbon based supercapacitor. The organic electrolyte consists in 1-ethyl-3-methylimidazolium 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. [ABSTRACT FROM AUTHOR]

The electrochemical CO 2 reduction reaction (eCO 2 RR) offers a pathway to produce valuable chemical fuels from CO 2 . However, its efficiency in aqueous electrolytes is hindered by the concurrent H 2 evolution reaction (HER), which takes place at similar potentials. While the influence of cations on this process has been extensively studied, the influence of anions remains largely unexplored. In this work, we study how eCO 2 RR selectivity and activity on a gold catalyst are affected by a wide range of inorganic and carboxylate anions. We utilize in situ differential electrochemical mass spectrometry (DEMS) for real-time product monitoring coupled with molecular dynamics (MD) simulations. We show that anions significantly impact eCO 2 RR kinetics and eCO 2 RR selectivity. MD simulations reveal a new descriptor─free energy of anion physisorption─where weakly adsorbing anions enable favorable CO 2 reduction kinetics, despite the negative charge carried by the electrode surface. By leveraging these fundamental insights, we identify propionate as the most promising anion, achieving nearly 100% Faradaic efficiency while showing high CO production rates that are comparable to those in bicarbonate. These insights underscore the vital role of anion selection in achieving a highly efficient eCO 2 RR in aqueous electrolytes.

Water-in-salt (WiS) electrolytes are promising systems for a variety of energy storage devices. Indeed, they represent a great alternative to conventional organic electrolytes thanks to their environmental friendliness, non-flammability, and good electrochemical stability. Understanding the behaviour of such systems and their local organisation is a key direction for their rational design and successful implementation at the industrial scale. In the present paper, we focus our investigation on the 21 m bis(trifluoromethanesulfonyl)imide (LiTFSI) WiS electrolyte, recently reported to have acidic pH values. We explore the speciation of an excess proton in this system and its dependence on the initial local environment using ab initio molecular dynamics simulations. In particular, we observe the formation of HTFSI acid in the WiS system, known to act as a superacid in water. This acid is stabilised in the WiS solution for several picoseconds thanks to the formation of a complex with water molecules and a neighboring TFSI - anion. We further investigate how the excess proton affects the microstructure of WiS, in particular, the recently observed oligomerisation of lithium cations, and we report possible notable perturbations of the lithium nanochain organisation. These two phenomena are particularly important when considering WiS as electrolytes in batteries and supercapacitors, and our results contribute to the comprehension of these systems at the molecular level.

Molecular dynamics (MD) 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, a simple parameterization method that allows to tune the metallicity of the electrode based on a quantum chemistry calculation of the density of states (DOS) is introduced. As a first illustration, the interface between graphite electrodes and two different liquid electrolytes, an aqueous solution of NaCl and a pure ionic liquid, at different applied potentials are studied. It is shown that the simulations reproduce qualitatively the experimentally-measured capacitance; in particular, they yield a minimum of capacitance at the point of zero charge (PZC), which is due to the quantum capacitance (QC) 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., (© 2024 The Author(s). Advanced Materials published by Wiley‐VCH GmbH.)

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis -[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

Capacitive storage devices allow for fast charge and discharge cycles, making them the perfect complements to batteries for high power applications. Many materials display interesting capacitive properties when they are put in contact with ionic solutions despite their very different structures and (surface) reactivity. Among them, nanocarbons are the most important for practical applications, but many nanomaterials have recently emerged, such as conductive metal-organic frameworks, 2D materials, and a wide variety of metal oxides. These heterogeneous and complex electrode materials are difficult to model with conventional approaches. However, the development of computational methods, the incorporation of machine learning techniques, and the increasing power in high performance computing now allow us to tackle these types of systems. In this Review, we summarize the current efforts in this direction. We show that depending on the nature of the materials and of the charging mechanisms, different methods, or combinations of them, can provide desirable atomic-scale insight on the interactions at play. We mainly focus on two important aspects: (i) the study of ion adsorption in complex nanoporous materials, which require the extension of constant potential molecular dynamics to multicomponent systems, and (ii) the characterization of Faradaic processes in pseudocapacitors, that involves the use of electronic structure-based methods. We also discuss how recently developed simulation methods will allow bridges to be made between double-layer capacitors and pseudocapacitors for future high power electricity storage devices.

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/Å, 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 nanosecond timescales in the accurate atomistic description of the electrical double layer., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Optimizing electrosynthetic reactions requires fine tuning of a vast chemical space, including charge transfer at electrocatalyst/electrode surfaces, engineering of mass transport limitations, and complex interactions of reactants and products with their environment. Hybrid electrolytes, in which supporting salt ions and substrates are dissolved in a binary mixture of organic solvent and water, represent a new piece of this complex puzzle as they offer a unique opportunity to harness water as the oxygen or proton source in electrosynthesis. In this work, we demonstrate that modulating water-organic solvent interactions drastically impacts the solvation properties of hybrid electrolytes. Combining various spectroscopies with synchrotron small-angle X-ray scattering (SAXS) and force field-based molecular dynamics (MD) simulations, we show that the size and composition of aqueous domains forming in hybrid electrolytes can be controlled. We demonstrate that water is more reactive for the hydrogen evolution reaction (HER) in aqueous domains than when strongly interacting with solvent molecules, which originates from a change in reaction kinetics rather than a thermodynamic effect. We exemplify novel opportunities arising from this new knowledge for optimizing electrosynthetic reactions in hybrid electrolytes. For reactions proceeding first via the activation of water, fine tuning of aqueous domains impacts the kinetics and potentially the selectivity of the reaction. Instead, for organic substrates reacting prior to water, aqueous domains have no impact on the reaction kinetics, while selectivity may be affected. We believe that such a fine comprehension of solvation properties of hybrid electrolytes can be transposed to numerous electrosynthetic reactions.

Using molecular dynamics simulations and graph-theory-based cluster analysis, we investigate the structure-transport properties of typical water-in-salt electrolytes. We demonstrate that ions exhibit distinct dynamics across different ionic clusters-namely, solvent-separated ion pairs (SSIPs), contact ion pairs (CIPs), and aggregates (AGGs). We assess the average proportions of various ionic species and their lifetimes. Our method reveals a dynamic decoupling of ion kinetics, with each species independently contributing to the overall molecular motion. This is evidenced by the fact that the total velocity autocorrelation function (VACF) and power spectrum can be expressed as a weighted sum of independent functions for each species. The experimental data on the ionic conductivity of the studied LiTFSI electrolytes align well with our theoretical predictions at various concentrations, based on the proportions and diffusion coefficients of free ions derived from our analysis. The insights gained into the solvation structures and dynamics of different ionic species enable us to elucidate the physical mechanisms driving ion transport in such superconcentrated electrolytes, providing a comprehensive framework for the future design and optimization of electrolytes., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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, investigation 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 coarse-grain simulations of all-atom systems and compute effective Lennard-Jones and Coulomb parameters, allowing subsequently 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 are dominant at a long range while van der Waals forces are dominant at a 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 that can be thought of as a data-driven generalization of the path collective variable concept. It consists of 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 reaching an accuracy far from the one achieved 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.

MXenes are 2D transition metal carbides, nitrides, and/or carbonitrides that can be intercalated with cations through chemical or electrochemical pathways. While the insertion of alkali and alkaline earth cations into Ti 3 C 2 T x MXenes is well studied, understanding of the intercalation of redox-active transition metal ions into MXenes and its impact on their electronic and electrochemical properties is lacking. In this work, we investigate the intercalation of Cu ions into Ti 3 C 2 T x MXene and its effect on its electronic and electrochemical properties. Using X-ray absorption spectroscopy (XAS) and ab initio molecular dynamics (AIMD), we observe an unusual phenomenon whereby Cu 2+ ions undergo partial reduction upon intercalation from the solution into the MXene. Furthermore, using in situ XAS, we reveal changes in the oxidation states of intercalated Cu ions and Ti atoms during charging. We show that the pseudocapacitive response of Cu-MXene originates from the redox of both the Cu intercalant and Ti 3 C 2 T x host. Despite highly reducing potentials, Cu ions inside the MXene show an excellent stability against full reduction upon charging. Our findings demonstrate how electronic coupling between Cu ions and Ti 3 C 2 T x modifies electrochemical and electronic properties of the latter, providing the framework for the rational design and utilization of transition metal intercalants for tuning the properties of MXenes for various electrochemical systems.

Water-in-salts (WiSs) have recently emerged as promising electrolytes for energy storage applications ranging from aqueous batteries to supercapacitors. Here, ab initio molecular dynamics is used to study the structure of a 21 m LiTFSI WiS. The simulation reveals a new feature, in which the lithium ions form polymer-like nanochains that involve up to 10 ions. Despite the strong Coulombic interaction between them, the ions in the chains are found at a distance of 2.5 Å. They show a drastically different solvation shell compared to that of the isolated ions, in which they share on average two water molecules. The nanochains have a highly transient character due to the low free energy barrier for forming/breaking them. Providing new insights into the nanostructure of WiS electrolytes, our work calls for reevaluating our current knowledge of highly concentrated electrolytes and the impact of the modification of the solvation of active species on their electrochemical performances.

When a surface is immersed in a solution, it usually acquires a charge, which attracts counterions and repels co-ions to form an electrical double layer. The ions directly adsorbed to the surface are referred to as the Stern layer. The structure of the Stern layer normal to the interface was described decades ago, but the lateral organization within the Stern layer has received scant attention. This is because instrumental limitations have prevented visualization of the ion arrangements except for atypical, model, crystalline surfaces. Here, we use high-resolution amplitude modulated atomic force microscopy (AFM) to visualize in situ the lateral structure of Stern layer ions adsorbed to polycrystalline gold, and amorphous silica and gallium nitride (GaN). For all three substrates, when the density of ions in the layer exceeds a system-dependent threshold, correlation effects induce the formation of close packed structures akin to Wigner crystals. Depending on the surface and the ions, the Wigner crystal-like structure can be hexagonally close packed, cubic, or worm-like. The influence of the electrolyte concentration, species, and valence, as well as the surface type and charge, on the Stern layer structures is described. When the system parameters are changed to reduce the Stern layer ion surface excess below the threshold value, Wigner crystal-like structures do not form and the Stern layer is unstructured. For gold surfaces, molecular dynamics (MD) simulations reveal that when sufficient potential is applied to the surface, ion clusters form with dimensions similar to the Wigner crystal-like structures in the AFM images. The lateral Stern layer structures presented, and in particular the Wigner crystal-like structures, will influence diverse applications in chemistry, energy storage, environmental science, nanotechnology, biology, and medicine.

Among a variety of promising cathode materials for Na-ion batteries, polyanionic Na-insertion compounds are among the preferred choices due to known fast sodium transfer through the ion channels along their framework structures. The most interesting representatives are Na 3 V 2 (PO 4 ) 3 (NVP) and Na 3 V 2 (PO 4 ) 2 F 3 (NVPF), which display large Na + diffusion coefficients (up to 10 -9 m 2 s -1 in NVP) and high voltage plateaux (up to 4.2 V for NVPF). While the diffusion in the solid material is well-known to be the rate-limiting step during charging, already being thoroughly discussed in the literature, interfacial transport of sodium ions from the liquid electrolyte toward the electrode was recently shown to be important due to complex ion desolvation effects at the surface. In order to fill the blanks in the description of the electrode/electrolyte interface in Na-ion batteries, we performed a molecular dynamics study of the local nanostructure of a series of carbonate-based sodium electrolytes at the NVP and the NVPF interfaces along with careful examination of the desolvation phenomenon. We show that the tightness of solvent packing at the electrode surface is a major factor determining the height of the free energy barrier associated with desolvation, which explains the differences between the NVP and the NVPF structures. To rationalize and emphasize the remarkable properties of this family of cathode materials, a complementary comparative analysis of the same electrolyte system at the carbon electrode interface was also performed.

Photoacids are organic molecules that release protons under illumination, providing spatiotemporal control of the pH. Such light-driven pH switches offer the ability to cyclically alter the pH of the medium and are highly attractive for a wide variety of applications, including CO 2 capture. Although photoacids such as protonated merocyanine can enable fully reversible pH cycling in water, they have a limited chemical stability against hydrolysis (<24 h). Moreover, these photoacids have low solubility, which limits the pH-switching ability in a buffered solution such as dissolved CO 2 . In this work, we introduce a simple pathway to dramatically increase stability and solubility of photoacids by tuning their solvation environment in binary solvent mixtures. We show that a preferential solvation of merocyanine by aprotic solvent molecules results in a 60% increase in pH modulation magnitude when compared to the behavior in pure water and can withstand stable cycling for >350 h. Our results suggest that a very high stability of merocyanine photoacids can be achieved in the right solvent mixtures, offering a way to bypass complex structural modifications of photoacid molecules and serving as the key milestone toward their application in a photodriven CO 2 capture process., Competing Interests: The authors declare no competing financial interest., (© 2023 American Chemical Society.)

BaSnF4 is a prospective solid state electrolyte for fluoride ion batteries. However, the diffusion mechanism of the fluoride ions remains difficult to study, both in experiments and in simulations. In principle, ab initio molecular dynamics could allow to fill this gap, but this method remains very costly from the computational point of view. Using machine learning potentials is a promising method that can potentially address the accuracy issues of classical empirical potentials while maintaining high efficiency. In this work, we fitted a dipole polarizable ion model and trained machine learning potential for BaSnF4 and made comprehensive comparisons on the ease of training, accuracy and efficiency. We also compared the results with the case of a simpler ionic system (NaF). We show that contrarily to the latter, for BaSnF4 the machine learning potential offers much higher versatility. The current work lays foundations for the investigation of fluoride ion mobility in BaSnF4 and provides insight on the choice of methods for atomistic simulations., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Redox-active organic species play an important role in catalysis, energy storage, and biotechnology. One of the representatives is the 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) radical, used as a mediator in organic synthesis and considered a safe alternative to heavy metals. In order to develop a TEMPO-based system with well-controlled electrochemical and catalytic properties, a reaction medium should be carefully chosen. Being highly conductive, stable, and low flammability fluids, ionic liquids (ILs) seem to be promising solvents with easily adjustable physical and solvation properties. In this work, we give an insight into the local structure of ILs around TEMPO and its oxidized form, TEMPO + , underlining striking differences in the solvation of these two species. The analysis is coupled with a study of thermodynamics and kinetics of oxidation in the frame of Marcus theory. Our systematic investigation includes imidazolium, pyrrolydinium, and phosphonium families combined with anions of different size, polarity, and flexibility, opting to provide a clear and comprehensive picture of the impact of the nature of IL ions on the behavior of radical/cation redox pairs. The obtained results will help to explain experimentally observed effects and to rationalize the design of TEMPO/IL systems.

Ionic liquids generally display peculiar structural features that impact their physical properties, such as the formation of polar and apolar domains. Recently, ionic liquids functionalized with anthraquinone and TEMPO redox groups were shown to increase the energy storage performance of supercapacitors, but their structure has not yet been characterized. In this work, we use polarizable molecular dynamics to study the nanostructuration of such biredox ionic liquids. We show that TEMPO nitroxyl functions tend to aggregate, while the anthraquinone groups favor stacked arrangements. The latter eventually percolate through the whole liquid, which sheds some light on the mechanisms at play within biredox ionic liquid-based supercapacitors.

The electrocatalytic epoxidation of alkenes at heterogeneous catalysts using water as the sole oxygen source is a promising safe route toward the sustainable synthesis of epoxides, which are essential building blocks in organic chemistry. However, the physicochemical parameters governing the oxygen-atom transfer to the alkene and the impact of the electrolyte structure on the epoxidation reaction are yet to be understood. Here, we study the electrocatalytic epoxidation of cyclooctene at the surface of gold in hybrid organic/aqueous mixtures using acetonitrile (ACN) solvent. Gold was selected, as in ACN/water electrolytes gold oxide is formed by reactivity with water at potentials less anodic than the oxygen evolution reaction (OER). This unique property allows us to demonstrate that a sacrificial mechanism is responsible for cyclooctene epoxidation at metallic gold surfaces, proceeding through cyclooctene activation, while epoxidation at gold oxide shares similar reaction intermediates with the OER and proceeds via the activation of water. More importantly, we show that the hydrophilicity of the electrode/electrolyte interface can be tuned by changing the nature of the supporting salt cation, hence affecting the reaction selectivity. At low overpotential, hydrophilic interfaces formed using strong Lewis acid cations are found to favor gold passivation. Instead, hydrophobic interfaces created by the use of large organic cations favor the oxidation of cyclooctene and the formation of epoxide. Our study directly demonstrates how tuning the hydrophilicity of electrochemical interfaces can improve both the yield and selectivity of anodic reactions at the surface of heterogeneous catalysts.

Metallic 1T-MoS 2 is a promising electrode material for supercapacitor applications. Its layered structure allows the efficient intercalation of ions, leading to experimental volumetric capacitance as high as 140 F/cm 3 . Molecular dynamics could in principle be used to characterize its charging mechanism; however, unlike conventional nanoporous carbon, 1T-MoS 2 is a multicomponent electrode. The Mo and S atoms have very different electronegativities so that 1T-MoS 2 cannot be simulated accurately using the conventional constant potential method. In this work, we show that controlling the electrochemical potential of the atoms allows one to recover average partial charges for the elements in agreement with electronic structure calculations for the material at rest, without compromising the ability to simulate systems under an applied voltage. The simulations yield volumetric capacitances in agreement with experiments. We show that due to the large electronegativity of S, the co-ion desorption is the main charging mechanism at play during the charging process. This contrasts drastically with carbon materials for which ion exchange and counterion adsorption usually dominate. In the future, our method can be extended to the study of a wide range of families of 2D layered materials such as MXenes.

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.

Quinones are redox active organic molecules that have been proposed as an alternative choice to metal-based materials in electrochemical energy storage devices. Functionalization allows one to fine tune not only their chemical stability but also the redox potential and kinetics of the electron transfer reaction. However, the reaction rate constant is not only determined by the redox species but also impacted by solvent effects. In this work, we show how the functionalization of benzoquinone with different functional groups impacts the solvent reorganization free energies of electron transfer half-reactions in acetonitrile. The use of molecular density functional theory, whose computational cost for studying the electron transfer reaction is considerably reduced compared to the state-of-the-art molecular dynamics simulations, enables us to perform a systematic study. We validate the method by comparing the predictions of the solvation shell structure and the free energy profiles for electron transfer reaction to the reference classical molecular dynamics simulations in the case of anthraquinone solvated in acetonitrile. We show that all the studied electron transfer half-reactions follow the Marcus theory, regardless of functional groups. Consequently, the solvent reorganization free energy decreases as the molecular size increases.

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

Li-ion batteries are the electrochemical energy storage technology of choice of today's electrical vehicles and grid applications with a growing interest for Na-ion and K-ion systems based on either aqueous or non-aqueous electrolyte for power, cost, and sustainable reasons. The rate capability of alkali-metal-ion batteries is influenced by ion transport properties in the bulk of the electrolyte, as well as by diverse effects occurring at the vicinity of the electrode and electrolyte interface. Therefore, identification of the predominant factor affecting the rate capability of electrodes still remains a challenge and requires suitable experimental and computational methods. Herein, we investigate the mechanistic of the K + insertion process in the Prussian blue phase, F e 4 I I I [ F e I I ( C N ) 6 ] 3 in both aqueous and non-aqueous electrolytes, which reveals drastic differences. Through combined electrochemical characterizations, electrochemical-quartz-crystal-microbalance and ac -electrogravimetric analyses, we provide evidences that what matters the most for fast ion transport is the positioning of the partially solvated cations adsorbed at the material surface in aqueous as opposed to non-aqueous electrolytes. We rationalized such findings by molecular dynamics simulations that establish the K + repartition profile within the electrochemical double layer. A similar trend was earlier reported by our group for the aqueous versus non-aqueous insertion of Li + into LiFePO 4 . Such a study unveils the critical but overlooked role of the electrode-electrolyte interface in ruling ion transport and insertion processes. Tailoring this interface structuring via the proper salt-solvent interaction is the key to enabling the best power performances in alkali-metal-ion batteries.

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 aqueous solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer.