Search NASA⌕ Search

SEARCH · Search NASA

Results for “solvates”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Molecular simulation of the separation of toluene and p-xylene with the thermally-robust ionic liquid triphenyl-p-phenyl sulfonyl phenyl phosphonium

In this work, we use molecular dynamics simulations to evaluate the behavior of toluene (TOL) and p-xylene (XYL) within three different solvents: n -hexane (HEX), n -heptane (HEP), and a thermally-robust ionic liquid (triphenyl-p-phenyl sulfonyl phenyl phosphonium + bis(trifluoromethylsulfonyl)imide (TPSP + Tf 2 N)). Several different temperatures are explored, while the solvation structures, energetics, and dynamics of TOL and XYL are evaluated. The solvation free energy (ΔG solv ) and transfer free energy (ΔG transfer ) of TOL and XYL are calculated using thermodynamic integration (TI). The transfer free energy of XYL from HEP/HEX to TPSP + Tf 2 N is more favorable than that of TOL over the entire temperature range, while the transfer of both TOL and XYL becomes less favorable as the temperature increases. The diffusion rates of TOL and XYL are lower in the ionic liquid when compared to HEP and HEX, but the thermal stability of TPSP + Tf 2 N may allow for higher operating temperatures and accelerated mass transport rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Polyacrylonitrile-based electrolytes: How processing and residual solvent affect ion transport and stability

Polyacrylonitrile (PAN) is one of the alternative candidate polymer hosts to form solid polymer electrolytes (SPEs) besides the widely used poly(ethylene oxide). In this study, we systematically investigate the processing of PAN based SPEs containing lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt, using dimethylformamide (DMF) as the solvent. The effects of PAN processing procedure including solution mixing, casting, and drying on the morphology, ion transport, solvation structure, Li and oxidative stability of PAN electrolytes are thoroughly examined. In particular, four drying conditions are investigated and the amount of residual DMF is accurately determined using infrared (IR) spectroscopy. Varying the drying conditions can lead to five orders of magnitude decrease in the ionic conductivity. As DMF content decreases, the SPE's stability again Li metal dramatically improves. The practical oxidative stability is also strongly affected by the residual DMF content, ranging from 2.5 V to 3.5 V, much lower than reported values. Finally, the role of the residual DMF solvent is elucidated. DMF content vitally influences ion solvation structure at different concentration regimes, which ultimately dictates the ion conduction mechanism and oxidative stability of PAN based SPEs. In conclusion, this thorough study lays the groundwork for future development of PAN based electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Calorimetry can detect the early onset of hydrolysis in hybrid supercapacitors with aqueous electrolytes

This study investigates the effect of cation species on the onset of electrolyte hydrolysis in hybrid supercapacitors with aqueous electrolytes using isothermal operando calorimetry. The cells consisted of a positive α-MnO 2 cryptomelane electrode and a negative activated carbon (AC) electrode with either 0.5 M K 2 SO 4 or 0.5 M Cs 2 SO 4 aqueous electrolytes. They were characterized using cyclic voltammetry and galvanostatic cycling. In addition, the instantaneous heat generation rate at each electrode was measured using a custom isothermal operando calorimeter. Heat generation associated with resistive losses (Joule heating) and reversible ion adsorption/desorption was clearly identified. For larger potential windows, an endothermic dip, attributed to the onset of hydrolysis, was observed at the positive α-MnO 2 electrode where K + and Cs + ions engaged in fast surface redox reactions. Interestingly, this endothermic dip appeared at 1.8 V and 2.0 V for K 2 SO 4 and Cs 2 SO 4 aqueous electrolytes, respectively. The difference in the stable operating potential window was attributed to thinner solvation shell around Cs + cation than for K + thus reducing the amount of water present near the electrodes as ions partially shed their solvation shells during adsorption. The early onset of hydrolysis could be observed by isothermal operando calorimetry before it could be observed with conventional electrochemical methods.

42 ENGINEERING↗
DOE OSTI2020

Ultrahigh coulombic efficiency electrolyte enables Li||SPAN batteries with superior cycling performance

