Polarizable Force Field Development and Molecular Dynamics Simulations of Ionic Liquids

Modeling of Nanomaterials for Supercapacitors: Beyond Carbon Electrodes

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.

Inhibition of silica scaling with functional polymers: Role of ionic strength, divalent ions, and temperature

High concentrations of dissolved silica in saline industrial wastewaters and brines cause silica scale formation, significantly hampering the efficacy of diverse engineered systems. Applying functional polymers as scale inhibitors in process feedwater is a common strategy to mitigate silica scaling. However, feedwater characteristics often vary widely, depending on the specific processes, making the inhibition of silica scaling challenging and complex. In this study, we systematically investigate the role of ionic composition, specifically ionic strength and divalent ions, and solution temperature, in inhibiting silica scaling using molecularly designed amine/amide polymers. The inhibitor demonstrates effective stabilization of silicic acid, with inhibition efficiency of 74 and 55 % in the absence and presence of 20,000 ppm NaCl, respectively. However, further increasing the ionic strength of oversaturated silicic acid solutions significantly diminishes inhibition performance, rendering it ineffective at 180,000 ppm NaCl. Divalent inorganic cations exhibit a stronger impact on reducing inhibition efficiency compared to sodium ions. Molecular dynamics simulations reveal a competition mechanism between anionic silicic acid reactants (i.e., H3SiO4-) and chlorides for binding to ammonium groups within the polymeric inhibitor. Additionally, cations form clusters with H3SiO4- ions, hindering their stabilization with polymeric inhibitor. Notably, at elevated temperatures, the inhibitor achieves near-perfect inhibition for 500 ppm silicic acid solutions. This comprehensive assessment provides important insights into the effectiveness of silica scaling inhibitors under solution conditions relevant to real-world applications, addressing the challenges posed by varying solution parameters in diverse industrial processes.

Sodium Triflate Water-in-Salt Electrolytes in Advanced Battery Applications: A First-Principles-Based Molecular Dynamics Study

Offering a compelling combination of safety and cost-effectiveness, water-in-salt (WiS) electrolytes have emerged as promising frontiers in energy storage technology. Still, there is a strong demand for research and development efforts to make these electrolytes ripe for commercialization. Here, we present a first-principles-based molecular dynamics (MD) study addressing in detail the properties of a sodium triflate WiS electrolyte for Na-ion batteries. We have developed a workflow based on a machine learning (ML) potential derived from ab initio MD simulations. As ML potentials are typically restricted to the interpolation of the data points of the training set and have hardly any predictive properties, we subsequently optimize a classical force field based on physics principles to ensure broad applicability and high performance. Performing and analyzing detailed MD simulations, we identify several very promising properties of the sodium triflate as a WiS electrolyte but also indicate some potential stability challenges associated with its use as a battery electrolyte.

Vehicular Motions Dominate the Ion Transport in Concentrated LiTFSI Aqueous Solutions?

Water-in-salt electrolytes (WiSEs) show great promise for applications in grid-scale energy storage. The design of high-performance WiSEs requires a comprehensive understanding of their microstructures and ion transport properties. In the present work, based on the CL&Pol force field, we have developed a polarizable force field (PFF) tailored for high-concentration LiTFSI aqueous solutions, which accurately reproduces the structural and dynamical properties. Unlike the literature, we do not observe the presence of bulk-like water in LiTFSI solutions exceeding 19 mol/kg. Furthermore, we find that the vast majority of Li(H20)n+ are short-lived, and thus, the structural motion rather than the vehicular motion is the main mode of ion transport. Our results have significant implications for understanding the ion dynamics in WiSEs. Additionally, further in-depth experimental analyses are imperative.

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB)

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. We present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.

System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields

The accuracy of classical force fields (FFs) has been shown to be limited for the simulation of cation-protein systems despite their importance in understanding the processes of life. Improvements can result from optimizing the parameters of classical FFs or by extending the FF formulation by terms describing charge transfer (CT) and polarization (POL) effects. In this work, we introduce our implementation of the CTPOL model in OpenMM, which extends the classical additive FF formula by adding CT and POL. Furthermore, we present an open-source parametrization tool, called FFAFFURR, that enables the (system-specific) parametrization of OPLS-AA and CTPOL models. The performance of our workflow was evaluated by its ability to reproduce quantum chemistry energies and by molecular dynamics simulations of a zinc-finger protein.

High voltage electrolytes for lithium-ion batteries with micro-sized silicon anodes

Micro-sized silicon anodes can significantly increase the energy density of lithium-ion batteries with low cost. However, the large silicon volume changes during cycling cause cracks for both organic-inorganic interphases and silicon particles. The liquid electrolytes further penetrate the cracked silicon particles and reform the interphases, resulting in huge electrode swelling and quick capacity decay. Here we resolve these challenges by designing a high-voltage electrolyte that forms silicon-phobic interphases with weak bonding to lithium-silicon alloys. The designed electrolyte enables micro-sized silicon anodes (5 µm, 4.1 mAh cm-2) to achieve a Coulombic efficiency of 99.8% and capacity of 2175 mAh g-1 for >250 cycles and enable 100 mAh LiNi0.8Co0.15Al0.05O2 pouch full cells to deliver a high capacity of 172 mAh g-1 for 120 cycles with Coulombic efficiency of >99.9%. The high-voltage electrolytes that are capable of forming silicon-phobic interphases pave new ways for the commercialization of lithium-ion batteries using micro-sized silicon anodes.

