Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium potential”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗
DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

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

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2020

Pre-equilibrium neutron emission in fission or fragmentation

Fissioning nuclei and fission fragments, nuclear fragments emerging from energetic collisions, or nuclei probed with various external fields can emit one or more pre-equilibrium neutrons, protons, and potentially other heavier nuclear fragments. Here, I describe a formalism which can be used to evaluate the pre-equilibrium neutron emission probabilities and the excitation energies of the remnant fragments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1984

Some minimization calculations concerning the effect of horizontal eddy momentum fluxes on the equilibrium zonal mean motions

A simplified linear two-level zonally averaged atmospheric model is used to determine the minimum zonal kinetic energy (ZKE) plus the zonal available potential energy (ZAPE) of equilibrium states for any distribution of horizontal eddy momentum fluxes. The model assumes fixed thermal forcing except for coefficients of the longwave radiative response to temperature. Some distinctive features of the model solutions are described. It is found that: the states of minimum ZKE or minimum ZKE + ZAPE contain a jet stream; the structure of the zonal winds is dependent on thermal forcing; the horizontal eddy momentum fluxes which minimize ZKE or ZKE + ZAPE can be up-gradient everywhere with respect to the upper level zonal mean zonal winds, or down-gradient depending on the thermal forcing; ZKE minimization results in a jet stream at 54 deg latitude independent of external parameters. The minimization of ZKE + ZAPE results in a jet stream at 35 deg latitude for a planet of the earth's radius and similar static stability.

Schneider, E. K.↗
NASA NTRS1990

IPC two-color analysis of x ray galaxy clusters

The mass distributions were determined of several clusters of galaxies by using X ray surface brightness data from the Einstein Observatory Imaging Proportional Counter (IPC). Determining cluster mass distributions is important for constraining the nature of the dark matter which dominates the mass of galaxies, galaxy clusters, and the Universe. Galaxy clusters are permeated with hot gas in hydrostatic equilibrium with the gravitational potentials of the clusters. Cluster mass distributions can be determined from x ray observations of cluster gas by using the equation of hydrostatic equilibrium and knowledge of the density and temperature structure of the gas. The x ray surface brightness at some distance from the cluster is the result of the volume x ray emissivity being integrated along the line of sight in the cluster.

White, Raymond E., III↗
DOE OSTI2021

Toward nonthermal control of excited quantum materials: framework and investigations by ultrafast electron scattering and imaging

Quantum material systems upon applying ultrashort laser pulses provide a rich platform to access excited material phases and their transformations that are not entirely like their equilibrium counterparts. The addressability and potential controls of metastable or long-trapped out-of-equilibrium phases have motivated interests both for the purposes of understanding the nonequilibrium physics and advancing the quantum technologies. Thus far, the dynamical spectroscopic probes eminently focus on microscopic electronic and phonon responses. For characterizing the long-range dynamics, such as order parameter fields and fluctuation effects, the ultrafast scattering probes offer direct sensitivity. Bridging the connections between the microscopic dynamics and macroscopic responses is central toward establishing the nonequilibrium physics behind the light-induced phases. Here, we present a path toward such understanding by cross-examining the structure factors associated with different dynamical states obtained from ultrafast electron scattering, imaging, and modeling. We give the basic theoretical framework on describing the non-equilibrium scattering problems and briefly describe how such framework relates to the out-of-equilibrium phenomena. We give effective models outlining the emergences of nonthermal critical points, hidden phases, and non-equilibrium relaxational responses from vacuum-suspended rare-earth tritellurides, tantalum disulfides thin films, and vanadium dioxide nanocrystalline materials upon light excitations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Ab Initio Characterization of Triatomic Bromine Molecules of Potential Interest in Stratospheric Chemistry

The equilibrium structures, harmonic vibrational frequencies, quadratic force fields, dipole moments, and IR intensities of several triatomic bromine compounds of known or potential importance in stratospheric ozone depletion chemistry have been determined using the CCSD(T) electron correlation method in conjunction with a basis set of triple zeta double polarized (TZ2P) quality. Specifically, the molecules included in the present study are HOBr, HBrO, FOBr, FBrO, BrNO, BrON, Br2O, BrBrO, BrCN, BrNC, ClOBr, ClBrO, and BrClO. Very accurate isomeric energy differences have also been determined at the CCSD(T) level with atomic natural orbital basis sets that include through g-type functions. In most cases, the isomer with a normal neutral Lewis dot structure is the lowest energy form, with the single exception that FBRO is predicted to be 11.1 kcal/mol (0 K) lower in energy than FOBr. In all cases, however, the hypervalent isomer is more stable relative to the isomer with a normal Lewis dot structure as compared to the chlorine analogs. Consistent with this observation, the energy of the last three molecules given above increases in the order ClOBr less than ClBrO less than BrClO. The CCSD(T)/TZ2P geometries and vibrational frequencies are in good agreement with the available experimental data. Heats of formation are determined for all species using a combination of theoretical isomeric, homodesmic, and isodesmic reaction energies. The accuracy of these quantities is ultimately dependent on the reliability of the experimental heat of formation of HOBr.

