Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal & statistical models”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Energy dependence of light hypernuclei production in heavy-ion collisions from a coalescence and statistical-thermal model perspective

A comparison of light hypernuclei production, from Ultra-relativistic Quantum Molecular Dynamics (UrQMD)+coalescence and the thermal model, in heavy ion collisions over a wide range of beam energies and system sizes is presented. We find that both approaches provide generally similar results, with differences in specific details. Especially the ratios of hypertriton to Λ are affected by both the source radius Δr of the coalescence procedure as well as canonical effects. On the other hand, the double ratio S 3 is almost independent of canonical effects, which is in contrast to coalescence. Thus, both the beam energy dependence and centrality dependence of S 3 can be used to constrain the hypertriton source radius. To do so the currently available data are not yet sufficient. Elliptic flow is shown to be unaffected by the source size of the nuclei and an almost perfect mass scaling of the elliptic flow is observed. Furthermore, our predictions further suggest that the existence of the H-dibaryon(ΛΛ) seems ruled out by ALICE data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Towards absolutely stable ergodicity breaking in two and three dimensions

We propose physically reasonable systems capable of avoiding ergodicity at infinite time in the thermodynamic limit, even with generic perturbations and when coupled to a heat bath. In two dimensions, the rainbow loop soup has (stretched) exponentially numerous absolutely stable nonergodic states with diverging energy but vanishing energy density. In three dimensions the rainbow membrane soup has (stretched) exponentially numerous nonergodic states with diverging energy barriers, leading to infinite-time robust ergodicity breaking that even survives coupling to a nonzero temperature heat bath. We describe our results in the language of exact emergent symmetries and demonstrate how the systems avoid common instabilities. Furthermore, our construction naturally connects to quantum dimer models, topologically ordered systems, the group word construction, and Hamiltonians whose low-energy eigenstates exhibit anomalous entanglement entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

First Measurement of Antideuteron Number Fluctuations at Energies Available at the Large Hadron Collider

The first measurement of event-by-event antideuteron number fluctuations in high energy heavy-ion collisions is presented. The measurements are carried out at midrapidity (|η|<0.8) as a function of collision centrality in Pb-Pb collisions at $\sqrt{^sNN}$=5.02 TeV using the ALICE detector. A significant negative correlation between the produced antiprotons and antideuterons is observed in all collision centralities. The results are compared with a state-of-the-art coalescence calculation. While it describes the ratio of higher order cumulants of the antideuteron multiplicity distribution, it fails to describe quantitatively the magnitude of the correlation between antiproton and antideuteron production. On the other hand, thermal-statistical model calculations describe all the measured observables within uncertainties only for correlation volumes that are different with respect to those describing proton yields and a similar measurement of net-proton number fluctuations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Observation of the Ω⁡(2012) baryon at the LHC

A signal consistent with the Ω⁡(2012) baryon has been observed with a significance of 15⁢𝜎 in pp collisions at $\sqrt{𝑠}$ = 13 TeV at the LHC. In this paper, the analysis technique is described and measurements of the mass and width of the Ω⁡(2012) are reported, along with the first measurement of its transverse-momentum spectrum and yield. This paper corroborates the observation by the Belle Collaboration of this excited Ω state and the observation that the Ω⁡(2012) has a rather narrow width for a strongly decaying resonance. The yield measurement is combined with a statistical thermal model calculation of strange baryon yield ratios to obtain estimates of the Ω⁢(2012) − → $Ξ⁢\bar{\textrm{K}}$ branching ratios. These results will improve our understanding of the internal structure and mass spectrum of excited baryon states and serve as a baseline for searches regarding modifications of these properties in high-temperature media.

Acharya, S. [INFN] (ORCID:0000000292135329)↗
DOE OSTI2026

Multiplicity dependence of f 0 (980) production in pp collisions at $\mathbf{\sqrt{s}}$ = 13 TeV

