Search

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

We revisit the statistical mechanics of charge fluctuations in capacitors. In constant-potential classical molecular simulations, the atomic charges 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 continuum 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. [ABSTRACT FROM AUTHOR]

Machine learning (ML) approaches are investigated for the prediction of nuclear magnetic resonance (NMR) parameters in aluminosilicate glasses, for which NMR has proven to be a cutting-edge method over the last decade. DFT computations have emerged as a new dimension for complementing these NMR methods although suffering from severe limitations in terms of size, time and computational resources consumption. While previous approaches tend to use DFT-GIPAW calculations for the prediction of NMR parameters in glassy systems, we propose to employ ML methods, characterized by a speed similar to that of classical molecular dynamics while the accuracy of ab initio methods can be reached. We design ML procedures to predict the isotropic magnetic shielding (σiso) for different multicomponent relevant glass compositions. The ML predictions of σiso deviate from DFT-GIPAW calculations, when including relaxed and room-temperature structures, by 0.7 ppm for 29Si (1.0% of the total span of the calculated) and 1.5 ppm for 17O (1.9%) in SiO2 glasses, 1.4 ppm for 23Na (1.5%) in Na2O–SiO2 and 1.5 ppm for 27Al (2.1%) in Al2O3–Na2O–SiO2 systems. We compare the performances obtained for a set of three descriptors suitable for encoding atomic local environments information (atom-centered representations) together with seven popular ML algorithms with a focus on the simple (but robust) linear ridge regression (LRR) and the popular smooth overlap of atomic positions (SOAP) descriptor. [ABSTRACT FROM AUTHOR]

Mixtures of trivalent metal halides with alkali halides are involved in many technologies but, from a more fundamental and general perspective, are worthy of study as interesting systems in which to examine the relationship between atomic-scale structure and physical properties. Here we examine the relationship between the viscosity and local and longer range structural measures in such mixtures where the trivalent metal cations span a significant size range and exhibit different behaviours in the dependence of their viscosity on the mixture composition. We characterise the structure and dynamics of the first coordination shell and the relationship between its structural relaxation time and the shear relaxation time of the mixture (the Maxwell relaxation time). We are then led to an examination of the structure of the networks which progressively form between the trivalent metal cations as their concentration increases in the mixtures. Here we find significant differences between small and larger cations, sufficient to explain the different behaviour of their viscosities. We draw attention to the similarities and differences of these networks with those which form in highly viscous, glass-forming materials like BeF2:LiF. [ABSTRACT FROM AUTHOR]

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

Ceria (CeO2) 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. [ABSTRACT FROM AUTHOR]

Ionic liquids are widely used as electrolytes in electrochemical devices. In this context, many experimental and theoretical approaches have been recently developed for characterizing their interface with electrodes. In this perspective article, we review the most recent advances in the field of computer simulations (mainly molecular dynamics). A methodology for simulating electrodes at constant electrical potential is presented. Several types of electrode geometries have been investigated by many groups in order to model planar, corrugated and porous materials and we summarize the results obtained in terms of the structure of the liquids. This structure governs the quantity of charge which can be stored at the surface of the electrode for a given applied potential, which is the relevant quantity for the highly topical use of ionic liquids in supercapacitors (also known as electrochemical double-layer capacitors). A key feature, which was also shown by atomic force microscopy and surface force apparatus experiments, is the formation of a layered structure for all ionic liquids at the surface of planar electrodes. This organization cannot take place inside nanoporous electrodes, which results in a much better performance for the latter in supercapacitors. The agreement between simulations and electrochemical experiments remains qualitative only though, and we outline future directions which should enhance the predictive power of computer simulations. In the longer term, atomistic simulations will also be applied to the case of electron transfer reactions at the interface, enabling the application to a broader area of problems in electrochemistry, and the few recent works in this field are also commented upon. [ABSTRACT FROM AUTHOR]