Lee. Timothy J.↗
NASA NTRS1995

Ab Initio Characterization of Triatomic Bromine Molecules of Potential Interest in Stratospheric Chemistry

The equilibrium structures, harmonic vibrational frequencies, quadratic force fields, dipole moments, and IR intensities of several triatomic bromine compounds of known or potential importance in stratospheric ozone depletion chemistry have been determined using the CCSD(T) electron correlation method in conjunction with a basis set of triple zeta double polarized (TZ2P) quality. Specifically, the molecules included in the present study are HOBr, HBrO, FOBr, FBrO, BrNO, BrON, Br2O, BrBrO, BrCN, BrNC, ClOBr, ClBrO, and BrClO. Very accurate isomeric energy differences have also been determined at the CCSD(T) level with atomic natural orbital basis sets that include through g-type functions. In most cases, the isomer with a normal neutral Lewis dot structure is the lowest energy form, with the single exception that FBrO is predicted to be 11.1 kcal/mol (0 K) lower in energy than FOBr. In all cases, however, the hypervalent isomer is more stable relative to the isomer with a normal Lewis dot structure as compared to the chlorine analogs. Consistent with this observation, the energy of the last three molecules given above increases in the order ClOBr less than ClBrO less than BrClO. The CCSD(T)/TZ2P geometries and vibrational frequencies are in good agreement with the available experimental data. Heats of formation are determined for all species using a combination of theoretical isomeric, homodesmic, and isodesmic reaction energies. The accuracy of these quantities is ultimately dependent on the reliability of the experimental heat of formation of HOBr.

Lee, Timothy J.↗
NASA NTRS1997

Response of the water status of soybean to changes in soil water potentials controlled by the water pressure in microporous tubes

Water transport through a microporous tube-soil-plant system was investigated by measuring the response of soil and plant water status to step change reductions in the water pressure within the tubes. Soybeans were germinated and grown in a porous ceramic 'soil' at a porous tube water pressure of -0.5 kpa for 28 d. During this time, the soil matric potential was nearly in equilibrium with tube water pressure. Water pressure in the porous tubes was then reduced to either -1.0, -1.5 or -2.0 kPa. Sap flow rates, leaf conductance and soil, root and leaf water potentials were measured before and after this change. A reduction in porous tube water pressure from -0.5 to -1.0 or -1.5 kPa did not result in any significant change in soil or plant water status. A reduction in porous tube water pressure to -2.0 kPa resulted in significant reductions in sap flow, leaf conductance, and soil, root and leaf water potentials. Hydraulic conductance, calculated as the transpiration rate/delta psi between two points in the water transport pathway, was used to analyse water transport through the tube-soil-plant continuum. At porous tube water pressures of -0.5 to-1.5 kPa soil moisture was readily available and hydraulic conductance of the plant limited water transport. At -2.0 kPa, hydraulic conductance of the bulk soil was the dominant factor in water movement.

NASA Center JSC↗
DOE OSTI2020

Effective pair potential for random fcc CoCrFeMnNi alloys

The single-phase equiatomic CoCrFeMnNi alloy is a random solid solution of five elements on the face-centered cubic lattice, whose pure constituents crystallize in very different structures and exhibit diverse magnetic properties. Due to the randomness of the alloy, 80% of nearest neighbor bonds are between unlike elements and thus the details of bonding in pure structures are less important. The elastic moduli of this alloy give rise to small Cauchy pressure C 12 - C 44 , which suggests that the dominant part of bonding may be described by a simple pair potential. We test this hypothesis by developing a long-range Lennard-Jones potential in which the equilibrium crystal structures of pure constituents are taken as reference. The standard mixing rules for regular solutions are then adopted to obtain parameters for bonds between unlike elements in the quinary system. The transferability of this potential to quaternary CoCrFeNi, ternary CoCrNi, and binary FeNi alloys is investigated and the predictions compared with experiments and density functional theory calculations. By sampling over a large number of random configurations, we investigate the effect of compositional randomness on misfit volumes, energies of point defects and stacking faults, and the dislocation friction stresses experienced by moving edge and screw dislocations.

Materials Science↗
NASA NTRS2012

Formation of Minor Phases in a Nickel-Based Disk Superalloy

