Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation structure”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗
DOE OSTI2024

Elucidating the Structure of the Eu‐EDTA Complex in Solution at Various Protonation States

Abstract Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu 3+ ion complexed with EDTA in all its protonation states in aqueous solution. Eu‐EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X‐ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu 3+ aqua ion and Eu‐EDTA complexes at pH 3 and 11. This work details how Eu‐EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu 3+ ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation.

Chemistry↗
DOE OSTI2020

Supercritical CO2 sterilization of N95 Masks

A preliminary investigation of the use of supercritical carbon dioxide for treating of 3M 1860 N95 masks was undertaken to evaluate a potential route to low-cost, scalable, sterilization of personal protective equipment for multiple reuse in hospital settings. Upon entering the supercritical regime, the normally distinct liquid and gaseous phases of CO 2 merge into a single homogeneous phase that has density, short-range order, and solvation capacity of a liquid, but the volume-filling and permeation properties that of a gas. This enables supercritical CO2 to function as a vehicle for delivery of biocidal agents such peracetic acid into microporous structures. The potentially adverse effect of a liquid-to-gas phase transition on mask filter media is avoided by conducting cleaning operations above 31 C, the critical temperature for carbon dioxide. A sample of fifteen 3M 1860 N95 masks was subjected to ten consecutive cycles of supercritical CO 2 cleaning to determine its effect on mask performance. These 15 masks, along with 5 control samples then underwent a battery of standardized tests at the CDC NIOSH NPPTL research facility in Pittsburgh, PA. The data from these tests strongly suggest (but do not prove) that supercritical carbon dioxide do not damage 3M 1860 N95 masks. Additional tests conducted during this project confirmed the compatibility of supercritical CO 2 with ventilator tubing that, like N95 masks, has been in short supply during portions of the COVID-19 pandemic and cannot be sterilized by conventional means. Finally, a control experiment was also conducted to examine the effect of supercritical CO 2 on a BSL-2 surrogate virus, vesicular stomatitis virus (VSV), Indiana serotype strain. In the absence of biocidal additives, supercritical CO 2 exhibited no measurable lethality against VSV. This surrogate virus experiment suggests that a biocidal additive such as peracetic acid will be necessary to achieve required sterilization metrics.

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

Critical role of solvation on CC13 porous organic cages for design of porous liquids

Efficient carbon capture requires the design of new materials with high CO 2 selectivity and gas adsorption capacity that can be incorporated into existing industrial processes. Porous liquids (PLs) are promising candidate materials that consist of a nanoporous host and a solvent forming a liquid with permanent porosity based on exclusion of the solvent from the interior of the nanoporous host. Stable PLs are based on solvent-nanoporous host interactions, which can be evaluated through molecular simulations. Here, time- and temperature-dependent density functional theory simulations were performed between four solvents, 2-bromophenol, 4-methylphenol, 2,4-dimethylphenol, and cyclohexanone and the CC13 porous organic cage (POC) as a prototypical PL composition. Overall, minimal reactions occurred in the PL including no changes in the POC structure. Additionally, POC-solvent coordination occurred through interactions of neighboring functional groups such as methyl/bromide and hydroxyl on the solvent molecules with the POC surface. Therefore, the location rather than the number of functional groups on the solvent molecule controls the POC-solvent interactions. Additionally, the POC pore window contracted or expanded up to 8% during solvation, which correlates with the experimental solubility and static solvent-POC binding, where solvents that caused less contraction of the POC pore window increased POC solubility. Finally, these results allow for the design of optimized POC-based PL compositions based on solvent-nanoporous host binding and variation in the pore window during solvation.

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

Changes in nitrate binding with lanthanides in BLPhen complexes

Preorganized ligands such as bis-lactam-1,10-phenanthroline (BLPhen) offer unique selectivity trends in solvent extraction of rare-earth ions, particularly from aqueous nitrate solutions. However, due to the ligand’s lipophilicity it is experimentally challenging to obtain crystal structures of the corresponding lanthanide (Ln):ligand complexes to describe the local structure of the first coordination shell and the role of the nitrate ions in the organic phase. In this work, we report first principles molecular dynamics (FPMD) simulations of the complexation of trivalent lanthanide ions (La, Nd, and Eu) with the BLPhen ligand in the presence of nitrate ions in the dichloroethane (DCE) solvent. We find that two nitrate anions are present in the first solvation shell and the third nitrate anion is far from the first solvation shell in all three [Ln(BLPhen) 2 ] 3+ complexes examined. Moving along the lanthanide series from La to Nd to Eu, the binding pocket formed by the two BLPhen ligands shrinks in size; as a result, the two nitrate configuration changes from one-monodentate/one-bidentate binding mode for La and Nd to all monodentate for Eu. Meanwhile, the total coordination number drops from 11 for La and Nd to 10 for Eu. More interestingly, the bidentate binding in [La(BLPhen) 2 ] 3+ and [Nd(BLPhen) 2 ] 3+ complexes is highly dynamic, frequently switching to monodentate and back. The FPMD insights into the differing modes and dynamics of nitrates in the first coordination shell will be useful for further atomistic understanding of the complex structure, formation, and stability in the organic phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2025

