Search

The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$\alpha$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$\alpha$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected through the [OIII]($\lambda4969,5008$) doublet at $5.35M_{\rm UV}>-19.5$)., Comment: 15 pages, 13 figures. Accepted for publication in A&A. Main figs. 1, 2, 10 (SFR surface density), 11 (spectra and positions of galaxies in the bubble). Corrected error in escape fraction estimate from SFR surface density

Satellite galaxies within the Milky Way's (MW) virial radius $R_{\mathrm{vir}}$ are typically devoid of cold gas due to ram pressure stripping by the MW's corona. The density of this corona is poorly constrained today and essentially unconstrained in the past, but can be estimated using ram pressure stripping. In this paper, we probe the MW corona at $z\approx 1.6$ using the Draco dwarf spheroidal galaxy. We assume that i) Draco's orbit is determined by its interaction with the MW, whose dark matter halo we evolve in time following cosmologically-motivated prescriptions, ii) Draco's star formation was quenched by ram pressure stripping and iii) the MW's corona is approximately smooth, spherical and in hydrostatic equilibrium. We used GAIA proper motions to set the initial conditions and Draco's star formation history to estimate its past gas content. We found indications that Draco was stripped of its gas during the first pericentric passage. Using 3D hydrodynamical simulations at a resolution that enables us to resolve individual supernovae and assuming no tidal stripping, which we estimate to be a minor effect, we find a density of the MW corona $\geq 8\times 10^{-4}$ cm$^{-3}$ at a radius $\approx 0.72R_{\mathrm{vir}}$. This provides evidence that the MW's corona was already in place at $z\approx 1.6$ and with a higher density than today. If isothermal, this corona would have contained all the baryons expected by the cosmological baryon fraction. Extrapolating to today shows good agreement with literature constraints if feedback has removed $\lesssim 30$% of baryons accreted onto the halo., Comment: 20 pages, 11 figures. Accepted for publication in MNRAS

Ly$\alpha$ nebulae ubiquitously found around z>2 quasars can supply unique constraints on the properties of the Circumgalactic Medium, such as its density distribution, provided the quasar halo mass is known. We present a new method to constrain quasar halo masses based on the line-of-sight velocity dispersion maps of Ly$\alpha$ nebulae. By using MUSE-like mock observations obtained from cosmological hydrodynamic simulations under the assumption of maximal quasar fluorescence, we show that the velocity dispersion radial profiles of Ly$\alpha$-emitting gas are strongly determined by gravity and that they are thus self-similar with respect to halo mass when rescaled by the virial radius. Through simple analytical arguments and by exploiting the kinematics of HeII1640\.A emission for a set of observed nebulae, we show that Ly$\alpha$ radiative transfer effects plausibly do not change the shape of the velocity dispersion profiles but only their normalisation without breaking their self-similarity. Taking advantage of these results, we define the variable $\eta^{140-200}_{40-100}$ as the ratio of the median velocity dispersion in two specifically selected annuli and derive an analytical relation between $\eta^{140-200}_{40-100}$ and the halo mass which can be directly applied to observations. We apply our method to 37 observed quasar Ly$\alpha$ nebulae at 3

Galactic halos accrete material from the intergalactic medium (IGM) and part of this accretion is expected to be in the form of cool ($T\sim10^4$ K) gas. A signature of this process could reside in the detection of numerous clouds in the circumgalactic medium (CGM). However, whether this material is able to accrete onto the galaxies and feed their star formation or, instead, evaporates into the CGM hot phase (corona, $T\sim10^6$ K), is not yet understood. Here, we investigate the evolution of cool CGM clouds accreted from the IGM and falling through the hot corona of low-redshift disc galaxies, using 3D high-resolution hydrodynamical simulations. We include the effects of gravity due to the dark matter halo, isotropic thermal conduction, radiative cooling and an ionizing UV background. We explored different values of parameters such as the halo mass, coronal mass, initial cloud velocity and strength of the thermal conduction. We find that the clouds lose the vast majority of their mass at distances larger than half of the galaxy virial radius and are completely dissolved in the corona before reaching the central galaxy. Resolving the Field length with at least 5-7 cells is crucial to correctly capture the evolution of the infalling cool gas. Our results indicate that cool IGM accretion can not feed star formation in $z\sim0$ star-forming galaxies in halos with masses of $10^{11.9}\ M_{\odot}$ or above. This suggests that present-day massive star-forming galaxies can sustain their star formation only via the spontaneous or induced cooling of their hot corona., Comment: 18 pages, 10 Figures. Accepted for publication in MNRAS

Observationally mapping the relation between galaxies and the intergalactic medium (IGM) is of key interest for studies of cosmic reionization. Diffuse hydrogen gas has typically been observed in HI Lyman-$\alpha$ (Ly$\alpha$) absorption in the spectra of bright background quasars. However, it is important to extend these measurements to background galaxies as quasars become increasingly rare at high redshift and rarely probe closely separated sight-lines. Here we use deep integral field spectroscopy in the MUSE eXtremely Deep Field to demonstrate the measurement of the Ly$\alpha$ transmission at $z\approx4$ in absorption to a background galaxy at $z=4.77$. The HI transmission is consistent with independent quasar sight-lines at similar redshifts. Exploiting the high number of spectroscopic redshifts of faint galaxies (500 between $z=4.0-4.7$ within a radius of 8 arcmin) that are tracers of the density field, we show that Ly$\alpha$ transmission is inversely correlated with galaxy density, i.e. transparent regions in the Ly$\alpha$ forest mark under-dense regions at $z\approx4$. Due to large-scale clustering, galaxies are surrounded by excess HI absorption over the cosmic mean out to 4 cMpc/h. We also find that redshifts from the peak of the Ly$\alpha$ line are typically offset from the systemic redshift by +170 km/s. This work extends results from $z\approx 2 - 3$ to higher redshifts and demonstrates the power of deep integral field spectroscopy to simultaneously measure the ionization structure of the IGM and the large-scale density field in the early Universe., Comment: Accepted for publication in MNRAS. Main text 10 pages, 9 figures. Key results in Fig 4 (Lya forest transmission in the MXDF field) and Fig 9 (transmission - galaxy distance cross-correlation)

