Early galaxy formation in warm dark matter cosmologies

We present a framework for high-redshift ($z \geq 7$) galaxy formation that traces their dark matter (DM) and baryonic assembly in four cosmologies: Cold Dark Matter (CDM) and Warm Dark Matter (WDM) with particle masses of $m_x =$ 1.5, 3 and 5 ${\rm keV}$. We use the same astrophysical parameters regulating star formation and feedback, chosen to match current observations of the evolving ultra violet luminosity function (UV LF). We find that the assembly of observable (with current and upcoming instruments) galaxies in CDM and $m_x \geq 3 {\rm keV}$ WDM results in similar halo mass to light ratios (M/L), stellar mass densities (SMDs) and UV LFs. However the suppression of small-scale structure leads to a notably delayed and subsequently more rapid stellar assembly in the $1.5 {\rm keV}$ WDM model. Thus galaxy assembly in $m_x \leq 2 {\rm keV}$ WDM cosmologies is characterized by: (i) a dearth of small-mass halos hosting faint galaxies; and (ii) a younger, more UV bright stellar population, for a given stellar mass. The higher M/L ratio (effect ii) partially compensates for the dearth of small-mass halos (effect i), making the resulting UV LFs closer to CDM than expected from simple estimates of halo abundances. We find that the redshift evolution of the SMD is a powerful probe of the nature of DM. Integrating down to a limit of $M_{UV} =-16.5$ for the James Webb Space Telescope (JWST), the SMD evolves as $\log$(SMD)$\propto -0.63 (1+z)$ in $m_x = 1.5 {\rm keV}$ WDM, as compared to $\log$(SMD)$\propto -0.44 (1+z)$ in CDM. Thus high-redshift stellar assembly provides a powerful testbed for WDM models, accessible with the upcoming JWST.Comment: Accepted for publication in Ap

Clustering and lifetime of Lyman Alpha Emitters in the Epoch of Reionization

We calculate Lyman Alpha Emitter (LAE) angular correlation functions (ACFs) at $z \simeq 6.6$ and the fraction of lifetime (for the 100 Myrs preceding $z\simeq6.6$) galaxies spend as Lyman Break Galaxies (LBGs) or as LBGs with Lyman Alpha (Ly$\alpha$) emission using a model that combines SPH cosmological simulations (GADGET-2), dust attenuation and a radiative transfer code (pCRASH). The ACFs are a powerful tool that significantly narrows the 3D parameter space allowed by LAE Ly$\alpha$ and UV luminosity functions (LFs) alone. With this work, we simultaneously constrain the escape fraction of ionizing photons $f_{esc}=0.05-0.5$, the mean fraction of neutral hydrogen in the intergalactic medium (IGM) $\langle \chi_{HI} \rangle \leq 0.01$ and the dust-dependent ratio of the escape fractions of Ly$\alpha$ and UV continuum photons $f_{\alpha}/f_c=0.6-1.2$. Our results show that reionization has the largest impact on the amplitude of the ACFs, and its imprints are clearly distinguishable from those of $f_{esc}$ and $f_\alpha/f_c$. We also show that galaxies with a critical stellar mass of $M_* = 10^{8.5} (10^{9.5})M_{\odot}$ produce enough luminosity to stay visible as LBGs (LAEs). Finally, the fraction of time during the past 100 Myrs prior to z=6.6 a galaxy spends as a LBG or as a LBG with Ly$\alpha$ emission increases with the UV magnitude (and the stellar mass $M_*$): considering observed (dust and IGM attenuated) luminosities, the fraction of time a galaxy spends as a LBG (LAE) increases from 65% to 100% (0-100%) as $M_{UV}$ decreases from $M_{UV} = -18.0$ to $-23.5$ ($M_*$ increases from $10^8-10^{10.5} M_{\odot}$). Thus in our model the brightest (most massive) LBGs most often show Ly$\alpha$ emission.Comment: 11 pages, 7 figures, accepted for publication in MNRAS, comments welcom

The habitability of the Universe through 13 billion years of cosmic time

