Source author record

Pratika Dayal

Pratika Dayal 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

33works
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

33 published item(s)

preprint2026arXiv

Probing General Relativity on Cosmological Scales in the 2040s

General relativity is exquisitely tested in strong-field regimes, yet its validity on cosmological scales remains largely unexplored. Upcoming wide and deep large-scale structure surveys will access the ultra-large, linear scales where relativistic effects - Doppler terms, gravitational redshift, lensing magnification, and potential evolution - leave significant imprints in the clustering of galaxies. These signatures represent unique probes of spacetime that are inaccessible to standard Newtonian analyses but increasingly important as survey volumes grow. We outline the scientific potential of next-generation facilities, such as those envisioned within ESO's Expanding Horizons programme, to deliver the first robust measurements of relativistic effects in large-scale structure through multi-tracer power spectra and the single-tracer bispectrum of high-redshift Lyman-break galaxies. Detecting these contributions would open a new window on gravity, enabling precision tests of general relativity and its alternatives on cosmological scales in the 2040s.

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2023arXiv

UNCOVER: Illuminating the Early Universe -- JWST/NIRSpec Confirmation of $z > 12$ Galaxies

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies.

preprint2022arXiv

A Search for H-Dropout Lyman Break Galaxies at z~12-16

We present two bright galaxy candidates at z~12-13 identified in our H-dropout Lyman break selection with 2.3 deg2 near-infrared deep imaging data. These galaxy candidates, selected after careful screening of foreground interlopers, have spectral energy distributions showing a sharp discontinuity around 1.7 um, a flat continuum at 2-5 um, and non-detections at <1.2 um in the available photometric datasets, all of which are consistent with z>12 galaxy. An ALMA program targeting one of the candidates shows a tentative 4sigma [OIII]88um line at z=13.27, in agreement with its photometric redshift estimate. The number density of the z~12-13 candidates is comparable to that of bright z~10 galaxies, and is consistent with a recently proposed double power-law luminosity function rather than the Schechter function, indicating little evolution in the abundance of bright galaxies from z~4 to 13. Comparisons with theoretical models show that the models cannot reproduce the bright end of rest-frame ultraviolet luminosity functions at z~10-13. Combined with recent studies reporting similarly bright galaxies at z~9-11 and mature stellar populations at z~6-9, our results indicate the existence of a number of star-forming galaxies at z>10, which will be detected with upcoming space missions such as James Webb Space Telescope, Nancy Grace Roman Space Telescope, and GREX-PLUS.

preprint2022arXiv

Are the Newly-Discovered $z \sim 13$ Drop-out Sources Starburst Galaxies or Quasars?

The detection of two $z\sim 13$ galaxy candidates has opened a new window on galaxy formation at an era only $330$ Myr after the Big Bang. Here, we investigate the physical nature of these sources: are we witnessing star forming galaxies or quasars at such early epochs? If powered by star formation, the observed ultraviolet (UV) luminosities and number densities can be jointly explained if: (i) these galaxies are extreme star-formers with star formation rates $5-24\times$ higher than those expected from extrapolations of average lower-redshift relations; (ii) the star formation efficiency increases with halo mass and is countered by increasing dust attenuation from $z \sim 10-5$; (iii) they form stars with an extremely top-heavy initial mass function. The quasar hypothesis is also plausible, with the UV luminosity produced by black holes of $\sim 10^8 \, \mathrm{M_{\odot}}$ accreting at or slightly above the Eddington rate ($f_{\rm Edd}\sim 1.0$). This black hole mass at $z\sim13$ would require very challenging, but not implausible, growth parameters. If spectroscopically confirmed, these two sources will represent a remarkable laboratory to study the Universe at previously inaccessible redshifts.

preprint2022arXiv

Astraeus VII: The environmental-dependent assembly of galaxies in the Epoch of Reionization

