SEARCH · Search NASA
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.
The order parameter profile between coexisting confined and plasma regions at the quantum chromodynamic (QCD) phase transition is constructed. The dimensionless combination of the surface energy (Sigma) and the correlation length (Zeta) is estimated to be Sigma Zeta 3 approximately equals 0.8.
We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the nrpmw waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f$$^{peak}_{2}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f$$^{peak}_{2}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f$$^{peak}_{2}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed 1.6σ. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. Furthermore, unless the phase transition is “strong,” disentangling it from other hadronic physics uncertainties will require significant theory improvements.
We reveal that the criticality of the chiral phase transition in QCD at the macroscale arises from the microscopic energy levels of its fundamental constituents, the quarks. We establish a novel relation between cumulants of the chiral order parameter (i.e., chiral condensate) and correlations among the energy levels of quarks (i.e., eigenspectra of the massless Dirac operator), which naturally leads to a generalization of the Banks-Casher relation. Based on this novel relation and through (2+1)-flavor lattice QCD calculations using the HISQ action with varying light quark masses in the vicinity of the chiral phase transition, we demonstrate that the correlations among the infrared part of the Dirac eigenspectra exhibit same universal scaling behaviors as expected of the cumulants of the chiral condensate. We find that these universal scaling behaviors extend up to the physical values of the up and down quark masses.
We discuss a universal nonequilibrium enhancement of long-wavelength Goldstone bosons induced by quenches to the broken phase in Model G—the dynamical universality class of an 𝑂(4) antiferromagnet and the chiral phase transition in QCD. Scaling arguments for the coarsening dynamics describing the formation of the chiral condensate predict a parametric enhancement in the infrared spectra of Goldstone bosons, a prediction confirmed by stochastic simulations of the transition. The details of the enhancement are determined by the nonlinear dynamics of a superfluid effective theory, which is a limit of Model G reflecting the broken 𝑂(4) symmetry. Our results translate to a parametric enhancement of low-momentum pions in heavy-ion collisions at the LHC, which are underpredicted in current hydrodynamic models without critical dynamics.
To date numerical simulations of lattice QCD have not found a chiral phase transition of first order that is expected to occur for sufficiently light pions. We show how the restoration of an exact global chiral symmetry can strongly decrease the breaking of the approximate, anomalous U A ( 1 ) symmetry. This is testable on the lattice through simulations for one through four flavors. In QCD a small breaking of the U A ( 1 ) symmetry in the chirally symmetric phase generates novel experimental signals. Published by the American Physical Society 2024
We study the cosmological evolution of a light QCD axion and identify the parameter space to obtain the correct relic dark matter abundance. The axion potential is flattened at the origin, corresponding to the only minimum, while it is unsuppressed at π. These potential features arise by assuming a mirror sector with the strong CP phase $\overline{θ}$ shifted by π compared to the SM sector, which allows the mirror axion potential to be tuned against the usual QCD axion potential. Before the QCD phase transition, assuming the mirror sector is decoupled and much colder than the SM thermal bath, the mirror sector potential dominates, causing the axion to initially roll to a temporary minimum at π. However, after the QCD phase transition, the potential minimum changes, and the axion relaxes from the newly created “hilltop” near π to the CP-conserving minimum at the origin. As the axion adiabatically tracks this shift in the potential minimum through the QCD phase transition, with non-adiabatic evolution near π and 0, it alters the usual prediction of the dark matter abundance. Consequently, this “hilltop” misalignment mechanism opens new regions of axion parameter space, with the correct relic abundance while still solving the strong CP problem, that could be explored in future experiments.
The quantum chromodynamics (QCD) axion arises as the pseudo-Goldstone mode of a spontaneously broken Abelian Peccei-Quinn (PQ) symmetry. If the scale of PQ symmetry breaking occurs below the inflationary reheat temperature and the domain wall number is unity, then there is a unique axion mass that gives the observed dark matter (DM) abundance. Computing this mass has been the subject of intensive numerical simulations for decades since the mass prediction informs laboratory experiments. Axion strings develop below the PQ symmetry-breaking temperature, and as the string network evolves, it emits axions that go on to become the DM. A key ingredient in the axion mass prediction is the spectral index of axion radiation emitted by the axion strings. We compute this index in this Letter using the most precise and accurate large-scale simulations to date of the axion-string network leveraging adaptive mesh refinement to achieve the precision that would, otherwise, require a static lattice with 262,144 3 lattice sites. We find a scale-invariant axion radiation spectrum to within 1% precision and find no evidence that the spectral index of radiation evolves with time. Accounting for axion production from strings prior to the QCD phase transition leads us to predict that the axion mass should be approximately 𝑚 𝑎 ∈ (45, 65) μ eV. However, we provide preliminary evidence that axions are produced in greater quantities from the string-domain-wall network collapse during the QCD phase transition, potentially increasing the mass prediction to as much as 300 μ eV.
Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.
We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.
Large primordial lepton flavor asymmetries with almost vanishing total baryon-minus-lepton number can evade the usual BBN and CMB constraints if neutrino oscillations lead to perfect flavor equilibration. Solving the momentum averaged quantum kinetic equations (QKEs) describing neutrino oscillations and interactions, we perform the first systematic investigation of this scenario, uncovering a rich flavor structure in stark contradiction to the assumption of simple flavor equilibration. We find (i) a particular direction in flavor space, ∆ne ≃ – 2/3 (– 1)∆n μ for normal (inverted) neutrino mass hierarchy, in which the flavor equilibration is efficient and primordial asymmetries are essentially unconstrained, (ii) a minimal washout factor, ∆$n_{e}^{2}$| BBN ≤ 0.03 (0.016) ∑ α ∆$n_{α}^{2}$| ini yielding a conservative estimate for the allowed primordial asymmetries in a generic flavor direction, and (iii) particularly strong or weak washout if one of the initial flavor asymmetries vanishes due to non-adiabatic muon- or electron-driven MSW transitions. These results open up the possibility of a first-order QCD phase transition facilitated by large lepton asymmetries as well as baryogenesis from large and compensated ∆n e = ∆n μ asymmetries. Our systematic approach of deriving momentum averaged QKEs includes collision terms beyond the damping approximation, energy transfer between the neutrino and electron-photon plasma, and provides a fast and reliable way to investigate the impact of primordial lepton asymmetries at the time of BBN. We publicly release the Mathematica code COFLASY-M on https://github.com/mariofnavarro/COFLASY which solves the QKEs numerically.
After presenting the current view of the processes taking place during the cosmological transition from 'quark soup' to normal hadron matter, attention is given to what happens to cosmological nucleosynthesis in the presence of small-scale baryon inhomogeneities. The QCD phase transition is among the plausible sources of this inhomogeneity. It is concluded that the formation of primordial 'quark nuggets' and other cold exotica requires very low entropy regions at the outset, and that even the more modest nonlinearities perturbing nucleosynthesis probably require some ingredient in addition to a quiescent, mildly supercooled transition.
At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.
Dense quantum chromodynamics (QCD) matter may exhibit crystalline phases. Their existence is reflected in a moat regime, where mesonic correlations feature spatial modulations. We study the real-time properties of pions at finite temperature and density in QCD in order to elucidate the nature of this regime. We show that the moat regime arises from particle-hole-like fluctuations near the Fermi surface. This gives rise to a characteristic peak in the spectral function of the pion at nonzero spacelike momentum. This peak can be interpreted as a new quasi particle, the moaton. In addition, our framework also allows us to directly test the stability of the homogeneous chiral phase against the formation of an inhomogeneous condensate in QCD. We find that an inhomogeneous instability is highly unlikely for baryon chemical potentials μ B ≤ 630 MeV . Published by the American Physical Society 2025
A nearby supernova will carry an unprecedented wealth of information about astrophysics, nuclear physics, and particle physics. Because supernova are fundamentally neutrino driven phenomenon, our knowledge about neutrinos—particles that remain quite elusive—will increase dramatically with such a detection. One of the biggest open questions in particle physics is related to the masses of neutrinos. Here, we show how a galactic supernova provides information about the masses of each of the three mass eigenstates , at some precision, and is well probed at JUNO. This information comes from several effects including time delay and the Mikheyev-Smirnov-Wolfenstein effect within the supernova. The time delay feature is strongest during a sharp change in the flux such as the neutronization burst; additional information may also come from a QCD phase transition in the supernova or if the supernova forms a black hole. We consider both standard cases as dictated by local oscillation experiments as well as new physics motivated scenarios where neutrino masses may differ across the galaxy. Published by the American Physical Society 2025
In gauge theories with fundamental matter there is typically no sharp way to distinguish confining and Higgs regimes, e.g. using generalized global symmetries acting on loop order parameters. It is standard lore that these two regimes are continuously connected, as has been explicitly demonstrated in certain lattice and continuum models. We point out that Higgsing and confinement sometimes lead to distinct symmetry protected topological (SPT) phases -- necessarily separated by a phase transition -- for ordinary global symmetries. We present explicit examples in 3+1 dimensions, obtained by adding elementary Higgs fields and Yukawa couplings to QCD while preserving parity P and time reversal T. In a suitable scheme, the confining phases of these theories are trivial SPTs, while their Higgs phases are characterized by non-trivial P- and T-invariant theta-angles θ f , θ g = π for flavor or gravity background gauge fields, i.e. they are topological insulators or superconductors. Finally, we consider conventional three-flavor QCD (without elementary Higgs fields) at finite U(1) B baryon-number chemical potential μ B , which preserves P and T. At very large μ B , three-flavor QCD is known to be a completely Higgsed color superconductor that also spontaneously breaks U(1) B . We argue that this high-density phase is in fact a gapless SPT, with a gravitational theta-angle θ g = π that safely co-exists with the U(1) B Nambu-Goldstone boson. We explain why this SPT motivates unexpected transitions in the QCD phase diagram, as well as anomalous surface modes at the boundary of quark-matter cores inside neutron stars.
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.
Abstract It is difficult to construct a post-inflation QCD axion model that solves the axion quality problem (and hence the Strong CP problem) without introducing a cosmological disaster. In a post-inflation axion model, the axion field value is randomized during the Peccei-Quinn phase transition, and axion domain walls form at the QCD phase transition. We emphasize that the gauge equivalence of all minima of the axion potential (i.e., domain wall number equals one) is insufficient to solve the cosmological domain wall problem. The axion string on which a domain wall ends must exist as an individual object (as opposed to a multi-string state), and it must be produced in the early universe. These conditions are often not satisfied in concrete models. Post-inflation axion models also face a potential problem from fractionally charged relics; solving this problem often leads to low-energy Landau poles for Standard Model gauge couplings, reintroducing the quality problem. We study several examples, finding that models that solve the quality problem face cosmological problems, and vice versa. This is not a no-go theorem; nonetheless, we argue that it is much more difficult than generally appreciated to find a viable post-inflation QCD axion model. Successful examples may have a nonstandard cosmological history (e.g., multiple types of cosmic axion strings of different tensions), undermining the widespread expectation that the post-inflation QCD axion scenario predicts a unique mass for axion dark matter.
When the quark masses are lighter than those in QCD, the standard lore is that a chiral transition of first order must emerge for three light flavors. Recently, however, numerical simulations on the lattice suggest that the chiral transition is of second order in the chiral limit. Using an extended linear sigma model in the mean field approximation, we study the relation between terms which break the anomalous, U ( 1 ) A symmetry and the order of the chiral phase transition, especially how a chiral transition of second order can arise for three massless flavors. We note that in an (unphysical) region of the “Columbia” phase diagram, when the strange quark mass is light and negative, corresponding to topological angle θ = π , the C P symmetry is spontaneously broken. Published by the American Physical Society 2025