Search NASA⌕ Search

SEARCH · Search NASA

Results for “solvation interactions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Comparative Study of Electrical Double Layer Effects for CO Reduction Reaction Kinetics

Solvation models describe how the interactions between the solutes and solvents affect the reactivity and selectivity in electrochemical processes. In this study, we developed a framework for evaluating the effects of applied potential and electrical double layer on CO reduction (COR), comparing fully explicit, implicit, and hybrid solvation models at the standard hydrogen electrode (SHE) scale. We analyzed all crucial intermediates leading to the production of C 1 and C 2 products and found good agreement across these models. Some notable differences were observed in the implicit description of *C and *CO adsorption at higher overpotentials and overall trends in the adsorption energies within the hybrid model. Using this unified SHE framework, we built comparable microkinetic models for COR kinetics and rates. Despite small differences in thermodynamics, solvation-model-based microkinetic simulations showed good agreement for onset potentials against the benchmark experiment. Only qualitative difference was observed for C 1+ versus hydrogen evolution at high overpotentials for the implicit model. Finally, we constructed a generalized C 2 selectivity map in descriptor space (U SHE , ΔG electrolyte CO ), which highlights the limitations of a copper-based COR catalyst and guides the search for optimal descriptor parameters to maximize C 2 selectivity. In conclusion, these findings demonstrate the importance of considering electrical double layer effects in reduction reactions and offer a useful framework for comparing solvation models and predicting optimal electrochemical conditions for specific applications.

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

Anti-Coulomb ion-ion interactions: A theoretical and computational study

The free energy of ion solvation can be decomposed into enthalpic and entropic contributions. This helps one to understand the connection between the dielectric properties and the underlying forces. We present a simple linear-response model of screened charge interactions that provides an alternative understanding of solvation barriers. Moreover, it explains the “anti-Coulomb” interactions (attraction between like-charged ions and repulsion between opposite-charged ions) observed in both simulations and experiments. We show that this is a universal behavior associated to the nonlocal response function of any dielectric or metallic systems. Published by the American Physical Society 2024

Wills, Alec (ORCID:0000000280184692)↗
DOE OSTI2023

Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear–Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics

Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear–electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. Furthermore, these calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Deciphering the Dynamic Balance Between Solvation Strength and Polysulfides Reaction Heterogeneity in Practical Lithium‐Sulfur Batteries

Achieving stable interfacial chemistry in lithium–sulfur batteries under practical conditions remains a key barrier to commercialization. Here, we demonstrate that interfacial dynamics can be effectively regulated by coupling solvation-power control with intrinsic heterogeneity of sulfur redox chemistry through the introduction of a weakly solvating fluorinated cosolvent, LIB 1200ET (1200ET). Compared with conventional fluorinated ethers, 1200ET efficiently shifts Li + solvation environment toward a more non-coordinated configuration at low volume fractions, enabling substantial solvation modulation without significantly impairing sulfur redox kinetics. This solvation transition weakens Li + –solvent interactions while strengthening Li + –anion and Li + –lithium polysulfide (LPS) coordination, suppressing LPS solubility and promoting reconstruction of solid–electrolyte interphase (SEI). Regulated LPS chemistry, together with 1200ET, leads to formation of a S 4+ -rich, LiF-reinforced SEI with enhanced ionic conductivity and mechanical robustness. Spatially resolved sulfur K-edge X-ray absorption spectroscopy on pouch cells reveals pronounced current-density-dependent chemical heterogeneity, distinguishing kinetically dominated and solvation-controlled regions. Under practical conditions (3.7 mg cm −2 sulfur loading, E/S = 6 µL mg −1 ), a single-layer pouch cell delivers 527 mAh g −1 over 200 cycles at C/3, while an Ah-level multilayer pouch cell achieves an energy density of 358 Wh kg −1 . These results establish non-coordinating cosolvent-driven solvation engineering as a scalable strategy for practical Li–S batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

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

Insights into the Solvation Characteristics of Zwitterionic Deep Eutectic Solvents Using Multiple Polarity Scales

