Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite temperature”

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

Fast and flexible neutrino decoupling. Part I. The Standard Model

Cosmological determinations of the number of relativistic neutrino species, N eff , are becoming increasingly accurate, and further improvements are expected both from CMB and BBN data. Given this context, we update the evaluation of N eff and the current entropy density via the momentum-averaged approach. This allows for a numerically fast description of neutrino decoupling, easily portable to an array of new physics scenarios. We revisit all aspects of this approach, including collision terms with full electron mass dependence, finite temperature QED corrections to the equation of state, neutrino oscillations, and the modelling of neutrino ensembles with effective chemical potentials. For integrated observables, our results differ by less than 0.04% from the solution of the momentum-dependent evolution equation. We outline how to extend the approach to BSM settings, and will highlight its power in Part II. To facilitate the practical implementation, we release a Mathematica and python code within nudec_BSM_v2, easily linkable to BBN codes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Toward machine learning interatomic potentials for modeling uranium mononitride

Uranium mononitride (UN) is a promising accident-tolerant fuel because of its high fissile density and high thermal conductivity. In this study, we developed the first machine learning interatomic potentials for reliable atomic-scale modeling of UN at finite temperatures. We constructed a training set using density functional theory (DFT) calculations that was enriched through an active learning procedure, and two neural network potentials were generated. Both potentials successfully reproduce key thermophysical properties of interest, such as temperature-dependent lattice parameter, specific heat capacity, and bulk modulus. We also evaluated the energy of stoichiometric defect reactions and defect migration barriers and found close agreement with DFT predictions, demonstrating that our potentials can be used for modeling defects in UN. Additional tests provide evidence that our potentials are reliable for simulating diffusion, noble gas impurities, and radiation damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Numerical relativity simulations of the neutron star merger GW190425: microphysics and mass ratio effects

ABSTRACT GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

r -process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants

ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magnetohydrodynamic stability of magnetars in the ultrastrong field regime – II. The crust

ABSTRACT We study the stability of Hall magnetohydrodynamic with strong magnetic fields in which Landau quantization of electrons is important. We find that the strong-field Hall modes can be destabilized by the dependence of the differential magnetic susceptibility on magnetic field strength. This hydrodynamic instability, thermodynamic in origin and stabilized by magnetic domain formation, is studied using linear perturbation theory. It is found to have typical growth time of order ≲103 yr, with the growth time decreasing as a function of wavelength of the perturbation. The instability is self-limiting, turning off following a period of local field growth by a few per cent of the initial value. Finite temperature is also shown to limit the instability, with sufficiently high temperatures eliminating it altogether. Alfvén waves can show similar unstable behaviour on shorter time-scales. We find that Ohmic heating due to the large fields developed via the instability and magnetic domain formation is not large enough to account for observed magnetar surface temperatures. However, Ohmic heating is enhanced by the oscillatory differential magnetic susceptibility of Landau-quantized electrons, which could be important to magnetothermal simulations of neutron star crusts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Theoretical description of atomtronic Josephson junctions in an optical lattice

Experimental realizations of “atomtronic” Josephson junctions have recently been created in annular traps in relative rotation with respect to potential barriers that generate the weak links. If these devices are additionally subjected to an optical lattice potential, then they can incorporate strong-coupling Mott physics within the design, which can modify the behavior and can allow for interesting new configurations of barriers and of superfluid flow patterns. We examine theoretically the behavior of a Bose superfluid in an optical lattice in the presence of an annular trap and a barrier across the annular region which acts as a Josephson junction. As the superfluid is rotated, circulating super-currents appear. Beyond a threshold superfluid velocity, phase slips develop, which generate vortices. We use a finite temperature strong-coupling expansion about the mean-field solution of the Bose Hubbard model to calculate various properties of such devices. Additionally, we discuss some of the rich behavior that can result when there are Mott regions within the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Breakdown of superdiffusion in perturbed quantum integrable spin chains and ladders

Superdiffusive transport with dynamical exponent z=3/2 has been firmly established at finite temperature for a class of integrable systems with a non-Abelian global symmetry G. On the inclusion of integrability-breaking perturbations, diffusive transport with z=2 is generically expected to hold in the limit of late time. Recent studies of the classical Haldane-Ishimori-Skylanin model have found that perturbations that preserve the global symmetry lead to a much slower time scale for the onset of diffusion, albeit with uncertainty over the exact scaling exponent. That is, for perturbations of strength λ, the characteristic timescale for diffusion goes as t*∼λ−α for some α. Using large-scale matrix product state simulations, we investigate this behavior for perturbations to the canonical quantum model showing superdiffusion: the S=1/2 quantum Heisenberg chain. We consider a ladder configuration and look at various perturbations that either break or preserve the SU(2) symmetry, leading to scaling exponents consistent with those observed in one classical study [McCarthy , Phys. Rev. B 110, L180301 (2024)2469-995010.1103/PhysRevB.110.L180301]: α=2 for symmetry-breaking terms and α=6 for symmetry-preserving terms. We also consider perturbations from another integrable point of the ladder model with G=SU(4) and find consistent results. Finally, we consider a generalization to an SU(3) ladder and find that the α=6 scaling appears to be universal across superdiffusive systems when the perturbations preserve the non-Abelian symmetry G.