Liquid Structure and Hydrogen Bonding in Aqueous Hydroxylammonium Nitrate

Hydroxylammonium nitrate (HAN) has emerged as a promising component in ionic liquid-based spacecraft propellants. However, the physicochemical and structural properties of aqueous HAN have been largely overlooked. The purpose of this study is to investigate the hydrogen bonding in aqueous HAN and understand its implications on these properties and the proton transfer mechanism as a function of concentration. Classical polarizable molecular dynamics simulations have been employed with the APPLE&P force field to analyze the geometry of individual hydrogen bonds and the overall hydrogen-bonding network in various concentrations of aqueous HAN. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) indicate the structural arrangement of the species and their hydrogen bonds. Projections of water density and the orientation of its electric dipole moment near the ions provide insight into the hydrogen-bonding network. The incorporation of water into the hydrogen-bonding network at high ion concentrations occurs via interstitial accommodation around the ions immediately outside the first solvation shell. While ion pairs are observed at all concentrations considered, the frequency of Ha···On hydrogen bonds increases substantially with the ion concentration. The findings contribute to a better fundamental understanding of HAN and the precursors of reactivity, crucial to the development of "green" spacecraft propellants.

GAFF-Based Polarizable Force Field Development and Validation for Ionic Liquids

Ionic liquids (ILs) have been used in many applications, including gas separations, electrochemistry, lubrication, and catalysis. Understanding how the different properties of ILs are related to their chemical structure and composition is crucial for these applications. Experimental investigations often provide limited insights and can be tedious in exploring a range of state points. Therefore, molecular simulations have emerged as a powerful tool that not only offers a microscopic perspective but also enables rapid screening and prediction of physical properties. The accuracy of these predictions, however, depends on the quality of the intermolecular potentials (force fields) used. The widely used classical fixed charge models, such as GAFF, OPLS, and CL&P, are popular due to their simplicity and computational efficiency. However, it has been shown that the use of integer charges with these classical models leads to sluggish dynamics. The use of scaled charge models can improve the dynamics, but these mean-field approaches are unable to account for polarization effects explicitly. Several different approaches have been proposed to include polarizability in IL force fields. In this work, we follow the protocol of the CL&Pol model to develop a Drude oscillator model based on the GAFF force field (Goloviznina, K., et al. J. Chem. Theory Comput. 2019, 15, 5858). We compare the performance of the model for eight imidazolium- and pyrrolidinium-based ILs against that of other models. We find that the new model provides reasonable estimations of density, self-diffusivity, and structural properties for these ILs and suggests a relatively simple way of extending the general GAFF model to more ILs.

Implementation and Validation of an OpenMM Plugin for the Deep Potential Representation of Potential Energy

Machine learning potentials, particularly the deep potential (DP) model, have revolutionized molecular dynamics (MD) simulations, striking a balance between accuracy and computational efficiency. To facilitate the DP model's integration with the popular MD engine OpenMM, we have developed a versatile OpenMM plugin. This plugin supports a range of applications, from conventional MD simulations to alchemical free energy calculations and hybrid DP/MM simulations. Our extensive validation tests encompassed energy conservation in microcanonical ensemble simulations, fidelity in canonical ensemble generation, and the evaluation of the structural, transport, and thermodynamic properties of bulk water. The introduction of this plugin is expected to significantly expand the application scope of DP models within the MD simulation community, representing a major advancement in the field.

Anomalous diffusion of lithium-anion clusters in ionic liquids

Lithium-ion transport is significantly retarded in ionic liquids (ILs). In this work, we performed extensive molecular dynamics simulations to mimic the kinetics of lithium ions in ILs using [N-methyl-N-propylpyrrolidium (pyr[Formula: see text])][bis(trifluoromethanesulfonyl)imide (Ntf[Formula: see text])] with added LiNtf[Formula: see text] salt. And we analyzed their transport, developing a two-state model and comparing it to the machine learning-identified states. The transport of lithium ions involves local shell exchanges of the Ntf[Formula: see text] in the medium. We calculated train size distributions over various time scales. The train size distribution decays as a power law, representing non-Poissonian bursty shell exchanges. We analyzed the non-Poissonian processes of lithium ions transport as a two-state (soft and hard) model. We analytically calculated the transition probability of the two-state model, which fits well to the lifetime autocorrelation functions of LiNtf[Formula: see text] shells. To identify two states, we introduced the graph neutral network incorporating local molecular structure. The results reveal that the shell-soft state mainly contributes to the transport of the lithium ions, and their contribution is more important in low temperatures. Hence, it is the key for enhanced lithium ion transport to increase the fraction of the shell-soft state.

Structural Properties of HEHN- and HAN-Based Ionic Liquid Mixtures: A Polarizable Molecular Dynamics Study

Molecular dynamics simulations of binary mixtures comprising 2-hydroxyethylhydrazinium nitrate (HEHN) and hydroxylammonium nitrate (HAN) were conducted using the polarizable APPLE&P force field to investigate fundamental properties of multimode propulsion (MMP) propellants. Calculated densities as a function of temperature were in good agreement with experiments and similar simulations. The structural properties of neat HEHN and HAN-HEHN provided insights into their inherent, protic nature. Radial distribution functions (RDFs) identified key hydrogen bonding sites located at N-H···O and O-H···O within a first solvation shell of approximately 2 Å. Angular distribution functions further affirmed the relatively strong nature of the hydrogen bonds with nearly linear directionality. The increased hydroxylammonium cation (HA+) mole fraction shows the influence of competitively strong hydrogen bonds on the overall hydrogen bond network. Dominant spatial motifs via three-dimensional distribution functions along with nearly nanosecond-long hydrogen bond lifetimes highlight the local bonding environment that may precede proton transfer reactions.

