Search NASA⌕ Search

SEARCH · Search NASA

Results for “local 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Comparison of excitation temperature of a laser-produced plasma by combining emission and absorption spectroscopy

Measurement of the temporal evolution of laser-produced plasma temperature is very important for many of its applications, and several plasma diagnostic tools are routinely used by researchers. However, it is very challenging to measure the properties of the plasma at the early and late times of its evolution using a single diagnostic tool. In this study, we combined emission and laser absorption spectroscopy to compare the excitation temperatures of a laser-produced uranium plasma system. Several U I transitions in the near-infrared spectral range (775–800 nm) were considered, and the Boltzmann plot method was used to measure the excitation temperatures using both emission and absorption spectroscopy. Emission spectroscopy provided early-time temperature measurements of the plasma up to times 2–20 µs, while absorption spectroscopy provided temperature measurements at late times of plasma evolution (for times 5–80 µs). The emission and absorbance of U I transitions were found to follow the Boltzmann distribution, indicating the plasma is likely in the state of local thermodynamic equilibrium even at late times of its lifetime. The emission and absorption-based time-resolved excitation temperatures demonstrated good agreement at earlier times (≤15 µs) in the overlapped temporal region, while a deviation in the measured values was seen at times (≥15 µs), and potential reasons for such a disagreement are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Key advancements toward eliminating the “drive deficit” in ICF hohlraum simulations

This paper, extended from a Letter [Chen et al., Phys. Rev. E 110, L013201 (2024)] published in Phys. Rev. E in 2024, presents detailed results from a set of experiments studying the x-ray drive produced in laser-heated ignition-scale hohlraums at the National Ignition Facility. The results reveal that deficiencies in non-local thermodynamic equilibrium atomics physics calculations of gold emissivity lead to the overestimation of the “M-band” radiation and are responsible for discrepancies in drive predictions vs measurements in these hohlraums. These experiments were performed on ViewFactor targets with the same gas fill and laser pulse shapes as currently used in ongoing ignition experiments. The common hohlraum model with a magnetohydrodynamics heat transport model and a 10%–20% reduction to the M-band absorption and emission opacity can bring the simulations into agreement with several independent measurements, including the absolute radiation flux, the radiation flux in the 2–4 keV spectral range, and spectroscopic measurements of both K-shell dopant and L-shell gold. This new model is applicable to ICF and high-energy density experiments that use gold hohlraums as drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resolving discrepancies in bang-time predictions for indirect-drive ICF experiments on the NIF: Insights from the Build-A-Hohlraum campaign

This study investigated discrepancies between measured and simulated x-ray drive in Indirect-Drive Inertial Confinement Fusion (ID-ICF) hohlraums at the National Ignition Facility. Despite advances in radiation-hydrodynamic simulations, a consistent “drive deficit” remains. Experimentally measured ID-ICF capsule bang-times are systematically 400–700 ps later than simulations predict. The Build-A-Hohlraum (BAH) campaign explored potential causes for this discrepancy by systematically varying hohlraum features, including laser entrance hole (LEH) windows, capsules, and gas fills. Overall, the agreement between simulated and experimental x-ray drive was found to be largely unaffected by these changes. The data allow us to exclude some hypotheses put forward to potentially explain the discrepancy. Errors in the local thermodynamic equilibrium (LTE) atomic modeling, errors in the modeling of LEH closure, and errors due to a lack of plasma species mix physics in simulations are shown to be inconsistent with our measurements. Instead, the data support the hypothesis that errors in NLTE emission modeling are a significant contributor to the discrepancy. X-ray emission in the 2–4 keV range is found to be approximately 30% lower than in simulations. This is accompanied by higher than predicted electron temperatures in the gold bubble region, pointing to errors in non-LTE modeling. Introducing an opacity multiplier of 0.87 on energy groups above 1.8 keV improves agreement with experimental data, reducing the bang-time discrepancy from 300 to 100 ps. These results underscore the need for refined NLTE opacity models to enhance the predictive power of hohlraum simulations.

Band emission↗
DOE OSTI2025

X-ray spectroscopy of multi-temperature plasmas using the differential emission measure formalism

We present a theoretical construct that nominally underlies spectroscopic data analysis of multi-temperature plasmas, known as the differential emission measure (DEM). From a data analytic perspective, the DEM formalism is used to derive temperature distributions from line spectra that are formed in the presence of temperature gradients and by time integrations of evolving plasmas. From a modeling perspective, DEMs are convenient intermediaries between radiation hydrodynamics simulations and spectroscopic measurements acquired in the laboratory. The DEM concept and its associated methodologies were originally developed by spectroscopists working with astrophysical data. We borrow from these earlier investigations. In this manuscript, intended primarily as a tutorial, we discuss the basic concepts, but also augment various aspects of the theory by the way of extension and example, including a detailed treatment of various weighting and averaging schemes, intended to mitigate ambiguities that often arise when reporting temperature information. We focus on high-temperature plasmas that are not in local thermodynamic equilibrium and the x-ray spectra that they produce, although the core ideas presented here are applicable to spectroscopy in other energy bands. A few examples involving the derivation and manipulation of model DEMs in simple geometries are provided.

