Source author record

A. Ferrara

A. Ferrara 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

50works
6topics
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

50 published item(s)

preprint2025arXiv

Lyman-$α$ radiation pressure regulates star formation efficiency

Order-unity star formation efficiencies (SFE) in early galaxies may explain the overabundance of bright galaxies observed by JWST at high redshift. Here we show that Lyman-$α$ (Ly$α$) radiation pressure limits the gas mass converted into stars, particularly in primordial environments. We develop a shell model including Ly$α$ feedback, and validate it with one-dimensional hydrodynamical simulations. To account for Ly$α$ resonant scattering, we adopt the most recent force multiplier fits, including the effect of Ly$α$ photon destruction by dust grains. We find that, independently of their gas surface density $Σ_g$, clouds are disrupted on a timescale shorter than a free-fall time, and even before supernova explosions if $Σ_g \gtrsim 10^3\,M_{\odot}\ \rm pc^{-2}$. At $\log(Z/Z_{\odot}) = -2$, relevant for high-redshift galaxies, the SFE is $0.01 \lesssim \hatε_{*} \lesssim 0.66$ for $10^3 \lesssimΣ_g [M_{\odot}\ \rm pc^{-2}] \lesssim 10^5$. The SFE is even lower for decreasing metallicity. Near-unity SFEs are possible only for extreme surface densities, $Σ_{g} \gtrsim 10^5\;M_{\odot}\ \rm pc^{-2}$, and near-solar metallicities. We conclude that Ly$α$ radiation pressure severely limits a possible extremely efficient, feedback-free phase of star formation in dense, metal-poor clouds.

preprint2025arXiv

Pushing JWST to the extremes: search and scrutiny of bright galaxy candidates at z$\simeq$15-30

We designed customized Lyman-break color selection techniques to identify galaxy candidates in the redshift ranges $15 \leq z \leq 20$ and $20 \leq z \leq 28$. The selection was performed on the ASTRODEEP-JWST multi-band catalogs of the CEERS, Abell-2744, JADES, NGDEEP, and PRIMER survey fields, covering a total area of $\sim0.2$ sq. deg. We identify five candidates at $15 \leq z \leq 20$, while no objects are found based on the $z\gtrsim20$ color selection criteria. Despite exhibiting a $>$1.5 mag break, all the objects display multimodal redshift probability distributions across different SED-fitting codes and methodologies. The alternative solutions correspond to poorly understood populations of low-mass quiescent or dusty galaxies at z$\sim$3-7. This conclusion is supported by the analysis of five F200W-dropout objects that we find to be interlopers on the basis of NIRSpec PRISM spectra: four dusty star-forming galaxies at z$\sim$2.2-6.6, and a passive galaxy at z=4.91 with log$(M_{\rm star}/{\rm M}_{\odot}) \lesssim$ 9. We measured the UV luminosity function under different assumptions on the contamination level within our sample. We find that if even a fraction of the candidates is indeed at $z\gtrsim15$, the resulting UV LF points to a very mild evolution compared to estimates at $z<15$, implying a significant tension with existing theoretical models. In particular, confirming our bright ($M_{\text{UV}}<-21$) candidates would require substantial revisions to the theoretical framework. In turn, if all these candidates will be confirmed to be interlopers, we conclude that future surveys may need ten times wider areas to select $M_{\text{UV}}\lesssim-20$ galaxies at $z>15$. Observations in the F150W and F200W filters at depths comparable to those in the NIRCam LW bands are also required to mitigate contamination from rare red objects at z$\lesssim$8.

preprint2024arXiv

CONCERTO: Extracting the power spectrum of the [C II ] emission line

CONCERTO is the first experiment to perform a [CII] line intensity mapping survey to target $z>5.2$. Measuring the [CII] power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds. We evaluate our ability to retrieve the [CII] signal in mock observations using the Simulated Infrared Dusty Extragalactic Sky. We compared two methods for dealing with the dust continuum emission from galaxies: the standard PCA and the arPLS method. For line interlopers, the strategy relies on masking low-redshift galaxies using external catalogues. As we do not have observations of CO or classical CO proxies ,we relied on the COSMOS stellar mass catalogue. To measure the power spectrum of masked data, we adapted the P of K EstimatoR and discuss its use on LIM data. The arPLS method achieves a reduction of the continuum background to a sub-dominant level of the [CII] at z=7 by a factor of>70. When using PCA, this factor is only 0.7. The masking lowers the power amplitude of line contamination down to $2 \times 10^2 Jy^2/sr$ This residual level is dominated by faint undetected sources. For our [CII] model, this results in a detection at z = 5.2 with a power ratio [CII]/(residual interlopers) = $62 \pm 32$ for a 22 % area survey loss. However, at z = 7, [C II ] / (residual interlopers)$=2.0 \pm 1.4$. Thanks to the large area covered by SIDES-Uchuu, we show that the power amplitude of line residuals varies by 12-15% for z=5.2-7. We present an end-to-end simulation of the extragalactic foreground removal that we ran to detect the [CII] at high redshift via its power spectrum. We show that dust continuum emission are not a limiting foreground for [CII] LIM. Residual CO and [CI] limits our ability to measure the [CII] power spectrum at z>7.

preprint2022arXiv

[CII] halos in ALPINE galaxies: smoking-gun of galactic outflows?

ALMA observations have revealed that many high redshift galaxies are surrounded by extended (10-15 kpc) [CII]-emitting halos which are not predicted by even the most advanced zoom-in simulations. Using a semi-analytical model, in a previous work we suggested that such halos are produced by starburst-driven, catastrophically cooling outflows. Here, we further improve the model and compare its predictions with data from 7 star-forming ($10\lesssim \rm SFR/ M_\odot \rm yr^{-1}<100$) galaxies at z=4-6, observed in the ALPINE survey. We find that (a) detected [CII] halos are a natural by-product of starburst-driven outflows; (b) the outflow mass loading factors are in the range $4\lesssimη\lesssim 7$, with higher $η$ values for lower-mass, lower-SFR systems, and scale with stellar mass as $η\propto M_*^{-0.43}$, consistently with the momentum-driven hypothesis. Our model suggests that outflows are widespread phenomena in high-z galaxies. However, in low-mass systems the halo extended [CII] emission is likely too faint to be detected with the current levels of sensitivity.

