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.

We investigate the disc-halo connection in massive (Mstar/Msun>5e10) disc galaxies from the cosmological hydrodynamical simulations EAGLE and IllustrisTNG, and compare it with that inferred from the study of HI rotation curves in nearby massive spirals from the Spitzer Photometry and Accurate Rotation Curves (SPARC) dataset. We find that discrepancies between the the simulated and observed discs arise both on global and on local scales. Globally, the simulated discs inhabit halos that are a factor ~4 (in EAGLE) and ~2 (in IllustrisTNG) more massive than those derived from the rotation curve analysis of the observed dataset. We also use synthetic rotation curves of the simulated discs to demonstrate that the recovery of the halo masses from rotation curves are not systematically biased. We find that the simulations predict dark-matter dominated systems with stellar-to-total enclosed mass ratios that are a factor of 1.5-2 smaller than real galaxies at all radii. This is an alternative manifestation of the `failed feedback problem', since it indicates that simulated halos hosting massive discs have been too inefficient at converting their baryons into stars, possibly due to an overly efficient stellar and/or AGN feedback implementation.

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.

Point spread function (PSF) reconstruction (PSF-R) is a well established technique to determine reliably and accurately the PSF from Adaptive Optics (AO) control loop data. We have successfully applied this technique to improve the precision on photometry and astrometry to observation of NGC6121 obtained with SPHERE/ZIMPOL as it will be presented in a forthcoming letter. Firstly, we present the methodology we followed to reconstruct the PSF combining pupil-plane and focal-plane measurements using using our PSF-R method PRIME (Beltramo-Martin et al. 2019), with upgrade of both the model and best-fitting steps compared to previous papers. Secondly, we highlight that PRIME allows to maintain the PSF fitting residual below 0.2% over 2 hours of observation and using only 30 s of AO telemetry, which may have important consequences for telemetry storage for PSF-R purpose on future 30-40 m class telescopes. Finally, we deploy PRIME in a more realistic regime using faint stars so as to identify the precision needed on the initial guess parameters to ensure the convergence towards the optimal solution.

Precise photometric and astrometric measurements on astronomical images require an accurate knowledge of the Point Spread Function (PSF). When the PSF cannot be modelled directly from the image, PSF-reconstruction techniques become the only viable solution. So far, however, their performance on real observations has rarely been quantified. Aims. In this Letter, we test the performance of a novel hybrid technique, called PRIME, on Adaptive Optics-assisted SPHERE/ZIMPOL observations of the Galactic globular cluster NGC6121. Methods. PRIME couples PSF-reconstruction techniques, based on control-loop data and direct image fitting performed on the only bright point-like source available in the field of view of the ZIMPOL exposures, with the aim of building the PSF model. Results. By exploiting this model, the magnitudes and positions of the stars in the field can be measured with an unprecedented precision, which surpasses that obtained by more standard methods by at least a factor of four for on-axis stars and by up to a factor of two on fainter, off-axis stars. Conclusions. Our results demonstrate the power of PRIME in recovering precise magnitudes and positions when the information directly coming from astronomical images is limited to only a few point-like sources and, thus, paving the way for a proper analysis of future Extremely Large Telescope observations of sparse stellar fields or individual extragalactic objects.

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.

We present radial velocities for five member stars of the recently discovered young (age= 100-150 Myr) stellar system Price-Whelan 1 (PW 1), that is located far away in the Galactic Halo (D~ 29 kpc, Z~ 15 kpc), and that is probably associated to the Leading Arm (LA) of the Magellanic Stream. We measure the systemic radial velocity of PW 1, V_r=275 +/- 10 km/s, significantly larger than the velocity of the LA gas in the same direction. We re-discuss the main properties and the origin of this system in the light of these new observations, computing the orbit of the system and comparing its velocity with that of the HI in its surroundings. We show that the bulk of the gas at the velocity of the stars is more than 10 deg (5~kpc) away from PW 1 and the velocity difference between the gas and the stars become larger as gas closer to the stars is considered. We discuss the possibilities that (a) the parent gas cloud was dissolved by the interaction with the Galactic gas, and (b) that the parent cloud is the high velocity cloud HVC287.5+22.5+240, lagging behind the stellar system by ~25 km/s and ~10 deg.

We present a systematic study of the extraplanar gas (EPG) in a sample of 15 nearby late-type galaxies at intermediate inclinations using publicly available, deep interferometric HI data from the HALOGAS survey. For each system we mask the HI emission coming from the regularly rotating disc and use synthetic datacubes to model the leftover "anomalous" HI flux. Our model consists of a smooth, axisymmetric thick component described by 3 structural and 4 kinematical parameters, which are fit to the data via a Bayesian MCMC approach. We find that extraplanar HI is nearly ubiquitous in disc galaxies, as we fail to detect it in only two of the systems with the poorest spatial resolution. The EPG component encloses ~5-25% of the total HI mass, with a mean value of 14%, and has a typical thickness of a few kpc, incompatible with expectations based on hydrostatic equilibrium models. The EPG kinematics is remarkably similar throughout the sample, and consists of a lagging rotation with typical vertical gradients of about -10 km/s/kpc, a velocity dispersion of 15-30 km/s and, for most galaxies, a global inflow in both the vertical and radial directions with speeds of 20-30 km/s. The EPG HI masses are in excellent agreement with predictions from simple models of the galactic fountain powered by stellar feedback. The combined effect of photo-ionisation and interaction of the fountain material with the circumgalactic medium can qualitatively explain the kinematics of the EPG, but dynamical models of the galactic fountain are required to fully test this framework.

We study the gas kinematics traced by the 21-cm emission of a sample of six HI$-$rich low surface brightness galaxies classified as ultra-diffuse galaxies (UDGs). Using the 3D kinematic modelling code $\mathrm{^{3D}}$Barolo we derive robust circular velocities, revealing a startling feature: HI$-$rich UDGs are clear outliers from the baryonic Tully-Fisher relation, with circular velocities much lower than galaxies with similar baryonic mass. Notably, the baryon fraction of our UDG sample is consistent with the cosmological value: these UDGs are compatible with having no "missing baryons" within their virial radii. Moreover, the gravitational potential provided by the baryons is sufficient to account for the amplitude of the rotation curve out to the outermost measured point, contrary to other galaxies with similar circular velocities. We speculate that any formation scenario for these objects will require very inefficient feedback and a broad diversity in their inner dark matter content.

In $\Lambda$CDM cosmology, to first order, galaxies form out of the cooling of baryons within the virial radius of their dark matter halo. The fractions of mass and angular momentum retained in the baryonic and stellar components of disc galaxies put strong constraints on our understanding of galaxy formation. In this work, we derive the fraction of angular momentum retained in the stellar component of spirals, $f_j$, the global star formation efficiency $f_M$, and the ratio of the asymptotic circular velocity ($V_{\rm flat}$) to the virial velocity $f_V$, and their scatter, by fitting simultaneously the observed stellar mass-velocity (Tully-Fisher), size-mass, and mass-angular momentum (Fall) relations. We compare the goodness of fit of three models: (i) where the logarithm of $f_j$, $f_M$, and $f_V$ vary linearly with the logarithm of the observable $V_{\rm flat}$; (ii) where these values vary as a double power law; and (iii) where these values also vary as a double power law but with a prior imposed on $f_M$ such that it follows the expectations from widely used abundance matching models. We conclude that the scatter in these fractions is particularly small ($\sim 0.07$ dex) and that the linear model is by far statistically preferred to that with abundance matching priors. This indicates that the fundamental galaxy formation parameters are small-scatter single-slope monotonic functions of mass, instead of being complicated non-monotonic functions. This incidentally confirms that the most massive spiral galaxies should have turned nearly all the baryons associated with their haloes into stars. We call this the failed feedback problem.

It is commonly believed that galaxies use, throughout the Hubble time, a very small fraction of the baryons associated to their dark matter halos to form stars. This so-called low "star formation efficiency" $f_\star\equiv M_\star/f_{\rm b} M_{\rm halo}$, where $f_{\rm b}\equiv\Omega_{\rm b}/\Omega_{\rm c}$ is the cosmological baryon fraction, is expected to reach its peak at nearly $L^\ast$ (at efficiency $\approx 20\%$) and decline steeply at lower and higher masses. We have tested this using a sample of nearby star-forming galaxies, from dwarfs ($M_\star\simeq 10^7 M_\odot$) to high-mass spirals ($M_\star\simeq 10^{11} M_\odot$) with HI rotation curves and 3.6$\mu$m photometry. We fit the observed rotation curves with a Bayesian approach by varying three parameters, stellar mass-to-light ratio, halo concentration and mass. We found two surprising results: 1) the star formation efficiency is a monotonically increasing function of $M_\star$ with no sign of a decline at high masses, and 2) the most massive spirals ($M_\star\simeq 1-3 \times 10^{11} M_\odot$) have $f_\star\approx 0.3-1$, i.e. they have turned nearly all the baryons associated to their haloes into stars. These results imply that the most efficient galaxies at forming stars are massive spirals (not $L^\ast$ galaxies), they reach nearly 100% efficiency and thus, once both their cold and hot gas is considered into the baryon budget, they have virtually no missing baryons. Moreover, there is no evidence of mass quenching of the star formation occurring in galaxies up to halo masses of $M_{\rm halo}\approx {\rm a\, few}\times 10^{12} M_\odot$.

We investigate the relation between stellar mass and specific stellar angular momentum, or `Fall relation', for a sample of 17 isolated, regularly rotating disc galaxies at z=1. All galaxies have a) rotation curves determined from Halpha emission-line data; b) HST imaging in optical and infrared filters; c) robust determinations of their stellar masses. We use HST images in f814w and f160w filters, roughly corresponding to rest-frames B and I bands, to extract surface brightness profiles for our systems. We robustly bracket the specific angular momentum by assuming that rotation curves beyond the outermost Halpha rotation point stay either flat or follow a Keplerian fall-off. By comparing our measurements with those determined for disc galaxies in the local Universe, we find no evolution in the Fall relation in the redshift range 0<z<1, regardless of the band used and despite the uncertainties in the stellar rotation curves at large radii. This result holds unless stellar masses at z=1 are systematically underestimated by more than 50%. Our findings are compatible with expectations based on a LCDM cosmological framework and support a scenario where both the stellar Tully-Fisher and mass-size relations for spirals do not evolve significantly in this redshift range.

