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.

The recent discovery of high redshift dusty galaxies implies a rapid dust enrichment of their interstellar medium (ISM). To interpret these observations, we run a cosmological simulation in a 30$h^{-1}$ cMpc/size volume down to $z \approx 4$. We use the hydrodynamical code dustyGadget, which accounts for the production of dust by stellar populations and its evolution in the ISM. We find that the cosmic dust density parameter ($\Omega_{\rm d}$) is mainly driven by stellar dust at $z \gtrsim 10$, so that mass- and metallicity-dependent yields are required to assess the dust content in the first galaxies. At $z \lesssim 9$ the growth of grains in the ISM of evolved systems (Log$(M_{\star}/M_{\odot})>8.5$) significantly increases their dust mass, in agreement with observations in the redshift range $4 \lesssim z < 8$. Our simulation shows that the variety of high redshift galaxies observed with ALMA can naturally be accounted for by modeling the grain-growth timescale as a function of the physical conditions in the gas cold phase. In addition, the trends of dust-to-metal (DTM) and dust-to-gas (${\cal D}$) ratios are compatible with the available data. A qualitative investigation of the inhomogeneous dust distribution in a representative massive halo at $z \approx 4$ shows that dust is found from the central galaxy up to the closest satellites along polluted filaments with $\rm Log({\cal D}) \leq -2.4$, but sharply declines at distances $d \gtrsim 30$ kpc along many lines of sight, where $\rm Log({\cal D}) \lesssim -4.0$.

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.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

We present cosmological hydrodynamical simulations including atomic and molecular non-equilibrium chemistry, multi-frequency radiative transfer (0.7-100 eV sampled over 150 frequency bins) and stellar population evolution to investigate the host candidates of the seeds of supermassive black holes coming from direct collapse of gas in primordial haloes (direct-collapse black holes, DCBHs). We consistently address the role played by atomic and molecular cooling, stellar radiation and metal spreading of C, N, O, Ne, Mg, Si, S, Ca, Fe, etc. from primordial sources, as well as their implications for nearby quiescent proto-galaxies under different assumptions for early source emissivity, initial mass function and metal yields. We find that putative DCBH host candidates need powerful primordial stellar generations, since common solar-like stars and hot OB-type stars are neither able to determine the conditions for direct collapse nor capable of building up a dissociating Lyman-Werner background radiation field. Thermal and molecular features of the identified DCBH host candidates in the scenario with very massive primordial stars seem favourable, with illuminating Lyman-Werner intensities featuring values of 1-50 J21. Nevertheless, additional non-linear processes, such as merger events, substructure formation, rotational motions and photo-evaporation, should inhibit pure DCBH formation in 2/3 of the cases. Local turbulence may delay gas direct collapse almost irrespectively from other environmental conditions. The impact of large Lyman-Werner fluxes at distances smaller than 5 kpc is severely limited by metal pollution.

Here we investigate how LyC-opaque systems present in the intergalactic medium at $z\approx3$ can distort the spectral shape of a uniform UV background (UVB) through radiative transfer (RT) effects. With this aim in mind, we perform a multi-frequency RT simulation through a cosmic volume of $10h^{-1}$~cMpc scale polluted by metals, and self-consistently derive the ions of all the species. The UVB spatial fluctuations are traced by the ratio of He$\, \rm \scriptstyle II\ $ and H$\, \rm \scriptstyle I\ $ column density, $\eta$, and the ratio of C$\,{\rm {\scriptstyle IV\ }}$ and Si$\,{\rm {\scriptstyle IV\ }}$ optical depths, $\zeta$. We find that: (i) $\eta$ spatially fluctuates through over-dense systems ($\Delta$) with statistically significant deviations $\delta\eta >25$\% in 18\% of the volume ; (ii) same fluctuations in $\zeta$ are also present in $34$\% of the enriched domain (only 8\% of the total volume) and derive from a combination of RT induced effects and in-homogeneous metal enrichment, both effective in systems with $\Delta > 1.5$.

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

