Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anion dynamics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗
DOE OSTI2026

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗
DOE OSTI2026

A Solid State Zwitterionic Plastic Crystal With High Static Dielectric Constant

Solid materials with a high dielectric constant have a wide range of applications in the energy storage field. In this research, an imidazolium-based zwitterion is designed, synthesized, and confirmed to have a plastic crystal phase based on the following experimental and computational evidence: (i) the presence of long-range order with weak intermolecular forces and competing attractive-repulsive interactions along different crystallographic directions; (ii) the observation of more than one endotherm on heating including a solid-solid phase transition at T s-s = −26 °C and melting of the plastic crystal at T m = 72°C; (iii) a low entropy of fusion at melting (2.1 JK −1 mol −1 ); (iv) a strongly anisotropic morphology; (v) relatively fast dynamics originating from short-range degrees of freedom. Furthermore, it exhibits a very high dielectric constant in the plastic crystal solid state (147 at −10°C and 103 at 70°C) due to the rotational degrees of freedom of plastic crystals that arise from weak net intermolecular interactions of zwitterions due to only having two carbons between the anion and cation. This material conveniently remains in the plastic crystal phase within 50 K of ambient temperature. This discovery opens new opportunities in the search for solid-state high dielectric constant materials.

dielectric relaxation↗
DOE OSTI2020

Onsager Transport Coefficients and Transference Numbers in Polyelectrolyte Solutions and Polymerized Ionic Liquids

Electrolytes featuring negatively charged polymers, such as nonaqueous polyelectrolyte solutions and polymerized ionic liquids, are currently under investigation as potential high cation transference number (t+) electrolytes for lithium-ion batteries. Herein, we use coarse-grained molecular dynamics simulations to characterize the Onsager transport coefficients of polyelectrolyte solutions as a function of chain length and concentration. For all the systems studied, we find that the rigorously computed transference number is substantially lower than that approximated by the ideal solution (Nernst-Einstein) equations typically used to characterize these systems due to the presence of strong anion-anion and cation-anion correlations. None of the polyelectrolyte solutions achieve t+ greater than that of the conventional binary salt electrolyte, with some solutions having negative t+. Overall, this work demonstrates that the Nernst-Einstein assumption does not provide a physically meaningful estimate of the transference number in these solutions and calls into question the expectation of polyelectrolytes to exhibit a high cation transference number.

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

Characterization of Contaminants and Mobility in WMA A-AX Boreholes D0006 and D0008

This study was initiated to determine if subsurface contaminants in the A-AX vadose zone at the Hanford site (boreholes D0006 and D0008) caused corrosion of nearby stainless steel well casings (299-E-24-19, 299-E25-46, and 299-E25-236) and, based on the type and distribution of contaminants, determine if contaminants are from one or more sources including a) the A-104 and/or A-105 tanks, b) the 242A evaporator, and c) the 200-E-286 ditch. A tiered approach was used for characterizing contaminants in vadose zone sediments starting with mobile and total contaminant concentrations (Tier 1), characterizing contaminant speciation (Tier 2), and measuring contaminant leaching (Tier 3). Total and mobile contaminant measurements at seven depths showed elevated sulfate and nitrate in all samples (130 to 280 ft depth), elevated chloride in five depths, and no elevated radionuclides except I-129 in D0006 at 206 to 208 ft depth. If these contaminants were initially acidic, stainless-steel corrosion would occur. However, because sediments likely have a high acid neutralization capacity, an acidic spill would have been neutralized within tens of feet, so casing corrosion would likely occur in shallow sediments. GEL labs data also shows low (but above natural) chromate from shallow to 170 ft depth, then higher chromate to 280 ft, suggesting a surface Cr spill or casing corrosion. GEL labs data also showed evidence of more than one plume, with deep nitrate, sulfate, and chloride migration, but only shallow tritium migration. Tritium should migrate nearly unretarded in the subsurface like nitrate and chloride. Elevated aqueous cations and anions in the D0006 206 to 208' depth pore water (SO 4 2- , NO 3 - , Na + , K + , Ca 2+ , Mg 2+ , Si) could be the result of acid disposal and calcite dissolution at shallow depth, with deeper migration of SO 4 2- , NO 3 - from acids and Ca 2+ and Mg 2+ from calcite. Elevated Si (160 mg/L) indicates dissolution of clays or other silicates, which can also occur in acids. Acid extractable metals did show elevated Cr (but not Fe, Ni, Mo, and Mn) above natural levels, which may indicate stainless-steel well corrosion or transport of metals from a shallow spill. The elevated I-129 in D0006 at 206 to 208' depth was low with 0.34 μg/g in pore water and 1.5 μg/g total extractable I-129 in the sediment, which indicates significant fraction of I-129 was bound in one or more precipitates. Although total I-129 was below detection limits in most extractions due to high molybdate interference, I-127 analysis provided some insight into potential precipitate phases. Total I-127 was mainly extracted by acetate or acetic acid (67%), which may be iodine in calcite. Inorganic carbon measurement identified 1.65% calcite, and total organic carbon (TOC) was below detection limits, so iodine was not associated with organic matter. Nearly all the I-127 was present as iodide. In 1-D leach experiments, I-127 as iodide leached from the sediment within 2 pore volumes. Due to high and dynamic molybdate leaching from the sediment and interference with I-129 analysis, there were few quantifiable I-129 measurements. The calculated molybdate release rates from the sediment suggested Mo was being released from the same precipitate phase, such as molybdate incorporated into calcite. Overall, because there was little radioactive contamination measured in D0006 and D0008 boreholes, potential releases from A-104 and A-105 did not reach this location south of the A tank farm. Elevated cations and anions found in D0006 and D0008 at depth as the result of acidic disposal could cause shallow stainless steel casing corrosion. The nearby 242A evaporator generally contained alkaline waste and tritium, but the nearby 200-E-286 ditch generally contained high Cl - liquid effluent, so these are unlikely to be acidic sources.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Cations impact radical reaction dynamics in concentrated multicomponent aqueous solutions

