SEARCH · Search NASA

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Gravothermal phase transition, black holes and space dimensionality

In the framework of gravothermal evolution of an ideal monatomic fluid, I examine the dynamical instability of the fluid sphere in (N+1) dimensions by exploiting Chandrasekhar’s criterion to each quasi-static equilibrium along the sequence of the evolution. Once the instability is triggered, it would probably collapse into a black hole if no other interaction halts the process. From this viewpoint, the privilege of (3+1)-dimensional spacetime is manifest, as it is the marginal dimensionality in which the ideal monatomic fluid is stable but not too stable. Moreover, it is the unique dimensionality that allows stable hydrostatic equilibrium with positive cosmological constant. While all higher dimensional (N > 3) spheres are genuinely unstable. In contrast, in (2+1)-dimensional spacetime it is too stable either in the context of Newton’s theory of Gravity or Einstein’s General Relativity. It is well-known that the role of negative cosmological constant is crucial to have the Banados-Teitelboim-Zanelli (BTZ) black hole solution and the equilibrium configurations of a fluid disk. Owing to the negativeness of the cosmological constant, there is no unstable configuration for a homogeneous fluid disk to collapse into a naked singularity, which supports the Cosmic Censorship Conjecture. However, BTZ holes of mass M BTZ > 0 could emerge from collapsing fluid disks. The implications of spacetime dimensionality are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Unraveling Hidden Order and Dynamics in a Heterogeneous Ferroelectric System Using Machine Learning

We uploaded two batches of datasets: Batch-a) Equilibrium dynamics at a constant temperature at zero electric-field, with one trajectory each for 4-different defect structures (i.e. SETs 1 to 4). Each trajectory has a time-step of 0.25 fs, and is run for 7775000 time-steps, with snapshots written out every 4 time-steps (i.e. every 1fs); and Batch-b) Non-Equilibrium dynamics at a constant temperature with the same 0.25 fs time-step, but data dumped every 500 times-steps (i.e. every 125fs) for each of the SETs. The total trajectory of each defect structure (i.e. each SETs 1 to 4) is 2800000 time-steps (5600 snapshots), with stepping of electric-field by 0.01 V/Ã…, after every 100,000 time-steps (i.e. every 200 snapshots), from E=0 to E=0.05 V/Ã… to E = -0.05 V/Ã… to E=0.05 V/Ã….

36 MATERIALS SCIENCE↗

DOE OSTI2024

EF-Hand Battle Royale: Hetero-ion Complexation in Lanmodulin

The lanmodulin (LanM) protein has emerged as an effective means for rare earth element (REE) extraction and separation from complex feedstocks without the use of organic solvents. Whereas the binding of LanM to individual REEs has been well characterized, little is known about the thermodynamics of mixed metal binding complexes (i.e., heterogeneous ion complexes), which limits the ability to accurately predict separation performance for a given metal ion mixture. In this paper, we employ the law of mass action to establish a theory of perfect cooperativity for LanM-REE complexation at the two highest-affinity binding sites. The theory is then used to derive an equation that explains the nonintuitive REE binding behavior of LanM, where separation factors for binary pairs of ions vary widely based on the ratio of ions in the aqueous phase, a phenomenon that is distinct from single-ion-binding chemical chelators. We then experimentally validate this theory and perform the first quantitative characterization of LanM complexation with heterogeneous ion pairs using resin-immobilized LanM. Importantly, the resulting homogeneous and heterogeneous constants enable accurate prediction of the equilibrium state of LanM in the presence of mixtures of up to 10 REEs, confirming that the perfect cooperativity model is an accurate mechanistic description of REE complexation by LanM. We further employ the model to simulate separation performance over a range of homogeneous and heterogeneous binding constants, revealing important insights into how mixed binding differentially impacts REE separations based on the relative positioning of the ion pairs within the lanthanide series. In addition to informing REE separation process optimization, these results provide mathematical and experimental insight into competition dynamics in other ubiquitous and medically relevant, cooperative binding proteins, such as calmodulin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Heating and ionization of stellar chromospheres by nonthermal proton beams: Implications for impulsive phase, redshifted Lyman-alpha radiation in stellar flares

