Source author record

Jordan Mirocha

Jordan Mirocha appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

A galaxy-free phenomenological model for the 21-cm power spectrum during reionization

Upper limits from the current generation of interferometers targeting the 21-cm signal from high redshifts have recently begun to rule out physically realistic, though still extreme, models of the Epoch of Reionization (EoR). While inferring the detailed properties of the first galaxies is one of the most important motivations for measuring the high-$z$ 21-cm signal, they can also provide useful constraints on the properties of the intergalactic medium (IGM). Motivated by this, we build a simple, phenomenological model for the 21-cm power spectrum that works directly in terms of IGM properties, which bypasses the computationally expensive 3-D semi-numerical modeling generally employed in inference pipelines and avoids explicit assumptions about galaxy properties. The key simplifying assumptions are that (i) the ionization field is binary, and composed of spherical bubbles with an abundance described well by a parametric bubble size distribution, and (ii) that the spin temperature of the ``bulk'' IGM outside bubbles is uniform. Despite the simplicity of the model, the mean ionized fraction and spin temperature of the IGM recovered from mock 21-cm power spectra generated with \textsc{21cmfast} are generally in good agreement with the true input values. This suggests that it is possible to obtain comparable constraints on the IGM using models with very different assumptions, parameters, and priors. Our approach will thus be complementary to semi-numerical models as upper limits continue to improve in the coming years.

preprint2022arXiv

Constraining Warm Dark Matter and Pop III stars with the Global 21-cm Signal

Upcoming ground and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21-cm signal of the intergalactic medium. In the early universe, when the relative abundance of low-mass DM halos is important, measuring the global signal would place constraints on the damping of structure formation caused by DM having a higher relic velocity (warm dark matter, or WDM) than in cold dark matter (CDM). Such damping, however, can be mimicked by altering the star formation efficiency (SFE) and difficult to detect because of the presence of Pop III stars with unknown properties. We study these various cases and their degeneracies with the WDM mass parameter $m_X$ using a Fisher matrix analysis. We study the $m_X = 7$ keV case and a star-formation model that parametrizes the SFE as a strong function of halo mass and include several variations of this model along with three different input noise levels for the likelihood; we also use a minimum halo virial temperature for collapse near the molecular cooling threshold. We find that when the likelihood includes only Pop II stars, $m_X$ is constrained to an uncertainty of $\sim 0.4$ keV for all models and noise levels at 68$\%$ CI. When the likelihood includes weak Pop III stars, $m_X \sim 0.3$ keV, and if Pop III star formation is relatively efficient, $m_X \sim 0.1$ keV uncertainty, with tight Pop III star-formation parameter constraints. Our results show that the global 21-cm signal is a promising test-bed for WDM models, even in the presence of strong degeneracies with astrophysical parameters.

preprint2022arXiv

Signatures of reionization feedback in the near-infrared background