Raising the coulombic efficiency of lithium metal anode cycling is the deciding step in realizing long-life rechargeable lithium batteries. Here, we designed a highly concentrated salt/ether electrolyte diluted in a fluorinated ether: 1.8 M LiFSI in DEE/BTFE (diethyl ether/bis(2,2,2-trifluoroethyl)ether), which realized an average coulombic efficiency of 99.37% at 0.5 mA cm -2 and 1 mAh cm -2 for more than 900 cycles. This electrolyte also maintained a record coulombic efficiency of 98.7% at 10 mA cm -2 , indicative of its ability to provide fast-charging with high cathode loadings. Morphological studies reveal dense, dendrite free Li depositions after prolonged cycling, while surface analyses confirmed the formation of a robust LiF-rich SEI layer on the cycled Li surface. Moreover, we discovered that this ether-based electrolyte is highly compatible with the low-cost, high-capacity SPAN (Sulfurized polyacrylonitrile) cathode, where the constructed Li||SPAN cell exhibited reversible cathode capacity of 579 mAh g -1 and no capacity decay after 1200 cycles. A cell where a high areal loading SPAN electrode (>3.5 mAh cm -2 ) is paired with only onefold excess Li was constructed and cycled at 1.75 mA cm -2 , maintaining a coulombic efficiency of 99.30% for the lithium metal. Computational simulations revealed that at saturation, the Li-FSI complex forms contact ion pairs, with a first solvation shell comprising DEE molecules, and a second solvation shell with a mix of DEE/BTFE. This study provides a path to enable high energy density Li||SPAN batteries with stable cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Unexpected Hydroxide Ion Structure and Properties at Low Hydration

Understanding the behavior of hydroxide ions in aqueous and non-aqueous media is fundamental to many chemical, biological, and electrochemical processes. Research has primarily focused on a single fully solvated hydroxide ion, either as an isolated cluster or in bulk. This work presents the first computational study to consider hydroxide under low hydration levels in detail, where the anion may not be fully solvated. Under such conditions, we find that the anions are predominantly present as unique water-bridged hydroxide pair complexes, distinct from previously reported structures under fully hydrated conditions. Although similar hydroxide pair structures were previously reported, we analyze these structures for the first time in the disordered liquid state where they are found to be unusually stable in the presence of bulky quaternary ammonium cations. Our findings help explain the unusual diffusion behavior as well as the higher reactivity of hydroxide anions observed under low hydration conditions.

computational studies↗
DOE OSTI2022

Electronic Structure Calculations of the Fundamental Interactions in Solvent Extraction Desalination

Solvent extraction desalination (SED) is a promising membrane-free desalination technology, but new solvents are still needed that can improve the molecular-level water selectivity, while minimizing solvent contamination in the aqueous phase. Improvements in solvent design for SED require a fundamental understanding of the molecular interactions governing solvation behavior. Herein, quantum chemical calculations are used to delineate different fundamental interactions in water-solvent systems involving 38 different solvents, with a particular focus on the role of the electrostatic potential characteristics. Furthermore, these interactions are rigorously analyzed using the independent gradient model based on Hirshfeld partition (IGMH), Symmetry-Adapted Perturbation Theory (SAPT), as well as the energy decomposition analysis based on force field (EDA-FF) to compare against a standard force field model. The solvation free energies and partition coefficients are also calculated to estimate water extraction performance in a bulk solvent system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Striking stability of a mixed-valence thallium(III)-thallium(I) complex in some solvents

At the dissolution of solid anhydrous thallium(III) trifluoromethanesulfonate, Tl(CF 3 SO 3 ) 3 , or thallium(III) trifluoroacetate, Tl(CF 3 COO) 3 , in dimethylsulfoxide (dmso) or N,N,N’,N’-tetramethylurea (tmu), intensely red-colored complexes are formed. This red thallium complex is stable for years in dmso, while it is reduced fairly rapidly to thallium(I) in tmu with a half-life time of an hour. At the dissolution of Tl(CF 3 SO 3 ) 3 in N,N-dimethylpropyleneurea (dmpu) an immediate reduction to thallium(I) takes place. A stable colorless aqueous thallium(III) solution is obtained at the dissolution in acidic water. Stable dmso solutions and solid dmso solvates of thallium(III) perchlorate, nitrate and trifluormethanesulfonate can be prepared by adding dmso to concentrated acidic aqueous thallium(III) solutions. These experimental observations conclude that the pure solids Tl(CF 3 SO 3 ) 3 and Tl(CF 3 COO) 3 play an essential role in the formation of the red-colored thallium complexes. 205 Tl NMR data show that the red thallium complex contains equal amounts of thallium(III) and thallium(I). The structure of the red thallium complex in dmso, as determined by EXAFS, has Tl–O bond distances of 2.216(3) and 2.80(2) Å, which are in very close agreement with the bond distances obtained in the pure dmso solvates of the thallium(III) and thallium(I) ions, respectively, and a Tl···Tl distance of 3.49(1) Å bridged by oxygen atoms. From the EXAFS data it is impossible to distinguish if dmso molecules and/or trifluoromethanesulfonate ions act as bridges. DFT calculations could eliminate some structures due to the irrelevant structural parameters or the energetics of the proposed reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Choline Fluoride-Ethylene glycol deep eutectic solvent mixture – Synthesis and physicochemical properties