Ultraviolet (UV) photolysis of nitrite ions (NO 2 – ) in aqueous solutions produces a suite of radicals, viz., NO · , O – , · OH, and · NO 2 . The O – and NO · radicals are initially formed from the dissociation of photoexcited NO 2 – . The O – radical undergoes reversible proton transfer with water to generate · OH. Both · OH and O – oxidize the NO 2 – to · NO 2 radicals. The reactions of · OH occur at solution diffusion limits, which are influenced by the nature of the dissolved cations and anions. Here, we systematically varied the alkali metal cation, spanning the range from strongly to weakly hydrating ions, and measured the production of NO · , · OH, and · NO 2 radicals during UV photolysis of alkaline nitrite solutions using electron paramagnetic resonance spectroscopy with nitromethane spin trapping. Comparing the data for the different alkali cations revealed that the nature of the cation had a significant effect on production of all three radical species. Radical production was inhibited in solutions with high charge density cations, e.g., lithium, and promoted in solutions containing low charge density cations, e.g., cesium. Through complementary investigations with multinuclear single pulse direct excitation nuclear magnetic resonance (NMR) spectroscopy and pulsed field gradient NMR diffusometry, cation-controlled solution structures and extent of NO 2 – solvation were determined to alter the initial yields of · NO and · OH radicals as well as alter the reactivity of NO 2 – toward · OH, impacting the production of · NO 2 . The implications of these results for the retrieval and processing of low-water, highly alkaline solutions that comprise legacy radioactive waste are discussed.

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

NMR Relaxometry and Diffusometry Analysis of Dynamics in Ionic Liquids and Ionogels for Use in Lithium-Ion Batteries

Solvation and transport properties of methlyviologen dichloride in mixtures of choline chloride (ChCl) and ethylene glycol (EG) including the deep eutectic solvent (DES) ethaline were studied through physical property measurements, cyclic voltammetry, and Raman spectroscopy. The transport properties of the solvate structures in ethaline and ChCl:EG mixtures were compared to a well-studied ionic liquid (IL), 1-methyl-1-propylpyrrolidinium bis(trifluoromethylsulfonyl)imide, [PYR13][TFSI] through the lens of hole theory. Exchanging the chloride anion of ethaline with bis(trifuoromethylsufonyl)imide) anion, [TFSI], is found to increase the fluidity of the solvent and promoted stronger solute-solvent interactions. Raman spectroscopy shows that methylviologen is solvated by EG in solutions with ChTFSI but interstitially accommodated in void spaces in both ChCl solutions and [PYR13][TFSI]. Complex solvents such as DESs are regarded as “designer solvents” and we demonstrate that the physical properties and solvation characteristics of these fluids can be varied by anion exchange of hydrogen bond acceptors.

25 ENERGY STORAGE↗
DOE OSTI2022

Ion Pairing and Molecular Orientation at Liquid/Liquid Interfaces: Self-Assembly and Function

