SEARCH · Engineering Papers

Results for “Equilibrium constant”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Hydrodynamic Models of Line-Driven Accretion Disk Winds III: Local Ionization Equilibrium

We present time-dependent numerical hydrodynamic models of line-driven accretion disk winds in cataclysmic variable systems and calculate wind mass-loss rates and terminal velocities. The models are 2.5-dimensional, include an energy balance condition with radiative heating and cooling processes, and includes local ionization equilibrium introducing time dependence and spatial dependence on the line radiation force parameters. The radiation field is assumed to originate in an optically thick accretion disk. Wind ion populations are calculated under the assumption that local ionization equilibrium is determined by photoionization and radiative recombination, similar to a photoionized nebula. We find a steady wind flowing from the accretion disk. Radiative heating tends to maintain the temperature in the higher density wind regions near the disk surface, rather than cooling adiabatically. For a disk luminosity L (sub disk) = solar luminosity, white dwarf mass M(sub wd) = 0.6 solar mass, and white dwarf radii R(sub wd) = 0.01 solar radius, we obtain a wind mass-loss rate of M(sub wind) = 4 x 10(exp -12) solar mass yr(exp -1) and a terminal velocity of approximately 3000 km per second. These results confirm the general velocity and density structures found in our earlier constant ionization equilibrium adiabatic CV wind models. Further we establish here 2.5D numerical models that can be extended to QSO/AGN winds where the local ionization equilibrium will play a crucial role in the overall dynamics.

Pereyra, Nicolas Antonio↗

NASA NTRS1928

The Gaseous Explosive Reaction : The Effect of Inert Gases

Attention is called in this report to previous investigations of gaseous explosive reactions carried out under constant volume conditions, where the effect of inert gases on the thermodynamic equilibrium was determined. The advantage of constant pressure methods over those of constant volume as applied to studies of the gaseous explosive reaction is pointed out and the possibility of realizing for this purpose a constant pressure bomb mentioned. The application of constant pressure methods to the study of gaseous explosive reactions, made possible by the use of a constant pressure bomb, led to the discovery of an important kinetic relation connecting the rate of propagation of the zone of explosive reaction within the active gases, with the initial concentrations of those gases: s = K(sub 1)(A)(sup n1)(B)(sup n2)(C)(sup n3)------. By a method analogous to that followed in determining the effect of inert gases on the equilibrium constant K, the present paper records an attempt to determine their kinetic effect upon the expression given above.

Stevens, F W↗

DOE OSTI2024

Coupled feldspar dissolution and secondary mineral precipitation in batch systems: 6. Labradorite dissolution, calcite growth, and clay precipitation at 60 °C and pH 8.2–8.4

Here, we conducted experiments on concurrent labradorite dissolution, calcite precipitation, and clay precipitation in batch reactor systems and tracked reaction processes using multiple isotope tracers. Labradorite was chosen for its role as a major and reactive component in basalt; the experiments thus directly impact our understanding of CO 2 storage in basalt aquifers and enhanced rock weathering. We doped initial solutions with 29 Si, 43 Ca, and Ca 13 CO 3 (s). Experiments were conducted at 60 °C and pH ~ 8.3 for up to 840 h, with isotope ratios in the experimental aqueous solutions measured using MC-ICP-MS. Unidirectional rates of labradorite dissolution near equilibrium were approximately two orders of magnitude slower than far-from-equilibrium rates reported in the literature. Calcite growth occurred near equilibrium and the rates were limited by the labradorite dissolution rates. In the steady state phase, the interplay of these three heterogeneous reactions—labradorite dissolution, calcite growth, and clay precipitation—results in a coupled system that approaches a near-equilibrium state. The system does not reach true equilibrium because labradorite continues to dissolve, albeit at a much slower rate near equilibrium. The overall reaction can be approximated as, Na 0.4 Ca 0.6 Al 1.6 Si 2.4 O 8 + 0.6HCO 3 - + 1·.7H 2 O + 0.4H + → 0.4Na + + 0.6CaCO 3(s) + 0.5Al 2 Si 2 O 5 (OH) 4(s) + 0.6Al(OH) 4 - + 1.4SiO 2 o (aq). The experimental results show that using short-term far-from-equilibrium rate constants would lead to an overestimation of feldspar weathering rates at the Earth’s surface (e.g., basalt weathering and enhanced rock weathering) and CO 2 mineralization in basalt aquifers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Coordination-Induced Weakening of N–H Bonds Driven by Bimetallic Cooperativity in Zr/Co Compounds

