Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

We provide an end-to-end exploration of a distinct modified gravitational theory in Jordan-Brans-Dicke (JBD) gravity, from an analytical and numerical description of the background expansion and linear perturbations, to the nonlinear regime captured with a hybrid suite of N-body simulations, to the cosmological constraints from existing probes of the expansion history, the large-scale structure, and the cosmic microwave background (CMB). We have focused on JBD gravity as it both approximates a wider class of Horndeski scalar-tensor theories on cosmological scales and allows us to adequately model the nonlinear corrections to the matter power spectrum. In a combined analysis of the Planck 2018 CMB temperature, polarization, and lensing reconstruction, together with Pantheon supernova distances and the Baryon Oscillation Spectroscopic Survey (BOSS) measurements of baryon acoustic oscillation distances, the Alcock-Paczynski effect, and the growth rate, we constrain the JBD coupling constant to ωBD>970 (95% confidence level; C.L.) in agreement with the General Relativistic expectation given by ωBD→∞. In the unrestricted JBD model, where the effective gravitational constant at present, Gmatter/G, is additionally varied, increased dataset concordance (e.g., within 1σ agreement in S8=σ8ωm/0.3) enables us to further include the combined ("3×2pt") dataset of cosmic shear, galaxy-galaxy lensing, and overlapping redshift-space galaxy clustering from the Kilo Degree Survey and the 2-degree Field Lensing Survey (KiDS×2dFLenS). In analyzing the weak lensing measurements, the nonlinear corrections due to baryons, massive neutrinos, and modified gravity are simultaneously modeled and propagated in the cosmological analysis for the first time. In the joint analysis of all datasets, we constrain ωBD>1540 (95% C.L.), Gmatter/G=0.997±0.029, the sum of neutrino masses, mν<0.12 eV (95% C.L.), and the baryonic feedback amplitude, B<2.8 (95% CL), all in agreement with the standard model expectation. In fixing the sum of neutrino masses, the lower bound on the coupling constant strengthens to ωBD>1460 and ωBD>2230 (both at 95% C.L.) in the restricted and unrestricted JBD models, respectively. We explore the impact of the JBD modeling choices, and show that a more restrictive parametrization of the coupling constant degrades the neutrino mass bound by up to a factor of three. In addition to the improved concordance between KiDS×2dFLenS and Planck, the tension in the Hubble constant between Planck and the direct measurement of Riess et al. (2019) is reduced to ∼3σ; however, we find no substantial model selection preference for JBD gravity relative to ΛCDM. We further show that a positive shift in the effective gravitational constant suppresses the CMB damping tail, which might complicate future inferences of small-scale physics, given its degeneracy with the primordial helium abundance, the effective number of neutrinos, and the running of the spectral index.

Abstract:

CMB-S4—the next-generation ground-based cosmic microwave background (CMB) experiment—is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the universe. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semianalytic projection tool, targeted explicitly toward optimizing constraints on the tensor-to-scalar ratio, r, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2–3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments, given a desired scientific goal. To form a closed-loop process, we couple this semianalytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for r > 0.003 at greater than 5σ, or in the absence of a detection, of reaching an upper limit of r < 0.001 at 95% CL.

Abstract:

The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields, and fractions of supernova energy deposited as CRs (χCR), accounting for CR anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: pSNR=2.87×105(χCR+1)4.82(ncm−3)−0.196M⊙ km s−1. Assuming a 10 per cent fraction of SN energy in CRs leads to a 50 per cent boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and CR outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarization.