sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Dayal, Pratika

Normalized to: Dayal, P.

53 article(s) in total. 412 co-authors, from 1 to 24 common article(s). Median position in authors list is 2,0.

[1]  [pdf] - 2063373
A hint on the metal-free star formation rate density from 21cm-EDGES data
Comments: Submitted to MNRAS
Submitted: 2020-03-12
We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dot{\rho}_{*}^{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$\alpha$), 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$\alpha$ 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{\rho}_{*}^{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$\alpha$ 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.
[2]  [pdf] - 2049097
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
Comments: 18 pages, 7 figures, v.2: Accepted for publications in MNRAS
Submitted: 2019-05-01, last modified: 2020-02-16
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.
[3]  [pdf] - 2032988
Reionization with galaxies and active galactic nuclei
Comments: Submitted to MNRAS; comments welcome
Submitted: 2020-01-16
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 positively 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$.
[4]  [pdf] - 2025610
The 21cm bispectrum during reionization: a tracer of the ionization topology
Comments: 16 pages, 7 figures. Accepted for publication in MNRAS
Submitted: 2019-07-09, last modified: 2019-11-18
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.
[5]  [pdf] - 1894233
Ruling out 3 keV warm dark matter using 21 cm-EDGES data
Comments: Accepted to MNRAS
Submitted: 2019-02-25, last modified: 2019-05-22
Weakly interacting cold dark matter (CDM) particles, which are otherwise extremely successful in explaining various cosmological observations, exhibit a number of problems on small scales. One possible way of solving these problems is to invoke (so-called) warm dark matter (WDM) particles with masses $m_x \sim$ keV. Since the formation of structure is delayed in such WDM models, it is natural to expect that they can be constrained using observations related to the first stars, e.g., the 21 cm signal from cosmic dawn. In this work, we use a detailed galaxy formation model, Delphi, to calculate the 21 cm signal at high-redshifts and compare this to the recent EDGES observations. We find that while CDM and 5 keV WDM models can obtain a 21 cm signal within the observed redshift range, reproducing the amplitude of the observations requires the introduction of an excess radio background. On the other hand, WDM models with $m_x \leq 3$ keV can be ruled out since they are unable to match either the redshift range or the amplitude of the EDGES signal, irrespective of the parameters used. Comparable to values obtained from the low-redshift Lyman Alpha forest, our results extend constraints on the WDM particle to an era inaccessible by any other means; additional forthcoming 21 cm data from the era of cosmic dawn will be crucial in refining such constraints.
[6]  [pdf] - 1863871
The hierarchical assembly of galaxies and black holes in the first billion years: predictions for the era of gravitational wave astronomy
Comments: Accepted to MNRAS
Submitted: 2018-10-25, last modified: 2019-03-26
In this work we include black hole (BH) seeding, growth and feedback into our semi-analytic galaxy formation model, Delphi. Our model now fully tracks the, accretion- and merger-driven, hierarchical assembly of the dark matter halo, baryonic and BH masses of high-redshift ($z>5$) galaxies. We use a minimal set of mass- and $z$-independent free parameters associated with star formation and BH growth (and feedback) and include suppressed BH growth in low-mass galaxies to explore a number of physical scenarios including: (i) two types of BH seeds (stellar and those from Direct Collapse BH; DCBH); (ii) the impact of reionization feedback; and (iii) the impact of instantaneous versus delayed galaxy mergers on the baryonic growth. While both reionization feedback and delayed galaxy mergers have no sensible impact on the evolving ultra-violet luminosity function, the latter limits the maximum BH masses achieved at these high-$z$. We then use this model, baselined against all available high-$z$ galaxy and BH data-sets, to predict the LISA detectability of merger events at $z > 5$. As expected, the merger rate is dominated by stellar BH mergers for all scenarios and our model predicts an expected upper limit of about 20 mergers in the case of instantaneous merging and no reionization feedback over the 4-year mission duration. Including the impact of delayed mergers and reionization feedback reduces this to about 12 events over the same observational time-scale.
[7]  [pdf] - 1846964
The MUSE Ultra Deep Field (MUDF). I. Discovery of a group of Ly$\alpha$ nebulae associated with a bright $z\approx 3.23$ quasar pair
Comments: MNRAS Letters in press, 3 figures, revised version: minor changes following proofs
Submitted: 2019-03-01, last modified: 2019-03-12
We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations.
[8]  [pdf] - 1845984
Cosmic Dawn and Reionization: Astrophysics in the Final Frontier
Comments: Science White paper submitted to Astro2020 Decadal Survey
Submitted: 2019-03-08
The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.
[9]  [pdf] - 1845983
Astro2020 Science White Paper: A proposal to exploit galaxy-21cm synergies to shed light on the Epoch of Reionization
Comments: 6 pages, 3 figures
Submitted: 2019-03-08
This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.
[10]  [pdf] - 1832668
The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s
Akeson, Rachel; Armus, Lee; Bachelet, Etienne; Bailey, Vanessa; Bartusek, Lisa; Bellini, Andrea; Benford, Dominic; Bennett, David; Bhattacharya, Aparna; Bohlin, Ralph; Boyer, Martha; Bozza, Valerio; Bryden, Geoffrey; Novati, Sebastiano Calchi; Carpenter, Kenneth; Casertano, Stefano; Choi, Ami; Content, David; Dayal, Pratika; Dressler, Alan; Doré, Olivier; Fall, S. Michael; Fan, Xiaohui; Fang, Xiao; Filippenko, Alexei; Finkelstein, Steven; Foley, Ryan; Furlanetto, Steven; Kalirai, Jason; Gaudi, B. Scott; Gilbert, Karoline; Girard, Julien; Grady, Kevin; Greene, Jenny; Guhathakurta, Puragra; Heinrich, Chen; Hemmati, Shoubaneh; Hendel, David; Henderson, Calen; Henning, Thomas; Hirata, Christopher; Ho, Shirley; Huff, Eric; Hutter, Anne; Jansen, Rolf; Jha, Saurabh; Johnson, Samson; Jones, David; Kasdin, Jeremy; Kelly, Patrick; Kirshner, Robert; Koekemoer, Anton; Kruk, Jeffrey; Lewis, Nikole; Macintosh, Bruce; Madau, Piero; Malhotra, Sangeeta; Mandel, Kaisey; Massara, Elena; Masters, Daniel; McEnery, Julie; McQuinn, Kristen; Melchior, Peter; Melton, Mark; Mennesson, Bertrand; Peeples, Molly; Penny, Matthew; Perlmutter, Saul; Pisani, Alice; Plazas, Andrés; Poleski, Radek; Postman, Marc; Ranc, Clément; Rauscher, Bernard; Rest, Armin; Roberge, Aki; Robertson, Brant; Rodney, Steven; Rhoads, James; Rhodes, Jason; Ryan, Russell; Sahu, Kailash; Sand, David; Scolnic, Dan; Seth, Anil; Shvartzvald, Yossi; Siellez, Karelle; Smith, Arfon; Spergel, David; Stassun, Keivan; Street, Rachel; Strolger, Louis-Gregory; Szalay, Alexander; Trauger, John; Troxel, M. A.; Turnbull, Margaret; van der Marel, Roeland; von der Linden, Anja; Wang, Yun; Weinberg, David; Williams, Benjamin; Windhorst, Rogier; Wollack, Edward; Wu, Hao-Yi; Yee, Jennifer; Zimmerman, Neil
Comments: 14 pages, 3 figures, 2 tables
Submitted: 2019-02-14
The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.
[11]  [pdf] - 1791061
Early galaxy formation and its large-scale effects
Comments: Solicited review published in Physics Reports. Typeset PDF available (till mid-January) at
Submitted: 2018-09-24, last modified: 2018-11-23
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.
[12]  [pdf] - 1769130
Science with an ngVLA: [CII] 158$\mu$m Emission from $z \ge 10$ Galaxies
Comments: 8 pages, 3 figures; to be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: substantial text overlap with arXiv:1709.02001
Submitted: 2018-10-16
We consider the capabilities of ALMA and the ngVLA to detect and image the[CII] 158\,$\mu$m line from galaxies into the cosmic `dark ages' ($z \sim 10$ to 20). The [CII] line may prove to be a powerful tool in determining spectroscopic redshifts, and galaxy dynamics, for the first galaxies. In 40\,hr, ALMA has the sensitivity to detect the integrated [CII] line emission from a moderate metallicity, active star-forming galaxy [$Z_A = 0.2\,Z_{\odot}$; star formation rate (SFR)= 5\,$M_\odot$\,yr$^{-1}$], at $z = 10$ at a significance of 6$\sigma$. The ngVLA will detect the integrated [CII] line emission from a Milky-Way like star formation rate galaxy ($Z_{A} = 0.2\,Z_{\odot}$, SFR = 1\,$M_\odot$\,yr$^{-1}$), at $z = 15$ at a significance of 6$\sigma$. Imaging simulations show that the ngVLA can determine rotation dynamics for active star-forming galaxies at $z \sim 15$, if they exist. The [CII] detection rate in blind surveys will be slow (of order unity per 40\,hr pointing.
[13]  [pdf] - 1767642
Probing the fluctuating Ultra-violet background using the Hubble Frontier Fields
Comments: Accepted for publication in MNRAS Letters
Submitted: 2018-09-05, last modified: 2018-09-30
In recent years, the rise in the number of Lyman Break Galaxies detected at high redshifts z >= 6 has opened up the possibility of understanding early galaxy formation physics in great detail. In particular, the faint-end slope (alpha) of the Ultra-violet luminosity function (UV LF) of these galaxies is a potential probe of feedback effects that suppress star formation in low mass haloes. In this work, we propose a proof-of-concept calculation for constraining the fluctuating UV background during reionization by constraining alpha in different volumes of the Universe. Because of patchy reionization, different volumes will experience different amount of photo-heating which should lead to a scatter in the measured alpha. Our approach is based on a simple model of the UV LF that is a scaled version of the halo mass function combined with an exponential suppression in the galaxy luminosity at the faint-end because of UV feedback. Although current data is not sufficient to constrain alpha in different fields, we expect that, in the near future, observations of the six lensed Hubble Frontier Fields with the James Webb Space Telescope (JWST) will offer an ideal test of our concept.
[14]  [pdf] - 1732812
Spectroscopic investigation of a reionized galaxy overdensity at z=7
Comments: Accepted for publication in The Astrophysical Journal Letters. 6 pages, 3 figures
Submitted: 2018-07-24
We present deep spectroscopic follow-up observations of the Bremer Deep Field (BDF) where the two $z\sim$7 bright Ly$\alpha$ emitters (LAE) BDF521 and BDF3299 were previously discovered by Vanzella et al. (2011) and where a factor of $\sim$3-4 overdensity of faint LBGs has been found by Castellano et al. (2016). We confirm a new bright Ly$\alpha$ emitter, BDF2195, at the same redshift of BDF521, $z=7.008$, and at only $\sim$90 kpc physical distance from it, confirming that the BDF area is likely an overdense, reionized region. A quantitative assessment of the Ly$\alpha$ fraction shows that the number of detected bright emitters is much higher than the average found at z$\sim$7, suggesting a high Ly$\alpha$ transmission through the inter-galactic medium (IGM). However, the line visibility from fainter galaxies is at odds with this finding, since no Ly$\alpha$ emission is found in any of the observed candidates with $M_{UV}>$-20.25. This discrepancy can be understood either if some mechanism prevents Ly$\alpha$ emission from fainter galaxies within the ionized bubbles from reaching the observer, or if faint galaxies are located outside the reionized area and bright LAEs are solely responsible for the creation of their own HII regions. A thorough assessment of the nature of the BDF region and of its sources of re-ionizing radiation will be made possible by JWST spectroscopic capabilities.
[15]  [pdf] - 1709552
Survey parameters for detecting 21cm - Ly$\alpha$ emitter cross correlations with the Square Kilometre Array
Comments: 5 pages, 3 figures, accepted for publication in MNRAS
Submitted: 2018-06-20
Detections of the cross correlation signal between the 21cm signal during reionization and high-redshift Lyman Alpha emitters (LAEs) are subject to observational uncertainties which mainly include systematics associated with radio interferometers and LAE selection. These uncertainties can be reduced by increasing the survey volume and/or the survey luminosity limit, i.e. the faintest detectable Lyman Alpha (Ly$\alpha$) luminosity. We use our model of high-redshift LAEs and the underlying reionization state to compute the uncertainties of the 21cm-LAE cross correlation function at $z\simeq6.6$ for observations with SKA1-Low and LAE surveys with $\Delta z=0.1$ for three different values of the average IGM ionization state ($\langle\chi_\mathrm{HI}\rangle$=0.1, 0.25, 0.5). At $z\simeq6.6$, we find SILVERRUSH type surveys, with a field of view of 21 deg$^2$ and survey luminosity limits of $L_\alpha\geq7.9\times10^{42}$erg~s$^{-1}$, to be optimal to distinguish between an inter-galactic medium (IGM) that is 50%, 25% and 10% neutral, while surveys with smaller fields of view and lower survey luminosity limits, such as the 5 and 10 deg$^2$ surveys with WFIRST, can only discriminate between a 50% and 10% neutral IGM.
[16]  [pdf] - 1659764
Probing the nature of Dark Matter through the metal enrichment of the intergalactic medium
Comments: Accepted to MNRAS
Submitted: 2018-03-14, last modified: 2018-03-20
We focus on exploring the metal enrichment of the intergalactic medium (IGM) in Cold and Warm (1.5 and 3 keV) Dark Matter (DM) cosmologies, and the constraints this yields on the DM particle mass, using a semi-analytic model, Delphi, that jointly tracks the Dark Matter and baryonic assembly of galaxies at $z \simeq 4-20$ including both Supernova and (a range of) reionization feedback (models). We find that while ${\rm M_{UV}} \geq -15$ galaxies contribute half of all IGM metals in the Cold Dark Matter model by $z \simeq 4.5$, given the suppression of low-mass halos, larger halos with ${\rm M_{UV}} \leq -15$ provide about 80\% of the IGM metal budget in 1.5 keV Warm Dark Matter models using two different models for the metallicity of the interstellar medium. Our results also show that the only models compatible with two different high-redshift data sets, provided by the evolving Ultra-Violet luminosity function at $z \simeq 6-10$ and IGM metal density (e.g. Simcoe et al. 2011), are standard Cold Dark Matter and 3 keV Warm DM that do not include any reionization feedback; a combination of the UV LF and the Diaz et al. (2016) points provides a weaker constraint, allowing Cold and 3 keV and 1.5 keV Warm DM models with SN feedback only, as well as CDM with complete gas suppression of all halos with $v_{circ} \leq 30\, {\rm km\, s^{-1}}$. Tightening the error bars on the IGM metal enrichment, future observations, at $z \geq 5.5$, could therefore represent an alternative way of shedding light on the nature of Dark Matter.
[17]  [pdf] - 1588041
Galaxies into the Dark Ages
Comments: 11pages, 8 figures, Accepted for the Astrophysical Journal
Submitted: 2017-09-06
We consider the capabilities of current and future large facilities operating at 2\,mm to 3\,mm wavelength to detect and image the [CII] 158\,$\mu$m line from galaxies into the cosmic "dark ages" ($z \sim 10$ to 20). The [CII] line may prove to be a powerful tool in determining spectroscopic redshifts, and galaxy dynamics, for the first galaxies. We emphasize that the nature, and even existence, of such extreme redshift galaxies, remains at the frontier of open questions in galaxy formation. In 40\,hr, ALMA has the sensitivity to detect the integrated [CII] line emission from a moderate metallicity, active star-forming galaxy [$Z_A = 0.2\,Z_{\odot}$; star formation rate (SFR) = 5\,$M_\odot$\,yr$^{-1}$], at $z = 10$ at a significance of 6$\sigma$. The next-generation Very Large Array (ngVLA) will detect the integrated [CII] line emission from a Milky-Way like star formation rate galaxy ($Z_{A} = 0.2\,Z_{\odot}$, SFR = 1\,$M_\odot$\,yr$^{-1}$), at $z = 15$ at a significance of 6$\sigma$. Imaging simulations show that the ngVLA can determine rotation dynamics for active star-forming galaxies at $z \sim 15$, if they exist. Based on our very limited knowledge of the extreme redshift Universe, we calculate the count rate in blind, volumetric surveys for [CII] emission at $z \sim 10$ to 20. The detection rates in blind surveys will be slow (of order unity per 40\,hr pointing). However, the observations are well suited to commensal searches. We compare [CII] with the [OIII] 88$\mu$m line, and other ancillary information in high $z$ galaxies that would aid these studies.
[18]  [pdf] - 1582998
Warm dark matter constraints from high-$z$ Direct Collapse Black Holes using the JWST
Comments: Accepted to MNRAS with minor revision
Submitted: 2017-05-01, last modified: 2017-09-02
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.
[19]  [pdf] - 1585537
Theia: Faint objects in motion or the new astrometry frontier
The Theia Collaboration; Boehm, Celine; Krone-Martins, Alberto; Amorim, Antonio; Anglada-Escude, Guillem; Brandeker, Alexis; Courbin, Frederic; Ensslin, Torsten; Falcao, Antonio; Freese, Katherine; Holl, Berry; Labadie, Lucas; Leger, Alain; Malbet, Fabien; Mamon, Gary; McArthur, Barbara; Mora, Alcione; Shao, Michael; Sozzetti, Alessandro; Spolyar, Douglas; Villaver, Eva; Albertus, Conrado; Bertone, Stefano; Bouy, Herve; Boylan-Kolchin, Michael; Brown, Anthony; Brown, Warren; Cardoso, Vitor; Chemin, Laurent; Claudi, Riccardo; Correia, Alexandre C. M.; Crosta, Mariateresa; Crouzier, Antoine; Cyr-Racine, Francis-Yan; Damasso, Mario; da Silva, Antonio; Davies, Melvyn; Das, Payel; Dayal, Pratika; de Val-Borro, Miguel; Diaferio, Antonaldo; Erickcek, Adrienne; Fairbairn, Malcolm; Fortin, Morgane; Fridlund, Malcolm; Garcia, Paulo; Gnedin, Oleg; Goobar, Ariel; Gordo, Paulo; Goullioud, Renaud; Hambly, Nigel; Hara, Nathan; Hobbs, David; Hog, Erik; Holland, Andrew; Ibata, Rodrigo; Jordi, Carme; Klioner, Sergei; Kopeikin, Sergei; Lacroix, Thomas; Laskar, Jacques; Poncin-Lafitte, Christophe Le; Luri, Xavier; Majumdar, Subhabrata; Makarov, Valeri; Massey, Richard; Mennesson, Bertrand; Michalik, Daniel; de Almeida, Andre Moitinho; Mourao, Ana; Moustakas, Leonidas; Murray, Neil; Muterspaugh, Matthew; Oertel, Micaela; Ostorero, Luisa; Perez-Garcia, Angeles; Platais, Imants; de Mora, Jordi Portell i; Quirrenbach, Andreas; Randall, Lisa; Read, Justin; Regos, Eniko; Rory, Barnes; Rybicki, Krzysztof; Scott, Pat; Schneider, Jean; Scholtz, Jakub; Siebert, Arnaud; Tereno, Ismael; Tomsick, John; Traub, Wesley; Valluri, Monica; Walker, Matt; Walton, Nicholas; Watkins, Laura; White, Glenn; Evans, Dafydd Wyn; Wyrzykowski, Lukasz; Wyse, Rosemary
Submitted: 2017-07-02
In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.
[20]  [pdf] - 1530222
Reionization and Galaxy Formation in Warm Dark Matter Cosmologies
Comments: Accepted to ApJ
Submitted: 2015-01-12, last modified: 2017-02-19
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.
[21]  [pdf] - 1530670
Exploring 21cm - Lyman Alpha emitter synergies for SKA
Comments: 20 pages, 9 figures, accepted in ApJ, comments welcome
Submitted: 2016-05-05, last modified: 2017-01-10
We study the signatures of reionization and ionizing properties of the early galaxies in the cross-correlations between the 21cm emission from the spin-flip transition of neutral hydrogen (HI) and the underlying galaxy population, in particular a sub-population of galaxies visible as Lyman Alpha Emitters (LAEs). With both observables simultaneously derived from a $z\simeq6.6$ hydrodynamical simulation (GADGET-2) snapshot post-processed with a radiative transfer code (pCRASH) and a dust model, we perform a parameter study and aim to constrain both the average intergalactic medium (IGM) ionization state ($1-\langle \chi_{HI} \rangle$) and the reionization topology (outside-in versus inside-out). We find that in our model LAEs occupy the densest and most-ionized regions resulting in a very strong anti-correlation between the LAEs and the 21cm emission. A 1000h SKA-LOW1 - Subaru Hyper Suprime Cam experiment can provide constraints on $\langle \chi_{HI} \rangle$, allowing us to distinguish between IGM ionization levels of 50%, 25%, 10% and fully ionized at scales $r<10$ comoving Mpc (assuming foreground avoidance for SKA). Our results support the inside-out reionization scenario where the densest knots (under-dense voids) are ionized first (last) for $\langle \chi_{HI} \rangle >= 0.1$. Further, 1000h SKA-LOW1 observations should be able to confirm the inside-out scenario by detecting a lower 21cm brightness temperature (by about 2-10 mK) in the densest regions ($> 2$ arcminute scales) hosting LAEs compared to lower-density regions devoid of them.
[22]  [pdf] - 1470835
Coevolution of metallicity and star formation in galaxies to z=3.7: I. A fundamental plane
Comments: 19 pages, 10 figures, 4 tables, accepted for publication in MNRAS
Submitted: 2016-08-18
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).
[23]  [pdf] - 1470836
Coevolution of metallicity and star formation in galaxies to z=3.7: II. A theoretical model
Comments: 13 pages, 6 figures, 1 table, accepted for publication in MNRAS
Submitted: 2016-08-18
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.
[24]  [pdf] - 1441106
Interpreting the evolution of galaxy colours from $z = 8$ to $z = 5$
Comments: accepted in MNRAS
Submitted: 2016-07-14, last modified: 2016-07-20
We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, $\beta$, 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$.
[25]  [pdf] - 1430765
The habitability of the Universe through 13 billion years of cosmic time
Comments: Under refereeing in ApJ
Submitted: 2016-06-29
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.
[26]  [pdf] - 1359220
First observational support for overlapping reionized bubbles generated by a galaxy overdensity
Comments: Accepted for publication in The Astrophysical Journal Letters. 6 pages, 5 figures
Submitted: 2016-01-13
We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.
[27]  [pdf] - 1570877
Evaluating Galactic Habitability Using High Resolution Cosmological Simulations of Galaxy Formation
Comments: 24 pages, 33 figures, accepted for publication in the International Journal of Astrobiology
Submitted: 2015-11-05
We present the first model that couples high-resolution simulations of the formation of Local Group galaxies with calculations of the galactic habitable zone (GHZ), a region of space which has sufficient metallicity to form terrestrial planets without being subject to hazardous radiation. These simulations allow us to make substantial progress in mapping out the asymmetric three-dimensional GHZ and its time evolution for the Milky Way (MW) and Triangulum (M33) galaxies, as opposed to works that generally assume an azimuthally symmetric GHZ. Applying typical habitability metrics to MW and M33, we find that while a large number of habitable planets exist as close as a few kiloparsecs from the galactic centre, the probability of individual planetary systems being habitable rises as one approaches the edge of the stellar disc. Tidal streams and satellite galaxies also appear to be fertile grounds for habitable planet formation. In short, we find that both galaxies arrive at similar GHZs by different evolutionary paths, as measured by the first and third quartiles of surviving biospheres. For the Milky Way, this interquartile range begins as a narrow band at large radii, expanding to encompass much of the galaxy at intermediate times before settling at a range of 2-13kpc. In the case of M33, the opposite behaviour occurs - the initial and final interquartile ranges are quite similar, showing gradual evolution. This suggests that galaxy assembly history strongly influences the time evolution of the GHZ, which will affect the relative time lag between biospheres in different galactic locations. We end by noting the caveats involved in such studies and demonstrate that high resolution cosmological simulations will play a vital role in understanding habitability on galactic scales, provided that these simulations accurately resolve chemical evolution.
[28]  [pdf] - 1273216
The quest for cradles of life: using the fundamental metallicity relation to hunt for the most habitable type of galaxy
Comments: Accepted to APJ letters. Comments welcome
Submitted: 2015-07-15
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.
[29]  [pdf] - 1450563
The Cosmic Dawn and Epoch of Reionization with the Square Kilometre Array
Comments: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015. PoS(AASKA14)001
Submitted: 2015-05-28
Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).
[30]  [pdf] - 1182789
Early galaxy formation in warm dark matter cosmologies
Comments: Accepted for publication in ApJ
Submitted: 2014-08-05, last modified: 2015-05-07
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.
[31]  [pdf] - 1810100
The dust mass in z > 6 normal star forming galaxies
Comments: 5 pages, 2 figures, accepted for publication as a Letter to MNRAS
Submitted: 2015-05-07
We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with Mstar > 10^9 Msun allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population
[32]  [pdf] - 989942
Clustering and lifetime of Lyman Alpha Emitters in the Epoch of Reionization
Comments: 11 pages, 7 figures, accepted for publication in MNRAS, comments welcome
Submitted: 2015-03-17, last modified: 2015-05-07
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.
[33]  [pdf] - 1209622
Essential physics of early galaxy formation
Comments: Accepted to MNRAS with minor changes
Submitted: 2014-05-19, last modified: 2014-09-05
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 $\alpha = -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$.
[34]  [pdf] - 1215233
Constraining dust formation in high-redshift young galaxies
Comments: 10 pages, 3 figures, accepted for publication in MNRAS
Submitted: 2014-06-25
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.
[35]  [pdf] - 818774
The visibility of Lyman Alpha Emitters: constraining reionization, ionizing photons and dust
Comments: 20 pages, 13 figures, 1 table, comments welcome, accepted for publication in MNRAS
Submitted: 2014-01-13, last modified: 2014-05-06
We build a physical model for high-redshift Lyman Alpha emitters (LAEs) by coupling state of the art cosmological simulations (GADGET-2) with a dust model and a radiative transfer code (pCRASH). We post-process the cosmological simulation with pCRASH using five different values of the escape fraction of hydrogen ionizing photons (f_esc=0.05,0.25,0.5,0.75,0.95) until reionization is complete, i.e. the average neutral hydrogen fraction drops to <X_HI>~10^-4. Then, the only free-parameter left to match model results to the observed Lya and UV luminosity functions of LAEs at z~6.6 is the relative escape of Lyman Alpha (Lya) and continuum photons from the galactic environment (f_alpha/f_c). We find a three-dimensional degeneracy such that the theoretical model can be reconciled with observations for an IGM Lya transmission <T_alpha>_LAE~38-50% (which translates to <X_HI>~0.5-10^-4 for Gaussian emission lines), f_esc~0.05-0.50 and f_alpha/f_c~0.6-1.8.
[36]  [pdf] - 802582
The colour distribution of galaxies at redshift five
Comments: 13 pages, 8 colour figures. Updated in response to referee report; accepted in MNRAS
Submitted: 2013-12-17, last modified: 2014-03-24
We present the results of a study investigating the rest-frame ultra-violet (UV) spectral slopes of redshift z~5 Lyman-break galaxies (LBGs). By combining deep Hubble Space Telescope imaging of the CANDELS and HUDF fields with ground-based imaging from the UKIDSS Ultra Deep Survey (UDS), we have produced a large sample of z~5 LBGs spanning an unprecedented factor of >100 in UV luminosity. Based on this sample we find a clear colour-magnitude relation (CMR) at z~5, such that the rest-frame UV slopes (beta) of brighter galaxies are notably redder than their fainter counterparts. We determine that the z~5 CMR is well described by a linear relationship of the form: d beta = (-0.12 +/- 0.02) d Muv, with no clear evidence for a change in CMR slope at faint magnitudes (i.e. Muv > -18.9). Using the results of detailed simulations we are able, for the first time, to infer the intrinsic (i.e. free from noise) variation of galaxy colours around the CMR at z~5. We find significant (12 sigma) evidence for intrinsic colour variation in the sample as a whole. Our results also demonstrate that the width of the intrinsic UV slope distribution of z~5 galaxies increases from Delta(beta)=0.1 at Muv=-18 to Delta(beta)=0.4 at Muv=-21. We suggest that the increasing width of the intrinsic galaxy colour distribution and the CMR itself are both plausibly explained by a luminosity independent lower limit of beta=-2.1, combined with an increase in the fraction of red galaxies in brighter UV-luminosity bins.
[37]  [pdf] - 1157533
Simulating the assembly of galaxies at redshifts z = 6 - 12
Comments: Accepted for publication in MNRAS
Submitted: 2012-11-05, last modified: 2013-06-18
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.
[38]  [pdf] - 1158270
The UV continua and inferred stellar populations of galaxies at z ~ 7 - 9 revealed by the Hubble Ultra Deep Field 2012 campaign
Comments: 15 pages, 7 figures, updated to match accepted version, now including an appendix tabulating the individual beta values for the objects used in this study
Submitted: 2012-12-04, last modified: 2013-04-22
We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope.
[39]  [pdf] - 1164907
CLUES to the past: Local Group progenitors amongst high-redshift Lyman Break Galaxies
Comments: Accepted for publication in MNRAS
Submitted: 2013-02-27, last modified: 2013-03-11
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.
[40]  [pdf] - 1116806
The physics of the fundamental metallicity relation
Comments: Accepted for publication in MNRAS
Submitted: 2012-02-21, last modified: 2013-01-18
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.
[41]  [pdf] - 1093050
Molecular hydrogen in Lyman Alpha Emitters
Comments: Accepted for publication in MNRAS
Submitted: 2012-01-17
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.
[42]  [pdf] - 1083772
Lyman Alpha Emitters and Lyman Break Galaxies: dichotomous twins?
Comments: Accepted for publication in MNRAS; figures updated with latest observational data sets
Submitted: 2011-09-01, last modified: 2012-01-04
(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.
[43]  [pdf] - 1078284
Local Group progenitors: Lyman Alpha bright?
Comments: Accepted for publication in MNRAS Letters
Submitted: 2011-07-28, last modified: 2011-10-05
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.
[44]  [pdf] - 1051985
Ancient giants: on the farthest galaxy at z=8.6
Comments: Accepted to MNRAS letters
Submitted: 2011-02-08, last modified: 2011-07-11
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.
[45]  [pdf] - 956181
Cosmic Lighthouses : Unveiling the nature of high-redshift galaxies
Comments: PhD Thesis
Submitted: 2011-01-20
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).
[46]  [pdf] - 1025833
Simulating high-redshift galaxies
Comments: 14 pages, 10 figures, MNRAS in press
Submitted: 2010-03-19, last modified: 2011-01-12
Recent observations have gathered a considerable sample of high redshift galaxy candidates and determined the evolution of their luminosity function (LF). To interpret these findings, we use cosmological SPH simulations including, in addition to standard physical processes, a detailed treatment of the Pop III-Pop II transition in early objects. The simulated high-z galaxies match remarkably well the amplitude and slope of the observed LF in the redshift range 5<z<10. The LF shifts towards fainter luminosities with increasing redshift, while its faint-end slope keeps an almost constant value, \alpha ~-2. The stellar populations of high-z galaxies have ages of 100-300 (40-130) Myr at z=5 (z=7-8), implying an early (z>9.4) start of their star formation activity; the specific star formation rate is almost independent of galactic stellar mass. These objects are enriched rapidly with metals and galaxies identified by HST/WFC3 (M_UV < -18) show metallicities ~0.1 Zsun even at z=7-8. Most of the simulated galaxies at z~7 (noticeably the smallest ones) are virtually dust-free, and none of them has an extinction larger than E(B-V) = 0.01. The bulk (50%) of the ionizing photons is produced by objects populating the faint-end of the LF (M_UV < -16), which JWST will resolve up to z=7.3. PopIII stars continue to form essentially at all redshifts; however, at z=6 (z=10) the contribution of Pop III stars to the total galactic luminosity is always less than 5% for M_UV < -17 (M_UV < -16). The typical high-z galaxies closely resemble the GRB host galaxy population observed at lower redshifts, strongly encouraging the use of GRBs to detect the first galaxies.
[47]  [pdf] - 1032708
High redshift Lya emitters: clues on the Milky Way infancy
Comments: Accepted for publication in MNRAS letter
Submitted: 2010-05-24
With the aim of determining if Milky Way (MW) progenitors could be identified as high redshift Lyman Alpha Emitters (LAEs) we have derived the intrinsic properties of z ~ 5.7 MW progenitors, which are then used to compute their observed Lyman-alpha luminosity, L_alpha, and equivalent width, EW. MW progenitors visible as LAEs are selected according to the canonical observational criterion, L_alpha > 10^42 erg/s and EW > 20 A. Progenitors of MW-like galaxies have L_alpha = 10^(39-43.25) erg/s, making some of them visible as LAEs. In any single MW merger tree realization, typically only 1 (out of ~ 50) progenitor meets the LAE selection criterion, but the probability to have at least one LAE is very high, P = 68%. The identified LAE stars have ages, t_* ~ 150-400 Myr at z ~ 5.7 with the exception of five small progenitors with t_* < 5 Myr and large EW = 60-130 A. LAE MW progenitors provide > 10% of the halo very metal-poor stars [Fe/H] < -2, thus establishing a potentially fruitful link between high-z galaxies and the Local Universe.
[48]  [pdf] - 1024986
The visibility of Lyman Alpha Emitters during reionization
Comments: Submitted to MNRAS, comments welcome
Submitted: 2010-02-04
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.
[49]  [pdf] - 31747
Constraining cosmic reionization with quasar, gamma ray burst, and Lya emitter observations
Comments: 4 pages, 4 figures, to appear in the Proceedings `L'Universo quattro secoli dopo Galileo', 53rd Congresso SAIt, Pisa, May 2009
Submitted: 2009-12-14
We investigate the cosmic reionization history by comparing semi-analytical models of the Lya forest with observations of high-z quasars and gamma ray bursts absorption spectra. In order to constrain the reionization epoch z_rei, we consider two physically motivated scenarios in which reionization ends either early (ERM, z_rei>= 7) or late (LRM, z_rei~6). We analyze the transmitted flux in a sample of 17 quasars spectra at 5.7<z<6.4 and in the spectrum of the gamma ray burst 050904 at z=6.3, studying the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that current observational data do not require any sudden change in the ionization state of the IGM at z~6, favouring indeed a highly ionized Universe at these epochs, as predicted by the ERM. Moreover, we test the predictions of this model through Lya emitters observations, finding that the ERM provide a good fit to the evolution of the luminosity function of Lya emitting galaxies in the redshift range z=5.7-6.5. The overall result points towards an extended reionization process which starts at z>=11 and completes at z_rei>=7, in agreement with the recent WMAP5 data.
[50]  [pdf] - 1003394
Detecting Lyman Alpha Emitters in the Sub-millimeter
Comments: 6 pages, 3 figures, Accepted for publication in MNRAS with minor additions
Submitted: 2009-08-11, last modified: 2009-12-07
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_\alpha > 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.
[51]  [pdf] - 1003147
The cool side of Lyman Alpha Emitters
Comments: 9 pages, accepted for publication in MNRAS
Submitted: 2009-07-29, last modified: 2009-11-06
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.
[52]  [pdf] - 1002776
Lyman Alpha Emitter Evolution in the Reionization Epoch
Comments: 13 pages, 9 figures, accepted to MNRAS
Submitted: 2009-07-02, last modified: 2009-11-02
Combining cosmological SPH simulations with a previously developed Lyman Alpha production/transmission model and the Early Reionization Model (ERM, reionization ends at redshift z~7), we obtain Lyman Alpha and UV Luminosity Functions (LFs) for Lyman Alpha Emitters (LAEs) for redshifts between 5.7 and 7.6. Matching model results to observations at z~5.7 requires escape fractions of Lyman Alpha, f_alpha=0.3, and UV (non-ionizing) continuum photons, f_c=0.22, corresponding to a color excess, E(B-V)=0.15. We find that (i) f_c increases towards higher redshifts, due the decreasing mean dust content of galaxies, (ii) the evolution of f_alpha/f_c hints at the dust content of the ISM becoming progressively inhomogeneous/clumped with decreasing redshift. The clustering photoionization boost is important during the initial reionization phases but has little effect on the Lyman Alpha LF for a highly ionized IGM. Halo (stellar) masses are in the range 10.0 < \log M_h < 11.8 (8.1 < \log M_* < 10.4) with M_h \propto M_*^{0.64}. The star formation rates are between 3-120 solar masses per year, mass-weighted mean ages are greater than 20 Myr at all redshifts, while the mean stellar metallicity increases from Z=0.12 to 0.22 solar metallicity from z~7.6 to z~5.7; both age and metallicity positively correlate with stellar mass. The brightest LAEs are all characterized by large star formation rates and intermediate ages (~200 Myr), while objects in the faint end of the Lyman Alpha LF show large age and star formation rate spreads. With no more free parameters, the Spectral Energy Distributions of three LAE at z~5.7 observed by Lai et al. (2007) are well reproduced by an intermediate age (182-220 Myr) stellar population and the above E(B-V) value.
[53]  [pdf] - 14654
Signatures of reionization on Lyman alpha emitters
Comments: 15 pages, 10 figures, Accepted to MNRAS
Submitted: 2008-07-18, last modified: 2008-07-23
We use a semi-analytic model of Lyman alpha emitters (LAEs) to constrain the reionization history. By considering two physically motivated scenarios in which reionization ends either early (ERM, z_i ~ 7) or late (LRM, z_i ~ 6), we fix the global value of the IGM neutral fraction (e.g. chi_{HI}=3 times 10^{-4}, 0.15 at z=6.56 for the ERM and LRM, respectively) leaving only the star formation efficiency and the effective escape fraction of Lya photons as free parameters. The ERM fits the observed LAE luminosity function (LF) at z=5.7 and 6.56 requiring no redshift evolution or mass dependence of the star formation efficiency, and LAE star formation rates (SFR) of 3-103 solar masses/year, contributing approximately 8% of the cosmic SFR density at z=5.7. The LRM requires a physically uncomfortable drop of approximately 4.5 times in the SFR of the emitters from z=6.5 to 5.7. Thus, the data seem to imply that the Universe was already highly ionized at z=6.56. The mass-dependent Lya transmissivity is between 0.36-0.51 (ERM) and less than 0.26 (LRM) at z=6.56. The LF data at z=4.5 imply an extra Lya line damping factor of approximately 0.25 possibly due to dust; the presence of a (clumpy) dust component with E(B-V) ~ 0.28 is also required to reproduce the observed large Lya equivalent widths at the same redshift. Additional useful information can be extracted from the line profile (weighted) skewness, found to be S_W=10-17 Angstrom for the two reionization models, which shows an interesting L_alpha-chi_{HI} anti-correlation, holding under the model assumptions. The shortcomings of the model and strategies to overcome them are discussed.