Source author record

Umberto Maio

Umberto Maio 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

26works
7topics
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

26 published item(s)

preprint2025arXiv

High$-z$ [OI] emission lines: ColdSIM simulations and ALMA observations

Neutral-oxygen [OI] far-infrared emission lines at $63μ$m and $145μ$m are powerful probes of the physical conditions in the interstellar medium, although they have not been fully exploited in high-redshift studies. We investigate the connection between [OI] emission lines and key galaxy properties, such as star formation rate (SFR) and H$_2$ content. Our predictions are compared with existing observations and new data analysed in this work. We post-process the outputs of the ColdSIM cosmological simulations with the DESPOTIC model, taking into account [OI]$63μ$m self-absorption by cold foreground material. A Random Forest algorithm is employed to accelerate computations and new observational ALMA data for galaxies at redshift $z\simeq 5-7$ are used to validate our model. Our predictions show significant [OI]$63μ$m luminosities ($\approx 10^8\,\rm L_\odot$) for galaxies with SFRs of $\approx 10^2\,\rm M_\odot\,{\rm yr}^{-1}$. The $145μ$m line luminosity is typically $15 \%$ the [OI]$63μ$m one and is a factor $\approx 2-20$ below high-$z$ observations. Both [OI] lines correlate with SFR and molecular mass, but exhibit flattening in scaling relations with metallicity and stellar mass. Foreground self-absorption reduces the [OI] flux by a factor of $2-4$, consistent with empirical corrections in observational studies. We find typical line ratios of [OI]$63μ$m / [CII]$158μ$m $\approx 1$ and [OI]$145μ$m / [CII]$158μ$m $ \approx 0.2 \, -$ consistent with $z\gtrsim 6$ observations, but only when [OI]$63μ$m self-absorption is included. Both [OI]$63μ$m and [OI]$145μ$m lines serve as tracers of star formation and molecular gas at high redshift. Their joint detection can provide constraints on the properties of the early interstellar medium and self-absorption of the [OI]$63μ$m line.

preprint2021arXiv

Atomic and molecular gas from the epoch of reionization down to redshift 2

Cosmic gas makes up about 90% of baryonic matter in the Universe and H$_2$ is the closest molecule to star formation. In this work we study cold neutral gas and its H$_2$ component at different epochs, exploiting state-of-the-art hydrodynamic simulations that include time-dependent atomic and molecular non-equilibrium chemistry coupled to star formation, feedback effects, different UV backgrounds presented in the recent literature and a number of additional processes - such as gas self-shielding, H$_2$ dust grain catalysis, photoelectric and cosmic-ray heating - occurring during structure formation (ColdSIM). We find neutral-gas mass density parameters $ Ω_{\rm neutral} \simeq $10$^{-3}$ and increasing from lower to higher redshift, in agreement with available HI data. Resulting H$_2$ fractions can be as high as $\sim $50% at $z\sim $4-8, in line with the latest high-$z$ measurements. Albeit dependent on the adopted UV background, derived $ Ω_{\rm H_2} $ values agree with observations up to $z\sim$7 and both HI and H$_2$ trends are better reproduced by our non-equilibrium H$_2$-based star formation modelling. The predicted gas depletion timescales decrease towards lower $z$, with H$_2$ depletion times remaining below the Hubble time and comparable to the dynamical time at all considered redshifts. This implies that non-equilibrium molecular cooling is efficient at driving cold-gas collapse in a broad variety of environments and since the very early cosmic epochs. In appendix, we show detailed analyses of individual processes, as well as simple numerical parameterizations and fits to account for them. Our findings suggest that, in addition to HI, non-equilibrium H$_2$ observations are pivotal probes for assessing cold-gas abundances and the role of UV background radiation - Abridged

preprint2020arXiv

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

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

preprint2020arXiv

Early black-hole seeds in the first billion years