We use numerical N-body hydrodynamical simulations with varying PopIII stellar models to investigate the possibility of detecting first star signatures with observations of high-redshift damped Ly$\alpha$ absorbers (DLAs). The simulations include atomic and molecular cooling, star formation, energy feedback and metal spreading due to the evolution of stars with a range of masses and metallicities. Different initial mass functions (IMFs) and corresponding metal-dependent yields and lifetimes are adopted to model primordial stellar populations. The DLAs in the simulations are selected according to either the local gas temperature (temperature selected) or the host mass (mass selected). We find that 3\% (40\%) of mass (temperature) selected high-$z$ ($z\ge5.5$) DLAs retain signatures of pollution from PopIII stars, independently from the first star model. Such DLAs have low halo mass ($<10^{9.6}\,\rm M_{\odot}$), metallicity ($<10^{-3}\,\rm Z_{\odot}$) and star formation rate ($<10^{-1.5}\,\rm M_{\odot}\,yr^{-1}$). { Metal abundance ratios of DLAs imprinted in the spectra of QSO} can be useful tools to infer the properties of the polluting stellar generation and to constrain the first star mass ranges. Comparing the abundance ratios derived from our simulations to those observed in DLAs at $z\ge5$, we find that most of these DLAs are consistent within errors with PopII stars dominated enrichment and strongly disfavor the pollution pattern of very massive first stars (i.e. 100~$\rm M_{\odot}$-500~$\rm M_{\odot}$). However, some of them could still result from the pollution of first stars in the mass range [0.1, 100]~$\rm M_{\odot}$. In particular, we find that the abundance ratios from SDSS J1202+3235 are consistent with those expected from PopIII enrichment dominated by massive (but not extreme) first stars.

By using N-body hydrodynamical cosmological simulations in which the chemistry of major metals and molecules is consistently solved for, we study the interaction of metallic fine-structure lines with the CMB. Our analysis shows that the collisional induced emissions in the OI 145 $\mu$m and CII 158 $\mu$m lines during reionization introduce a distortion of the CMB spectrum at low frequencies ($\nu < 300$ GHz) with amplitudes up to $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-8}$-$10^{-7}$, i.e., at the $\sim 0.1$ percent level of FIRAS upper limits. Shorter wavelength fine-structure transitions (OI 63 $\mu$m, FeII 26 $\mu$m, and SiII 35 $\mu$m) typically sample the reionization epoch at higher observing frequencies ($\nu > 400$ GHz). This corresponds to the Wien tail of the CMB spectrum and the distortion level induced by those lines may be as high as $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-4}$. The angular anisotropy produced by these lines should be more relevant at higher frequencies: while practically negligible at $\nu=145 $GHz, signatures from CII 158 $\mu$m and OI 145 $\mu$m should amount to 1%-5% of the anisotropy power measured at $l \sim 5000$ and $\nu=220 $GHz by the ACT and SPT collaborations (after assuming $\Delta \nu_{\rm obs}/\nu_{\rm obs}\simeq 0.005$ for the line observations). Our simulations show that anisotropy maps from different lines (e.g., OI 145 $\mu$m and CII 158 $\mu$m) at the same redshift show a very high degree ($>0.8$) of spatial correlation, allowing for the use of observations at different frequencies to unveil the same snapshot of the reionization epoch. Finally, our simulations demonstrate that line-emission anisotropies extracted in narrow frequency/redshift shells are practically uncorrelated in frequency space, thus enabling standard methods for removal of foregrounds that vary smoothly in frequency, just as in HI 21 cm studies.

We study the possibility to detect and distinguish signatures of enrichment from PopIII stars in observations of PopII GRBs (GRBIIs) at high redshift by using numerical N-body/hydrodynamical simulations including atomic and molecular cooling, star formation and metal spreading from stellar populations with different initial mass functions (IMFs), yields and lifetimes. PopIII and PopII star formation regimes are followed simultaneously and both a top-heavy and a Salpeter-like IMF for pristine PopIII star formation are adopted. We find that the fraction of GRBIIs hosted in a medium previously enriched by PopIII stars (PopIII-dominated) is model independent. Typical abundance ratios, such as [Si/O] vs [C/O] and [Fe/C] vs [Si/C], can help to disentangle enrichment from massive and intermediate PopIII stars, while low-mass first stars are degenerate with regular PopII generations. The properties of galaxies hosting PopIII-dominated GRBIIs are not very sensitive to the particular assumption on the mass of the first stars.

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