Wang, Kevin↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2026

Krylov winding and emergent coherence in operator growth dynamics

The operator wavefunction provides a fine-grained description of quantum chaos and of the irreversible growth of simple operators into increasingly complex ones. Remarkably, at finite temperature this wavefunction can acquire a phase that increases linearly with the operator’s size, a phenomenon called . Although size winding occurs naturally in a holographic setting, the emergence of a coherent phase in a scrambled operator remains mysterious from the standpoint of a thermalizing quantum many-body system. Here, in this article, we elucidate this phenomenon by introducing the related concept of , whereby the operator wavefunction acquires a phase which winds linearly with the Krylov index. We show that Krylov winding is a generic feature of quantum chaotic systems and is a direct consequence of the universal operator growth bound hypothesis. It gives rise to size winding under two additional conditions: (i) a low-rank mapping between the Krylov and size bases, which ensures phase alignment among operators of the same size, and (ii) the saturation of the "chaos-operator growth" bound 𝜆 𝐿 ≤ 2⁢𝛼 (with 𝜆 𝐿 the Lyapunov exponent and 𝛼 the growth rate), which ensures a linear phase dependence on size. For systems which do not saturate this bound, with ℎ = 𝜆 𝐿 /2⁢𝛼 < 1, the winding with Pauli size ℓ becomes superliner, behaving as ℓ 1/ℎ . We illustrate these results with two classes of microscopic models: the Sachdev-Ye-Kitaev (SYK) model and its variants, and a disordered 𝑘-local spin model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Delocalized Excitation Transfer in Open Quantum Systems with Long-Range Interactions

The interplay between coherence and system-environment interactions is at the basis of a wide range of phenomena, from quantum information processing to charge and energy transfer in molecular systems, biomolecules, and photochemical materials. In this work, we use a Frenkel exciton model with long-range interacting qubits coupled to a damped collective bosonic mode to investigate vibrationally assisted transfer processes in donor-acceptor systems featuring internal substructures analogous to light-harvesting complexes. We find that certain delocalized excitonic states maximize the transfer rate and that the entanglement is preserved during the dissipative transfer over a wide range of parameters. We investigate the reduction in transfer caused by static disorder, white noise, and finite temperature and study how transfer efficiency scales as a function of the number of dimerized monomers and the component number of each monomer, finding which excitonic states lead to optimal transfer. Finally, we provide a realistic experimental setting to realize this model in analog trapped-ion quantum simulators. Analog quantum simulation of systems comprising many and increasingly complex monomers could offer valuable insights into the design of light-harvesting materials, particularly in the nonperturbative intermediate parameter regime examined in this study, where classical simulation methods are resource intensive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
DOE OSTI2025

Monomer-dimer tensor-network basis for qubit-regularized lattice gauge theories

Traditional SU⁡(𝑁) lattice gauge theories (LGTs) can be formulated using an orthonormal basis constructed from the irreducible representations (irreps) 𝑉 𝜆 of the SU⁡(𝑁) gauge symmetry. On a lattice, the elements of this basis are tensor networks comprising dimer tensors on the links labeled by a set of irreps {𝜆 ℓ } and monomer tensors on sites labeled by {𝜆 𝑠 }. These tensors naturally define a local site Hilbert space, ℋ$^𝑔_𝑠$, on which gauge transformations act. Gauss’s law introduces an additional index 𝛼 𝑠 =1,2,…,𝒟⁡(ℋ$^𝑔_𝑠$) that labels an orthonormal basis of the gauge-invariant subspace of ℋ$^𝑔_𝑠$. This monomer-dimer tensor-network (MDTN) basis, |{𝜆 𝑠 },{𝜆 ℓ },{𝛼 𝑠 }⟩, of the physical Hilbert space enables the construction of new qubit-regularized SU⁡(𝑁) gauge theories that are free of sign problems while preserving key features of traditional LGTs. Here, we investigate finite-temperature confinement-deconfinement transitions in a simple qubit-regularized SU(2) and SU(3) gauge theory in 𝑑 =2 and 𝑑 =3 spatial dimensions, formulated using the MDTN basis, and show that they reproduce the universal results of traditional LGTs at these transitions. Additionally, in 𝑑 =1, we demonstrate using a plaquette chain that the string tension at zero temperature can be continuously tuned to zero by adjusting a model parameter that plays the role of the gauge coupling in traditional LGTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Relaxation time approximation for a multispecies relativistic gas