Molecular Modeling of Water-in-Salt Electrolytes: A Comprehensive Analysis of Polarization Effects and Force Field Parameters in Molecular Dynamics Simulations

Accurate modeling of highly concentrated aqueous solutions, such as water-in-salt (WiS) electrolytes in battery applications, requires proper consideration of polarization contributions to atomic interactions. Within the force field molecular dynamics (MD) simulations, the atomic polarization can be accounted for at various levels. Nonpolarizable force fields implicitly account for polarization effects by incorporating them into their van der Waals interaction parameters. They can additionally mimic electron polarization within a mean-field approximation through ionic charge scaling. Alternatively, explicit polarization description methods, such as the Drude oscillator model, can be selectively applied to either a subset of polarizable atoms or all polarizable atoms to enhance simulation accuracy. The trade-off between simulation accuracy and computational efficiency highlights the importance of determining an optimal level of accounting for atomic polarization. In this study, we analyze different approaches to include polarization effects in MD simulations of WiS electrolytes, with an example of a Na-OTF solution. These approaches range from a nonpolarizable to a fully polarizable force field. After careful examination of computational costs, simulation stability, and feasibility of controlling the electrolyte properties, we identify an efficient combination of force fields: the Drude polarizable force field for salt ions and non-polarizable models for water. This cost-effective combination is sufficiently flexible to reproduce a broad range of electrolyte properties, while ensuring simulation stability over a relatively wide range of force field parameters. Furthermore, we conduct a thorough evaluation of the influence of various force field parameters on both the simulation results and technical requirements, with the aim of establishing a general framework for force field optimization and facilitating parametrization of similar systems.

Mechanistic understanding of the correlation between structure and dynamics of liquid carbonate electrolytes: impact of polarization

Liquid electrolyte design and modelling is an essential part of the development of improved lithium ion batteries. For mixed organic carbonates (ethylene carbonate (EC) and ethyl-methyl carbonate (EMC) mixtures)-based electrolytes with LiPF₆ as salt, we have compared a polarizable force field with the standard non-polarizable force field with and without charge rescaling to model the structural and dynamic properties. The result of our molecular dynamics simulations shows that both polarizable and non-polarizable force fields have similar structural factors, which are also in agreement with X-ray diffraction experimental results. In contrast, structural differences are observed for the lithium neighborhood, while the lithium-anion neighbourhood is much more pronounced for the polarizable force field. Comparison of EC/EMC coordination statistics with Fourier transformed infrared spectroscopy (FTIR) shows the best agreement for the polarizable force field. Also for transport quantities such as ionic conductivities, transference numbers, and viscosities, the agreement with the polarizable force field is by far better for a large range of salt concentrations and EC : EMC ratios. In contrast, for the non-polarizable variants, the dynamics are largely underestimated. The excellent performance of the polarizable force field is explored in different ways to pave the way to a realistic description of the structure-dynamics relationships for a wide range of salt and solvent compositions for this standard electrolyte. In particular, we can characterize the distinct correlation terms between like and unlike ions, relate them to structural properties, and explore to which degree the transport in this electrolyte is mass or charge limited.

Molecular Dynamics Simulations of High-Concentration Li[TFSA] Sulfone Solution: Effect of Easy Conformation Change of Sulfolane on Fast Diffusion of Li Ion

The parameters of the polarizable force field used for molecular dynamics simulations of Li diffusion in high-concentration lithium bis(trifluoromethanesulfonyl)amide (Li[TFSA]) sulfone (sulfolane, dimethylsulfone, ethylmethylsulfone, and ethyl-i-propylsulfone) solutions were refined. The densities of the solutions obtained by molecular dynamics simulations reproduced well the experimental values. The calculated concentration, temperature, and solvent dependencies of self-diffusion coefficients of ions and solvents in the mixtures well reproduce the experimentally observed dependencies. Ab initio calculations show that the intermolecular interactions between Li ions and four sulfones are not largely different. Conformational analyses show that sulfolane can change the conformation more easily owing to lower barrier height for pseudorotation compared to the rotational barrier heights of diethylsulfone and ethylmethylsulfone. Molecular dynamics simulations indicate that the easy conformation change of solvent affects the rotational relaxation of the solvent and the diffusion of Li ion in the mixture. The easy conformation change of sulfolane is one of the causes of faster diffusion of Li ion in the mixture of Li[TFSA] and sulfolane compared to the mixtures of smaller dimethylsulfone and ethylmethylsulfone.

Simulation of polymeric mixed ionic and electronic conductors with a combined classical and quantum mechanical model