Deep eutectic solvents (DESs) formed by mixing zwitterionic hydrogen bond acceptors (HBAs) with diol- and organic acid-based hydrogen bond donors (HBDs) are popular sustainable solvents. In addition to possessing an amphiphilic character, zwitterionic DESs offer unique solvation properties, unlike those composed of halide-containing HBAs. While all other popular classes of DESs have been investigated, no study has attempted to characterize the solvation characteristics of zwitterionic DESs. In this study, the solvation properties of betaine [Bet]-, sulfobetaine-, and l-carnitine [l-Car]-based DESs were thoroughly investigated for the first time using solvatochromic dye methods and inverse gas chromatography (IGC). Herein, the Nile red polarity scale revealed that DESs composed of organic acids were generally more polar than diol-based DESs. Kamlet–Taft solvent parameters revealed that glycerol-based DESs offered exceptionally high hydrogen bond donating ability compared to those composed of 1,4-butane diol and triethylene glycol. The Abraham solvation parameter model revealed hydrogen bond basicity to be the most dominant interaction, and it was generally higher for diol-based HBDs compared to organic acid-based HBDs. Additionally, sulfobetaine HBAs possessing longer alkyl chains generally offered higher dispersive-type interactions and hydrogen bond basicity compared to [Bet] and [l-Car] salts. This is the first study to combine and compare results from solvatochromic dye methods and IGC in an attempt to comprehensively explain and predict the solvation characteristics of zwitterionic DESs.

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

Desorption lifetimes and activation energies influencing gas–surface interactions and multiphase chemical kinetics

Abstract. Adsorption and desorption of gases on liquid or solid substrates are involved in multiphase processes and heterogeneous chemical reactions. The desorption energy (Edes0), which depends on the intermolecular forces between adsorbate and substrate, determines the residence time of chemical species at interfaces. We show how Edes0 and temperature influence the net uptake or release of gas species, the rates of surface–bulk exchange and surface or bulk reactions, and the equilibration timescales of gas–particle partitioning. Using literature data, we derive a parameterization to estimate Edes0 for a wide range of chemical species based on the molecular mass, polarizability, and oxygen-to-carbon ratio of the desorbing species independent of substrate-specific properties, which is possible because of the dominant role of the desorbing species' properties. Correlations between Edes0 and the enthalpies of vaporization and solvation are rooted in molecular interactions. The relation between Edes0 and desorption kinetics reflects the key role of interfacial exchange in multiphase processes. For small molecules and semi-volatile organics (VOC, IVOC, SVOC), Edes0 values around 10–100 kJ mol−1 correspond to desorption lifetimes around nanoseconds to days at room temperature. Even higher values up to years are obtained at low temperatures and for low volatile organic compounds (LVOC, ELVOC/ULVOC) relevant for secondary organic aerosols (SOA). Implications are discussed for SOA formation, gas–particle partitioning, organic phase changes, and indoor surface chemistry. We expect these insights to advance the mechanistic and kinetic understanding of multiphase processes in atmospheric and environmental physical chemistry, aerosol science, materials science, and chemical engineering.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Moving beyond the Solvent-Tip Approximation to Determine Site-Specific Variations of Interfacial Water Structure through 3D Force Microscopy

Although interfacial solution structure impacts environmental, biological, and technological phenomena, including colloidal stability, protein assembly, heterogeneous nucleation, and water desalination, its molecular details remain poorly understood. Here, we visualize the three-dimensional (3D) hydration structure at the boehmite(010)-water interface using fast force mapping (FFM). Using a self-consistent scheme to decouple long-range tip-sample interactions from short-range solvation forces, we obtain the solution structure with lattice resolution. The results are benchmarked against molecular dynamics simulations that explicitly include the effects of the tip with different levels of approximation and systematically account for tip size, chemistry, and confinement effects. We find four laterally structured water layers within 1 nm of the surface, with the highest water densities at sites adjacent to hydroxyl groups. The key features beyond the first two layers can only be predicted using a full-scale simulation of the boehmite-water-silica system. Our findings further reveal a complex relationship between site-specific chemistry, water density, and long-range particle interactions; and present important advances toward quantitative data interpretation in 3D FFM.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electrolyte Reactivity on the MgV 2 O 4 Cathode Surface

