SEARCH · Engineering Papers

Results for “Solvation models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Solvation Environment and Interface Dynamics of Li 2 B 12 H 12 and Li 2 B 12 F 12 Electrolytes Uncovered by Theory and Operando Optical and FTIR Spectroelectrochemistry

In this work, we evaluated two closo-borate salts (Li 2 B 12 H 12 and Li 2 B 12 F 12 ) in propylene carbonate from theoretical and experimental perspectives to understand how the coordination environment influences their spectroscopic and electrochemical properties. The coordination environments of the closo-borate salts were modeled via density functional theory (DFT) and molecular dynamics (MD). Vibrational spectra calculated from the predicted coordination environments are in agreement with experimentally measured steady-state FTIR data. This theoretical investigation also suggested that Li 2 B 12 F 12 would possess a higher ionic conductivity than Li 2 B 12 H 12 , which was corroborated experimentally. Additionally, an electrochemical cell was designed and fabricated that enabled operando optical and FTIR spectroelectrochemical (OP-IR-SEC) measurements. This allowed for the simultaneous measurement of the relative changes of species at a lithium electrode-liquid electrolyte interface and the visualization of lithium plating at the electrode surface. This technique could provide new chemical insights and potentially link optical changes at the electrode–electrolyte interface to specific chemical species in similar electrochemical systems. The Li 2 B 12 F 12 electrolyte was found to have a higher thermal stability, which may find utility in applications for batteries that are subject to high-temperature conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗

DOE OSTI2024

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction

Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗

DOE OSTI2025

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗

DOE OSTI2026

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗

DOE OSTI2025

Self‐Assembly Methods Induce Different Individual‐Polymer‐Block Solvation Responses

An understanding of block-specific responses to stimuli in self-assembled copolymers is essential for the effective design of responsive materials. Here, materials formed from different post-reaction processing methods of the same ROMP block copolymer exhibit different block-specific solvent responses, as characterized in situ by fluorescence lifetime imaging microscopy (FLIM). FLIM data were acquired by tagging the polar or nonpolar block separately with a covalently incorporated viscosity-sensitive fluorescent molecular rotor to measure the changing tightness or looseness of assembly upon changes in solvent composition. Rapid precipitation from an insoluble solvent to form a film results in tighter assembly and less difference in tightness/looseness between the two blocks. Further, it interrupts individual block-solvent response behavior. Both results indicate an apparent disruption of the core–shell assembly. This model from FLIM is further supported by differential scanning calorimetry (DSC) and small-angle X-ray scattering (SAXS) data from isolated solids, which show tighter long-range interactions and more long-range order in these kinetically trapped states. Here, the experiments develop FLIM as a method for pinpointing stimuli-responsive behaviors from different processing methods to individual blocks. Together these outcomes provide a future handle for tailoring solvent-triggered assembly/disassembly behavior.