The field of astrobiology has made tremendous progress in modelling galactic-scale habitable zones which offer a stable environment for life to form and evolve in complexity. Recently, this idea has been extended to cosmological scales by studies modelling the habitability of the local Universe in its entirety (e.g. Dayal et al. 2015; Li & Zhang 2015). However, all of these studies have solely focused on estimating the potentially detrimental effects of either Type II supernovae (SNII) or Gamma Ray Bursts (GRBs), ignoring the contributions from Type Ia supernovae (SNIa) and active galactic nuclei (AGN). In this study we follow two different approaches, based on (i) the amplitude of deleterious radiation and (ii) the total planet-hosting volume irradiated by deleterious radiation. We simultaneously track the contributions from the key astrophysical sources (SNII, SNIa, AGN and GRBs) for the entire Universe, for both scenarios, to determine its habitability through 13.8 billion years of cosmic time. We find that SNII dominate the total radiation budget and the volume irradiated by deleterious radiation at any cosmic epoch closely followed by SNIa (that contribute half as much as SNII), with GRBs and AGN making up a negligible portion (<1%). Secondly, as a result of the total mass in stars (or the total number of planets) slowly building-up with time and the total deleterious radiation density, and volume affected, falling-off after the first 3 billion years, we find that the Universe has steadily increased in habitability through cosmic time. We find that, depending on the exact model assumptions, the Universe is 2.5 to 20 times more habitable today compared to when life first appeared on the Earth 4 billion years ago. We find that this increase in habitability will persist until the final stars die out over the next hundreds of billions of years.Comment: Under refereeing in Ap

Coevolution of metallicity and star formation in galaxies to z=3.7: I. A fundamental plane

With the aim of understanding the coevolution of star formation rate (SFR), stellar mass (M*), and oxygen abundance (O/H) in galaxies up to redshift z=3.7, we have compiled the largest available dataset for studying Metallicity Evolution and Galaxy Assembly (MEGA); it comprises roughly 1000 galaxies with a common O/H calibration and spans almost two orders of magnitude in metallicity, a factor of 10^6 in SFR, and a factor of 10^5 in stellar mass. From a Principal Component Analysis, we find that the 3-dimensional parameter space reduces to a Fundamental Plane of Metallicity (FPZ) given by 12+log(O/H) = -0.14 log (SFR) + 0.37 log (M*) + 4.82. The mean O/H FPZ residuals are small (0.16 dex) and consistent with trends found in smaller galaxy samples with more limited ranges in M*, SFR, and O/H. Importantly, the FPZ is found to be redshift-invariant within the uncertainties. In a companion paper, these results are interpreted with an updated version of the model presented by Dayal et al. (2013).Comment: 19 pages, 10 figures, 4 tables, accepted for publication in MNRA

Constraining dust formation in high-redshift young galaxies

Core-collapse supernovae (SNe) are believed to be the first significant source of dust in the Universe. Such SNe are expected to be the main dust producers in young high-redshift Lyman $\alpha$ emitters (LAEs) given their young ages, providing an excellent testbed of SN dust formation models during the early stages of galaxy evolution. We focus on the dust enrichment of a specific, luminous LAE (Himiko, $z\simeq 6.6$) for which a stringent upper limit of $52.1~\mu$Jy ($3\sigma$) has recently been obtained from ALMA continuum observations at 1.2 mm. We predict its submillimetre dust emission using detailed models that follow SN dust enrichment and destruction and the equilibrium dust temperature, and obtain a plausible upper limit to the dust mass produced by a single SN: $m_\mathrm{d,SN} < 0.15$--0.45 M$_\odot$, depending on the adopted dust optical properties. These upper limits are smaller than the dust mass deduced for SN 1987A and that predicted by dust condensation theories, implying that dust produced in SNe are likely to be subject to reverse shock destruction before being injected into the interstellar medium. Finally, we provide a recipe for deriving $m_\mathrm{d,SN}$ from submillimetre observations of young, metal poor objects wherein condensation in SN ejecta is the dominant dust formation channel.Comment: 10 pages, 3 figures, accepted for publication in MNRA

Warm dark matter constraints from high-zz Direct Collapse Black Holes using the JWST

