Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Charge equilibration models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

35 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Probing electronic dead layers in homoepitaxial n-SrTiO 3 (001) films

We combine state-of-the-art oxide epitaxial growth by hybrid molecular beam epitaxy with transport, x-ray photoemission, and surface diffraction, along with classical and first-principles quantum mechanical modeling to investigate the nuances of insulating layer formation in otherwise high-mobility homoepitaxial n-SrTiO 3 (001) films. Our analysis points to charge immobilization at the buried n-SrTiO 3 /undoped SrTiO 3 (001) interface as well as within the surface contamination layer resulting from air exposure as the drivers of electronic dead-layer formation. As Fermi level equilibration occurs at the surface and the buried interface, charge trapping reduces the sheet carrier density (n 2D ) and renders the n-STO film insulating if n 2D falls below the critical value for the metal-to-insulator transition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Charged particle transport coefficient challenges in high energy density plasmas

High energy density physics (HEDP) and inertial confinement fusion (ICF) research typically relies on computational modeling using radiation-hydrodynamics codes in order to design experiments and understand their results. These tools, in turn, rely on numerous charged particle transport and relaxation coefficients to account for laser energy absorption, viscous dissipation, mass transport, thermal conduction, electrical conduction, non-local ion (including charged fusion product) transport, non-local electron transport, magnetohydrodynamics, multi-ion-species thermalization, and electron-ion equilibration. In many situations, these coefficients couple to other physics, such as imposed or self-generated magnetic fields. Furthermore, how these coefficients combine are sensitive to plasma conditions as well as how materials are distributed within a computational cell. Uncertainties in these coefficients and how they couple to other physics could explain many of the discrepancies between simulation predictions and experimental results that persist in even the most detailed calculations. This paper reviews the challenges faced by radiation-hydrodynamics in predicting the results of HEDP and ICF experiments with regard to these and other physics models typically included in simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Local spin polarizations in relativistic heavy-ion collisions

Here, based on a generalized side-jump formalism for massless chiral fermions, which naturally takes into account the spin-orbit coupling in the scattering of two chiral fermions and the chiral vortical effect in a rotating chiral fermion matter, we have developed a covariant and total-angular-momentum-conserved chiral transport model to study both the global and local polarizations of this matter. For a system of massless quarks of random spin orientations and finite vorticity in a box, we have demonstrated that the model can exactly conserve the total angular momentum of the system and dynamically generate the quark spin polarization expected from a thermally equilibrated quark matter. Using this model to study the spin polarization in relativistic heavy-ion collision, we have found that the local quark spin polarizations depend strongly on the reference frame where they are evaluated as a result of the nontrivial axial charge distribution caused by the chiral vortical effect. We have further shown that because of the anomalous orbital or side-jump contribution to the quark spin polarization, the local quark polarizations calculated in the medium rest frame are qualitatively consistent with the local polarizations of Lambda hyperons measured in experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Liquid Sorption-Enhanced Haber–Bosch Process

The use of a liquid sorbent in a traditional Haber-Bosch process enables significant improvements in energy efficiency and potential cost savings for arguably the most important chemical process on the planet. The approach presented in this report employs an incompressible liquid sorbent that absorbs and releases ammonia (NH 3 ) under specific conditions. To achieve this, we investigate reactions of ammonia and pure phosphoric acid (H 3 PO 4 , PA), which rapidly neutralize to form an equilibrated solution of monoammonium phosphate (MAP) and diammonium phosphate (DAP) that functions as a reversible and regenerable sorbent. Through intimate contact of the gas-phase Haber-Bosch reaction mixture with this liquid absorbent, complete equilibrium uptake may be achieved in an appropriately sized separator, and facile separation occurs through the use of independent liquid and gas phases. Following depressurization and release of the ammonia product, only the incompressible fluid needs to be repressurized and returned to the reactor. This study documents proof-of-concept absorption and desorption experiments carried out in 75 mL batch reactors, predominantly charged with precise MAP and DAP mixtures that equilibrate at process-relevant temperatures and pressures. We then assemble the first thermodynamic relationships that underlie this advantaged separation strategy, validated by reactive force field (ReaxFF) interatomic potential simulations, and benchmarked with traditional separation routes via process modeling and technoeconomic analysis. The scale of energy consumption in the century-old Haber-Bosch process is massive, and the elegant liquid sorption approach reported here offers opportunities to enhance its energy efficiency for the next frontier of ammonia synthesis.

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

Teaching a neural network to attach and detach electrons from molecules