We report that molecular orientation plays a pivotal role in defining the functionality and chemistry of interfaces, yet accurate measurements probing this important feature are few, due, in part, to technical and analytical limitations in extracting information from molecular monolayers. For example, buried liquid/liquid interfaces, where a complex and poorly understood balance of inter- and intramolecular interactions impart structural constraints that facilitate the formation of supramolecular assemblies capable of new functions, are difficult to probe experimentally. Here, we use vibrational sum-frequency generation spectroscopy, numerical polarization analysis, and atomistic molecular dynamics simulations to probe molecular orientations at buried oil/aqueous interfaces decorated with amphiphilic oligomers. We show that the orientation of self-assembled oligomers changes upon the addition of salts in the aqueous phase. The evolution of these structures can be described by competitive ion effects in the aqueous phase altering the orientations of the tails extending into the oil phase. These specific anionic effects occur via interfacial ion pairing and associated changes in interfacial solvation and hydrogen-bonding networks. These findings provide more quantitative insight into orientational changes encountered during self-assembly and pave the way for the design of functional interfaces for chemical separations, neuromorphic computing applications, and related biomimetic systems.

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

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.

Coarse-grained simulations↗
DOE OSTI2023

Polymer electrolytes based on protic ionic liquids with perfluorinated anions for safe lithium-ion batteries

Here, the quest for safe and high-performance polymer electrolytes in lithium-ion batteries (LIBs) has led researchers to explore protic ionic liquids (PILs) as potential candidates to be entrapped in polymer matrices. In this context, we present an investigation into solid polymeric systems based on poly(methyl methacrylate) (PMMA) as a host for PILs, featuring 1,8-diazabicyclo-[5,4,0]-undec-7-ene (DBU) cation paired with three different anions: bis(trifluoromethanesulfonyl)imide (TFSI – ), trifluoromethanesulfonate (TFO – ), and (trifluoromethanesulfonyl-nonafluorobutylsulfonyl)imide (IM14 – ). Additionally, we explore the lithium-doped IM14-gel-like system to broaden our understanding of these intriguing materials. Through comprehensive thermal analysis, solid-state NMR, and diffusion NMR techniques, we delve into the interactions and structural features of these binary and ternary polymeric systems. Our investigation reveals unique dynamics and ion interactions within the PMMA matrix, shedding light on the potential of these materials for advanced energy storage technologies. Particularly, we highlight the distinctive features of DBUH-IM14 and its specific interaction with the polymeric matrix and the lithium ions, underscoring its significance in advancing safer and more efficient energy storage devices.

25 ENERGY STORAGE↗
DOE OSTI2020

Substituent Effects on Exchange Coupling and Magnetic Relaxation in 2,2'-Bipyrimidine Radical-Bridged Dilanthanide Complexes