preprint2022arXiv

A survey of high-$z$ galaxies: SERRA simulations

We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 202 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $κ_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios can be reproduced by a part of the SERRA galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$).

preprint2022arXiv

Interpreting ALMA non-detections of JWST super-early galaxies

Recent attempts to detect [OIII] 88$μ$m emission from super-early ($z>10$) galaxy candidates observed by JWST have been unsuccessful. By using zoom-in simulations, we show that these galaxies are faint, and mostly fall below the local metal-poor $\rm [OIII]-SFR$ relation as a result of their low ionization parameter, $U_{\rm ion}\lesssim 10^{-3}$. Such low $U_{\rm ion}$ values are found in galaxies that are in an early assembly stage, and whose stars are still embedded in high-density natal clouds. However, the most luminous galaxy in our sample ($\rm{log}[L_{\rm{[OIII]}}/L_\odot] = 8.4$, $U_{\rm ion} \approx 0.1$) could be detected by ALMA in only $2.8$ hrs.

preprint2022arXiv

Reionization Era Bright Emission Line Survey: Selection and Characterization of Luminous Interstellar Medium Reservoirs in the z>6.5 Universe

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~79% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 28 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization.

preprint2022arXiv

The ALMA REBELS Survey: Average [CII] $158\,\rm{μm}$ sizes of Star-Forming Galaxies from $z\sim 7$ to $z\sim 4$

We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm μm}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.

preprint2022arXiv

The ALMA REBELS Survey: The Cosmic HI Gas Mass Density in Galaxies at $z\approx 7$

The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine HI 21-cm transition. Here we estimate the HI gas mass of a sample of main-sequence star-forming galaxies at $z\sim 6.5 - 7.8$ surveyed for [CII]$-158μ$m emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [CII]-to-HI conversion factor. We find that the HI gas mass excess in galaxies increases as a function of redshift, with an average of $M_{\rm HI} / M_\star \approx 10$, corresponding to HI gas mass fractions of $f_{\rm HI} = M_{\rm HI} / (M_\star + M_{\rm HI}) = 90\%$, at $z\approx 7$. Based on the [CII]-$158μ$m luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic HI gas mass density in galaxies ($ρ_{\rm HI}$) at this redshift, which we measure to be $ρ_{\rm HI} = 7.1^{+6.4}_{-3.0} \times 10^{6}\,M_{\odot}\,{\rm Mpc^{-3}}$. This estimate is substantially lower by a factor of $\approx 10$ than that inferred from an extrapolation of damped Lyman-$α$ absorber (DLA) measurements, and largely depend on the exact [CII] LF adopted. However, we find this decrease in $ρ_{\rm HI}$ to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sightlines predominantly probing the substantial fraction of HI gas in high-$z$ galactic halos, whereas [CII] traces the HI in the ISM associated with star formation. We make predictions for this build-up of neutral gas in galaxies as a function of redshift, showing that at $z\gtrsim 5$ only $\approx 10\%$ of the cosmic HI gas content is confined in galaxies and associated with the star-forming ISM.

preprint2022arXiv

The ALMA REBELS Survey: the dust content of $z \sim 7$ Lyman Break Galaxies

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $τ_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (ψ_{\rm UV}) = -0.05 ~[log (ψ)]^{2} + 0.86 ~log(ψ) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $τ_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions.

preprint2022arXiv

The ALMA REBELS Survey. Epoch of Reionization giants: properties of dusty galaxies at $z \approx 7$

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($β, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(β-β_{\rm int})$, where $β_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$β$ relation.

preprint2022arXiv

The ionizing properties of two bright Ly$α$ emitters in the BDF reionized bubble at z=7

We investigate the ionizing properties of the pair of bright Ly$α$ emitting galaxies BDF521 and BDF2195 at z=7.012 in order to constrain their contribution to the formation of the BDF "reionized bubble". We obtain constraints on UV emission lines (CIV$λ1548$ doublet, HeII$λ1640$, OIII]$λ1660$ doublet, and CIII]$λ1909$ doublet) from deep VLT-XSHOOTER observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and AGN emission. We use this spectroscopic information together with the photometry available in the field to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. We do not detect any significant emission at the expected position of UV lines, with 3$σ$ upper limits of EW$\lesssim$2-7 AA rest-frame. We find that the two objects have lower CIII] emission than expected on the basis of the correlation between the Ly$α$ and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at $\sim2-3σ$ significance, and only models with a $\lesssim$40% AGN contribution are compatible with the observations. The two objects are found to be relatively young ($\sim$20-30 Myrs) and metal-poor ($\lesssim 0.3 Z_{\odot}$) with stellar masses of a few $10^9M_{\odot}$. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log($ξ_{ion}^*$/Hz erg$^{-1}$)=25.02-25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in $z>$7 galaxies with strong CIII] and/or optical line emission ($\simeq$25.6-25.7). The two BDF emitters have no evidence of higher than average ionizing capabilities and are not capable of reionizing their surroundings by their own means under realistic assumptions on the escape fraction of ionizing photons. (abridged)

preprint2022arXiv

The REBELS ALMA Survey: cosmic dust temperature evolution out to z $\sim$ 7

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$μ$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures.

preprint2022arXiv

Widespread, strong outflows in XQR-30 quasars at the Reionisation epoch

Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dynamical mass correlation, thus following the black hole-dominance growth path. It is unknown what slowed down the black hole growth, leading towards the symbiotic growth observed in the local Universe, and when this process started, although black hole feedback is a likely driver. This deadlock is due to the lack of large, homogeneous samples of high-redshift quasars with high-quality, broad-band spectroscopic information. Here we report results from a VLT/X-shooter survey of 30 quasars at redshift 5.8$\le$z$\le$6.6 (XQR-30). About 50% of their spectra reveal broad blue-shifted absorption line (BAL) throughs, tracing powerful ionised winds. The BAL fraction in z$\gtrsim$6 quasars is 2-3 times higher than in quasars at z~2-4.5. XQR-30 BAL quasars exhibit extreme outflow velocities, up to 17% of the light speed, rarely observed at lower redshift. These outflows inject large amounts of energy into the galaxy interstellar medium, which can contrast nuclear gas accretion, slowing down the black-hole growth. The star-formation rate in high-z quasar hosts is generally $>$100 M$_\odot$/yr, so these galaxies are growing at a fast rate. The BAL phase may then mark the beginning of significant feedback, acting first on black hole growth and possibly later on galaxy growth. The red optical colors of BAL quasars at z$\gtrsim$6 indeed suggest that these systems are dusty and may be caught during an initial quenching phase of obscured accretion.