We use mock interferometric HI measurements and a conventional tilted-ring modelling procedure to estimate circular velocity curves of dwarf galaxy discs from the APOSTLE suite of {\Lambda}CDM cosmological hydrodynamical simulations. The modelling yields a large diversity of rotation curves for an individual galaxy at fixed inclination, depending on the line-of-sight orientation. The diversity is driven by non-circular motions in the gas; in particular, by strong bisymmetric fluctuations in the azimuthal velocities that the tilted-ring model is ill-suited to account for and that are difficult to detect in model residuals. Large misestimates of the circular velocity arise when the kinematic major axis coincides with the extrema of the fluctuation pattern, in some cases mimicking the presence of kiloparsec-scale density 'cores', when none are actually present. The thickness of APOSTLE discs compounds this effect: more slowly-rotating extra-planar gas systematically reduces the average line-of-sight speeds. The recovered rotation curves thus tend to underestimate the true circular velocity of APOSTLE galaxies in the inner regions. Non-circular motions provide an appealing explanation for the large apparent cores observed in galaxies such as DDO 47 and DDO 87, where the model residuals suggest that such motions might have affected estimates of the inner circular velocities. Although residuals from tilted ring models in the simulations appear larger than in observed galaxies, our results suggest that non-circular motions should be carefully taken into account when considering the evidence for dark matter cores in individual galaxies.