In organic polymeric materials with mixed ionic and electronic conduction (OMIEC), the excess charge in doped polymers is very mobile and the dynamics of the polymer chain cannot be accurately described with a model including only fixed point charges. Ions and polymer are comparatively slower and a methodology to capture the correlated motions of excess charge and ions is currently unavailable. Considering a prototypical interface encountered in this type of materials, we constructed a scheme based on the combination of MD and QM/MM to evaluate the classical dynamics of polymer, water and ions, while allowing the excess charge of the polymer chains to rearrange following the external electrostatic potential. We find that the location of the excess charge varies substantially between chains. The excess charge changes across multiple timescales as a result of fast structural fluctuations and slow rearrangement of the polymeric chains. Our results indicate that such effects are likely important to describe the phenomenology of OMIEC, but additional features should be added to the model to enable the study of processes such as electrochemical doping.

Theory and Simulations of Ionic Liquids in Nanoconfinement

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL-water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]⁺ cation and nitrate [NO₃]^- anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV^E, and liquid structures gr. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid-water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid-water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ∼0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. The ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

AMOEBA Polarizable Force Field for Molecular Dynamics Simulations of Glyme Solvents

Classical molecular dynamics (MD) simulations of electrolyte systems are important to gain insight into the atom-scale properties that determine the battery-relevant performance. The recent Tinker-HP software release enables efficient and accurate MD simulations with the AMOEBA polarizable force field. In this work, we developed a procedure to construct a universal AMOEBA model for the solvent family of glymes (glycol methyl ethers), which involves a refinement scheme for valence parameters by fitting the AMOEBA-derived atomic forces to those computed at the DFT level. The refined AMOEBA model provides a good description of both local and nonlocal properties in terms of the spectroscopic response of glyme molecules, as well as the liquid glyme density and dielectric constant. In addition, the complexation energies of alkali and alkaline-earth metal cations with tetraglyme molecules obtained from AMOEBA calculations are in good agreement with DFT results, demonstrating the suitability of the developed AMOEBA model for an accurate simulation of glyme-based battery electrolytes. We also expect the procedure to be transferable to the development of AMOEBA models for other battery electrolyte systems.

Effect of the Electric Double Layer (EDL) in Multicomponent Electrolyte Reduction and Solid Electrolyte Interphase (SEI) Formation in Lithium Batteries

Electrolytes, consisting of salts, solvents, and additives, must form a stable solid electrolyte interphase (SEI) to ensure the performance and durability of lithium(Li)-ion batteries. However, the electric double layer (EDL) structure near charged surfaces is still unsolved, despite its importance in dictating the species being reduced for SEI formation near a negative electrode. In this work, a newly developed model was used to illustrate the effect of EDL on SEI formation in two essential electrolytes, the carbonate-based electrolyte for Li-ion batteries and the ether-based electrolyte for batteries with Li-metal anodes. Both electrolytes have fluoroethylene carbonate (FEC) as a common additive to form the beneficial F-containing SEI component (e.g., LiF). However, the role of FEC drastically differs in these electrolytes. FEC is an effective SEI modifier for the carbonate-based electrolyte by being the only F-containing species entering the EDL and being reduced, as the anion (PF₆^-) will not enter the EDL. For the ether-based electrolyte, both the anion (TFSI^-) and FEC can enter the EDL and be reduced. The competition of the two species within the EDL due to the surface charge and temperature leads to a unique temperature effect observed in prior experiments: the FEC additive is more effective in modulating SEI components at a low temperature (-40 °C) than at room temperature (20 °C) in the ether-based electrolyte. These collective quantitative agreements with experiments emphasize the importance of incorporating the effect of the EDL in multicomponent electrolyte reduction reactions in simulations/experiments to predict/control the formation of the SEI layer.

Ion-Solvent Interplay in Concentrated Electrolytes Enables Subzero Temperature Li-Ion Battery Operations

Despite the essential role of ethylene carbonate (EC) in solid electrolyte interphase (SEI) formation, the high Li⁺ desolvation barrier and melting point (36 °C) of EC impede lithium-ion battery operation at low temperatures and induce sluggish Li⁺ reaction kinetics. Here, we demonstrate an EC-free high salt concentration electrolyte (HSCE) composed of lithium bis(fluorosulfonyl)imide salt and tetrahydrofuran solvent with enhanced subzero temperature operation originating from unusually rapid low-temperature Li⁺ transport. Experimental and theoretical characterizations reveal the dominance of intra-aggregate ion transport in the HSCE that enables efficient low-temperature transport by increasing the exchange rate of solvating counterions relative to that of solvent molecules. This electrolyte also produces a <5 nm thick anion-derived LiF-rich SEI layer with excellent graphite electrode compatibility and electrochemical performance at subzero temperature in half-cells. Full cells based on LiNi_0.6Co_0.2Mn_0.2O₂||graphite with tailored HSCE electrolytes outperform state-of-the-art cells comprising conventional EC electrolytes during charge-discharge operation at an extreme temperature of -40 °C. These results demonstrate the opportunities for creating intrinsically robust low-temperature Li⁺ technology.

Molecular Dynamics of Ionic Liquids from Fast-Field Cycling NMR and Molecular Dynamics Simulations

Understanding the connection between the molecular structure of ionic liquids and their properties is of paramount importance for practical applications. However, this connection can only be established if a broad range of physicochemical properties on different length and time scales is already available. Even then, the interpretation of the results often remains ambiguous due to the natural limits of experimental approaches. Here we use fast-field cycling (FFC) to access both translational and rotational dynamics of ionic liquids. These combined with a comprehensive physicochemical characterization and MD simulations provide a toolkit to give insight into the mechanisms of molecular mechanics. The FFC results are consistent with the computer simulation and conventional physicochemical approaches. We show that curling of the side chains around the positively charged cationic core is essential for the properties of ether-functionalized ionic liquids, and we demonstrate that neither geometry nor polarity alone are sufficient to explain the macroscopic properties.