The discovery of giant quasar Ly$\alpha$ nebulae at $z>2$ has opened up the possibility to directly study in emission the Circumgalactic and Intergalactic Medium (CGM/IGM). However, the resonant nature of the Ly$\alpha$ line and its different emission mechanisms hamper the ability to constrain both the kinematics and physical properties of the CGM/IGM. Here, we present results of a pilot project aiming at the detection of CGM H$\alpha$ emission, a line which does not suffer from these limitations. To this end, we first used KCWI to detect Ly$\alpha$ emission around three bright quasars with $2.2510^{12}$M$_{\odot}$)., Comment: 15 pages, 12 figures, 2 tables

We use a dynamical model of galactic fountain to study the neutral extraplanar gas (EPG) in the nearby spiral galaxy NGC 2403. We have modelled the EPG as a combination of material ejected from the disc by stellar feedback (i.e. galactic fountain) and gas accreting from the inner circumgalactic medium (CGM). This accretion is expected to occur because of cooling/condensation of the hot CGM (corona) triggered by the fountain. Our dynamical model reproduces the distribution and kinematics of the EPG H$\mathrm{\scriptsize{I}}$ emission in NGC 2403 remarkably well and suggests a total EPG mass of $4.7^{+1.2}_{-0.9}\times10^8\mathrm{M}_\odot$, with a typical scale height of around 1 kpc and a vertical gradient of the rotation velocity of $-10.0\pm2.7\,\mathrm{km\,s^{-1}\,kpc^{-1}}$. The best-fitting model requires a characteristic outflow velocity of $50\pm10\,\mathrm{km\,s^{-1}}$. The outflowing gas starts out mostly ionised and only becomes neutral later in the trajectory. The accretion rate from the condensation of the inner hot CGM inferred by the model is 0.8$\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$, approximately equal to the star formation rate in this galaxy (0.6$\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$). We show that the accretion profile, which peaks at a radius of about 4.5$\,$kpc, predicts a disc growth rate compatible with the observed value. Our results indicate that fountain-driven corona condensation is a likely mechanism to sustain star formation as well as the disc inside-out growth in local disc galaxies., Comment: 15 pages, 9 figures, accepted for publication in MNRAS

As the closest $L^{\ast}$ galaxy to our own Milky Way, the Andromeda galaxy (M31) is an ideal laboratory for studies of galaxy evolution. The AMIGA project has recently provided observations of the cool ($T\sim10^4$ K) phase of the circumgalactic medium (CGM) of M31, using HST/COS absorption spectra along $\sim40$ background QSO sightlines, located up to and beyond the galaxy virial radius. Based on these data, and by the means of semi-analytic models and Bayesian inference, we provide here a physical description of the origin and dynamics of the cool CGM of M31. We investigate two competing scenarios, in which (i) the cool gas is mostly produced by supernova(SN)-driven galactic outflows or (ii) it mostly originates from infall of gas from the intergalactic medium. In both cases, we take into account the effect of gravity and hydrodynamical interactions with a hot corona, which has a cosmologically motivated angular momentum. We compare the outputs of our models to the observed covering factor, silicon column density and velocity distribution of the AMIGA absorbers. We find that, to explain the observations, the outflow scenario requires an unphysically large (> 100\%) efficiency for SN feedback. Our infall models, on the other hand, can consistently account for the AMIGA observations and the predicted accretion rate, angular momentum and metallicity are consistent with a cosmological infall from the intergalactic medium., Comment: Accepted for publication in MNRAS

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$\alpha$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$\alpha$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Ly\alpha}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$\alpha$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$\alpha$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $\xi_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background. We conclude proposing tests to further develop our suggested Ly$\alpha$-anchored formalism., Comment: Accepted for publication in MNRAS. Figure 2 shows the main result -- the comoving emissivity due to bright LAEs. Our fiducial model is based on results presented in our companion paper -- Naidu & Matthee et al. 2022, arXiv: 2110.11961

The ionizing photon escape fraction (LyC $f_{\rm{esc}}$) of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshift. Here we circumvent this uncertainty by inferring $f_{\rm{esc}}$ with resolved (R>4000) LyA profiles from the X-SHOOTER LyA survey at z=2 (XLS-z2). We select leakers ($f_{\rm{esc}}>20$%) and non-leakers ($f_{\rm{esc}}<5$%) from a representative sample of $>0.2 L^{*}$ LyA emitters (LAEs). With median stacked spectra of these subsets covering 1000-8000 {\AA} (rest-frame) we investigate the conditions for LyC $f_{\rm{esc}}$. We find the following differences between leakers vs. non-leakers: (i) strong nebular CIV and HeII emission vs. non-detections, (ii) O32~8.5 vs. ~3, (iii) Ha/Hb indicating no dust vs. E(B-V)~0.3, (iv) MgII emission close to the systemic velocity vs. redshifted, optically thick MgII, (v) LyA $f_{\rm{esc}}$ of ~50% vs. ~10%. The extreme EWs in leakers (O3+Hb~1100 {\AA}) constrain the characteristic timescale of LyC escape to ~3-10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers -- extremely ionizing stellar populations, low column densities, a dust-free, high ionization state ISM -- occur simultaneously in the $f_{\rm{esc}}>20\%$ stack, suggesting they are causally connected, and motivating why indicators like O32 may suffice to constrain $f_{\rm{esc}}$ at z>6 with JWST. The leakers comprise half our sample, have a median LyC $f_{\rm{esc}}$~50%, and an ionising production efficiency $\log({\xi_{\rm{ion}}/\rm{Hz\ erg^{-1}}})$~25.9. These results show LAEs -- the type of galaxies rare at z=2, but that become the norm at higher redshift -- are highly efficient ionizers, with extreme $\xi_{\rm{ion}}$ and prolific $f_{\rm{esc}}$ occurring in sync. (ABRIDGED), Comment: Submitted to MNRAS. Figure 4 and Table 2 show key results from the stacks. Figure 8 presents a unifying scheme for the LyC duty cycle. Our companion paper (Matthee & Naidu et al. 2021) works out implications of these results for the cosmic UV background at $z\approx2-8$. Comments warmly welcomed and greatly appreciated!