In an inorganic material, a measure of the Lewis basicity of a solvent is commonly provided by the optical basicity, i.e.the ability of the solvent molecules to donate their electrons to an acidic species. This quantity is known to vary with the polarizability of the Lewis base. In protic solvents, the Lewis definition of basicity is barely used; it is replaced by the more purpose-built Brønsted–Lowry scale, and its generalized variant proposed by Hammett. In this study, individual molecular polarizabilities were computed from first-principles for a series of protic solvents: pure water, hydrogen fluoride, ammonia and their mixtures. From these calculations optical basicity scales were set up for each Lewis base. It was shown that these scales correlate with the Hammett acidity. It is therefore possible to build a common optical basicity scale, in which any material (protic solvents, inorganic materials) can unambiguously be classified. [ABSTRACT FROM AUTHOR]

The internal mobility gives the rate at which one ionic species moves relative to the other species present in an ionic mixture, it mirrors the differential strength of the interactions between different ionic species. In this work we examine the dependence of the internal mobilities of the Li+and K+ions on the composition in molten mixtures of LiF and KF. We compare them to the behaviour of the individual diffusion coefficients and the self-exchange velocities, which measure the rate at which an ion separates from its nearest-neighbour coordination shell. The examination is made using molecular dynamics simulations with polarizable, first-principles parameterised interaction potentials which are shown to reproduce the limited available experimental data on the transport properties of these mixtures extremely well. The results confirm that the composition-dependence of the internal mobilities in LiF/KF follows the unusual type-II behaviour, which is not reflected in that of the diffusion coefficients or the self-exchange velocities. [ABSTRACT FROM AUTHOR]

The local structure and the dynamics of molten LiF–KF mixtures have been studied by nuclear magnetic resonance (NMR) and molecular dynamics simulations. We have measured and calculated the self-diffusion coefficients of fluorine, lithium and potassium across the full composition range around the liquidus temperature and at 1123 K. Close to the liquidus temperature, DF, DLiand DKchange with composition in a way that mimics the phase diagram shape. At 1123 K DF, DLiand DKdepend linearly on the LiF molar fraction. These results show that the composition affects the self-diffusion of anions and cations more weakly than the temperature. The activation energy for diffusion was also determined and its value can be correlated with the strength of the anion–cation interaction in molten fluoride salts. [ABSTRACT FROM AUTHOR]

Chemical doping and other surface modifications have been used to engineer the bulk properties of materials, but their influence on the surface structure and consequently the surface chemistry are often unknown. Previous work has been successful in fluorinating anatase TiO 2 with charge balance achieved via the introduction of Ti vacancies rather than the reduction of Ti. Our work here investigates the interface between this fluorinated titanate with cationic vacancies and a monolayer of water via density functional theory based molecular dynamics. We compute the projected density of states for only those atoms at the interface and for those states that fall within 1 eV of the Fermi level for various steps throughout the simulation, and we determine that the variation in this visualization of the density of states serves as a reasonable tool to anticipate where surfaces are most likely to be reactive. In particular, we conclude that water dissociation at the surface is the main mechanism that influences the anatase (001) surface whereas the change in the density of states at the surface of the fluorinated structure is influenced primarily through the adsorption of water molecules., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Beyond the dielectric continuum description initiated by Marcus theory, the 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 present 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 within 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 → Cl + and Cl 0 → 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 a solid-solvent interface on the ET, by investigating the evolution of the reorganization free energy of the Cl 0 → Cl + reaction when the atom approaches an atomistically resolved wall.

The local structure and the dynamics of molten LiF-KF mixtures have been studied by nuclear magnetic resonance (NMR) and molecular dynamics simulations. We have measured and calculated the self-diffusion coefficients of fluorine, lithium and potassium across the full composition range around the liquidus temperature and at 1123 K. Close to the liquidus temperature, D(F), D(Li) and D(K) change with composition in a way that mimics the phase diagram shape. At 1123 K D(F), D(Li) and D(K) depend linearly on the LiF molar fraction. These results show that the composition affects the self-diffusion of anions and cations more weakly than the temperature. The activation energy for diffusion was also determined and its value can be correlated with the strength of the anion-cation interaction in molten fluoride salts.