preprint2021arXiv

Infrared emission of $z \sim 6$ galaxies: AGN imprints

We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. In dusty (${M_{\rm d}\gtrsim 3\times 10^7 M_{\rm \odot}}$) galaxies, $\gtrsim 50-90 \%$ of the UV radiation is obscured by dust inhomogeneities on scales ${\gtrsim 100}$ pc. In runs with AGN, a clumpy, warm ($\approx 250$ K) dust component co-exists with a colder ($\approx 60$ K) and more diffuse one, heated by stars. Warm dust provides up to ${50 \%}$ of the total IR luminosity, but only $\lesssim 0.1 \%$ of the total mass content. The AGN boosts the MIR flux by ${10-100 \times}$ with respect to star forming galaxies, without significantly affecting the FIR. Our simulations successfully reproduce the observed SED of bright (${M_{\rm UV}\sim -26}$) ${z\sim 6}$ quasars, and show that these objects are part of complex, dust-rich merging systems, containing multiple sources (accreting BHs and/or star forming galaxies) in agreement with recent HST and ALMA observations. Our results show that the proposed ORIGINS missions will be able to investigate the MIR properties of dusty star forming galaxies and to obtain good quality spectra of bright quasars at $z\sim 6$. Finally, the MIR-to-FIR flux ratio of faint (${M_{\rm UV}\sim -24}$) AGN is ${>10\times}$ higher than for normal star forming galaxies. This implies that combined JWST/ORIGINS/ALMA observations will be crucial to identify faint and/or dust-obscured AGN in the distant Universe.

preprint2020arXiv

Shaping the structure of a GMC with radiation and winds

We study the effect of stellar feedback (photodissociation/ionization, radiation pressure and winds) on the evolution of a Giant Molecular Cloud (GMC), by means of a 3D radiative transfer, hydro-simulation implementing a complex chemical network featuring ${\rm H}_2$ formation and destruction. We track the formation of individual stars with mass $M>1\,{\rm M}_\odot$ with a stochastic recipe. Each star emits radiation according to its spectrum, sampled with 10 photon bins from near-infrared to extreme ultra-violet bands; winds are implemented by energy injection in the neighbouring cells. We run a simulation of a GMC with mass $M=10^5\,{\rm M}_\odot$, following the evolution of different gas phases. Thanks to the simultaneous inclusion of different stellar feedback mechanisms, we identify two stages in the cloud evolution: (1) radiation and winds carve ionized, low-density bubbles around massive stars, while FUV radiation dissociates most ${\rm H}_2$ in the cloud, apart from dense, self-shielded clumps; (2) rapid star formation (SFR$\simeq 0.1\,{\rm M}_\odot\,{\rm yr}^{-1}$) consumes molecular gas in the dense clumps, so that UV radiation escapes and ionizes the remaining HI gas in the GMC. ${\rm H}_2$ is exhausted in $1.6$ Myr, yielding a final star formation efficiency of 36 per cent. The average intensity of FUV and ionizing fields increases almost steadily with time; by the end of the simulation ($t=2.5$ Myr) we find $\langle G_0 \rangle \simeq 10^3$ (in Habing units), and a ionization parameter $\langle U_{\rm ion} \rangle \simeq 10^2$, respectively. The ionization field has also a more patchy distribution than the FUV one within the GMC. Throughout the evolution, the escape fraction of ionizing photons from the cloud is $f_{\rm ion, esc} < 0.03$.

preprint2020arXiv

Star formation law in the EoR from [CII] and CIII] lines

We present a novel method to simultaneously characterise the star formation law and the interstellar medium properties of galaxies in the Epoch of Reionization (EoR) through the combination of [CII] 158$μ$m (and its known relation with star formation rate) and CIII]$λ$1909Å emission line data. The method, based on a Markov Chain Monte Carlo algorithm, allows to determine the target galaxy average density, $n$, gas metallicity, $Z$, and burstiness parameter, $κ_s$, quantifying deviations from the Kennicutt-Schmidt relation. As an application, we consider COS-3018 (z=6.854), the only EoR Lyman Break Galaxy so far detected in both [CII] and CIII]. We show that COS-3018 is a moderate starburst ($κ_s \approx 3$), with $Z\approx 0.4\, Z_{\odot}$, and $n \approx 500\, {\rm cm^{-3}}$. Our method will be optimally applied to joint ALMA and JWST targets.

preprint2020arXiv

The stellar populations of high-redshift dwarf galaxies

We use high-resolution ($\approx 10$ pc), zoom-in simulations of a typical (stellar mass $M_\star\simeq10^{10}M_\odot$) Lyman Break Galaxy (LBG) at $z\simeq 6$ to investigate the stellar populations of its six dwarf galaxy satellites, whose stellar [gas] masses are in the range $\log (M_\star/M_\odot) \simeq 6-9$ [$\log (M_{gas}/M_\odot) \simeq4.3-7.75$]. The properties and evolution of satellites show no dependence on the distance from the central massive LBG ($< 11.5$ kpc). Instead, their star formation and chemical enrichment histories are tightly connected their stellar (and sub-halo) mass. High-mass dwarf galaxies ($\rm M_\star \gtrsim 5\times 10^8 M_\odot$) experience a long history of star formation, characterised by many merger events. Lower-mass systems go through a series of short star formation episodes, with no signs of mergers; their star formation activity starts relatively late ($z\approx 7$), and it is rapidly quenched by internal stellar feedback. In spite of the different evolutionary patterns, all satellites show a spherical morphology, with ancient and more metal-poor stars located towards the inner regions. All six dwarf satellites experienced high star formation rate ($\rm >5\,M_\odot yr ^{-1}$) bursts, which can be detected by JWST while targeting high-$z$ LBGs.

preprint2020arXiv

Warm dust in high-z galaxies: origin and implications

