Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NEUTRON STARS”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Symmetry-Energy Dependence of the Bulk Viscosity of Nuclear Matter

We clarify how the weak-interaction-driven bulk viscosity ζ and the bulk relaxation time τ Π of neutrino-transparent n p e matter depend on the nuclear symmetry energy. We show that, at saturation density, the equation-of-state dependence of these transport quantities is fully determined by the experimentally constrained nuclear symmetry energy S and its slope L . Variations of L can change the bulk viscosity by orders of magnitude, which can affect both the dissipative and the conservative tidal response of neutron stars. Furthermore, this suggests that both conservative and dissipative effects encoded in the gravitational-wave signatures of binary neutron star inspirals may help constrain nuclear symmetry-energy properties.

Neutron stars & pulsars↗
DOE OSTI2021

Structure of Stellar Remnants with Coupling to a Light Scalar

In this paper we study how a Yukawa coupling of the Standard Model fermions to a light scalar field effects the stellar structure of cold stellar remnants such as neutron stars. We elucidate the stellar structure phenomenology using a simple model of a massive scalar coupled to a single dominant fermion with no other interactions. For a broad scalar mass range ($10^{-10}\,\mathrm{eV}\ll m_\phi\ll10^3\,\mathrm{eV}$ for neutron stars) we show that the equation-of-state and stellar structure depends only the effective coupling $\mathfrak{g}=\frac{g_f\,m_f}{m_\phi}$, where $g_f$ is the Yukawa coupling, $m_f$ the fermion mass and $m_\phi$ is the scalar kinematic mass at nuclear densities. If $\mathfrak{g}>\mathcal{O}(1)$ the Yukawa coupled matter exhibits various anomalous behaviors including hydrodynamic instability, negative pressure, distinct phases (soft and hard) of matter with sharp phase boundaries between them and with vacuum. These anomalies can lead to stars consisting of only soft, only hard or hybrid of soft and hard matter. These stars can have varying signs of the slope of the mass-radius relation, anomalously large and small masses, gaps in allowed radii, multiple radii for the same mass, very thin crusts and radiate anomalously large amounts of energy when they form (in the form of neutrinos for neutron stars). To the extent that these anomalies have not and/or will not be observed limits the effective coupling to $\mathfrak{g}<\mathcal{O}(1)$. We argue this phenomenology is generic to realistic models of stars with Yukawa coupled matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Constraining ΛCDM cosmological parameters with Einstein Telescope mock data

ABSTRACT We investigate the capability of Einstein Telescope to constrain the cosmological parameters of the non-flat ΛCDM cosmological model. Two types of mock data sets are considered depending on whether or not a short gamma-ray burst is detected, and associated with the gravitational wave emitted by binary neutron stars merger, using the THESEUS satellite. Depending on the mock data set, two statistical estimators are applied: one assumes that the redshift is known, while the other marginalizes over it assuming a specific redshift prior distribution. We demonstrate that (i) using mock catalogues collecting gravitational wave signals emitted by binary neutron stars systems to which a short gamma-ray burst has been associated, Einstein Telescope may achieve an accuracy on the cosmological parameters of $\sigma _{H_0}\approx 0.40$ km s−1 Mpc−1, $\sigma _{\Omega _{k,0}}\approx 0.09$, and $\sigma _{\Omega _{\Lambda ,0}}\approx 0.07$; while (ii) using mock catalogues collecting all gravitational wave signals emitted by binary neutron stars systems for which an electromagnetic counterpart has not been detected, Einstein Telescope may achieve an accuracy on the cosmological parameters of $\sigma _{H_0}\approx 0.04$ km s−1 Mpc−1, $\sigma _{\Omega _{k,0}}\approx 0.01$, and $\sigma _{\Omega _{\Lambda ,0}}\approx 0.01$, once the redshift probability distribution of GW events is known from from population synthesis simulations and/or the measure of the tidal deformability parameter. These results show an improvement of a factor 2–75 with respect to earlier results using complementary data sets.

Califano, Matteo↗
DOE OSTI2021

Detailed examination of astrophysical constraints on the symmetry energy and the neutron skin of 208 Pb with minimal modeling assumptions

We report the symmetry energy and its density dependence are pivotal for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, $R^{^{208}Pb}_{skin}$, implying a symmetry-energy slope parameter L of 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with chiral effective field theory (χEFT) and PREX-II constraints yields S 0 = 33.0$^{+2.0}_{-1.8}$ MeV, L = 53$^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin}$ = 0.17 $^{+0.04}_{-0.04}$ fm. We also examine the consistency of several individual χEFT calculations with astrophysical observations and terrestrial experiments. We find that there is only mild tension between χ EFT, astrophysical data, and PREX-II's $R^{^{208}Pb}_{skin}$ measurement (p value = 12.3%) and that there is excellent agreement between χEFT, astrophysical data, and other nuclear experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of 208 Pb with Minimal Modeling Assumptions

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, implying a slope parameter L = 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Finally, combining astrophysical data with PREX-II and chiral effective field theory constraints yields S 0 = $33.0^{+2.0}_{-1.8}$ MeV , L = $53^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin} = 0.17^{+0.04}_{-0.04}$ fm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Finding Structure in the Speed of Sound of Supranuclear Matter from Binary Love Relations

