Search

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

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

A tight relation between the [CII]158$\mu$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($\Sigma_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [CII]-deficiency we have developed an analytical model that fits local [CII] data, and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of $\Sigma_{\rm CII}$ with the surface star formation rate ($\Sigma_*$). However, for $\Sigma_* > 1 M_\odot~{\rm yr}^{-1}~{\rm kpc}^{-2}$, $\Sigma_{\rm CII}$ saturates. We conclude that underluminous [CII] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation ($\kappa_s \gg 1$), parameterized by the "burstiness" parameter $\kappa_s$; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [CII] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [OIII]88$\mu$m, CIII]1909A, CO lines). Simple formulae are given to interpret available data for low and high-$z$ galaxies., Comment: Accepted for publication in MNRAS. 12 pages, 6 Figures. Comments welcome

We study the photoevaporation of Jeans-unstable molecular clumps by isotropic FUV (6 eV $< {\rm h}\nu$ < 13.6 eV) radiation, through 3D radiative transfer hydrodynamical simulations implementing a non-equilibrium chemical network that includes the formation and dissociation of H$_2$. We run a set of simulations considering different clump masses ($M=10-200$ M$_\odot$) and impinging fluxes ($G_0=2\times 10^3-8\times 10^4$ in Habing units). In the initial phase, the radiation sweeps the clump as an R-type dissociation front, reducing the H$_2$ mass by a factor 40-90%. Then, a weak ($\mathcal{M}\simeq 2$) shock develops and travels towards the centre of the clump, which collapses while loosing mass from its surface. All considered clumps remain gravitationally unstable even if radiation rips off most of the clump mass, showing that external FUV radiation is not able to stop clump collapse. However, the FUV intensity regulates the final H$_2$ mass available for star formation: for example, for $G_0 < 10^4$ more than 10% of the initial clump mass survives. Finally, for massive clumps ($\sim 100$ M$_\odot$) the H$_2$ mass increases by 25-50% during the collapse, mostly because of the rapid density growth that implies a more efficient H$_2$ self-shielding., Comment: 17 pages, 16 figures; added references, corrected typos

We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([CII], [NII], and [OIII]). At $z=8$, the most massive galaxy, `Freesia', has an age $t_\star \simeq 409\,{\rm Myr}$, stellar mass $M_{\star} \simeq 4.2\times 10^9 {\rm M}_{\odot}$, and a star formation rate ${\rm SFR} \simeq 11.5\,{\rm M}_{\odot}{\rm yr}^{-1}$, due to a recent burst. Freesia has two stellar components (A and B) separated by $\simeq 2.5\, {\rm kpc}$; other 11 galaxies are found within $56.9 \pm 21.6 \, {\rm kpc}$. The mean ISRF in the Habing band is $G = 7.9\, G_0$ and is spatially uniform; in contrast, the ionisation parameter is $U = 2^{+20}_{-2} \times 10^{-3}$, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionising escape fraction from Freesia is $f_{\rm esc}\simeq 2\%$. While [CII] emission is extended (radius 1.54 kpc), [OIII] is concentrated in Freesia-A (0.85 kpc), where the ratio $\Sigma_{\rm [OIII]}/\Sigma_{\rm [CII]} \simeq 10$. As many high-$z$ galaxies, Freesia lies below the local [CII]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) which disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [CII]-SFR deviations., Comment: 22 pages, 14 figures, accepted by MNRAS