Supermassive black holes with billion solar masses are in place already within the first Gyr, however, their origin and growth in such a short lapse of time is extremely challenging to understand. Here, we discuss the formation paths of early black-hole seeds, showing the limits of light black-hole seeds from stellar origin and the expected characteristics of heavy/massive black-hole seeds originated by gas direct collapse in peculiar primordial conditions. To draw conclusions on the possible candidates and the role of the ambient medium, we use results from N-body hydrodynamic simulations including atomic and molecular non-equilibrium abundance calculations, cooling, star formation, feedback mechanisms, stellar evolution, metal spreading of several heavy elements from SNII, AGB and SNIa, and multifrequency radiative transfer over 150 frequencies coupled to chemistry and SED emission for popII-I and popIII stellar sources. Standard stellar-origin light black holes are unlikely to be reliable seeds of early supermassive black holes, because, under realistic assumptions, they cannot grow significantly in less than a billion years. Alternatively, massive black-hole seeds might originate from direct collapse of pristine gas in primordial quiescent mini-haloes that are exposed to stellar radiation from nearby star forming regions. The necessary conditions required to form these heavy seeds must be complemented with information on the complex features of local environments and the fine balance between chemistry evolution and radiative transfer.

preprint2016arXiv

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

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

preprint2015arXiv

Bulk Flows and End of the Dark Ages with the SKA

The early Universe is a precious probe of the birth of primordial objects, first star formation events and consequent production of photons and heavy elements. Higher-order corrections to the cosmological linear perturbation theory predicts the formation of coherent supersonic gaseous streaming motions at decoupling time. These bulk flows impact the gas cooling process and determine a cascade effect on the whole baryon evolution. By analytical estimates and N-body hydrodynamical chemistry numerical simulations including atomic and molecular evolution, gas cooling, star formation, feedback effects and metal spreading for individual species from different stellar populations according to the proper yields and lifetimes, we discuss the role of these primordial bulk flows at the end of the dark ages and their detectable impacts during the first Gyr in view of the upcoming SKA mission. Early bulk flows can inhibit molecular gas cooling capabilities, suppressing star formation, metal spreading and the abundance of small primordial galaxies in the infant Universe. This can determine a delay in the re-ionization process and in the heating of neutral hydrogen making the observable HI signal during cosmic evolution patchier and noisier. The planned SKA mission will represent a major advance over existing instruments, since it will be able to probe the effects on HI 21cm at z ~ 6-20 and on molecular line emissions from first collapsing sites at z ~ 20-40. Therefore, it will be optimal to address the effects of primordial streaming motions on early baryon evolution and to give constraints on structure formation in the first Gyr.

preprint2015arXiv

Origin of Cosmic Chemical Abundances

Cosmological N-body hydrodynamic computations following atomic and molecular chemistry (e$^-$, H, H$^+$, H$^-$, He, He$^+$, He$^{++}$, D, D$^+$, H$_2$, H$_2^+$, HD, HeH$^+$), gas cooling, star formation and production of heavy elements (C, N, O, Ne, Mg, Si, S, Ca, Fe, etc.) from stars covering a range of mass and metallicity are used to explore the origin of several chemical abundance patterns and to study both the metal and molecular content during simulated galaxy assembly. The resulting trends show a remarkable similarity to up-to-date observations of the most metal-poor damped Lyman-$α$ absorbers at redshift $z\gtrsim 2$. These exhibit a transient nature and represent collapsing gaseous structures captured while cooling is becoming effective in lowering the temperature below $\sim 10^4\,\rm K$, before they are disrupted by episodes of star formation or tidal effects. Our theoretical results agree with the available data for typical elemental ratios, such as [C/O], [Si/Fe], [O/Fe], [Si/O], [Fe/H], [O/H] at redshifts $z\sim 2-7$. Correlations between HI and H$_2$ abundances show temporal and local variations and large spreads as a result of the increasing cosmic star formation activity from $z\sim 6$ to $z\sim 3$. The scatter we find in the abundance ratios is compatible with the observational data and is explained by simultaneous enrichment by sources from different stellar phases or belonging to different stellar populations. Simulated synthetic spectra support the existence of metal-poor cold clumps with large optical depth at $z\sim 6$ that could be potential population~III sites at low or intermediate redshift. The expected dust content is in line with recent determinations.

preprint2014arXiv

Signatures of very massive stars: supercollapsars and their cosmological rate

