Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMODYNAMIC EQUILIBRIUM”

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

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

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

Laser-heated diamond anvil cell synthesis and recovery of metastable MnSb 2 and YbZn 2 for post-synthesis transport studies

The creation and exploration of new materials under extreme pressure–temperature conditions has become increasingly reliant on laser-heated diamond anvil cell (LHDAC) techniques, which provide direct access to previously unexplored regions of multinary phase diagrams. Whereas numerous high-pressure phases have been identified in situ, systematic recovery and post-synthesis physical property characterization of these materials remain significant challenges. In this work, we describe the setup and implementation of an LHDAC-based synthesis and recovery workflow and demonstrate its application to metastable MnSb 2 and YbZn 2 phases. Synchrotron x-ray diffraction and spatial mapping confirm dominant formation of the targeted phases, whereas laboratory-based refinement quantifies phase fractions despite intrinsic microstrain and minor secondary phases. High-pressure transport measurements on recovered samples reveal pressure-tunable electronic instabilities in both systems. In MnSb 2 , pressure suppresses two high-temperature magnetic ordering anomalies, observed in transport, by ∼5 GPa and, for higher pressures, induces a new low-temperature feature that increases with further pressure increase. In hexagonal high-pressure YbZn 2 , an electronic reconstruction emerges at ∼11 GPa, characterized by semiconducting-like behavior from ∼30 to 300 K and a broad low-temperature coherence crossover near 30 K. Our results establish LHDAC synthesis not only as a structural discovery tool but also as an experimental platform for investigating correlated quantum states stabilized far from equilibrium thermodynamic conditions.

