Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MOLECULAR ION”

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

Nanoconfined Grain Boundaries Increase the Conductivity of Polycrystalline Molecular Crystals

Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.

Crystals↗
DOE OSTI2020

Semiconducting Small Molecules as Active Materials for p-Type Accumulation Mode Organic Electrochemical Transistors

A series of semiconducting small molecules with bithiophene or bis-3,4-ethylenedioxythiophene cores are designed and synthesized. The molecules display stable reversible oxidation in solution and can be reversibly oxidized in the solid state with aqueous electrolyte when functionalized with polar triethylene glycol side chains. Evidence of promising ion injection properties observed with cyclic voltammetry is complemented by strong electrochromism probed by spectroelectrochemistry. Blending these molecules with high molecular weight polyethylene oxide (PEO) is found to improve both ion injection and thin film stability. The molecules and their corresponding PEO blends are investigated as active layers in organic electrochemical transistors (OECTs). For the most promising molecule:polymer blend (P4E4:PEO), p-type accumulation mode OECTs with µA drain currents, μS peak transconductances, and a µC* figure-of-merit value of 0.81 F V -1 cm -1 s -1 are obtained.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular origin of negative lithium transference in electrolytes with star-shaped multivalent anions

Large-scale molecular dynamics simulations illustrate that highly correlated cation–anion motion leads to negative t 0+ on the order of 1 in lithium electrolytes with star-shaped multivalent anions. Large multivalent anions have gained increasing attention for their potential to improve lithium transference in electrolytes. We employ large-scale molecular dynamics simulations based on the Onsager transport framework to investigate ion transport in a lithium electrolyte with star-shaped multivalent anions. The simulations show that t 0+, the cation transference number with respect to solvent velocity, is negative over a wide range of concentration. This is consistent with experimental data reported previously. The simulation-based Onsager transport coefficients reveal that the magnitudes of the cation–cation, anion–anion, and cation–anion correlations are comparable, a signature of highly correlated motion in the electrolyte. Examination of the cation solvation environment indicates the presence of strong cation–anion association across the entire concentration range, which leads to negative t 0+ on the order of −1. Both simulation and experiment also show that the maximum value of t 0+ reaches 0 when the cation concentration is c + = 0.4 M. This is the concentration at which the anions begin to spatially overlap, and lithium ions serve as dynamic linkers to balance cation–cation and cation–anion correlations. Our results provide molecular-level insights into the origin of transference in multivalent electrolytes.

Fang, Chao↗
DOE OSTI2026

Molecular insights into CO 2 -to-bicarbonate transformation in functionalized anion exchange ionomers for electrochemical separations

Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.

Bipolar membranes, Reactive transformation of CO2,↗
DOE OSTI2023

Quantum simulation of conical intersections using trapped ions

Conical intersections often control the reaction products of photochemical processes and occur when two electronic potential energy surfaces intersect. Theory predicts that the conical intersection will result in a geometric phase for a wavepacket on the ground potential energy surface, and although conical intersections have been observed experimentally, the geometric phase has not been directly observed in a molecular system. Herein we use a trapped atomic ion system to perform a quantum simulation of a conical intersection. The ion’s internal state serves as the electronic state, and the motion of the atomic nuclei is encoded into the motion of the ions. The simulated electronic potential is constructed by applying state-dependent optical forces to the ion. We experimentally observe a clear manifestation of the geometric phase using adiabatic state preparation followed by motional state measurement. Our experiment shows the advantage of combining spin and motion degrees for quantum simulation of chemical reactions.

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

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

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

Joll, Kit↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

Machine learning predictions of diffusion in bulk and confined ionic liquids using simple descriptors

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350–500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. Here, we also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Salts in Hot Water: Developing a Scientific Basis for Supercritical Desalination and Strategic Metal Recovery