We compute the rate of supercollapsars by using cosmological, N-body, hydro, chemistry simulations of structure formation, following detailed stellar evolution according to proper yields (for He, C, N, O, Si, S, Fe, Mg, Ca, Ne, etc.) and lifetimes for stars having different masses and metallicities, and for different stellar populations (population III and population II-I). We find that supercollapsars are usually associated to dense, collapsing gas with little metal pollution and with abundances dominated by oxygen. The resulting supercollapsar rate is about $10^{-2}\,\rm yr^{-1} sr^{-1}$ at redshift $z=0$, and their contribution to the total rate is $ < 0.1 $ per cent, which explains why they have never been detected so far. Expected rates at redshift $z\simeq 6$ are of the order of $\sim 10^{-3}\,\rm yr^{-1} sr^{-1}$ and decrease further at higher $z$. Because of the strong metal enrichment by massive, short-lived stars, only $\sim 1$ supercollapsar generation is possible in the same star forming region. Given their sensitivity to the high-mass end of the primordial stellar mass function, they are suitable candidates to probe pristine population III star formation and stellar evolution at low metallicities.

preprint2014arXiv

The First Billion Years of a Warm Dark Matter Universe

We present results of cosmological N-body hydrodynamic chemistry simulations of primordial structure growth and evolution in a scenario with warm dark matter (WDM) having a mass of 3 keV (thermal relic) and compare with a model consisting of standard cold dark matter (CDM). We focus on the high-redshift universe ($z>6$), where the structure formation process should better reflect the primordial (linear) differences in terms of matter power spectrum. We find that early epochs can be exceptional probes of the dark-matter nature. Non-linear WDM power spectra and mass functions are up to 2 dex lower than in CDM and show spreads of factor of a few persisting in the whole first Gyr. Runaway molecular cooling in WDM haloes results severely inhibited because of the damping of power at large $k$ modes and hence cosmic (population III and II-I) star formation rate (SFR) is usually suppressed with respect to CDM predictions. Luminous objects formed in a WDM background are very rare at $z>10$, due to the sparser and retarded evolution of early WDM mini-haloes during the dark ages and their lack can be fitted with a simple analytical formula depending only on magnitude and redshift. Future high-$z$ observations of faint galaxies have the potential to discriminate between CDM and WDM scenarios by means of cosmic stellar mass density (SMD) and specific SFR, as well. When compared to the effects of alternative cosmologies (e.g. non-Gaussian or dark-energy models) or of high-order corrections at large $z$ (e.g. primordial streaming motions or changes in the pristine IMF) the ones caused by WDM are definitely more dramatic.

preprint2013arXiv

Effects Of Circum-nuclear Disk Gas Evolution And The Spin Of Central Black Holes

Mass and spin are the only two parameters needed to completely characterize black holes in General Relativity. However, the interaction between black holes and their environment is where complexity lies, as the relevant physical processes occur over a large range of scales. That is particularly relevant in the case of super-massive black holes (SMBHs), hosted in galaxy centers, and surrounded by swirling gas and various generations of stars. These compete with the SMBH for gas consumption and affect both dynamics and thermodynamics of the gas itself. How the behavior of such fiery environment influence the angular momentum of the gas accreted onto SMBHs, and, hence, black-hole spins is uncertain. We explore the interaction between SMBHs and their environment via first 3D sub-parsec resolution simulations (ranging from 0.1 pc to 1 kpc scales) that study the evolution of the SMBH spin by including the effects of star formation, stellar feedback, radiative transfer, and metal pollution according to the proper stellar yields and lifetimes. This approach is crucial to investigate the impact of star formation processes and feedback effects on the angular momentum of the material that could accrete on the central hole. We find that star formation and feedback mechanisms can locally inject significant amounts of entropy in the surrounding medium, and impact on the inflow inclination angles and Eddington fractions. As a consequence, the resulting trends show upper-intermediate equilibrium values for the spin parameter, a, of about 0.6 - 0.9, corresponding to radiative efficiencies ε= 9% - 15%. These results suggest that star formation feedback taking place in the circum-nuclear disk during the in-fall cannot induce alone very strong chaotic trends in the gas flow, quite independently from the different numerical parameters.

preprint2013arXiv

Hydrodynamical chemistry simulations of the SZ effect and the impacts from primordial non-Gaussianities