We study the shape and kinematics of simulated dwarf galaxy discs in the APOSTLE suite of $\Lambda$CDM cosmological hydrodynamical simulations. We find that a large fraction of these gas-rich, star-forming discs show weak bars in their stellar component, despite being dark matter-dominated systems. The bar pattern shape and orientation reflect the ellipticity of the dark matter potential, and its rotation is locked to the slow figure rotation of the triaxial dark halo. The bar-like nature of the potential induces non-circular motions in the gas component, including strong bisymmetric flows that can be readily seen as m=3 harmonic perturbations in the HI line-of-sight velocity fields. Similar bisymmetric flows are seen in many galaxies of the THINGS and LITTLE THINGS surveys, although on average their amplitudes are a factor of ~2 weaker than in our simulated discs. Our results indicate that bar-like patterns may arise even when baryons are not dominant, and that they are common enough to warrant careful consideration when analyzing the gas kinematics of dwarf galaxy discs.

The detailed distribution and kinematics of the atomic and the CO-bright molecular hydrogen in the disc of the Milky Way inside the Solar circle are derived under the assumptions of axisymmetry and pure circular motions. We divide the Galactic disc into a series of rings, and assume that the gas in each ring is described by four parameters: its rotation velocity, velocity dispersion, midplane density and its scale height. We fit these parameters to the Galactic HI and CO (J=1-0) data by producing artificial HI and CO line-profiles and comparing them with the observations. Our approach allows us to fit all parameters to the data simultaneously without assuming a-priori a radial profile for one of the parameters. We present the distribution and kinematics of the HI and H2 in both the approaching (QIV) and the receding (QI) regions of the Galaxy. Our best-fit models reproduces remarkably well the observed HI and CO longitude-velocity diagrams up to a few degrees of distance from the midplane. With the exception of the innermost 2.5 kpc, QI and QIV show very similar kinematics. The rotation curves traced by the HI and H2 follow closely each other, flattening beyond R=6.5 kpc. Both the HI and the H2 surface densities show a) a deep depression at 0.5<R<2.5 kpc, analogous to that shown by some nearby barred galaxies, b) local overdensities that can be interpreted in terms of spiral arms or ring-like features in the disk. The HI (H2) properties are fairly constant in the region outside the depression, with typical velocity dispersion of 8.9+/-1.1 (4.4+/-1.2) km/s, density of 0.43+/-0.11 (0.42+/-0.22) cm-3 and HWHM scale height of 202+/-28 (64+/-12) pc. We also show that the HI opacity in the LAB data can be accounted for by using an `effective' spin temperature of about 150 K: assuming an optically thin regime leads to underestimate the HI mass by about 30%.

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

The recent discovery of an enriched metallicity for the Smith high-velocity HI cloud (SC) lends support to a Galactic origin for this system. We use a dynamical model of the galactic fountain to reproduce the observed properties of the SC. In our model, fountain clouds are ejected from the region of the disc spiral arms and move through the halo interacting with a pre-existing hot corona. We find that a simple model where cold gas outflows vertically from the Perseus spiral arm reproduces the kinematics and the distance of the SC, but is in disagreement with the cloud's cometary morphology, if this is produced by ram-pressure stripping by the ambient gas. To explain the cloud morphology we explore two scenarios: a) the outflow is inclined with respect to the vertical direction; b) the cloud is entrained by a fast wind that escapes an underlying superbubble. Solutions in agreement with all observational constraints can be found for both cases, the former requires outflow angles >40 deg while the latter requires >1000 km/s winds. All scenarios predict that the SC is in the ascending phase of its trajectory and have large - but not implausible - energy requirements.

We use the EAGLE suite of cosmological hydrodynamical simulations to study how the HI content of present-day galaxies depends on their environment. We show that EAGLE reproduces observed HI mass-environment trends very well, while semi-analytic models typically overpredict the average HI masses in dense environments. The environmental processes act primarily as an on/off switch for the HI content of satellites with stellar mass Mstar>10^9 Msun. At a fixed Mstar, the fraction of HI-depleted satellites increases with increasing host halo mass M200 in response to stronger environmental effects, while at a fixed M200 it decreases with increasing satellite Mstar as the gas is confined by deeper gravitational potentials. HI-depleted satellites reside mostly, but not exclusively, within the virial radius r200 of their host halo. We investigate the origin of these trends by focussing on three environmental mechanisms: ram pressure stripping by the intra-group medium, tidal stripping by the host halo, and satellite-satellite encounters. By tracking back in time the evolution of the HI-depleted satellites, we find that the most common cause of HI removal is satellite encounters. The timescale for HI removal is typically less than 0.5 Gyr. Tidal stripping occurs in halos of M200<10^{14} Msun within half r200, while the other processes act also in more massive halos, generally within r200. Conversely, we find that ram pressure stripping is the most common mechanism that disturbs the HI morphology of galaxies at redshift z=0. This implies that HI removal due to satellite-satellite interactions occurs on shorter timescales than the other processes.

We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.

High-velocity clouds consist of cold gas that appears to be raining down from the halo to the disc of the Milky Way. Over the past fifty years, two competing scenarios have attributed their origin either to gas accretion from outside the Galaxy or to circulation of gas from the Galactic disc powered by supernova feedback (galactic fountain). Here we show that both mechanisms are simultaneously at work. We use a new galactic fountain model combined with high-resolution hydrodynamical simulations. We focus on the prototypical cloud complex C and show that it was produced by an explosion that occurred in the Cygnus-Outer spiral arm about 150 million years ago. The ejected material has triggered the condensation of a large portion of the circumgalactic medium and caused its subsequent accretion onto the disc. This fountain-driven cooling of the lower Galactic corona provides the low-metallicity gas required by chemical evolution models of the Milky Way's disc.

Disc galaxies like the Milky Way are expected to be surrounded by massive coronae of hot plasma that may contain a significant fraction of the so-called missing baryons. We investigate whether the local (|vLSR|<400 km/s) warm-hot absorption features observed towards extra-Galactic sources or halo stars are consistent with being produced by the cooling of the Milky Way's corona. In our scheme, cooling occurs at the interface between the disc and the corona and it is triggered by positive supernova feedback. We combine hydrodynamical simulations with a dynamical 3D model of the galactic fountain to predict the all-sky distribution of this cooling material, and we compare it with the observed distribution of detections for different `warm' (SiIII, SiIV, CII, CIV) and `hot' (OVI) ionised species. The model reproduces the position-velocity distribution and the column densities of the vast majority of warm absorbers and about half of OVI absorbers. We conclude that the warm-hot gas responsible for most of the detections lies within a few kiloparsecs from the Galactic plane, where high-metallicity material from the disc mixes efficiently with the hot corona. This process provides an accretion of a few Mo/yr of fresh gas that can easily feed the star formation in the disc of the Galaxy. The remaining OVI detections are likely to be a different population of absorbers, located in the outskirts of the Galactic corona and/or in the circumgalactic medium of nearby galaxies.