Ferroelectric nematic liquid-crystalline phases

Recent experimental realization of ferroelectric nematic liquid crystalline phases stimulated material development and numerous experimental studies of these phases, guided by their fundamental and applicative interest. In this Perspective, we give an overview of this emerging field by linking history and theoretical predictions to a general outlook of the development and properties of the materials exhibiting ferroelectric nematic phases. We will highlight the most relevant observations to date, e.g., giant dielectric permittivity values, polarization values an order of magnitude larger than in classical ferroelectric liquid crystals, and nonlinear optical coefficients comparable with several ferroelectric solid materials. Key observations of anchoring and electro-optic behavior will also be examined. The collected contributions lead to a final discussion on open challenges in materials development, theoretical description, experimental explorations, and possible applications of the ferroelectric phases.

Toward Bottom-Up Understanding of Transport in Concentrated Battery Electrolytes

Bottom-up understanding of transport describes how molecular changes alter species concentrations and electrolyte voltage drops in operating batteries. Such an understanding is essential to predictively design electrolytes for desired transport behavior. We herein advocate building a structure-property-performance relationship as a systematic approach to accurate bottom-up understanding. To ensure generalization across salt concentrations as well as different electrolyte types and cell configurations, the property-performance relation must be described using Newman's concentrated solution theory. It uses Stefan-Maxwell diffusivity, _ij , to describe the role of molecular motions at the continuum scale. The key challenge is to connect _ij to the structure. We discuss existing methods for making such a connection, their peculiarities, and future directions to advance our understanding of electrolyte transport.

Coarse-Grained Water Model Development for Accurate Dynamics and Structure Prediction

Several coarse-graining (CG) methods have been combined to develop a CG model of water capable of the accurate prediction of structure and dynamics properties. The multiscale coarse-graining (MS-CG) method based on force matching and the PDF-based coarse-graining method were used for accurate dynamics prediction. The iterative Boltzmann inversion (IBI) method was added for accurate structure representation. The approach is applied to bulk water, and the results show close reproduction of the CG structure when compared with the reference atomistic data. The combination of MS-CG and IBI methods facilitates the development of CG force fields at different temperatures based on a single MS-CG coarse-graining procedure. The dynamic properties of the CG water model closely match those obtained from the reference atomistic system. The general application of this approach to any existing coarse-graining methods is discussed.

Nonpolarizable Force Fields through the Self-Consistent Modeling Scheme with MD and DFT Methods: From Ionic Liquids to Self-Assembled Ionic Liquid Crystals

A key to achieve the accuracy of molecular dynamics (MD) simulation is the set of force fields used to express the atomistic interactions. In particular, the electrostatic interaction remains the main issue for the precise simulation of various ionic soft materials from ionic liquids to their supramolecular compounds. In this study, we test the nonpolarizable force fields of ionic liquids (ILs) and self-assembled ionic liquid crystals (ILCs) for which the intermolecular charge transfer and intramolecular polarization are significant. The self-consistent modeling scheme is adopted to refine the atomic charges of ionic species in a condensed state through the use of density functional theory (DFT) under the periodic boundary condition. The atomic charges of the generalized amber force field (GAFF) are effectively updated to express the electrostatic properties of ionic molecules obtained by the DFT calculation in condensed phase, which improves the prediction accuracy of ionic conductivity with the obtained force field (GAFF-DFT). The derived DFT charges then suggest that the substitution of a hydrophobic liquid-crystalline moiety into IL-based cations enhances the charge localization of ionic groups in the amphiphilic molecules, leading to the amplification of the electrostatic interactions among the hydrophilic/ionic groups in the presence of hydrophobic moieties. In addition, we focus on an ion-conductive pathway hidden in the self-assembled nanostructure. The MD results indicate that the ionic groups of cation and anion interact strongly for keeping the bicontinuous nanosegregation of ionic nanochannel. The partial fractions of hydrophilic/ionic and hydrophobic nanodomains are then quantified with the volume difference from referenced IL systems, while the calculated ionic conductivity decreases in the self-assembled ILCs more than the occupied volume of ionic nanodomains. These analyses suggest that the mobility of ions in the self-assembled ILCs remains quite restricted even with small tetrafluoroborate anions because of strong attractive interaction among ionic moieties.

Aqueous Electrolytes Reinforced by Mg and Ca Ions for Highly Reversible Fe Metal Batteries

Iron (Fe) metal batteries, such as Fe-ion batteries and all Fe flow batteries, are promising energy storage technologies for grid applications due to the extremely low cost of Fe and Fe salts. Nonetheless, the cycle life of Fe metal batteries is poor primarily due to the low Coulombic efficiency of the Fe deposition/stripping reaction. Current aqueous electrolytes based on Fe chloride or sulfate salts can only operate at a Coulombic efficiency of <91% under mild operation conditions (<5 mA/cm²), largely due to undesired hydrogen evolution reaction (HER). This work reports a series of novel Fe electrolytes, Fe electrolytes reinforced with Mg ions (FERMI) and Ca ions (FERCI), which have remarkably better Coulombic efficiency, higher conductivity, and faster deposition/stripping kinetics. By the addition of 4.5 M MgCl₂ or CaCl₂ into the baseline FeCl₂ electrolyte, the Fe deposition/stripping efficiency can be significantly improved to 99.1%, which greatly boosts the cycling performance of Fe metal batteries in both half-cells and full-cells. Mechanistic studies reveal that the remarkably improved efficiency is due to a reduced amount of "dead Fe" as well as suppressed HER. By the combination of experiments and molecular dynamics and density functional theory computation, the electrolyte structure is revealed, and the mechanism for enhanced water reduction resistance is elucidated. These novel electrolytes not only enable a highly reversible Fe metal anode for low-cost energy storage technologies but also have the potential to address the HER side reaction problem in other electrochemical technologies based on aqueous electrolytes, such as CO₂ reduction, NH₃ synthesis, etc.