Predictive understanding of the solvation-dependent reactivity and molecular interaction of electrolyte ions and solvent molecules on reactive electrodes has been a major challenge but is essential for addressing instabilities and surface passivation that occur at electrode-electrolyte interface (EEI) of multivalent Mg batteries. In this work, the isolated intrinsic reactivities of prominent chemical species present in magnesium bis(trifluoromethanesulfonimide) (Mg(TFSI) 2 ) in diglyme (G2) electrolytes, including ionic (TFSI - , [Mg(TFSI)] + , [Mg(TFSI):G2] + , [Mg(TFSI):2G2] + ) as well as neutral molecules (G2) on magnesium vanadate cathode (MgV 2 O 4 ) surface has been studied using a combination of first-principles calculations and multimodal analysis of well-defined cathode electrolyte interphase (CEI) layers. Here, our calculations show that non-solvated [Mg(TFSI)] + is the strongest adsorbing species on the MgV 2 O 4 surface compared to all other ions while fully solvated [Mg(TFSI):2G2] + are least favorable to decomposition. The cleavage of C-S bonds in TFSI - to form CF 3 - is predicted to be most desired pathway for all ionic species, which is followed by the cleavage of C-O bonds of G2 to yield CH 3 + or OCH 3 - species. The strong stabilization and electron transfer between ionic electrolyte species and MgV 2 O 4 is found to significantly favor these decomposition reactions on the surface compared to intrinsic gas phase dissociation. Experimentally, we used state-of-the-art ion soft landing to selectively deposit mass-selected TFSI - , [Mg(TFSI):G2] + and [Mg(TFSI):2G2] + on MgV 2 O 4 thin film to form well-defined electrolyte-MgV 2 O 4 interface. Analysis of soft-landed interphase using X-ray photoelectron, X-ray absorption near edge structure, electron energy-loss spectroscopies as well as transmission electron microscopy confirmed the presence of decomposition species (e.g., MgFx, carbonates) formed in the interfacial region and the higher amount of MgFx with [Mg(TFSI):G2] + , which corroborates the theoretical observation. Overall, we established the mechanistic pathway for the electrolyte-induced formation of passivating fluorides on MgV 2 O 4 cathode facilitated by the surface adsorption and charge transfer, which provided essential knowledge for rational design of stable electrolytes for multivalent cathodes.

Cathode-electrolyte interphase formation↗
DOE OSTI2024

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

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

HCOO - aq Degradation in Droplets by OH aq in an Atmospheric Pressure Glow Discharge

Plasmas in contact with liquids can degrade organic molecules in solution, as reactive oxygen and nitrogen species produced in the plasma solvate into the liquid. Immersing small droplets (tens of microns in diameter) in the plasma can more rapidly activate the liquid compared to treating a large volume of liquid with a smaller surface-to-volume ratio. The interactions between a radio frequency glow discharge sustained in He/H 2 O and a water droplet containing formate (HCOO - aq ) immersed in and flowing through the plasma were modeled using a 0-dimensional global plasma chemistry model to investigate these activation processes. HCOO - aq interacts with OH aq , produced from solvation of OH from the gas phase. HCOO - aq concentrations were benchmarked with previously reported experimental measurements. The diameter of the droplet, initial HCOO - aq concentration, and gas flow rate affect only the HCOO - aq concentration and OH aq density, leaving the OH density in the gas phase unaffected. Power deposition and gas mixture (e.g., percentage of H 2 O) change both the gas and liquid phase chemistry. A general trend was observed – during the first portion of droplet exposure to the plasma, OH aq primarily consumes HCOO - aq . However, O2 - aq , a byproduct of HCOO - aq consumption, consumes OH aq once O2 - aq reaches a critically large density. Using HCOO - aq as a surrogate for OH aq sensitive contaminants, combinations of residence time, droplet diameter, water vapor density, and power will determine the optimum remediation strategy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

