Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

It has only recently become possible to simulate the full nonlinear dynamics of binary black holes in scalar-Gauss-Bonnet theories of gravity. The simulations remain technically challenging and evolutions of unequal mass binaries in particular have been difficult to follow through the merger. Even when the merger is successful, accurately quantifying the physical dephasing, as opposed to contributions from transients in the initial data and gauge adjustments, remains difficult. We show the first full simulations of 2∶1 and 3∶1 binaries through merger, and we discuss how specific choices in the setup affect the dephasing observed and our ability to obtain reliable results. In cases with weaker couplings, we match the expected PN value for the dephasing, whereas for larger couplings, eccentricity introduced by the initial data transients can lead to artificial deviations. Our work highlights the need for improvements in the initial data methods used, to ensure reliable waveforms are obtained for data analysis in beyond-GR models.

Abstract:

Cosmic ray (CR) feedback has been proposed as a powerful mechanism for driving warm gas outflows in galaxies. We use cosmological magnetohydrodynamic simulations to investigate the impact of CR feedback on neutral hydrogen (H i) in a 1011 M⊙ dark matter halo at 2 < z < 4. To this end, we postprocess the simulations with ionizing radiative transfer and perform Monte Carlo Lyman-α (Lyα) transfer calculations. CR feedback reduces H i column densities around young stars, thereby allowing more Lyα photons to escape and consequently offering a better match to the Lyα luminosities of observed Lyα emitters. Although galaxies with CR-driven outflows have more extended H i in the circumgalactic medium, two Lyα line properties sensitive to optical depth and gas kinematic—the location of the red peak relative to the Lyα line center in velocity space (vred) and relative strength of the blue-to-red peaks (B/R)—cannot distinguish between the CR-driven and non-CR simulations. This is because Lyα photons propagate preferentially along low H i density channels created by the ionizing radiation, thereby limiting the scattering with volume-filling H i. In contrast, the observed low flux ratios between the valley and peak and the surface brightness profiles are better reproduced in the model with CR-driven outflows because the Lyα photons interact more before escaping, rather than being destroyed by dust as is the case in the non-CR simulation. We discuss the potential cause of the paucity of sightlines in simulations that exhibit prominent red peaks and large vred, which may require the presence of more volume-filling H i.