Green Polymer Electrolytes Based on Polycaprolactones for Solid-State High-Voltage Lithium Metal Batteries

Solid polymer electrolytes (SPEs) have attracted considerable attention for high energy solid-state lithium metal batteries (LMBs). In this work, potentially ecofriendly, solid-state poly(ε-caprolactone) (PCL)-based star polymer electrolytes with cross-linked structures (xBt-PCL) are introduced that robustly cycle against LiNi_0.6 Mn_0.2 Co_0.2 O₂ (NMC622) composite cathodes, affording long-term stability even at higher current densities. Their superior features allow for sufficient suppression of dendritic lithium deposits, as monitored by ⁷ Li solid-state NMR. Advantageous electrolyte|electrode interfacial properties derived from cathode impregnation with 1.5 wt% PCL enable decent cell performance until up to 500 cycles at rates of 1C (60 °C), illustrating the high potential of PCL-based SPEs for application in high-voltage LMBs.

A Unified Approach to Derive Atomic Partial Charges and Polarizabilities of Ionic Liquids

We propose a unified approach to fit simultaneously a set of atomic partial charges and polarizabilities of the polarizable model against the ab initio electrostatic potential (ESP) and polarizability. The polarizable model is represented with interactive atomic dipoles with distance-dependent attenuation. For the polarizable model employed in this study, the internal electric field on the polarization sites is fully turned on, and thus allows self-induced dipoles, which persist even for an isolated molecule/ion. By such treatment, the contribution of ESP stems not only from the partial charges but also from the self-induced dipoles, and the atomic partial charges and polarizabilities can be fitted simultaneously against ESP in a unified manner. The fitting with 1-ethyl-3-methylimidazolium (EMIM⁺) and nitrate (NO₃^-), a prototypical organic cation and inorganic anion, respectively, that can form ionic liquid, demonstrates that allowance of the self-induced dipoles gives much better fitness. Moreover, test on the total dipole of an EMIM⁺/NO₃^- ion pair shows that the agreement with the ab initio dipole is also much improved for the polarizable model, which highlights the importance of the polarization effects of ionic liquids.

Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions

Polarizable molecular dynamics simulations are a fast progressing field in the scientific research of ionic liquids. The fundamentals of polarizable simulations, as well as their application to ionic liquids, were summarized in a review [Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, Jr., A. D.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem. Rev. 2019, 119, 7940–7995] in 2019. Since then, new methods to treat intermolecular interaction of induced dipoles in these highly charged systems were developed. This concerns the damping of these interactions and additional charge transfer as well as the prediction of ionic materials with ultrahigh refractive indices. In addition to the progress of the polarizable force fields, also thermostats and barostats for polarizable simulations evolved recently.

Applying Classical, Ab Initio, and Machine-Learning Molecular Dynamics Simulations to the Liquid Electrolyte for Rechargeable Batteries

Rechargeable batteries have become indispensable implements in our daily life and are considered a promising technology to construct sustainable energy systems in the future. The liquid electrolyte is one of the most important parts of a battery and is extremely critical in stabilizing the electrode-electrolyte interfaces and constructing safe and long-life-span batteries. Tremendous efforts have been devoted to developing new electrolyte solvents, salts, additives, and recipes, where molecular dynamics (MD) simulations play an increasingly important role in exploring electrolyte structures, physicochemical properties such as ionic conductivity, and interfacial reaction mechanisms. This review affords an overview of applying MD simulations in the study of liquid electrolytes for rechargeable batteries. First, the fundamentals and recent theoretical progress in three-class MD simulations are summarized, including classical, ab initio, and machine-learning MD simulations (section 2). Next, the application of MD simulations to the exploration of liquid electrolytes, including probing bulk and interfacial structures (section 3), deriving macroscopic properties such as ionic conductivity and dielectric constant of electrolytes (section 4), and revealing the electrode-electrolyte interfacial reaction mechanisms (section 5), are sequentially presented. Finally, a general conclusion and an insightful perspective on current challenges and future directions in applying MD simulations to liquid electrolytes are provided. Machine-learning technologies are highlighted to figure out these challenging issues facing MD simulations and electrolyte research and promote the rational design of advanced electrolytes for next-generation rechargeable batteries.

Structure of water-in-salt and water-in-bisalt electrolytes

We report a systematic diffraction study of two "water-in-salt" electrolytes and a "water-in-bisalt" electrolyte combining high-energy X-ray diffraction (HEXRD) with polarized and unpolarized neutron diffraction (ND) on both H₂O and D₂O solutions. The measurements provide three independent combinations of correlations between the different pairs of atom types that reveal the short- and intermediate-range order in considerable detail. The ND interference functions show pronounced peaks around a scattering vector Q ∼ 0.5 Å^-1 that change dramatically with composition, indicating significant rearrangements of the water network on a length scale around 12 Å. The experimental results are compared with two sets of Molecular Dynamics (MD) simulations, one including polarization effects and the other based on a non-polarizable force field. The two simulations reproduce the general shapes of the experimental structure factors and their changes with concentration, but differ in many detailed respects, suggesting ways in which their force fields might be modified to better represent the actual systems.