We generalize a recent prescription for the relaxation time approximation for the relativistic Boltzmann equation for systems with multiple particle species at finite temperature. This is performed by adding counter-terms to the traditional Anderson-Witting ansatz for each particle species. Our approach allows for the use of momentum-dependent relaxation times and the obedience of local conservation laws regardless of the definition of the local equilibrium state. As an application, we derive the first order Chapman-Enskog corrections to the equilibrium distribution and display results for the hadron-resonance gas. We also demonstrate that our collision term ansatz obeys the second law of thermodynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗
DOE OSTI2025

Prototypes of Nonrelativistic Spin Splitting and Polarization in Symmetry Broken Antiferromagnets

Antiferromagnets that break both space-time reversal and translation-spin-rotation symmetries were recently predicted [L.-D. Yuan, Z. Wang, J.-W. Luo, E. I. Rashba, and A. Zunger, Phys. Rev. B 102, 014422 (2020)] to possess splitting between the otherwise spin-degenerate energy bands even without the relativistic spin-orbit coupling (SOC). Here, we point out that such nonrelativistic spin splitting (NRSS)—in particular, “spin splitting type 4” (SST-4) symmetry-broken antiferromagnets—can be divided into subgroups having distinct patterns of spin splitting and spin textures, depending on additional auxiliary symmetries of spin interconversion and polarity. These SST-4 subgroups include the 𝛼-type (no spin-interconverting symmetry) having spin splitting at the Brillouin zone center, as well as the 𝛽 subgroup in which a rotation symmetry is applied and determines the alternating spin texture and the 𝛾 subgroup having exclusively reflection spin-interconverting symmetry. Unlike ferrimagnets, the 𝛼-type compounds are shown to have tiny net magnetization at finite temperature and thus avoid the adverse effect of the stray field. The 𝛼 and 𝛽 subgroups can be either polar or nonpolar, whereas the 𝛾 subgroup is polar only, providing a basis for possible switching by external fields. The combination of NRSS-enabling and auxiliary symmetries is used here as a filter for identifying previously synthesized compounds as specific prototypes. Their characteristic splitting and spin polarization are calculated by density functional theory to the benefit of potential future experimental testing. Interesting results are as follows: (i) SOC-independent NRSS can exceed the magnitude of the SOC-induced Rashba and Dresselhaus spin splitting in semiconductors. (ii) Examples of predicted 𝛼-type insulating compounds include BiCrO 3 (nonpolar) and Mn 2 ⁢ScSbO 6 (polar), the latter having spin splitting of 158 meV and 160 meV in the valence and conduction bands, respectively. (iii) The 𝛽-type (Cu 2 ⁢Y 2 ⁢O 5 and FeF 2 ) and 𝛾-type compounds (Mn 4 ⁢Nb 2 ⁢O 9 and FeScO 3 ) are distinguished both by their auxiliary symmetries and polarity. The spin textures of 𝛾-type compounds are mirror reflected with spin degeneracy of the wave vectors on that mirror. These observations will likely broaden the experimental playing field of NRSS physics significantly.

antiferromagnets↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2025

Distinguishing Thermal Fluctuations from Polaron Formation in Halide Perovskites

Recent angle-resolved photoelectron spectroscopy (ARPES) measurements of the hole effective mass in CsPbBr$_3$ revealed an enhancement of $\sim$50 % compared to the bare mass computed from first principles for CsPbBr$_3$ at $T = 0 K$. This large enhancement was interpreted as evidence of polaron formation. Employing accurate finite-temperature first-principles calculations, we show that the calculated hole effective mass of CsPbBr$_3$ at $T = 300 K$ can explain experimental results without invoking polarons. Thermal fluctuations are particularly strong in halide perovskites compared to conventional semiconductors such as Si and GaAs, and cannot be ignored when comparing with experiment. We not only resolve the debate on polaron formation in halide perovskites, but also demonstrate the general importance of including thermal fluctuations in first-principles calculations for strongly anharmonic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Universal thermodynamics of an SU ( N ) Fermi-Hubbard model

The SU(2) symmetric Fermi-Hubbard model (FHM) plays an essential role in strongly correlated fermionic many-body systems. In the one particle per site and strongly interacting limit U/t $\gg$ 1, it is effectively described by the Heisenberg Hamiltonian. In this limit, enlarging the spin and extending the typical SU(2) symmetry to SU(N) has been predicted to give exotic phases of matter in the ground state, with a complicated dependence on N. This raises the question of what—if any—are the finite-temperature signatures of these phases, especially in the currently experimentally relevant regime near or above the superexchange energy. We explore this question for thermodynamic observables by numerically calculating the thermodynamics of the SU(N) FHM in the two-dimensional square lattice near densities of one particle per site, using determinant quantum Monte Carlo and numerical linked cluster expansion. Interestingly, we find that for temperatures above the superexchange energy, where the correlation length is short, the energy, number of on-site pairs, and kinetic energy are universal functions of N. Although the physics in the regime studied is well beyond what can be captured by low-order high-temperature series, we show that an analytic description of the scaling is possible in terms of only one- and two-site calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