Given the steady growth in world population and the ever-increasing fresh water demand, the only sustainable water supply option is desalination. Conceptually, desalination is very simple - just separate dissolved salts from water. While desalination research programs have existed since the 1960s and have resulted in a multitude of approaches, realizing affordable desalination has proven to be a challenge. Our goal is to meet this challenge using low cost heat sources to drive a liquid-discharge-free desalination process that can be co-mingled with the extraction of economically valuable co-products. High temperatures and pressures will be used to manipulate water’s properties, such as the dielectric constant, and hence its ability to solvate ions. Selective precipitation/recovery of valuable metal co-products creates new revenue streams (i.e., cost offsets). We will develop the scientific understanding necessary to control salt precipitation from supercritical seawater, inland brines, and water commonly co-produced with oil and gas. Atomistic simulations will be used to understand at a molecular level the changes in hydrogen bonding and ion solvation that occur as we increase the temperature and/or pressure and change the density and ion concentration. Our simulation results will also aid in the interpretation of experimental data and guide the development of applied thermodynamic models for engineering design calculations.

58 GEOSCIENCES↗
DOE OSTI2023

Solution Structures of Europium Terpyridyl Complexes with Nitrate and Triflate Counterions in Acetonitrile

Lanthanide–ligand complexes are key components of technological applications, and their properties depend on their structures in the solution phase, which are challenging to resolve experimentally or computationally. The coordination structure of the Eu 3+ ion in different coordination environments in acetonitrile is examined using ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. AIMD simulations are conducted for the solvated Eu 3+ ion in acetonitrile, both with or without a terpyridyl ligand, and in the presence of either triflate or nitrate counterions. EXAFS spectra are calculated directly from AIMD simulations and then compared to experimentally measured EXAFS spectra. In acetonitrile solution, both nitrate and triflate anions are shown to coordinate directly to the Eu 3+ ion forming either ten- or eight-coordinate solvent complexes where the counterions are binding as bidentate or monodentate structures, respectively. Coordination of a terpyridyl ligand to the Eu 3+ ion limits the available binding sites for the solvent and anions. In certain cases, the terpyridyl ligand excludes any solvent binding and limits the number of coordinated anions. The solution structure of the Eu-terpyridyl complex with nitrate counterions is shown to have a similar arrangement of Eu 3+ coordinating molecules as the crystal structure. Finally, this study illustrates how a combination of AIMD and EXAFS can be used to determine how ligands, solvent, and counterions coordinate with the lanthanide ions in solution.

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

Metadynamics simulations reveal mechanisms of Na + and Ca 2+ transport in two open states of the channelrhodopsin chimera, C1C2

Cation conducting channelrhodopsins (ChRs) are a popular tool used in optogenetics to control the activity of excitable cells and tissues using light. ChRs with altered ion selectivity are in high demand for use in different cell types and for other specialized applications. However, a detailed mechanism of ion permeation in ChRs is not fully resolved. Here, we use complementary experimental and computational methods to uncover the mechanisms of cation transport and valence selectivity through the channelrhodopsin chimera, C1C2, in the high- and low-conducting open states. Electrophysiology measurements identified a single-residue substitution within the central gate, N297D, that increased Ca 2+ permeability vs. Na + by nearly two-fold at peak current, but less so at stationary current. We then developed molecular models of dimeric wild-type C1C2 and N297D mutant channels in both open states and calculated the PMF profiles for Na + and Ca 2+ permeation through each protein using well-tempered/multiple-walker metadynamics. Results of these studies agree well with experimental measurements and demonstrate that the pore entrance on the extracellular side differs from original predictions and is actually located in a gap between helices I and II. Cation transport occurs via a relay mechanism where cations are passed between flexible carboxylate sidechains lining the full length of the pore by sidechain swinging, like a monkey swinging on vines. In the mutant channel, residue D297 enhances Ca 2+ permeability by mediating the handoff between the central and cytosolic binding sites via direct coordination and sidechain swinging. We also found that altered cation binding affinities at both the extracellular entrance and central binding sites underly the distinct transport properties of the low-conducting open state. This work significantly advances our understanding of ion selectivity and permeation in cation channelrhodopsins and provides the insights needed for successful development of new ion-selective optogenetic tools.

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

Collision-Induced Dissociation at TRIUMF's Ion Trap for Atomic and Nuclear science