Abstract Interatomic potentials derived with Machine Learning algorithms such as Deep-Neural Networks (DNNs), achieve the accuracy of high-fidelity quantum mechanical (QM) methods in areas traditionally dominated by empirical force fields and allow performing massive simulations. Most DNN potentials were parametrized for neutral molecules or closed-shell ions due to architectural limitations. In this work, we propose an improved machine learning framework for simulating open-shell anions and cations. We introduce the AIMNet-NSE (Neural Spin Equilibration) architecture, which can predict molecular energies for an arbitrary combination of molecular charge and spin multiplicity with errors of about 2–3 kcal/mol and spin-charges with error errors ~0.01e for small and medium-sized organic molecules, compared to the reference QM simulations. The AIMNet-NSE model allows to fully bypass QM calculations and derive the ionization potential, electron affinity, and conceptual Density Functional Theory quantities like electronegativity, hardness, and condensed Fukui functions. We show that these descriptors, along with learned atomic representations, could be used to model chemical reactivity through an example of regioselectivity in electrophilic aromatic substitution reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quark flavor equilibration of the quark-gluon plasma

The early stage of a heavy-ion collision is marked by rapid entropy production and the transition from a gluon saturated initial condition to a plasma of quarks and gluons that evolves hydrodynamically. However, during the early times of the hydrodynamic evolution, the chemical composition of the QCD medium is still largely unknown. We present a study of quark chemical equilibration in the (Q)GP using a novel model of viscous hydrodynamic evolution in partial chemical equilibrium. Motivated by the success of gluon saturated initial condition models, we initialize the QCD medium as a completely gluon dominated state. Local quark production during the hydrodynamic phase is then simulated through the evolution of time-dependent fugacities for each independent quark flavor, with the timescales set as free parameters to compare different rates of equilibration. We present the results of complete heavy-ion collision simulations using this partial chemical equilibrium model, and show the effects on hadronic and electromagnetic observables. In particular, we show that the development of flow is sensitive to the equilibration timescale, providing an empirical way to probe the chemical equilibration of the QCD medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Relaxation times for disoriented isospin condensates in high energy heavy ion collisions

Fluctuations between charged and neutral kaons measured by the ALICE Collaboration in Pb-Pb collisions at the CERN Large Hadron Collider (LHC) exceed conventional explanations. Previously it was shown that if the scalar condensate is accompanied by an electrically neutral isospin-1 field then the combination can produce large equilibrium fluctuations where ⟨$\overline{𝑢}$ ⁢𝑢⟩≠⟨ $\overline{𝑑}$𝑑⟩. Hadronizing strange and antistrange quarks might then strongly fluctuate between charged ($\overline{𝑢}$⁢𝑠 or $\overline{𝑠}$⁢𝑢) and neutral (𝑑⁢$\overline{𝑠}$⁢ or 𝑠⁢$\overline{𝑑}$) kaons. Here, we estimate the times for the condensates to achieve their equilibrium probability distributions within causal volumes in high energy heavy ion collisions. This is achieved by modeling the temperature dependence of the condensates, mesonic collective excitations, decay rates of the associated fields, and employing the Langevin and Fokker-Planck equations. Within this model, we find that the equilibration times are short compared with the expansion time, making disoriented isospin condensates a viable explanation for the anomalous fluctuations observed at the LHC.

Quark-gluon plasma↗
DOE OSTI2024

Review of the second charged-particle transport coefficient code comparison workshop

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24–27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron–ion properties such as the electron–ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Aqueous ion partitioning in Nafion: Applicability of Manning's counter-ion condensation theory

Nafion is a well-studied polymer of significant interest as a fuel cell membrane due to rapid transport of protons, which is often ascribed, at least in part, to its phase-separated morphology. The literature detailing ion sorption and transport in Nafion diverges considerably from that of water purification membranes, with contrasting analyses and models. To help bridge these gaps and facilitate comparisons of membranes with various structures, this study presents measurements of ion activity coefficients and sorption in Nafion 117 equilibrated with aqueous solutions of NaCl, MgCl 2 , and Na 2 SO 4 , which are relevant to water purification. The experimental data are compared to a theoretical framework based on Manning's counter-ion condensation theory and Donnan theory. Good agreement between the model and data is observed for NaCl and MgCl 2 , using no adjustable parameters. However, the model significantly underpredicts Na 2 SO 4 solubility, which is attributed to ion pairing effects. Furthermore, the results highlight the general applicability of counter-ion condensation theory to highly charged polymers, including those exhibiting nanoscale ordering, while revealing fundamental limitations related to ion pairing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1979