We show that the ultraviolet absorption features, newly discovered in HST spectra, are consistent with being formed in a layer that extends a few kpc above the disk of the Milky Way. In this interface between the disk and the Galactic corona, high-metallicity gas ejected from the disk by supernova feedback can mix efficiently with the virial-temperature coronal material. The mixing process triggers the cooling of the lower corona down to temperatures encompassing the characteristic range of the observed absorption features, producing a net supernova-driven gas accretion onto the disk at a rate of a few Msun/yr. We speculate that this mechanism explains how the hot-mode of cosmological accretion feeds star formation in galactic disks.

We modelled the distribution and the kinematics of HI in the inner Milky Way (R<Rsun) at latitude b=0{\deg} assuming axi-symmetry. We fitted the line profiles of the LAB 21-cm survey using an iterative approach based on the tangent-point method. The resulting model reproduces the HI data remarkably well, with significant differences arising only for R<2 kpc. This suggests that, despite the presence of a barred potential, the neutral gas in the inner Milky Way is distributed in a fairly axi-symmetric disc.

We use a model of the Galactic fountain to simulate the neutral-hydrogen emission of the Milky Way Galaxy. The model was developed to account for data on external galaxies with sensitive HI data. For appropriate parameter values, the model reproduces well the HI emission observed at Intermediate Velocities. The optimal parameters imply that cool gas is ionised as it is blasted out of the disc, but becomes neutral when its vertical velocity has been reduced by ~30 per cent. The parameters also imply that cooling of coronal gas in the wakes of fountain clouds transfers gas from the virial-temperature corona to the disc at ~2 Mo/yr. This rate agrees, to within the uncertainties with the accretion rate required to sustain the Galaxy's star formation without depleting the supply of interstellar gas. We predict the radial profile of accretion, which is an important input for models of Galactic chemical evolution. The parameter values required for the model to fit the Galaxy's HI data are in excellent agreement with values estimated from external galaxies and hydrodynamical studies of cloud-corona interaction. Our model does not reproduce the observed HI emission at High Velocities, consistent with High Velocity Clouds being extragalactic in origin. If our model is correct, the structure of the Galaxy's outer HI disc differs materially from that used previously to infer the distribution of dark matter on the Galaxy's outskirts.

Aims: we studied the global distribution and kinematics of the extra-planar neutral gas in the Milky Way. Methods: we built 3D models for a series of Galactic HI layers, projected them for an inside view, and compared them with the Leiden-Argentina-Bonn 21-cm observations. Results: we show that the Milky Way disk is surrounded by an extended halo of neutral gas with a vertical scale-height of 1.6[+0.6/-0.4] kpc and an HI mass of 3.2[+1.0/-0.9]x10^8 solar masses, which is 5-10% of the total Galactic HI. This HI halo rotates more slowly than the disk with a vertical velocity gradient of -15[+/-4] km/s/kpc. We found evidence for a global infall motion, both vertical (20[+5/-7] km/s) and radial (30[+7/-5]km/s). Conclusions: the Milky Way HI halo shows properties similar to the halos of external galaxies and is compatible with being predominantly produced by supernova explosions in the disk. It is most likely composed of distinct gas complexes with masses of 10^4-10^5 solar masses of which the Intermediate Velocity Clouds are the local manifestations. The classical High Velocity Clouds appear to be a separate population.