Two currently debated problems in galaxy evolution, the fundamentally local or global nature of the main sequence of star formation and the evolution of the mass-size relation of star forming galaxies (SFGs), are shown to be intimately related to each other. As a preliminary step, a growth function $g$ is defined, which quantifies the differential change in half-mass radius per unit increase in stellar mass ($g = d \log R_{1/2}/d \log M_\star$) due to star formation. A general derivation shows that $g = K \Delta(sSFR)/sSFR$, meaning that $g$ is proportional to the relative difference in specific star formation rate between the outer and inner half of a galaxy, with $K$ a dimensionless structural factor for which handy expressions are provided. As an application, it is shown that galaxies obeying a fundamentally local main sequence also obey, to a good approximation, $g \simeq \gamma n$, where $\gamma$ is the slope of the normalized local main sequence ($sSFR \propto \Sigma_\star^{-\gamma}$) and $n$ the Sersic index. An exact expression is also provided. Quantitatively, a fundamentally local main sequence is consistent with SFGs growing along a stationary mass-size relation, but inconsistent with the continuation at $z=0$ of evolutionary laws derived at higher $z$. This demonstrates that either the main sequence is not fundamentally local, or the mass-size relation of SFGs has converged to an equilibrium state some finite time in the past, or both., Comment: Accepted for publication in MNRAS

The relations between the specific angular momenta ($j$) and masses ($M$) of galaxies are often used as a benchmark in analytic models and hydrodynamical simulations as they are considered to be amongst the most fundamental scaling relations. Using accurate measurements of the stellar ($j_\ast$), gas ($j_{\rm gas}$), and baryonic ($j_{\rm bar}$) specific angular momenta for a large sample of disc galaxies, we report the discovery of tight correlations between $j$, $M$, and the cold gas fraction of the interstellar medium ($f_{\rm gas}$). At fixed $f_{\rm gas}$, galaxies follow parallel power laws in 2D $(j,M)$ spaces, with gas-rich galaxies having a larger $j_\ast$ and $j_{\rm bar}$ (but a lower $j_{\rm gas}$) than gas-poor ones. The slopes of the relations have a value around 0.7. These new relations are amongst the tightest known scaling laws for galaxies. In particular, the baryonic relation ($j_{\rm bar}-M_{\rm bar}-f_{\rm gas}$), arguably the most fundamental of the three, is followed not only by typical discs but also by galaxies with extreme properties, such as size and gas content, and by galaxies previously claimed to be outliers of the standard 2D $j-M$ relations. The stellar relation ($j_{\ast}-M_{\ast}-f_{\rm gas}$) may be connected to the known $j_\ast-M_\ast-$bulge fraction relation; however, we argue that the $j_{\rm bar}-M_{\rm bar}-f_{\rm gas}$ relation can originate from the radial variation in the star formation efficiency in galaxies, although it is not explained by current disc instability models., Comment: A&A Letters, in press. Data catalogue will be available via CDS and at this link https://unishare.nl/index.php/s/NMQRYfrrDpj8iaj

The intensity of the Cosmic UV background (UVB), coming from all sources of ionising photons such as star-forming galaxies and quasars, determines the thermal evolution and ionization state of the intergalactic medium (IGM) and is, therefore, a critical ingredient for models of cosmic structure formation. Most of the previous estimates are based on the comparison between observed and simulated Lyman-$\alpha$ forest. We present the results of an independent method to constrain the product of the UVB photoionisation rate and the covering fraction of Lyman limit systems (LLSs) by searching for the fluorescent Lyman-$\alpha$ emission produced by self-shielded clouds. Because the expected surface brightness is well below current sensitivity limits for direct imaging, we developed a new method based on three-dimensional stacking of the IGM around Lyman-$\alpha$ emitting galaxies (LAEs) between 2.93 that are consistent with previous measurements, with a preference for relatively low UVB intensities at z=3, and which suggest a non-monotonic decrease of $\Gamma$HI with increasing redshift between 3

We present the first results from the X-SHOOTER Lyman-$\alpha$ survey at $z=2$ (XLS-$z2$). XLS-$z2$ is a deep spectroscopic survey of 35 Lyman-$\alpha$ emitters (LAEs) utilising $\approx90$ hours of exposure time with VLT/X-SHOOTER and covers rest-frame Ly$\alpha$ to H$\alpha$ emission with R$\approx4000$. We present the sample selection, the observations and the data reduction. Systemic redshifts are measured from rest-frame optical lines for 33/35 sources. In the stacked spectrum, our LAEs are characterised by an interstellar medium with little dust, a low metallicity and a high ionisation state. The ionising sources are young hot stars that power strong emission-lines in the optical and high ionisation lines in the UV. The LAEs exhibit clumpy UV morphologies and have outflowing kinematics with blue-shifted SiII absorption, a broad [OIII] component and a red-skewed Ly$\alpha$ line. Typically 30 % of the Ly$\alpha$ photons escape, of which one quarter on the blue side of the systemic velocity. A fraction of Ly$\alpha$ photons escapes directly at the systemic suggesting clear channels enabling a $\approx10$ % escape of ionising photons, consistent with an inference based on MgII. A combination of a low effective HI column density, a low dust content and young star-burst determine whether a star forming galaxy is observed as a LAE. The first is possibly related to outflows and/or a fortunate viewing angle, while we find that the latter two in LAEs are typical for their stellar mass of 10$^9$ M$_{\odot}$., Comment: Main text 26 pages. Accepted for publication in MNRAS