Systematic analysis of related compounds is crucial to the design of single-molecule magnets with improved properties, yet such studies on multinuclear lanthanide complexes with strong magnetic coupling remain rare. In this work, we present the synthesis and magnetic characterization of the series of radical-bridged dilanthanide complex salts [(Cp* 2 Ln) 2 (μ-5,5'-R 2 bpym)](BPh 4 ) (Ln = Gd, Dy; R = NMe 2 (1), OEt (2), Me (3), F (4); bpym = 2,2'-bipyrimidine). Modification of the substituent on the bridging 5,5'-R 2 bpym radical anion allows the magnetic exchange coupling constant, J Gd-rad , for the gadolinium compounds in this series to be tuned over a range from -2.7 cm -1 (1) to -11.1 cm -1 (4), with electron-withdrawing or -donating substituents increasing or decreasing the strength of exchange coupling, respectively. Modulation of the exchange coupling interaction has a significant impact on the magnetic relaxation dynamics of the single-molecule magnets 1-Dy through 4-Dy, where stronger J Gd-rad for the corresponding Gd 3+ compounds is associated with larger thermal barriers to magnetic relaxation ( U eff ), open magnetic hysteresis at higher temperatures, and slower magnetic relaxation rates for through-barrier processes. Further, we derive an empirical linear correlation between the experimental U eff values for 1-Dy through 4-Dy and the magnitude of J Gd-rad for the corresponding gadolinium derivatives that provides insight into the electronic structure of these complexes. This simple model applies to other organic radical-bridged dysprosium complexes in the literature, and it establishes clear design criteria for increasing magnetic operating temperatures in radical-bridged molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Surfactant-Aided Low-Salinity Waterflooding for Low-Temperature Carbonate Reservoirs

Carbonate reservoirs tend to be oil-wet/mixed-wet and heterogeneous because of mineralogy and diagenesis. The objective of this study is to improve oil recovery in low-temperature dolomite reservoirs using low-salinity and surfactant-aided spontaneous imbibition. The low-salinity brine composition was optimized using ζ-potential measurements, contact-angle (CA) experiments, and a novel wettability-alteration measure. Significant wettability alteration was observed on dolomite rocks at a salinity of 2,500 ppm. We evaluated 37 surfactants by performing CA, interfacial-tension (IFT), and spontaneous-imbibition experiments. Here, three (quaternary ammonium) cationic and one (sulfonate) anionic surfactants showed significant wettability alteration and produced 43–63% of original oil in place (OOIP) by spontaneous imbibition. At a low temperature (35°C), oil recovery by low-salinity effect is small compared with that by wettability-altering surfactants. Coreflood tests were performed with a selected low-salinity cationic surfactant solution. A novel coreflood was proposed that modeled heterogeneity and dynamic imbibition into low-permeability regions. The results of the “heterogeneous” coreflood were consistent with that of spontaneous-imbibition tests. These experiments demonstrated that a combination of low-salinity brine and surfactants can make originally oil-wet dolomite rocks more water-wet and improve oil recovery from regions bypassed by waterflood at a low temperature of 35°C.

02 PETROLEUM↗
DOE OSTI2021

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

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

Molecular Insights Into the Ionic Assembly of Poly-Galacturonic Acid Oligomers - Impact of Charge, Ionic Radius, and Polymer Functionalization

Pectin, a major class of matrix polysaccharides present in plant cell walls (PCW), contains widespread anionic saccharides that cross-link in the presence of cations. It modulates important functions such as cell-cell adhesion and determines the PCW's biomechanical properties. It is known that mono-, di-, and tri-valent cations facilitate cross-linking; however, significant knowledge gaps remain in understanding the structure and mechanism of pectin cross-linking. In this study, replica-exchange molecular dynamics (REMD) simulations were employed to elucidate the role of ionic charge, ionic radii, and functional groups on the cross-linking of homogalacturonan (HG), the most abundant pectin molecule. Our enhanced sampling approach in fully solvated environments suggests more effective cross-linking with higher-valent and smaller ions, and that the "zipper" conformation is more favorable than the prevalent "egg-box" conformation. These findings advance our fundamental understanding of pectin matrix structure in PCWs and provide a solid foundation to probe structure-property relationships in pectic polysaccharides.

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

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

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

Hidden features in the OH-stretching spectra of amino acid decorated air–water interfaces

Chemical reactivity at the air–water interface is governed by the interfacial solvation of reactive species. For instance, during aqueous amino acid-based CO 2 absorption, water reorganizes around the reactive sites and couples dynamically with reaction pathways, facilitating the reaction. In this context, surface-sensitive vibrational sum-frequency generation (vSFG) spectroscopy can probe the OH stretch vibrations of interfacial water and determine the solvation structures around reactants and products, thereby furthering our understanding of the role of interfacial solvation. However, vSFG spectra of the air–water interface in the presence of charged species can be remarkably complex; key species-bound local water structures with distinct orientations may be hidden beneath prominent vSFG peaks arising from water–water hydrogen bonds and remain difficult to resolve. Here, we measure and compute vSFG spectra of the water OH stretch at air–water interfaces decorated with amino acids in their zwitterionic and anionic forms, as well as equimolar mixtures of these forms with bicarbonate. The latter represents post-CO 2 -absorption conditions. We find that computing depth- and frequency-dependent spectral densities—decomposed into contributions from water molecules hydrogen-bonded exclusively to other water molecules, exclusively to amines, exclusively to carboxylates, or shared between these polar/charged groups—is indispensable for accurate interpretation of the vSFG spectra. Key findings include orientational flip-flop in water sub-layers, strong carboxylate-water H-bonding, and water orientational ordering extending into the bulk aqueous phase induced by anionic amino acids. Here, this study provides a computational spectroscopic platform for improved understanding of interfacial solvation relevant to interfacial reactivity.

Air-water interface↗
DOE OSTI2021

Understanding Gas Absorption in Multivalent Ionic Liquids via Solute-Solvent Interaction Analyses

Solute-solvent interactions are often evaluated to understand and explain gas solubilities. This information is regularly used to interpret gas solubility and separation within ionic liquids (ILs). Here, we use several different quantum chemical and previous molecular dynamics simulation results to show that strong solute-solvent interactions neither lead to high solubility nor to high selectivity, with respect to other solutes. Instead, the solute-solvent interaction strength with different isolated IL ions is affected by the electrostatic extremes. Also, the solute-solvent interaction with different IL molecules is strongly correlated to the molecular volumes or polarizabilities (within a series of same-anion based ILs).

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

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