ALMA observations have revealed the presence of dust in galaxies in the Epoch of Reionization (redshift $z>6$). However, the dust temperature, $T_d$, remains unconstrained, and this introduces large uncertainties, particularly in the dust mass determinations. Using an analytical and physically-motivated model, we show that dust in high-$z$, star-forming giant molecular clouds (GMC), largely dominating the observed far-infrared luminosity, is warmer ($T_d > 60\ \mathrm{K}$) than locally. This is due to the more compact GMC structure induced by the higher gas pressure and turbulence characterizing early galaxies. The compactness also delays GMC dispersal by stellar feedback, thus $\sim 40\%$ of the total UV radiation emitted by newly born stars remains obscured. A higher $T_d$ has additional implications: it (a) reduces the tension between local and high-$z$ IRX-$β$ relation, (b) alleviates the problem of the uncomfortably large dust masses deduced from observations of some EoR galaxies.

preprint2016arXiv

Constraints on photoionization feedback from number counts of ultra-faint high-redshift galaxies in the Frontier Fields

We exploit a sample of ultra-faint high-redshift galaxies (demagnified HST $H_{160}$ magnitude $>30$) in the Frontier Fields clusters A2744 and M0416 to constrain a theoretical model for the UV luminosity function (LF) in the presence of photoionization feedback. The objects have been selected on the basis of accurate photometric redshifts computed from multi-band photometry including 7 HST bands and deep $K_s$ and IRAC observations. Magnification is computed on an object-by-object basis from all available lensing models of the two clusters. We take into account source detection completeness as a function of luminosity and size, magnification effects and systematics in the lens modeling of the clusters under investigation. We find that our sample of high-$z$ galaxies constrain the cut-off halo circular velocity below which star-formation is suppressed by photo-ionization feedback to $v_c^{\rm cut} < 50$ km s$^{-1}$. This circular velocity corresponds to a halo mass of $\approx5.6\times10^9~M_\odot$ and $\approx2.3\times10^9~M_\odot$ at $z=5$ and 10 respectively: higher mass halos can thus sustain continuous star formation activity without being quenched by external ionizing flux. More stringent constraints are prevented by the uncertainty in the modeling of the cluster lens, as embodied by systematic differences among the lens models available.

preprint2016arXiv

First observational support for overlapping reionized bubbles generated by a galaxy overdensity

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.

preprint2015arXiv

ALMA constraints on the faint millimetre source number counts and their contribution to the cosmic infrared background

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$μ$Jy, to derive differential number counts down to 60$μ$Jy and 100$μ$Jy at $λ=$1.3 mm and $λ=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$μ$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{λ= 1.1~mm} = 130~μ$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{λ= 1.3~mm} = 60~μ$Jy. At the faintest flux limits, i.e. 30$μ$Jy and 40$μ$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

preprint2015arXiv