N -Arylimide Molecular Balances: A Comprehensive Platform for Studying Aromatic Interactions in Solution

Noncovalent interactions of aromatic surfaces play a key role in many biological processes and in determining the properties and utility of synthetic materials, sensors, and catalysts. However, the study of aromatic interactions has been challenging because these interactions are usually very weak and their trends are modulated by many factors such as structural, electronic, steric, and solvent effects. Recently, N-arylimide molecular balances have emerged as highly versatile and effective platforms for studying aromatic interactions in solution. These molecular balances can accurately measure weak noncovalent interactions in solution via their influence on the folded–unfolded conformational equilibrium. The structure (i.e., size, shape, π-conjugation, and substitution) and nature (i.e., element, charge, and polarity) of the π-surfaces and interacting groups can be readily varied, enabling the study of a wide range of aromatic interactions. These include aromatic stacking, heterocyclic aromatic stacking, and alkyl–π, chalcogen–π, silver–π, halogen–π, substituent–π, and solvent–π interactions. The ability to measure a diverse array of aromatic interactions within a single model system provides a unique perspective and insights as the interaction energies, stability trends, and solvent effects for different types of interactions can be directly compared. Some broad conclusions that have emerged from this comprehensive analysis include: (1) The strongest aromatic interactions involve groups with positive charges such as pyridinium and metal ions which interact with the electrostatically negative π-face of the aromatic surface via cation–π or metal–π interactions. Attractive electrostatic interactions can also form between aromatic surfaces and groups with partial positive charges. (2) Electrostatic interactions involving aromatic surfaces can be switched from repulsive to attractive using electron-withdrawing substituents or heterocycles. These electrostatic trends appear to span many types of aromatic interactions involving a polar group interacting with a π-surface such as halogen–π, chalcogen–π, and carbonyl–π. (3) Nonpolar groups form weak but measurable stabilizing interactions with aromatic surfaces in organic solvents due to favorable dispersion and/or solvophobic effects. Furthermore, a good predictor of the interaction strength is provided by the change in solvent-accessible surface area. (4) Solvent effects modulate the aromatic interactions in the forms of solvophobic effects and competitive solvation, which can be modeled using solvent cohesion density and specific solvent–solute interactions.

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

Using a Chromatographic Pseudophase Model To Elucidate the Mechanism of Olefin Separation by Silver(I) Ions in Ionic Liquids

Silver(I) ions undergo selective olefin complexation and have been utilized in various olefin/paraffin separation techniques such as argentation chromatography and facilitated transport membranes. Ionic liquids (ILs) are solvents known for their low vapor pressure, high thermal stability, low melting points, and ability to promote a favorable solvation environment for silver(I) ion–olefin interactions. To develop highly selective separation systems, a fundamental understanding of analyte partitioning to the stationary phase and the thermodynamic driving forces behind solvation is required. In this study, a chromatographic model treating silver(I) ions as a pseudophase is constructed and employed for the first time to investigate the olefin separation mechanism in silver(I) salt/IL mixtures. Stationary phases containing varying amounts of noncoordinated silver(I) salt ([Ag + ][NTf 2 – ]) dissolved in the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 – ]) IL are utilized to determine the partition coefficients of various analytes including alkanes, alkenes, alkynes, aromatics, aldehyde, esters, and ketones. As ligand coordination to silver(I) ions is known to lower its olefin complexation capability, this study also examines two different types of coordinated silver(I) ion pseudophases, namely, monocoordinated silver(I) salt ([Ag + (1-decyl-2-methylimidazole, DMIM)][NTf 2 – ]) and dicoordinated silver(I) salt ([Ag + (1-methylimidazole, MIM)(DMIM)][NTf 2 – ]). The extent of olefin partitioning to the coordinated silver(I) ion pseudophases over the carrier gas and IL decreased by up to two orders of magnitude. Values for enthalpy, entropy, and free energy of solvation were determined for the three silver(I) ion-containing systems. Olefin retention was observed to be enthalpically dominated, while ligand coordination to the silver(I) ion pseudophase resulted in variations for both enthalpic and entropic contributions to the free energy of solvation. Furthermore, the developed model can be used to study chemical changes that occur in silver(I) ions over time as well as identify optimal silver(I) salt/IL mixtures that yield high olefin selectivity.

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