The minor phases of powder metallurgy disk superalloy LSHR were studied. Samples were consistently heat treated at three different temperatures for long times to approximate equilibrium. Additional heat treatments were also performed for shorter times, to then assess non-equilibrium conditions. Minor phases including MC carbides, M23C6 carbides, M3B2 borides, and sigma were identified. Their transformation temperatures, lattice parameters, compositions, average sizes and total area fractions were determined, and compared to estimates of an existing phase prediction software package. Parameters measured at equilibrium sometimes agreed reasonably well with software model estimates, with potential for further improvements. Results for shorter times representing non-equilibrium indicated significant potential for further extension of the software to such conditions, which are more commonly observed during heat treatments and service at high temperatures for disk applications.

Gabb, T. P.↗
DOE OSTI2023

First-Principles-Based Machine Learning Models for Phase Behavior and Transport Properties of CO 2

In this work, we construct distinct first-principles-based machine-learning models of CO 2 , reproducing the potential energy surface of the PBE-D3, BLYP-D3, SCAN and SCAN-rvv10 approximations of density functional theory. Here we employ the Deep Potential methodology to develop the models and consequently achieve a significant computational efficiency over ab initio molecular dynamics (AIMD) that allows for larger system sizes and time scales to be explored. Although our models are trained only with liquid phase configurations, they are able to simulate a stable interfacial system and predict vapor-liquid equilibrium properties, in good agreement with results from the literature. Because of the computational efficiency of the models, we are also able to obtain transport properties, such as viscosity and diffusion coefficients. We find that the SCAN-based model presents a temperature shift in the position of the critical point, while the SCAN-rvv10-based model shows improvement but still exhibits a temperature shift that remains approximately constant for all properties investigated in this work. We find that the BLYP-D3-based model generally performs better for liquid phase and vapor-liquid equilibrium properties, but the PBE-D3-based model is better suited for predicting transport properties.

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

Kinetic and Thermodynamic Characterization of the Cobalt and Manganese Oxyhydroxide Cores Formed in Horse Spleen Ferritin

Horse spleen ferritin (HoSF) containing 800-1500 cobalt or 250-1200 manganese atoms as Co(O)OH and Mn(O)OH mineral cores within the HoSF interior (Co-HoSF and Mn-HoSF) was synthesized, and the chemical reactivity, kinetics of reduction, and the reduction potentials were measured. Microcoulometric and chemical reduction of HoSF containing the M(O)OH mineral core (M = Co or Mn) was rapid and quantitative with a reduction stoichiometry of 1.05+/-0.10 e/M forming a stable M(OH)2 mineral core. At pH 9.0, ascorbic acid (AH2), a two-electron reductant, effectively reduced the mineral cores; however, the reaction was incomplete and rapidly reached equilibrium. The addition of excess AH2 shifted the reaction to completion with a M(3+)/AH2 stoichiometry of 1.9-2.1, consistent with a single electron per metal atom reduction. The rate of reaction between M(0)OH and excess AH2 was measured by monitoring the decrease in mineral core absorbance with time. The reaction was first order in each reactant with second-order rate constants of 0.53 and 4.74/M/min, respectively, for Co- and Mn-HoSF at pH 9.0. From the variation of absorbance with increasing AH2 concentration, equilibrium constants at pH 9.0 of 5.0+/-1.9 for Co-HoSF and 2.9+/-0.9 for Mn-HoSF were calculated for 2M(O)OH + AH2 = 2M(OH)2 f D, where AH2 and D are ascorbic acid and dehydroascorbic acid, respectively. Consistent with these equilibrium constants, the standard potential for the reduction of Co(III)-HoSF is 42 mV more positive than that of the ascorbic acid reaction, while the standard potential of Mn(III)-HoSF is 27 mV positive relative to AH2. Fe(2+) in solution with Co- and Mn-HoSF under anaerobic conditions was oxidized to form Fe(O)OH within the HoSF interior, resulting in partial displacement of the Co or Mn by iron.

Zhang, Bo↗
DOE OSTI2024

How Precisely Can Individual Molecules Be Analyzed? A Case Study on Locally Quantifying Forces and Energies Using Scanning Probe Microscopy

Recent advances in scanning probe microscopy methodology have enabled the measurement of tip-sample interactions with picometer accuracy in all three spatial dimensions, thereby providing a detailed site-specific and distance-dependent picture of the related properties. This paper explores the degree of detail and accuracy that can be achieved in locally quantifying probe-molecule interaction forces and energies for adsorbed molecules. Toward this end, cobalt phthalocyanine (CoPc), a promising CO 2 reduction catalyst, was studied on Ag(111) as a model system using low-temperature, ultrahigh vacuum noncontact atomic force microscopy. Data were recorded as a function of distance from the surface, from which detailed three-dimensional maps of the molecule's interaction with the tip for normal and lateral forces as well as the tip-molecule interaction potential were constructed. The data were collected with a CO molecule at the tip apex, which enabled a detailed visualization of the atomic structure. Determination of the tip-substrate interaction as a function of distance allowed isolation of the molecule-tip interactions; when analyzing these in terms of a Lennard-Jones-type potential, the atomically resolved equilibrium interaction energies between the CO tethered to the tip and the CoPc molecule could be recovered. Interaction energies peaked at less than 160 meV, indicating a physisorption interaction. As expected, the interaction was weakest at the aromatic hydrogens around the periphery of the molecule and strongest surrounding the metal center. In conclusion, the interaction, however, did not peak directly above the Co atom but rather in pockets surrounding it.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atmospheric-pressure ammonia synthesis on AuRu catalysts enabled by plasmon-controlled hydrogenation and nitrogen-species desorption

