Source author record

Andrea Pallottini

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

14works
2topics
4close collaborators

Published work

14 publications

preprint2026arXiv

Amaryllis: a digital twin of the earliest galaxies in the Universe

Synergies between JWST and ALMA are unveiling a population of bright, super-early (z>10z>10) galaxies, including systems like GS-z14-0 (z=14.2z=14.2) and GHZ2 (z=12.3z=12.3) with extreme FIR line ratios ([OIII] 88 \,um / [CII] 158 \,um >3>3) that challenge galaxy formation models. To address this, we identify a synthetic analog of these sources, "Amaryllis", within the SERRA zoom-in simulations, and track its evolution from z=16z=16 to z=7z=7. During this period, Amaryllis grows from log⁡(M⋆/M⊙)∼7.4\log(M_\star/M_{\odot}) \sim 7.4 to 10.310.3, linking super-early progenitors to the massive galaxy population at the end of reionization. At z∼11.3z \sim 11.3, Amaryllis closely matches the observed properties of GS-z14-0, including M⋆M_\star, SFR, and the luminosity of FIR ([OIII] 88 \,um) and UV (e.g. CIII] 1908\,1908) lines. We find that high [OIII]/[CII] ratios appear during short, merger-driven starburst episodes, when low metallicity (Z∼0.1 Z⊙Z \sim 0.1\,Z_{\odot}) and high ionization conditions (Uion∼0.3U_{\mathrm{ion}} \sim 0.3) push the ISM far from equilibrium. These extreme FIR line ratios are thus transient and linked to major mergers that ignite strong ionized gas outflows. Strikingly, despite this dynamical activity, Amaryllis develops a rotation-supported gaseous disk (V/σ∼4V/σ\sim 4-6) by z∼11z \sim 11, while stars remain dispersion-dominated. This coexistence of ordered gas rotation and merger-driven disturbances occurs within a massive yet typical ΛΛCDM halo, enabling disk formation even at cosmic dawn.

preprint2022arXiv

Hydrogen reionisation ends by z=5.3z=5.3: Lyman-αα optical depth measured by the XQR-30 sample

The presence of excess scatter in the Ly-αα forest at z∼5.5z\sim 5.5, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-αα transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at z>5.5z>5.5 with high signal-to-noise ratios of >10>10 per ≤15\leq 15 km s−1^{-1} spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR >20>20 per pixel) spectra of 30 quasars at z>5.7z>5.7. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-αα systems. We present improved measurements of the mean Ly-αα transmission over 4.9<z<6.14.9<z<6.1. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-αα transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to z≤5.2z\leq5.2 (<1σ<1 σ), and mild tension (∼2.5σ\sim 2.5 σ) at z=5.3z=5.3. Homogeneous UVB models are ruled out by excess Ly-αα transmission scatter at z≥5.4z\geq5.4 with high confidence (>3.5σ>3.5 σ). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least z=5.3z=5.3 (t=1.1t=1.1 Gyr after the Big Bang). This is further evidence for a late end to reionisation.

preprint2022arXiv

Long Dark Gaps in the Lyββ Forest at z<6z<6: Evidence of Ultra Late Reionization from XQR-30 Spectra

We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-ββ (Lyββ) forest. With its lower optical depth, Lyββ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Lyαα line. We identify dark gaps in the Lyββ forest using spectra of 42 QSOs at zem>5.5z_{\rm em}>5.5, including new data from the XQR-30 VLT Large Programme. Approximately 40%40\% of these QSO spectra exhibit dark gaps longer than 10h−1Mpc10h^{-1}{\rm Mpc} at z≃5.8z\simeq5.8. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Lyαα forest are caused solely by ionizing ultraviolet background (UVB) fluctuations and with models that include large neutral hydrogen patches at z<6z<6 due to a late end to reionization. Of particular interest is a very long (L=28h−1MpcL=28h^{-1}{\rm Mpc}) and dark (τeff≳6τ_{\rm eff} \gtrsim 6) gap persisting down to z≃5.5z\simeq 5.5 in the Lyββ forest of the z=5.85z_{\rm}=5.85 QSO PSO J025−-11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of ⟨xHI⟩≳5% \langle x_{\rm HI}\rangle \gtrsim 5\% by z=5.6z=5.6. Finally, we infer constraints on ⟨xHI⟩\langle x_{\rm HI}\rangle over 5.5≲z≲6.05.5 \lesssim z \lesssim 6.0 based on the observed Lyββ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find ⟨xHI⟩≤0.05\langle x_{\rm HI}\rangle \leq 0.05, 0.17, and 0.29 at z≃5.55z\simeq 5.55, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than z=6z = 6.