López, Pía A. [University of California, Irvine, C↗

DOE OSTI2024

Modeling prebiotic chemistries with quantum accuracy at classical costs

Molecular Dynamics (MD) simulations using classical force-fields are commonly employed in numerous scientific investigations. However, many natural processes involve bond breaking and quantum forces. This complexity is compounded by the presence of multiple competing length and timescales. For example, accurately modeling the thermodynamics and dynamics of a chemical reaction requires accounting for the concerted movements of numerous solvent molecules and ions with their own fast or slow timescales. While widely used static Density Functional Theory (DFT) calculations at 0 temperature can be beneficial for such investigations, they do not account for dynamics, and lack precision in describing the molecular environments. They particularly fail at correct, rigorous treatments of finite-temperature fluctuations, and thus generalization to experimentally relevant conditions. In PNAS Benayad et al develop a scalable, generalizable approach for designing Neural Network Potentials (NNPs) that can handle chemical reactivity in solvated systems with quantum accuracy at classical costs. Specifically, they study phosphoester bond formation and rupture, which is fundamentally relevant to the Phosphorus-Oxygen bond formation central to life, and especially for the RNA world hypothesis. The framework developed here has the potential to generalize to different chemical reactions of energy and biological relevance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Solvation Structure of 237 Np 4+ in a Noncomplexing Environment

Here, the solvation structure of an Np 4+ ion in an aqueous, noncomplexing and nonoxidizing environment of trifluoromethanesulfonic (triflic) acid was investigated with X-ray absorption spectroscopy (XAS) combined with ab initio molecular dynamics (AIMD) and time-dependent density functional theory (TDDFT) calculations. Np L III -edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data were collected for Np 4+ in 1, 3, and 7 M triflic acid using a laboratory-scale spectrometer and separately at a synchrotron facility, producing data sets in excellent agreement. TDDFT calculations revealed a weak pre-edge feature not previously reported for Np L III -edge XANES. AIMD modeling results showed differences in the hydration shell of the Np 4+ ion at different concentrations of triflic acid; these results are supported by the experiment. EXAFS fit models to the experiment resulted in similar coordination of Np 4+ in noncomplexing aqueous media as reported in the literature for 1 M perchloric acid but, together with calculations, revealed more than one distance between Np and O atoms in 7 M triflic acid. These results imply monodentate coordination with sulfonate groups in 7 M triflic acid and suggest the possibility of proto-neptunyl species in relatively low-concentration Np 4+ acid solutions.

Boglaienko, Daria V. [Pacific Northwest National L↗

DOE OSTI2025

The Influence of Charge Correlation and Ion Solvation on the Phase Behavior of Single-Ion Conducting Polymer Blend Electrolytes Using SAXS/SANS

Single-ion conducting polymer blends (SICPBs) have demonstrated exceptional electrochemical performance as solid-state battery electrolytes; however, their nanoscale morphology and thermodynamic behavior remain unexplored. In this work, we investigate blends composed of deuterated poly(ethylene oxide) and poly[lithium sulfonyl(trifluoromethane sulfonyl)imide methacrylate], dPEO/P(LiMTFSI), and report the first experimental study of the nanostructures of charge-neutral polymer blends using small-angle neutron scattering (SANS) and small-angle X-ray scattering (SAXS). Despite the macroscopic miscibility indicated by a single glass-transition temperature, SANS and SAXS results reveal disordered, charge-correlated nanostructures that are strongly influenced by blend composition and temperature. At low concentrations of charge polymer, the scattering is dominated by concentration fluctuations, and the random phase approximation is applied to extract values of the Flory–Huggins interaction parameter, χ SC . At higher charged polymer content, concentration fluctuations are suppressed, and a correlation model is used to characterize the nanostructures of the charge correlations. We find that the structures of the charge correlations are highly dependent on blend composition─consistent with predictions from Sing’s self-consistent field theory-liquid state models. Understanding these features is essential for uncovering the ion transport mechanism that leads to improved electrochemical performance previously reported in SICPB systems.

25 ENERGY STORAGE↗

DOE OSTI2026

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data provided valuable insights into another type of non-DLVO interaction, namely the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. In conclusion, our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1997

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗

NASA NTRS1996

Mechanism of unassisted ion transport across membrane bilayers

To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.

NASA Discipline Exobiology↗

NASA NTRS1990

Cavities in molecular liquids and the theory of hydrophobic solubilities

Thermal configurational data on neat liquids are used to obtain the work of formation of hard spherical cavities of atomic size in six molecular solvents: n-hexane, n-dodecane, n-undecyl alcohol, chloroform, carbon tetrachloride, and water. These results are used to test a recent suggestion that the differences between nonaqueous solvents and liquid water in solvation of inert gases are not principally due to the hydrogen-bonded structure of liquid water but rather to the comparatively small size of the water molecule. The frequencies of occurrence of cavities in liquid water can be meaningfully distinguished from those in the organic solvents. Liquid water has a larger fractional free volume, but that free volume is distributed in smaller packets. With respect to cavity work, water is compared to a solvent of the same molecular density and composed of hard spheres of the same size as the water molecule. That comparison indicates that the hard-sphere liquid finds more ways to configure its free volume in order to accommodate an atomic solute of substantial size and thus, would be more favorable solvent for inert gases. The scaled particle model of inert gas solubility in liquid water predicts cavity works 20% below the numerical data for TIP4P water at 300 K and 1.0 g/cm3 for cavity radii near 2.0 angstroms. It is argued that the sign of this difference is just the sign that ought to be expected and that the magnitude of this difference measures structural differences between water and the directly comparable hard-sphere liquid. In conjunction with previous data, these results indicate that atomic sized cavities should be considered submacroscopic.

NASA Discipline Exobiology↗

DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗

DOE OSTI2025

Revealing the roles of the solid–electrolyte interphase in designing stable, fast-charging, low-temperature Li-ion batteries

Designing the solid–electrolyte interphase (SEI) is critical for stable, fast-charging, low-temperature Li-ion batteries. Fostering a “fluorinated interphase,” SEI enriched with LiF, has become a popular design strategy. Although LiF possesses low Li-ion conductivity, many studies have reported favorable battery performance with fluorinated SEIs. Such a contradiction suggests that optimizing SEI must extend beyond chemical composition design to consider spatial distributions of different chemical species. In this work, we demonstrate that the impact of a fluorinated SEI on battery performance should be evaluated on a case-by-case basis. Sufficiently passivating the anode surface without impeding Li-ion transport is key. We reveal that a fluorinated SEI containing excessive and dense LiF severely impedes Li-ion transport. In contrast, a fluorinated SEI with well-dispersed LiF (i.e., small LiF aggregates well mixed with other SEI components) is advantageous, presumably due to the enhanced Li-ion transport across heterointerfaces between LiF and other SEI components. An electrolyte, 1 M LiPF 6 in 2-methyl tetrahydrofuran (2MeTHF), yields a fluorinated SEI with dispersed LiF. This electrolyte allows anodes of graphite, μSi/graphite composite, and pure Si to all deliver a stable Coulombic efficiency of 99.9% and excellent rate capability at low temperatures. Pouch cells containing layered cathodes also demonstrate impressive cycling stability over 1,000 cycles and exceptional rate capability down to −20 °C. Through experiments and theoretical modeling, we have identified a balanced SEI-based approach that achieves stable, fast-charging, low-temperature Li-ion batteries.

25 ENERGY STORAGE↗