The bond dissociation free energy (BDFE) of the element-hydrogen bonds of protic substrates have been found to decrease upon metal coordination. Herein, an early/late heterobimetallic complex is used to examine the impact on the BDFE N−H when the substrate binding site and the redox-active site are two different metals that are spatially separated. A tris- (phosphinoamide) framework is used to link a d 0 ZrIV center with an accessible substrate binding site to a coordinatively saturated redox-active Co center, which serves as an appended electron reservoir. A series of aniline, amido, and imido Zr/Co model compounds were synthesized starting from the Zr IV /Co −I aniline adduct PhH 2 N−Zr(MesNP i Pr 2 ) 3 CoCN t Bu (2). 2,4,6-tristert- butylphenoxyl radical ( t Bu 3 ArO • ) was used to abstract one or two H atoms and produce the amido and imido complexes PhHNZr( MesNP i Pr 2 ) 3 CoCN t Bu (3) and PhN≡Zr(MesNP i Pr 2 ) 3 CoCNtBu (4), respectively. Using open-circuit potential measurements, the BDFE N−H within 2 and 3 were determined to be 37 kcal/mol (2) and 55 kcal/mol (3). Cyclic voltammetry measurements were conducted to determine the Co I/0 and Co 0/−I redox potentials. The pK a s were then estimated using the Bordwell equation to provide further insight into the thermochemical aspects of the observed proton coupled electron transfer (PCET) reactions.

Bond dissociation free energy↗

DOE OSTI2026

Radiation Damage Mechanisms in Oxide-Based Solid Electrolytes

Ion implantation is a promising approach to mitigate dendrite formation in metal anode solid-state batteries by modifying the chemical and mechanical properties of oxide-based electrolyte surfaces such as Li 7 La 3 Zr 2 O 12 (LLZO). Yet, the detailed atomistic mechanisms and potentially adverse side effects of ion implantation remain poorly understood. In this study, we used molecular dynamics simulations to elucidate the evolution of radiation-induced damage in LLZO. The results reveal that radiation damage is primarily driven by recoils from heavy ions in the lattice, producing antisite defect clusters whose density decreases as recoil energy increases, a counterintuitive finding that highlights complex cascade fragmentation. Radiation-induced disruption to the crystal sublattice connectivity, especially in the La–O network critical for Li-ion transport, leads to degradation in ionic pathways at low energies, while higher energy cascades can form new conduction pathways via enhanced Zr–O network connectivity. These findings provide a mechanistic foundation for optimizing ion implantation strategies that balance disruption of the lattice for surface modification for dendrite resistance with preservation of essential ionic conductivity.

Defects↗

DOE OSTI2025

Reliable p K a Prediction through Efficient Incorporation of Anharmonicity within the Nuclear–Electronic Orbital Framework

Accurate pK a prediction is critical for understanding chemical reactivity and molecular properties across a wide range of applications. Computational approaches usually invoke a harmonic treatment of the vibrational modes for zero-point energies, as well as thermal and entropic contributions. Herein, we present a general protocol for relative pK a prediction that incorporates the significant anharmonic effects using nuclear–electronic orbital (NEO) theory. This protocol is validated against experimental data for a range of molecules in acetonitrile, including protonated nitrogen bases, nitrophenols, anilines, and diamines, as well as cobalt electrocatalysts. For simple acids, the NEO approach offers only a slight improvement over conventional density functional theory with the standard harmonic vibrational treatment, whereas for hydrogen-bonded acids, the NEO approach offers more significantly improved performance at a comparable computational cost. This accessible methodology provides a practical route for accurate pKa prediction in challenging systems and is extendable to related thermodynamic properties such as hydricities and proton-coupled redox potentials.

Density functional theory↗

DOE OSTI2025

The impact of curation errors in the PDBBind Database on machine learning predictions of protein–protein binding affinity