We investigate the physical basis for the timescale of impulsive-phase, redshifted Lyman-alpha emission in stellar flares on the assumption that it is determined by energy losses in a nonthermal proton beam that is penetrating the chromosphere from above. The temporal evolution of ionization and heating in representative model chromospheres subjected to such beams is calculated. The treatment of 'stopping' of beam protons takes into account their interactions with (1) electrons bound in neutral hydrogen, (2) nuclei of neutral hydrogen, (3) free electrons, and (4) ambient thermal protons. We find that, for constant incident beam flux, the system attains an equilibrium with the beam energy input to the chromosphere balanced by radiative losses. In equilibrium, the beam penetration depth is constant, and erosion of the chromosphere ceases. If the redshifted, impulsive-phase stellar flare Lyman-alpha emission is produced by downstreaming hydrogen formed through charge exchange between beam protons and ambient hydrogen, then the emission should end when the beam no longer reaches neutral hydrogen. The durations of representative emission events calculated on this assumption range from 0.1 to 14 s. The stronger the beam, the shorter the timescale over which the redshifted Lyman-alpha emission can be observed.

Brosius, Jeffrey W.↗

DOE OSTI2024

Reductive Dynamic and Static Excited State Quenching of a Homoleptic Ruthenium Complex Bearing Aldehyde Groups

A new homoleptic Ru polypyridyl complex bearing two aldehyde groups on each bipyridine ligand, [Ru(dab) 3 ](PF 6 ) 2 , where dab is 4,4′-dicarbaldehyde-2,2′-bipyridine, was synthesized, characterized, and utilized for iodide photo-oxidation studies. In acetonitrile (CH 3 CN) solution, the complex displayed an intense metal-to-ligand charge transfer (MLCT) absorbance maximum at 475 nm (ε = 22,000 M –1 cm –1 ) and an infrared (IR) band at 1712 cm –1 assigned to the pendent aldehyde groups. Visible light excitation in air-saturated solution resulted in room temperature photoluminescence (PL) with a maximum at 675 nm, a quantum yield, ϕ PL = 0.048, and an excited state lifetime, τ ο = 440 ns, from which radiative and nonradiative relaxation rate constants were extracted, k r = 9.1 × 10 4 s –1 and knr = 1.8 × 10 6 s –1 . Pulsed visible light excitation yielded transient UV–vis and IR absorption spectra consistent with an MLCT excited state; relaxation occurred with the maintenance of two isosbestic points in the visible region, and a lifetime that agreed with that measured by time-resolved PL. Cyclic voltammetry studies in a CH 3 CN solution with 0.1 M TBAPF 6 electrolyte revealed a quasi-reversible oxidation, E°(Ru III/II ) = +1.25 V vs. Fc +/0 , and three sequential one-electron reductions at −1.10, −1.25, and −1.54 V vs. Fc +/0 . Here, an excited state reduction potential of E°(Ru *2+/+ ) = +0.89 V vs. Fc +/0 was estimated with the Rehm–Weller expression. Titration of tetrabutylammonium iodide, TBAI, into a CD 3 CN solution of [Ru(dab) 3 ](PF 6 ) 2 resulted in significant shifts in the aldehyde H atom and 3,3′-biypridyl resonances that were analyzed with a 1:1 equilibrium model, from which K eq = 460 M –1 was extracted, increasing to 5800 M –1 when the solvent was changed to acetone-d 6 . Iodide titrations resulted in a significant quenching of the [Ru(dab) 3 ] *2+ lifetime and quantum yield in both CH 3 CN and acetone solvents. In CH 3 CN, the quenching was mainly dynamic and well described by the Stern–Volmer model, from which a quenching rate constant, k q , of 4.5 × 10 10 M –1 s –1 and an equilibrium constant, K eq , of 8.3 × 10 3 M –1 were obtained. In acetone, the static quenching pathway by iodide was greatly enhanced, with a K eq of 1.2 × 10 4 M –1 and a higher k q of 9.2 × 10 10 M –1 s –1 .

Aldehydes↗

NASA NTRS1993

Direct numerical simulation of turbulent H2-O2 combustion using reduced chemistry

Results of direct numerical simulations of hydrogen-oxygen combustion using a partial-equilibrium chemistry scheme in constant density, decaying, isotropic turbulence are reported. The simulations qualitatively reproduce many features of experimental results, such as superequilibrium radical species mole fractions, with temperature and major species mole fractions closer to chemical equilibrium. It was also observed that the peak reaction rates occur in narrow zones where the stoichiometric surface intersects regions of high scalar dissipation, as might be expected for combustion conditions close to chemical equilibrium. Another finding was that high OH mole fraction correspond more closely to the stoichiometric surface than to areas of high reaction rate for conditions of the simulations. Simulation results were compared to predictions of the Conditional Moment Closure model. This model was found to give good results for all quantities of interest when the conditionally averaged scalar dissipation was used in the prediction. When the nonconditioned average dissipation was used, the predictions compared well to the simulations for most of the species and temperature, but not for the reaction rate. The comparison would be expected to improve for higher Reynolds number flows, however.