The performance of high-precision mass spectrometry of radioactive isotopes can often be hindered by large amounts of contamination, including molecular species, stemming from the production of the radioactive beam. In this paper, we report on the development of Collision-Induced Dissociation (CID) as a means of background reduction for experiments at TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). This study was conducted to characterize the quality and purity of radioactive ion beams and the reduction of molecular contaminants to allow for mass measurements of radioactive isotopes to be done further from nuclear stability. Furthermore, this is the first demonstration of CID at an ISOL-type radioactive ion beam facility, and it is shown that molecular contamination can be reduced up to an order of magnitude.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ion Conductivity and Correlations in Model Salt-Doped Polymers: Effects of Interaction Strength and Concentration

Correlated anion and cation motion can greatly reduce the overall ion conductivity in electrolytes versus the ideal conductivity calculated based on the diffusion constants alone. Using coarse-grained molecular dynamics simulations, we calculate conductivity and the degree of uncorrelated ion motion in salt-doped homopolymers and block copolymers as a function of concentration and interaction strengths. Calculating conductivity from ion mobility under an applied electric field increases accuracy versus the typical use of fluctuation dissipation relationships in equilibrium simulations. In typical electrolytes, correlation in cation-anion motion is often expected to be reduced at low ion concentrations. Yet, for these polymer electrolytes with strong ion-polymer and ion-ion interactions, we find cation-anion motion is more correlated at lower concentrations when other variables are held constant. We show this phenomenon is related to the slower ion cluster relaxation rate at low concentrations rather than the static spatial state of ion aggregation or the fraction of free ions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electronic Structure Distortions in Chromium Chelates Impair Redox Kinetics in Flow Batteries

Aminopolycarboxylate chelates are emerging as a promising class of electrolyte materials for aqueous redox flow batteries, offering tunable redox potentials, solubility, and pH stability through careful selection of ligands and transition metal ions. Despite their potential, the impact of molecular structure modifications on the electronic and electrochemical properties of these chelates remains underexplored. Here, in this study, we examine how introducing a hydroxyl group, often employed for its solubilizing properties, to the backbone of CrPDTA, a reference chelate material, significantly changes the thermodynamics and kinetics of the chelate's redox process. We correlate changes in molecular and electronic structures to different electrochemical responses resulting from the hydroxyl addition and show that the introduction of this functional group leads to a distortion in the octahedral coordination of chromium. Furthermore, increased anisotropic spin density and nonintegral oxidation state changes in the Cr metal center result in a larger barrier for electron transfer in CrPDTA‐OH. It is demonstrated that preserving a hexacoordinate chelate structure across a broad pH range is crucial for efficient flow battery application and it is emphasized that ligand modifications must avoid distorting the octahedral coordination of the transition metal.

25 ENERGY STORAGE↗
DOE OSTI2023

Chemical Kinetic Method for Active-Site Quantification in Fe-N-C Catalysts and Correlation with Molecular Probe and Spectroscopic Site-Counting Methods

Mononuclear Fe ions ligated by nitrogen (FeN x ) dispersed on nitrogen-doped carbon (Fe-N-C) serve as active centers for electrocatalytic O 2 reduction and thermocatalytic aerobic oxidations. Despite their promise as replacements for precious metals in a variety of practical applications, such as fuel cells, the discovery of new Fe-N-C catalysts has relied primarily on empirical approaches. In this context, the development of quantitative structure–reactivity relationships and benchmarking of catalysts prepared by different synthetic routes and by different laboratories would be facilitated by the broader adoption of methods to quantify atomically dispersed FeN x active centers. In the present study, we develop a kinetic probe reaction method that uses the aerobic oxidation of a model hydroquinone substrate to quantify the density of FeN x centers in Fe- N-C catalysts. The kinetic method is compared with low-temperature Mössbauer spectroscopy, CO pulse chemisorption, and electrochemical reductive stripping of NO derived from NO 2 – on a suite of Fe-N-C catalysts prepared by diverse routes and featuring either the exclusive presence of Fe as FeN x sites or the coexistence of aggregated Fe species in addition to FeN x . The FeN x site densities derived from the kinetic method correlate well with those obtained from CO pulse chemisorption and Mössbauer spectroscopy. The broad survey of Fe-N-C materials also reveals the presence of outliers and challenges associated with each site quantification approach. Furthermore, the kinetic method developed here does not require pretreatments that may alter active-site distributions nor specialized equipment beyond reaction vessels and analytical instrumentation (e.g., NMR).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Directing polymorph specific calcium carbonate formation with de novo protein templates