The characterization of the large amount of gas residing in the galaxy halos, the so called circumgalactic medium (CGM), is crucial to understand galaxy evolution across cosmic time. We focus here on the the cool ($T\sim10^4$ K) phase of this medium around star-forming galaxies in the local universe, whose properties and dynamics are poorly understood. We developed semi-analytical parametric models to describe the cool CGM as an outflow of gas clouds from the central galaxy, as a result of supernova explosions in the disc (galactic wind). The cloud motion is driven by the galaxy gravitational pull and by the interactions with the hot ($T\sim10^6$ K) coronal gas. Through a bayesian analysis, we compare the predictions of our models with the data of the COS-Halos and COS-GASS surveys, which provide accurate kinematic information of the cool CGM around more than 40 low-redshift star-forming galaxies, probing distances up to the galaxy virial radii. Our findings clearly show that a supernova-driven outflow model is not suitable to describe the dynamics of the cool circumgalactic gas. Indeed, to reproduce the data, we need extreme scenarios, with initial outflow velocities and mass loading factors that would lead to unphysically high energy coupling from the supernovae to the gas and with supernova efficiencies largely exceeding unity. This strongly suggests that, since the outflows cannot reproduce most of the cool gas absorbers, the latter are likely the result of cosmological inflow in the outer galaxy halos, in analogy to what we have previously found for early-type galaxies., Comment: Accepted for publication on MNRAS

