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.

A prediction of the classic active galactic nuclei (AGN) unification model is the presence of ionisation cones with different orientations depending on the AGN type. Confirmations of this model exist for present times, but it is less clear in the early Universe. Here, we use the morphology of giant Ly$\alpha$ nebulae around AGNs at redshift z$\sim$3 to probe AGN emission and therefore the validity of the AGN unification model at this redshift. We compare the spatial morphology of 19 nebulae previously found around type I AGNs with a new sample of 4 Ly$\alpha$ nebulae detected around type II AGNs. Using two independent techniques, we find that nebulae around type II AGNs are more asymmetric than around type I, at least at radial distances $r>30$~physical kpc (pkpc) from the ionizing source. We conclude that the type I and type II AGNs in our sample show evidence of different surrounding ionising geometries. This suggests that the classical AGN unification model is also valid for high-redshift sources. Finally, we discuss how the lack of asymmetry in the inner parts (r$\lesssim$30 pkpc) and the associated high values of the HeII to Ly$\alpha$ ratios in these regions could indicate additional sources of (hard) ionizing radiation originating within or in proximity of the AGN host galaxies. This work demonstrates that the morphologies of giant Ly$\alpha$ nebulae can be used to understand and study the geometry of high redshift AGNs on circum-nuclear scales and it lays the foundation for future studies using much larger statistical samples.

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.

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.

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.

Several open questions on galaxy formation and evolution have their roots in the lack of a universal star formation law, that could univocally link the gas properties, e.g. its density, to the star formation rate (SFR) density. In a recent paper, we used a sample of nearby disc galaxies to infer the volumetric star formation (VSF) law, a tight correlation between the gas and the SFR volume densities derived under the assumption of hydrostatic equilibrium for the gas disc. However, due to the dearth of information about the vertical distribution of the SFR in these galaxies, we could not find a unique slope for the VSF law, but two alternative values. In this paper, we use the scale height of the SFR density distribution in our Galaxy adopting classical Cepheids (age$\lesssim 200$ Myr) as tracers of star formation. We show that this latter is fully compatible with the flaring scale height expected from gas in hydrostatic equilibrium. These scale heights allowed us to convert the observed surface densities of gas and SFR into the corresponding volume densities. Our results indicate that the VSF law $\rho_\mathrm{SFR} \propto \rho_\mathrm{gas}^\alpha$ with $\alpha \approx 2$ is valid in the Milky Way as well as in nearby disc galaxies.

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermal, kinematic, and chemical composition of the circum-galactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying up to half of the sky over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view (~1 deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5" in the soft X-ray band, and a Low-Earth orbit that ensures a stable instrumental background and high sensitivity to faint, extended emission.

The well-known quasar SDSS J095253.83+011421.9 (J0952+0114) at z=3.02 has one of the most peculiar spectra discovered so far, showing the presence of narrow Ly$\alpha$ and broad metal emission lines. Although recent studies have suggested that a Proximate Damped Ly$\alpha$ system (PDLA) causes this peculiar spectrum, the origin of the gas associated with the PDLA is unknown. Here we report the results of MUSE observations that reveal a new giant ($\approx$ 100 physical kpc) Lyman $\alpha$ nebula. The detailed analysis of the Ly$\alpha$ velocity, velocity dispersion, and surface brightness profiles suggests that the J0952+0114 Ly$\alpha$ nebula shares similar properties of other QSO nebulae previously detected with MUSE, implying that the PDLA in J0952+0144 is covering only a small fraction of the QSO emission solid angle. We also detected bright and spectrally narrow CIV$\lambda$1550 and HeII$\lambda$1640 extended emission around J0952+0114 with velocity centroids similar to the peak of the extended and central narrow Ly$\alpha$ emission. The presence of a peculiarly bright, unresolved, and relatively broad HeII$\lambda$1640 emission in the central region at exactly the same PDLA redshift hints at the possibility that the PDLA originates in a clumpy outflow with a bulk velocity of about 500 km/s. The smaller velocity dispersion of the large scale Ly$\alpha$ emission suggests that the high-speed outflow is confined to the central region. Lastly, the derived spatially resolved HeII/Ly$\alpha$ and CIV/Ly$\alpha$ maps show a positive gradient with the distance to the QSO hinting at a non-homogeneous ionization parameter distribution.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We present a study of the [OIII]5007/[OII]3727 (O32) ratios of star-forming galaxies drawn from MUSE data spanning a redshift range 0.28<z<0.85. Recently discovered Lyman continuum (LyC) emitters have extremely high oxygen line ratios: O32>4. Here we aim to understand the properties and the occurrences of galaxies with such high line ratios. Combining data from several MUSE GTO programmes, we select a population of star-forming galaxies with bright emission lines, from which we draw 406 galaxies for our analysis based on their position in the z-dependent star formation rate (SFR) - stellar mass (M*) plane. Out of this sample 15 are identified as extreme oxygen emitters based on their O32 ratios (3.7%) and 104 galaxies have O32>1 (26%). Our analysis shows no significant correlation between M*, SFR, and the distance from the SFR-M* relation with O32. We find a decrease in the fraction of galaxies with O32>1 with increasing M*, however, this is most likely a result of the relationship between O32 and metallicity, rather than between O32 and M*. We draw a comparison sample of local analogues with <z>~0.03 from SDSS, and find similar incidence rates for this sample. In order to investigate the evolution in the fraction of high O32 emitters with redshift, we bin the sample into three redshift subsamples of equal number, but find no evidence for a dependence on redshift. Furthermore, we compare the observed line ratios with those predicted by nebular models with no LyC escape and find that most of the extreme oxygen emitters can be reproduced by low metallicity models. The remaining galaxies are likely LyC emitter candidates. Finally, based on a comparison between electron temperature estimates from the [OIII4363]/[OIII]5007 ratio of the extreme oxygen emitters and nebular models, we argue that the galaxies with the most extreme O32 ratios have young light-weighted ages.

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).

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(<R)$ radial profile). We find the relation to be well described by a single, unbroken power-law $j_\ast\propto M_\ast^\alpha$ over the entire mass range, with $\alpha=0.55\pm 0.02$ and orthogonal intrinsic scatter of $0.17\pm 0.01$ dex. We finally discuss some implications for galaxy formation models of this fundamental scaling law and, in particular, the fact that it excludes models in which discs of all masses retain the same fraction of the halo angular momentum.

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<z<4. Our method should increase the probability of detecting filamentary Lyman alpha emission, provided that these structures are Lyman Limit Systems (LLSs). By stacking 390 oriented subcubes we reach a 2 sigma sensitivity level of SB ~ 0.44e-20 erg/s/cm^2/arcsec^2 in an aperture of 1 arcsec^2 x 6.25 Angstrom, which is three times below the expected fluorescent Lyman alpha signal from the Haardt-Madau 2012 (HM12) UVB at z~3.5. No detectable emission is found on intergalactic scales, implying that at least two thirds of our subcubes do not contain oriented LLSs for a HM12 UVB. On the other hand, significant emission is detected in the circum-galactic medium (CGM) of galaxies in the direction of the neighbours. The signal is stronger for galaxies with a larger number of neighbours and appears to be independent of any other galaxy properties such as luminosity, redshift and neighbour distance. We estimate that preferentially oriented satellite galaxies cannot contribute significantly to this signal, suggesting instead that gas densities in the CGM are typically larger in the direction of neighbouring galaxies on cosmological scales.

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.

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.

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.

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.