skip to main content

Title: Dynamics of hot galactic winds launched from spherically-stratified starburst cores

The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes uniform energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as r−Δ within the starburst volume R, providing solutions for Δ = 0, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model (Δ = 0), which predicts zero velocity at the centre, for a singular isothermal sphere profile (Δ = 2), we find that the flow maintains a constant Mach number of $\mathcal {M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be written as $v_{r \lt R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where v∞ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r\lt R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the Δ = 0 (CC85) and Δ = 2 models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions.

more » « less
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society: Letters
Medium: X Size: p. L87-L91
["p. L87-L91"]
Sponsoring Org:
National Science Foundation
More Like this

    Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

    more » « less

    We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lyα emitters (LAEs) at 2.36 < z < 2.82 selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as LLyα ≃ 30 × 1042 erg s−1 and MUV ≃ −23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lyα rest-frame equivalent widths, $EW_{0}\,\rm (Ly\alpha)$, ranging from 16 to 50 Å and Lyα escape fractions, $f_{\rm esc}\, \rm (Ly\alpha)$, from 10 per cent to 40 per cent. Large $EW_{0}\, \rm (Ly\alpha)$ and $f_{\rm esc}\, \rm (Ly\alpha)$ are found predominantly in LAEs showing weak low-ionization ISM absorption (EW0 ≲ 1 Å) and narrow Lyα profiles (≲300 km s−1 FWHM) with their peak close (≲80 km s−1) to their systemic redshifts, suggestive of less scatter from low H i column densities that favours the escape of Lyα photons. We infer stellar metallicities of Z/Z⊙ ≃ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines N v1240 Å and C iv1549 Å with those from stellar synthesis models. We also find a trend between MUV and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to ≃180 M⊙ yr−1, yet they appear relatively blue (βUV ≃ −1.8 to −2.0) showing evidence of little dust attenuation [E(B − V) = 0.10–0.14]. These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

    more » « less
  3. Abstract

    M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [Cii] 158μm line in the central disk and the southern outflow of M82 using the upGREAT instrument on board SOFIA. We also report the first detections of velocity-resolved (ΔV= 10 km s−1) [Cii] emission in the outflow of M82 at projected distances of ≈1–2 kpc south of the galaxy center. We compare the [Cii] line profiles to observations of CO and Hiand find that likely the majority (>55%) of the [Cii] emission in the outflow is associated with the neutral atomic medium. We find that the fraction of [Cii] actually outflowing from M82 is small compared to the bulk gas outside the midplane (which may be in a halo or tidal streamers), which has important implications for observations of [Cii] outflows at higher redshift. Finally, by comparing the observed ratio of the [Cii] and CO intensities to models of photodissociation regions, we estimate that the far-ultraviolet (FUV) radiation field in the disk is ∼103.5G0, in agreement with previous estimates. In the outflow, however, the FUV radiation field is 2–3 orders of magnitudes lower, which may explain the high fraction of [Cii] arising from the neutral medium in the wind.

    more » « less
  4. ABSTRACT A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a ‘velocity kick’ parameter. The simplest 2cDM model with cross-section $0.01\lesssim \sigma /m\lt 1\, \textrm {cm}^{2}{ \rm g}^{-1}$ and the kick velocity $V_{\mathrm{ k}}\simeq 100\, \rm {km\, s}^{-1}$ have been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in N-body cosmological simulations tested on Milky Way (MW)-like haloes of a virial mass ${\sim}5 \times 10^{11}\, {\rm M_{\odot }}$ (Papers I & II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of ∼107−108 and ${\sim}10^{13} - 10^{14}\, {\rm M_{\odot }}$, respectively. We find that σ0/m ≳ 0.1 cm2g−1 is preferentially disfavoured for both dwarfs and galaxy cluster haloes in comparison with observations, while σ0/m = 0.001 cm2g−1 causes little perceptible difference from that of the CDM counterpart for most of the cross-section’s velocity dependence studied in this work. Our main result is that within the reasonable set of parameters, the 2cDM model can successfully explain the observational trends seen in dwarf galaxy and galaxy cluster haloes, and the model leaves us an open window for other possible alternative DM models. 
    more » « less
  5. 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