The Ionization and Dynamics of the Makani Galactic Wind

The Makani galaxy hosts the poster child of a galactic wind on scales of the circumgalactic medium. It consists of a two-episode wind in which the slow, outer wind originated 400 Myr ago (Episode I; R I = 20 − 50 kpc) and the fast, inner wind is 7 Myr old (Episode II; R II = 0 − 20 kpc). While this wind contains ionized, neutral, and molecular gas, the physical state and mass of the most extended phase—the warm, ionized gas—are unknown. Here we present Keck optical spectra of the Makani outflow. These allow us to detect hydrogen lines out to r = 30–40 kpc and thus constrain the mass, momentum, and energy in the wind. Many collisionally excited lines are detected throughout the wind, and their line ratios are consistent with 200–400 km s−1 shocks that power the ionized gas, with v shock = σ wind. Combining shock models, density-sensitive line ratios, and mass and velocity measurements, we estimate that the ionized mass and outflow rate in the Episode II wind could be as high as those of the molecular gas: MIIHII∼MIIH2=(1−2)×109M⊙ and dM/dtIIHII∼dM/dtIIH2=170−250M⊙ yr−1. The outer wind has slowed, so that dM/dtIHII∼10M⊙ yr−1, but it contains more ionized gas, MIHII=5×109 M ⊙. The momentum and energy in the recent Episode II wind imply a momentum-driven flow (p “boost” ∼7) driven by the hot ejecta and radiation pressure from the Eddington-limited, compact starburst. Much of the energy and momentum in the older Episode I wind may reside in a hotter phase, or lie further into the circumgalactic medium.


INTRODUCTION
The circumgalactic medium (CGM) of galaxies contains at least 80% of the baryonic mass within their dark matter halos (Shull et al. 2012;Werk et al. 2014;Tum-linson et al. 2017).The CGM is also home to >50% of the metals within a halo (Peeples et al. 2014;Bordoloi et al. 2014;Oppenheimer et al. 2016).The origin of this gas is likely diverse, but an important source is metalenriched gas ejected from galaxies by galactic winds.
Catching galactic winds in the act of depositing gas in the circumgalactic medium has proven challenging, as the observed sizes of these winds are tpyically com-parable to the scales of the galaxies themselves (of order 10 kpc).We reported in Rupke et al. (2019) a 100 kpc nebula surrounding the Makani galaxy at z = 0.459.The galaxy Makani (SDSS J211824.06+001729.4) is a massive (M = 10 11.1 ) star-forming galaxy with r e = 2.5 kpc (Sell et al. 2014).The nebula, observed in [O II] 3726, 3729 Å, is consistent with two galactic wind episodes over the past ∼400 Myr, based on analysis of its morphology, kinematics, and stellar populations.Episode I was powered by a star formation episode 400 Myr in the past, and includes most of the outer 20-50 kpc of the wind.This wind has slow projected speeds ∼100 km s −1 and linewidths σ = 200 km −1 , with a shape characteristic of a giant, bipolar outflow.The star formation timescale, projected ballistic flow speed, and radius of this flow are consistent: v I ∼ R I /t * ,I = 50 kpc/400 Myr = 120 km s −1 .Episode II was powered by star formation a mere 7 Myr ago, and consists of a fast wind with maximum speeds exceeding 2000 km s −1 , similar to the extended, high-velocity flow seen in other compact starbursts (Geach et al. 2014).Most of the Episode II wind is within 20 kpc of the host galaxy, though there is a faint southern extension to 40 kpc.As with Episode I, the approximate timescale, speed, and size line up: v II ∼ R II /t * ,II = 10 kpc/7 Myr = 1400 km s −1 .The overall size of the Episode I+II nebula relative to the parent galaxy (r wind /r e, * 20) is direct evidence that the wind has moved into the galaxy's CGM.
The KCWI observations covered only blue wavelengths and emission lines from [O II], Mg II 2796, 2803 Å, and [Ne V] 3345, 3426 Å.They thus leave unanswered two important questions regarding the ionized gas in the Makani wind.First, the ionized mass in the wind cannot be determined without recombination-line measurements.Second, the ionization state of the wind is unknown.The former is critical for constraining the impact of the wind on its host and surroundings; the second is important for understanding its interaction with the CGM.An understanding of how the nebula is powered may also inform how it is driven.
The spatially-resolved molecular component of the Episode II wind in Makani is outflowing with similar velocities to the ionized gas at a rate of 245 M yr −1 (Rupke et al. 2019).This suggests a rapid exhaustion of star formation fuel, at a level comparable to the star formation rate of several hundred M yr −1 (Petter et al. 2020).The ionized mass of either wind episode is unknown, however, and since no other gas phases have yet been observed in the extended Episode I wind, we do not know its outflow rate.In its time-resolved structure, Makani presents a unique opportunity to constrain the evolution of outflow rates with time.
Due to its compact size, most of the host galaxy itself is unresolved in ground-based observations.Based on Keck/NIRSPEC, MMT longslit, and SDSS data, we know the unresolved host emission-which includes a narrow component near systemic and a broad, blueshifted component that dominates the flux-lies in the AGN area of excitation diagrams (Rupke et al. 2019;Perrotta et al. 2021).The outflowing component in the galaxy also shows weak [Ne V] 3324, 3426 Å emission (Rupke et al. 2019).Normally these would be indicative of an AGN, but Makani shows no other evidence of an AGN (Sell et al. 2014;Rupke et al. 2019).We have tentatively interpreted these as evidence of high-velocity shocks as the Episode II wind propagates through the gas in the center of the galaxy.However, it is also possible that the ionization of the Makani host galaxy and that of its parent sample (Tremonti et al. 2007) reflect leakage of Lyman continuum (LyC) photons (Perrotta et al. 2021).Further constraining its ionization may distinguish between these scenarios.
To address these questions, here we present Keck/ESI long-slit spectra of Makani in the optical.In Section 2, we present the observations, data reduction, and data analysis.We describe the results in Section 3, and discuss them in the light of shock and outflow models in Section 4. We summarize in the final Section 5.