Insight Into Sulfur‐Containing Additive to Boost Anti‐Oxidation Ability of the Ether‐Based Electrolyte for Sodium‐Ion Full Batteries

Ether-based electrolytes are considered as promising candidates for sodium-ion batteries (SIBs) due to their high ionic conductivity and good compatibility with hard carbon (HC) anode. However, they suffer from poor anti-oxidative ability with unstable cathode electrolyte interface, inducing successive electrolyte decomposition and rapid capacity fading. Herein, the sulfur-containing additive (1,3-propanediol cyclic sulfate, PCS) is introduced to boost the stability of the electrode-electrolyte interface (EEI) in ether-based electrolyte. PCS largely participates in the inner Na + sheath, causing more diglyme (G2) molecules and fewer PF 6 - anions to occupy the inner Na + solvation sheath, avoiding the decomposition of G2 molecules at high operation voltage. Moreover, PCS is beneficial for simultaneously constructing thin and robust sulfur-containing EEI on both Prussian blue (PB) cathode and HC anode, maintaining the structural stability of PB and alleviating the dissolution of transition metals. These enable PB||HC pouch cells to deliver a capacity retention of 60.3% after 400 cycles, significantly higher than the 32.1% in PCS-free electrolytes. In conclusion, this work discloses the underlying mechanism of PCS to evoke the anti-oxidation ability of ether-based electrolyte and provides a promising method for realizing advanced sodium-ion full cells.

anti-oxidation ability↗
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 OSTI2026

Transfer Learning Meets Embedded Correlated Wavefunction Theory for Chemically Accurate Molecular Simulations: Application to Calcium Carbonate Ion Pairing

Achieving chemical accuracy for molecular simulations remains a central challenge in computational chemistry. Here, we present an embedded correlated wavefunction transfer learning (ECW-TL) framework for accurately simulating molecular dynamics in the condensed phase. ECW-TL incorporates high-level electron exchange and correlation effects in ECW theory while preserving the training and computational efficiency of machine-learned interatomic potentials. We demonstrate the framework on Ca 2+ –CO 3 2– ion pairing in aqueous solution, a key process underlying CO 2 mineralization in seawater. As proof of principle, we first show that fine-tuning a DFT-revPBE-D3(BJ) baseline model with embedded-DFT-SCAN data reproduces the DFT-SCAN free-energy surface within 1 kcal/mol across all solvation states. Extending the framework to embedded MP2 and localized natural-orbital CCSD(T) further refines the free-energy profile, revealing the crucial role of exact electron exchange and correlation in determining ion-pair stability and structure. The computed ion-pair association free energy is in quantitative agreement with experimental measurements, further validating the accuracy of the ECW-TL framework. ECW-TL thus provides a general, data-efficient route for transferring CW accuracy to efficient simulations of complex aqueous and interfacial chemical processes.

cluster chemistry↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗
DOE OSTI2026

Cyclodextrin-Derived Porous Liquids Enabled by In Situ Solvation Shell Formation

Porous liquids (PLs) represent a unique platform for molecular separations by combining permanent porosity with liquid-phase mobility. However, it remains a formidable challenge to construct and stabilize PLs with sub-5 Å pores using readily available porous host and liquid media. Here, we report the construction of cyclodextrin (CD)-derived PLs enabled by in situ solvation shell formation. The acid–base neutralization reaction between CD and an organic base was leveraged to generate a thin ionic solvation shell around the CD host, effectively liquefying CD and preventing its segregation in the liquid base medium while preserving accessible molecular-scale cavities. Spectroscopic analysis, neutron scattering, density functional theory calculations, and molecular dynamics simulations collectively confirm the structural evolution and existence of abundant internal porosity in PLs. The unique architectures of CD-derived PLs enable highly selective encapsulation of fluorinated alkanes and significantly enhanced uptake of inert gases. This facile and generalizable strategy enables construction of high-quality PLs with engineered ultramicroporosity to facilitate molecular separations.

