Far-infrared (FIR) emission lines are a powerful tool to investigate the properties of the interstellar medium, especially in high-redshift galaxies, where ALMA observations have provided unprecedented information. Interpreting such data with state-o
f-the-art cosmological simulations post-processed with CLOUDY, has provided insights on the internal structure and gas dynamics of these systems. However, no detailed investigation of the consistency and uncertainties of this kind of analysis has been performed to date. Here, we compare different approaches to estimate FIR line emission from state-of-the-art cosmological simulations, either with CLOUDY or with on-the-fly non-equilibrium chemistry. We find that [CII]$_{158mu}$ predictions are robust to the model variations we explored. [OI] emission lines, that typically trace colder and denser gas relative to [CII]$_{158mu}$, are instead model-dependent, as these lines are strongly affected by the thermodynamic state of the gas and non-equilibrium photoionisation effects. For the same reasons, [OI] lines represent an excellent tool to constrain emission models, hence future observations targeting these lines will be crucial.
Recent stacked ALMA observations have revealed that normal, star-forming galaxies at $zapprox 6$ are surrounded by extended ($approx 10,mathrm{kpc}$) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a
model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading factor, $eta$, and the parent galaxy dark matter halo circular velocity, $v_c$. The outflow model successfully matches the observed [CII] surface brightness profile if $eta = 3.20 pm 0.10$ and $v_c = 170 pm 10{,rm km,s^{-1}}$, corresponding to a dynamical mass of $approx 10^{11}, mathrm{M}_odot$. The predicted outflow rate and velocity range are $128 pm 5 ,mathrm{M}_odot {rm yr}^{-1}$ and $300-500 {,rm km,s^{-1}}$, respectively. We conclude that: (a) extended halos can be produced by cooling outflows; (b) the large $eta$ value is marginally consistent with starburst-driven outflows, but it might indicate additional energy input from AGN; (c) the presence of [CII] halos requires an ionizing photon escape fraction from galaxies $f_{rm esc} ll 1$. The model can be readily applied also to individual high-$z$ galaxies, as those observed, e.g., by the ALMA ALPINE survey now becoming available.
Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of Walking the Line, a 3-day conference held in 2
018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.
We investigate quasar outflows at $z geq 6$ by performing zoom-in cosmological hydrodynamical simulations. By employing the SPH code GADGET-3, we zoom in the $2 R_{200}$ region around a $2 times 10^{12} M_{odot}$ halo at $z = 6$, inside a $(500 ~ {rm
Mpc})^3$ comoving volume. We compare the results of our AGN runs with a control simulation in which only stellar/SN feedback is considered. Seeding $10^5 M_{odot}$ BHs at the centers of $10^{9} M_{odot}$ halos, we find the following results. BHs accrete gas at the Eddington rate over $z = 9 - 6$. At $z = 6$, our most-massive BH has grown to $M_{rm BH} = 4 times 10^9 M_{odot}$. Fast ($v_{r} > 1000$ km/s), powerful ($dot{M}_{rm out} sim 2000 M_{odot}$/yr) outflows of shock-heated low-density gas form at $z sim 7$, and propagate up to hundreds kpc. Star-formation is quenched over $z = 8 - 6$, and the total SFR (SFR surface density near the galaxy center) is reduced by a factor of $5$ ($1000$). We analyse the relative contribution of multiple physical process: (i) disrupting cosmic filamentary cold gas inflows, (ii) reducing central gas density, (iii) ejecting gas outside the galaxy; and find that AGN feedback has the following effects at $z = 6$. The inflowing gas mass fraction is reduced by $sim 12 %$, the high-density gas fraction is lowered by $sim 13 %$, and $sim 20 %$ of the gas outflows at a speed larger than the escape velocity ($500$ km/s). We conclude that quasar-host galaxies at $z geq 6$ are accreting non-negligible amount of cosmic gas, nevertheless AGN feedback quenches their star formation dominantly by powerful outflows ejecting gas out of the host galaxy halo.