Analyses that connect observations of neutron stars with nuclear-matter properties can rely on equation-of-state insensitive relations. We show that the slope of the binary Love relations (between the tidal deformabilities of binary neutron stars) encodes the baryon density at which the speed of sound rapidly changes. Twin stars lead to relations that present a signature “hill,” “drop,” and “swoosh” due to the second (mass-radius) stable branch, requiring a new description of the binary Love relations. Together, these features can reveal new properties and phases of nuclear matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Direct astrophysical tests of chiral effective field theory at supranuclear densities

Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS equipped with reliable uncertainty estimates from chiral effective field theory ( x EFT) allow us to bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on x EFT to directly constrain the underlying physical properties of the compact objects without introducing modeling systematics. We discuss how the data alone constrain the EoS at high densities when we condition on x EFT at low densities. Here, we also demonstrate how to exploit astrophysical data to directly test the predictions of x EFT for the EoS up to twice nuclear saturation density, in order to estimate the density at which these predictions might break down. We nd that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor x EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While x EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at or slightly above nuclear saturation density. Additionally, for these QMC calculations, we marginalize over the uncertainty in the density at which x EFT begins to break down, constraining the radius of a 1.4M $\bigodot$ neutron star to R 1.4 = $11.40^{+1.38}_{–1.04}$ ($12.54^{+0.71}_{–0.63}$) km and the pressure at twice nuclear saturation density to p(2n sat ) = $14.2^{+18.1}_{–8.4}$ ($28.7^{+15.3}_{–15.0}$) MeV=fm 3 with massive pulsar and gravitational-wave (and NICER) data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Particle-in-cell Simulations of Mildly Relativistic Outflows in Kilonova Emissions

Abstract The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet and thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation, and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL’s vector particle-in-cell (VPIC) code, we have run simulations of such mildly relativistic, collisionless, weakly magnetized plasmas and computed the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton-to-electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
DOE OSTI2021

The Challenges Ahead for Multimessenger Analyses of Gravitational Waves and Kilonova: A Case Study on GW190425

In recent years, there have been significant advances in multimessenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, in this work, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so-called kilonova, and particularly on black hole–neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole or a neutron star. We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Polarization formalism for ALP-induced X-ray emission from magnetars

Missions like NASA's Imaging X-ray Polarimetry Explorer (IXPE) are poised to provide an unprecedented view of the Universe in polarized X-rays. Polarization probes physical anisotropies, a fact exploited by particle physicists to look for the anisotropic a E∙B operator in the axion-like-particle (ALP) Lagrangian. Such studies have typically focused on polarization in the radio and microwaves, through local or cosmic birefringence effects. To such polarization studies we add X-rays emanating from magnetars — a class of neutron stars with near-critical strength magnetic fields — that are important targets for IXPE. ALPs produced in the neutron star core convert to X-rays in the magnetosphere; such X-rays are polarized along the direction parallel to the dipolar magnetic field at the point of conversion. Here, we develop the full theoretical formalism for ALP-induced polarization in the presence of dipolar magnetic fields. For uncorrelated photon and ALP production mechanisms, we completely disentangle the ALP contributions to the Stokes parameters in terms of the ALP intensity, the ALP-to-photon conversion probability, and the ALP-induced birefringence. In the proper limit, our results demonstrate that the inclusion of ALPs suppresses the observed degree of circular polarization compared to its pure astrophysical value. Our results can also be used to impose limits on ALP couplings with IXPE polarization data from magnetars 4U 0142+61 and 1RXS J170849.0-400910, the subject of upcoming work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Supernova Axions Convert to Gamma Rays in Magnetic Fields of Progenitor Stars

It has long been established that axions could have been produced within the nascent proto-neutron star formed following the type II supernova SN1987A, escaped the star due to their weak interactions, and then converted to gamma rays in the Galactic magnetic fields; the nonobservation of a gamma-ray flash coincident with the neutrino burst leads to strong constraints on the axion-photon coupling for axion masses m a ≲10 –10 eV. In this Letter, we use SN1987A to constrain higher mass axions, all the way to m a ~10 –3 eV, by accounting for axion production from the Primakoff process, nucleon bremsstrahlung, and pion conversion along with axion-photon conversion on the still-intact magnetic fields of the progenitor star. Moreover, we show that gamma-ray observations of the next Galactic supernova, leveraging the magnetic fields of the progenitor star, could detect quantum chromodynamics axions for masses above roughly 50 μ⁢eV, depending on the supernova. We propose a new full-sky gamma-ray satellite constellation that we call the GALactic AXion Instrument for Supernova (GALAXIS) to search for such future signals along with related signals from extragalactic neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Revealing the strength of three-nucleon interactions with the proposed Einstein Telescope