The dependence of f 0 (980) production on the final-state charged-particle multiplicity is reported for proton–proton (pp) collisions at the centre-of-mass energy, $\mathbf{\sqrt{s}}$ = 13 TeV. The production of (980) is measured with the ALICE detector via the decay channel in a midrapidity region of |y| < 0.5. The evolution of the integrated yields and mean transverse momentum of f 0 (980) as a function of charged-particle multiplicity measured in pp at $\mathbf{\sqrt{s}}$ = 13 TeV follows the trends observed in pp at $\mathbf{\sqrt{s}}$ = 5.02 TeV and in proton–lead (p–Pb) collisions at$\mathbf{\sqrt{s_{NN}}}$ = 5.02 TeV. Particle yield ratios of f 0 (980) to π ± and K* (892) are found to decrease with increasing charged-particle multiplicity. These particle ratios are compared with calculations from the canonical statistical thermal model as a function of charged-particle multiplicity. The thermal model calculations provide a better description of the decreasing trend of particle ratios when no strange or antistrange quark composition for f 0 (980) is assumed, which suggests that the data do not support significant hidden strangeness in the f 0 (980).

Abualrob, I. J. [University of Houston] (ORCID:000↗
DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗
DOE OSTI2025

Hypertriton puzzle in relativistic heavy-ion collisions

The yields of hadrons and light nuclei in relativistic collisions of heavy-nuclei at a center of mass energy of $\sqrt{s_{NN}}$ = 2.6 TeV can be described remarkably well by a thermal distribution of an ideal gas of hadrons and light nuclei interacting only via the decay of resonances. Given the particularly small binding energy of hypertritons relative to the temperature describing the yields (about 156 MeV), one might naturally expect hypertrions to dissociate in medium, making the agreement of hypertriton yields with thermal predictions highly puzzling. The puzzle is compounded by the fact that small binding energy is associated with the large size of the hypertriton. This size is on a similar scale to the overall size of the fireball and much larger than the length scale over which temperatures in the fireball vary over phenomenologically relevant amounts. Here, this paper quantifies the tension this effect causes and shows that it is sufficiently large to render the thermal model inconsistent: its natural assumptions are in conflict with its outputs. The possibility that hypertritons are formed at freeze out as compact objects, quark droplets, that subsequently evolve into hypertritons is considered as a way to resolve the puzzle. It is noted that beyond making the assumption that compact quark droplets form, additional detailed dynamical assumptions which have not been justified are needed to make the thermal model work. The issue of why, despite these issues, the hypertriton is well described by a simple statistical description at freeze out is unresolved. Resolving the hypertriton puzzle is important as it may clarify whether the phenomenological success of the simple thermal model for yields accurately reflects the simple picture of the underlying physics on which it is based.

hydrodynamic models↗
DOE OSTI2022

Coupled fission fragment angular momenta

Here, nuclear fission produces fragments whose spins are coupled to the relative angular motion via angular momentum conservation. It is shown how ensembles of such spins can readily be obtained by either direct microcanonical sampling or by sampling of the associated normal modes of rotation. The resulting distribution of the spin-spin opening angle is illustrated in various three- and two-dimensional scenarios and it is demonstrated how recent mutually conflicting model calculations can be well reproduced with different assumptions about the scission geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

First simultaneous measurement of the γ-ray and neutron emission probabilities in inverse kinematics at a heavy-ion storage ring

The probabilities for γ-ray and particle emission as a function of the excitation energy of a decaying nucleus are valuable observables for constraining the ingredients of the models that describe the deexcitation of nuclei near the particle emission threshold. These models are essential in nuclear astrophysics and applications. In this paper, we have for the first time simultaneously measured the γ-ray and neutron emission probabilities of 208 Pb. The measurement was performed in inverse kinematics at the Experimental Storage Ring (ESR) of the GSI/FAIR facility, where a 208 Pb beam interacted through the 208 Pb (p,p') reaction with a hydrogen gas jet target. Instead of detecting the γ rays and neutrons emitted by 208 Pb, we detected the heavy beamlike residues produced after γ and neutron emission. These heavy residues were fully separated by a dipole magnet of the ESR and were detected with outstanding efficiencies. The comparison of the measured probabilities with model calculations has allowed us to test and select different descriptions of the γ-ray strength function and the nuclear level density available in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Spin cutoff factor and level density for Ni 59 from an analysis of compound nuclear reactions

The spin cutoff parameter for 59 Ni has been studied from different types of experimental data, including neutron angular distributions from the 56 Fe(α,n) 59 Ni reaction, spin of discrete levels, the level density from the proton evaporation spectrum of the 54 Fe( 6 Li,p) 59 Ni reaction and neutron resonance spacing. Experimental data points were compared with calculations using models widely used in literature. It was found that the available empirical models overestimate data points in the energy region below the neutron separation energy, while microscopic calculations which take into account pairing correlations within Hartree-Fock +Bardeen-Cooper-Schrieffer approach are consistent with data. As a result, it confirmed earlier findings that pairing correlations play an important role and need to be taken into account when the spin cutoff parameter is calculated below the neutron separation energy.

(a,n) reaction↗
DOE OSTI2022

β — -delayed fission in the coupled quasiparticle random-phase approximation plus Hauser-Feshbach approach

Beta-delayed neutron emission and β-delayed fission (βdf) probabilities were calculated for heavy, neutronrich nuclei using the Los Alamos coupled quasiparticle random-phase approximation plus Hauser-Feshbach (QRPA + HF) approach. In this model, the compound nucleus is initially populated by β decay and is followed through subsequent statistical decays taking into account competition among neutrons, γ rays, and fission. The primary output of these calculations includes branching ratios along with neutron and γ -ray spectra. Here, we find a relatively large region of heavy nuclides where the probability of βdf is near 100%. For a subset of nuclei near the neutron drip line, delayed neutron emission and the probability of fission are both large, which leads to the possibility of multichance βdf (mc-βdf). We comment on prospective neutron-rich nuclei that could be probed by future experimental campaigns and provide a full table of branching ratios in ASCII format in the Supplemental Material for use in various applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Finite density QCD equation of state: Critical point and lattice-based 𝑇′ expansion

Here, we present a novel construction of the QCD equation of state (EoS) at finite baryon density. Our work combines a recently proposed resummation scheme for lattice QCD results with the universal critical behavior at the QCD critical point. This allows us to obtain a family of equations of state in the range 0 ≤ 𝜇 𝐵 ≤ 700 MeV and 25 MeV ≤ 𝑇 ≤ 800 MeV, which match lattice QCD results near 𝜇 𝐵 =0 while featuring a critical point in the 3D Ising model universality class. The position of the critical point can be chosen within the range accessible to beam-energy scan heavy-ion collision experiments. The strength of the singularity and the shape of the critical region are parametrized using a standard parameter set. We impose stability and causality constraints and discuss the available ranges of critical point parameter choices, finding that they extend beyond earlier parametric QCD EoS proposals. We present thermodynamic observables, including baryon density, pressure, entropy density, energy density, baryon susceptibility and speed of sound, that cover a wide range in the QCD phase diagram relevant for experimental exploration.

Astronomy & Astrophysics↗
DOE OSTI2023

Entropy and Seebeck signals meet on the edges

Here, we explore the electronic entropy per particle s and Seebeck coefficient S in zigzag graphene ribbons. Pristine and edge-doped ribbons are considered using tight-binding models to inspect the role of edge states in the observed thermal transport properties. As a band gap opens when the ribbons are doped at one or both edges, due to asymmetric edge potentials, we find that s and S signals are closely related to each other: both develop sharp dip-peak line shapes as the chemical potential lies in the gap, while the ratio s/S exhibits a near-constant value equal to the elementary charge e at low temperatures. This constant ratio suggests that S can be seen as the transport differential entropy per charge, as suggested by some authors. Our calculations also indicate that measurement of s and S may be useful as a spectroscopic probe of different electronic energy scales involved in such quantities in gapped materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Low-energy enhancement of the magnetic dipole radiation in odd-mass lanthanides

We compute the magnetic dipole (M1) $\gamma$-ray strength functions ($\gamma$SF) for the odd-mass lanthanides $^{\textrm{143-151}}$Nd and $^{\textrm{147-153}}$Sm using the shell-model Monte Carlo method in combination with the static-path approximation and the maximum-entropy method. In particular, we quantify the statistical uncertainties in the calculated M1 $\gamma$SFs and show that they are under control for the excitation energies relevant to the experiments despite a Monte Carlo sign problem that originates in the projection onto an odd number of neutrons. We identify a low-energy enhancement (LEE) in the M1 $\gamma$SFs of these odd-mass lanthanides, which was recently observed experimentally in some of them. We also find a scissors mode resonance (SR) in the strongly deformed isotopes. We observe that the decrease in the LEE strength with neutron number along an isotopic chain is compensated for by an increase in the SR strength in the deformed nuclei. Furthermore, we compare our results with recent experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Measurement of fifth- and sixth-order fluctuations of (net-)proton number in Au + Au collisions from phase II of the beam energy scan program at RHIC

We report high-statistics measurements of fifth- and sixth-order factorial cumulants and cumulant ratios of (net-)proton multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7–27GeV, using data from the STAR experiment collected during the Beam Energy Scan Phase II at RHIC. Protons and antiprotons are identified at midrapidity (|𝑦| < 0.5) with transverse momentum 0.4 < 𝑝 𝑇 < 2.0GeV/𝑐. The proton factorial cumulants 𝜅 4 , 𝜅 5 , and 𝜅 6 increase with order but exhibit no sign alternation within current uncertainties, offering no evidence for a two-component structure in the proton multiplicity distribution, as might be expected near a first-order phase transition. The cumulant ratios 𝐶 5 /𝐶 1 and 𝐶 6 /𝐶 2 fluctuate around zero in collisions at 0–40% centrality. The results are consistent with both the negative predictions from lattice QCD and the positive trends obtained from the ultrarelativistic quantum molecular dynamics (UrQMD) model. At $\sqrt{s_{NN}}$ ≳ 27GeV, the 𝐶 4 /𝐶 2 and 𝐶 5 /𝐶 1 results are compatible with predictions from lattice QCD, functional renormalization group (FRG), and hadron resonance gas (HRG) models, while UrQMD describes the data better at lower energies. Here, these measurements place constraints on baryon number fluctuations and offer valuable insights into the QCD phase structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

First 𝛽-Delayed Two-Neutron Spectroscopy of the 𝑟-Process Nucleus 134 In and Observation of the 𝑖 13/2 Single-Particle Neutron State in 133 Sn

This manuscript reports on the direct observation of a 𝛽-delayed two-neutron emission in a study of 134 In at the ISOLDE Decay Station using neutron spectroscopy. We also report on the first measurement in 𝛽 − decay of the long-sought 13/2 + excited state in 133 Sn, attributed to be the neutron single-particle 𝑖 13/2 orbital. The observation of sequential neutron emission is used to extract the relative population of the 𝑖 13/2 state, which was found to be much smaller than the predictions of the statistical model. The experiment was possible because of the innovative use of a neutron array with neutron discrimination and interaction tracking capabilities. This is the first study of the details of the two-neutron emission for a nucleus, which belongs to the 𝑟-process path. Understanding 𝛽-delayed two-neutron emission probabilities is essential to validate models used in astrophysical 𝑟-process nucleosynthesis calculations. Observing two-neutron emissions in 𝛽 − decay paves the way for new experiments to study energy and angular correlations for 𝛽-delayed multineutron emitters.

Beta decay↗
DOE OSTI2025

Deviations from the Porter-Thomas Distribution due to Nonstatistical 𝛾 Decay below the 150 Nd Neutron Separation Threshold

We introduce a new method for the study of fluctuations of partial transition widths based on nuclear resonance fluorescence experiments with quasimonochromatic linearly polarized photon beams below particle separation thresholds. It is based on the average branching of decays of 𝐽=1 states of an even-even nucleus to the 2$^{+}_{1}$ state in comparison to the ground state. Between 5 and 7 MeV, a constant average branching ratio for 𝛾 decays from 1 − states of 0.490(16) is observed for the nuclide 150 Nd. Assuming 𝜒 2 -distributed partial transition widths, this average branching ratio is related to a degree of freedom of 𝜈 = 1.93⁢(12), rejecting the validity of the Porter-Thomas distribution, requiring 𝜈 = 1. The observed deviation can be explained by nonstatistical effects in the 𝛾-decay behavior with contributions in the range of 9.4(10)% up to 94(10)%.

150 ≤ A ≤ 189↗