-
Inferring the Ionizing Photon Contributions of High-Redshift Galaxies to Reionization with JWST NIRCam Photometry
Authors:
Nicholas Choustikov,
Richard Stiskalek,
Aayush Saxena,
Harley Katz,
Julien Devrient,
Adrianne Slyz
Abstract:
JWST observations have the potential to provide unprecedented constraints on the history of reionization and the sources responsible for the ionizing photons due to the detection of large populations of faint galaxies at $z\gg6$. Modelling reionization requires knowing both the number of ionizing photons that are produced by galaxies and the fraction of those photons that escape into the intergala…
▽ More
JWST observations have the potential to provide unprecedented constraints on the history of reionization and the sources responsible for the ionizing photons due to the detection of large populations of faint galaxies at $z\gg6$. Modelling reionization requires knowing both the number of ionizing photons that are produced by galaxies and the fraction of those photons that escape into the intergalactic medium. Observational estimates of these values generally rely on spectroscopy for which large samples with well-defined selection functions remain limited. To overcome this challenge, we present an implicit likelihood inference (ILI) pipeline trained on mock photometry to predict the escaped ionizing luminosity of individual galaxies ($\dot{n}_{\rm ion}$) based on photometric magnitudes and redshifts. Compared to traditional SED-fitting methods, the new ILI pipeline is consistently more accurate and significantly faster. We deploy the method on a sample of 4,559 high-redshift galaxies from the JADES Deep survey, finding a gentle redshift evolution of $\log_{10}(\dot{n}_{\rm ion}) = (0.08\pm0.01)z + (51.60\pm0.06)$, with late-time values for $\dot{N}_{\rm ion}$ consistent with theoretical models and observations. We measure the evolution of the volume-averaged ionized fraction and optical depth to find that observed populations of star-forming galaxies are capable of driving reionization to completion at $z\sim 5.3$ without the need for AGN or other exotic sources. The $20\%$ of UV-brightest galaxies ($M_{\rm UV}<-18.5$) in our sample can reionize only $\sim30\%$ of the survey volume, demonstrating that faint LyC emitters are crucial for reionization.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
Evaluating the reconstruction of individual haloes in constrained cosmological simulations
Authors:
Richard Stiskalek,
Harry Desmond,
Julien Devriendt,
Adrianne Slyz
Abstract:
Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation histories. We introduce a method for determining the extent to which individual haloes are reliably reconstructed between constrained simulations, and apply it to the Constrained Simulations in BORG (CSiBORG) suite of…
▽ More
Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation histories. We introduce a method for determining the extent to which individual haloes are reliably reconstructed between constrained simulations, and apply it to the Constrained Simulations in BORG (CSiBORG) suite of $101$ high-resolution realisations across the posterior probability distribution of initial conditions from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. The method is based on the overlap of the initial Lagrangian patch of a halo in one simulation with those in another, and therefore measures the degree to which the haloes' particles are initially coincident. By this metric we find consistent reconstructions of $M\gtrsim10^{14}~M_\odot / h$ haloes across the CSiBORG simulations, indicating that the constraints from the BORG algorithm are sufficient to pin down the masses, positions and peculiar velocities of clusters to high precision. The effect of the constraints tapers off towards lower mass however, and the halo spins and concentrations are largely unconstrained at all masses. We document the advantages of evaluating halo consistency in the initial conditions, describe how the method may be used to quantify our knowledge of the halo field given galaxy survey data analysed through the lens of probabilistic inference machines such as BORG, and describe applications to matched but unconstrained simulations.
△ Less
Submitted 31 October, 2023;
originally announced October 2023.
-
Probing general relativistic spin-orbit coupling with gravitational waves from hierarchical triple systems
Authors:
Marius A. Oancea,
Richard Stiskalek,
Miguel Zumalacárregui
Abstract:
Wave packets propagating in inhomogeneous media experience a coupling between internal and external degrees of freedom and, as a consequence, follow spin-dependent trajectories. These are known as spin Hall effects, which are well known in optics and condensed matter physics. Similarly, the gravitational spin Hall effect is expected to affect the propagation of gravitational waves on curved spacet…
▽ More
Wave packets propagating in inhomogeneous media experience a coupling between internal and external degrees of freedom and, as a consequence, follow spin-dependent trajectories. These are known as spin Hall effects, which are well known in optics and condensed matter physics. Similarly, the gravitational spin Hall effect is expected to affect the propagation of gravitational waves on curved spacetimes. In this general-relativistic setup, the curvature of spacetime acts as impurities in a semiconductor or inhomogeneities in an optical medium, leading to a frequency- and polarization-dependent propagation of wave packets. In this letter, we study this effect for strong-field lensed gravitational waves generated in hierarchical triple black hole systems in which a stellar-mass binary merges near a more massive black hole. We calculate how the gravitational spin Hall effect modifies the gravitational waveforms and show its potential for experimental observation. If detected, these effects will bear profound implications for astrophysics and tests of general relativity.
△ Less
Submitted 4 July, 2023;
originally announced July 2023.
-
On the fundamentality of the radial acceleration relation for late-type galaxy dynamics
Authors:
Richard Stiskalek,
Harry Desmond
Abstract:
Galaxies have been observed to exhibit a level of simplicity unexpected in the complex galaxy formation scenario posited by standard cosmology. This is particularly apparent in their dynamics, where scaling relations display much regularity and little intrinsic scatter. However, the parameters responsible for this simplicity have not been identified. Using the Spitzer Photometry & Accurate Rotatio…
▽ More
Galaxies have been observed to exhibit a level of simplicity unexpected in the complex galaxy formation scenario posited by standard cosmology. This is particularly apparent in their dynamics, where scaling relations display much regularity and little intrinsic scatter. However, the parameters responsible for this simplicity have not been identified. Using the Spitzer Photometry & Accurate Rotation Curves galaxy catalogue, we argue that the radial acceleration relation (RAR) between galaxies' baryonic and total dynamical accelerations is the fundamental $1$-dimensional correlation governing the radial (in-disk) dynamics of late-type galaxies. In particular, we show that the RAR cannot be tightened by the inclusion of any other available galaxy property, that it is the strongest projection of galaxies' radial dynamical parameter space, and that all other statistical radial dynamical correlations stem from the RAR plus the non-dynamical correlations present in our sample. We further provide evidence that the RAR's fundamentality is unique in that the second most significant dynamical relation does not possess any of these features. Our analysis reveals the root cause of the correlations present in galaxies' radial dynamics: they are nothing but facets of the RAR. These results have important ramifications for galaxy formation theory because they imply that to explain statistically late-type galaxy dynamics within the disk it is necessary and sufficient to explain the RAR and lack of any significant, partially independent correlation. While simple in some modified dynamics models, this poses a challenge to standard cosmology.
△ Less
Submitted 18 September, 2023; v1 submitted 31 May, 2023;
originally announced May 2023.
-
Frequency- and polarization-dependent lensing of gravitational waves in strong gravitational fields
Authors:
Marius A. Oancea,
Richard Stiskalek,
Miguel Zumalacárregui
Abstract:
The propagation of gravitational waves can be described in terms of null geodesics by using the geometrical optics approximation. However, at large but finite frequencies the propagation is affected by the spin-orbit coupling corrections to geometrical optics, known as the gravitational spin Hall effect. Consequently, gravitational waves follow slightly different frequency- and polarization-depend…
▽ More
The propagation of gravitational waves can be described in terms of null geodesics by using the geometrical optics approximation. However, at large but finite frequencies the propagation is affected by the spin-orbit coupling corrections to geometrical optics, known as the gravitational spin Hall effect. Consequently, gravitational waves follow slightly different frequency- and polarization-dependent trajectories, leading to dispersive and birefringent phenomena. We study the potential for detecting the gravitational spin Hall effect in hierarchical triple black hole systems, consisting of an emitting binary orbiting a more massive body acting as a gravitational lens. We calculate the difference in time of arrival with respect to the geodesic propagation and find that it follows a simple power-law dependence on frequency with a fixed exponent. We calculate the gravitational spin Hall-corrected waveform and its mismatch with respect to the original waveform. The waveform carries a measurable imprint of the strong gravitational field if the source, lens, and observer are sufficiently aligned, or for generic observers if the source is close enough to the lens. We present constraints on dispersive time delays from GWTC-3, translated from limits on Lorentz invariance violation. Finally, we address the sensitivity of current and future ground detectors to dispersive lensing. Our results demonstrate that the gravitational spin Hall effect can be detected, providing a novel probe of general relativity and the environments of compact binary systems.
△ Less
Submitted 17 May, 2024; v1 submitted 14 September, 2022;
originally announced September 2022.
-
The scatter in the galaxy-halo connection: a machine learning analysis
Authors:
Richard Stiskalek,
Deaglan J. Bartlett,
Harry Desmond,
Dhayaa Anbajagane
Abstract:
We apply machine learning, a powerful method for uncovering complex correlations in high-dimensional data, to the galaxy-halo connection of cosmological hydrodynamical simulations. The mapping between galaxy and halo variables is stochastic in the absence of perfect information, but conventional machine learning models are deterministic and hence cannot capture its intrinsic scatter. To overcome t…
▽ More
We apply machine learning, a powerful method for uncovering complex correlations in high-dimensional data, to the galaxy-halo connection of cosmological hydrodynamical simulations. The mapping between galaxy and halo variables is stochastic in the absence of perfect information, but conventional machine learning models are deterministic and hence cannot capture its intrinsic scatter. To overcome this limitation, we design an ensemble of neural networks with a Gaussian loss function that predict probability distributions, allowing us to model statistical uncertainties in the galaxy-halo connection as well as its best-fit trends. We extract a number of galaxy and halo variables from the Horizon-AGN and IllustrisTNG100-1 simulations and quantify the extent to which knowledge of some subset of one enables prediction of the other. This allows us to identify the key features of the galaxy-halo connection and investigate the origin of its scatter in various projections. We find that while halo properties beyond mass account for up to 50 per cent of the scatter in the halo-to-stellar mass relation, the prediction of stellar half-mass radius or total gas mass is not substantially improved by adding further halo properties. We also use these results to investigate semi-analytic models for galaxy size in the two simulations, finding that assumptions relating galaxy size to halo size or spin are not successful.
△ Less
Submitted 8 June, 2022; v1 submitted 28 February, 2022;
originally announced February 2022.
-
The dependence of subhalo abundance matching on galaxy photometry and selection criteria
Authors:
Richard Stiskalek,
Harry Desmond,
Thomas Holvey,
Michael G. Jones
Abstract:
Subhalo abundance matching (SHAM) is a popular technique for assigning galaxy mass or luminosity to haloes produced in N-body simulations. The method works by matching the cumulative number functions of the galaxy and halo properties, and is therefore sensitive both to the precise definitions of those properties and to the selection criteria used to define the samples. Further dependence follows w…
▽ More
Subhalo abundance matching (SHAM) is a popular technique for assigning galaxy mass or luminosity to haloes produced in N-body simulations. The method works by matching the cumulative number functions of the galaxy and halo properties, and is therefore sensitive both to the precise definitions of those properties and to the selection criteria used to define the samples. Further dependence follows when SHAM parameters are calibrated with galaxy clustering, which is known to depend strongly on the manner in which galaxies are selected. In this paper we introduce a new parametrisation for SHAM and derive the best-fit SHAM parameters as a function of various properties of the selection of the galaxy sample and of the photometric definition, including Sérsic vs Petrosian magnitudes, stellar masses vs r-band magnitudes and optical (SDSS) vs HI (ALFALFA) selection. In each case we calculate the models' goodness-of-fit to measurements of the projected two-point galaxy correlation function. In the optically-selected samples we find strong evidence that the scatter in the galaxy-halo connection increases towards the faint end, and that AM performs better with luminosity than stellar mass. The SHAM parameters of optically- and HI-selected galaxies are mutually exclusive, with the latter suggesting the importance of properties beyond halo mass. We provide best-fit parameters for the SHAM galaxy-halo connection as a function of each of our input choices, extending the domain of validity of the model while reducing potential systematic error in its use.
△ Less
Submitted 21 July, 2021; v1 submitted 7 January, 2021;
originally announced January 2021.
-
Are stellar--mass binary black hole mergers isotropically distributed?
Authors:
Richard Stiskalek,
John Veitch,
Chris Messenger
Abstract:
The Advanced LIGO and Advanced Virgo gravitational wave detectors have detected a population of binary black hole mergers in their first two observing runs. For each of these events we have been able to associate a potential sky location region represented as a probability distribution on the sky. Thus, at this point we may begin to ask the question of whether this distribution agrees with the iso…
▽ More
The Advanced LIGO and Advanced Virgo gravitational wave detectors have detected a population of binary black hole mergers in their first two observing runs. For each of these events we have been able to associate a potential sky location region represented as a probability distribution on the sky. Thus, at this point we may begin to ask the question of whether this distribution agrees with the isotropic model of the Universe, or if there is any evidence of anisotropy. We perform Bayesian model selection between an isotropic and a simple anisotropic model, taking into account the anisotropic selection function caused by the underlying antenna patterns and sensitivity of the interferometers over the sidereal day. We find an inconclusive Bayes factor of $1.3:1$, suggesting that the data from the first two observing runs is insufficient to pick a preferred model. However, the first detections were mostly poorly localised in the sky (before the Advanced Virgo joined the network), spanning large portions of the sky and hampering detection of potential anisotropy. It will be appropriate to repeat this analysis with events from the recent third LIGO observational run and a more sophisticated cosmological model.
△ Less
Submitted 19 November, 2020; v1 submitted 5 March, 2020;
originally announced March 2020.