Huyan, S. [Iowa State University, Ames, IA (United↗
DOE OSTI2021

Multiply-excited states and their contribution to opacity in CO 2 laser-driven tin-plasma conditions

A recent study (2020 Nat. Commun. 11 2334) has found that transitions between multiply-excited configurations in open 4d-subshell tin ions are the dominant contributors to intense EUV emission from dense, Nd:YAG-driven (laser wavelength λ = 1.064 μm) tin plasmas. In the present study, we employ the Los Alamos Atomic code to investigate the spectral contribution from these transitions under industrially-relevant, CO 2 laser-driven (λ = 10.6 μm) tin plasma conditions. First, we employ Busquet's ionisation temperature method to match the average charge state < Z > of a non-local-thermodynamic equilibrium (non-LTE) plasma with an LTE one. This is done by varying the temperature of the LTE calculations until a so-called ionisation temperature T Z is established. Importantly, this approach generates LTE-computed configuration populations in excellent agreement with the non-LTE populations. A corollary of this observation is that the non-LTE populations are well-described by Boltzmann-type exponential distributions having effective temperatures T eff ≈ T Z . In the second part of this work, we perform extensive level-resolved LTE opacity calculations at T Z . It is found that 66% of the opacity in the industrially-relevant 2% bandwidth centred at 13.5 nm arises from transitions between multiply-excited states. These results reinforce the need for the consideration of complex, multiply-excited states in modelling the radiative properties of laser-driven plasma sources of EUV light.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Carrier density oscillation in the photoexcited semiconductor

Abstract The perturbation of a semiconductor from the thermodynamic equilibrium often leads to the display of non-linear dynamics and formation of spatiotemporal patterns due to the spontaneous generation of competing processes. Here, we use scanning ultrafast electron microscopy to show that the transport of hot carriers in the strongly excited semiconductor slows down by turning into an oscillatory process; this is evidenced by the expansions and contractions in the second moment of the distribution. We attribute such a response to the electric field generated by the spatial separation of photo-excited electrons and holes under intrinsic and photo-induced fields; we then introduce a transport model that mimics the experimental observation. Our finding provides a direct imaging evidence for the electrostatic oscillation of hot carriers in the highly excited semiconductor and offers new insights into their dynamics in space and time.

Physics↗
DOE OSTI2020

Non-LTE modelling of cyanoacetylene: evidence for isomer-specific excitation

ABSTRACT Cyanoacetylene molecules are widespread in the interstellar medium (ISM) and several of its isomers have been detected in cold molecular clouds and circumstellar gas. Accurate estimates of the abundance ratio of cyanoacetylene isomers may provide deep insight into their environment. Such knowledge requires rigorous modelling of the emission spectra based on non-local thermodynamic equilibrium (LTE) radiative transfer calculations. To this end, we computed excitation cross-sections of HC2NC and HNC3 induced by collision with para- and ortho-H2, using a quantum mechanical close-coupling method. Then, by thermally averaging these data, we derived rate coefficients for the first 31 low-lying rotational levels of each isomer for temperatures up to 80 K. For the para-H2 collider, the propensity rules are in favour of rotational transitions involving Δj1 = 2 for both isomers, while for the ortho-H2 collider, Δj1 = 2 and Δj1 = 1 rotational transitions are favoured for HC2NC and HNC3, respectively. A comparison of rate coefficients for the HC3N isomers shows differences up to an order of magnitude, especially at low temperatures. Finally, we performed non-LTE radiative transfer calculations to assess the impact of such variations in the analysis of observations. Our simulation suggests that the lack of collisional data specific to each isomer could lead to errors up to a factor of 2–3 in the excitation temperatures. We expect that these data could help in better understanding the cyanoacetylene chemistry and constraining the nitrogen chemistry in the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACT We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Can we observe the ion-neutral drift velocity in prestellar cores?

Given the low-ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer, and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ‘survives’ radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that HCN and HCNH + ⁠, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Outflowing photoionized plasma in Circinus X-1 using the high-resolution X-ray spectrometer Resolve onboard XRISM and the radiative transfer code cloudy

High-resolution X-ray spectroscopy is a key to understanding the mass inflow and outflow of compact objects. Spectral lines carry information about the ionization, density, and velocity structures through their intensity ratios and profiles. They are formed in non-local thermodynamic equilibrium conditions under the intense radiation field from the compact objects; thus radiative transfer (RT) calculation is a requisite for proper interpretations. We present such a study for a low-mass X-ray binary, Circinus X-1, from which the P Cygni profile was discovered using the X-ray grating spectrometer onboard Chandra. We observed the source using the X-ray microcalorimeter onboard XRISM at an orbital phase of 0.93–0.97 and revealed many spectral features unidentified before: the higher series transitions ($n{\rightarrow }1$; $n>2$) of highly ionized (H- and He-like) S, Ca, Ar, and Fe in emission and absorption, the Fe K$\alpha$ and K$\beta$ inner-shell excitation absorption of mildly ionized (O- to Li-like) Fe, and resolved fine-structure level transitions in the Fe Ly$\alpha$ and He$\alpha$ complexes. They blend with each other at different velocity shifts on top of apparently variable continuum emission that changed its flux by an order of magnitude within a 70 ks telescope time. Despite such complexity in the observed spectra, most of them can be explained by a simple model consisting of the photoionized plasma outflowing at $\sim$300 km s$^{-1}$ and the variable blocking material in the line of sight of the incident continuum emission from the accretion disk. We demonstrate this with the aid of the RT code cloudy for the line ratio diagnostics and spectral fitting. We further constrain the physical parameters of the outflow and argue that the outflow is launched close to the outer edge of the accretion disk and can be driven radiatively by being assisted by the line force calculated using the RT simulation.

X-rays: binaries↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2020

Negative capacitance regime in ferroelectrics demystified from nonequilibrium molecular dynamics

Negative capacitance regime in ferroelectrics is associated with the negative slope in the dependence of the polarization on electric field and is prohibited by the laws of equilibrium thermodynamics. Very recently, it developed into a controversial but very productive area with a promise for next generation low-power transistors. At least, in part, the controversies arise from the lack of direct experimental insight into the mechanism for such phenomena. Although the experimental insight is, indeed, challenging at present, we probe the effect in ferroelectric PbTiO 3 bulk and a thin film with realistic first-principles-based nonequilibrium molecular dynamics simulations which model most closely experimental conditions. Here, the following insights emerge from the simulations: (i) the ferroelectric never enters the thermodynamically forbidden region even under extreme rates of electric-field application; (ii) the negative slope naturally arises when the compensating charge is either scarce or slowed down by the external circuitry; (iii) the negative slope is well controlled by the interplay between the applied step pulse field and the amount of the screening charge. The latter can be used to design the negative capacitance regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

High-pressure synthesis and thermodynamic stability of PdH 1±ε up to 8 GPa

Palladium hydride alloys are superconductors and hydrogen storage materials. One synthesis route is compression of Pd to high pressure in a hydrogen-rich environment. Here we report the evolution of the unit cell volume of PdH x synthesized by compressing Pd in a pure H 2 medium to pressures from 0.2 to 8 GPa in a diamond anvil cell at room temperature. The volume of the face-centered cubic unit cell changes nonmonotonically with pressure, increasing upon compression from 0.2 to 1 GPa and decreasing upon compression from 1 to 8 GPa. Volume is reversible upon decompression and is independent of whether the sample was heated to 600 K at low pressure (P < 2 GPa). The x-ray diffraction data show no evidence for a phase transition between 0.2 and 8 GPa. The volume maximum at 1 GPa must be caused by progressive hydrogenation from 0 to 1 GPa. Assuming a pressure-volume-composition equation of state derived from previously published data, the [H]:[Pd] ratio in this study increases to a maximum value of x = 1 ± 0.02 at 2 ± 0.5 GPa and remains stable upon further compression to and from 8 GPa. These results add to a mounting body of evidence that PdH 1±ϵ is in thermodynamic equilibrium with pure H 2 at room temperature from 2 GPa to at least 8 GPa. Here, the simplest interpretation is that H atoms occupy all octahedral sites and no tetrahedral sites in face-centered cubic PdH 1.0 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Laser-induced quenching of metastability at the Mott insulator to metal transition

There is growing interest in strongly correlated insulator thin films because the intricate interplay of their intrinsic and extrinsic state variables causes memristive behavior that might be used for biomimetic devices in the emerging field of neuromorphic computing. Here, in this study, we find that laser irradiation tends to drive V 2 ⁢O 3 from supercooled/superheated metastable states toward thermodynamic equilibrium, most likely in a nonthermal way. We study thin films of the prototypical Mott-insulator V2⁢O3, which show spontaneous phase separation into metal-insulator herringbone domains during the Mott transition. Here, we use low-temperature microscopy to investigate how these metal-insulator domains can be modified by scanning a focused laser beam across the thin film surface. We find that the response depends on the thermal history: When the thin film is heated from below the Mott transition temperature, the laser beam predominantly induces metallic domains. On the contrary, when the thin film is cooled from a temperature above the transition, the laser beam predominantly induces insulating domains. Very likely, the V 2 ⁢O 3 thin film is in a superheated or supercooled state, respectively, during the first-order phase transition, and the perturbation by a laser beam drives these metastable states into stable ones. This way, the thermal history is locally erased. Our findings are supported by a phenomenological model with a laser-induced lowering of the energy barrier between the metastable and equilibrium states.

Materials Science↗
DOE OSTI2023

Stochastic fluctuations in relativistic fluids: Causality, stability, and the information current

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments in the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. Furthermore, the latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Understanding the deficiency in inertial confinement fusion hohlraum x-ray flux predictions using experiments at the National Ignition Facility

The predicted implosion performance of deuterium-tritium fuel capsules in indirect-drive inertial confinement fusion experiments relies on precise calculations of the x-ray drive in laser-heated cavities (hohlraums). This requires accurate, spectrally dependent simulations of laser to x-ray conversion efficiencies and x-ray absorption losses to the hohlraum wall. A set of National Ignition Facility experiments have identified a cause for the long-standing hohlraum “drive deficit” as the overprediction of gold emission at ~2.5 keV in nonlocal thermodynamic equilibrium coronal plasma regions within the hohlraum. Finally, reducing the emission and absorption opacity in this spectral region by ~20% brings simulations into agreement with measured x-ray fluxes and spectra.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermal energy transport in laser-driven high x-ray conversion efficiency metallic silver nanowire foams

Here, we report on record brightness from Ag x-ray emission obtained using a novel laser-produced plasma source. The reported K-shell conversion efficiency of nearly 1% with a radiant energy of ≈ 0.6 kJ/sr from Ag ions is the highest presently recorded and is about twofold greater than more conventional metal-lined cavity targets. He-like Ag ions are the dominant radiators at ≈ 22.7 keV, which contrast the x-ray sources from other pulse power facilities that produce K-shell x-rays from Auger processes in near-neutral ions driven by the nonthermal hot electrons produced in the Z-pinch implosion. The reported x-ray emission was produced from a 4-mm-diameter, 4-mm-long underdense silver nanowire target at 10⁢ mg/⁢cm 3 . The National Ignition Facility laser beams deposited ≈ 1000 kJ of 3ω light into the target with a ≈ 400 TW, 2.5 ns square pulse. Finally, we show via comparisons of data with radiation-hydrodynamics simulations that flux-limited classical thermal energy transport models are inadequate in correctly modeling the behavior of these non-local thermodynamic equilibrium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