A critical factor in achieving tunability in DESs is varying the type of hydrogen bond acceptors (HBAs) and donors (HBDs). Choline halides (ChX), such as ChCl and ChBr, have been studied to investigate the effect of the choline anion on DESs. However, ChF-based DESs have been underrepresented in scientific research. ChF has excellent potential, such as high electron affinity, ionization potential, and strong hydrogen bond interactions. Moreover, ChF-based DES can lead to a deeper understanding of DESs due to the NMR-active 19F nuclear spin nature. We report the first synthesis and characterization of choline fluoride-EG (1:2, ChF:EG) DES, which we named “EthalineF.” We studied this mixture using various spectroscopic techniques, including NMR, FT-IR, Raman, and UV–vis spectroscopies. Polarity, density, viscosity, and conductivity measurements were collected as well. Transient absorption spectroscopy (TAS) was used to study the solvation dynamics of EthalineF. 19F NMR indicated that some fluoride anions (~15%) are coordinated to the choline cation. In EthalineF, the new type of observed fluoride anions is proposed to be a tight-binding choline-fluoride ion pair. Further, EthalineF has a larger solvodynamic radius than Ethaline. This can explain the lower conductivity for EthalineF compared to Ethaline despite the higher fluidity and solvation dynamics for EthalineF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

CO 2 adsorption mechanisms at the ZIF-8 interface in a Type 3 porous liquid

Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO 2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO 2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO 2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO 2 molecules were unable to diffuse into the ZIF-8 pore, with CO 2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO 2 adsorption should be based on the strength of gas binding at the solvated particle surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Boosting the sodium storage behaviors of carbon materials in ether-based electrolyte through the artificial manipulation of microstructure

The porous carbon blacks rationally designed by a facile yet efficient NH 3 thermal etching route have been investigated as anode materials in an ether-based electrolyte for sodium-ion batteries. The as-synthesized CBN35 carbon black with a 35% weight loss after NH3 thermal etching exhibited a large specific charge capacity of 352 mAh g -1 at 50 mA g -1 and a superior rate capability of 101 mAh g -1 at 16000 mA g -1 , due to its highest microporosity, an appropriate surface area, a desirable microstructure, and a promising hybrid intercalation mechanism. Impressively, even cycled at 1600 mA g -1 over 3200 cycles, an outstanding reversible capacity of 103 mAh g -1 with a negligible 0.0162% capacity loss per cycle can still be achieved. Based on the multimodal characterizations including the structural probes of phase evolution for carbon materials, the electrochemical techniques, and the surface-sensitive XAS measurements, the exceptional electrochemical properties should stem from several merits of modified carbon black system. While the particular microporous structure provides relatively more accessible sodium storage sites, a novel hybrid intercalation mechanism in ether-based electrolyte would incorporate the sodium ion insertion into the disordered structure with the solvated sodium ion species co-intercalation into the graphitic phase. In addition to the diffusion-controlled redox reactions, the noticeable surface-induced pseudocapacitive reactions also significantly contribute to the charge storage upon sodiation and guarantee the rapid migrations of sodium ions/solvated compounds. In conclusion, this system further features a controlled emergence of a robust but thin solid electrolyte interphase layer, which could suppress the side reactions of active electrode with reactive electrolyte, maintain the fragile porous structure upon cycling, and facilitate the migrations of sodium ions and solvated sodium ion compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Superior long-term cycling of high-voltage lithium-ion batteries enabled by single-solvent electrolyte

A new single-solvent electrolyte system comprising lithium bis(fluorosuflonyl) imide (LiFSI) and beta-fluorinated sulfone (TFPMS) was designed to enable very stable long-term cycling of high-voltage lithium-ion batteries. Compared to other fluorinated solvents such as alpha-fluorinated sulfone (FMES) and fluorinated carbonate (FEMC), which are prone to reduction on the graphite anode, the LiFSI-TFPMS system displayed outstanding compatibility with graphite. While regular carbonate and sulfone from the LiFSI electrolyte system are compatible with the graphite anode, their high solvating power not only induces severe corrosion on the aluminum cathode current collector at high voltage, but also renders a low aggregation level at a normal salt concentration (about 1.0 M), resulting in the formation of an unstable solid-electrolyte interphase (SEI) on the graphite anode. Owing to the low solvating power of TFPMS, the aggregation level of the LiFSI-TFPMS system is relatively high even at normal salt concentration, which not only facilitates the formation of a robust SEI by the sacrificial decomposition of LiFSI, but also suppresses the aluminum corrosion of the LiFSI electrolyte system at high voltage. Together with the high intrinsic anodic stability of TFPMS, the superior cycling performance of graphite parallel to LiNi 0.6 Co 0.2 Mn 0.2 O 2 cells was achieved by employing the non-flammable LiFSI-TFPMS single-solvent electrolyte system.