The PDBBind database has been widely utilized for the computational prediction of protein–protein binding affinities. While the accuracy of the PDBBind-curated equilibrium dissociation constants (K D ) has been reported for the protein–ligand subset of the PDBBind database, the curation accuracy has not been reported for the protein–protein subset. Here, we present a detailed manual analysis for the subset of PDBBind records with PubMed Central Open Access primary publications and find that ~19% of these records had K D values that were not supported by their primary publications. The impact of these putative curation errors on the machine learning-based prediction of K D from experimental protein–protein 3D structures was evaluated and correcting the curation errors improved the Pearson correlation coefficient between measured and random forest-predicted log 10 (K D ) values by ~8 percentage points. This finding underscores the importance of dataset accuracy for computational modelling and highlights the need for more stringent curation processes when extracting information from the scientific literature.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Mapping Interfacial Solution Forces that Drive Fluorescent Nanoparticle Incorporation into Crystals

Interfacial solution structure governs processes ranging from catalytic and electrochemical reactions to particle aggregation and crystallization. Here, we design bio-inspired nanocomposites by mapping interfacial forces and controlling nanoparticle incorporation during crystal growth. Our analytical model contains measured kinetic barriers for surface approaches and equilibrium binding constants with the crystal surface. We validate this model using fluorescent silica nanoparticles and calcite. Our results show that surface chemistry dictates incorporation pathways: methoxy groups are least favorable, hydroxyls reduce kinetic barriers by interacting favorably with hydration layers, carboxylates bind strongly but must overcome a kinetic barrier, and amines outperform all kinetically and thermodynamically. This framework resolves the interplay between interfacial solution structures, particle forces and dynamics, and growth kinetics, enabling predicative control of nanocomposite formation and multiplexed mixed-particle systems.