Mapping metals at high redshift with far-infrared lines

Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $μ$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($ΔI/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($ν_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/μ{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4σ$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,μ$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.

preprint2015arXiv

Missing cosmic metals revealed by X-ray absorption towards distant sources

The census of heavy elements (metals) produced by all stars through cosmic times up to present-day is limited to ~50%; of these only half are still found within their parent galaxy. The majority of metals is expelled from galaxies into the circumgalactic (or even more distant, intergalactic) space by powerful galactic winds, leaving unpleasant uncertainty on the amount, thermal properties and distribution of these key chemical species. These dispersed metals unavoidably absorb soft X-ray photons from distant sources. We show that their integrated contribution can be detected in the form of increasing X-ray absorption with distance, for all kinds of high-energy cosmic sources. Based on extensive cosmological simulations, we assess that $\sim$ 10\% of all cosmic metals reside in the intergalactic medium. Most of the X-ray absorption arises instead from a few discrete structures along the line of sight. These extended structures, possibly pin-pointing galaxy groups, contain million degree, metal-enriched gas, 100-1,000 times denser than the cosmic mean. An additional ~10% of cosmic metals could reside in this phase.

preprint2015arXiv

On the [CII]-SFR relation in high redshift galaxies

After two ALMA observing cycles, only a handful of [CII] $158\,μm$ emission line searches in z>6 galaxies have reported a positive detection, questioning the applicability of the local [CII]-SFR relation to high-z systems. To investigate this issue we use the Vallini et al. 2013 (V13) model, based on high-resolution, radiative transfer cosmological simulations to predict the [CII] emission from the interstellar medium of a z~7 (halo mass $M_h=1.17\times10^{11}M_{\odot}$) galaxy. We improve the V13 model by including (a) a physically-motivated metallicity (Z) distribution of the gas, (b) the contribution of Photo-Dissociation Regions (PDRs), (c) the effects of Cosmic Microwave Background on the [CII] line luminosity. We study the relative contribution of diffuse neutral gas to the total [CII] emission ($F _{diff}/F_{tot}$) for different SFR and Z values. We find that the [CII] emission arises predominantly from PDRs: regardless of the galaxy properties, $F _{diff}/F_{tot}\leq 10$% since, at these early epochs, the CMB temperature approaches the spin temperature of the [CII] transition in the cold neutral medium ($T_{CMB}\sim T_s^{CNM}\sim 20$ K). Our model predicts a high-z [CII]-SFR relation consistent with observations of local dwarf galaxies ($0.02<Z/Z_{\odot}<0.5$). The [CII] deficit suggested by actual data ($L_{CII}<2.0\times 10^7 L_{\odot}$ in BDF3299 at z~7.1) if confirmed by deeper ALMA observations, can be ascribed to negative stellar feedback disrupting molecular clouds around star formation sites. The deviation from the local [CII]-SFR would then imply a modified Kennicutt-Schmidt relation in z>6 galaxies. Alternatively/in addition, the deficit might be explained by low gas metallicities ($Z<0.1 Z_{\odot}$).

preprint2015arXiv

Probing high-redshift galaxies with Ly$α$ intensity mapping

We present a study of the cosmological Ly$α$ emission signal at $z > 4$. Our goal is to predict the power spectrum of the spatial fluctuations that could be observed by an intensity mapping survey. The model uses the latest data from the HST legacy fields and the abundance matching technique to associate UV emission and dust properties with the halos, computing the emission from the interstellar medium (ISM) of galaxies and the intergalactic medium (IGM), including the effects of reionization, self-consistently. The Ly$α$ intensity from the diffuse IGM emission is 1.3 (2.0) times more intense than the ISM emission at $z = 4(7)$; both components are fair tracers of the star-forming galaxy distribution. However the power spectrum is dominated by ISM emission on small scales ($k > 0.01 h{\rm Mpc}^{-1}$) with shot noise being significant only above $k = 1 h{\rm Mpc}^{-1}$. At very lange scales ($k < 0.01h{\rm Mpc}^{-1}$) diffuse IGM emission becomes important. The comoving Ly$α$ luminosity density from IGM and galaxies, $\dot ρ_{{\rm Ly}α}^{\rm IGM} = 8.73(6.51) \times 10^{40} {\rm erg}{\rm s}^{-1}{\rm Mpc}^{-3}$ and $\dot ρ_{{\rm Ly}α}^{\rm ISM} = 6.62(3.21) \times 10^{40} {\rm erg}{\rm s}^{-1}{\rm Mpc}^{-3}$ at $z = 4(7)$, is consistent with recent SDSS determinations. We predict a power $k^3 P^{{\rm Ly}α}(k, z)/2π^2 = 9.76\times 10^{-4}(2.09\times 10^{-5}){\rm nW}^2{\rm m}^{-4}{\rm sr}^{-2}$ at $z = 4(7)$ for $k = 0.1 h {\rm Mpc}^{-1}$.

preprint2015arXiv

Probing the stellar initial mass function with high-$z$ supernovae

The first supernovae will soon be visible at the edge of the observable universe, revealing the birthplaces of Population III stars. With upcoming near-infrared missions, a broad analysis of the detectability of high-$z$ supernovae is paramount. We combine cosmological and radiation transport simulations, instrument specifications, and survey strategies to create synthetic observations of primeval core-collapse, Type IIn and pair-instability supernovae with the James Webb Space Telescope ($JWST$). We show that a dedicated observational campaign with the $JWST$ can detect up to $\sim 15$ pair-instability explosions, $\sim 300$ core-collapse supernovae, but less than one Type IIn explosion per year, depending on the Population III star formation history. Our synthetic survey also shows that $\approx 1-2 \times10^2$ supernovae detections, depending on the accuracy of the classification, are sufficient to discriminate between a Salpeter and flat mass distribution for high redshift stars with a confidence level greater than 99.5 per cent. We discuss how the purity of the sample affects our results and how supervised learning methods may help to discriminate between CC and PI SNe.

preprint2015arXiv

The assembly of "normal" galaxies at z=7 probed by ALMA

We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5-15 Msun/yr these systems are much more representative of the high-z galaxy population than other systems targeted in the past by millimeter observations. For the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected in any of the three galaxies, with 3sigma upper limits significantly lower than the [CII] emission observed in lower reshift galaxies. These results suggest that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. These findings are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.

preprint2015arXiv

The Brightest Ly$α$ Emitter: Pop III or Black Hole?

CR7 is the brightest $z=6.6 \, {\rm Ly}α$ emitter (LAE) known to date, and spectroscopic follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III stars. We examine this interpretation using cosmological hydrodynamical simulations. Several simulated galaxies show the same "Pop III wave" pattern observed in CR7. However, to reproduce the extreme CR7 ${\rm Ly}α$/HeII1640 line luminosities ($L_{\rm α/He II}$) a top-heavy IMF and a massive ($>10^{7}{\rm M}_{\odot}$) PopIII burst with age $<2$ Myr are required. Assuming that the observed properties of ${\rm Ly}α$ and HeII emission are typical for Pop III, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs at $z=6.6$ with $L_α >10^{43.3}{\rm erg}\,{\rm s}^{-1}$ should also host Pop III stars producing an observable $L_{\rm He II}>10^{42.7}{\rm erg}\,{\rm s}^{-1}$. As an alternate explanation, we explore the possibility that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). Our model predicts $L_α$, $L_{\rm He II}$, and X-ray luminosities that are in agreement with the observations. In any case, the observed properties of CR7 indicate that this galaxy is most likely powered by sources formed from pristine gas. We propose that further X-ray observations can distinguish between the two above scenarios.

preprint2014arXiv

Initial mass function of intermediate mass black hole seeds

We study the Initial Mass Function (IMF) and host halo properties of Intermediate Mass Black Holes (IMBH, 10^{4-6} Msun) formed inside metal-free, UV illuminated atomic cooling haloes (virial temperature T_vir > 10^4 K) either via the direct collapse of the gas or via an intermediate Super Massive Star (SMS) stage. We achieve this goal in three steps: (a) we derive the gas accretion rate for a proto-SMS to undergo General Relativity instability and produce a direct collapse black hole (DCBH) or to enter the ZAMS and later collapse into a IMBH; (b) we use merger-tree simulations to select atomic cooling halos in which either a DCBH or SMS can form and grow, accounting for metal enrichment and major mergers that halt the growth of the proto-SMS by gas fragmentation. We derive the properties of the host halos and the mass distribution of black holes at this stage, and dub it the "Birth Mass Function"; (c) we follow the further growth of the DCBH due to accretion of leftover gas in the parent halo and compute the final IMBH mass.We consider two extreme cases in which minihalos (T_vir < 10^4 K) can (fertile) or cannot (sterile) form stars and pollute their gas leading to a different IMBH IMF. In the (fiducial) fertile case the IMF is bimodal extending over a broad range of masses, M= (0.5-20)x10^5 Msun, and the DCBH accretion phase lasts from 10 to 100 Myr. If minihalos are sterile, the IMF spans the narrower mass range M= (1-2.8)x10^6 Msun, and the DCBH accretion phase is more extended (70-120 Myr). We conclude that a good seeding prescription is to populate halos (a) of mass 7.5 < log (M_h/Msun) < 8, (b) in the redshift range 8 < z < 17, (c) with IMBH in the mass range 4.75 < log (M_BH/Msun) < 6.25.

preprint2014arXiv

Simulating cosmic metal enrichment by the first galaxies

We study cosmic metal enrichment via AMR hydrodynamical simulations in a (10 Mpc/h)$^3$ volume following the Pop III-Pop II transition and for different Pop III IMFs. We have analyzed the joint evolution of metal enrichment on galactic and intergalactic scales at z=6 and z=4. Galaxies account for <9% of the baryonic mass; the remaining gas resides in the diffuse phases: (a) voids, i.e. regions with extremely low density ($Δ$<1), (b) the true intergalactic medium (IGM, 1<$Δ$<10) and (c) the circumgalactic medium (CGM, 10<$Δ<10^{2.5}$), the interface between the IGM and galaxies. By z=6 a galactic mass-metallicity relation is established. At z=4, galaxies with a stellar mass $M_*=10^{8.5}M_\odot$ show log(O/H)+12=8.19, consistent with observations. The total amount of heavy elements rises from $Ω^{SFH}_Z=1.52\, 10^{-6}$ at z=6 to 8.05 $10^{-6}$ at z=4. Metals in galaxies make up to ~0.89 of such budget at z=6; this fraction increases to ~0.95 at z=4. At z=6 (z=4) the remaining metals are distributed in CGM/IGM/voids with the following mass fractions: 0.06/0.04/0.01 (0.03/0.02/0.01). Analogously to galaxies, at z=4 a density-metallicity ($Δ$-Z) relation is in place for the diffuse phases: the IGM/voids have a spatially uniform metallicity, Z~$10^{-3.5}$Zsun; in the CGM Z steeply rises with density up to ~$10^{-2}$Zsun. In all diffuse phases a considerable fraction of metals is in a warm/hot (T>$10^{4.5}$K) state. Due to these physical conditions, CIV absorption line experiments can probe only ~2% of the total carbon present in the IGM/CGM; however, metal absorption line spectra are very effective tools to study reionization. Finally, the Pop III star formation history is almost insensitive to the chosen Pop III IMF. Pop III stars are preferentially formed in truly pristine (Z=0) gas pockets, well outside polluted regions created by previous star formation episodes.

preprint2014arXiv

Very extended cold gas, star formation and outflows in the halo of a bright QSO at z>6

Past observations of QSO host galaxies at z >6 have found cold gas and star formation on compact scales of a few kiloparsecs. We present new high sensitivity IRAM PdBI follow-up observations of the [CII] 158micron emission line and FIR continuum in the host galaxy of SDSS J1148+5152, a luminous QSO at redshift 6.4189. We find that a large fraction of the gas traced by [CII] is at high velocities, up to ~1400 km/s relative to the systemic velocity, confirming the presence of a major quasar-driven outflow indicated by previous observations. The outflow has a complex morphology and reaches a maximum projected radius of ~30 kpc. The extreme spatial extent of the outflow allows us, for the first time in an external galaxy, to estimate mass-loss rate, kinetic power and momentum rate of the outflow as a function of the projected distance from the nucleus and the dynamical time-scale. These trends reveal multiple outflow events during the past 100 Myr, although the bulk of the mass, energy and momentum appear to have been released more recently, within the past ~20 Myr. Surprisingly, we discover that also the quiescent gas at systemic velocity is extremely extended. More specifically, we find that, while 30% of the [CII] within v\in(-200, 200) km/s traces a compact component that is not resolved by our observations, 70% of the [CII] emission in this velocity range is extended, with a projected FWHM size of 17.4+-1.4 kpc. We detect FIR continuum emission associated with both the compact and the extended [CII] components, although the extended FIR emission has a FWHM of 11+-3 kpc, thus smaller than the extended [CII] source. Overall, our results indicate that the cold gas traced by [CII] is distributed up to r~30 kpc. A large fraction of extended [CII] is likely associated with star formation on large scales, but the [CII] source extends well beyond the FIR continuum.

preprint2013arXiv

Can supermassive black hole seeds form in galaxy mergers?

It has been recently suggested that supermassive black holes at z = 5-6 might form from super-fast (\dot M > 10^4 Msun/yr) accretion occurring in unstable, massive nuclear gas disks produced by mergers of Milky-Way size galaxies. Interestingly, such mechanism is claimed to work also for gas enriched to solar metallicity. These results are based on an idealized polytropic equation of state assumption, essentially preventing the gas from cooling. We show that under more realistic conditions, the disk rapidly (< 1 yr) cools, the accretion rate drops, and the central core can grow only to \approx 100 Msun. In addition, most of the disk becomes gravitationally unstable in about 100 yr, further quenching the accretion. We conclude that this scenario encounters a number of difficulties that possibly make it untenable.

preprint2013arXiv

Light from the Cosmic Frontier: Gamma-Ray Bursts

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

preprint2013arXiv

Metals and ionizing photons from dwarf galaxies

We estimate the potential contribution of M < 10^9 Msun dwarf galaxies to the reionization and early metal-enrichment of the Milky Way environment, or circum-Galactic Medium. Our approach is to use the observed properties of ancient stars (> 12 Gyr old) measured in nearby dwarf galaxies to characterize the star-formation at high-z. We use a merger-tree model for the build-up of the Milky Way, which self-consistently accounts for feedback processes, and which is calibrated to match the present-day properties of the Galaxy and its dwarf satellites. We show that the high-z analogues of nearby dwarf galaxies can produce the bulk of ionizing radiation (>80%) required to reionize the Milky Way environment. Our fiducial model shows that the gaseous environment can be 50% reionized at z ~ 8 by galaxies with 10^7 Msun < M < 10^8 Msun. At later times, radiative feedback stops the star-formation in these small systems, and reionization is completed by more massive dwarf galaxies by z_rei = 6.4\pm 0.5. The metals ejected by supernova-driven outflows from M < 10^9 Msun dwarf galaxies almost uniformly fill the Milky Way environment by z ~ 5, enriching it to Z ~ 2 10^-2 Zsun. At z ~ 2 these early metals are still found to represent ~ 50% of the total mass of heavy elements in the circum-Galactic Medium.

preprint2013arXiv

Simulating intergalactic quasar scintillation

Intergalactic scintillation of distant quasars is sensitive to free electrons and therefore complements Ly$α$ absorption line experiments probing the neutral intergalactic medium (IGM). We present a new scheme to compute IGM refractive scintillation effects on distant sources in combination with Adaptive Mesh Refinement cosmological simulations. First we validate our model by reproducing the well-known interstellar scintillation (ISS) of Galactic sources. The simulated cosmic density field is then used to infer the statistical properties of intergalactic scintillation. Contrary to previous claims, we find that the scattering measure of the simulated IGM at $z<2$ is $\langle \mbox{SM}_{\equ}\rangle=3.879$, i.e. almost 40 times larger than for the usually assumed smooth IGM. This yield an average modulation index ranging from 0.01 ($ν_s=5$ GHz) up to 0.2 ($ν_s=50$ GHz); above $ν_{s}\gsim30$ GHz the IGM contribution dominates over ISS modulation. We compare our model with data from a $0.3\leq z\leq 2$ quasar sample observed at $ν_{\obs}=8.4$ GHz. For this high frequency ($10.92\leq ν_s \leq 25.2$), high galactic latitude sample ISS is negligible, and IGM scintillation can reproduce the observed modulation with a 4% accuracy, without invoking intrinsic source variability. We conclude by discussing the possibility of using IGM scintillation as a tool to pinpoint the presence of intervening high-$z$ groups/clusters along the line of sight, thus making it a probe suitably complementing Sunyaev-Zeldovich data recently obtained by \textit{Planck}.

preprint2012arXiv

Escape Fraction of Ionizing Radiation from Starburst Galaxies at High Redshifts

Recent data indicates that the cosmic UV emissivity decreased with decreasing redshift z near the end of reionization. Lacking evidence for very massive early stars, this could signal a decline with time in the mass-averaged escape fraction of ionizing radiation from galaxies <fesc> at z > 6. We calculate the evolution of ionization fronts in dark matter halos which host gas in hydrostatic equilibrium at its cooling temperature floor (T~10^4 K for atomic hydrogen). We find a high escape fraction only for the lowest mass halos (with M< 10^8.7 Msun at (1+z)=10) provided their star formation efficiency f_star > 10^-3. Since the low-mass galaxy population is depleted by radiative feedback, we find that indeed <fesc> decreases with time during reionization.

preprint2012arXiv

Evidence of strong quasar feedback in the early Universe

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

preprint2012arXiv

Joint Lyα emitters - quasars reionization constraints

We present a novel method to investigate cosmic reionization, using joint spectral information on high redshift Lyman Alpha Emitters (LAE) and quasars (QSOs). Although LAEs have been proposed as reionization probes, their use is hampered by the fact their Lyα line is damped not only by intergalactic HI but also internally by dust. Our method allows to overcome such degeneracy. First, we carefully calibrate a reionization simulation with QSO absorption line experiments. Then we identify LAEs in two simulation boxes at z=5.7 and z=6.6 and we build synthetic images/spectra of a prototypical LAE. At redshift 5.7, we find that the Lyα transmissivity (T_LAE) ~ 0.25, almost independent of the halo mass. This constancy arises from the conspiracy of two effects: (i) the intrinsic Lyα line width and (ii) the infall peculiar velocity. At higher redshift, z=6.6, where the transmissivity is instead largely set by the local HI abundance and LAE transmissivity consequently increases with halo mass from 0.15 to 0.3. Although outflows are present, they are efficiently pressure-confined by infall in a small region around the LAE; hence they only marginally affect transmissivity. Finally, we cast LOS originating from background QSOs passing through foreground LAEs at different impact parameters, and compute the quasar transmissivity (T_QSO). At smaller impact parameters, d < 1 cMpc, a positive correlation between T_QSO and halo mass is found at z = 5.7, which tends to become less pronounced (i.e. flatter) at larger distances. By cross-correlating T_LAE and T_QSO, we can obtain a HI density estimate unaffected by dust. At z= 5.7, the cross-correlation is relatively weak,whereas at z = 6.6 we find a clear positive correlation. We conclude by briefly discussing the perspectives for the application of the method to existing and forthcoming data.

preprint2011arXiv

Simulating high-redshift galaxies

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, α~-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.

preprint2010arXiv

First star formation with dark matter annihilation

We include an energy term based on Dark Matter (DM) self-annihilation during the cooling and subsequent collapse of the metal-free gas, in halos hosting the formation of the first stars in the Universe. We have found that the feedback induced on the chemistry of the cloud does modify the properties of the gas throughout the collapse. However, the modifications are not dramatic, and the typical Jeans mass within the halo is conserved throughout the collapse, for all the DM parameters we have considered. This result implies that the presence of Dark Matter annihilations does not substantially modify the Initial Mass Function of the First Stars, with respect to the standard case in which such additional energy term is not taken into account. We have also found that when the rate of energy produced by the DM annihilations and absorbed by the gas equals the chemical cooling (at densities yet far from the actual formation of a proto-stellar core) the structure does not halt its collapse, although that proceeds more slowly by a factor smaller than few per cent of the total collapse time.

preprint2010arXiv

High redshift Lya emitters: clues on the Milky Way infancy

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.

preprint2010arXiv

Mining the Galactic Halo for Very Metal-Poor Stars

We study the age and metallicity distribution function (MDF) of metal-poor stars in the Milky Way halo as a function of galactocentric radius by combining N-body simulations and semi-analytical methods. We find that the oldest stars populate the innermost region, while extremely metal-poor stars are more concentrated within r < 60 kpc. The MDF of [Fe/H] < -2 stars varies only very weakly within the central 50 kpc, while the relative contribution of [Fe/H] < -2 stars strongly increases with r, varying from 16% within 7 kpc < r < 20 kpc up to > 40% for r > 20 kpc. This is due to the faster descent of the spatial distribution (as seen from Earth) of the more enriched population. This implies that the outer halo < 40 kpc is the best region to search for very metal-poor stars. Beyond ~ 60 kpc the density of [Fe/H] < -2 stars is maximum within dwarf galaxies. All these features are imprinted by a combination of (i) the virialization epoch of the star-forming haloes, and (ii) the metal enrichment history of the Milky Way environment.

preprint2009arXiv

Lyman Alpha Emitter Evolution in the Reionization Epoch

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.

preprint2009arXiv

WIMP annihilation effects on primordial star formation

We study the effects of WIMP dark matter (DM) annihilations on the thermal and chemical evolution of the gaseous clouds where the first generation of stars in the Universe is formed. We follow the collapse of the gas inside a typical halo virializing at very high redshift, from well before virialization until a stage where the heating from DM annihilations exceeds the gas cooling rate. The DM energy input is estimated by inserting the energy released by DM annihilations (as predicted by an adiabatic contraction of the original DM profile) in a spherically symmetric radiative transfer scheme. In addition to the heating effects of the energy absorbed, we include its feedback upon the chemical properties of the gas, which is critical to determine the cooling rate in the halo, and hence the fragmentation scale and Jeans mass of the first stars. We find that DM annihilation does alter the free electron and especially the H2 fraction when the gas density is n>~ 10^4 cm^-3, for our fiducial parameter values. However, even if the change in the H2 abundance and the cooling efficiency of the gas is large (sometimes exceeding a factor 100), the effects on the temperature of the collapsing gas are far smaller (a reduction by a factor <~1.5), since the gas cooling rate depends very strongly on temperature: then, the fragmentation mass scale is reduced only slightly, hinting towards no dramatic change in the initial mass function of the first stars.

preprint2008arXiv

Dark matter annihilation effects on the first stars

We study the effects of WIMP dark matter (DM) on the collapse and evolution of the first stars in the Universe. Using a stellar evolution code, we follow the pre-Main Sequence (MS) phase of a grid of metal-free stars with masses in the range 5-600 solar mass forming in the centre of a 1e6 solar mass halo at redhisft z=20. DM particles of the parent halo are accreted in the proto-stellar interior by adiabatic contraction and scattering/capture processes, reaching central densities of order 1e12 GeV/cm3 at radii of the order of 10 AU. Energy release from annihilation reactions can effectively counteract the gravitational collapse, in agreement with results from other groups. We find this stalling phase (known as "dark" star) is transients and lasts from 2.1e3 yr (M=600 solar mass) to 1.8e4 yr (M=9 solar mass). Later in the evolution, DM scattering/capture rate becomes high enough that energy deposition from annihilations significantly alters the pre-MS evolution of the star in a way that depends on DM (i) velocity dispersion, (ii) density, (iii) elastic scattering cross section with baryons. For our fiducial set of parameters (10 km/s, 1e11 GeV/cm3, 1e-38 cm2) we find that the evolution of stars of mass lower than 40 solar masses "freezes" on the HR diagram before reaching the ZAMS. Stars with bigger masses manage to ignite nuclear reactions; however, DM "burning" prolonges their lifetimes by a factor 2 (5) for a 600 (40) solar mass star.

preprint2008arXiv

Was the Universe neutral beyond redshift six?

We provide measurements of the neutral hydrogen fraction xHI at z~6, by comparing semi-analytical models of the Lyalpha forest with observations of high-z quasars and Gamma Ray Bursts absorption spectra. We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.74<zem<6.42 studying separately the narrow transmission windows (peaks) and the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that xHI evolves smoothly from 10^{-4.4} at z=5.3 to 10^{-4.2} at z=5.6, with a robust upper limit xHI<0.36 at z=6.3. We show the results of the first-ever detected transverse proximity effect in the HI Lyalpha forest, produced by the HII region of the faint quasar RD J1148+5253 at z=5.70 intervening along the LOS of SDSS J1148+5251 at z=6.42. Moreover, we propose a novel method to study cosmic reionization using absorption line spectra of high-redshift GRBs afterglows. We show that the time evolution and the statistics of gaps in the observed spectra represent exquisite tools to discriminate among different reionization models. By applying our methods to GRB 050904 detected at z=6.29, we show that the observation of this burst provides strong indications of a highly ionized intergalactic medium at z~6, with an estimated mean neutral hydrogen fraction xHI=6.4\pm 0.3\times 10^{-5} along that line of sight.

preprint2007arXiv

CMB polarization constraints on radiative feedback

We compute the imprints left on the CMB by two cosmic reionization models consistent with current observations but characterized by alternative radiative feedback prescriptions (suppression and filtering) resulting in a different suppression of star formation in low-mass halos. The models imply different ionization and thermal histories and 21 cm background signals. The derived Comptonization, u, and free-free distortion, y_B, parameters are below current observational limits for both models. However, the value of u = 1.69 * 10^-7 (9.65 * 10^-8) for the suppression (filtering) model is in the detectability range of the next generation of CMB spectrum experiments. Through the dedicated Boltzmann code CMBFAST, modified to include the above ionization histories, we compute the CMB angular power spectrum (APS) of the TT, TE, and EE modes. For the EE mode the differences between these models are significantly larger than the cosmic and sampling variance over the multipole range l = 5-15, leaving a good chance of discriminating between these feedback mechanisms with forthcoming/future CMB polarization experiments. The main limitations come from foreground contamination: it should be subtracted at per cent level in terms of APS, a result potentially achievable by novel component separation techniques and mapping of Galactic foreground.

preprint2003arXiv

Low-Mass Relics of Early Star Formation

The earliest stars to form in the Universe were the first sources of light, heat and metals after the Big Bang. The products of their evolution will have had a profound impact on subsequent generations of stars. Recent studies of primordial star formation have shown that, in the absence of metals (elements heavier than helium), the formation of stars with masses 100 times that of the Sun would have been strongly favoured, and that low-mass stars could not have formed before a minimum level of metal enrichment had been reached. The value of this minimum level is very uncertain, but is likely to be between 10^{-6} and 10^{-4} that of the Sun. Here we show that the recent discovery of the most iron-poor star known indicates the presence of dust in extremely low-metallicity gas, and that this dust is crucial for the formation of lower-mass second-generation stars that could survive until today. The dust provides a pathway for cooling the gas that leads to fragmentation of the precursor molecular cloud into smaller clumps, which become the lower-mass stars.

preprint1996arXiv

Properties of the Lyman Alpha Clouds from non-equilibrium Photoionization Models

We investigate the thermal and ionization history of Lyman alpha clouds photoionized by a time--dependent UV background, including non equilibrium effects. The results show that it is possible to obtain temperatures as low as T~15000 K (or, equivalently, Doppler parameters b~15 km/s) at z=3 for cloud total densities n~10^-4 cm^-3, if (i) the reionization epoch occurred at z_i~10, and (ii) the UV background has a factor 70-100 decrease at the HeII edge. A trend towards smaller b with increasing redshift is present in the redshift interval z=1-5. Higher densities lead to higher values of b and smaller hydrogen correction factors, n_{HII}/n_{HI}. The correction factors for helium are also given. For a hydrogen column density N_{HI}=3 10^{14} cm^{-2}, cloud sizes are larger than 100 kpc, consistent with recent observations of quasar pairs. Pressure confined models, instead, yield implausibly low cloud densities at low redshift, and too small sizes at intermediate redshift. The implications of the model are confronted with the available observational data.