Liedahl, Duane A. [Lawrence Livermore National Lab↗
DOE OSTI2020

Soft X-ray and EUV emission spectra of beryllium plasma produced by neodymium-glass laser radiation with broad frequency and angular spectra

We present the results of an experimental study of soft X-ray (SXR) and extreme ultraviolet emission spectra of the plasma produced by exposing a plane solid beryllium target to laser radiation with broad frequency and angular spectra. SXR lines up to 1s – 9p of Be IV as well as the plasma continuum are recorded for a laser focal-spot intensity of 5.3 × 10{sup 13} W cm{sup −2}. To model the SXR beryllium plasma spectra, simulations are carried out using the INDHAUS programme and the FLYCHK code in the framework of local thermodynamic equilibrium model, which agree nicely with experimentally obtained results. (laser plasma)

36 MATERIALS SCIENCE↗
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 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 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 OSTI2024

Actinide signatures in low electron fraction kilonova ejecta

ABSTRACT Neutron star (NS) mergers are known to produce heavy elements through rapid neutron capture (r-process) nucleosynthesis. Actinides are expected to be created solely by the r-process in the most neutron-rich environments. Confirming if NS mergers provide the requisite conditions for actinide creation is therefore central to determining their origin in the Universe. Actinide signatures in kilonova (KN) spectra may yield an answer, provided adequate models are available in order to interpret observational data. In this study, we investigate actinide signatures in neutron-rich merger ejecta. We use three ejecta models with different compositions and radioactive power, generated by nucleosynthesis calculations using the same initial electron fraction ($Y_e = 0.15$) but with different nuclear physics inputs and thermodynamic expansion history. These are evolved from 10 to 100 d after merger using the sumo non-local thermodynamic equilibrium (NLTE) radiative transfer code. We highlight how uncertainties in nuclear properties, as well as choices in thermodynamic trajectory, may yield entirely different outputs for equal values of $Y_e$. We consider an actinide-free model and two actinide-rich models, and find that the emergent spectra and light-curve evolution are significantly different depending on the amount of actinides present, and the overall decay properties of the models. We also present potential key actinide spectral signatures, of which doubly ionized $_{89}$Ac and $_{90}$Th may be particularly interesting as spectral indicators of actinide presence in KN ejecta.

Pognan, Quentin (ORCID:0000000152550782)↗
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 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 OSTI2020

Impact of thermodynamic fluctuations and pattern size on the nucleation behavior during area selective deposition

In this work, the authors explore the impact that thermodynamic fluctuations have on the spontaneous appearance of defects during area selective deposition on patterned surfaces that are fully covered by adsorbates under local thermodynamic equilibrium, such as self-assembled monolayers. By using a simple lattice gas model for the adsorbed monolayer, the authors were able to track the spontaneous formation of defects in the monolayer as a function of the pattern width. The results indicate that, for pattern widths of the order of tens of nanometers, roughening effects at the pattern edge can be the leading source for the spontaneous appearance of nucleation defects. This leads to an enhancement of the density of defects that can be up to three orders of magnitude higher than those expected in uniform (not patterned) surfaces. The model also predicts a density of defects that is inversely proportional to the pattern width. Finally, if the dynamic nucleation of defects during area selective deposition is driven by thermodynamic fluctuations, the model predicts that the nucleation rate should be proportional to the total precursor fluency and independent of purge times. Moreover, a tight confinement of the monolayer through a high quality smooth interface in the patterned substrate and strong cohesive interactions between adsorbates should each contribute to a reduction of the overall defect density.

36 MATERIALS SCIENCE↗
DOE OSTI2024

$$\Lambda $$ polarization in very high energy heavy ion collisions as a probe of the quark–gluon plasma formation and properties

Abstract We have studied the spin polarization of $$\Lambda $$ Λ hyperons in heavy ion collisions at center-of-mass energies $$\sqrt{s_\textrm{NN}}= 200$$ s NN = 200 GeV and $$\sqrt{s_\textrm{NN}}= 5.02$$ s NN = 5.02 TeV carried out at RHIC and LHC colliders. We have calculated the mean spin vector at local thermodynamic equilibrium, including all known first-order terms in the gradients of the thermo-hydrodynamic fields, assuming that the hadronization hypersurface has a uniform temperature. We have also included the feed-down contributions to the polarization of $$\Lambda $$ Λ stemming from the decays of polarized $$\Sigma ^*$$ Σ ∗ and $$\Sigma ^0$$ Σ 0 hyperons. The obtained results are in good agreement with the data. In general, the component of the spin vector along the global angular momentum, orthogonal to the reaction plane, shows strong sensitivity to the initial longitudinal flow velocity. Furthermore, the longitudinal component of the spin vector turns out to be very sensitive to the bulk viscosity of the plasma at the highest LHC energy. Therefore, the azimuthal dependence of spin polarization can effectively constrain the initial hydrodynamic conditions and the transport coefficients of the quark gluon plasma.

Physics↗
DOE OSTI2024

SuperNu Version 4.x

We seek to release SuperNu, Version 4.x, as a continuation of development for the open source SuperNu software. The SuperNu, Version 3.x Monte Carlo radiative transfer code for astrophysical transients was released with GPLv3 copyright, asserted by LANL in 2015. For the next release we have features planned for development , including: opacity implementation (including non-local thermodynamic equilibrium effects), generalized source implementation (e.g. for emulating shock heating in Type II supernovae), 3T (electron, ion, radiation) internal energy update, special relativity corrections through O(v^2/c^2), light polarization (e.g. for comparison to spectroplarimetry observations of supernovae and kilonovae), and infrastructure features (checkpoint and restart of simulations, and tools including setup and Slurm batch scripts and simulation post-processing/analysis scripts). These features are intended to improve the fidelity and/or better understand uncertainty in supernova and kilonova light curve calculations.

Wollaeger, Ryan↗