Cavities↗
DOE OSTI2023

Non-polar ether-based electrolyte solutions for stable high-voltage non-aqueous lithium metal batteries

The electrochemical instability of ether-based electrolyte solutions hinders their practical applications in high-voltage Li metal batteries. To circumvent this issue, here, we propose a dilution strategy to lose the Li + /solvent interaction and use the dilute non-aqueous electrolyte solution in high-voltage lithium metal batteries. We demonstrate that in a non-polar dipropyl ether (DPE)-based electrolyte solution with lithium bis(fluorosulfonyl) imide salt, the decomposition order of solvated species can be adjusted to promote the Li + /salt-derived anion clusters decomposition over free ether solvent molecules. This selective mechanism favors the formation of a robust cathode electrolyte interphase (CEI) and a solvent-deficient electric double-layer structure at the positive electrode interface. When the DPE-based electrolyte is tested in combination with a Li metal negative electrode (50 μm thick) and a LiNi 0.8 Co 0.1 Mn 0.1 O 2 -based positive electrode (3.3 mAh/cm 2 ) in pouch cell configuration at 25 °C, a specific discharge capacity retention of about 74% after 150 cycles (0.33 and 1 mA/cm 2 charge and discharge, respectively) is obtained.

25 ENERGY STORAGE↗
DOE OSTI2025

Computationally Guided and Experimentally Validated Design of Custom Chelators for Critical Mineral Recovery

Selective, high throughput separation of target critical metals from complex environments such as fly ash leachates and mining process streams presents a significant challenge for economical production. Custom chelators and sorbents are an attractive technology for selective metal extraction, however it can be difficult to predict their performance, and significant experimental efforts are often required to develop chelating technologies. Here, we present a computational strategy focused on modelling chelator-metal binding interactions and benchmark these results versus experimental data. A computational pipeline combining forcefield, semiempirical, and meta-GGA methods with a thermodynamic framework optimized for error cancellation has been developed to predict binding energies of chelator complexes towards critical mineral recovery applications. This approach, originally validated on [2.2.2] cryptates binding mono- and divalent cations, demonstrated robust predictive capabilities with an R2 of 0.850 against experimental aqueous binding energies. The workflow includes metadynamics for exploring high-dimensional potential energy surfaces and a cluster-continuum model for accurate yet computationally efficient solvation modeling. Error cancellation between solvation energies of free and chelator-coordinated ions enables faster convergence, even with finite cluster sizes. Initial studies on the cryptates revealed consistent metal-ligand coordination patterns, with systematic variations influenced by ion size and charge, highlighting key structural features linked to binding selectivity. Further studies of a proprietary chelator have resulted in identification of previously unreported selectivity towards economically significant metals, which in-house experiments have confirmed, demonstrating the feasibility of this approach. By applying this methodology to new chelators targeting critical minerals such as lithium, cobalt, nickel and other strategic metals, we aim to accelerate the discovery of next-generation chelators for efficient recovery, recycling, and separation processes. This computational framework serves as the backbone of a high-throughput design pipeline tailored for sustainable resource utilization and may be applied to a wide range of systems to meet experimental needs.

computational materials↗
DOE OSTI2024

Theoretical Investigation of the Adsorbate and Potential–Induced Stability of Cu Facets During Electrochemical CO 2 and CO Reduction

The activity and product selectivity of electrocatalysts for reactions like the carbon dioxide reduction reaction (CO 2 RR) are intimately dependent on the catalyst's structure and composition. While engineering catalytic surfaces can improve performance, discovering the key sets of rational design principles remains challenging due to limitations in modeling catalyst stability under operating conditions. Herein, we perform first-principles density functional calculations adopting implicit solvation methods with potential control to study the influence of adsorbates and applied potential on the stability of different facets of model Cu electrocatalysts. Using coverage dependencies extracted from microkinetic models, we describe an approach for calculating potential and adsorbate-dependent contributions to surface energies under reaction conditions, where Wulff constructions are used to understand the morphological evolution of Cu electrocatalysts under CO 2 RR conditions. Here we identify that CO*, a key reaction intermediate, exhibits higher kinetically and thermodynamically accessible coverages on (100) relative to (111) facets, which can translate into an increased relative stabilization of the (100) facet during CO 2 RR. Our results support the known tendency for increased (111) faceting of Cu nanoparticles under more reducing conditions and that the relative increase in (100) faceting observed under CO 2 RR conditions is likely attributed to differences in CO* coverage between these facets.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

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

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