Using the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) framework, we explore the impact of environmental density and radiative feedback on the assembly of galaxies and their host halos during the Epoch of Reionization. The ASTRAEUS framework allows us to study the evolution of galaxies with masses ($\rm 10^{8.2}M_\odot < M_{\rm h} < 10^{13}M_\odot$) in wide variety of environment ($-0.5 < {\rm log}(1+δ) < 1.3$ averaged over $(2~{\rm cMpc})^3$). We find that : (i) there exists a mass- and redshift- dependent "characteristic" environment (${\rm log} (1+δ_a(M_{\rm h}, z)) = 0.021\times (M_{\rm h}/M_\odot)^{0.16} + 0.07 z -1.12$, up to $z\sim 10$) at which galaxies are most efficient at accreting dark matter, e.g at a rate of $0.2\%$ of their mass every Myr at $z=5$; (ii) the number of minor and major mergers and their contributions to the dark matter assembly increases with halo mass at all redshifts and is mostly independent of the environment; (iii) at $z=5$ minor mergers contribute slightly more (by up to $\sim 10\%$) to the dark matter assembly while for the stellar assembly, major mergers dominate the contribution from minor mergers for $M_{\rm h}\lesssim 10^{11.5}M_\odot$ galaxies; (iv) radiative feedback quenches star formation more in low-mass galaxies ($M_{\rm h} \lesssim 10^{9.5}M_\odot$) in over-dense environments (${\rm log}(1+δ) > 0.5$); dominated by their major branch, this yields star formation histories biased towards older ages with a slower redshift evolution.

preprint2022arXiv

Schrodinger's Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $Λ$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2μ$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($Δz\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical.

preprint2022arXiv

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates.

preprint2020arXiv

A hint on the metal-free star formation rate density from 21cm-EDGES data

We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dotρ_{*}^{III}$ required at high-redshifts ($z \gtrsim 16$) in order to reproduce both the amplitude and the redshift of the EDGES 21-cm global signal. Our model accounts for the Lyman Alpha (Ly$α$), radio and X-ray backgrounds from both Pop III and metal-enriched Population II (Pop II) stars. For the latter, we use the star formation rate density estimates (and the Ly$α$ background) from the {\it Delphi} semi-analytic model that has been shown to reproduce all key observables for galaxies at $z \gtrsim 5$; the radio and X-ray backgrounds are fixed using low-$z$ values. The constraints on the free parameters characterizing the properties of the Pop III stars are obtained using a Markov Chain Monte Carlo analysis. Our results yield a $\dotρ_{*}^{III}$ that whilst increasing from $z \sim 21-16$ thereafter shows a sharp decline which is in excellent agreement with the results found by \citet{valiante2016} to simulate the growth of $z \sim 6 - 7$ quasars and their host galaxies, suggesting that the bulk of Pop III star formation occurs in the rarest and most massive metal-poor halos at $z \lesssim 20$. This allows Pop III stars to produce a rapidly growing Ly$α$ background between $z \sim 21-15$. Further, Pop III stars are required to provide a radio background that is about $3-4$ orders of magnitude higher than that provided by Pop II stars although Pop II stars dominate the X-ray background.

preprint2020arXiv

Bubble mapping with the Square Kilometer Array -- I. Detecting galaxies with Euclid, JWST, WFIRST and ELT within ionized bubbles in the intergalactic medium at z>6

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

preprint2020arXiv

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

preprint2019arXiv

The 21cm bispectrum during reionization: a tracer of the ionization topology

We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.

preprint2018arXiv

Early galaxy formation and its large-scale effects

Galaxy formation is at the heart of our understanding of cosmic evolution. Although there is a consensus that galaxies emerged from the expanding matter background by gravitational instability of primordial fluctuations, a number of additional physical processes must be understood and implemented in theoretical models before these can be reliably used to interpret observations. In parallel, the astonishing recent progresses made in detecting galaxies that formed only a few hundreds of million years after the Big Bang is pushing the quest for more sophisticated and detailed studies of early structures. In this review, we combine the information gleaned from different theoretical models/studies to build a coherent picture of the Universe in its early stages which includes the physics of galaxy formation along with the impact that early structures had on large-scale processes as cosmic reionization and metal enrichment of the intergalactic medium.

preprint2016arXiv

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).

preprint2016arXiv

Coevolution of metallicity and star formation in galaxies to z=3.7: II. A theoretical model

Recent work suggests that galaxy evolution, and the build-up of stellar mass (M*) over cosmic time, is characterized by changes with redshift of star formation rate (SFR) and oxygen abundance (O/H). In a companion paper, we have compiled a large dataset to study Metallicity Evolution and Galaxy Assembly (MEGA), consisting of roughly 1000 galaxies to z=3.7 with a common O/H calibration. Here we interpret the MEGA scaling relations of M*, SFR, and O/H with an updated version of the model presented by Dayal et al. (2013). This model successfully reproduces the observed O/H ratio of 80,000 galaxies selected from the Sloan Digital Sky Survey to within 0.05-0.06 dex. By extending the model to the higher redshift MEGA sample, we find that although the specific mass loading of outflows does not change measurably during the evolution, the accretion rate and gas content of galaxies increase significantly with redshift. These two effects can explain, either separately or possibly in tandem, the observed lower metal abundance of high-z galaxies.