Dynamic Interconversion of Metal Active Site Ensembles in Zeolite Catalysis

Catalysis science is founded on understanding the structure, number, and reactivity of active sites. Kinetic models that consider active sites to be static and noninteracting entities are routinely successful in describing the behavior of heterogeneous catalysts. Yet, active site ensembles often restructure in response to their external environment and even during steady-state catalytic turnover, sometimes requiring non-mean-field kinetic treatments to describe distance-dependent interactions among sites. Such behavior is being recognized more frequently in modern catalysis research, with the advent of experimental methods to quantify turnover rates with increasing precision, an expanding arsenal of operando characterization tools, and computational descriptions of atomic structure and motion at chemical potentials and timescales increasingly relevant to reaction conditions. This review focuses on dynamic changes to metal active site ensembles on zeolite supports, which are silica-based crystalline materials substituted with Al that generate binding sites for isolated and low-nuclearity metal site ensembles. Metal sites can become solvated and mobilized during reaction, facilitating interactions among sites that change their nuclearity and function. Such intersite communication can be regulated by the zeolite support, resulting in non-single-site and potentially non-mean-field kinetic behavior arising from mechanisms of catalytic action that combine elements of those canonically associated with homogeneous and heterogeneous catalysis. Here, we discuss recent literature examples that document dynamic active site behavior in metal-zeolites and outline methodologies to identify and interpret such behavior. We conclude with our outlook on future research directions to develop this evolving branch of catalysis science and harness it for practical applications.

operando characterization↗
DOE OSTI2020

Dependency of solvation effects on metal identity in surface reactions

Solvent interactions with adsorbed moieties involved in surface reactions are often believed to be similar for different metal surfaces. However, solvents alter the electronic structures of surface atoms, which in turn affects their interaction with adsorbed moieties. To reveal the importance of metal identity on aqueous solvent effects in heterogeneous catalysis, we studied solvent effects on the activation free energies of the O–H and C–H bond cleavages of ethylene glycol over the (111) facet of six transition metals (Ni, Pd, Pt, Cu, Ag, Au) using an explicit solvation approach based on a hybrid quantum mechanical/molecular mechanical (QM/MM) description of the potential energy surface. A significant metal dependence on aqueous solvation effects was observed that suggests solvation effects must be studied in detail for every reaction system. The main reason for this dependence could be traced back to a different amount of charge-transfer between the adsorbed moieties and metals in the reactant and transition states for the different metal surfaces.

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

Trace LiBF 4 Enabling Robust LiF-Rich Interphase for Durable Low-Temperature Lithium-Ion Pouch Cells

Lithium-ion batteries (LIBs) with electrolytes containing lithium tetrafluoroborate (LiBF 4 ) can achieve large capacity retention under low temperature, but the accompanying rapid capacity decay inhibits commercialization. Here, in this study, the impact of LiBF 4 as a supplemental salt to LiPF 6 is systematically investigated using low ethylene carbonate (EC)-content electrolytes, along with a low-melting-point cosolvent. It is found that rational adjustment of the amount of LiBF 4 could not only regulate the interactions of anions and solvents in Li + solvation sheaths but also tune the composition and morphology of solid electrolyte interphase (SEI). It is worth noting that electrolytes with trace amount of LiBF 4 (0.05 M) show synergetic interaction between PF 6 - and Li + and decreased interaction between EC and Li + , achieving a dense and LiF-rich SEI, which enables a 200 mAh pouch cell with less gas generation, long-lived cycling, and higher low-temperature capacity, simultaneously. This work provides new insight into utilizing trace LiBF 4 for stable interface construction of durable low-temperature LIBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