The basic processes of the formation of the first stars in the primordial Universe are outlined and the implications for cosmological structure formation discussed. By employing theoretical and numerical models of cosmic structure evolution embedded within N-body hydrodynamical chemistry simulations, predictions for the production of the first heavy elements in the Universe are given. These results are then compared against measured data of UV luminosities and metal abundances in different kinds of observations in order to draw conclusions on the chemical and thermal state of the cosmic medium at different cosmological epochs.

We present results from multifrequency radiative hydrodynamical chemistry simulations addressing primordial star formation and related stellar feedback from various populations of stars, stellar energy distributions (SEDs) and initial mass functions. Spectra for massive stars, intermediate-mass stars and regular solar-like stars are adopted over a grid of 150 frequency bins and consistently coupled with hydrodynamics, heavy-element pollution and non-equilibrium species calculations. Powerful massive population III stars are found to be able to largely ionize H and, subsequently, He and He$^+$, causing an inversion of the equation of state and a boost of the Jeans masses in the early intergalactic medium. Radiative effects on star formation rates are between a factor of a few and 1 dex, depending on the SED. Radiative processes are responsible for gas heating and photoevaporation, although emission from soft SEDs has minor impacts. These findings have implications for cosmic gas preheating, primordial direct-collapse black holes, the build-up of "cosmic fossils" such as low-mass dwarf galaxies, the role of AGNi during reionization, the early formation of extended disks and angular-momentum catastrophe.

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-$\alpha$ 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.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with Mstar > 10^9 Msun allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population

Revealing hidden patterns in astronomical data is often the path to fundamental scientific breakthroughs; meanwhile the complexity of scientific inquiry increases as more subtle relationships are sought. Contemporary data analysis problems often elude the capabilities of classical statistical techniques, suggesting the use of cutting edge statistical methods. In this light, astronomers have overlooked a whole family of statistical techniques for exploratory data analysis and robust regression, the so-called Generalized Linear Models (GLMs). In this paper -- the first in a series aimed at illustrating the power of these methods in astronomical applications -- we elucidate the potential of a particular class of GLMs for handling binary/binomial data, the so-called logit and probit regression techniques, from both a maximum likelihood and a Bayesian perspective. As a case in point, we present the use of these GLMs to explore the conditions of star formation activity and metal enrichment in primordial minihaloes from cosmological hydro-simulations including detailed chemistry, gas physics, and stellar feedback. We predict that for a dark mini-halo with metallicity $\approx 1.3 \times 10^{-4} Z_{\bigodot}$, an increase of $1.2 \times 10^{-2}$ in the gas molecular fraction, increases the probability of star formation occurrence by a factor of 75%. Finally, we highlight the use of receiver operating characteristic curves as a diagnostic for binary classifiers, and ultimately we use these to demonstrate the competitive predictive performance of GLMs against the popular technique of artificial neural networks.

We investigate signatures of population III (PopIII) stars in the metal-enriched environment of GRBs originating from population II-I (PopII/I) stars by using abundance ratios derived from numerical simulations that follow stellar evolution and chemical enrichment. We find that at $z>10$ more than $10%$ of PopII/I GRBs explode in a medium previously enriched by PopIII stars (we refer to them as GRBII$\rightarrow$III). Although the formation of GRBII$\rightarrow$III is more frequent than that of pristine PopIII GRBs (GRBIIIs), we find that the expected GRBII$\rightarrow$III observed rate is comparable to that of GRBIIIs, due to the usually larger luminosities of these latter. GRBII$\rightarrow$III events take place preferentially in small proto-galaxies with stellar masses $\rm M_\star \sim 10^{4.5} - 10^7\,\rm M_\odot$, star formation rates $\rm SFR \sim 10^{-3}-10^{-1}\,\rm M_\odot/yr$ and metallicities $Z \sim 10^{-4}-10^{-2}\,\rm Z_\odot$. On the other hand, galaxies with $Z < 10^{-2.8}\,\rm Z_\odot$ are dominated by metal enrichment from PopIII stars and should preferentially host GRBII$\rightarrow$III. Hence, measured GRB metal content below this limit could represent a strong evidence of enrichment by pristine stellar populations. We discuss how to discriminate PopIII metal enrichment on the basis of various abundance ratios observable in the spectra of GRBs' afterglows. By employing such analysis, we conclude that the currently known candidates at redshift $z\simeq 6$ -- i.e. GRB 050904 \cite[][]{2006Natur.440..184K} and GRB 130606A \cite[][]{2013arXiv1312.5631C} -- are likely not originated in environments pre-enriched by PopIII stars.

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.

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.