OBSERVATIONS AND DATA PROCESSING
We observed Makani on 12 Sep 2020 UT with the Echellette Spectrograph and Imager (ESI; Sheinis et al. 2002) on Keck II.We chose the 1. 0 slit in echellette mode, which yields moderate velocity resolution (R = 4000).The broad, simultaneous wavelength coverage allows us to observe many rest-frame optical strong lines ([O II] through [S II] 6716, 6731 Å), while the 20 slit length is well-matched to the extent of the 17 nebula.
Conditions were photometric, with seeing 0. 8.We took 3×30 minute exposures centered on the galaxy with PA = 45 • E of N, and 4×30 minute exposures centered on a point offset 3. 3 W at PA = 0 • (Figure 1).These slit positions were designed to probe the highest-surface brightness regions of the nebula outside of the nucleus while catching significant areas of low surface-brightness features.To calibrate, we took bias, arc lamp, dome flat, twilight flat, and pinhole flat observations.We dithered ±1 along the slit between exposures.We observed the flux standard BD+28 4211 for photometric and telluric corrections.Airmasses were 1.06-1.07for the flux standard and PA = 45 • slit; and 1.1-1.4 for the PA = 0 • slit.We reduce the data using the ESIRedux package (Prochaska et al. 2003) in XIDL.We use all calibrations but the twilight flats.We additionally telluric-correct the spectra using the flux standard after normalizing it by a 82,000 K blackbody.The final spectra are in vacuum wavelengths.
Along each slit, we extract five apertures (Figure 1 and Table 1, labelled in order of increasing projected galactocentric radius).The lengths of the inner apertures are 2×FWHM seeing (10 pixels, or 1. 54), and those of the outer apertures are 4×FWHM seeing (20 pixels, or 3. 08).We coadd the two overlapping apertures (ap7).We searched for emission in apertures outside the previously detected nebula but found none that surpassed the detection threshold.
We extract a 3 ×1 rectangular aperture centered on the galaxy to compare to the existing 3 diameter SDSS spectrum.The SDSS spectrum is 26% higher in flux over 6500-9000 Å, but the difference rises smoothly with decreasing wavelength to 80% at 4000 Å.This flux difference is also observed in a 3 circular extraction from the KCWI data cube.We apply this upward correction to each extracted spectrum.The effect on flux ratios is minimal, but may be 5% for [O II]/Hα, as the total upward correction at [O II] is 33%.
We fit the spectra using IFSFIT (Rupke 2014).We fit the continuum with pPXF (Cappellari 2012(Cappellari , 2017) ) and the same C3K Solar-metallicity stellar models (Conroy  et al., in prep.)used in Rupke et al. (2019).We use the best-fit stellar model to the host spectrum (ap1) as a template for the two r ∼ 10 kpc apertures (ap2 and ap3), which also contain some stellar emission.In these fits, we allow for non-zero stellar attenuation and include an order-4 additive polynomial to account for uncertainties in data reduction (flux calibration, scattered light, etc.).The only aperture with a non-zero best-fit stellar attenuation is ap1, for which E(B − V ) * = 0.3.There is no visible stellar continua in other apertures, though very low levels of continuum flux may indicate imperfect sky subtraction in some cases.To remove this residual signal, we fit an order-10 polynomial continuum.
After subtracting the best-fit continuum model, we fit two velocity components-one narrow and near systemic, the other broad and blueshifted-to the emission lines in ap1, ap2, and ap3, and one component otherwise (Figure 2).(Hereafter we refer to the narrow, near-systemic component of ap1 as ap1.1 and the broad, blueshifted component of ap1 as ap1.2.)These components are robustly detected in the host galaxy and we treat them separately in some parts of the analysis.While multiple components (including the same broad line seen in ap1) are clearly required in ap2 and ap3, the lower S/N in some lines means that these components are not significantly detected in every line.For instance, [S II] suffers from lower throughput and sky line contamination.Thus, in our analysis we treat only  the total flux for ap2 and ap3.For these apertures, we also fix the [O II] 3729/3726 flux ratios to match that of ap1.
For the most part, the velocities and linewidths of all emission lines are tied together.However, we find that the fit to [O I] in apertures ap1, ap4, and ap5 is improved if we fit its velocity and linewidth separately.(In other apertures, the line is too faint to fit separately.)The most significant difference is that in ap1, the broad component in [O I] is redshifted compared to the broad component in other lines by ∆v = 350 km/s and narrower by ∆σ = 200 km/s (40% smaller).In the narrow component of ap1 and in ap4 and ap5, [O I] is only slightly redshifted ( 50 km/s) and the linewidth difference is small ( 30 km/s).
In several apertures (ap2, ap7, and ap9) we first fit [O II] only and use the resulting line centers and linewidths as fixed priors for the other lines.As this line is the brightest and highest S/N, this serves to improve the fit for faint lines ([O I]), lines impacted by sky contamination (in the red), and lines with some degeneracy due to blending (Hα and [N II]).
In ap9, Hα is only marginally detected, as it is impacted by sky lines.For this aperture, we instead estimate Hα from Hβ assuming no extinction and the Case B ratio of 2.86.
There is residual interstellar Na I D absorption and emission in two apertures, ap1 and ap2.We fit these features also using IFSFIT (Figure 3), and report the results in Section 3.2.

RESULTS
The high throughput and long wavelength coverage of ESI allow us to detect numerous strong emission lines across the nebula.The ap1 spectrum also includes weak emission lines such as [Ne V] 3426 Å, [O III] 4363 Å, and [N II] 5755 Å.We report observed and extinctioncorrected line fluxes of all emission lines detected at >2σ in Table 2, as well as Hα fluxes and luminosities.Observed fluxes have been corrected for a Galactic extinction of E(B −V ) = 0.0684 (Schlafly & Finkbeiner 2011).
The ESI spectra cover the full range of Balmer lines.In Figure 4, we show observed [O II] and Hα surface brightnesses as a function of radius.The ratio steadily rises from about 1/3 in the center, to unity at 20 kpc, and finally to a maximum of 2 at r = 30 − 40 kpc.This increase in observed [O II]/Hα is due to a combination of decreasing extinction and increasing collisional excitation with increasing radius.We compute E(B − V ) from Hα/Hβ for Case B conditions using the Cardelli et al. (1989) extinction curve and R V = 3.1.E(B − V ) decreases with increasing radius from a peak of 1.0 at r = 10 kpc to 0 at r > 25 kpc (Figure 4).The intrinsic [O II]/Hα flux ratio is near unity in ap1 and ap3, rising to a value of 2-5 in larger-radius apertures.
By combining the ESI Hα and KCWI [O II] measurements, we can estimate the spatially-resolved and integrated Hα luminosity.We do so by bootstrapping from the ESI measurements to estimate how the ratio ([O II], observed)/(Hα, intrinsic) changes with radius.We fit a linear relationship between 0 and 35 kpc (Figure 4) with an RMS of 0.3 dex.We then apply this relationship to the Voronoi-binned KCWI data to infer the spatiallyresolved Hα flux.The resulting total Hα luminosity is L tot Hα = (9.4+6.0 −3.0 ) × 10 42 erg s −1 , where we calculate the error by applying the RMS dispersion at r > 5.5 kpc (the edge of ap1).(Using instead a parabolic fit increases the luminosity by 15%, but does not change the RMS and results in line ratios at the edge of the windr ∼ 50 kpc-that appear unphysical.) We calculate the electron density using the flux ratios [O II]3729/3726 and [S II]6716/6731, which are allowed to freely vary between physical limits (Sanders et al. 2016).Both lines are observed throughout the nebula, though [S II] is not detected in all apertures.The [O II] doublet is one of the brightest lines in the rest-frame optical spectrum of the nebula and is unaffected by sky lines.However, the lines are separated by only 3 Å, leading to covariance among fit parameters.Overlap of multiple velocity components adds further degeneracy.While the [S II] lines are better separated in wavelength (15 Å), they are impacted by sky lines and fall in a lowersensitivity part of the data.Thus, we report density estimates only in ap1, ap4, and ap5; the significance of the individual line detections of the [S II] doublet, and thus our ability to reliably decompose them, is too low in other apertures, even if the total doublet is detected.