Montgomery, Christopher J.↗

NASA NTRS1990

Acoustic positioning and orientation prediction

A method is described for use with an acoustic positioner, which enables a determination of the equilibrium position and orientation which an object assumes in a zero gravity environment, as well as restoring forces and torques of an object in an acoustic standing wave field. An acoustic standing wave field is established in the chamber, and the object is held at several different positions near the expected equilibrium position. While the object is held at each position, the center resonant frequency of the chamber is determined, by noting which frequency results in the greatest pressure of the acoustic field. The object position which results in the lowest center resonant frequency is the equilibrium position. The orientation of a nonspherical object is similarly determined, by holding the object in a plurality of different orientations at its equilibrium position, and noting the center resonant frequency for each orientation. The orientation which results in the lowest center resonant frequency is the equilibrium orientation. Where the acoustic frequency is constant, but the chamber length is variable, the equilibrium position or orientation is that which results in the greatest chamber length at the center resonant frequency.

Barmatz, Martin B.↗

NASA NTRS1986

Thermal stability analysis of the fine structure of solar prominences

The linear thermal stability of a 2D periodic structure (alternatively hot and cold) in a uniform magnetic field is analyzed. The energy equation includes wave heating (assumed proportional to density), radiative cooling and both conduction parallel and orthogonal to magnetic lines. The equilibrium is perturbed at constant gas pressure. With parallel conduction only, it is found to be unstable when the length scale 1// is greater than 45 Mn. In that case, orthogonal conduction becomes important and stabilizes the structure when the length scale is smaller than 5 km. On the other hand, when the length scale is greater than 5 km, the thermal equilibrium is unstable, and the corresponding time scale is about 10,000 s: this result may be compared to observations showing that the lifetime of the fine structure of solar prominences is about one hour; consequently, our computations suggest that the size of the unresolved threads could be of the order of 10 km only.

Demoulin, Pascal↗

NASA NTRS1982

Electric field observations of time constants related to charging and charge neutralization processes in the ionosphere

The Polar 5 electric field results are reviewed, and the transients from Polar 3 are presented. The phenomena are discussed from the standpoint of space charge. On the basis of the Polar 5 results, the large magnitude of the electric field from Polar 3 is seen as indicating that the observed space charge was probably within a few km or less of the payload. Reference is made to Cole's prediction (1960) that charges in the ionosphere would reach equilibrium with a time constant of the order of a few microsec. The processes involved in the two cases presented here require time constants of the order of ms. If the sheath dimensions are taken to be between 50 and 100 m, which is not considered unreasonable in view of the electric field measurements, then a qualitative estimate of the neutralization time would be the transit time for ions across the sheath. Since the kinetic velocity of a 1-eV proton is approximately 14 m/s, it would traverse the distance in 4 to 8 ms, assuming freedom of movement across magnetic field lines. This is the order of the decay times observed on Polar 5.

Maynard, N. C.↗

DOE OSTI2023

A general predictive methodology for fuel-mixture properties up to supercritical conditions

A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME 3–4 ), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME 3–4 . The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.

33 ADVANCED PROPULSION SYSTEMS↗

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↗

DOE OSTI2021

Ultrafast X-ray scattering offers a structural view of excited-state charge transfer

Intramolecular charge transfer and the associated changes in molecular structure in N,N'-dimethylpiperazine are tracked using femtosecond gas-phase X-ray scattering. The molecules are optically excited to the 3p state at 200 nm. Following rapid relaxation to the 3s state, distinct charge-localized and charge-delocalized species related by charge transfer are observed. The experiment determines the molecular structure of the two species, with the redistribution of electron density accounted for by a scattering correction factor. The initially dominant charge-localized state has a weakened carbon–carbon bond and reorients one methyl group compared with the ground state. Subsequent charge transfer to the charge-delocalized state elongates the carbon–carbon bond further, creating an extended 1.634 Å bond, and also reorients the second methyl group. At the same time, the bond lengths between the nitrogen and the ring-carbon atoms contract from an average of 1.505 to 1.465 Å. The experiment determines the overall charge transfer time constant for approaching the equilibrium between charge-localized and charge-delocalized species to 3.0 ps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Initial study on cross-section generation requirements for a PBR equilibrium core