Zhang, Mingyi [Pacific Northwest National Laborato↗

NASA NTRS1971

A system debugging model.

Consideration of the nature of the 'debugging' process applied to a new complex system during the initial period of its life. During this period failures and errors are corrected as they occur, with resulting improvement in the subsequent performance of the system. One mathematical idealization of this process leads to the assumption that system failure rate is decreasing with time. In practice, the debugging phase is considered completed when the failure rate reaches an equilibrium or constant value. Models are formulated for this phenomenon. Maximum likelihood estimates are obtained for relevant failure rate functions and for the end of the debugging period. A conservative upper confidence bound on the stable failure rate is obtained.

Barlow, R. E.↗

NASA NTRS1982

The numerical analysis of the rotational theory for the formation of lunar globules

The morphology of lunar globules is studied through the application of a numerical analysis of their rotation in space during cooling. It is assumed that molten rock is shot from the surface of the moon, solidifies in space above the moon and then falls back to the surface. The rotational theory studied makes the following assumptions: the volume of the molten rock does not change during cooling; the angular momentum is conserved; there are no internal motions because of the high viscosity of the molten rock, i.e., in equilibrium the globule is rotating as a rigid body; finally, the kinetic reaction of the globule to the forces is fast relative to the rate of cooling, i.e., the globule reaches equilibrium at constant energy. These assumptions are subjected to numerical analysis yielding good agreement between the actual globule shapes and the numerical results, but leaving some doubt as to the validity of the rotational theory due to the failure to establish the existence of true local minima and an incomplete understanding of the thermokentics.

Ross, J.↗

NASA NTRS1981

Effects of diffusion and mass flows on C IV and Si IV lines formed in the solar atmosphere

A model for the transition region is derived from an initial interpretation of EUV observations, assuming ionization equilibrium and constant elemental abundance with height. The effects of diffusion and mass flows are then included in the initial model and the emergent profiles of several C IV and Si IV lines are computed. It is found that diffusion and mass flows have a strong effect on both the emergent intensity and spectral shape of these lines. Diffusion acts to deplete the transition region of heavy ions to an extent which depends on the detailed temperature and density structure. The net effect is a weighting of the coronal emission relative to that in the transition. On the basis of the results, it is suggested that the downflows observed in the network in lines of C IV and Si IV could be due to gravitational settling of the ions following their injection, via diffusion, into the corona from spicules. When flows are superposed on the basic diffusion model, the rapid change in elemental abundance characteristic of the model is virtually eliminated and the coronal contribution to the emission in C IV and Si IV lines becomes negligible relative to that in the transition region. Flows have a strong effect on the computed line intensities and introduce large asymmetries into the line shapes.

Roussel-Dupre, R.↗

NASA NTRS1989

Periodic-disturbance accommodating control of the Space Station for asymptotic momentum management

Periodic-disturbance accommodating control is investigated for asymptotic momentum management of control moment gyros used as primary actuating devices for the Space Station. The proposed controller utilizes the concepts of quaternion feedback control and periodic-disturbance accommodation to achieve oscillations about the constant torque equilibrium attitude, while minimizing the control effort required. Three-axis coupled equations of motion, written in terms of quaternions, are derived for roll/yaw controller design and stability analysis. The quaternion feedback controller designed using the linear-quadratic regulator synthesis technique is shown to be robust for a wide range of pitch angles. It is also shown that the proposed controller tunes the open-loop unstable vehicle to a stable oscillatory motion which minimizes the control effort needed for steady-state operations.

Warren, Wayne↗

NASA NTRS1989

Periodic-disturbance accommodating control of the space station for asymptotic momentum management

Warren, Wayne↗

NASA NTRS1989

Periodic-disturbance accommodating control of the space station for asymptotic momentum management

Periodic maneuvering control is developed for asymptotic momentum management of control gyros used as primary actuating devices for the Space Station. The proposed controller utilizes the concepts of quaternion feedback control and periodic-disturbance accommodation to achieve oscillations about the constant torque equilibrium attitude, while minimizing the control effort required. Three-axis coupled equations of motion, written in terms of quaternions, are derived for roll/yaw controller design and stability analysis. It is shown that the quaternion feedback controller is very robust for a wide range of pitch angles. It is also shown that the proposed controller tunes the open-loop unstable vehicle to a stable oscillatory motion which minimizes the control effort needed for steady-state operations.

Warren, Wayne↗

NASA NTRS1990

Application of the embedded atom method to Pb and Be

We have derived the embedding energy functional and two-body potential of the embedded atom method using decreasing exponentials for both the electron density and the two body potential. The embedding function was obtained from the equation of state given by Rose et al. (1984). Because of the form of the embedding function, the equilibrium lattice constant, cohesive energy, and bulk modulus are automatically satisfied. The two parameters phi(e) and gamma of the two-body potential were determined by fitting to shear modulus and the single vacancy formation energy. Contributions of up to the third nearest neighbors were included in the evaluation of the charge density rho and the two-body potential phi. The stability and anisotropy of each structure were estimated and compared with the available experimental data.

Karimi, M.↗

NASA NTRS1995

An Accurate ab initio Quartic Force Field and Vibrational Frequencies for CH4 and Isotopomers

A very accurate ab initio quartic force field for CH4 and its isotopomers is presented. The quartic force field was determined with the singles and doubles coupled-cluster procedure that includes a quasiperturbative estimate of the effects of connected triple excitations, CCSD(T), using the correlation consistent polarized valence triple zeta, cc-pVTZ, basis set. Improved quadratic force constants were evaluated with the correlation consistent polarized valence quadruple zeta, cc-pVQZ, basis set. Fundamental vibrational frequencies are determined using second-order perturbation theory anharmonic analyses. All fundamentals of CH4 and isotopomers for which accurate experimental values exist and for which there is not a large Fermi resonance, are predicted to within +/- 6 cm(exp -1). It is thus concluded that our predictions for the harmonic frequencies and the anharmonic constants are the most accurate estimates available. It is also shown that using cubic and quartic force constants determined with the correlation consistent polarized double zeta, cc-pVDZ, basis set in conjunction with the cc-pVQZ quadratic force constants and equilibrium geometry leads to accurate predictions for the fundamental vibrational frequencies of methane, suggesting that this approach may be a viable alternative for larger molecules. Using CCSD(T), core correlation is found to reduce the CH4 r(e), by 0.0015 A. Our best estimate for r, is 1.0862 +/- 0.0005 A.

Lee, Timothy J.↗