Excitation of the wind
The detection of multiple emission lines throughout the Makani nebula allows us to constrain the spatiallyresolved excitation of the wind using standard line ratios (Figure 5; Baldwin et al. 1981;Veilleux & Osterbrock 1987;Cid Fernandes et al. 2010).
Beginning with the nucleus, ap1, we that find that the broad component dominates its line flux and is securely in the AGN region of these diagrams.This is consistent with previous results (Perrotta et al. 2021) Perrotta et al. (2021) interpreted these line ratios as ionization from fast shocks; we discuss  [NeV]3426 -1.02±0.12-1.20±0.10-0.43±0.17this in detail in Section 4. Makani lacks of evidence for an AGN in the rest-frame UV, X-ray, radio, or infrared (Diamond-Stanic et al. 2012;Sell et al. 2014;Petter et al. 2020;Whalen et al. 2022).The rest-frame optical-UV continuum of Makani shows no indication of non-stellar emission, and a stellar population fit to the SED is perfectly consistent with the data (Rupke et al. 2019).Though [Ne V] emission is present, its luminosity is low compared to AGN (Rupke et al. 2019).We thus discount AGN photoionization as the origin of these line ratios.
The narrow component of ap1 is instead likely ionized by stars, with some possible contribution from shocks based on its location in the composite region of the Taken as a whole, the apertures other than ap1 lie largely in the composite region of the [N II]/Hα diagram and in the low-ionization nuclear emission-line region (LINER) area in the other three diagrams.An exception is aperture ap2, which lies in between the properties of ap1 and ap3.This aperture is in the direction of the most highly blueshifted molecular and ionized gas (Rupke et al. 2019).It must therefore include more flux from this higher-ionization wind component (mixed in with a lower-ionization component) than ap3, which is at the same radius as ap2.
We detect the weak emission lines  (Nagao et al. 2001;Molina et al. 2018).Similarly, the ratio [N II] 5755/6548 is in the range -1.3 to -1.4, also consistent with T ∼ 2 × 10 4 K. Makani may be the highest-redshift system in which interstellar Na I D absorption and/or emission has been seen in a galactic wind.Resonant Na I D absorption is detected in the host galaxy (ap1) with a low equivalent width (rest-frame W eq of 0.37 +0.07 −0.05 Å), but blueshifted by -342 km s −1 (Figure 3).The outflow was not previously detected in resonant absorption in Mg II 2796, 2803 Å or Mg I transitions in the near-UV (Rupke et al. 2019).However, there is residual, weak Fe II 2585 Å absorption (Rupke et al. 2019;Perrotta et al. 2021).
Several lines of evidence suggest that resonant absorption in ap1 is filling in the emission, reducing the observed |W eq | and increasing |v − v sys | of both absorp-tion and emission.First, a fit to Fe II 2585 Å gives a velocity at maximum depth of -100 km s −1 and a covering factor of 0.12, assuming it is optically thick.Thus v NaI − v FeII ∼ −200 km s −1 .Furthermore, if the Na I D and Fe II absorption lines are optically thick, the ratio of the observed covering factors is C f,FeII /C f,NaI ∼ 3 − 4, suggesting that Na I D is too shallow.
Second, neutral gas absorption in outflows is closely connected with the foreground dust column (Rupke & Veilleux 2015;Rupke et al. 2021).Though there is significant scatter, the observed E(B−V ) = 0.6−1.0 in ap1 and ap2 (Figure 4) are typically associated with 5-10× higher W eq on average.
Finally, a comparison to the radial surface-brightness profile of Mg II (Figure 6) shows that the Na I D surface brightness declines much more slowly than that of Mg II.If this is due to infilling of emission by foreground absorption, Na I D should intrinsically be brighter in ap1 by a factor of a few (0.6 ± 0.2 dex).This in turn would raise the absorption equivalent width by the same factor and the covering factor by some (smaller) amount.This reduction of Na I D flux due to infilling is observed in other spatially-resolved observations (Rupke & Veilleux 2015).
Due to this infilling, we do not attempt to estimate the column density, and thus the mass and dynamics, of the neutral phase of the wind.Makani contains a compact (r ∼ 400 pc) stellar core (Sell et al. 2014;Rupke et al. 2019;Diamond-Stanic et al. 2021) with a young stellar population (Rupke et al. 2019).The low-velocity component of the innermost (r 5.5 kpc) aperture, ap1.1, thus contains most or all of the star formation activity.Based on the observed line ratios (Section 3.1 and Perrotta et al. 2021), this component is photoionized by the young stellar population.The Hα luminosity of this component is 2.1 × 10 42 erg s −1 (Table 2), which is 40% of the total luminosity in ap1 (i.e., the host galaxy).We note that this luminosity, however, is larger than the difference between the total nebular luminosity (Section 3) and wind luminosity (Section 4.4).This difference is L starburst Makani has a star formation rate (SFR) of 224-300 M yr −1 , as inferred from radio and infrared data (Petter et al. 2020).The radio emission is concentrated in a 1 point source coinciding with ap1.Using a standard calibration (Moustakas et al. 2006), we compute from L starburst Hα a SFR of 14 M yr −1 .This is a factor of 16 smaller than the lowest previous estimate from IR data.A similar discrepancy is seen in other galaxies in the parent galaxy sample of Makani (Edmonds et al. 2022).It may in fact be evidence of the stellar feedback in Makani shutting off star formation, as Hα traces star formation over the shortest timescales 6 Myr (Calzetti 2013).This is just below the estimated 7 Myr estimate of the most recent burst of star formation.An earlier burst occurred 400 Myr ago (Rupke et al. 2019), which is unlikely to be reflected in current SFR estimates.The observed discrepancy could also result in part from LyC leakage (Moustakas et al. 2006), to which compact starbursts with strong feedback like Makani may be susceptible (Perrotta et al. 2021).
Finally, it may also be the case that the bulk of the star formation is obscured by larger columns than traced by optical light, as seen in nearby compact mergers (Spoon et al. 2007).Future rest-frame far-ultraviolet and mid-infrared spectra of Makani and its parent sample could distinguish among these possibilities.
The oxygen abundance of this component is 12+log(O/H) ∼ 9.1 (or ∼2× Solar; Asplund et al. 2021), based on standard [N II]/[O II] and R 23 strong-line calibrations (Kewley & Dopita 2002).For component 1 of ap1, [N II]/[O II] = −0.12 and R 23 = −0.13.This is consistent with the high mass of the galaxy (Rupke et al. 2019;Diamond-Stanic et al. 2021).The actual abundance could be closer to Solar, as strong-line calibrations may overestimate the true gas abundance by a factor of ∼2 (Kewley & Ellison 2008).The weak lines [O III] 4363 Å and [N II] 5755 Å could in principle also provide a weak-line measurement based on electron temperature.However, the inferred temperatures are much too high for high-metallicity H II regions, suggesting that component blending or simply an overestimate of the ap1.1 component in these weak lines is producing these high ratios.