Porous Liquids: Computational Design for Targeted Gas Adsorption

In this Perspective, we present the unique gas adsorption capabilities of porous liquids (PLs) and the value of complex computational methods in the design of PL compositions. Traditionally, liquids only contain transient pore space between molecules that limit long-term gas capture. However, PLs are stable fluids that that contain permanent porosity due to the combination of a rigid porous host structure and a solvent. PLs exhibit remarkable adsorption and separation properties, including increased solubility and selectivity. The unique gas adsorption properties of PLs are based on their structure, which exhibits multiple gas binding sites in the pore and on the cage surface, varying binding mechanisms including hydrogen-bonding and π-π interactions, and selective diffusion in the solvent. Tunable PL compositions will require fundamental investigations of competitive gas binding mechanisms, thermal effects on binding site stability, and the role of nanoconfinement on gas and solvent diffusion that can be accelerated through molecular modeling. With these new insights PLs promise to be an exceptional material class with tunable properties for targeted gas adsorption.

Transference Number in Polymer Electrolytes: Mind the Reference-Frame Gap

The transport coefficients, in particular the transference number, of electrolyte solutions are important design parameters for electrochemical energy storage devices. The recent observation of negative transference numbers in PEO–LiTFSI under certain conditions has generated much discussion about its molecular origins, by both experimental and theoretical means. However, one overlooked factor in these efforts is the importance of the reference frame (RF). This creates a non-negligible gap when comparing experiment and simulation because the fluxes in the experimental measurements of transport coefficients and in the linear response theory used in the molecular dynamics simulation are defined in different RFs. In this work, we show that, by applying a proper RF transformation, a much improved agreement between experimental and simulation results can be achieved. Moreover, it is revealed that the anion mass and the anion–anion correlation, rather than ion aggregates, play a crucial role for the reported negative transference numbers.

Emulating proton transfer reactions in the pseudo-protic ionic liquid 1-methylimidazolium acetate

Proton transfer reactions can enhance conductivity in protic ionic liquids. However, several proton reactions are possible in a multicomponent system of charged and neutral species, resulting in a complex reaction network. Probabilities and equilibrium concentrations of the involved species are modeled by the combination of reducible Markov chains and quantum-mechanical calculations.

Proton transfer reactions can enhance conductivity in protic ionic liquids.

GenX in water: Interactions and self-assembly

2,3,3,3-tetrafluoro-2-(heptafluoropropoxy) propanoate, a.k.a. "GenX", is a surfactant introduced as a safer alternative to replace perfluorooctanoate (PFOA) in the manufacturing of fluorinated polymers, however, GenX is shown to cause adverse health effects similar to, or even worse than, those of the legacy PFOA. With an overarching goal to understand the behavior of GenX molecules in aqueous media, we report here on GenX micelle formation and structure in aqueous solutions, on the basis of results obtained from a combination of experimental techniques such as surface tension, fluorescence, viscosity, and small-angle neutron scattering (SANS), and molecular dynamics (MD) simulations. To our best knowledge, this is the first report on GenX micelles. The critical micelle concentration (CMC) of GenX ammonium salt in water is 175 mM. GenX forms small micelles with association number 6-8 and 10 Å radius. GenX molecules prefer to align along the micelle surface, and the ether oxygen of GenX has very little interaction with and exposure to water. Information on the surfactant and interfacial properties of GenX is crucial, since such properties are manifestations of interactions between GenX molecules and between GenX and water molecules and, in turn, the amphiphilic character of GenX dictates its fate and transport in the aqueous environment, its interactions with various biomolecules, and its binding to adsorbent materials.

Ionic Conductivity Enhancement in UHMW PEO Gel Electrolytes Based on Room-Temperature Ionic Liquids and Deep Eutectic Solvents

Physical gels made of poly(ethylene oxide) (PEO) and deep eutectic solvents urea-Li bis(trifluoromethanesulfonyl)imide (TFSI) and ethylene glycol/LiTFSI, or pyrrolidinium ionic liquid solutions PYR13TFSI-LiTFSI and PYR14TFSI-LiTFSI, are prepared by a fast, single-step process, which involves no auxiliary solvents or intermediates and is reproducible and scalable. The properties of these gels are studied as a function of the PEO content and its molecular weight and the nature of the liquid electrolyte. The gels prepared with a low concentration (1-5 wt %) of ultrahigh molecular weight (UHMW) PEO are tough, stretchable materials which resemble soft elastomers and are also self-healing and transparent. Their rheology shows the conventional behavior of physical polymer gels, so that the higher the molecular weight of PEO, the lower the polymer concentration needed to produce the gel. However, the ion conductivities and diffusivities of the gels are striking, in many cases being equal to or significantly higher than those of pure liquid electrolytes. This ion conductivity enhancement is the highest for the lowest PEO concentration with the highest molecular weight. This unprecedented molecular weight dependence of conductivity and diffusivity is the result of two combined effects: the liquid electrolyte chemical structure modification as a consequence of the addition of PEO and the development of elastic networks, where ion mobility and rheology are uncoupled when the PEO added is of UHMW.

