An official website of the United States government
Abstract The observed prevalence of galaxies exhibiting bursty star formation histories (SFHs) atz≳ 6 has created new challenges and opportunities for understanding their formation pathways. The degenerate effects of the efficiency and burstiness of star formation on the observed UV luminosity function are separable by galaxy clustering. However, quantifying the timescales of burstiness requires more than just the continuum UV measurements. Here we develop a flexible semi-analytic framework for modeling both the amplitude of star formation rate (SFR) variations and their temporal correlation, from which the luminosity function and clustering can be derived for SFR indicators tracing different characteristic timescales (e.g., UV continuum and Hα luminosities). Based on this framework, we study the prospect of using galaxy summary statistics to distinguish models where SFR fluctuations are prescribed by different power spectral density (PSD) forms. Using the Fisher matrix approach, we forecast the constraints on parameters in our PSD-based model that can be extracted from mock JWST observations of the UV and Hαluminosity functions and clustering bias factors atz∼ 6. If potential confusion due to e.g., dust attenuation and stellar population effects can be properly quantified, these results imply the possibility of probing the burstiness of high-zgalaxies with one-point and two-point statistics and highlight the benefits of combining long-term and short-term SFR tracers. Our flexible framework can be readily extended to characterize the SFH of high-redshift galaxies with a wider range of observational diagnostics.
more »
« less
A Fisher matrix for gravitational-wave population inference
ABSTRACT We derive a Fisher matrix for the parameters characterizing a population of gravitational-wave events. This provides a guide to the precision with which population parameters can be estimated with multiple observations, which becomes increasingly accurate as the number of events and the signal-to-noise ratio of the sampled events increase. The formalism takes into account individual event measurement uncertainties and selection effects, and can be applied to arbitrary population models. We illustrate the framework with two examples: an analytical calculation of the Fisher matrix for the mean and variance of a Gaussian model describing a population affected by selection effects, and an estimation of the precision with which the slope of a power-law distribution of supermassive black hole masses can be measured using extreme-mass-ratio inspiral observations. We compare the Fisher predictions to results from Monte Carlo analyses, finding very good agreement.
more »
« less
- Award ID(s):
- 2207502
- PAR ID:
- 10390221
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 519
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 2736-2753
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract It is well known that the power spectrum is not able to fully characterize the statistical properties of non-Gaussian density fields. Recently, many different statistics have been proposed to extract information from non-Gaussian cosmological fields that perform better than the power spectrum. The Fisher matrix formalism is commonly used to quantify the accuracy with which a given statistic can constrain the value of the cosmological parameters. However, these calculations typically rely on the assumption that the sampling distribution of the considered statistic follows a multivariate Gaussian distribution. In this work, we follow Sellentin & Heavens and use two different statistical tests to identify non-Gaussianities in different statistics such as the power spectrum, bispectrum, marked power spectrum, and wavelet scattering transform (WST). We remove the non-Gaussian components of the different statistics and perform Fisher matrix calculations with theGaussianizedstatistics using Quijote simulations. We show that constraints on the parameters can change by a factor of ∼2 in some cases. We show with simple examples how statistics that do not follow a multivariate Gaussian distribution can achieve artificially tight bounds on the cosmological parameters when using the Fisher matrix formalism. We think that the non-Gaussian tests used in this work represent a powerful tool to quantify the robustness of Fisher matrix calculations and their underlying assumptions. We release the code used to compute the power spectra, bispectra, and WST that can be run on both CPUs and GPUs.more » « less
-
Abstract Gravitational-wave (GW) radiation from a coalescing compact binary is a standard siren, as the luminosity distance of each event can be directly measured from the amplitude of the signal. One possibility to constrain cosmology using the GW siren is to perform statistical inference on a population of binary black hole (BBH) events. In essence, this statistical method can be viewed as follows. We can modify the shape of the distribution of observed BBH events by changing the cosmological parameters until it eventually matches the distribution constructed from an astrophysical population model, thereby allowing us to determine the cosmological parameters. In this work, we derive the Cramér–Rao bound for both cosmological parameters and those governing the astrophysical population model from this statistical dark siren method by examining the Fisher information contained in the event distribution. Our study provides analytical insights and enables fast yet accurate estimations of the statistical accuracy of dark siren cosmology. Furthermore, we consider the bias in cosmology due to unmodeled substructures in the merger rate and mass distribution. We find that a 1% deviation in the astrophysical model can lead to a more than 1% error in the Hubble constant. This could limit the accuracy of dark siren cosmology when there are more than 104BBH events detected.more » « less
-
Abstract Gravitational-wave observations of neutron star mergers can probe the nuclear equation of state by measuring the imprint of the neutron star’s tidal deformability on the signal. We investigate the ability of future gravitational-wave observations to produce a precise measurement of the equation of state from binary neutron star inspirals. Because measurability of the tidal effect depends on the equation of state, we explore several equations of state that span current observational constraints. We generate a population of binary neutron stars as seen by a simulated Advanced LIGO–Virgo network, as well as by a planned Cosmic Explorer observatory. We perform Bayesian inference to measure the parameters of each signal, and we combine measurements across each population to determineR1.4, the radius of a 1.4M⊙neutron star. We find that, with 321 signals, the LIGO–Virgo network is able to measureR1.4to better than 2% precision for all equations of state we consider; however, we also find that achieving this precision could take decades of observation, depending on the equation of state and the merger rate. On the other hand, we find that with one year of observation, Cosmic Explorer will measureR1.4to better than 0.6% precision. In both cases, we find that systematic biases, such as from an incorrect mass prior, can significantly impact measurement accuracy, and efforts will be required to mitigate these effects.more » « less
-
ABSTRACT Line intensity mapping (LIM) experiments probing the nearby Universe can expect a considerable amount of cosmic infrared background (CIB) continuum emission from near and far-infrared galaxies. For the purpose of using LIM to constrain the star formation rate (SFR), we argue that the CIB continuum – traditionally treated as contamination – can be combined with the LIM signal to enhance the SFR constraints achievable. We first present a power spectrum model that combines continuum and line emissions assuming a common SFR model. We subsequently analyse the effectiveness of the joint model in the context of the EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM), which utilizes the $$[{\rm C\, \small {II}}]$$ molecular line to study the SFR. We numerically compute the theoretical power spectra according to our model and the EXCLAIM survey specifics, and perform Fisher analysis to forecast the SFR constraints. We find that although the joint model has no considerable advantage over LIM alone assuming the current survey level of EXCLAIM, its effects become significant when we consider more optimistic values of survey resolution and angular span that are expected of future LIM experiments. We show that the CIB is not only an additional SFR sensitive signal, but also serves to break the SFR parameter degeneracy that naturally emerges from the $$[{\rm C\, \small {II}}]$$ Fisher matrix. For this reason, addition of the CIB will allow improvements in the survey parameters to be better reflected in the SFR constraints, and can be effectively utilized by future LIM experiments.more » « less