Comparison to shock models
Line ratios are commonly used to distinguish between stellar and shock ionization in galactic winds (Lehnert & Heckman 1996;Veilleux & Rupke 2002;Sharp & Bland-Hawthorn 2010).Aside from ap1.1, the line ratios observed in Makani (Figure 5) are inconsistent with those predicted by stellar photoionization models (Kewley et al. 2001;Byler et al. 2017).An obvious place to turn instead is shock models, both due to their typical consistency with LINER-like line ratios (e.g., Dopita et al. 1997, among numerous examples) and because of the extreme outflows we detect in the Makani nebula.
We first compare the position of apertures ap3, ap4, and ap5 in Figure 5 to the low-velocity range of fast shock models (Allen et al. 2008).We assume solar metallicity, and the benchmark model pre-shock density n H = 1 cm −3 (their model series M n1).In the [N II]/Hα diagram, shock velocities v s ∼ 200 km s −1 and low magnetic field strengths B/n 1/2 < 0.5 µG cm −3/2 can reproduce the results, whether the emission comes from the shock only or from the shock and precursor.The precursor is the gas that is about to be impacted by the shock, and which in this case is fully pre-ionized by the shock radiation.For [S II]/Hα, the same low v s and low B fit best in the shock-only grid, but the shock+precursor grid does not overlap the observed values.The [O I]/Hα grids allow for shock or shock+precursor solutions, though with slightly higher velocity, v s ∼250 km s −1 .
The fast shock models cannot match the low observed ionization parameters from [O III]/[O II], even in the lower-ionization shocked region.This may point to only partial pre-ionization of the precursor, as assumed in slow shock models (Rich et al. 2010;Dopita & Sutherland 2017), or additional post-shock physics.The slow-shock models are consistent with the observed [S II]/Hα and [O I]/Hα at the highest model velocities, v s ∼ 200 km s −1 , but produce values of [N II]/Hα that are too high at the same velocities.
The fast shock models are mostly consistent with the observed values of [O II]/Hα = 1 − 5 (Figure 4).For shocked gas only, the model ratio is 2-4 for v s = 100 − 300 km s −1 , with a minimum around 200 km s −1 .Including precursor, the model ratio is slightly smaller, in the 1.8-2.8range.The observed values are much larger than those in nearby star-forming galaxies, for which the maximum observed [O II]/Hα is about 2 (Moustakas et al. 2006).It is also larger than observed in the diffuse ionized gas in extraplanar regions of nearby galaxies (Jones et al. 2017).
We estimate the emitted Hβ fluxes, F shock Hβ , across the shock in the context of these models by dividing the luminosity in each aperture (Table 2) by the projected physical area of each aperture in cm −2 (Table 1).This assumes the shock is plane-parallel to the line-of-sight, and so may be an underestimate if the aperture is tracing the edge of the wind.However, given the large size of the wind, projection effects are probably not significant for the inner apertures (ap1-ap5).The resulting shock fluxes F shock Hβ in ap3-ap5 are approximately 10 −4.0 to 10 −3.7 erg s −1 cm −2 .We then compare these observed fluxes to those predicted by the the Z = Z models M n1 and V n10 from Allen et al. (2008).We use the shock velocity and magnetic field parameter estimated from comparing the observed line ratios to the same models, v s = 200 km s −1 and B/n 1/2 < 0.5 µG cm −3/2 .(The result is not sensitive to B.) Finally, interpolating between the model grids, the extrapolated pre-shock density from fast models is then n H = 2 − 3 cm −3 , since the flux scales linearly with n H (Dopita & Sutherland 1995).At larger radius, the Hβ fluxes decrease to the range 10 −4.3 to 10 −5.7 erg s −1 cm −1 , which suggests an increase in projection effects that mean these are lower limits; or a lower velocity and/or density at these radii.Lower velocity would be consistent, for instance, with the lower line ratios in ap6 and ap8 compared to ap3-ap5, but it seems probable that the gas density would also decrease with radius.
At r = 0-5 kpc, the high-velocity outflow component ap1.The Hβ flux is also higher at 10 −2.65 erg s −1 cm −1 .The line ratios place this broad component in the AGN regions of the excitation diagrams.However, n H ∼ 10 cm −3 shock+precursor models from Allen et al. (2008) (their models V n10) are entirely consistent with the observed line ratios if v s is in the range 300-400 km s −1 and B/n 1/2 < 0.5 µG cm −3/2 .These models also reproduce the observed ratios [Ne V] 3426/Hβ = (−1.33±0.20)dex and [Ne V] 3426/[Ne III] 3869 = (−0.45± 0.22) dex.The higher ionization state is due to the harder radiation field produced by the higher-velocity shocks.This harder radiation also produces high temperatures in the precursor and post-shock regions.The observed [O III] 4363/5007 ratio in ap1.2 is -1.1 dex, consistent with these high temperature (T 2 × 10 4 K; Section 3).Such ratios are commonly observed in LIN-ERs (Nagao et al. 2001;Molina et al. 2018).However, the models cannot perfectly reproduce the combined [O III] 4363/5007 and [O III]/Hβ ratios, an ongoing problem (Dopita & Sutherland 1995;Allen et al. 2008).
In summary, the fast shock models of Allen et al. (2008) are very consistent with the observed strong line ratios and fluxes throughout the nebula.In ap1.2 (the outflow at radii 0-5 kpc), high shock velocities (300-400 km s −1 ) and an ionized precursor are needed.These shock velocities align with the observed σ II = 400 km s −1 (Rupke et al. 2019).Farther from the host, the required shock velocities are lower, 200 km s −1 , but these still align with the observed linewidths σ I = 200 km s −1 .The precursor emission contributes fractionally less at these velocities.In the next section, we apply the shock models to constrain the wind density.2008) shock models to constrain the ambient, pre-shock and post-shock densities in the wind and extended nebula.The high-velocity Episode II wind seen in ap1-ap3 has a high electron density n e = 3000 cm −3 (1σ range 1200-10,000 km s −1 ) from the [S II] ratio of ap1.2.This density is also consistent with the [O II] ratio.The shock model that best fits the observed line ratios and shock fluxes in ap1.2-ap3 has a pre-shock density of n H ∼ 10 cm −1 , as we discuss above in Section 4.2.However, compression in the shock raises the density in the post-shock region up to a factor of n post /n pre = (8πµm H ) 1/2 v shock /(B/n 1/2 ) (Dopita & Sutherland 1996).For the estimated velocities of 300-400 km s −1 and low magnetic field parameter 0.5 µG cm −3/2 , the maximum value of log(n post /n pre ) is 2.75.This translates to a post-shock density of at most n post ∼ 6000 cm −3 for n pre = 10 cm −3 .(The compression can also be seen in Figures 3-6 of Allen et al. (2008).) The correspondence between the observed density and predicted post-shock density suggests that the observed [S II] and [O II] emission is arising primarily in the compressed post-shock region.The line ratios, however, appear to require an ionized precursor.That the postshock region seems to dominate the flux of these lines could point to additional physics occurring behind the shock front, perhaps due to the accumulation of sweptup interstellar gas or to the driving mechanism of the wind.
At larger radii, in the extended, Episode I nebula, the line ratios and fluxes point to much lower ionized gas densities n e 10 cm −3 .The shock models are consistent with low pre-shock densities of n H 2 − 3 cm −3 , with this value possibly decreasing with increasing radius.The corresponding maximum post-shock density is n H 600 cm −3 .In contrast to the high-velocity, highdensity Episode II gas, this implies that the [S II] ratio is primarily tracing the lower-density precursor gas.The line ratios are, for the most part, consistent with either shock-only or shock+precursor models.
The 200 and 400 km s −1 models with n H = 1 − 10 cm −3 (models M n1 and V n10) predict that just a little over half of the Hβ luminosity is produced in the post-shock gas, while the other fraction arises in the precursor, with little dependence on magnetic parameter B/n 1/2 (Allen et al. 2008).Thus, in the calculations that follow, we assume complete ionization in each shock region and divide the flux evenly between these low-and high-density regions.Using the estimated preshock densities and calculated maximum compression as a guideline, we assume n e = 10 and 1000 cm −3 preand post-shock for Episode I, and n e = 1 and 100 cm −3 for Episode II.The post-shock densities could be even higher, but higher values do not significantly change the results, as the lower-density pre-shock gas dominates the mass.

