skip to main content


Title: Intermediate- and high-velocity clouds in the Milky Way – I. Covering factors and vertical heights
ABSTRACT

Intermediate- and high-velocity clouds (IVCs, HVCs) are a potential source of fuel for star formation in the Milky Way (MW), but their origins and fates depend sensitively on their distances. We search for IVCs and HVCs in HST high-resolution ultraviolet spectra of 55 halo stars at vertical heights $|z|\gtrsim \,1$ kpc. We show that IVCs (40 ≤ |$v$LSR| < 90 ${\rm km\, s}^{-1}$) have a high detection rate – the covering factor, fc – that is about constant (fc = 0.90 ± 0.04) from $z$ = 1.5 to 14 kpc, implying IVCs are essentially confined to |$z$| ≲ 1.5 kpc. For the HVCs (90 ≤ |$v$LSR| ≲ 170 ${\rm km\, s}^{-1}$), we find fc increases from fc ≃ 0.14 ± 0.10 at |$z$| ≲ 2–3 kpc to fc = 0.60 ± 0.15 at 6 ≲ |$z$| ≲ 14 kpc, the latter being similar to that found towards QSOs. In contrast, the covering factor of very high-velocity clouds (VHVCs; |$v$LSR| ≳ 170 ${\rm km\, s}^{-1}$) is $f_c \lt 0.04$ in the stellar sample compared to 20 per cent towards QSOs, implying these clouds must be at d ≳ 10–15 kpc (|$z$| ≳ 10 kpc). Gas clouds with |$v$LSR| > 40 ${\rm km\, s}^{-1}$ at |b| ≳ 15° have therefore |$v$LSR| decreasing with decreasing |$z$|. Our findings are consistent with a Galactic rain and/or fountain origin for these clouds. In the latter scenario, VHVCs may mostly serve as fuel for the MW halo. In view of their high covering factors and since all the IVCs and some HVCs are found in the thick disc, they appear good candidates as gas reservoirs to help sustain star formation in the MW.

 
more » « less
PAR ID:
10366910
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
513
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
p. 3228-3240
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We use FIRE-2 simulations to examine 3D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 MW and M31-mass galaxies across their formation histories at z ≤ 1.5 ($t_{\rm lookback} \le 9.4 \, \rm {Gyr}$), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within $1 \, \rm {kpc}$ of the disc mid-plane is vertically homogeneous to $\lesssim 0.008 \, \rm {dex}$ at all z ≤ 1.5. We find negative radial gradients (metallicity decreases with galactocentric radius) at all times, which steepen over time from $\approx \! -0.01 \, \rm {dex}\, \rm {kpc}^{-1}$ at z = 1 ($t_{\rm lookback} = 7.8 \, \rm {Gyr}$) to $\approx \! -0.03 \, \rm {dex}\, \rm {kpc}^{-1}$ at z = 0, and which broadly agree with observations of the MW, M31, and nearby MW/M31-mass galaxies. Azimuthal variations at fixed radius are typically $0.14 \, \rm {dex}$ at z = 1, reducing to $0.05 \, \rm {dex}$ at z = 0. Thus, over time radial gradients become steeper while azimuthal variations become weaker (more homogeneous). As a result, azimuthal variations were larger than radial variations at z ≳ 0.8 ($t_{\rm lookback} \gtrsim 6.9 \, \rm {Gyr}$). Furthermore, elemental abundances are measurably homogeneous (to ≲0.05 dex) across a radial range of $\Delta R \approx 3.5 \, \rm {kpc}$ at z ≳ 1 and $\Delta R \approx 1.7 \, \rm {kpc}$ at z = 0. We also measure full distributions of elemental abundances, finding typically negatively skewed normal distributions at z ≳ 1 that evolve to typically Gaussian distributions by z = 0. Our results on gas abundances inform the initial conditions for stars, including the spatial and temporal scales for applying chemical tagging to understand stellar birth in the MW. 
    more » « less
  2. Abstract We report the first direct detection of molecular hydrogen associated with the Galactic nuclear wind. The Far-Ultraviolet Spectroscopic Explorer spectrum of LS 4825, a B1 Ib–II star at l , b = 1.67°,−6.63° lying d = 9.9 − 0.8 + 1.4 kpc from the Sun, ∼1 kpc below the Galactic plane near the Galactic center, shows two high-velocity H 2 components at v LSR = −79 and −108 km s −1 . In contrast, the FUSE spectrum of the nearby (∼0.6° away) foreground star HD 167402 at d = 4.9 − 0.7 + 0.8 kpc reveals no H 2 absorption at these velocities. Over 60 lines of H 2 from rotational levels J = 0 to 5 are identified in the high-velocity clouds. For the v LSR = −79 km s −1 cloud we measure total log N (H 2 ) ≥ 16.75 cm −2 , molecular fraction f H 2 ≥ 0.8%, and T 01 ≥ 97 and T 25 ≤ 439 K for the ground- and excited-state rotational excitation temperatures. At v LSR = −108 km s −1 , we measure log N (H 2 ) = 16.13 ± 0.10 cm −2 , f H 2 ≥ 0.5%, and T 01 = 77 − 18 + 34 and T 25 = 1092 − 117 + 149 K, for which the excited-state ortho- to para-H 2 is 1.0 − 0.1 + 0.3 , much less than the equilibrium value of 3 expected for gas at this temperature. This nonequilibrium ratio suggests that the −108 km s −1 cloud has been recently excited and has not yet had time to equilibrate. As the LS 4825 sight line passes close by a tilted section of the Galactic disk, we propose that we are probing a boundary region where the nuclear wind is removing gas from the disk. 
    more » « less
  3. ABSTRACT

    As part of our program to identify host galaxies of known z = 2–3 Mg ii absorbers with the Keck Cosmic Web Imager (KCWI), we discovered a compact group giving rise to a z = 2.431 DLA with ultrastrong Mg ii absorption in quasar field J234628+124859. The group consists of four star-forming galaxies within 8–28 kpc and v ∼ 40–340 km s−1 of each other, where tidal streams are weakly visible in deep HST imaging. The group geometric centre is D = 25 kpc from the quasar (D = 20–40 kpc for each galaxy). Galaxy G1 dominates the group (1.66L*, SFRFUV = 11.6 M⊙ yr−1) while G2, G3, and G4 are less massive (0.1–0.3L*, SFRFUV = 1.4–2.0 M⊙ yr−1). Using a VLT/UVES quasar spectrum covering the H i Lyman series and metal lines such as Mg ii, Si iii, and C iv, we characterized the kinematic structure and physical conditions along the line of sight with cloud-by-cloud multiphase Bayesian modelling. The absorption system has a total $\log (N({{{\rm H}\,\rm{\small I}}})/{\rm cm}^{-2})=20.53$ and an $N({{{\rm H}\,\rm{\small I}}})$-weighted mean metallicity of log (Z/Z⊙) = −0.68, with a very large Mg ii linewidth of Δv ∼ 700 km s−1. The highly kinematically complex profile is well modelled with 30 clouds across low- and intermediate-ionization phases with values ${13\lesssim \log (N({{{\rm H}\,\rm{\small I}}})/{\rm cm}^{-2})\lesssim 20}$ and −3 ≲ log (Z/Z⊙) ≲ 1. Comparing these properties to the galaxy properties, we infer a wide range of gaseous environments, including metal-rich outflows, metal-poor IGM accretion, and tidal streams from galaxy–galaxy interactions. This diversity of structures forms the intragroup medium around a complex compact group environment at the epoch of peak star formation activity. Surveys of low-redshift compact groups would benefit from obtaining a more complete census of this medium for characterizing evolutionary pathways.

     
    more » « less
  4. ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $\sim 10^{12}\, {\rm M_{\odot }}$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of ${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $\gtrsim {\dot{M}}_{\rm crit}$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare ${\dot{M}}_{\rm crit}$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that ${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $\sim 10^{12}\, {\rm M_{\odot }}$. 
    more » « less
  5. null (Ed.)
    ABSTRACT It has been proposed that gravothermal collapse due to dark matter self-interactions (i.e. self-interacting dark matter, SIDM) can explain the observed diversity of the Milky Way (MW) satellites’ central dynamical masses. We investigate the process behind this hypothesis using an N-body simulation of a MW-analogue halo with velocity-dependent SIDM (vdSIDM) in which the low-velocity self-scattering cross-section, $\sigma _{\rm T}/m_{\rm x}$, reaches 100 cm2 g−1; we dub this model the vd100 model. We compare the results of this simulation to simulations of the same halo that employ different dark models, including cold dark matter (CDM) and other, less extreme SIDM models. The masses of the vd100 haloes are very similar to their CDM counterparts, but the values of their maximum circular velocities, Vmax, are significantly higher. We determine that these high Vmax subhaloes were objects in the mass range [5 × 106, 1 × 108] M⊙ at z = 1 that undergo gravothermal core collapse. These collapsed haloes have density profiles that are described by single power laws down to the resolution limit of the simulation, and the inner slope of this density profile is approximately −3. Resolving the ever decreasing collapsed region is challenging, and tailored simulations will be required to model the runaway instability accurately at scales <1 kpc. 
    more » « less