preprint2016arXiv

Interpreting the evolution of galaxy colours from $z = 8$ to $z = 5$

We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, $β$, of galaxies at $5 \leq z \leq 8$, to improve our understanding of their dust content and ISM properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess - UV color relation we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. (1999) one, usually adopted to correct the observed UV luminosities of high-$z$ galaxies for the effects of dust extinction. As a result, their total star formation rates might have been over-estimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-$z$ galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at $z \lesssim 7$.

preprint2016arXiv

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.

preprint2015arXiv

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$α$) 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$α$ 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 χ_{HI} \rangle \leq 0.01$ and the dust-dependent ratio of the escape fractions of Ly$α$ and UV continuum photons $f_α/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_α/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$α$ 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$α$ emission.

preprint2015arXiv

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.

preprint2015arXiv

The quest for cradles of life: using the fundamental metallicity relation to hunt for the most habitable type of galaxy

The field of astrobiology has made huge strides in understanding the habitable zones around stars (Stellar Habitable Zones) where life can begin, sustain its existence and evolve into complex forms. A few studies have extended this idea by modelling galactic-scale habitable zones (Galactic Habitable Zones) for our Milky Way and specific elliptical galaxies. However, estimating the habitability for galaxies spanning a wide range of physical properties has so far remained an outstanding issue. Here, we present a "cosmobiological" framework that allows us to sift through the entire galaxy population in the local Universe and answer the question "Which type of galaxy is most likely to host complex life in the cosmos"? Interestingly, the three key astrophysical criteria governing habitability (total mass in stars, total metal mass and ongoing star formation rate) are found to be intricately linked through the "fundamental metallicity relation" as shown by SDSS (Sloan Digital Sky Survey) observations of more than a hundred thousand galaxies in the local Universe. Using this relation we show that metal-rich, shapeless giant elliptical galaxies at least twice as massive as the Milky Way (with a tenth of its star formation rate) can potentially host ten thousand times as many habitable (earth-like) planets, making them the most probable "cradles of life" in the Universe.

preprint2014arXiv

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 $α$ 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~μ$Jy ($3σ$) 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.

preprint2014arXiv

Essential physics of early galaxy formation

We present a theoretical model embedding the essential physics of early galaxy formation (z = 5-12) based on the single premise that any galaxy can form stars with a maximal limiting efficiency that provides enough energy to expel all the remaining gas, quenching further star formation. This simple idea is implemented into a merger-tree based semi-analytical model that utilises two mass and redshift-independent parameters to capture the key physics of supernova feedback in ejecting gas from low-mass halos, and tracks the resulting impact on the subsequent growth of more massive systems via halo mergers and gas accretion. Our model shows that: (i) the smallest halos (halo mass $M_h \leq 10^{10} M_\odot$) build up their gas mass by accretion from the intergalactic medium; (ii) the bulk of the gas powering star formation in larger halos ($M_h \geq 10^{11.5} M_\odot$) is brought in by merging progenitors; (iii) the faint-end UV luminosity function slope evolves according to $α= -1.75 \log \,z -0.52$. In addition, (iv) the stellar mass-to-light ratio is well fit by the functional form $\log\, M_* = -0.38 M_{UV} -0.13\, z + 2.4$, which we use to build the evolving stellar mass function to compare to observations. We end with a census of the cosmic stellar mass density (SMD) across galaxies with UV magnitudes over the range $-23 \leq M_{UV} \leq -11$ spanning redshifts $5 < z < 12$: (v) while currently detected LBGs contain $\approx 50$% (10%) of the total SMD at $z=5$ (8), the JWST will detect up to 25% of the SMD at $z \simeq 9.5$.

preprint2013arXiv

CLUES to the past: Local Group progenitors amongst high-redshift Lyman Break Galaxies

We use state-of-the-art numerical simulations to explore the observability and the expected physical properties of the progenitors of the local group galaxies at z=6-8, within 1 billion years of the big bang. We find that the most massive progenitors of the Milky Way (MW) and Andromeda (M31) at z=6 and 7 are predicted to have absolute UV continuum magnitudes between -17 and -18, suggesting that their analogues lie close to the detection limits of the deepest near-infrared surveys conducted to date (i.e. HST WFC3/IR UDF12). This in turn confirms that the majority of currently known z=6-8 galaxies are expected to be the seeds of present-day galaxies which are more massive than L* spirals. We also discuss the properties of the local-group progenitors at these early epochs, extending down to absolute magnitudes M_UV = -13. The most massive MW/M31 progenitors at z=7 have stellar masses of about 10^7.5-8 solar masses, stellar metallicities between 3 and 6% of the solar value, and predicted observed UV continuum slopes between -2.4 and -2.5.