Nuclear systems, ranging from atomic nuclei to dense matter probed in neutron stars, are governed by strong interactions. Three-nucleon forces have been found to be a crucial ingredient for the reliable description of these systems. Here, in this work, we explore how astrophysical data on neutron stars and their mergers from current and next-generation observatories will enable us to distinguish nuclear Hamiltonians. In particular, we investigate two different nuclear Hamiltonians that have been adjusted to reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. We find that no significant constraints can be obtained from current data, but that the proposed Einstein Telescope could provide strong evidence to distinguish among these Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Investigating effects of relativistic kinematics, dimensionality, interactions, and short-range correlations on the ratio of quartic over quadratic nuclear symmetry energies

While ample evidence for the so-called empirical parabolic law of the Equation of State (EOS) of isospin asymmetric nuclear matter (ANM) has been obtained in many studies within both non-relativistic and relativistic nuclear many-body theories using various interactions, it has been unclear if there is any fundamental physics reason for the small quartic symmetry energy compared to the quadratic one even as the ANM approaches pure neutron matter. Within both relativistic and non-relativistic Free Fermi Gas (FFG) models in coordinate spaces of arbitrary dimension d with and without considering Short-Range Correlations (SRC) as well as non-linear Relativistic Mean Field (RMF) models, we study effects of relativistic kinematics, dimensionality, interactions and SRC on the ratio Ψ(ρ) of quartic over quadratic symmetry energies in ANM EOSs. We found that the ratio Ψ(ρ) in the FFG model depends strongly on the dimension d. While it is very small already in the normal 3D space, it could be even smaller in spaces with reduced dimensions for sub-systems of particles in heavy-ion reactions and/or whole neutron stars due to constraints, collectivities and/or symmetries. Here, we also found that the ratio Ψ(ρ) could theoretically become very large only at the ultra-relativistic limit far above the density reachable in neutron stars. On the other hand, nuclear interaction directly and/or indirectly through SRC-induced high-momentum nucleons affect significantly the density dependence of Ψ(ρ) compared to the relativistic FFG model prediction. The SRC affects significantly not only the kinetic energy of symmetric nuclear matter but also the ratio Ψ(ρ) while the relativistic corrections are found negligible. The results may help better understand the EOS of dense neutron-rich matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effect of the Nuclear Equation of State on Relativistic Turbulence-induced Core-collapse Supernovae

Abstract The nuclear equation of state (EOS) is an important component in the evolution of core-collapse supernovae. In this paper we make a survey of various EOSs in the literature and analyze their effect on spherical core-collapse models in which the effects of three-dimensional turbulence is modeled by a general relativistic formulation of Supernova Turbulence In Reduced-dimensionality (STIR). We show that the viability of the explosion is quite EOS dependent and that it best correlates with the early-time interior entropy density of the proto–neutron star. We check that this result is not progenitor dependent, although the lowest-mass progenitors show different explosion properties, due to the different pre-collapse nuclear composition. Larger central entropies also induce more vigorous proto–neutron star convection in our one-dimensional turbulence model, as well as a wider convective layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Covariant Density Functional Theory in Nuclear Physics and Astrophysics

How does subatomic matter organize itself? Neutron stars are cosmic laboratories uniquely poised to answer this fundamental question that lies at the heart of nuclear science. Newly commissioned rare isotope facilities, telescopes operating across the entire electromagnetic spectrum, and ever more sensitive gravitational wave detectors will probe the properties of neutron-rich matter with unprecedented precision over an enormous range of densities. A coordinated effort between observation, experiment, and theoretical research is of paramount importance for realizing the full potential of these investments. Theoretical nuclear physics provides valuable insights into the properties of neutron-rich matter in regimes that are not presently accessible to experiment or observation. In particular, nuclear density functional theory is likely the only tractable framework that can bridge the entire nuclear landscape by connecting finite nuclei to neutron stars. This compelling connection is the main scope of the present review.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Bulk viscosity of nuclear matter with pions in the neutrino-trapped regime

Recent work has shown that a population of thermal pions could modify the equation of state and transport properties of hot and dense neutron-rich matter and introduce new reaction pathways to change the proton fraction. In this article we study their impact on the bulk viscosity of dense matter, focusing on the neutrino-trapped regime that would be realized in neutron star mergers and supernovae. We find that the presence of a thermal population of pions alters the bulk viscosity by modifying the EoS (via the susceptibilities) and by providing new reaction pathways to achieve beta equilibrium. In neutron star merger conditions, the bulk viscosity in neutrino-trapped npeμ matter (without pions) has its peak at temperatures of at most a couple MeV and is quite small at temperatures of tens of MeV. We find that thermal pions enhance the low-temperature peak of the bulk viscosity by a factor of a few and shift it to slightly lower temperatures. At higher temperatures, where the pion abundance is large but the bulk viscosity is traditionally small, pions can increase the bulk viscosity by an order of magnitude or more, although it is still orders of magnitude smaller than its peak value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