preprint2022arXiv

The 21-cm signal from the Cosmic Dawn: metallicity dependence of high mass X-ray binaries

X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; z∼10−15z\sim 10-15) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs (LX/SFRL_{\rm X}/{\rm SFR}) should increase in metal-poor environments, typical of early galaxies. Here we study the impact of the metallicity (ZZ) dependence of LX/SFRL_{\rm X}/{\rm SFR} during the EoH. For our fiducial models, galaxies with star formation rates of order 10−3−10−110^{-3} - 10^{-1} M⊙M_\odot yr−1^{-1} and metallicities of order 10−3−10−210^{-3} - 10^{-2} Z⊙Z_\odot are the dominant contributors to the X-ray background (XRB) during this period. Different LX/SFRL_{\rm X}/{\rm SFR}-ZZ relations result in factors of ∼\sim 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000h integration with the upcoming Square Kilometer Array (SKA), for two different LX/SFRL_{\rm X}/{\rm SFR}-ZZ relations. We perform inference on these mock observations using the common simplification of a constant LX/SFRL_{\rm X}/{\rm SFR}, finding that constant LX/SFRL_{\rm X}/{\rm SFR} models can recover the IGM evolution of the more complicated LX/SFRL_{\rm X}/{\rm SFR}-ZZ simulations only during the EoH. At z<10z<10, where the typical galaxies are more polluted, constant LX/SFRL_{\rm X}/{\rm SFR} models over-predict the XRB and its relative contribution to the early stages of the reionization.

preprint2022arXiv

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

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

preprint2022arXiv

The interstellar medium of high-redshift galaxies: Gathering clues from C III] and [C II] lines

A tight relation between [C II] line luminosity and the star formation rate (SFR) has been observed for local galaxies. At high redshift (z > 5), galaxies instead deviate downwards from the local ΣΣ_[C II] - ΣΣ_SFR relation. This deviation might be caused by different interstellar medium (ISM) properties in galaxies at early epochs. To test this hypothesis, we combined the [C II] and SFR data with C III] line observations and our physical models. We additionally investigated how ISM properties, such as burstiness, κsκ_s, total gas density, nn, and metallicity, ZZ, affect the deviation from the ΣΣ_[C II] - ΣΣ_SFR relation in these sources. We present the VLT/X-SHOOTER observations targeting the C III] λ1909λ1909 line emission in three galaxies at 5.5 < z < 7. We include X-SHOOTER and VLT/MUSE archival data of eight galaxies at 2 < z < 7, and eleven star-forming systems at 6 < z < 7.5, with either C III] or [C II] detection reported in the literature. We detected C III] λλ1907,1909λλ1907, 1909 line emission in HZ10 and we derived the intrinsic, integrated flux of the C III] λ1909λ1909 line. We constrained the ISM properties for our sample of galaxies, κsκ_s , nn, and ZZ, by applying our physically motivated model based on the MCMC algorithm. For the most part, high-z star-forming galaxies show subsolar metallicities. The majority of the sources have log(κs)>1log(κ_s) > 1, that is, they overshoot the Kennicutt-Schmidt (KS) relation by about one order of magnitude, implying that the whole KS relation might be shifted upwards at early times. Furthermore, all the high-z galaxies of our sample lie below the ΣΣ_[C II] - ΣΣ_SFR local relation. The total gas density, nn, shows the strongest correlation with the deviation from the local ΣΣ_[C II] - ΣΣ_SFR relation.