We present a novel approach, based on robust principal components analysis (RPCA) and maximal information coefficient (MIC), to study the redshift dependence of halo baryonic properties. Our data are composed of a set of different physical quantities for primordial minihaloes: dark-matter mass ($M_{\mathrm{dm}}$), gas mass ($M_{\mathrm{gas}}$), stellar mass ($M_{\mathrm{star}}$), molecular fraction (${\mathrm{x_{mol}}}$), metallicity ($Z$), star formation rate (SFR) and temperature. We find that $M_{\mathrm{dm}}$ and $M_{\mathrm{gas}}$ are dominant factors for variance, particularly at high redshift. Nonetheless, with the emergence of the first stars and subsequent feedback mechanisms, ${\mathrm{x_{mol}}}$, SFR and $Z$ start to have a more dominant role. Standard PCA gives three principal components (PCs) capable to explain more than 97 per cent of the data variance at any redshift (two PCs usually accounting for no less than 92 per cent), whilst the first PC from the RPCA analysis explains no less than 84 per cent of the total variance in the entire redshift range (with two PCs explaining $ \gtrsim 95$ per cent anytime). Our analysis also suggests that all the gaseous properties have a stronger correlation with $M_{\mathrm{gas}}$ than with $M_{\mathrm{dm}}$, while $M_{\mathrm{gas}}$ has a deeper correlation with ${\mathrm{x_{mol}}}$ than with $Z$ or SFR. This indicates the crucial role of gas molecular content to initiate star formation and consequent metal pollution from Population III and Population II/I regimes in primordial galaxies. Finally, a comparison between MIC and Spearman correlation coefficient shows that the former is a more reliable indicator when halo properties are weakly correlated.

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.

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.

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.

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.

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.

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 \epsilon = 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.

We investigate the nature of high-z host galaxies of long Gamma-Ray Bursts (LGRBs) by means of state-of-the-art numerical simulations of cosmic structure formation and evolution of galaxies. We combine results from different runs with various box sizes and resolutions. By assigning to each simulated galaxy the probability to host a LGRB, assumed to be proportional to the mass of young stars, we provide a full description of the physical properties of high-z LGRB host galaxy population. We find that LGRBs at z>6 are hosted in galaxies with typical star formation rates SFR \sim 0.03-0.3 Msun yr^{-1}, stellar masses M \sim 10^{6-8} Msun, and metallicities Z \sim 0.01-0.1 Zsun. Furthermore, the ratio between their doubling time and the corresponding cosmic time seems to be universally equal to ~0.1-0.3, independently from the redshift. The distribution of their UV luminosity places LGRB hosts in the faint-end of the galaxy luminosity function, well below the current capabilities of space- or ground-based optical facilities. This is in line with recent reports of non-detection of LGRB hosts using extremely deep HST and VLT observations. In conclusion, high-z LGRBs are found to trace the position of those faint galaxies that are thought to be the major actors in the re-ionization of the Universe.

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^{\alpha_s}$ and $T \propto M_h^{\alpha_T}$, with $\alpha_s \approx 0.128$ and $\alpha_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 $<\lambda> =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.

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.

We propose to use high-redshift long $\gamma$-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.

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.

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.

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.

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.