Mass, Momentum, and Energy
As the KCWI observations did not cover strong H recombination lines, we were previously unable to measure the mass of the ionized Makani nebula.Here we use our bootstrapped estimate of the spatially-resolved Hα luminosity of the wind (Section 3) to estimate the mass and dynamics of both the Episode I and II winds.
We start with the Voronoi-binned, KCWI [O II] map and spatially delineate the Episode I and II winds using a velocity cut at v 98% = −700 km s −1 (Rupke et al. 2019).To do this, we characterize multi-component line fits in the KCWI data with the cumulative velocity distribution.We calculate velocities that enclose a particular percentage p of flux integrated from the red side of the line, v p% .(E.g., v 98% is the velocity at which 98% of the line flux is redshifted from this velocity, and so characterizes the maximum blueshift observed.)We compute the Hα flux and luminosity in each bin, applying our model for ([O II], observed)/(Hα, intrinsic) vs. radius and using the mean radius in each Voronoi bin.From the inner r < 5.5 kpc, we also conservatively remove 40% of the flux in each bin for the contribution from star-forming gas (Section 4.1).
As discussed in Section 4.3, we use the observed line ratios and shock models to infer how much of the observed luminosity arises from regions of different density.Using these luminosities and densities, we then calculate the ionized gas mass in each bin.This ionized density model differs from many estimates in the literature, which typically use only the density computed from line ratios.However, the shock models provide further infor- Note-Properties of the ionized wind for each episode are listed for three different methods of calculation, as described in Section 4.4.Methods 1a and 1b make use of the outflow dynamical timescale to calculate the mass, momentum, and energy outflow rates, while method 2 is based on the stellar population ages for episodes I and II.For Episode I, we assume a wind radius RI = 40 kpc, electron density ne = 1 − 100 cm −3 , and time when the wind left the starburst (for Method 2) equal to the stellar population age t * ,I = 400 Myr.For Episode II, we assume RII = 15 kpc, ne = 10 − 1000 cm −3 , and t * ,II = 7 Myr.We also remove 40% of the flux at r < 5.5 kpc to conservatively account for star formation in the inner nebula (Section 4.1), which is equivalent to removing 1.4 × 10 42 erg s −1 from Episode II.Our preferred methods are 1a and 2 for Episode I; and 1b and 2 for Episode II.The 1σ errors propagated from our bootstrapping method (Section 3) are included for LHα, and similar fractional errors apply to other quantities calculated from LHα. Systematic errors due to our choice of density model, etc., are not included.
mation on the density structure that could yield a more accurate estimate.Given that M HII ∼ n −1 e , the masses are largely dependent on the lower, pre-shock density in our model.This lower density is mostly consistent with the measured densities (Section 3).However, the preshock density in the Episode II shock model is lower than the density measured in the densest part of the wind, in ap1.2, by a factor of ∼100.Assuming a constant density wind with the higher density measured from [S II] and [O II] in ap1.2 would result in a mass outflow rate that is lower by a similar factor.
The resulting Hα luminosities and gas masses are listed in Table 3.The Episode I gas is fainter but contains more mass because of our density model.Together, the total ionized gas mass of 6.4 +5.2 −2.6 × 10 9 M is comparable to the integrated CO mass of 1 × 10 10 M (Rupke et al. 2019).The detected CO is much less extended than the ionized gas (reaching a radius of only 20 kpc; Figure 6).The Episode II ionized gas mass is also comparable to the CO mass of (2.4 ± 0.6) × 10 9 M in a similar velocity range (Rupke et al. 2019).
To compute outflow rates (Table 3), we use two methods.First, we assume a single-radius wind for each episode, using R I = 40 kpc and R II = 15 kpc (Method 1).We use both the central velocity v 50% (Method 1a) and maximum velocity v 98% (Method 1b) to compute the wind properties.We deproject the velocity (either v 50% or v 98% ) in each bin to compute the 3D v rad , using the projected galactocentric distance of the bin R obs and assuming that it lies on the wind surface at the assumed 3D radius R wind : We then compute the ballistic flow time for each bin, t flow = v rad /R wind and divide the mass in that bin by t flow .We follow a similar procedure for momentum and energy, adding in velocity dispersion for energy (following Rupke et al. 2005).
As an alternative method we assume the stellar population ages t * for Episodes I and II, 400 and 7 Myr, represent the times when all of the gas in each episode left the starburst (Method 2).We then divide the ionized mass by these times for each episode to compute dM/dt = M/t * .This method is effectively a timeaveraged outflow rate, where the time average is over the lifetime of the outflow.We use it as a sanity check on Method I.
For Episode I, Method 1a (using v 50% ) yields a result consistent with Method 2, as the mean deprojected radial velocity v rad = 90 km/s is close to the velocity inferred from R I /t * ,I = 40 kpc/400 Myr = 97 km/s.In other words, most of the gas is at a velocity consistent with the outer radius of the nebula and the ballistic flow velocity, if this gas left the starburst 400 Myr ago.The gas that began at higher velocity has either slowed or since escaped to much larger radius.The remaining detected amounts of high-velocity gas do not contribute significantly to the mass, or mass outflow rate, of the outflow.Using v 98% to trace the radial velocity of the bulk of the gas thus overestimates the flow rate.
In this case, while the fastest gas has flowed to the largest observed radius, the lower velocity gas has not had time to do so.Thus, the assumption of a single radius for the inner wind is probably incorrect, and a flow rate computed from the flow time in each spaxel, t * (i.e., Method 2), is more correct.This intepretation implies that the inferred radius, R II = v rad /t * ,II , is velocity-dependent.Lower-velocity spaxels will have smaller radius, and our single-radius deprojection will underestimate the radial velocity.For both of these reasons, Method 1a produces a significant underestimate for the Episode II wind.
We conclude that Methods 1a and 2 are closer to the correct flow rates for Episode I, while Methods 1b and 2 are closer to the correct flow rates for Episode II.The inferred mass outflow rates are thus in the range 10-13 M yr −1 for Episode I and 151-186 M yr −1 for Episode II.Again, dM/dt II for the ionized gas is comparable to the CO value, 245 M yr −1 .In most nearby compact starbursts, dM/dt in the ionized gas is much smaller than in the molecular gas, though with significant scatter (e.g., Rupke & Veilleux 2013;Fluetsch et al. 2021).However, the canonical density model in these calculations is inferred from optical line ratios; we combine the line ratio densities with the shock model density structure, which raises the inferred masses.In Perrotta et al. (2022, in prep.), we will show that measurements of dM/dt from Mg II and Fe II absorption lines in other compact starbursts from the Tremonti et al. (2007) sample are consistent with mass outflow rates of order 10 2−3 M yr −1 .
Note that we choose the single wind radius for each episode based on two factors.The first is the observed wind morphology (Rupke et al. 2019).We refine the radius to approximately match the average deprojected velocity v rad and the velocity R/t * for the two methods (1a or 1b and 2) that produce the best agreement in dM/dt for each episode.The resulting dM/dt is not sensitive to a 30% change in R wind .