preprint2021arXiv

Measuring the Density Fields around Bright Quasars at z∼6z\sim 6 with XQR-30 Spectra

Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-αα absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-αα absorption is not saturated. In this study, we use 1010 high-resolution (R≳10,000R\gtrsim 10,000) z∼6z\sim 6 quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within 33 pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the Cosmic Reionization on Computers\textit{Cosmic Reionization on Computers} simulation, and find a good agreement between 1.51.5 to 33pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter σ8=0.6±0.3σ_8=0.6\pm0.3. The uncertainty is mainly due to the limited number of high-quality quasar sightlines currently available. Utilizing the more than >200>200 known quasars at z≳6z\gtrsim 6, this method will allow us in the future to tighten the constraint on σ8σ_8 to the percent level. In the region closer to the quasar within 1.51.5 pMpc, we find the density is higher than predicted in the simulation by 1.23±0.171.23 \pm 0.17, suggesting the typical host dark matter halo mass of a bright quasar (M1450<−26.5M_{\rm 1450}<-26.5) at z∼6z\sim 6 is log⁡10(M/M⊙)=12.5−0.7+0.4\log_{\rm 10} (M/M_\odot)=12.5^{+0.4}_{-0.7}.

preprint2020arXiv

Dynamical Properties of Molecular-forming Gas Clumps in Galaxies at the Epoch of Reionization