We present the results from a MUSE survey of twelve $z\simeq3.15$ quasars, which were selected to be much fainter (20

Galaxy scaling laws, such as the Tully-Fisher, mass-size and Fall relations, can provide extremely useful clues on our understanding of galaxy formation in a cosmological context. Some of these relations are extremely tight and well described by one single parameter (mass), despite the theoretical existence of secondary parameters such as spin and concentration, which are believed to impact these relations. In fact, the residuals of these scaling laws appear to be almost uncorrelated with each other, posing significant constraints on models where secondary parameters play an important role. Here, we show that a possible solution is that such secondary parameters are correlated amongst themselves, in a way that removes correlations in observable space. In particular, we focus on how the existence of an anti-correlation between the dark matter halo spin and its concentration -- which is still debated in simulations -- can weaken the correlation of the residuals of the Tully-Fisher and mass-size relations. Interestingly, using simple analytic galaxy formation models, we find that this happens only for a relatively small portion of the parameter space that we explored, which suggests that this idea could be used to derive constraints to galaxy formation models that are still unexplored., Comment: Accepted for publication in A&A

In the last decades, much effort has been put into finding the star formation law which could unequivocally link the gas and the star formation rate (SFR) densities measured on sub-kiloparsec scale in star-forming galaxies. The conventional approach of using the observed surface densities to infer star formation laws has however revealed a major and well-known issue, as such relations are valid for the high-density regions of galaxies but break down in low-density and HI-dominated environments. Recently, an empirical correlation between the total gas (HI+H$_2$) and the star formation rate (SFR) volume densities was obtained for a sample of nearby disc galaxies and for the Milky Way. This volumetric star formation (VSF) law is a single power-law with no break and a smaller intrinsic scatter with respect to the star formation laws based on the surface density. In this work, we explore the VSF law in the regime of dwarf galaxies in order to test its validity in HI-dominated, low-density, and low-metallicity environments. In addition, we assess this relation in the outskirts of spiral galaxies, which are low-density and HI-dominated regions similar to dwarf galaxies. Remarkably, we find the VSF law, namely $\rho_\mathrm{SFR} \propto \rho_\mathrm{gas}^\alpha$ with $\alpha \approx 2$, is valid for both these regimes. This result indicates that the VSF law, which holds unbroken for a wide range of gas ($\approx 3$ dex) and SFR ($\approx 6$ dex) volume densities, is the empirical relation with the smallest intrinsic scatter and is likely more fundamental than surface-based star formation laws., Comment: 17 pages, 7 figures, accepted by A&A

CR7 is among the most luminous Lyman-$\alpha$ emitters (LAEs) known at $z = 6.6$ and consists of at least three UV components that are surrounded by Lyman-$\alpha$ (Ly$\alpha$) emission. Previous studies have suggested that it may host an extreme ionising source. Here, we present deep integral field spectroscopy of CR7 with VLT/MUSE. We measure extended emission with a similar halo scale length as typical LAEs at $z\approx5$. CR7's Ly$\alpha$ halo is clearly elongated along the direction connecting the multiple components, likely tracing the underlying gas distribution. The Ly$\alpha$ emission originates almost exclusively from the brightest UV component, but we also identify a faint kinematically distinct Ly$\alpha$ emitting region nearby a fainter component. Combined with new near-infrared data, the MUSE data show that the rest-frame Ly$\alpha$ equivalent width (EW) is $\approx100$ {\AA}. This is a factor four higher than the EW measured in low-redshift analogues with carefully matched Ly$\alpha$ profiles (and thus arguably HI column density), but this EW can plausibly be explained by star formation. Alternative scenarios requiring AGN powering are also disfavoured by the narrower and steeper Ly$\alpha$ spectrum and much smaller IR to UV ratio compared to obscured AGN in other Ly$\alpha$ blobs. CR7's Ly$\alpha$ emission, while extremely luminous, resembles the emission in more common LAEs at lower redshifts very well and is likely powered by a young metal poor starburst., Comment: Main text 16 pages, 9 figures. Accepted for publication by MNRAS

It is well known that gas in galaxy discs is highly turbulent, but there is much debate on which mechanism can energetically maintain this turbulence. Among the possible candidates, supernova (SN) explosions are likely the primary drivers but doubts remain on whether they can be sufficient in regions of moderate star formation activity, in particular in the outer parts of discs. In this paper, we measure the SN efficiency $\eta$, namely the fraction of the total SN energy needed to sustain turbulence in galaxies, and verify that SNe can indeed be the sole driving mechanism. The key novelty of our approach is that we take into account the increased turbulence dissipation timescale associated to the flaring in outer regions of gaseous discs. We analyse the distribution and kinematics of HI and CO in 10 nearby star-forming galaxies to obtain the radial profiles of the kinetic energy per unit area, for both the atomic gas and the molecular gas. We use a theoretical model to reproduce the observed energy with the sum of turbulent energy from SNe, as inferred from the observed star formation rate (SFR) surface density, and the gas thermal energy. We find that the observed kinetic energy is remarkably well reproduced by our model across the whole extent of the galactic discs, assuming $\eta$ constant with the galactocentric radius. Taking into account the uncertainties on the SFR surface density and on the atomic gas phase, we obtain that the median SN efficiencies for our sample of galaxies are $\langle \eta_\mathrm{atom} \rangle=0.015_{-0.008}^{+0.018}$ for the atomic gas and $\langle \eta_\mathrm{mol} \rangle = 0.003_{-0.002}^{+0.006}$ for the molecular gas. We conclude that SNe alone can sustain gas turbulence in nearby galaxies with only few percent of their energy and that there is essentially no need for any further source of energy., Comment: 18 pages, 6 figures, and 2 appendices; accepted by Astronomy & Astrophysics; abridged abstract

We study the gas kinematics of a sample of six isolated gas-rich low surface brightness galaxies, of the class called ultra-diffuse galaxies (UDGs). These galaxies have recently been shown to be outliers from the baryonic Tully-Fisher relation (BTFR), as they rotate much slower than expected given their baryonic mass, and to have baryon fractions similar to the cosmological mean. By means of a 3D kinematic modelling fitting technique, we show that the HI in our UDGs is distributed in "thin" regularly rotating discs and we determine their rotation velocity and gas velocity dispersion. We revisit the BTFR adding galaxies from other studies. We find a previously unknown trend between the deviation from the BTFR and the disc scale length valid for dwarf galaxies with circular speeds < 45 km/s, with our UDGs being at the extreme end. Based on our findings, we suggest that the high baryon fractions of our UDGs may originate due to the fact that they have experienced weak stellar feedback, likely due to their low star formation rate surface densities, and as a result they did not eject significant amounts of gas out of their discs. At the same time, we find indications that our UDGs may have higher-than-average stellar specific angular momentum, which can explain their large optical scale lengths., Comment: Accepted for publication in MNRAS. v2: a few typos have been corrected and a couple of references added

Galactic outflows are thought to eject baryons back out to the circum-galactic medium (CGM). Studies based on metal absorption lines (MgII in particular) in the spectra of background quasars indicate that the gas is ejected anisotropically, with galactic winds likely leaving the host in a bi-conical flow perpendicular to the galaxy disk. In this paper, we present a detailed analysis of an outflow from a z = 0.7 "green-valley" galaxy (log($M_*$/$\mathrm{M}_\odot$) = 9.9; SFR = 0.5 $\mathrm{M}_\odot\,\mathrm{yr}^{-1}$) probed by two background sources part of the MUSE Gas Flow and Wind (MEGAFLOW) survey. Thanks to a fortuitous configuration with a background quasar (SDSSJ1358+1145) and a bright background galaxy at $z = 1.4$, both at impact parameters of $\approx 15\,\mathrm{kpc}$, we can - for the first time - probe both the receding and approaching components of a putative galactic outflow around a distant galaxy. We measure a significant velocity shift between the MgII absorption from the two sightlines ($84\pm17\,\mathrm{km}\,\mathrm{s}^{-1}$), which is consistent with the expectation from our simple fiducial wind model, possibly combined with an extended disk contribution., Comment: Submitted to MNRAS. After addressing referee's comments

The observed properties of the Lyman-$\alpha$ (Ly$\alpha$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$\alpha$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$\alpha$ equivalent width (EW$_0\approx38$ {\AA}). These data are combined with deep resolved [CII]$_{\rm 158 \mu m}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$\alpha$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$\alpha$ and [CII] are similarly extended beyond the UV, with effective radius r$_{\rm eff} = 2.1\pm0.2$ kpc for a single exponential model or r$_{\rm eff, Ly\alpha, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$\alpha$ profile is broader and redshifted with respect to the [CII] line (by 213 km s$^{-1}$), but there are spatial variations that are qualitatively similar in both lines and coincide with resolved UV components. This suggests that the emission originates from two components with plausibly different HI column densities. We place VR7 in the context of other galaxies at similar and lower redshift. The Ly$\alpha$ halo scale length is similar at different redshifts and velocity shifts with respect to the systemic are typically smaller. Overall, we find little indications of a more neutral vicinity at higher redshift. This means that the local ($\sim 10$ kpc) neutral gas conditions that determine the observed Ly$\alpha$ properties in VR7 resemble the conditions in post-re-ionisation galaxies., Comment: 13 pages, 10 figures. Accepted for publication in MNRAS

Cosmic filaments play a crucial role in galaxy evolution transporting gas from the intergalactic medium into galaxies. However, little is known about the efficiency of this process and whether the gas is accreted in a homogenous or clumpy way. Recent observations suggest the presence of broad gas density distributions in the circumgalactic medium which could be related to the accretion of filaments. By means of high-resolution hydrodynamical simulations, we explore here the evolution of cold accreting filaments flowing through the hot circumgalactic medium (CGM) of high-z galaxies. In particular, we examine the nonlinear effects of Kelvin-Helmholtz instability (KHI) on the development of broad gas density distributions and on the formation of cold, dense clumps. We explore a large parameter space in filament and perturbation properties, such as, filament Mach number, initial perturbation wavelength, and thickness of the interface between the filament and the halo. We find that the time averaged density distribution of the cold gas is qualitatively consistent with a skewed log-normal probability distribution function (PDF) plus an additional component in form of a high density tail for high Mach-numbers. Our results suggest a tight correlation between the accreting velocity and the maximum densities developing in the filament which is consistent with the variance-Mach number relation for turbulence. Therefore, cosmological accretion could be a viable mechanism to produce turbulence and broad gas density distributions within the CGM., Comment: 12 pages, 14 figures, submitted to MNRAS on April 3rd 2019

We investigate the baryon content of the circumgalactic medium (CGM) within the virial radius of $M_h \sim 10^{12} \; M_\odot$ haloes at z ~ 3, by modelling the surface brightness profile of the giant Ly$\alpha$ nebulae recently discovered by MUSE around bright QSOs at this redshift. We initially assume fluorescent emission from cold photo-ionized gas confined by the pressure of a hot halo. Acceptable CGM baryon fractions (equal or smaller than the cosmological value) require that the cold gas occupies $\lesssim$ 1% of the volume, but is about as massive as the hot gas. CGM baryon fractions as low as 30% of the cosmic value, as predicted by some strongly ejective feedback models at this redshift, are not easy to reconcile with observations, under our assumptions, unless both the QSO-hosting haloes at $z\sim3$ are more massive than recent BOSS estimates based on clustering and the photo-ionized gas is colder than expected in a standard QSO ionizing radiation field. We also consider the option that the emission is dominated by photons scattered from the QSO broad line region. In this scenario, a very stringent lower limit to the baryon fraction can be obtained under the extreme assumption of optically thin scattering. We infer in this case a baryon fraction of at least 70% of the cosmic value, for fiducial parameters. Lower values require halo masses or gas temperatures different than expected, or that some mechanism keeps the cold gas systematically over-pressured with respect to the ambient medium., Comment: Accepted for publication in MNRAS. 22 pages, 10 Figures

The circumgalactic medium (CGM) of galaxies consists of a multiphase gas with components at very different temperatures, from $10^ {4}$ K to $10^ {7}$ K. One of the greatest puzzle about this medium is the presence of a large amount of low-temperature ($T\sim10^4$ K) gas around quiescent early-type galaxies (ETGs). Using semi-analytical parametric models, we describe the cool CGM around massive, low-redshift ETGs as the cosmological accretion of gas into their dark matter halos, resulting in an inflow of clouds from the external parts of the halos to the central galaxies. We compare our predictions with the observations of the COS-LRG collaboration. We find that inflow models can successfully reproduce the observed kinematics, the number of absorbers and the column densities of the cool gas. Our MCMC fit returns masses of the cool clouds of about $10^5\ \rm{M}_{\odot}$ and shows that they must evaporate during their journey due to hydrodynamic interactions with the hot gas. We conclude that the cool gas present in the halos of ETGs likely cannot reach the central regions and feed the galaxy star formation, thus explaining why these passive objects are no longer forming stars., Comment: 17 pages, 11 figures and 3 tables, accepted for publication in A&A

We use the MusE GAs FLOw and Wind (MEGAFLOW) survey to study the kinematics of extended disk-like structures of cold gas around $z\approx1$ star-forming galaxies. The combination of VLT/MUSE and VLT/UVES observations allows us to connect the kinematics of the gas measured through MgII quasar absorption spectroscopy to the kinematics and orientation of the associated galaxies constrained through integral field spectroscopy. Confirming previous results, we find that the galaxy-absorber pairs of the MEGAFLOW survey follow a strong bimodal distribution, consistent with a picture of MgII absorption being predominantly present in outflow cones and extended disk-like structures. This allows us to select a bona-fide sample of galaxy-absorber pairs probing these disks for impact parameters of 10-70 kpc. We test the hypothesis that the disk-like gas is co-rotating with the galaxy disks, and find that for 7 out of 9 pairs the absorption velocity shares the sign of the disk velocity, disfavouring random orbits. We further show that the data are roughly consistent with inflow velocities and angular momenta predicted by simulations, and that the corresponding mass accretion rates are sufficient to balance the star formation rates., Comment: 22 pages + Supplementary Appendix (15 pages); submitted to MNRAS; version after addressing referee's comments

To understand star formation in galaxies, we investigate the star formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) profiles for galaxies, based on a well-defined sample of 976 star-forming MaNGA galaxies. We find that the typical $\Sigma_{\rm SFR}$ profiles within 1.5Re of normal SF galaxies can be well described by an exponential function for different stellar mass intervals, while the sSFR profile shows positive gradients, especially for more massive SF galaxies. This is due to the more pronounced central cores or bulges rather than the onset of a `quenching' process. While galaxies that lie significantly above (or below) the star formation main sequence (SFMS) show overall an elevation (or suppression) of $\Sigma_{\rm SFR}$ at all radii, this central elevation (or suppression) is more pronounced in more massive galaxies. The degree of central enhancement and suppression is quite symmetric, suggesting that both the elevation and suppression of star formation are following the same physical processes. Furthermore, we find that the dispersion in $\Sigma_{\rm SFR}$ within and across the population is found to be tightly correlated with the inferred gas depletion time, whether based on the stellar surface mass density or the orbital dynamical time. This suggests that we are seeing the response of a simple gas-regulator system to variations in the accretion rate. This is explored using a heuristic model that can quantitatively explain the dependence of $\sigma(\Sigma_{\rm SFR})$ on gas depletion timescale. Variations in accretion rate are progressively more damped out in regions of low star-formation efficiency leading to a reduced amplitude of variations in star-formation., Comment: 24 pages, 16 figures, accepted for publication in ApJ

With a projected size of about 450 kpc at z~2.3, the Slug Ly-alpha nebula is a rare laboratory to study, in emission, the properties of the intergalactic gas in the Cosmic Web. Since its discovery, the Slug has been the subject of several spectroscopic follow-ups to constrain the properties of the emitting gas. Here we report the results of a deep MUSE integral-field spectroscopic search for non-resonant, extended HeII1640 and metal emission. Extended HeII radiation is detected on scales of about 100 kpc, but only in some regions associated with the bright Ly-alpha emission and a continuum-detected source, implying large and abrupt variations in the line ratios across adjacent regions in projected space. The recent detection of associated H-alpha emission and similar abrupt variations in the Ly-alpha kinematics, strongly suggest that the HeII/Ly-alpha gradient is due to large variations in the physical distances between the associated quasar and these regions. This implies that the overall length of the emitting structure could extend to physical Mpc scales and be mostly oriented along our line of sight. At the same time, the relatively low HeII/Ly-alpha values suggest that the emitting gas has a broad density distribution that - if expressed in terms of a lognormal - implies dispersions as high as those expected in the interstellar medium of galaxies. These results strengthen the possibility that the density distribution of intergalactic gas at high-redshift is extremely clumpy and multiphase on scales below our current observational spatial resolution of a few physical kpc., Comment: accepted for publication in MNRAS

Throughout the Hubble time, gas makes its way from the intergalactic medium into galaxies fuelling their star formation and promoting their growth. One of the key properties of the accreting gas is its angular momentum, which has profound implications for the evolution of, in particular, disc galaxies. Here, we discuss how to infer the angular momentum of the accreting gas using observations of present-day galaxy discs. We first summarize evidence for ongoing inside-out growth of star forming discs. We then focus on the chemistry of the discs and show how the observed metallicity gradients can be explained if gas accretes onto a disc rotating with a velocity 20-30% lower than the local circular speed. We also show that these gradients are incompatible with accretion occurring at the edge of the discs and flowing radially inward. Finally, we investigate gas accretion from a hot corona with a cosmological angular momentum distribution and describe how simple models of rotating coronae guarantee the inside-out growth of disc galaxies., Comment: 5 pages, 2 figures, proceeding of FM6 "Galactic Angular Momentum", XXX IAU General Assembly, Vienna 2018

Star formation (SF) laws are fundamental relations between the gas content of a galaxy and its star formation rate (SFR) and play key roles in galaxy evolution models. In this paper, we present new empirical SF laws of disc galaxies based on volume densities. Following the assumption of hydrostatic equilibrium, we calculated the radial growth of the thickness of the gaseous discs in the combined gravitational potential of dark matter, stars, and gas for 12 nearby star-forming galaxies. This allowed us to convert the observed surface densities of gas and SFR into the deprojected volume densities. We found a tight correlation with slope in the range 1.3-1.9 between the volume densities of gas (HI+H$_2$) and the SFR with a significantly smaller scatter than the surface-based (Kennicutt) law and no change in the slope over five orders of magnitude. This indicates that taking into account the radial increase of the thickness of galaxy discs is crucial to reconstruct their three-dimensional density profiles, in particular in their outskirts. Moreover, our result suggests that the break in the slope seen in the Kennicutt law is due to disc flaring rather than to a drop of the SF efficiency at low surface densities. Surprisingly, we discovered an unexpected correlation between the volume densities of HI and SFR, indicating that the atomic gas is a good tracer of the cold star-forming gas, especially in low density HI-dominated environments., Comment: 24 pages, 12 figures, accepted by A&A

Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literature are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearly unrealistic because the rotational velocity $v_\phi(R,z)$ does not depend on the distance $z$ from the galactic plane, while (ii) models obtained as a result of cosmological galaxy formation simulations are more realistic, but are impractical to fit to observations due to high computational cost. Here we bridge this gap by presenting a general method to construct axisymmetric baroclinic equilibrium models of rotating galactic coronae in arbitrary external potentials. We consider in particular a family of models whose equipressure surfaces in the $(R,z)$ plane are ellipses of varying axis ratio. These models are defined by two one-dimensional functions, the axial ratio of pressure $q_{\rm axis}(z)$ and the value of the pressure $P_{\rm axis}(z)$ along the galaxy's symmetry axis. These models can have a rotation speed $v_\phi(R,z)$ that realistically decreases as one moves away from the galactic plane, and can reproduce the angular momentum distribution found in cosmological simulations. The models are computationally cheap to construct and can thus be used in fitting algorithms. We provide a python code that given $q_{\rm axis}(z)$, $P_{\rm axis}(z)$ and $\Phi(R,z)$ returns $\rho(R,z)$, $T(R,z)$, $P(R,z)$, $v_\phi(R,z)$. We show a few examples of these models using the Milky Way as a case study., Comment: Accepted for publication in MNRAS

In a $\Lambda$CDM Universe, the specific stellar angular momentum ($j_\ast$) and stellar mass ($M_\ast$) of a galaxy are correlated as a consequence of the scaling existing for dark matter haloes ($j_{\rm h}\propto M_{\rm h}^{2/3}$). The shape of this law is crucial to test galaxy formation models, which are currently discrepant especially at the lowest masses, allowing to constrain fundamental parameters, e.g. the retained fraction of angular momentum. In this study, we accurately determine the empirical $j_\ast-M_\ast$ relation (Fall relation) for 92 nearby spiral galaxies (from S0 to Irr) selected from the Spitzer Photometry and Accurate Rotation Curves (SPARC) sample in the unprecedented mass range $7 \lesssim \log M_\ast/M_\odot \lesssim 11.5$. We significantly improve all previous estimates of the Fall relation by determining $j_\ast$ profiles homogeneously for all galaxies, using extended HI rotation curves, and selecting only galaxies for which a robust $j_\ast$ could be measured (converged $j_\ast(

We derive the stellar-to-halo specific angular momentum relation (SHSAMR) of galaxies at $z=0$ by combining i) the standard $\Lambda$CDM tidal torque theory ii) the observed relation between stellar mass and specific angular momentum (Fall relation) and iii) various determinations of the stellar-to-halo mass relation (SHMR). We find that the ratio $f_j = j_\ast/j_{\rm h}$ of the specific angular momentum of stars to that of the dark matter i) varies with mass as a double power-law, ii) it always has a peak in the mass range explored and iii) it is $3-5$ times larger for spirals than for ellipticals. The results have some dependence on the adopted SHMR and we provide fitting formulae in each case. For any choice of the SHMR, the peak of $f_j$ occurs at the same mass where the stellar-to-halo mass ratio $f_\ast = M_\ast/M_{\rm h}$ has a maximum. This is mostly driven by the straightness and tightness of the Fall relation, which requires $f_j$ and $f_\ast$ to be correlated with each other roughly as $f_j\propto f_\ast^{2/3}$, as expected if the outer and more angular momentum rich parts of a halo failed to accrete onto the central galaxy and form stars (biased collapse). We also confirm that the difference in the angular momentum of spirals and ellipticals at a given mass is too large to be ascribed only to different spins of the parent dark-matter haloes (spin bias)., Comment: matches MNRAS published version

Cosmological simulations suggest that most of the matter in the Universe is distributed along filaments connecting galaxies. Illuminated by the cosmic UV background (UVB), these structures are expected to glow in fluorescent Lyman alpha emission with a Surface Brightness (SB) that is well below current observational limits for individual detections. Here, we perform a stacking analysis of the deepest MUSE/VLT data using three-dimensional regions (subcubes) with orientations determined by the position of neighbouring Lyman alpha galaxies (LAEs) at 3

Massive and diffuse haloes of hot gas (coronae) are important intermediaries between cosmology and galaxy evolution, storing mass and angular momentum acquired from the cosmic web until eventual accretion on to star-forming discs. We introduce a method to reconstruct the rotation of a galactic corona, based on its angular momentum distribution (AMD). This allows us to investigate in what conditions the angular momentum acquired from tidal torques can be transferred to star forming discs and explain observed galaxy-scale processes, such as inside-out growth and the build-up of abundance gradients. We find that a simple model of an isothermal corona with a temperature slightly smaller than virial and a cosmologically motivated AMD is in good agreement with galaxy evolution requirements, supporting hot-mode accretion as a viable driver for the evolution of spiral galaxies in a cosmological context. We predict moderately sub-centrifugal rotation close to the disc and slow rotation close to the virial radius. Motivated by the observation that the Milky Way has a relatively hot corona (T ~ 2 x 10^6 K), we also explore models with a temperature larger than virial. To be able to drive inside-out growth, these models must be significantly affected by feedback, either mechanical (ejection of low angular momentum material) or thermal (heating of the central regions). However, the agreement with galaxy evolution constraints becomes, in these cases, only marginal, suggesting that our first and simpler model may apply to a larger fraction of galaxy evolution history., Comment: MNRAS, in press. 21 pages, 10 figures, 2 tables

Gas accretion and radial flows are key ingredients of the chemical evolution of spiral galaxies. They are also tightly linked to each other (accretion drives radial flows, due to angular momentum conservation) and should therefore be modelled simultaneously. We summarise an algorithm that can be used to consistently compute accretion profiles, radial flows and abundance gradients under quite general conditions and we describe illustrative applications to the Milky Way. We find that gas-phase abundance gradients strongly depend on the angular momentum of the accreting material and, in the outer regions, they are significantly affected by the choice of boundary conditions., Comment: 4 pages, 2 figures. Proceedings of the 592 WE-Heraeus Seminar. To appear in Astronomische Nachricthen, special issue "Reconstructing the Milky Way's history: spectroscopic surveys, asteroseismology and chemodynamical models", Guest Editors C. Chiappini, J. Montalban and M. Steffen

The metal-poor gas continuously accreting onto the discs of spiral galaxies is unlikely to arrive from the intergalactic medium (IGM) with exactly the same rotation velocity as the galaxy itself and even a small angular momentum mismatch inevitably drives radial gas flows within the disc, with significant consequences to galaxy evolution. Here we provide some general analytic tools to compute accretion profiles, radial gas flows and abundance gradients in spiral galaxies as a function of the angular momentum of accreting material. We generalize existing solutions for the decomposition of the gas flows, required to reproduce the structural properties of galaxy discs, into direct accretion from the IGM and a radial mass flux within the disc. We then solve the equation of metallicity evolution in the presence of radial gas flows with a novel method, based on characteristic lines, which greatly reduces the numerical demand on the computation and sheds light on the crucial role of boundary conditions on the abundance profiles predicted by theoretical models. We also discuss how structural and chemical constraints can be combined to disentangle the contributions of inside-out growth and radial flows in the development of abundance gradients in spiral galaxies. Illustrative examples are provided throughout with parameters plausible for the Milky Way. We find that the material accreting on the Milky Way should rotate at 70-80 per cent of the rotational velocity of the disc, in agreement with previous estimates., Comment: MNRAS, in press; 16 pages, 9 Figures

We present a new and simple method to measure the instantaneous mass and radial growth rates of the stellar discs of spiral galaxies, based on their star formation rate surface density (SFRD) profiles. Under the hypothesis that discs are exponential with time-varying scalelengths, we derive a universal theoretical profile for the SFRD, with a linear dependence on two parameters: the specific mass growth rate $\nu_\textrm{M} \equiv \dot{M_\star}/M_\star$ and the specific radial growth rate $\nu_\textrm{R} \equiv \dot{R}_\star/R_\star$ of the disc. We test our theory on a sample of 35 nearby spiral galaxies, for which we derive a measurement of $\nu_\textrm{M}$ and $\nu_\textrm{R}$. 32/35 galaxies show the signature of ongoing inside-out growth ($\nu_\textrm{R} > 0$). The typical derived e-folding timescales for mass and radial growth in our sample are ~ 10 Gyr and ~ 30 Gyr, respectively, with some systematic uncertainties. More massive discs have a larger scatter in $\nu_\textrm{M}$ and $\nu_\textrm{R}$, biased towards a slower growth, both in mass and size. We find a linear relation between the two growth rates, indicating that our galaxy discs grow in size at ~ 0.35 times the rate at which they grow in mass; this ratio is largely unaffected by systematics. Our results are in very good agreement with theoretical expectations if known scaling relations of disc galaxies are not evolving with time., Comment: MNRAS, accepted. 14 pages, 4 figures, 3 tables. Additional material (Atlas.pdf) available at http://www.filippofraternali.com/downloads/index.html