preprint2013arXiv

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.

preprint2013arXiv

The physics of the fundamental metallicity relation

We present a simple, redshift-independent analytic model that explains the local Fundamental Metallicity Relation (FMR), taking into account the physical processes of star formation, inflow of metal-poor intergalactic medium (IGM) gas, and the outflow of metal rich interstellar medium (ISM) gas. We show that the physics of the FMR can be summarised as follows: for massive galaxies with stellar mass M_* >= 10^11 solar masses, ISM metal enrichment due to star formation is compensated by inflow of metal poor IGM gas, leading to a constant value of the gas metallicity with star formation rate (SFR); outflows are rendered negligible as a result of the large potential wells of these galaxies. On the other hand, as a result of their smaller SFR, less massive galaxies produce less heavy elements that are also more efficiently ejected due to their shallow potential wells; as a result, for a given M_*, the gas metallicity decreases with SFR. For such galaxies, the outflow efficiency determines both the slope, and the knee of the metallicity-SFR relation. Without changing any parameters, this simple model is also successfully matched to the gas fraction - gas metallicity relation observed for a sample of about 260 nearby galaxies.

preprint2012arXiv

Lyman Alpha Emitters and Lyman Break Galaxies: dichotomous twins?

(Abridged) We extend our previous studies aimed at modeling Lyman Alpha Emitters (LAEs) to the second population of high redshift sources, Lyman Break Galaxies (LBGs), with the final goal of investigating the physical relationship between them at 6 < z <8. We use a set of large (~10^6 Mpc^3) cosmological SPH simulations that include a detailed treatment of star formation, feedback, metal enrichment and supernova dust production; these same simulations have already been shown to successfully reproduce a large number of observed properties of LAEs (Dayal et al. 2010). We identify LBGs as galaxies with an absolute ultraviolet (UV) magnitude M_{UV} <= -17, consistent with current observational criteria. With no further parameter tuning, the model reproduces the redshift evolution of the LBG UV luminosity function, stellar mass function, and star formation rate (SFR) density extremely well, compared to the data at z~6-8. It predicts a z-independent average beta value of about -2.2, consistent with the most recent observational results. The mean LBG specific SFR increases from 6.7 Gyr^-1 at z ~ 6 to 13.9 Gyr^-1 at z ~ 8, and is largest for the smallest LBGs. From a comparison of the simulated LAE and LBG populations, we find no appreciable differences between them in terms of the stellar masses, ages, SFR, and dust content; only the faintest LBGs with M_{UV} >= -18 (-19) at z ~ 6 (8) do not show an observable Lyman Alpha line. LAEs hence represent a luminous LBG subset, whose relative extent depends only on the adopted equivalent width (EW) selection criteria. Further, in line with current observations, we find that the fraction of LBGs showing a Lyman Alpha line decreases with increasing redshift due to the combined effects of dust and reionization. We conclude with a brief critical model discussion, which emphasizes the uncertainties inherent to theoretical EW determinations.

preprint2012arXiv

Molecular hydrogen in Lyman Alpha Emitters