25 ENERGY STORAGE↗
DOE OSTI2021

Synergistic Coupling of CO2 and H2O during Expansion of Clays in Supercritical CO2-CH4 Fluid Mixtures

A combination of in operando IR and XRD methods were used to investigate the interaction of variably hydrated supercritical CO2-CH4 fluids with Na/NH4/Cs substituted montmorillonites. Comparing the behavior of Na-clay exposed to CH4 versus CO2 dominant fluid phases demonstrated that CO2 disrupts the H-bond network of intercalated H2O. At relatively low potentials of H2O, CO2 facilitates H2O intercalation and expansion of the interlayer region of the Na-clay. In contrast, CO2 inhibited H2O intercalation but promoted expansion of the Cs-clay. The NH4-clay displayed intermediated behavior. Methane likely intercalated the clays but opportunistically, filling unoccupied space after CO2 and H2O actively expanded the interlayer region. By comparing and contrasting the behavior of the Cs and NH4 clays to the Na Clay we conclude that, whereas H2O is required for initial Na-clay expansion, CO2 synergistically facilitated expansion and hydration of Na-clay by mixing with H2O (the entropic contribution) and participating in the outer solvation sphere of Na (the enthalpic contribution). In contrast, CO2 can directly solvate interlayer Cs and NH4 and cause expansion in the absence of H2O; nonetheless, H2O and CO2 also act synergistically to facilitate expansion of the Cs and NH4 clays at low H2O concentrations.

Loring, John S.↗
DOE OSTI2024

Pairing a Global Optimization Algorithm with EXAFS to Characterize Lanthanide Structure in Solution

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. Here, to accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Understanding the Interplay between the Nonvalence and Valence States of the Uracil Anion upon Monohydration

In this work we present an ab initio investigation into the effect of monohydration on the interaction of uracil with low energy electrons. Electron attachment and photodetachment experimental studies have previously shown dramatic changes in uracil upon solvation with even a single water molecule, due to an inversion of the character of the ground state of the anion. Here we explore the interplay between the nonvalence and valence states of the uracil anion, as a function of geometry and site of solvation. Our model provides unambiguous interpretation of previous photoelectron studies, reproducing the binding energies and photoelectron images for bare uracil and a single isomer of the U•(H 2 O) 1 cluster. Furthermore, the results of this study provide insight into how electrons may attach to hydrated nucleobases. These results lay the foundations for further investigations into the effect of microhydration on the electronic structure and electron capture dynamics of nucleobases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗
DOE OSTI2021

Sub-picosecond Production of Solute Radical Cations in Tetrahydrofuran after Radiolysis

Ultrafast hole transfer from solvent radical cations produced by radiolysis with ~10 ps, 9 MeV electron pulses to solutes was investigated in tetrahydrofuran (THF). Because of rapid fragmentation of initially produced THF +• , solute radical cations are not expected and have not previously been reported. Using 9,9-Dihexyl-2,7-dibromofluorene (Br 2 F) at 5 to 1000 mM, Br 2 F +• was observed with radiation chemical yields up to G = 2.23 / 100 eV absorbed. While more than half of this was the result of direct solute ionization, the results highlight the importance of capture of holes from THF +• prior to solvation and fragmentation. The observed data show a time-resolution limited (15 ps) rise in transient absorption of Br 2 F +• , identical in form to reports of pre-solvated or dry electron capture in water and a few organic liquids, including THF. The results were thus interpreted with a similar formalism, finding C 37 = 1.7 M, the concentration at which 37% of holes escape capture. The yield of solvent hole capture can be accounted for by the formation of solvent holes adjacent to solute molecules reacting faster than they can fragment, however mechanisms such as delocalized holes or rapid hopping may play a role. Finally, low temperature results find over two times more capture, supporting the speculation that if THF +• was longer lived, the yield of capture in under 15 ps would have been at least 2x larger at 1 M Br 2 F, possibly capturing nearly all available holes from the solvent.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