Powering the Nebula and Driving the Wind
We have shown that the line emission from the wind is most consistent with shock models (Section 4.2).We nonetheless consider if the radiative energy from the starburst is capable of powering the Hα luminosity of the nebula.
We first assume the current SFR as given by IR/radio tracers, rather than Hα.A 6.5 Myr instantaneous burst or a 7 Myr continuous burst with star formation rate of 224-300 M yr −1 produces Q(H 0 ) = (1.9−2.6)× 10 55 ionizing photons s −1 (Leitherer et al. 1999).The ratio of Hα photons to ionizing photons is α eff Hα /α, or the ratio of the effective to total recombination coefficients.
For Case B at T = 10 4 K and n e = 100 cm −2 , α eff Hα /α B = α eff Hβ /α B × (j λ λ) Hα /(j λ λ) Hβ = 0.452, where j λ are the line emissivities and we use the emissivity and recombination coefficient tabulations of Hummer & Storey (1987).This yields an Hα luminosity predicted from the radio/IR SFR of L SFR Hα = (2.6 − 3.5) × 10 43 erg s −1 , which is 3 − 4× larger than the observed luminosity of the nebula, L tot Hα (Hα) = 9.4 × 10 42 erg s −1 (Section 3).Photoionizing the entire nebula would require a significant fraction of the Lyman continuum to leak outside of the inner starburst region.We can infer from the predicted Hα luminosity L SFR Hα and the Hα emission attributed to the starburst, L starburst Hα = 1.8 × 10 42 erg s −1 (Section 4.1), that a fraction L starburst Hα /L SFR Hα ∼ 0.05 − 0.07 of the LyC is required to ionize the gas consistent with stellar photoionization (ap1.1).To energize the Hα emisison in the wind (L wind Hα = 8.0 × 10 42 erg s −1 ; Table 3) would require an additional 23-31% of the LyC emitted by the starburst.The total escape fractions for the inner starburst and entire nebula, respectively, would then be of order 90% and 65%, which are uncomfortably high for a galaxy of this mass and ionization parameter (Izotov et al. 2022).Thus, if the nebula were photoionized by the starburst rather than shocks, significant absorption of LyC by dust or a very recent SFR that is much lower would be necessary (Section 4.1).
The shock models also make energetic predictions about the total flux radiated in the shocks to which we can compare the data.The mechanical energy flux that powers the shock is completely radiated by the cooling gas (Dopita & Sutherland 1996).Assuming Hα/Hβ = 2.86, the ratio of the total radiated luminosity L rad to the observed Hα luminosity is given by (Dopita & Sutherland 1995) L shock L Hα = 107 v 0.59 s,100 (1 + 1.32 v −0.13 s,100 where v s,100 is the shock velocity in units of 100 km s −1 .
Here we compare the energy produced in the shocks to the mechanical energy observed in the warm, ionized and molecular phases of the wind episodes.We also look at how the mass, momentum and energy in the wind could arise from the starburst itself.

Episode II ionized wind
The best-fitting model shock velocity is v s = 400 km s −1 , which yields L shock,II = 115L wind Hα,II = 6.6 × 10 44 erg s −1 .If we assume, based on their similar mass and dynamics, that the ionized and molecular components of the fast wind contribute similarly to the energy flow rate dE/dt, then dE/dt HII+H2 II = 1.8 × 10 45 erg s −1 = 2.7L shock,II .In other words, the inferred mechanical energy in the ionized and molecular components of the wind is capable of powering the total radiated luminosity in the shock models.
The mechanical energy produced by the starburst is probably less than both of these numbers.Assuming a continuous starburst of age t * ,II , Solar metallicity, and a Salpeter IMF with stellar masses M = 0.1 − 100 M , Starburst99 predicts a mechanical luminosity in the hot ejecta of dE/dt * = 2.4 × 10 43 erg s −1 (Leitherer et al. 1999).Because of the relatively young age of this starburst, some SNe have not yet released their energy.So older, continuous bursts would contain a greater fraction of the asymptotic value of dE/dt * , which is ∼8 × 10 43 erg s −1 for SFR = 300 M yr −1 .Even if we assume that 100% of this mechanical energy is thermalized into hot ejecta that then drives the outflow, the range of possible dE/dt * is 4% of dE/dt HII+H2 II .An alternative to being driven by the energy in the hot ejecta from the central starburst is that the Makani wind is driven by the mechanical and radiative momentum of the burst.Compact starbursts with strong winds, like Makani, are prime candidates for outflows driven by the radiation pressure from Eddington-limited star formation (Murray et al. 2005;Thompson et al. 2005;Diamond-Stanic et al. 2012).The momentum in the stellar winds and supernova ejecta from the starburst is dp/dt * = 2 dM/dt * dE/dt * (Veilleux et al. 2020).The input mass loss rate from the hot ejecta is dM/dt * ∼ 12 M yr −1 (Leitherer et al. 1999), so that dp/dt * = 1.9 × 10 35 dyn.However, observations of the hot wind eject in M82 suggest that dM/dt * = (0.1 − 0.3)SFR (Strickland & Heckman 2009), depending on the mass-loading of the hot wind by cold clouds within the starburst region.This would imply dM/dt * ∼ 30−90 M yr −1 , and in turn dp/dt * = (3−5)×10 35 dyn.For the calculations below, we take the midpoint of this range, 4 × 10 35 dyn.
We find dp/dt HII+H2 II = 7 × (dp/dt * + dp/dt rad ) for τ eff,IR = 0.A similar momentum "boost" is commonly observed in molecular AGN outflows (Sturm et al. 2011;Cicone et al. 2014;González-Alfonso et al. 2017); in that context, the boost is defined as (dp/dt)/(L/c) .Boosts (dp/dt)/(dp/dt * ) = 1−10 are common in the warm, ionized outflows of 10-100 M yr −1 starbursts (Heckman et al. 2015).As we discuss above, the energy delivered by the starburst hot ejecta may be insufficient to drive the wind.The bulk momentum dp/dt * of the hot wind may be transferred to the cold clouds, but additional momentum can come from work done by the expanding hot bubble during its energy-conserving phase of expansion (Lochhaas et al. 2018;Veilleux et al. 2020).The maximum momentum of a wind powered by a singlephase, radiatively cooling hot wind is also predicted to be several times higher than the estimate of dp/dt * based on M82, at 1.2×10 36 dyne (Lochhaas et al. 2021).When combined with the momentum from radiation pressure, this is only a factor of two below the observed momentum flow rate of the wind.Finally, τ eff,IR 1 (i.e., many photon scatterings) could in principle contribute significantly to driving the wind (Zhang & Davis 2017; though see also Menon et al. 2022).This would be consistent with the dusty, molecular+neutral phase observed in the Episode II flow (Section 4.6).
The gravitational potential of the galaxy will work to decelerate this flow.At a radius of 10 kpc, F grav ∼ 1 × 10 35 dyne for the 4 × 10 9 M of outflowing molecular and ionized gas we observe.Thus, gravity is not significantly decelerating the fast, Episode II wind.Since F grav ∝ M wind /r 2 , at smaller radius the wind may have experienced significant deceleration, but it also may have entrained much less material at earlier times.
Momentum is certainly flowing from the hot to the cool wind.A substantial fraction of the cool, outflowing phase thus arises from clouds that are accelerated by the hot phase.The cool clouds may also acquire mass directly from the hot wind and CGM surrounding the galaxy.The hot wind can transfer mass to the cold phase via radiative cooling (Thompson et al. 2016) and/or turbulent mixing at cloud boundaries (Gronke & Oh 2018, 2020;Fielding et al. 2020;Schneider et al. 2020).
The mass outflow rate in the Episode II wind implies significant entrainment of cool gas, in that the ratio of mass flowing out in the cool, ionized and molecular wind is many times the mass predicted to be in the hot wind originating in the starburst.Using the values of dM/dt * from Strickland & Heckman (2009) that assume some mass-loading of the hot ejecta from stellar winds and SNe within the starburst injection zone (see above), (dM/dt HII+H2 II )/(dM/dt * ) = 5 − 15.The mass outflow rate is also larger than the star formation rate: dM/dt HII+H2 II ∼ 1.4 SFR for SFR = 300 M yr −1 .This value of (dM/dt II )/SFR is consistent with other measurements of star-forming galaxies of similar mass (Heckman et al. 2015), which range from 0.3 to 3. Though the (dM/dt II )/(dM/dt * ) ratio is conceptually closer to representing the mass loading of the hot wind with cool clouds, the ratio (dM/dt II )/SFR is also called the "mass loading factor."The degree of cold cloud entrainment in analytic two-phase models then points to large initial cloud sizes M cl 10 5−6 M in the Makani wind, such that the cooling time is short compared to the mixing time (Fielding & Bryan 2022).It also implies that the large observed cool cloud "mass loading" of (dM/dt II )/SF R ∼ 1.4 is not, in fact, too large to cause the wind to stall, as suggested by some energydriven models (Thompson et al. 2016;Fielding & Bryan 2022).This is probably consistent with the need for significant radiation-pressure driving of the cool wind.
This mixing of mass and energy between the cool and hot phases could also be an alternative power source for the ionized nebular emission.In the single-phase analytic model of Thompson et al. (2016), the hot wind cools radiatively as it expands into the galaxy, and then again as it shocks with the surrounding CGM.By contrast, in the explicitly multiphase analytic model of Fielding & Bryan (2022), the cool phase exists even at the base of the wind.The multiphase wind carries millions of n 10 cm −3 clouds that are area-filling (though not volume-filling).The cooling gas in the cloud boundary layer between the hot and cool phases, with T = 0.2 − 1 × 10 5 K, then produces line emission in the wind.
Increasing n e in the ionized gas, so that all of the Hα emission arises in the compressed, post-shock region (Section 4.3), would lower the mass, momenta, and energy in the ionized gas by a factor ∼50.This would reduce the observed momentum boost and mass loading factor.However, the fact would remain that the mass, momentum, and energy in the molecular gas are still at the level shown in Table 3 for Episode II, so that the total flow rates would only decrease by a factor of order 2. Furthermore, the large flow rates we observe are consistent with those measured in the neutral and ionized phase of similar galaxies via other probes (Perrotta et al. 2022, in prep.).

Episode I ionized wind
The Episode I ionized wind, with v s = 200 km s −1 , has an estimated ionized mass M I ∼ 4M II .This mass may be distributed over a much larger volume, however: [(R I /R II ) 3 − 1] ∼ 15.The outflow rate in Episode I is also more than 10× smaller.As discussed above, high-velocity gas from Episode I may have carried significant amounts of mass beyond r wind = 50 kpc, further into the CGM.Alternatively, the wind may simply have slowed down in the extended halo/CGM of the galaxy, reducing dM/dt while still entraining more gas at larger radius (Lochhaas et al. 2018).Consistent with this possibility is the observation that the outflow velocity in the parent sample of Makani decreases with increasing light-weighted age (Davis et al. 2022, submitted).In Episode II, some or all of the cool wind arose from the host galaxy, originating in cool clouds or condensed from the hot wind.However, much more of the mass entrainment at larger radii may be due to radiative cooling and/or turbulent mixing from the hot wind and/or CGM (Thompson et al. 2016;Gronke & Oh 2018;Fielding et al. 2020).
The radiated energy in the Episode I shocks is L shock,I = 73L Hα,I = 1.6 × 10 44 erg s −1 .For Episode I, we include only the ionized gas, since the CO gas is contained in the Episode II footprint.Thus, dE/dt HII I = 0.003L shock,I , a significant discrepancy.We conclude that much of the energy in the shocks beyond r ∼ 20 kpc must come from the mechanical energy of a different phase of the outflow.This could be the hot phase traced by higher ionization states like O VI, as observed in the CGM in absorption at large impact parameters (e.g., Tumlinson et al. 2011).Shocks may in fact arise naturally in the interaction between a hot wind and surrounding CGM (Thompson et al. 2016).

Episode II dusty wind
Observations of molecular gas and Mg II 2796, 2803 Å emission to radii 20-25 kpc point to the presence of significant H I and H 2 in the inner, Episode II wind.The presence of these two tracers also indicates significant amounts of dust out to similar radii.However, these dusty gas phases appear only at the inner edge of the Episode I wind (Rupke et al. 2019).
Two new lines of evidence confirm this picture.First, Na I D is detected in two apertures out to r = 10−15 kpc (Figure 6).The presence of Na I D has long been closely associated with dust extinction (see, e.g., Rupke et al. 2021, and references therein), due to dust being required to keep Na in the neutral phase.Second, we detect substantial gas extinction, E(B − V ) ∼ 0.5, out to r = 20 − 25 kpc (Figure 4).
Such dusty, neutral outflows on ∼10 kpc scales are commonly observed in compact starbursts in nearby mergers (Rupke & Veilleux 2013).Frequently, high spatial resolution observations show filamentary dust structures associated with these outflows.Future highresolution observations of Makani with the James Webb Space Telescope may resolve these dust structures, as well as more directly determine the radial extent of dust in the wind.

CONCLUSION
The giant Makani galactic wind is driving ionized, neutral, and molecular gas and dust out of a compact starburst into the CGM.In the preceding sections we have presented a spatially-resolved analysis of the physical state and mass, momentum, and energy in the ionized phase of the wind using rest-frame optical spectra.These ionization and dynamical measurements are key to unlocking the driving force of the Makani wind, as well as its impact on and interaction with the surrounding CGM.
The fast, inner wind extending to R II = 20 kpc-Episode II-is powered by a young starburst of age 7 Myr (Rupke et al. 2019) and star formation rate 224-300 M yr −1 (Petter et al. 2020).Using line ratios and the luminosity of recombination lines, we conclude that the ionized gas in the inner wind is energized by fast, v s = 400 km s −1 shocks moving through a low-density n H ∼ 10 cm −3 medium (Allen et al. 2008), consistent with earlier results (Perrotta et al. 2021).The velocity dispersion of the gas is equal to this model shock velocity: σ II = 400 km s −1 .The shock compresses this gas to high density, yielding n e 10 3 cm −3 in the post-shock region.The hard radiation field of the shock produces high-ionization line emission commonly seen in AGN ([Ne V] 3426 Å), as well as high gas temperatures T > 2 × 10 4 K, as traced by [O III] 4363 Å and [N II] 5755 Å.
Molecular and neutral gas was previously detected throughout the Episode II wind (Rupke et al. 2019).We find new evidence for neutral, dusty gas in the wind in the form of Na I D absorption in the host galaxy and emission out to 10-15 kpc; and extinction in the ionized gas, as traced by the Balmer decrement, to 20-25 kpc.
A much older and larger wind, extending from R I = 20 − 50 kpc-Episode I-was powered by a star formation event 400 Myr ago (Rupke et al. 2019).Presently, only ionized gas is detected in the Episode I wind, and it is energized by 200 km s −1 shocks in n H 2 − 3 cm −3 gas.Again, this is consistent with the observed velocity dispersion of this outflow, σ I = 200 km s −1 (Rupke et al. 2019).
Using measurements of [O II] 3726, 3729 Å and Hα throughout the wind, we model the radial dependence of their ratio.We apply this model to the fully spatiallyresolved [O II] map of the nebula to compute L Hα across the full nebula.We then combine the shock structure predicted by models (Dopita & Sutherland 1996;Allen et al. 2008) with the density from line ratios to measure the mass in each wind episode.Finally, we employ three different methods for computing dM/dt, p, dp/dt, E, and dE/dt in the ionized wind from the 2D mass and velocity maps of the nebula.
Our preferred model of the Episode II nebula is a massive, powerful, ionized wind that is comparable to that of the molecular phase, with M HII II ∼ M H2 II = (1 − 2) × 10 9 M and dM/dt HII II ∼ dM/dt H2 II = 170−250 M yr −1 .Together, the ionized and molecular phases are depleting gas from the galaxy at 1.4× the star formation rate.These phases carry as much energy as is predicted to be fully radiated by the shocked gas (Dopita & Sutherland 1995).They need significantly more momentum than is initially contained in the hot ejecta from the recent starburst, with a required boost of ∼7.We suggest radiation pressure from an Eddington-limited starburst is a likely culprit (Diamond-Stanic et al. 2012;Thompson et al. 2015), on top of the momentum produced by the hot wind (Thompson et al. 2016;Lochhaas et al. 2021;Fielding & Bryan 2022).
The slower, Episode I wind has a much lower outflow rate of dM/dt HII I ∼ 10 M yr −1 .However, this ionized wind has moved into the galaxy's CGM and is depositing gas in a much larger volume.It contains even more ionized gas; in our density model, M HII I = 5 × 10 9 M .This increase in mass as the wind slows but moves to larger radius may be due to loading from radiative cooling of the hot wind or CGM and/or mixing of the hot gas into the cool phase through mixing layers (Lochhaas et al. 2018;Gronke & Oh 2020;Fielding et al. 2020).The gas may also be shocking on the surrounding CGM (Thompson et al. 2016).The energy contained in this flow is quite small, however, compared to the flux predicted to be radiated in the shocked gas, with dE/dt I = 0.003L rad .Thus, much of the energy and momentum in Episode I may be carried in a hotter phase.Alternatively, significant energy from the fastestmoving gas in the earlier wind may simply have escaped to large radius.
The clearest way to make progress on the Makani wind and its connection to the CGM will be to try to detect hot gas in the extended nebula and measure its mass and energy content.This could be in the form of deep Xray observations or UV observations of ionization states like O VI.On the other side of the energetic spectrum, the dusty, neutral phase is detected to 20-25 kpc, but its physical conditions are mostly unknown.Deep mid-infrared imaging and spectroscopy with the JWST would produce much stronger constraints on its presence and properties in both wind episodes.Finally, we could reduce uncertainties in the present work through integral-field maps of the recombination lines in Makani with wide-field, red-sensitive instruments like the Multi Unit Spectroscopic Explorer (MUSE) or the Keck Cosmic Reionization Mapper (KCRM).
Figure 1.The Makani [O II] nebula as observed by KCWI (Rupke et al. 2019), with ESI long slits overlaid.The peak [OII] KCWI spaxel is labeled with a cross.The extracted apertures are shown as colored boxes.North is up and East is left.
2: Aperture extraction size.Column 3: Radius of the aperture center, plus or minus the full range of radii represented by the boundaries of the aperture.Column 4: Observed [O II] 3726, 3729 Å flux in the aperture, with 1σ errors.

Figure 2 .
Figure 2. Observed spectrum of each aperture and model fits near strong emission lines.In each panel, the lines plotted are in the upper left-hand corner.The black line is the data, in units of 10 −16 erg s −1 cm −2 Å−1 ; the red line is the total model fit (lines+continuum); and dashed lines are the emission-line model.Different velocity components are shown as different colors.The best-fit redshift of each component is shown in the [O II] panel, as is an inset showing the velocity profile of [O II] 3729 Å.Because the velocity of [O I] is allowed to vary separately in some apertures, its best-fit redshifts are shown in the corresponding panel.The data are smoothed with a 3-pixel boxcar.The data and model residuals after subtraction of the emission-line model (black and red lines) and the 1σ error in the data (teal line) are shown in the bottom third of each panel.(The complete figure set (9 images) is available.) Figure 3. Fit to He I 5876 Å emission and Na I D 5890, 5896 Å absorption and emission in the central spectrum, ap1.Vertical dashed lines show the expected wavelengths based on our fit to the stellar continuum (magenta) and [O I] emission lines (teal) in ap1.Velocities with respect to these reference redshifts (zstar for Na I D and zem for He I) are shown in the upper right.
Figure 4. (Top) Radial dependence of observed [O II] and Hα surface brightness in each spectrum.The line shows the azimuthally-averaged radial surface brightness profile of [O II] measured with KCWI (Rupke et al. 2019).Each point is centered on the radius of the aperture center; horizontal grey lines are the radii spanned by each extraction aperture; and vertical grey lines are 1σ errors.Points in the same aperture are connected by vertical bars, and points at radius > 0 kpc are randomly offset slightly in radius to prevent visual overlap.(Upper middle) Color excess E(B − V ) vs. radius.E(B − V ) is computed from the Balmer decrement.Vertical error bars are 1σ.Point colors correpond to the aperture colors from Figure 1.Reddening is observed out to 20 kpc.(Lower middle) Variation of observed and intrinsic [O II]/Hα] flux ratio with distance from the host galaxy.The final point is a lower limit due to the large uncertainty in E(B − V ).Both ratios increase steadily from the host to the outskirts of the nebula.(Bottom) Correction from observed [O II] to intrinsic Hα vs. radius.The lines show a linear fit to the data and (fit±RMS).
[N II]/Hα vs. [O III]/Hβ diagram.The two apertures with r ∼ 10 kpc (ap2 and ap3) combine information from two velocity components, one broad and one narrow.Ap2 has line ratios similar to those from total fluxes in ap1, with higher low-ionization line ratios ([O I]/Hα) and lower ionization parameter [O III]/[O II].Apertures ap3, ap4, and ap5 have similar [O I]/Hα and [O III]/[O II] as ap2, but with lower highionization line ratios ([O III]/Hβ), larger [S II]/Hα, and slightly smaller [N II]/Hα.Two other apertures with significant detections across several lines (ap6 and ap8) lie in roughly the same locations as the other apertures, but these apertures have low S/N in the weakest lines so their positioning is more uncertain.The lines other than [O II] and Hα in ap7 and ap9 are not strong enough to accurately place them in these diagrams, though ap9 does appear in several with large error bars.

