An official website of the United States government
ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($$M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$$, $$M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($$\gtrsim 1\,$$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($$M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$$, $$M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($$M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $$\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $$\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant.
more »
« less
Inferring the neutron star equation of state with nuclear-physics informed semiparametric models
Abstract Over the past decade, an abundance of information from neutron-star observations, nuclear experiments and theory has transformed our efforts to elucidate the properties of dense matter. However, at high densities relevant to the cores of neutron stars, substantial uncertainty about the dense matter equation of state (EoS) remains. In this work, we present a semiparametric equation of state framework aimed at better integrating knowledge across these domains in astrophysical inference. We use a Meta-model and realistic crust at low densities, and Gaussian Process extensions at high densities. Comparisons between our semiparametric framework to fully nonparametric EoS representations show that imposing nuclear theoretical and experimental constraints through the Meta-model up to nuclear saturation density results in constraints on the pressure up to twice nuclear saturation density. We also show that our Gaussian Process trained on EoS models with nucleonic, hyperonic, and quark compositions extends the range of EoS explored at high density compared to a piecewise polytropic extension schema, under the requirements of causality of matter and of supporting the existence of heavy pulsars. We find that maximum TOV masses above $$3.2 M_{\odot}$$ can be supported by causal EoS compatible with nuclear constraints at low densities. We then combine information from existing observations of heavy pulsar masses, gravitational waves emitted from binary neutron star mergers, and X-ray pulse profile modeling of millisecond pulsars within a Bayesian inference scheme using our semiparametric EoS prior. With information from all public NICER pulsars (including PSR J0030$$+$$0451, PSR J0740$$+$$6620, PSR J0437-4715, and PSR J0614-3329), we find an astrophysically favored pressure at two times nuclear saturation density of $$P(2\rho_{\rm nuc}) = 1.98^{+2.13}_{-1.08}\times10^{34}$$ dyn/cm$$^{2}$$, a radius of a $$1.4 M_{\odot}$$ neutron star value of $$R_{1.4} = 11.4^{+0.98}_{-0.60}$$\;km, and $$M_{\rm max} = 2.31_{-0.23}^{+0.35} M_{\odot}$$ at the 90\% credible level.
more »
« less
- PAR ID:
- 10641641
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Classical and Quantum Gravity
- ISSN:
- 0264-9381
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
For densities beyond nuclear saturation, there is still a large uncertainty in the equations of state (EoSs) of dense matter that translate into uncertainties in the internal structure of neutron stars. The MUSES Calculation Engine provides a free and open-source composable workflow management system, which allows users to calculate the EoSs of dense and hot matter that can be used, e.g., to describe neutron stars. For this work, we make use of two MUSES EoS modules, i.e., Crust Density Functional Theory and Chiral Mean Field model, with beta-equilibrium with leptons enforced in the Lepton module, then connected by the Synthesis module using different functions: hyperbolic tangent, generalized Gaussian, bump, and smoothstep. We then calculate stellar structure using the QLIMR module and discuss how the different interpolating functions affect our results.more » « less
-
Neutron stars provide a window into the properties of dense nuclear matter. Several recent observational and theoretical developments provide powerful constraints on their structure and internal composition. Among these are the first observed binary neutron star merger, GW170817, whose gravitational radiation was accompanied by electromagnetic radiation from a short γ-ray burst and an optical afterglow believed to be due to the radioactive decay of newly minted heavy r-process nuclei. These observations give important constraints on the radii of typical neutron stars and on the upper limit to the neutron star maximum mass and complement recent pulsar observations that established a lower limit. Pulse-profile observations by the Neutron Star Interior Composition Explorer (NICER) X-ray telescope provide an independent, consistent measure of the neutron star radius. Theoretical many-body studies of neutron matter reinforce these estimates of neutron star radii. Studies using parameterized dense matter equations of state (EOSs) reveal several EOS-independent relations connecting global neutron star properties.more » « less
-
Abstract This study investigates the properties of symmetric and asymmetric nuclear matter using the relativistic Brueckner–Hartree–Fock formalism, examining both zero and finite temperatures up to 70 MeV. Employing the full Dirac space, we incorporate three Bonn potentials (A, B, and C), which account for meson masses, coupling strengths, cutoff parameters, and form factors. The calculated properties of asymmetric nuclear matter form the basis for constructing equation-of-state (EOS) models tailored for neutron stars. These models, in turn, enable the computation of bulk properties for nonrotating, uniformly rotating, and differentially rotating neutron stars. Notably, the EOS models studied in this paper are sufficiently versatile to accommodate the mass of the most massive neutron star ever detected, PSR J0952–0607, estimated to be 2.35 ± 0.17M⊙. Furthermore, they yield masses and radii for PSR J0030+451 that align with the confidence intervals established for this pulsar.more » « less
-
Recent experimental and ab initio theory investigations of the 208Pb neutron skin thickness have the potential to inform the neutron star equation of state. In particular, the strong correlation between the 208Pb neutron skin thickness and the pressure of neutron matter at normal nuclear densities leads to modified predictions for the radii, tidal deformabilities, and moments of inertia of typical 1.4M⊙ neutron stars. In the present work, we study the relative impact of these recent analyses of the 208Pb neutron skin thickness on bulk properties of neutron stars within a Bayesian statistical analysis. Two models for the equation of state prior are employed in order to highlight the role of the highly uncertain high-density equation of state. From our combined Bayesian analysis of nuclear theory, nuclear experiment, and observational constraints on the dense matter equation of state, we find at the 90% credibility level R1.4=12.36−0.73+0.38 km for the radius of a 1.4M⊙ neutron star, R2.0=11.96−0.71+0.94 km for the radius of a 2.0M⊙ neutron star, Λ1.4=440−144+103 for the tidal deformability of a 1.4M⊙ neutron star, and I1.338=1.425−0.146+0.074×1045gcm2 for the moment of inertia of PSR J0737-3039A whose mass is 1.338M⊙.more » « less