Influence of water on the microscopic dynamics of 1-butyl-3-methylimidazolium tetrafluoroborate studied by means of quasielastic neutron scattering

We present a systematic study on the effect of water on the microscopic dynamics of 1-butyl-3-methylimidazolium tetrafluoroborate by means of quasielastic neutron scattering. By mixing the ionic liquid with either heavy or light water, the different contributions to the quasielastic broadening could be identified and treated separately. This study was performed at room temperature, which is more than 15 °C above the demixing line. Our results show that even small amounts of water accelerate the diffusion mechanisms considerably. While samples with small water percentage reveal a diffusion process confined within ionic liquid nanodomains, an admixture of more than 15 wt. % water relieves the confinement. Furthermore, the presence of two water species was identified: one behaving as free water, whereas the other was interpreted as a component bound to the ionic liquid motion. Based on the fact that water preferentially binds to the BF4 anion, which itself has a negligible contribution to the scattered intensity, our experiments reveal unprecedented information about the microscopic anion dynamics.

Coarse-Grained Dynamically Accurate Simulations of Ionic Liquids: [pyr14][TFSI] and [EMIM][BF4]

In this work, coarse-grained (CG) models for two different sets of ionic liquids were developed from atomistic molecular dynamics (MD) reference systems, expanding their system size and time duration capabilities. The bonded force field of the CG systems was built using harmonic oscillator potential (HOP) fitting, while the nonbonded force field was generated with the multiscale coarse-graining (MS-CG) approach based on force matching. The dynamics of each system were corrected using the probability distribution function-based coarse-grained molecular dynamics (PDF-based CGMD) method. The structure and dynamics of each system were proven to match reference system data at two temperature scales. CG models and force fields for these liquids were developed to exemplify a general purpose methodology for producing MD results of ionic liquids and other fluids with accurate structural as well as dynamic properties. As an application, developed ionic liquids CG models were then applied to study vacuum-interface interaction. Density profile results of vacuum-interface exposure show significant deviation from bulk behavior. At the interface, multilayer ordering of ionic liquids is predicted to be similar to those observed from an experimental work. This ordering is intensified by decreasing temperature and use of the PDF-based CGMD method as opposed to conventional CG methods.

Charge transfer and polarisability in ionic liquids: a case study

The practical use of ionic liquids (ILs) is benefiting from a growing understanding of the underpinning structural and dynamic properties, facilitated through classical molecular dynamics (MD) simulations. The predictive and explanatory power of a classical MD simulation is inextricably linked to the underlying force field. A key aspect of the forcefield for ILs is the ability to recover charge based interactions. Our focus in this paper is on the description and recovery of charge transfer and polarisability effects, demonstrated through MD simulations of the widely used 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C₄C₁im][NTf₂] IL. We study the charge distributions generated by a range of ab initio methods, and present an interpolation method for determining atom-wise scaled partial charges. Two novel methods for determining the mean field (total) charge transfer from anion to cation are presented. The impact of using different charge models and different partial charge scaling (unscaled, uniformly scaled, atom-wise scaled) are compared to fully polarisable simulations. We study a range of Drude particle explicitly polarisable potentials and shed light on the performance of current approaches to counter known problems. It is demonstrated that small changes in the charge description and MD methodology can have a significant impact; biasing some properties, while leaving others unaffected within the structural and dynamic domains.

An automated framework for high-throughput predictions of NMR chemical shifts within liquid solutions

Identifying stable speciation in multi-component liquid solutions is fundamentally important to areas from electrochemistry to organic chemistry and biomolecular systems. Here we introduce a fully automated, high-throughput computational framework for the accurate prediction of stable species in liquid solutions by computing the nuclear magnetic resonance (NMR) chemical shifts. The framework automatically extracts and categorizes hundreds of thousands of atomic clusters from classical molecular dynamics simulations, identifies the most stable species in solution and calculates their NMR chemical shifts via density functional theory calculations. Additionally, the framework creates a database of computed chemical shifts for liquid solutions across a wide chemical and parameter space. We compare our computational results to experimental measurements for magnesium bis(trifluoromethanesulfonyl)imide Mg(TFSI)2 salt in dimethoxyethane solvent. Our analysis of the Mg2+ solvation structural evolutions reveals key factors that influence the accuracy of NMR chemical shift predictions in liquid solutions. Furthermore, we show how the framework reduces the performance of over 300 13C and 600 1H density functional theory chemical shift predictions to a single submission procedure.

Operando Characterization of Organic Mixed Ionic/Electronic Conducting Materials

Operando characterization plays an important role in revealing the structure-property relationships of organic mixed ionic/electronic conductors (OMIECs), enabling the direct observation of dynamic changes during device operation and thus guiding the development of new materials. This review focuses on the application of different operando characterization techniques in the study of OMIECs, highlighting the time-dependent and bias-dependent structure, composition, and morphology information extracted from these techniques. We first illustrate the needs, requirements, and challenges of operando characterization then provide an overview of relevant experimental techniques, including spectroscopy, scattering, microbalance, microprobe, and electron microscopy. We also compare different in silico methods and discuss the interplay of these computational methods with experimental techniques. Finally, we provide an outlook on the future development of operando for OMIEC-based devices and look toward multimodal operando techniques for more comprehensive and accurate description of OMIECs.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.