Figure 5 .
Figure 5. Excitation plots for the Makani nebula.The dashed and dotted lines divide nuclear classification regions of star-forming (H II), AGN, LINER, and composite (C) using canonical dividing lines (Kauffmann et al. 2003; Kewley et al. 2006; Cid Fernandes et al. 2010).The symbol size reflects projected distance from the host galaxy, with smaller symbols at larger radius.Filled circles are from total fluxes.The bowtie is the central aperture broad component (ap1.2), and the open circle is the corresponding narrow component (ap1.1).
Properties of the Star-Forming Host

[
Figure 6.Radial variation of neutral gas tracers Mg II 2796, 2803 Å and Na I D in Makani.The blue line shows the azimuthally-averaged radial dependence of Mg II surface brightness from KCWI data (Rupke et al. 2019), with half-light radius r Mg II e = 9 kpc and maximum extent 20-25 kpc.The slower decline of Na I D emission, based on the current ESI data, suggests emission-line infilling in the innermost region, reducing its flux by ∼0.6 dex (Section 3.2).The lack of neutral/molecular gas (based on Mg II and CO from Rupke et al. 2019 and Na I from the current work) and dust (based on E(B −V ) from the current work; Figure 4) beyond 20-25 kpc is indicated with the vertical, orange dotted line.The full extent of the ionized gas nebula from KCWI observations is shown as the solid green line, with r [O II] e = 17 kpc and r [O II] max = 50 kpc.

4. 3 .
Shock structure and wind density We can combine the observed line flux ratios [S II]6716/6731 and [O II]3729/3726 with the Allen et al. (

Table 2 .
Fluxes and Luminosities yields an escape fraction in the Mg II 2796 Å line of 25%.Using the line fluxes from ap1.1 in the current data and the method of Xu et al. (2022), we can make a second, independent measure of f esc (2796).We apply the highest-metallicity model relating the intrinsic values of [O III]5007/[O II] and Mg II2796/[O III]5007

Table 3 .
Ionized Wind Properties