Relaxation to photochemical equilibrium of charged species within displaced air parcels in the D-region

The effects of the relaxation to photochemical equilibrium are investigated for the electrons and NO(+) ions in the upper D-region and for the electrons and several species of positive and negative ions present in the lower D-region. In the upper D-region the photochemical model includes electrons, nitric oxide, and NO(+). A self-consisting model of the equilibrium concentrations is developed for the lower D-region, which includes electrons, nitric oxide, O2(-), heavy negative and positive ions, and the oxygen species O, O3, and O2(1 delta g). At altitudes of 75 and 80 km the time dependence of the oxygen species has a significant effect on the equilibration of the electron and O2(-) concentrations. The relaxation to photochemical equilibrium proceeds by five e-folding time scales in the lower D-region, where the time required for final equilibration due to recombination is very long.

Hill, R. J.↗
DOE OSTI2023

Tracking Photochemical and Photophysical Processes for Solar Energy Conversion Via Multidimensional Visible and Vibrational Spectroscopic Methods

This project was motivated by an overarching goal to elucidate the mechanism of energy and electron transfer that governs the efficient charge separation in photosystem I (PSI) complexes. PSI is a natural light harvesting complex that drives oxygenic photosynthesis in plants, algae, and cyanobacteria. It uses ~300 tightly packed chlorophylls (Chls) to absorb photons, transfer the excitation energy to the reaction center (RC), and generate a charge separated state with near unity quantum efficiency (QE). A better understanding of the mechanism of energy transfer and charge separation in PSI is required for understanding the high QE of natural light harvesting complexes, and it could lead to the further development of artificial photosynthetic systems for solar energy conversion. We applied two-dimensional electronic spectroscopy (2DES) to different cyanobacterial photosystem I complexes to map energy transfer pathways to gain insight into the efficient light harvesting of PSI. We interpreted the 2DES spectra through a global analysis procedure, finding two ultrafast equilibration processes involving red-shifted Chl states. The analysis procedure is general and can be readily applied to other natural and artificial light harvesting systems. In addition to investigating PSI complexes, we also performed measurements on model systems to establish general data analysis procedures for interpreting 2D spectra.

14 SOLAR ENERGY↗
DOE OSTI2022

The Donnan potential revealed

Selective transport of solutes across a membrane is critical for many biological, water treatment and energy conversion and storage systems. When a charged membrane is equilibrated with an electrolyte, an unequal distribution of ions arises between phases, generating the so-called Donnan electrical potential at the solution/membrane interface. The Donnan potential results in the partial exclusion of co-ion, providing the basis of permselectivity. Although there are well-established ways to indirectly estimate the Donnan potential, it has been widely reported that it cannot be measured directly. Here we report the first direct measurement of the Donnan potential of an ion exchange membrane equilibrated with salt solutions. Our results highlight the dependence of the Donnan potential on external salt concentration and counter-ion valence, and show a reasonable agreement with current theoretical models of IEMs, which incorporate ion activity coefficients. By directly measuring the Donnan potential, we eliminate ambiguities that arise from limitations inherent in current models.

25 ENERGY STORAGE↗
DOE OSTI2024

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

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

Deciphering the excited-state landscape of Cr(iii) tris(diimines) using [Cr(phen) 3 ] 3+

Earth-abundant Cr(III) polypyridines chromophores feature long-lived doublet ligand-field excited states. The thermally equilibrated 2 T 1 / 2 E ligand-field excited state manifold produces broad, structured UV-Vis-NIR absorption bands. Although ubiquitous across pseudo-octahedral Cr(III) diimine complexes, the dynamics and individual electronic transitions underlying these spectra have remained elusive. Here, we combine time-resolved spectroscopy and multiconfigurational calculations to assign ground- and excited-state absorption spectra of Cr(III) tris(1,10-phenanthroline), [Cr(phen) 3 ] 3+ , as a prototypical example for this class of complexes. Excited-state kinetics reveal a rapid biexponential decay of the initially prepared intraligand/ligand-to-metal charge transfer [ 4 ( 1 IL)/ 4 LMCT] excited state into the lowest 2 T 1 / 2 E state within 10 picoseconds, assigned to intersystem crossing and vibrational relaxation. Calculations identify broad excited-state absorption features associated with metal-centered ligand-field states ( 2 MC) mixed with intraligand ( 2 IL) transitions within the doublet manifold. These results provide detailed spectral assignments and highlight the need to include multiconfigurational ligand-based orbitals when modeling the spectral properties of Cr(III) coordination complexes.

Barth, Alexandra T. [North Carolina State Universi↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