We present a physically motivated model to estimate the molecular hydrogen (H2) content of high-redshift (z~5.7,6.6) Lyman Alpha Emitters (LAEs) extracted from a suite of cosmological simulations. We find that the H2 mass fraction, (f_H2), depends on three main LAE physical properties: (a) star formation rate, (b) dust mass, and (c) cold neutral gas mass. At z~5.7, the value of f_H2 peaks and ranges between 0.5-0.9 for intermediate mass LAEs with stellar mass M_* ~ 10^{9-10} solar mass, decreasing for both smaller and larger galaxies. However, the largest value of the H2 mass is found in the most luminous LAEs. These trends also hold at z\sim6.6, although, due to a lower dust content, f_H2(z=6.6)\sim0.5 f_H2(z=5.7) when averaged over all LAEs; they arise due to the interplay between the H2 formation/shielding controlled by dust and the intensity of the ultraviolet (UV) Lyman-Werner photo-dissociating radiation produced by stars. We then predict the carbon monoxide (CO) luminosities for such LAEs and check that they are consistent with the upper limits found by Wagg et al. (2009) for two z>6 LAEs. At z\sim(5.7, 6.6), the lowest CO rotational transition observable for both samples with the actual capabilities of Atacama Large Millimeter Array (ALMA) is the CO(6-5). We find that at z\sim5.7, about 1-2% of LAEs, i.e., those with an observed Lyman Alpha luminosity larger than 10^{43.2} erg/s would be detectable with an integration time of 5-10 hours (S/N=5); at z\sim6.6 none of the LAEs would be detectable in CO, even with an ALMA integration time of 10 hours. We also build the CO `flux function', i.e., the number density of LAEs as a function of the line-integrated CO flux, S_CO, and show that it peaks at S_CO = 0.1 mJy at z\sim5.7, progressively shifting to lower values at higher redshifts. We end by discussing the model uncertainties.

preprint2011arXiv

Ancient giants: on the farthest galaxy at z=8.6

The observational frontiers for the detection of high-redshift galaxies have recently been pushed to unimaginable distances with the record-holding Lyman Alpha Emitter (LAE) UDFy-38135539 discovered at redshift z=8.6. However, the physical nature and the implications of this discovery have yet to be assessed. By selecting galaxies with observed luminosities similar to UDFy-38135539 in state-of-the-art cosmological simulations tuned to reproduce the large scale properties of LAEs, we bracket the physical nature of UDFy-38135539: it has a star formation rate ~ 2.7-3.7 solar masses/yr, it contains ~ 10^{8.3-8.7} solar mass of stars 50-80 Myr old, with stellar metallicity ~ 0.03-0.12 of the solar value. For any of the simulated galaxies to be visible as a LAE in the observed range, the intergalactic neutral hydrogen fraction at z=8.6 must be <= 0.2 and extra ionizing radiation from sources clustered around UDFy-38135539 is necessary. Finally, we predict that there is a 70% (15%) probability of detecting at least 1 such source from JWST (HST/WFC3) observations in a physical radius ~ 0.4 Mpc around UDFy-38135539.

preprint2011arXiv

Cosmic Lighthouses : Unveiling the nature of high-redshift galaxies

We are in the golden age for the search for high-redshift galaxies, made possible by a combination of new instruments and innovative search techniques. One of the major aims of such searches is to constrain the epoch of reionization (EoR), which marks the second major change in the ionization state of the Universe. Understanding the EoR is difficult since whilst it is galaxy evolution which drives reionization, reionization itself influences galaxy evolution through feedback effects. Unraveling the interplay of reionization and galaxy evolution is further complicated by of a lack of understanding of the metal enrichment and dust distribution in high redshift galaxies. To this end, a class of galaxies called Lyman Alpha Emitters (LAEs) have been gaining enormous popularity as probes of all these three processes. In this thesis, we couple state of the art cosmological SPH simulations (GADGET-2) with a physically motivated, self-consistent model for LAEs, so as to be able to understand the importance of the intergalactic medium (IGM) ionization state, dust and peculiar velocities in shaping their observed properties. By doing so, the aim is to gain insight on the nature of LAEs, put precious constraints on their elusive physical properties and make predictions for future instruments such as the Atacama Large Millimeter Array (ALMA). Using our LAE model in conjunction with a code that builds the MW merger tree (GAMETE), we build a bridge between the high-redshift and the local Universe. We also use SPH simulations (GADGET-2) to study the nature of the earliest galaxies that have been detected as of yet, place constraints on their contribution to reionization, and predict their detectability using the next generation of instruments, such as the James Web Space Telescope (JWST).

preprint2011arXiv

Local Group progenitors: Lyman Alpha bright?

We present a novel approach of identifying the Milky Way (MW) and Andromeda (M31) progenitors that could be visible as LAEs at z~6: we couple a snapshot from the Constrained Local UniversE Simulations (CLUES) project, that successfully reproduces the MW and M31 galaxies situated in their correct environment, to a Lyman Alpha Emitter (LAE) model. Exploring intergalactic medium (IGM) ionization states ranging from an almost neutral to a fully ionized one, we find that including (excluding) the effects of clustered sources the first local group progenitor appears as a LAE for a neutral hydrogen fraction of 0.4 (0.1). This number increases to 5 progenitors each of the MW and M31 being visible as LAEs for a completely ionized IGM; the contribution from clustered sources is crucial in making many of the progenitors visible in the Lyman Alpha, for all the ionization states considered. We also put constrains on the physical properties of such progenitors. Finally, we find that the number density of these LAEs are higher than that of general field LAEs (observed in cosmological volumes) by about two (one) orders of magnitude for a neutral hydrogen fraction of 10^-5 (0.4). Detections of such high LAE number densities at z ~ 6 would be a clear signature of an over-dense region that could evolve and resemble the local group volume at z=0.

preprint2010arXiv

The visibility of Lyman Alpha Emitters during reionization

We present the first Lyman Alpha Emitter (LAE) study that combines: (i) cosmological SPH simulations run using GADGET-2, (ii) radiative transfer simulations (CRASH), and (iii) a previously developed LAE model. This complete LAE model accounts for the intrinsic LAE Lyman Alpha/continuum luminosity, dust enrichment and Lyman Alpha transmission through the intergalactic medium (IGM), to quantify the effects of reionization, dust and velocity fields on the Lyman Alpha and UV Luminosity Functions (LF). We find that a model neglecting dust sorely fails to reproduce either the slope or the magnitude of the observed Lyman Alpha and UV LFs. Clumped dust is required to simultaneously fit the observed UV and Lyman Alpha LFs, such that the intrinsic Lyman Alpha-to-continuum luminosity is enhanced by a factor f_alpha/f_c ~ 1.5 (3.7) excluding (including) peculiar velocities. The higher value including velocity fields arises since LAEs reside in large potential wells and inflows decrease their Lyman Alpha transmission. For the first time, a degeneracy is found between the the ionization state of the IGM and the clumping of dust inside high-redshift galaxies. The Lyman Alpha LF at z ~ 5.7 can be well reproduced (to within a 5 sigma error) by a wide range of IGM average neutral hydrogen fraction, 3.4/10^3 to 0.16, provided that the increase in the Lyman Alpha transmission through a more ionized IGM is compensated by a decrease in the Lyman Alpha escape fraction from the galaxy due to dust absorption. The physical properties of LAEs are presented, along with a discussion of the assumptions adopted.

preprint2009arXiv

Detecting Lyman Alpha Emitters in the Sub-millimeter

Using the results from a previously developed Lyman Alpha/continuum production/transmission and dust enrichment model for Lyman Alpha Emitters (LAEs), based on cosmological SPH simulations, we assess the detectability of their dust-reprocessed sub-millimeter (submm) radiation. As supernovae (rather than evolved stars) control dust formation and destruction processes, LAEs are relatively dust-poor with respect to local galaxies: they have low dust-to-gas ratios (0.05 times the dust-to-gas ratio of the Milky Way) in spite of their relatively high metallicity, which is approximately 10-50% of the solar value. Using the derived escape fraction of ultraviolet (UV) continuum photons we compute the UV luminosity absorbed by dust and re-emitted in the far infrared. The LAE submm fluxes correlate with their Lyman Alpha luminosity: about (3%, 1%) at z=(5.7, 6.6) of the LAEs in our simulated sample (those with with Log L_α> 43.1) would have fluxes at 850 micrometer (the optimal band for detection) in excess of 0.12 mJy and will be therefore detectable at 5 sigma with ALMA with an integration time of only 1 hour. Such detections would open a new window on the physical conditions prevailing in these most distant galaxies.

preprint2009arXiv

The cool side of Lyman Alpha Emitters

We extend a previous study of Lyman Alpha Emitters (LAEs) based on hydrodynamical cosmological simulations, by including two physical processes important for LAEs: (a) Lyman Alpha and continuum luminosities produced by cooling of collisionally excited HI in the galaxy, (b) dust formation and evolution; we follow these processes on a galaxy-by-galaxy basis. HI cooling on average contributes 16-18% of the Lyman Alpha radiation produced by stars, but this value can be much higher in low mass LAEs and further increased if the HI is clumpy. The continuum luminosity is instead almost completely dominated by stellar sources. The dust content of galaxies scales with their stellar mass, M_{dust} is proportional to M_*^0.7 and stellar metallicity, Z_*, such that M_{dust} is proportional to Z_*^1.7. As a result, the massive galaxies have Lyman Alpha escape fraction as low as f_alpha=0.1, with a LAE-averaged value decreasing with redshift such that it is (0.33,0.23) at z =(5.7,6.6). The UV continuum escape fraction shows the opposite trend with redshift, possibly resulting from clumpiness evolution. The model successfully reproduces the observed Lyman Alpha and UV luminosity functions at different redshifts and the Lyman Alpha equivalent width scatter to a large degree, although the observed distribution appears to be more more extended than the predicted one. We discuss possible reasons for such tension.