A Serpent model of the equilibrium core HTR-PM small modular nuclear reactor in China, was developed for use in cross-section preparation studies in order to guide methods development for the Griffin reactor multiphysics application. The model includes detailed isotopics for 10 distinct pebble burnup groups in 126 core zones with unique fuel and moderator temperatures obtained from a coupled neutronics-thermal-fluids equilibrium core calculation using Griffin-Pronghorn. A sensitivity study of the fuel and moderator temperatures for various core regions was performed with the MOOSE stochastic tools. The results show that the uncertainties are, not unexpectedly, dominated by the value of the fluid temperature and that the power level, heat transfer coefficient and effective conduction to neighboring pebbles and fluid constitute, at best, second order effects. The temperature uncertainty range varies from 28 K to 57 K between the core entry and exit planes, respectively, but these values are probably higher. We still have to quantify the significance of these uncertainties in the preparation of cross-sections in future work. In addition, we verified that the effective pebble approximation used in the PEBBED and V.S.O.P. computer codes works well for the preparation of region averaged cross-sections. Nevertheless, there are some discrepancies in the cross-sections when compared to the multi-pebble model, which could affect the prediction of peak values and the depletion calculation. We conclude that is highly desirable for future studies with Griffin to be able to handle both the 'effective' pebble approximation and the multi-pebble approach for various pebble burnup groups. This enables Griffin users with the flexibility to perform higher-fidelity studies. Finally, we initiated the preparation of cross-sections for various core regions from the full core Serpent reference model. We quantified the differences in 26 group cross-sections from infinite domain models versus the full core approach. These reference cross-sections will serve to verify the double heterogeneity, self-shielding, and spectrum-correction methods in Griffin. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1982

Sharp edges of planetary rings

The ring systems of Saturn and Uranus exhibit several sharp edges, across which the optical depth drops from order unity to essentially zero. At least two and perhaps all of these features are associated with the location of orbital resonances between a satellite and the ring particles. It is remarkable that the optical depth varies on a distance scale which is much finer than that over which angular momentum can be transferred between a satellite and the ring material. The important features of this phenomenon are: (1) a perturbed band of width Delta a/a of about (satellite mass/planetary mass) to the 1/2 power adjacent to the edge within which the angular momentum transfer occurs, (2) streamlines perturbed such that the angular momentum luminosity decreases smoothly across the band to zero at the edge even though the optical depth remains constant, and (3) dynamical equilibrium requires a relation between the random velocity, the rate of deformation and the optical depth.

Borderies, N.↗

NASA NTRS1990

Comparison of the quadratic configuration interaction and coupled cluster approaches to electron correlation including the effect of triple excitations

The recently proposed quadratic configuration interaction (QCI) method is compared with the more rigorous coupled cluster (CC) approach for a variety of chemical systems. Some of these systems are well represented by a single-determinant reference function and others are not. The finite order singles and doubles correlation energy, the perturbational triples correlation energy, and a recently devised diagnostic for estimating the importance of multireference effects are considered. The spectroscopic constants of CuH, the equilibrium structure of cis-(NO)2 and the binding energies of Be3, Be4, Mg3, and Mg4 were calculated using both approaches. The diagnostic for estimating multireference character clearly demonstrates that the QCI method becomes less satisfactory than the CC approach as non-dynamical correlation becomes more important, in agreement with a perturbational analysis of the two methods and the numerical estimates of the triple excitation energies they yield. The results for CuH show that the differences between the two methods become more apparent as the chemical systems under investigation becomes more multireference in nature and the QCI results consequently become less reliable. Nonetheless, when the system of interest is dominated by a single reference determinant both QCI and CC give very similar results.

Taylor, Peter R.↗

NASA NTRS2015

Mercury-T: A New Code to Study Tidally Evolving Multi-Planet Systems. Applications to Kepler-62

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass 0:1 M(sub ⨀), Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water..

Dynamical evolution and stability↗

NASA NTRS1994

Equilibration temperatures of EL chondrites: A major downward revision in the ferrosilite contents of enstatite

A study of enstatite grains separated from EL6 chondrites shows ferrosilite contents to be approx. 0.02-0.04 mol%, about an order of magnitude lower than reported in the literature. The higher values reported earlier can probably be attributed to the excitation of Fe atoms in Fe-Ni and FeS, in part from micrometer-sized blebs of FeS and Fe-Ni commonly present in enstatite grains but mainly from large neighboring opaque grains excited by bremsstrahlung or by doubly backscattered electrons. Our EL6 ferrosilite contents are similar to those observed in coarse enstatite grains in some aubrites. The equilibrium between ferrosilite and Si-bearing kamacite is highly temperature dependent. At constant Si concentration in kamacite, the equilibrium ferrosilite mole fraction drops by a factor of roughly 2.5 for each 100 K decrease in the equilibrium temperature. As a result, the compositions of these two most common EL6 phases can be used to estimate metamorphic equilibration temperatures. The inferred temperatures are near 1200 K.

Wasson, John T.↗