skip to main content

Title: Galactic Winds and Bubbles from Nuclear Starburst Rings

Galactic outflows from local starburst galaxies typically exhibit a layered geometry, with cool 104K flow sheathing a hotter 107K, cylindrically collimated, X-ray-emitting plasma. Here we argue that winds driven by energy injection in a ring-like geometry can produce this distinctive large-scale multiphase morphology. The ring configuration is motivated by the observation that massive young star clusters are often distributed in a ring at the host galaxy’s inner Lindblad resonance, where larger-scale spiral arm structure terminates. We present parameterized three-dimensional radiative hydrodynamical simulations that follow the emergence and dynamics of energy-driven hot winds from starburst rings. In this letter, we show that the flow shocks on itself within the inner ring hole, maintaining high 107K temperatures, while flows that emerge from the wind-driving ring unobstructed can undergo rapid bulk cooling down to 104K, producing a fast hot biconical outflow enclosed by a sheath of cooler nearly comoving material without ram pressure acceleration. The hot flow is collimated along the ring axis, even in the absence of pressure confinement from a galactic disk or magnetic fields. In the early stages of expansion, the emerging wind forms a bubble-like shape reminiscent of the Milky Way’s eROSITA and Fermi bubbles and can reach velocities usually associated with active-galactic-nucleus-driven winds. We discuss the physics of the ring configuration, the conditions for radiative bulk cooling, and the implications for future X-ray observations.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Medium: X Size: Article No. L24
["Article No. L24"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We analyze image and spectral data from ≈365 ks of observations from the Chandra X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the −1.1 to +0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central ∼250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and nonspherical outflow geometry. Our derived temperatures and densities yield cooling times in the nuclear region of a few million years, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne, which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates of the mass outflow rates, which are 2.8Myr−1in the northern outflow and 3.2Myr−1in the southern outflow. Additionally, we detect emission from charge exchange and find it makes a significant contribution (20%–42%) to the total broadband (0.5–7 keV) X-ray emission in the central and southern regions of the outflow.

    more » « less

    We study the properties of cosmic-ray (CR) driven galactic winds from the warm interstellar medium using idealized spherically symmetric time-dependent simulations. The key ingredients in the model are radiative cooling and CR-streaming-mediated heating of the gas. Cooling and CR heating balance near the base of the wind, but this equilibrium is thermally unstable, leading to a multiphase wind with large fluctuations in density and temperature. In most of our simulations, the heating eventually overwhelms cooling, leading to a rapid increase in temperature and a thermally driven wind; the exception to this is in galaxies with the shallowest potentials, which produce nearly isothermal $T \approx 10^4\,$ K winds driven by CR pressure. Many of the time-averaged wind solutions found here have a remarkable critical point structure, with two critical points. Scaled to real galaxies, we find mass outflow rates $\dot{M}$ somewhat larger than the observed star-formation rate in low-mass galaxies, and an approximately ‘energy-like’ scaling $\dot{M} \propto v_{\rm esc}^{-2}$. The winds accelerate slowly and reach asymptotic wind speeds of only ∼0.4vesc. The total wind power is $\sim 1~{{\ \rm per\ cent}}$ of the power from supernovae, suggesting inefficient preventive CR feedback for the physical conditions modelled here. We predict significant spatially extended emission and absorption lines from 104–105.5 K gas; this may correspond to extraplanar diffuse ionized gas seen in star-forming galaxies.

    more » « less

    Much progress has been made recently in the acceleration of ∼104 K clouds to explain absorption line measurements of the circumgalactic medium and the warm, atomic phase of galactic winds. However, the origin of the cold, molecular phase in galactic winds has received relatively little theoretical attention. Studies of the survival of ∼104 K clouds suggest efficient radiative cooling may enable the survival of expelled material from galactic discs. Alternatively, gas colder than 104 K may form within the outflow, including molecules if dust survives the acceleration process. We explore the survival of dusty clouds in a hot wind with three-dimensional hydrodynamic simulations including radiative cooling and dust modelled as tracer particles. We find that cold ∼103 K gas can be destroyed, survive, or transformed entirely to ${\sim}10^4\,$ K gas. We establish analytic criteria distinguishing these three outcomes that compare characteristic cooling times to the system’s ‘cloud crushing’ time. In contrast to typically studied ∼104 K clouds, colder clouds are entrained faster than the drag time as a result of efficient mixing. We find that while dust can in principle survive embedded in the accelerated clouds, the survival fraction depends critically on the time dust spends in the hot phase and on the effective threshold temperature for destruction. We discuss our results in the context of polluting the circumgalactic medium with dust and metals, as well as understanding observations suggesting rapid acceleration of molecular galactic winds and ram-pressure-stripped tails of jellyfish galaxies.

    more » « less
  4. Abstract Massive protostars attain high luminosities as they are actively accreting and the radiation pressure exerted on the gas in the star’s atmosphere may launch isotropic high-velocity winds. These winds will collide with the surrounding gas producing shock-heated ( T ∼ 10 7 K) tenuous gas that adiabatically expands and pushes on the dense gas that may otherwise be accreted. We present a suite of 3D radiation-magnetohydrodynamic simulations of the collapse of massive prestellar cores and include radiative feedback from the stellar and dust-reprocessed radiation fields, collimated outflows, and, for the first time, isotropic stellar winds to model how these processes affect the formation of massive stars. We find that winds are initially launched when the massive protostar is still accreting and its wind properties evolve as the protostar contracts to the main sequence. Wind feedback drives asymmetric adiabatic wind bubbles that have a bipolar morphology because the dense circumstellar material pinches the expansion of the hot shock-heated gas. We term this the “wind tunnel effect.” If the core is magnetized, wind feedback is less efficient at driving adiabatic wind bubbles initially because magnetic tension delays their growth. We find that wind feedback eventually quenches accretion onto ∼30 M ⊙ protostars that form from the collapse of the isolated cores simulated here. Hence, our results suggest that ≳30 M ⊙ stars likely require larger-scale dynamical inflows from their host cloud to overcome wind feedback. Additionally, we discuss the implications of observing adiabatic wind bubbles with Chandra while the massive protostars are still highly embedded. 
    more » « less
  5. null (Ed.)
    ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE winds may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes. 
    more » « less