The Haber–Bosch process for ammonia synthesis contributes up to ~3% of global greenhouse gas emissions. Plasmonic catalysts strongly concentrate light and can alter the reaction intermediates via out-of-equilibrium processes, providing the potential for an alternative, less-energy-intensive pathway to synthesize ammonia. Here, in this study, we show that gold-ruthenium (AuRu) bimetallic nanoparticles can synthesize ammonia at room temperature and pressure using visible light. We create AuRu alloys with varying compositions and achieve ammonia production rates of ~60 μmol per gram of catalyst bed per hour. In situ infrared spectroscopy reveals that light accelerates the hydrogenation of nitrogen intermediates compared to conventional thermal catalysis. Through computational modelling, we demonstrate that photo-excited electrons enable associative hydrogenation pathways for nitrogen activation rather than direct nitrogen–nitrogen bond breaking. This light-assisted mechanism requires both hydrogen and light working together to overcome the nitrogen activation barrier, mimicking how biological enzymes produce ammonia naturally and providing fundamental insights for developing sustainable, energy-efficient chemical synthesis.

Yuan, Lin [Stanford Univ., CA (United States)] (OR↗
DOE OSTI2021

Preventing leaching from lead water pipes with electrochemistry: an exploratory study

Toxic levels of lead leaching from ageing water distribution infrastructure affect over 5000 public drinking water systems in the US. Pipe replacement, the most effective solution to this problem, is prohibitively expensive. Additionally, chemical conditioning of drinking water using orthophosphates, although cost-effective, does not quickly stop lead leaching once it has started. We propose a novel approach to stop lead leaching: to rapidly form an insoluble scale within lead pipes using an external power supply. We report on the feasibility of this approach by first anodizing lead coupons and lead pipes reclaimed from a local water distribution using a phosphate electrolyte and different potentials, pH values, and phosphate concentrations. We subsequently exposed these anodized lead coupons and pipes to synthetic tap water to evaluate their lead leaching rates. We found that polarizing lead coupons in the presence of a 0.05 M phosphate solution decreased lead leaching by up to a 100-fold, relative to leaching from polished bare lead. Similarly, polarizing the reclaimed lead pipes (with a preexisting scale) decreased lead leaching from an average of 36 ppb to 7 ppb. These results were observed when applying potentials that favor Pb(IV) formation, which resulted in the buildup of PbO 2 and Pb 5 (PO 4 ) 3 OH in both lead coupons and pipes. Our findings indicate that this novel technology has the potential to rapidly decrease equilibrium lead levels in tap water below the EPA action limit of 15 ppb, and thus deserves further exploration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ground and excited state properties of ThB − and ThB: a theoretical study

In the present work, we studied a series of electronic and spin–orbit states of ThB − and ThB using high-level multireference and coupled-cluster theories. We report the potential energy curves (PECs), equilibrium electron configurations, spectroscopic constants, energetics, and spin–orbit coupling effects of 17 and 19 electronic states of ThB − and ThB, respectively. The ground state of ThB − is a single-reference 1 3 Π with a 1σ 2 2σ 2 3σ 1 1π 3 electron configuration. Detachment of an electron from the doubly occupied 1π orbital of ThB − (1 3 Π) produces the single-reference ground electronic state of ThB (1 4 Σ − ). The ground spin–orbit states of ThB − and of ThB are 1 3 Π + 0 and 1 4 Σ − 3/2 , respectively. The vertical electron detachment energy (VDE) of ThB − and the adiabatic electron attachment energy (AEA) of ThB at our largest CBS-C-CCSD(T)+δT(Q)+δSO (complete basis set effect, spin–orbit effect, and triple and perturbative quadruple electron correlation effect added coupled-cluster theory with single, double, and perturbative triple excitations) level are 1.473 eV and 1.459, respectively. The reaction of Th( 3 F) + B( 2 P o ) produces the ground state of ThB with a bond energy of 2.843 eV. Finally, we estimated a heat of formation, Δ H 0 f (298 K), of 891.01 kJ mol −1 for the ThB molecule. The high-level findings of this work are expected to aid and motivate future experimental spectroscopic investigations of ThB and ThB − species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