Search for: All records

Dust is a key component of galaxies, but its properties during the earliest eras of structure formation remain elusive. Here we present a simple semi-analytic model of the dust distribution in galaxies atz≳ 5. We calibrate the free parameters of this model to estimates of the UV attenuation (using the IRX-βrelation between infrared emission and the UV spectral slope) and to ALMA measurements of dust emission. We find that the observed dust emission requires that most of the dust expected in these galaxies is retained (assuming a similar yield to lower-redshift sources), but if the dust is spherically distributed, the modest attenuation requires that it be significantly more extended than the stars. Interestingly, the retention fraction is larger for less massive galaxies in our model. However, the required radius is a significant fraction of the host's virial radius and is larger than the estimated extent of dust emission from stacked high-zgalaxies. These can be reconciled if the dust is distributed anisotropically, with typical covering fractions of ∼ 0.2–0.7 in bright galaxies and ≲ 0.1 in fainter ones.

The first wave of observations with JWST has revealed a striking overabundance of luminous galaxies at early times (z> 10) compared to models of galaxies calibrated to pre-JWST data. Early observations have also uncovered a large population of supermassive black holes (SMBHs) atz> 6. Because many of the high-zobjects appear extended, the contribution of active galactic nuclei (AGNs) to the total luminosity has been assumed to be negligible. In this work, we use a semi-empirical model for assigning AGNs to galaxies to show that active galaxies can boost the stellar luminosity function (LF) enough to solve the overabundance problem while simultaneously remaining consistent with the observed morphologies of high-zsources. We construct a model for the composite AGN+galaxy LF by connecting dark matter halo masses to galaxy and SMBH masses and luminosities, accounting for dispersion in the mapping between host galaxy and SMBH mass and luminosity. By calibrating the model parameters — which characterize the M_∙-M_*relation — to a compilation ofz> 10 JWST UVLF data, we show that AGN emission can account for the excess luminosity under a variety of scenarios, including one where 10% of galaxies host BHs of comparable luminosities to their stellar components. Using a sample of simulated objects and real observations, we demonstrate that such low-luminosity AGNs can be `hidden' in their host galaxies and be missed in common morphological analyses. We find that for this explanation to be viable, our model requires a population of BHs that are overmassive (M_∙/M_*~ 10^-2) with respect to their host galaxies compared to the local relation and are more consistent with the observed relation atz= 4-8. We explore the implications of this model for BH seed properties and comment on observational diagnostics necessary to further investigate this explanation.

Early observations with JWST indicate an overabundance of bright galaxies at redshifts z ≳ 10 relative to Hubble-calibrated model predictions. More puzzling still is the apparent lack of evolution in the abundance of such objects between z ∼ 9 and the highest redshifts yet probed, z ∼ 13–17. In this study, we first show that, despite a poor match with JWST luminosity functions (LFs), semi-empirical models calibrated to rest-ultraviolet LFs and colours at 4 ≲ z ≲ 8 are largely consistent with constraints on the properties of individual JWST galaxies, including their stellar masses, ages, and spectral slopes. We then show that order-of-magnitude scatter in the star formation rate of galaxies (at fixed halo mass) can indeed boost the abundance of bright galaxies, provided that star formation is more efficient than expected in low-mass haloes. However, this solution to the abundance problem introduces tension elsewhere: because it relies on the upscattering of low-mass haloes into bright magnitude bins, one expects typical ages, masses, and spectral slopes to be much lower than constraints from galaxies observed thus far. This tension can be alleviated by non-negligible reddening, suggesting that – if the first batch of photometrically selected candidates are confirmed – star formation and dust production could be more efficient than expected in galaxies at z ≳ 10.

Observed scatter in the Lyαopacity of quasar sightlines atz< 6 has motivated measurements of the correlation between Lyαopacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Lyαtroughs atz∼ 5.7 with a deficit of Lyαemitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift toward thez= 6.02 quasar SDSS J1306+0356 and thez= 6.17 quasar PSO J359-06. We find that both fields are underdense in LAEs within 10h⁻¹Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Lyαtroughs. We combine our observations with three additional fields from the literature and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in intergalactic medium (IGM) Lyαopacities. There is tension at low opacities, however, as the models tend to associate lower IGM Lyαopacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent.

The epoch of reionization (EoR) offers a unique window into the dawn of galaxy formation, through which high-redshift galaxies can be studied by observations of both themselves and their impact on the intergalactic medium. Line intensity mapping (LIM) promises to explore cosmic reionization and its driving sources by measuring intensity fluctuations of emission lines tracing the cosmic gas in varying phases. Using LIMFAST, a novel seminumerical tool designed to self-consistently simulate LIM signals of multiple EoR probes, we investigate how building blocks of galaxy formation and evolution theory, such as feedback-regulated star formation and chemical enrichment, might be studied with multitracer LIM during the EoR. On galaxy scales, we show that the star formation law and the feedback associated with star formation can be indicated by both the shape and redshift evolution of LIM power spectra. For a baseline model of metal production that traces star formation, we find that lines highly sensitive to metallicity are generally better probes of galaxy formation models. On larger scales, we demonstrate that inferring ionized bubble sizes from cross-correlations between tracers of ionized and neutral gas requires a detailed understanding of the astrophysics that shape the line luminosity–halo mass relation. Despite various modeling and observational challenges, wide-area, multitracer LIM surveys will provide important high-redshift tests for the fundamentals of galaxy formation theory, especially the interplay between star formation and feedback by accessing statistically the entire low-mass population of galaxies as ideal laboratories, complementary to upcoming surveys of individual sources by new-generation telescopes.

One of the key processes driving galaxy evolution during the Cosmic Dawn is supernova feedback. This likely helps regulate star formation inside of galaxies, but it can also drive winds that influence the large-scale intergalactic medium. Here, we present a simple semi-analytic model of supernova-driven galactic winds and explore the contributions of different phases of galaxy evolution to cosmic metal enrichment in the high-redshift (z ≳ 6) Universe. We show that models calibrated to the observed galaxy luminosity function at z ∼ 6–8 have filling factors $\sim 1{{\%}}$ at z ∼ 6 and $\sim 0.1{{\%}}$ at z ∼ 12, with different star formation prescriptions providing about an order of magnitude uncertainty. Despite the small fraction of space filled by winds, these scenarios predict an upper limit to the abundance of metal-line absorbers in quasar spectra at $z \gtrsim 5$ which is comfortably above that currently observed. We also consider enrichment through winds driven by Pop III star formation in minihalos. We find that these can dominate the total filling factor at $z \gtrsim 10$ and even compete with winds from normal galaxies at z ∼ 6, at least in terms of the total enriched volume. But these regions have much lower overall metallicities, because each one is generated by a small burst of star formation. Finally, we show that Compton cooling of these supernova-driven winds at $z \gtrsim 6$ has only a small effect on the cosmic microwave background.

One of the most exciting advances of the current generation of telescopes has been the detection of galaxies during the epoch of reionization, using deep fields that have pushed these instruments to their limits. It is essential to optimize our analyses of these fields in order to extract as much information as possible from them. In particular, standard methods of measuring the galaxy luminosity function discard information on large-scale dark matter density fluctuations, even though this large-scale structure drives galaxy formation and reionization during the Cosmic Dawn. Measuring these densities would provide a bedrock observable, connecting galaxy surveys to theoretical models of the reionization process and structure formation. Here, we use existing Hubble deep field data to simultaneously fit the universal luminosity function and measure large-scale densities for each Hubble deep field at z = 6–8 by directly incorporating priors on the large-scale density field and galaxy bias. Our fit of the universal luminosity function is consistent with previous methods but differs in the details. For the first time, we measure the underlying densities of the survey fields, including the most over/underdense Hubble fields. We show that the distribution of densities is consistent with current predictions for cosmic variance. This analysis on just 17 fields is a small sample of what will be possible with the James Webb Space Telescope, which will measure hundreds of fields at comparable (or better) depths and at higher redshifts.

Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21-cm power spectrum. One of these parameters, Rmax, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum distance that ionizing photons can travel. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically motivated prescriptions of the mean free path (MFP) of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in seminumerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly defined MFP prescriptions more effectively suppress large-scale structures in the ionization field in seminumerical reionization simulations compared to the standard Rmax approach, and the magnitude of the MFP modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcoming 21-cm power spectrum observations. Finally, we show that in our model, the MFP plays a larger role in regulating the reionization history than in models using Rmax.

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 \gtrsim6$. 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 haloes, with halo masses up to ∼1010.5 M⊙ 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.

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 3D 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 21cm fast 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.