We study the properties of molecular-forming gas clumps (MGCs) at the epoch of reionization using cosmological zoom-in simulations. We identify MGCs in a z=6 prototypical galaxy ("Althaea") using an H2 density-based clump finder. We compare their mass, size, velocity dispersion, gas surface density, and virial parameter (alpha_vir) to observations. In Althaea, the typical MGC mass and size are Mgas=10^6.5 Msun and R=45-100 pc, which are comparable to those found in nearby spirals and starburst galaxies. MGCs are highly supersonic and supported by turbulence, with rms velocity dispersions of sigma_gas=20-100 km s^-1 and pressure of P/k_B=10^7.6 K cm^-3 (i.e., >1000x with respect to the Milky Way), similar to those found in nearby and z~2 gas-rich starburst galaxies. In addition, we perform stability analysis to understand the origin and dynamical properties of MGCs. We find that MGCs are globally stable in the main disk of Althaea. Densest regions where star formation is expected to take place in clouds and cores on even smaller scales instead have lower alpha_vir and Toomre-Q values. Detailed studies of the star-forming gas dynamics at the epoch of reionization thus require a spatial resolution of <40 pc (=0.01"), which is within reach with the Atacama Large (sub-)Millimeter Array and the Next Generation Very Large Array.

preprint2020arXiv

Large Population of ALMA Galaxies at z>6 with Very High [OIII]88um to [CII]158um Flux Ratios: Evidence of Extremely High Ionization Parameter or PDR Deficit?

We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

preprint2020arXiv

Outflows and extended [CII] halos in high redshift galaxies

Recent stacked ALMA observations have revealed that normal, star-forming galaxies at z≈6z\approx 6 are surrounded by extended (≈10 kpc\approx 10\,\mathrm{kpc}) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading factor, ηη, and the parent galaxy dark matter halo circular velocity, vcv_c. The outflow model successfully matches the observed [CII] surface brightness profile if η=3.20±0.10η= 3.20 \pm 0.10 and vc=170±10 km s−1v_c = 170 \pm 10{\,\rm km\,s^{-1}}, corresponding to a dynamical mass of ≈1011 M⊙\approx 10^{11}\, \mathrm{M}_\odot. The predicted outflow rate and velocity range are 128±5 M⊙yr−1128 \pm 5 \,\mathrm{M}_\odot {\rm yr}^{-1} and 300−500 km s−1300-500 {\,\rm km\,s^{-1}}, respectively. We conclude that: (a) extended halos can be produced by cooling outflows; (b) the large ηη value is marginally consistent with starburst-driven outflows, but it might indicate additional energy input from AGN; (c) the presence of [CII] halos requires an ionizing photon escape fraction from galaxies fesc≪1f_{\rm esc} \ll 1. The model can be readily applied also to individual high-zz galaxies, as those observed, e.g., by the ALMA ALPINE survey now becoming available.

preprint2019arXiv

First Identification of 10-kpc Scale [CII] 158um Halos around Star-Forming Galaxies at z=5-7

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

preprint2016arXiv

Inferring physical properties of galaxies from their emission line spectra

We present a new approach based on Supervised Machine Learning (SML) algorithms to infer key physical properties of galaxies (density, metallicity, column density and ionization parameter) from their emission line spectra. We introduce a numerical code (called GAME, GAlaxy Machine learning for Emission lines) implementing this method and test it extensively. GAME delivers excellent predictive performances, especially for estimates of metallicity and column densities. We compare GAME with the most widely used diagnostics (e.g. R23_{23}, [NII]λλ6584 / Hαα indicators) showing that it provides much better accuracy and wider applicability range. GAME is particularly suitable for use in combination with Integral Field Unit (IFU) spectroscopy, both for rest-frame optical/UV nebular lines and far-infrared/sub-mm lines arising from Photo-Dissociation Regions. Finally, GAME can also be applied to the analysis of synthetic galaxy maps built from numerical simulations.

preprint2015arXiv

Intensity mapping of [CII] from early galaxies

The intensity mapping of the [CII] 157.7 μ\rm μm fine-structure emission line represents an ideal experiment to probe star formation activity in galaxies, especially in those that are too faint to be individually detected. Here, we investigate the feasibility of such an experiment for z>5z > 5 galaxies. We construct the LCII−MhL_{\rm CII} - M_{\rm h} relation from observations and simulations, then generate mock [CII] intensity maps by applying this relation to halo catalogs built from large scale N-body simulations. Maps of the extragalactic far-infrared (FIR) continuum, referred to as "foreground", and CO rotational transition lines and [CI] fine-structure lines referred to as "contamination", are produced as well. We find that, at 316 GHz (corresponding to zCII=5z_{\rm CII} = 5), the mean intensities of the extragalactic FIR continuum, [CII] signal, all CO lines from J=1J=1 to 13 and two [CI] lines are ∼3×105\sim 3\times10^5 Jy sr−1^{-1}, ∼1200\sim 1200 Jy sr−1^{-1}, ∼800\sim 800 Jy sr−1^{-1} and ∼100\sim 100 Jy sr−1^{-1}, respectively. We discuss a method that allows us to subtract the FIR continuum foreground by removing a spectrally smooth component from each line of sight, and to suppress the CO/[CI] contamination by discarding pixels that are bright in contamination emission. The z>5z > 5 [CII] signal comes mainly from halos in the mass range 1011−12 M⊙10^{11-12} \,M_\odot; as this mass range is narrow, intensity mapping is an ideal experiment to investigate these early galaxies. In principle such signal is accessible to a ground-based telescope with a 6 m aperture, 150 K system temperature, a 128×128128\times128 pixels FIR camera in 5000 hr total integration time, however it is difficult to perform such an experiment by using currently available telescopes.

preprint2014arXiv

The circumgalactic medium of high redshift galaxies

We study the properties of the circumgalactic medium (CGM) of high-zz galaxies in the metal enrichment simulations presented in Pallottini et al. 2014. At z=4z=4, we find that the simulated CGM gas density profiles are self-similar, once scaled with the virial radius of the parent dark matter halo. We also find a simple analytical expression relating the neutral hydrogen equivalent width (EWHI{\rm EW}_{\rm HI}) of CGM absorbers as a function of the line of sight impact parameter (bb). We test our predictions against mock spectra extracted from the simulations, and show that the model reproduces the EWHI(b){\rm EW}_{\rm HI}(b) profile extracted from the synthetic spectra analysis. When compared with available data, our CGM model nicely predicts the observed EWHI(b){\rm EW}_{\rm HI}(b) in z≲2z\lesssim2 galaxies, and supports the idea that the CGM profile does not evolve with redshift.

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.