We use a semi-analytic model, ${\it Delphi}$, that jointly tracks the dark matter and baryonic assembly of high-redshift ($z \simeq 4-20$) galaxies to gain insight on the number density of Direct Collapse Black Hole (DCBH) hosts in three different cosmologies: the standard Cold Dark Matter (CDM) model and two Warm Dark Matter (WDM) models with particle masses of 3.5 and 1.5 keV. Obtaining the Lyman-Werner (LW) luminosity of each galaxy from ${\it Delphi}$, we use a clustering bias analysis to identify all, pristine halos with a virial temperature $T_{vir}>=10^4$ K that are irradiated by a LW background above a critical value as, DCBH hosts. In good agreement with previous studies, we find the DCBH number density rises from $\sim10^{-6.1}$ to $\sim 10^{-3.5}\, \mathrm{cMpc^{-3}}$ from $z\simeq 17.5$ to $8$ in the CDM model using a critical LW background value of $30 J_{21}$ (where $J_{21}= 10^{-21} \, {\rm erg\, s^{-1}\, Hz^{-1} \, cm^{-2} \, sr^{-1}}$). We find that a combination of delayed structure formation and an accelerated assembly of galaxies results in a later metal-enrichment and an accelerated build-up of the LW background in the 1.5 keV WDM model, resulting in DCBH hosts persisting down to much lower redshifts ($z \simeq 5$) as compared to CDM where DCBH hosts only exist down to $z \simeq 8$. We end by showing how the expected colours in three different bands of the Near Infrared Camera (NIRCam) onboard the forthcoming James Webb Space Telescope (${\it JWST}$) can be used to hunt for potential $z \simeq 5-9$ DCBHs, allowing hints on the WDM particle mass.Comment: Accepted to MNRAS with minor revisio

Reionization and Galaxy Formation in Warm Dark Matter Cosmologies

We compare model results from a semi-analytic (merger-tree based) framework for high-redshift (z ~ 5-20) galaxy formation against reionization indicators, including the Planck electron scattering optical depth and the ionizing photon emissivity, to shed light on the reionization history and sources in Cold (CDM) and Warm Dark Matter (WDM; particle masses of $m_x = 1.5,3$ and 5 keV) cosmologies. This model includes all the key processes of star formation, supernova feedback, the merger/accretion/ejection driven evolution of gas and stellar mass and the effect of the ultra-violet background (UVB) in photo-evaporating the gas content of low-mass galaxies. We find that the delay in the start of reionization in light (1.5 keV) WDM models can be compensated by a steeper redshift evolution of the ionizing photon escape fraction and a faster mass assembly, resulting in reionization ending at comparable redshifts (z~5.5) in all the DM models considered. We find the bulk of the reionization photons come from galaxies with a halo mass $M_h < 10^9 M_\odot$ and a UV magnitude $-15 < M_{UV} < -10$ in CDM. The progressive suppression of low-mass halos with decreasing $m_x$ leads to a shift in the reionization population to larger halo masses of $M_h > 10^9 M_\odot$ and $-17 < M_{UV} < -13$ for 1.5 keV WDM. We find that current observations of the electron scattering optical depth and the Ultra-violet luminosity function are equally compatible with all the (cold and warm) DM models considered in this work. We propose that global indicators including the redshift evolution of the stellar mass density and the stellar mass-halo mass relation, observable with the James Webb Space Telescope, can be used to distinguish between CDM and WDM (1.5 keV) cosmologies.Comment: Accepted to Ap

Simulating the assembly of galaxies at redshifts z = 6 - 12

We use state-of-the-art simulations to explore the physical evolution of galaxies in the first billion years of cosmic time. First, we demonstrate that our model reproduces the basic statistical properties of the observed Lyman-break galaxy (LBG) population at z = 6 - 8, including the evolving ultra-violet (UV) luminosity function (LF), the stellar-mass density (SMD), and the average specific star-formation rates (sSFR) of LBGs with M_{UV} < -18 (AB mag). Encouraged by this success we present predictions for the behaviour of fainter LBGs extending down to M_{UV} <= -15 (as will be probed with the James Webb Space Telescope) and have interrogated our simulations to try to gain insight into the physical drivers of the observed population evolution. We find that mass growth due to star formation in the mass-dominant progenitor builds up about 90% of the total z ~ 6 LBG stellar mass, dominating over the mass contributed by merging throughout this era. Our simulation suggests that the apparent "luminosity evolution" depends on the luminosity range probed: the steady brightening of the bright end of the LF is driven primarily by genuine physical luminosity evolution and arises due to a fairly steady increase in the UV luminosity (and hence star-formation rates) in the most massive LBGs. However, at fainter luminosities the situation is more complex, due in part to the more stochastic star-formation histories of lower-mass objects; at this end, the evolution of the UV LF involves a mix of positive and negative luminosity evolution (as low-mass galaxies temporarily brighten then fade) coupled with both positive and negative density evolution (as new low-mass galaxies form, and other low-mass galaxies are consumed by merging). We also predict the average sSFR of LBGs should rise from sSFR = 4.5 Gyr^-1 at z = 6 to about 11 Gyr^-1 by z = 9.Comment: Accepted for publication in MNRA