The reionization of the intergalactic medium at redshifts $z\gtrsim 6$ is expected to have a lasting impact on galaxies residing in low-mass dark matter halos. Unable to accrete or retain gas photo-heated to temperatures $T \gtrsim 10^4$ K, the star formation histories of faint galaxies in the early Universe are expected to decline as they exhaust their gas supply, resulting in a turn-over in the galaxy luminosity function (LF) and a potential solution to the missing satellites problem in the local group. Unfortunately, there are several challenges to constraining `reionization feedback' empirically, most notably that galaxies in low-mass halos are intrinsically faint, and that there are other physical mechanisms capable of inducing a turn-over in the LF. In this work, we investigate a new signature of reionization feedback that is in principle distinct from other processes: as faint galaxies are quenched by reionization, their stellar populations passively age and grow redder while the brighter galaxies nearby continue to form stars at an increasing rate and so remain relatively blue. We find that this contrast, between quenched and un-quenched galaxies induces a scale and colour-dependent signature in the present-day near-infrared background comparable to the expected sensitivity of NASA's upcoming SPHEREx mission. Whereas models with pure mass suppression largely affect the signal at wavelengths $\lesssim 2 μ\rm{m}$, $\sim 5$%-level differences in the background persist out to $\sim 5 μ\rm{m}$ for reionization feedback models. Finally, the power spectra of intensity ratio maps exhibit larger variations, and may thus be a promising target for future analyses.

preprint2021arXiv

Bursty star formation during the Cosmic Dawn driven by delayed stellar feedback

In recent years, several analytic models have demonstrated that simple assumptions about halo growth and feedback-regulated star formation can match the (limited) existing observational data on galaxies at z>6. By extending such models, we demonstrate that imposing a time delay on stellar feedback (as inevitably occurs in the case of supernova explosions) induces burstiness in small galaxies. Although supernova progenitors have short lifetimes (~5-30 Myr), the delay exceeds the dynamical time of galaxies at such high redshifts. As a result, star formation proceeds unimpeded by feedback for several cycles and "overshoots" the expectations of feedback-regulated star formation models. We show that such overshoot is expected even in atomic cooling halos, with masses up to ~10^10.5 Msun at z>6. However, these burst cycles damp out quickly in massive galaxies, because large haloes are more resistant to feedback so retain a continuous gas supply. Bursts in small galaxies - largely beyond the reach of existing observations - induce a scatter in the luminosity of these haloes (of ~1 mag) and increase the time-averaged star formation efficiency by up to an order of magnitude. This kind of burstiness can have substantial effects on the earliest phases of star formation and reionization.

preprint2020arXiv

Effects of self-consistent rest-ultraviolet colours in semi-empirical galaxy formation models

Connecting the observed rest-ultraviolet (UV) luminosities of high-$z$ galaxies to their intrinsic luminosities (and thus star formation rates) requires correcting for the presence of dust. We bypass a common dust-correction approach that uses empirical relationships between infrared (IR) emission and UV colours, and instead augment a semi-empirical model for galaxy formation with a simple -- but self-consistent -- dust model and use it to jointly fit high-$z$ rest-UV luminosity functions (LFs) and colour-magnitude relations ($M_{\mathrm{UV}}$-$β$). In doing so, we find that UV colours evolve with redshift (at fixed UV magnitude), as suggested by observations, even in cases without underlying evolution in dust production, destruction, absorption, or geometry. The observed evolution in our model arises due to the reduction in the mean stellar age and rise in specific star formation rates with increasing $z$. The UV extinction, $A_{\mathrm{UV}}$, evolves similarly with redshift, though we find a systematically shallower relation between $A_{\mathrm{UV}}$ and $M_{\mathrm{UV}}$ than that predicted by IRX-$β$ relationships derived from $z \sim 3$ galaxy samples. Finally, assuming that high $1600 Å$ transmission ($\gtrsim 0.6$) is a reliable LAE indicator, modest scatter in the effective dust surface density of galaxies can explain the evolution both in $M_{\mathrm{UV}}$-$β$ and LAE fractions. These predictions are readily testable by deep surveys with the James Webb Space Telescope.

preprint2020arXiv

High-mass X-ray binaries in nearby metal-poor galaxies: on the contribution to nebular HeII emission

Despite significant progress both observationally and theoretically, the origin of high-ionization nebular HeII emission in galaxies dominated by stellar photoionization remains unclear. Accretion-powered radiation from high-mass X-ray binaries (HMXBs) is still one of the leading proposed explanations for the missing $\mathrm{He^+}$-ionizing photons, but this scenario has yet to be conclusively tested. In this paper, we present nebular line predictions from a grid of photoionization models with input SEDs containing the joint contribution of both stellar atmospheres and a multi-color disk model for HMXBs. This grid demonstrates that HMXBs are inefficient producers of the photons necessary to power HeII, and can only boost this line substantially in galaxies with HMXB populations large enough to power X-ray luminosities of $10^{42}$ erg/s per unit star formation rate (SFR). To test this, we assemble a sample of eleven low-redshift star-forming galaxies with high-quality constraints on both X-ray emission from Chandra and HeII emission from deep optical spectra, including new observations with the MMT. These data reveal that the HMXB populations of these nearby systems are insufficient to account for the observed HeII strengths, with typical X-ray luminosities or upper limits thereon of only $10^{40}$-$10^{41}$ erg/s per SFR. This indicates that HMXBs are not the dominant source of $\mathrm{He^+}$ ionization in these metal-poor star-forming galaxies. We suggest that the solution may instead reside in revisions to stellar wind predictions, softer X-ray sources, or very hot products of binary evolution at low metallicity.

preprint2019arXiv

Canada and the SKA from 2020-2030

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

preprint2019arXiv

Global 21-cm Signal Extraction from Foreground and Instrumental Effects II: Efficient and Self-Consistent Technique for Constraining Nonlinear Signal Models

We present the completion of a data analysis pipeline that self-consistently separates global 21-cm signals from large systematics using a pattern recognition technique. In the first paper of this series, we obtain optimal basis vectors from signal and foreground training sets to linearly fit both components with the minimal number of terms that best extracts the signal given its overlap with the foreground. In this second paper, we utilize the spectral constraints derived in the first paper to calculate the full posterior probability distribution of any signal parameter space of choice. The spectral fit provides the starting point for a Markov Chain Monte Carlo (MCMC) engine that samples the signal without traversing the foreground parameter space. At each MCMC step, we marginalize over the weights of all linear foreground modes and suppress those with unimportant variations by applying priors gleaned from the training set. This method drastically reduces the number of MCMC parameters, augmenting the efficiency of exploration, circumvents the need for selecting a minimal number of foreground modes, and allows the complexity of the foreground model to be greatly increased to simultaneously describe many observed spectra without requiring extra MCMC parameters. Using two nonlinear signal models, one based on EDGES observations and the other on phenomenological frequencies and temperatures of theoretically expected extrema, we demonstrate the success of this methodology by recovering the input parameters from multiple randomly simulated signals at low radio frequencies (10-200 MHz), while rigorously accounting for realistically modeled beam-weighted foregrounds.

preprint2019arXiv

The Effects of Population III X-ray and Radio Backgrounds on the Cosmological 21-cm Signal

We investigate the effects of Population III (Pop III) stars and their remnants on the cosmological 21-cm global signal. By combining a semi-analytic model of Pop III star formation with a global 21-cm simulation code, we investigate how X-ray and radio emission from accreting Pop III black holes may affect both the timing and depth of the 21-cm absorption feature that follows the initial onset of star formation during the Cosmic Dawn. We compare our results to the findings of the EDGES experiment, which has reported the first detection of a cosmic 21-cm signal. In general, we find that our fiducial Pop III models, which have peak star formation rate densities of $\sim 10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ between $z \sim 10$ and $z \sim 15$, are able to match the timing of the EDGES signal quite well, in contrast to models that ignore Pop III stars. To match the unexpectedly large depth of the EDGES signal without recourse to exotic physics, we vary the parameters of emission from accreting black holes (formed as Pop III remnants) including the intrinsic strength of X-ray and radio emission as well as the local column density of neutral gas. We find that models with strong radio emission and relatively weak X-ray emission can self-consistently match the EDGES signal, though this solution requires fine-tuning. We are only able to produce signals with sharp features similar to the EDGES signal if the Pop~III IMF is peaked narrowly around $140 \, M_\odot$.

preprint2016arXiv

The Global 21-cm Signal in the Context of the High-z Galaxy Luminosity Function

Motivated by recent progress in studies of the high-$z$ Universe, we build a new model for the global 21-cm signal that is explicitly calibrated to measurements of the galaxy luminosity function (LF) and further tuned to match the Thomson scattering optical depth of the cosmic microwave background, $τ_e$. Assuming that the $z \lesssim 8$ galaxy population can be smoothly extrapolated to higher redshifts, the recent decline in best-fit values of $τ_e$ and the inefficient heating induced by X-ray binaries (HMXBs; the presumptive sources of the X-ray background at high-$z$) imply that the entirety of cosmic reionization and reheating occurs at redshifts $z \lesssim 12$. In contrast to past global 21-cm models, whose $z \sim 20$ ($ν\sim 70$ MHz) absorption features and strong $\sim 25$ mK emission features were driven largely by the assumption of efficient early star-formation and X-ray heating, our new fiducial model peaks in absorption at $ν\sim 110$ MHz at a depth of $\sim -160$ mK and has a negligible emission component. As a result, a strong emission signal would provide convincing evidence that HMXBs are not the only drivers of cosmic reheating. Shallow absorption troughs should accompany strong heating scenarios, but could also be caused by a low escape fraction of Lyman-Werner photons. Generating signals with troughs at $ν\lesssim 95$ MHz requires a floor in the star-formation efficiency in halos below $\sim 10^{9} M_{\odot}$, which is equivalent to steepening the faint-end of the galaxy LF. These findings demonstrate that the global 21-cm signal is a powerful complement to current and future galaxy surveys and efforts to better understand the interstellar medium in high-$z$ galaxies.

preprint2015arXiv

Interpreting the Global 21-cm Signal from High Redshifts. II. Parameter Estimation for Models of Galaxy Formation

Following our previous work, which related generic features in the sky-averaged (global) 21-cm signal to properties of the intergalactic medium, we now investigate the prospects for constraining a simple galaxy formation model with current and near-future experiments. Markov-Chain Monte Carlo fits to our synthetic dataset, which includes a realistic galactic foreground, a plausible model for the signal, and noise consistent with 100 hours of integration by an ideal instrument, suggest that a simple four-parameter model that links the production rate of Lyman-$α$, Lyman-continuum, and X-ray photons to the growth rate of dark matter halos can be well-constrained (to $\sim 0.1$ dex in each dimension) so long as all three spectral features expected to occur between $40 \lesssim ν/ \mathrm{MHz} \lesssim 120$ are detected. Several important conclusions follow naturally from this basic numerical result, namely that measurements of the global 21-cm signal can in principle (i) identify the characteristic halo mass threshold for star formation at all redshifts $z \gtrsim 15$, (ii) extend $z \lesssim 4$ upper limits on the normalization of the X-ray luminosity star-formation rate ($L_X$-SFR) relation out to $z \sim 20$, and (iii) provide joint constraints on stellar spectra and the escape fraction of ionizing radiation at $z \sim 12$. Though our approach is general, the importance of a broad-band measurement renders our findings most relevant to the proposed Dark Ages Radio Explorer, which will have a clean view of the global 21-cm signal from $\sim 40-120$ MHz from its vantage point above the radio-quiet, ionosphere-free lunar far-side.

preprint2015arXiv

Parametrizations of the global 21-cm signal and parameter estimation from single-dipole experiments

One approach to extracting the global 21-cm signal from total-power measurements at low radio frequencies is to parametrize the different contributions to the data and then fit for these parameters. We examine parametrizations of the 21-cm signal itself, and propose one based on modelling the Lyman-alpha background, IGM temperature and hydrogen ionized fraction using tanh functions. This captures the shape of the signal from a physical modelling code better than an earlier parametrization based on interpolating between maxima and minima of the signal, and imposes a greater level of physical plausibility. This allows less biased constraints on the turning points of the signal, even though these are not explicitly fit for. Biases can also be alleviated by discarding information which is less robustly described by the parametrization, for example by ignoring detailed shape information coming from the covariances between turning points or from the high-frequency parts of the signal, or by marginalizing over the high-frequency parts of the signal by fitting a more complex foreground model. The fits are sufficiently accurate to be usable for experiments gathering 1000 h of data, though in this case it may be important to choose observing windows which do not include the brightest areas of the foregrounds. Our assumption of pointed, single-antenna observations and very broad-band fitting makes these results particularly applicable to experiments such as the Dark Ages Radio Explorer, which would study the global 21-cm signal from the clean environment of a low lunar orbit, taking data from the far side.

preprint2014arXiv

Decoding the X-ray Properties of Pre-Reionization Era Sources

Evolution in the X-ray luminosity -- star formation rate ($L_X$-SFR) relation could provide the first evidence of a top-heavy stellar initial mass function in the early universe, as the abundance of high-mass stars and binary systems are both expected to increase with decreasing metallicity. The sky-averaged (global) 21-cm signal has the potential to test this prediction via constraints on the thermal history of the intergalactic medium, since X-rays can most easily escape galaxies and heat gas on large scales. A significant complication in the interpretation of upcoming 21-cm measurements is the unknown spectrum of accreting black holes at high-$z$, which depends on the mass of accreting objects and poorly constrained processes such as how accretion disk photons are processed by the disk atmosphere and host galaxy interstellar medium. Using a novel approach to solving the cosmological radiative transfer equation (RTE), we show that reasonable changes in the characteristic BH mass affects the amplitude of the 21-cm signal's minimum at the $\sim 10-20$ mK level --- comparable to errors induced by commonly used approximations to the RTE --- while modifications to the intrinsic disk spectrum due to Compton scattering (bound-free absorption) can shift the position of the minimum of the global signal by $Δz \approx 0.5$ ($Δz \approx 2$), and modify its amplitude by up to $\approx 10$ mK ($\approx 50$ mK) for a given accretion history. Such deviations are larger than the uncertainties expected of current global 21-cm signal extraction algorithms, and could easily be confused with evolution in the $L_X$-SFR relation.

preprint2013arXiv

Interpreting the Global 21-cm Signal from High Redshifts. I. Model Independent Constraints

The sky-averaged (global) 21-cm signal is a powerful probe of the intergalactic medium (IGM) prior to the completion of reionization. However, it has so far been unclear that even in the best case scenario, in which the signal is accurately extracted from the foregrounds, that it will provide more than crude estimates of when the Universe's first stars and black holes form. In contrast to previous work, which has focused on predicting the 21-cm signatures of the first luminous objects, we investigate an arbitrary realization of the signal, and attempt to translate its features to the physical properties of the IGM. Within a simplified global framework, the 21-cm signal yields quantitative constraints on the Lyman-alpha background intensity, net heat deposition, ionized fraction, and their time derivatives, without invoking models for the astrophysical sources themselves. The 21-cm absorption signal is most easily interpreted, setting strong limits on the heating rate density of the Universe with a measurement of its redshift alone, independent of the ionization history or details of the Lyman-alpha background evolution. In a companion paper we extend these results, focusing on the confidence with which one can infer source emissivities from IGM properties.

preprint2012arXiv

Optimized Multi-Frequency Spectra for Applications in Radiative Feedback and Cosmological Reionization

The recent implementation of radiative transfer algorithms in numerous hydrodynamics codes has led to a dramatic improvement in studies of feedback in various astrophysical environments. However, because of methodological limitations and computational expense, the spectra of radiation sources are generally sampled at only a few evenly-spaced discrete emission frequencies. Using one-dimensional radiative transfer calculations, we investigate the discrepancies in gas properties surrounding model stars and accreting black holes that arise solely due to spectral discretization. We find that even in the idealized case of a static and uniform density field, commonly used discretization schemes induce errors in the neutral fraction and temperature by factors of two to three on average, and by over an order of magnitude in certain column density regimes. The consequences are most severe for radiative feedback operating on large scales, dense clumps of gas, and media consisting of multiple chemical species. We have developed a method for optimally constructing discrete spectra, and show that for two test cases of interest, carefully chosen four-bin spectra can eliminate errors associated with frequency resolution to high precision. Applying these findings to a fully three-dimensional radiation-hydrodynamic simulation of the early universe, we find that the HII region around a primordial star is substantially altered in both size and morphology, corroborating the one-dimensional prediction that discrete spectral energy distributions can lead to sizable inaccuracies in the physical properties of a medium, and as a result, the subsequent evolution and observable signatures of objects embedded within it.