The impacts of Compton scattering of hot cosmic gas with the cosmic microwave background radiation (Sunyaev-Zel'dovich effect, SZ) are consistently quantified in Gaussian and non-Gaussian scenarios, by means of 3D numerical, N-body, hydrodynamic simulations, including cooling, star formation, stellar evolution and metal pollution (He, C, O, Si, Fe, S, Mg, etc.) from different stellar phases, according to proper yields for individual metal species and mass-dependent stellar lifetimes. Light cones are built through the simulation outputs and samples of one hundred maps for the resulting temperature fluctuations are derived for both Gaussian and non-Gaussian primordial perturbations. From them, we estimate the possible changes due to early non-Gaussianities on: SZ maps, probability distribution functions, angular power spectra and corresponding bispectra. We find that the different growth of structures in the different cases induces significant spectral distortions only in models with large non-Gaussian parameters, $f_{\rm NL}$. In general, the overall trends are covered by the non-linear, baryonic evolution, whose feedback mechanisms tend to randomize the gas behaviour and homogenize its statistical features, quite independently from the background matter distribution. Deviations due to non-Gaussianity are almost undistinguishable for $f_{\rm NL}\lesssim 100$, remaining always at few-per-cent level, within the error bars of the Gaussian scenario. Rather extreme models with $f_{\rm NL}\sim1000$ present more substantial deviations from the Gaussian case, overcoming baryon contaminations and showing discrepancies up to a factor of a few in the spectral properties.

preprint2013arXiv

Simulating extremely metal-poor gas and DLA metal content at redshift z=7

We present the first theoretical study of metals in damped Ly-alpha (DLA) systems at redshift z=7. The features of cold, primordial gas are studied by means of N-body, hydro, chemistry simulations, including atomic and molecular non-equilibrium chemistry, cooling, star formation for population III and population II-I regimes, stellar evolution, cosmic metal spreading according to proper yields (for He, C, O, Si, Fe, Mg, S, etc.) and lifetimes, and feedback effects. Theoretical expectations are then compared to recently available constraints from DLA observations. We find that DLA galaxies at z=7 account for 10 per cent of the whole galaxy population and for most of the metal-poor galaxies at these epochs. About 7 per cent of these DLA galaxies contain purely pristine material and 34 per cent of them consist of very weakly polluted gas, being, therefore, suitable candidate s as population III sites. The remaining 59 per cent are enriched above ~10^{-4} Zsun. Additionally, DLA candidates appear to have: gas masses <~2x10^8 Msun; very low star formation rate, ~10^{-3} - 10^{-2} Msun/yr (significantly weaker than late-time counterparts); mean molecular fractions covering a fairly wide range, x_mol~10^{-3}- 10^{-6}; typical metallicities Z<~3x10^{-3}Zsun and HI column densities NHI>~3x10^{20} cm^{-2}. They present no or weak correlations between their gas mass and Z, NHI, or x_mol; a moderate correlation between x_mol and Z, linked to the ongoing molecular-driven star formation and metal pollution processes; a mild anti-correlation between NHI and x_mol, due to H depletion into molecules; and a chemical content that is subject to environmental dependencies.

preprint2013arXiv

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

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

preprint2013arXiv

Statistical properties of mass, star formation, chemical content and rotational patterns in early z > 9 structures

We study the baryonic, chemical and dynamical properties of a significantly large sample of early proto-galaxies in the first 500 Myr of the Universe (redshift z>9), obtained from high-resolution numerical, N-body, hydrodynamical, chemistry simulations including atomic and molecular networks, gas cooling, star formation, stellar evolution and metal spreading for population III and population II-I regimes according to proper stellar yields and lifetimes. We find that first star formation events take place in halos with dark-matter mass M_{DM} > 2\times 10^6 M_{sun}. Early star forming objects have: molecular fractions from x_{mol} < 10^{-4} in quiescent structures up to x_{mol} > 0.1 in active regions; star formation rates SFR ~10^{-8}-10^{-3}M_{sun}/yr; and metallicities in the range ~10^{-8}-10^{-2}Z_{sun}. Roughly ~10% of high-z haloes host population II-I star formation and dominate the cosmic SFR density. They usually are bursty objects with mean specific SFR around ~10 Gyr^{-1} at z~9 and increasing with redshift up to ~10^2 Gyr^{-1}. Stellar feedback effects alter the baryonic content of the haloes and locally affect their chemical and thermo-dynamical properties, as reflected by the broadening of various physical relations. The establishment of gaseous rotationally-supported cores is quite uncommon, weakly related to the large-scale dark-matter behaviour and evolving in an intermittent fashion. The colder, molecular-rich phase tends to maintain any established rotational motion longer with respect to the hotter, metal-rich component, which is very sensitive to environmental processes. While the fraction of haloes featuring a significant amount of co-rotating, molecular-rich gas increases with cosmic time (from a few percent at z~20 up to ~5-15% at z~9), the chaotic nature of metal-enriched material does not lead to particular trends.

preprint2012arXiv

Counts of high-redshift GRBs as probe of primordial non-Gaussianities

We propose to use high-redshift long $γ$-ray bursts (GRBs) as cosmological tools to constrain the amount of primordial non-Gaussianity in the density field. By using numerical, N-body, hydrodynamic, chemistry simulations of different cosmological volumes with various Gaussian and non-Gaussian models, we self-consistently relate the cosmic star formation rate density to the corresponding GRB rate. Assuming that GRBs are fair tracers of cosmic star formation, we find that positive local non-Gaussianities, described in terms of the non-linear parameter, \fnl, might boost significantly the GRB rate at high redshift, $z \gg 6$. Deviations with respect to the Gaussian case account for a few orders of magnitude if \fnl$\sim 1000$, one order of magnitude for \fnl$\sim 100$, and a factor of $\sim 2$ for \fnl$\sim 50$. These differences are found only at large redshift, while at later times the rates tend to converge. Furthermore, a comparison between our predictions and the observed GRB data at $z > 6$ allows to exclude large negative \fnl, consistently with previous works. Future detections of any long GRB at extremely high redshift ($z\sim 15-20$) could favor non-Gaussian scenarios with positive \fnl. More stringent constraints require much larger high-$z$ GRB complete samples, currently not available in literature. By distinguishing the contributions to the GRB rate from the metal-poor population III regime, and the metal-enriched population II-I regime, we conclude that the latter is a more solid tracer of the underlying matter distribution, while the former is strongly dominated by feedback mechanisms from the first, massive, short-lived stars, rather than by possible non-Gaussian fluctuations. This holds quite independently of the assumed population III initial mass function.

preprint2012arXiv

Dark Matter Halo Environment for Primordial Star Formation

We study the statistical properties (such as shape and spin) of high-z halos likely hosting the first (PopIII) stars with cosmological simulations including detailed gas physics. In the redshift range considered ($11 < z < 16$) the average sphericity is $<s> = 0.3 \pm 0.1$, and for more than 90% of halos the triaxiality parameter is $T \lesssim 0.4$, showing a clear preference for oblateness over prolateness. Larger halos in the simulation tend to be both more spherical and prolate: we find $s \propto M_h^{α_s}$ and $T \propto M_h^{α_T}$, with $α_s \approx 0.128$ and $α_T= 0.276$ at z = 11. The spin distributions of dark matter and gas are considerably different at $z=16$, with the baryons rotating slower than the dark matter. At lower redshift, instead, the spin distributions of dark matter and gas track each other almost perfectly, as a consequence of a longer time interval available for momentum redistribution between the two components. The spin of both the gas and dark matter follows a lognormal distribution, with a mean value at z=16 of $<λ> =0.0184$, virtually independent of halo mass. This is in good agreement with previous studies. Using the results of two feedback models (MT1 and MT2) by McKee & Tan (2008) and mapping our halo spin distribution into a PopIII IMF, we find that at high-$z$ the IMF closely tracks the spin lognormal distribution. Depending on the feedback model, though, the distribution can be centered at $\approx 65 M_\odot$ (MT1) or $\approx 140 M_\odot$ (MT2). At later times, model MT1 evolves into a bimodal distribution with a second prominent peak located at $35-40 M_\odot$ as a result of the non-linear relation between rotation and halo mass. We conclude that the dark matter halo properties might be a key factor shaping the IMF of the first stars.

preprint2012arXiv

Effect of Intergalactic Medium on the Observability of Lyman Alpha Emitters during Cosmic Reionization

We perform a systematic study of how the inhomogeneities in the Inter-Galactic Medium (IGM) affect the observability of Lyman Alpha Emitters (LAEs) around the Epoch of Reionization. We focus on the IGM close to the galaxies as the detailed ionization distribution and velocity fields of this region could significantly influence the scattering of Ly-alpha photons off neutral H atoms as they traverse the IGM after escaping from the galaxy. We simulate the surface brightness (SB) maps and spectra of more than 100 LAEs at z=7.7 as seen by an observer at z=0. To achieve this, we extract the source properties of galaxies and their surrounding IGM from cosmological simulations of box sizes 5-30 Mpc/h and follow the coupled radiative transfer of ionizing and Ly-alpha radiation through the IGM using CRASH-alpha. We find that the simulated SB profiles are extended and their detailed structure is affected by inhomogeneities in the IGM, especially at high neutral fractions. The detectability of LAEs and the fraction of the flux observed depend heavily on the shape of the SB profile and the SB threshold (SB_th) of the observational campaign. Only ultradeep observations (e.g. SB_th ~ 10^-23 ergs/s/cm^2/arcsec^2) would be able to obtain the true underlying mass-luminosity relation and luminosity functions of LAEs. The details of our results depend on whether Ly-alpha photons are significantly shifted in the galaxy to longer wavelengths, the mean ionization fraction in the IGM and the clustering of ionizing sources. These effects can lead to an easier escape of Ly-alpha photons with less scattering in the IGM and a concentrated SB profile similar to the one of a point source. Finally, we show that the SB profiles are steeper at high ionization fraction for the same LAE sample which can potentially be observed from the stacked profile of a large number of LAEs.

preprint2012arXiv

Non-Gaussianities, early Universe, and GRBs. Tracing the primeval state of the Universe with number counts of Gamma-Ray Bursts

We investigate the effects of primordial non-Gaussianities in the primordial Universe on the baryonic structure formation process. By relating the cosmic star formation rate in Gaussian and non-Gaussian scenarios to the detectability of high-redshift sources of reionization, we derive the expected Gamma-Ray Burst rate in the different models. We find that counts of high-redshift Gamma-Ray Bursts can be used as cosmological probes of non-Gaussianities and that they are suitable candidates to distinguish non-Gaussian effects at early epochs.

preprint2012arXiv

Radiative feedback and cosmic molecular gas: numerical method

We present results from self-consistent 3D numerical simulations of cosmic structure formation with a multi-frequency radiative transfer scheme and non-equilibrium molecular chemistry of 13 primordial species (e-, H, H+, H-, He, He+, He++, H2, H2+, D, D+, HD, HeH+), performed by using the simulation code GADGET. We describe our implementation and show tests for ionized sphere expansion in a static and dynamic density field around a central radiative source, and for cosmological abundance evolution coupled with the cosmic microwave background radiation. As a demonstrative application of radiative feedback on molecular gas, we run also cosmological simulations of early structure formation in a ~1Mpc size box. Our tests agree well with analytical and numerical expectations. Consistently with other works, we find that ionization fronts from central sources can boost H2 fractions in shock-compressed gas. The tight dependence on H2 lead to a corresponding boost of HD fractions, as well. We see a strong lowering of the the typical molecular abundances up to several orders of magnitudes which partially hinders further gas collapse of pristine neutral gas, and clearly suggests the need of re-ionized gas or metal cooling for the formation of the following generation of structures.

preprint2011arXiv

Baryon history and cosmic star formation in non-Gaussian cosmological models: numerical simulations

We present the first numerical, N-body, hydrodynamical, chemical simulations of cosmic structure formation in the framework of non-Gaussian models. We study the impact of primordial non-Gaussianities on early chemistry (e, H, H+, H-, He, He+, He++, H2, H2+, D, D+, HD, HeH+), molecular and atomic gas cooling, star formation, metal (C, O, Si, Fe, Mg, S) enrichment, population III (popIII) and population II-I (popII) transition, and on the evolution of "visible" objects. We find that non-Gaussianities can have some consequences on baryonic structure formation at very early epochs, but the subsequent evolution at later times washes out any difference among the various models. When assuming reasonable values for primordial non-Gaussian perturbations, it turns out that they are responsible for: (i) altering early molecular fractions in the cold, dense gas phase of ~10 per cent; (ii) inducing small temperature fluctuations of <~10 per cent during the cosmic evolution of primordial objects; (iii) influencing the onset of the first star formation events, at z>~15, and of the popIII/popII transition of up to some 10^7yr; (iv) determining variations of <~10 per cent in the gas cloud and stellar mass distributions after the formation of the first structures; (v) causing only mild variations in the chemical history of the Universe. We stress, though, that purely non-Gaussian effects might be difficult to address, since they are strictly twisted with additional physical phenomena (e.g. primordial gas bulk flows, unknown primordial popIII stellar mass function, etc.) that have similar or stronger impact on the behaviour of the baryons.

preprint2011arXiv

Gas distribution, metal enrichment, and baryon fraction in Gaussian and non-Gaussian universes

We study the cosmological evolution of baryons in universes with and without primordial non-Gaussianities via (large scale) N-body/hydrodynamical simulations, including gas cooling, star formation, stellar evolution, chemical enrichment from both population III and population II regimes, and feedback effects. We find that large fnl values for non-Gaussianities can alter the gas probability distribution functions, the metal pollution history, the halo baryon, gas and stellar fractions, mostly at early times. More precisely: (i) non-Gaussianities lead to an earlier evolution of primordial gas, structures, and star formation; (ii) metal enrichment starts earlier (with respect to the Gaussian scenario) in non-Gaussian models with larger fnl; (iii) gas fractions within the haloes are not significantly affected by the different values of fnl, with deviations of ~1-10%; (iv) the stellar fraction is quite sensitive to non-Gaussianities at early times, with discrepancies reaching up to a factor of ~10 at very high z, and rapidly converging at low z; (v) the trends at low redshift are independent from fnl: they are mostly led by the ongoing baryonic evolution and by the feedback mechanisms, which determine a ~25%-30% discrepancy in the baryon fraction of galaxy groups/clusters with respect to the cosmic values; (vi) non-Gaussianity impacts on the cluster X-ray emission or on the SZ effect(s) are expected to be not very large and dominated by feedback mechanisms, whereas some effects on the 21-cm emission can be expected at early times; (vii) in order to address non-Gaussianities in the cosmological structure contest, high-redshift (z~10) investigations are required: first stars, galaxies, quasars, and GRBs may be potential cosmological probes of non-Gaussianities.

preprint2011arXiv

The imprint of cosmological non-Gaussianities on primordial structure formation

We study via numerical N-body/SPH chemistry simulations the effects of primordial non-Gaussianities on the formation of the first stars and galaxies, and investigate the impact of supernova feedback in cosmologies with different fnl. Density distributions are biased to higher values, so star formation and the consequent feedback processes take place earlier in high-fnl models and later in low-fnl ones. Mechanical feedback is responsible for shocking and evacuating the gas from star forming sites earlier in the highly non-Gaussian cases, because of the larger bias at high densities. Chemical feedback translates into high-redshift metal filling factors that are larger by some orders of magnitude for larger fnl, but that converge within one Gyr, for both population III and population II-I stellar regimes. The efficient enrichment process, though, leads to metallicities > 0.01 Zsun by redshift ~9, almost independently from fnl. The impact of non-Gaussianities on the formation of dark-matter haloes at high redshift is directly reflected in the properties of the gas in these haloes, as models with larger fnl show more concentrated gas profiles at early times. Non-Gaussian signatures in the gas behaviour are lost after the first feedback takes place and introduces a significant degree of turbulence and chaotic motions.

preprint2011arXiv

The interplay between chemical and mechanical feedback from the first generation of stars

We study cosmological simulations of early structure formation, including non-equilibrium molecular chemistry, metal pollution from stellar evolution, transition from population III (popIII) to population II (popII) star formation, regulated by a given critical metallicity, and feedback effects. We investigate the properties of early metal spreading from the different stellar populations and its interplay with primordial molecular gas. We find that, independently of the details about popIII modeling, after the onset of star formation, regions enriched below the critical level are mostly found in isolated environments, while popII star formation regions are much more clumped. Typical star forming haloes show average SN driven outflow rates of up to 10^{-4} Msun/yr in enriched gas, initially leaving the original star formation regions almost devoid of metals. The polluted material, which is gravitationally incorporated in over-dense environments on timescales of 10^7 yr, is mostly coming from external, nearby star forming sites ("gravitational enrichment"). In parallel, the pristine-gas inflow rates are between 10^{-3} - 10^{-1} Msun/yr. However, thermal feedback from SN generates turbulence and destroys molecules within the pristine gas, and only the polluted material, incorporated via gravitational enrichment, can continue to cool by atomic metal fine-structure transitions on time scales short enough to end the initial popIII regime within less than 10^8 yr.

preprint2010arXiv

The impact of primordial supersonic flows on early structure formation, reionization and the lowest-mass dwarf galaxies

Tseliakhovich & Hirata recently discovered that higher-order corrections to the cosmological linear-perturbation theory lead to supersonic coherent baryonic flows just after recombination (i.e. z~1020), with rms velocities of ~30 km/s relative to the underlying dark-matter distribution, on comoving scales of ~<3 Mpc/h. To study the impact of these coherent flows we performed high-resolution N-body plus SPH simulations in boxes of 5.0 and 0.7 Mpc/h, for bulk-flow velocities of 0 (as reference), 30 and 60 km/s. The simulations follow the evolution of cosmic structures by taking into account detailed, primordial, non-equilibrium gas chemistry (i.e. H, He, H2, HD, HeH, etc.), cooling, star formation, and feedback effects from stellar evolution. We find that these bulk flows suppress star formation in low-mass haloes (i.e. Mvir<~10^8Msun until z~13, lower the abundance of the first objects by ~1%-20%, and, as consequence, delay cosmic star formation history by ~2e7yr. The gas fractions in individual objects can change up to a factor of two at very early times. Coherent bulk flow, therefore, has implications for (i) the star-formation in the lowest-mass haloes (e.g. dSphs), (ii) the start of reionization by suppressing it in some patches of the Universe, and (iii) the heating (i.e. spin temperature) of neutral hydrogen. We speculate that the patchy nature of reionization and heating on several Mpc scales could lead to enhanced differences in the HI spin-temperature, giving rise to stronger variations in the HI brightness temperatures during the late dark ages.

preprint2009arXiv

Formation of Isolated Dwarf Galaxies with Feedback

We present results of high resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal-dependent cooling, star formation, feedback from type II and type Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky Way type galaxies. We find that a combination of feedback and the cosmic UV background results in the formation of galaxies with properties similar to the Local Group dwarf spheroidals, and that their effect is strongly moderated by the depth of the gravitational potential. Taking this into account, our models naturally reproduce the observed luminosities and metallicities. The final objects have halo masses between 2.3x10^8 and 1.1x10^9 solar masses, mean velocity dispersions between 6.5 and 9.7 kms-1, stellar masses ranging from 5x10^5 to 1.2x10^7 solar masses, median metallicities between [Fe/H] = -1.8 and -1.1, and half-light radii of the order of 200 to 300 pc, all comparable with Local Group dwarf spheroidals. Our simulations also indicate that the dwarf spheroidal galaxies observed today lie near a halo mass threshold around 10^9 solar masses, in agreement with stellar kinematic data, where supernova feedback not only suffices to completely expel the interstellar medium and leave the residual gas-free, but where the combination of feedback, UV radiation and self-shielding establishes a dichotomy of age distributions similar to that observed in the Milky Way and M31 satellites.

preprint2009arXiv

The onset of star formation in primordial haloes

Star formation remains an unsolved problem in astrophysics. Numerical studies of large-scale structure simulations cannot resolve the whole process and their approach usually assumes that only gas denser than a typical threshold can host and form stars. We investigate the onset of cosmological star formation and compare several very-high-resolution, three-dimensional, N-body/SPH simulations that include non-equilibrium, atomic and molecular chemistry, star formation prescriptions, and feedback effects. We study how primordial star formation depends on gas density thresholds, cosmological parameters and initial set-ups. For mean-density initial conditions, we find that standard low-density star-formation threshold (0.2 h^2/cm3) models predict the onset of star formation at z~25-31, depending on the adopted cosmology. In these models stars are formed regardless of the time between the moment when the threshold is reached and the effective runaway collapse. At high redshift, this time interval represents a significant fraction of the Hubble time and thus this assumption can induce large artificial off-sets to the onset of star formation. Choosing higher density thresholds (135 h^2/cm3) allows the entire cooling process to be followed, and the onset of star formation is then estimated to be at redshift z~12-16. When isolated, rare, high-density peaks are considered, the chemical evolution is much faster and the first star formation episodes occur at z > 40, almost regardless of the choice for the density threshold. These results could have implications for the formation redshift of the first cosmological objects, as inferred from direct numerical simulations of mean-density environments, and on the studies of the reionization history of the universe.