Because massive, low-metallicity population III (PopIII) stars may produce very powerful long gamma-ray bursts (LGRBs), high-redshift GRB observations could probe the properties of the first stars. We analyze the correlation between early PopIII stars and LGRBs by using cosmological N-body/hydrodynamical simulations, which include detailed chemical evolution, cooling, star formation, feedback effects and the transition between PopIII and more standard population I/II (PopII/I) stars. From the Swift observed rate of LGRBs, we estimate the fraction of black holes that will produce a GRB from PopII/I stars to be in the range 0.028<f_{GRB}<0.140, depending on the assumed upper metallicity of the progenitor. Assuming that as of today no GRB event has been associated to a PopIII star, we estimate the upper limit for the fraction of LGRBs produced by PopIII stars to be in the range 0.006<f_{GRB}<0.022. When we apply a detection threshold compatible with the BAT instrument, we find that the expected fraction of PopIII GRBs (GRB3) is ~10% of the full LGRB population at z>6, becoming as high has 40% at z>10. Finally, we study the properties of the galaxies hosting our sample of GRB3. We find that the average metallicity of the galaxies hosting a GRB3 is typically higher than the critical metallicity used to select the PopIII stars, due to the efficiency in polluting the gas above such low values. We also find that the highest probability of finding a GRB3 is within galaxies with a stellar mass <10^7 Msun, independently from the redshift.

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.

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.

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.

We present results from the first cosmological simulations which study the onset of primordial, metal-free (population III), cosmic star formation and the transition to the present-day, metal-rich star formation (population II-I), including molecular (H$_2$, HD, etc.) evolution, tracing the injection of metals by supernov{\ae} into the surrounding intergalactic medium and following the change in the initial stellar mass function (IMF) according to the metallicity of the corresponding stellar population. Our investigation addresses the role of a wide variety of parameters (critical metallicity for the transition, IMF slope and range, SN/pair-instability SN metal yields, star formation threshold, resolution, etc.) on the metal-enrichment history and the associated transition in the star formation mode. All simulations present common trends. Metal enrichment is very patchy, with rare, unpolluted regions surviving at all redshifts, inducing the simultaneous presence of metal-free and metal-rich star formation regimes. As a result of the rapid pollution within high-density regions due to the first SN/pair-instability SN, local metallicity is quickly boosted above the critical metallicity for the transition. The population III regime lasts for a very short period during the first stages of star formation ($\sim 10^7\,\rm yr$), and its average contribution to the total star formation rate density drops rapidly below $\sim 10^{-3}-10^{-2}$.

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.

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.

Cooling is the main process leading to the condensation of gas in the dark matter potential wells and consequently to star and structure formation. In a metal-free environment, the main available coolants are H, He, H$_2$ and HD; once the gas is enriched with metals, these also become important in defining the cooling properties of the gas. We discuss the implementation in Gadget-2 of molecular and metal cooling at temperatures lower that $\rm10^4 K$, following the time dependent properties of the gas and pollution from stellar evolution. We have checked the validity of our scheme comparing the results of some test runs with previous calculations of cosmic abundance evolution and structure formation, finding excellent agreement. We have also investigated the relevance of molecule and metal cooling in some specific cases, finding that inclusion of HD cooling results in a higher clumping factor of the gas at high redshifts, while metal cooling at low temperatures can have a significant impact on the formation and evolution of cold objects.

This submission has been withdrawn by arXiv administrators because it is a duplicate of 0704.2182.

We present the first hydrodynamic N-body simulations of primordial gas clouds responsible for the reionisation process in dark energy cosmologies. We compare the cosmological constant scenario with a SUGRA quintessence model with marked dynamics in order to highlight effects due to the different acceleration histories imposed by the dark energy. We show that both the number density of gas clouds and their clumpiness keep a record of the expansion rate during evolution, similar to the non-linear dark matter profile at virialisation, as was recently demonstrated by Dolag et al. (2004). Varying the shape of the primordial power spectrum, we show how this effect is mitigated by a running spectral index decreasing the power at small scales. Our results demonstrate that, in order to constrain the dark energy from large scale structures, one must track its effects down to the distribution of luminous matter.