Biomolecules modulate inorganic crystallization to generate hierarchically structured biominerals, but the atomic structure of the organic-inorganic interfaces that regulate mineralization remain largely unknown. We hypothesized that heterogeneous nucleation of calcium carbonate could be achieved by a structured flat molecular template that pre-organizes calcium ions on its surface. To test this hypothesis, we design helical repeat proteins (DHRs) displaying regularly spaced carboxylate arrays on their surfaces and find that both protein monomers and protein-C a2+ supramolecular assemblies directly nucleate nano-calcite with non-natural {110} or {202} faces while vaterite, which forms first in the absence of the proteins, is bypassed. These protein-stabilized nanocrystals then assemble by oriented attachment into calcite mesocrystals. We find further that nanocrystal size and polymorph can be tuned by varying the length and surface chemistry of the designed protein templates. Thus, bio-mineralization can be programmed using de novo protein design, providing a route to next-generation hybrid materials.

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

Bayesian inference of nuclear incompressibility from collective flow in mid-central Au+Au collisions at 400–1500 MeV/nucleon

The incompressibility K of symmetric nuclear matter (SNM) is determined through a Bayesian analysis of collective flow data from Au + Au collisions at beam energies $E = 400 -1500$ MeV/nucleon. This analysis utilizes a Gaussian process (GP) emulator applied to the isospin-dependent quantum molecular dynamics (IQMD) model for heavy-ion collisions, both with and without incorporating the momentum dependence of the single-nucleon potentials. Specifically, at the 68% confidence level, using rapidity and transverse velocity dependence of proton elliptic flow data with and without consideration of the momentum dependence, the inferred incompressibility values are $K=188.9^{+2.9}_{-4.5}$ MeV and $256.1^{+8.2}_{-8.7}$ MeV at $E = 400$ MeV/nucleon, respectively. When the transverse momentum dependence of proton-like directed flow data is included, the inferred incompressibility values become $K=222.3^{+9.0}_{-9.9}$ MeV and $K=285.5^{+6.7}_{-7.3}$ MeV, respectively. Furthermore, we found that the value of K derived from observables of proton elliptic flow increases with beam energy. Finally, this indicates that the equation of state (EoS) of nuclear matter hardens at higher densities and temperatures in reactions with higher beam energies.

Bayesian inference↗
DOE OSTI2021

Topical Report on Water Balance Model Development and Initial Application

This topical report summarizes results from developing a generic dynamic water balance tool which can be used by coal-fired power plant operators to evaluate FGD water management strategies. EPRI is using a statistical software model created on the GoldSim platform and data from coal-fired energy plants to determine the fraction of power plant water needs that can be satisfied by electrodialysis reversal treatment and, correspondingly, the discharge volume and/or byproducts that still need to be managed. The generic model is designed to accept a number of easily obtainable process variables and generate probabilities of resultant water characteristics. The model is also constructed as a universal process flow; one that can be adjusted to accommodate most major components of an FGD blowdown handling system. The generic model is dynamic and utilizes user-inputted operational variation in water quality and quantity to best provide a probabilistic representation of potential performance outcomes. In order to account for complexities of water composition and treatment, a separate chemical modeling software is used to evaluate water chemistry and unit operation performance. The thermodynamic modeling software, OLI Studio, is employed to convert ion data from laboratory analyses into balanced molecular concentrations that are classified as dissolved and suspended solids. The dynamic water model is designed to import a total of 11 species commonly found in FGD blowdown water as molecular concentrations.Water management unit operations are used in the model and categorized by either physical unit operation blocks or chemical unit operations blocks. As the name implies, the physical water treatment blocks are unit operations that are exclusively physical changes to the system (e.g., separations and volume changes). The chemical water treatment blocks are unit operations where chemical changes occur and require the use of OLI Studio to characterize subsequent streams. Water and material balances from three unique coal-fired power plants were obtained and preliminary modeling on the dynamic water balance tool was performed using the extracted data from these power plants. These three power plants were chosen due to their unique configurations and water profiles that represented a variety of operational changes. Results from using the dynamic tool with the three plant configurations are presented and discussed in the context of the existing tool and planned future capabilities.

01 COAL, LIGNITE, AND PEAT↗