The ALMA View of Positive Black Hole Feedback in the Dwarf Galaxy Henize 2–10 Hansung B. Gim1 and Amy E. Reines1,2 1 Department of Physics, Montana State University, Bozeman, MT 59717, USA; hansung.gim@montana.edu 2 eXtreme Gravity Institute, Department of Physics, Montana State University, Bozeman, MT 59717, USA Received 2023 May 25; revised 2023 December 13; accepted 2024 January 3; published 2024 March 4 Abstract Henize 2–10 is a dwarf starburst galaxy hosting a ∼106Me black hole (BH) that is driving an ionized outflow and triggering star formation within the central ∼100 pc of the galaxy. Here, we present Atacama Large Millimeter/ submillimeter Array continuum observations from 99 to 340 GHz, as well as spectral line observations of the molecules CO (1–0, 3–2), HCN (1–0, 3–2), and HCO+ (1–0, 3–2), with a focus on the BH and its vicinity. Incorporating centimeter-wave radio measurements from the literature, we show that the spectral energy distribution of the BH is dominated by synchrotron emission from 1.4 to 340 GHz, with a spectral index of α≈− 0.5. We analyze the spectral line data and identify an elongated molecular gas structure around the BH with a velocity distinct from the surrounding regions. The physical extent of this molecular gas structure is ≈130 pc× 30 pc and the molecular gas mass is ∼106Me. Despite an abundance of molecular gas in this general region, the position of the BH is significantly offset from the peak intensity, which may explain why the BH is radiating at a very low Eddington ratio. Our analysis of the spatially resolved line ratio between CO J= 3–2 and J= 1–0 implies that the CO gas in the vicinity of the BH is highly excited, particularly at the interface between the BH outflow and the regions of triggered star formation. This suggests that the cold molecular gas is being shocked by the bipolar outflow from the BH, supporting the case for positive BH feedback. Unified Astronomy Thesaurus concepts: Interstellar medium (847); Dwarf galaxies (416); Active galaxies (17) 1. Introduction Galaxy evolution is governed by the balance between star formation and quenching (Somerville & Davé 2015), and active galactic nuclei (AGNs) are thought to play a crucial role in the extinguishing of star formation in massive galaxies. Prior to the last decade, massive black holes (BHs) that sometimes accrete and shine as AGNs were almost exclusively found in giant galaxies (e.g., Kormendy & Ho 2013). However, we now know that AGNs in dwarf galaxies are much more common than previously thought and the BHs have typical masses of MBH< 106Me (for reviews, see Greene et al. 2020; Reines 2022). Searching for and studying these systems is important for our understanding of the role of AGN fueling and feedback in the low-mass regime and may also provide clues about the origin of the first BH “seeds” in the early Universe (e.g., Volonteri et al. 2021). Henize 2–10 (He 2–10) is a nearby (d= 9Mpc) dwarf starburst galaxy (e.g., Nguyen et al. 2014) that is particularly remarkable for displaying evidence of positive AGN feedback (i.e., enhancing star formation; Schutte & Reines 2022; see also Figure 1 in this paper). Multiwavelength studies have shown that this galaxy hosts a massive BH in its center (Reines et al. 2011) through coincident peaks at radio and X-ray wavelengths (Reines et al. 2016) detected with the Very Large Array (VLA) and the Chandra X-ray Observatory (Chandra), respectively. The BH is also detected as a strong nonthermal source by high-angular-resolution observations using the Long Baseline Array (LBA; Reines & Deller 2012). The mass of the central BH is estimated to be MBH∼ 106Me, based on stellar velocity measurements (Riffel 2020) and the scaling relation between the total stellar mass and BH mass (Reines & Volonteri 2015). Combined with the X-ray luminosity, this indicates the BH is accreting at a very low Eddington ratio (Reines et al. 2016). A recent study using Hubble Space Telescope (HST) spectroscopy revealed that the central BH is driving an outflow of ionized gas, which is triggering the formation of star clusters located ∼70 pc (∼1 5) from the BH (Schutte & Reines 2022). This is the first case of positive BH feedback found in a dwarf galaxy and may represent a low-power analog of radio-loud AGNs experiencing “jet-mode” feedback (Schutte & Reines 2022). Molecular gas and dust can both feed and obscure AGNs. While the cold gas distributions of He 2–10 have been extensively studied for the galaxy as a whole and with a focus on the star-forming regions, very little is known about the molecular gas in the immediate vicinity of the BH. Here we present and analyze high-resolution Atacama Large Milli- meter/submillimeter Array (ALMA)3 observations of He 2–10, with a specific focus on the central massive BH and its immediate surroundings (∼130 pc). A summary of previous studies of the molecular gas in He 2–10 is given below. Baas et al. (1994) first observed the three CO transitions (J= 1–0, J= 2–1, and J= 3–2) of He 2–10, showing that hot and cold molecular gas is distributed throughout the galaxy. The first interferometric observations of CO J= 1–0 with a synthesized beam size of 6 5× 5 5 (284 pc× 240 pc) revealed an elongated morphology, indicating that He 2–10 could be a post-merger system and/or experiencing tidal interactions (Kobulnicky et al. 1995). Meier et al. (2001) observed the CO J= 3–2 line to investigate the interstellar medium conditions in the central 13″ (567 pc) and found a gas temperature of TK∼ 5–10K and a number density of n(H2)> 103.5 cm−3. Higher transitions of CO lines, J= 3–2, 4–3, 6–5, and 7–6, were observed by Bayet The Astrophysical Journal, 963:103 (14pp), 2024 March 10 https://doi.org/10.3847/1538-4357/ad1b62 © 2024. The Author(s). Published by the American Astronomical Society. Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. 3 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 1 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X mailto:hansung.gim@montana.edu http://astrothesaurus.org/uat/847 http://astrothesaurus.org/uat/416 http://astrothesaurus.org/uat/17 https://doi.org/10.3847/1538-4357/ad1b62 https://crossmark.crossref.org/dialog/?doi=10.3847/1538-4357/ad1b62&domain=pdf&date_stamp=2024-03-04 https://crossmark.crossref.org/dialog/?doi=10.3847/1538-4357/ad1b62&domain=pdf&date_stamp=2024-03-04 http://creativecommons.org/licenses/by/4.0/ et al. (2004). They found that the central region hosts warm and dense molecular gas with TK≈ 50–100K and n(H2)> 104 cm−3, which is an order of magnitude higher than the temperature calculated with lower transitions of CO lines (J= 1–0, 2–1, and 3–2) by Meier et al. (2001). Beck et al. (2018) observed the CO J= 3–2 line with a synthesized beam size of 0 3 (13 pc) using ALMA and found several clumpy filaments separated spatially and kinematically. Their kinematic analysis led to the conclusion that the molecular filaments feed the star formation in the star clusters. Imara & Faesi (2019) traced the CO J= 1–0 line with a synthesized beam size of 0 67× 0 58 (29 pc× 25 pc) using ALMA and showed that He 2–10 has a more complex morphology with clumps, ring-shaped structures, and high-density peak emissions. Higher-critical-density molecular gas tracers have also been observed in He 2–10. Imanishi et al. (2007) observed the HCN J= 1–0 and HCO+ J= 1–0 lines with a beam size of 10 8× 5 5 (471 pc× 240 pc) and reported the ratio of HCN (1–0)/HCO+ (1–0) < 0.6, which is low compared to luminous infrared galaxies. The strong HCN J= 1–0 line was also recovered by the single-dish observations of Santangelo et al. (2009). Johnson et al. (2018) used ALMA observations tracing HCN J= 1–0, HCO+ J= 1–0, HNC J= 1–0, and CCH J= 1–0 with a synthesized beam size of 1 7× 1 6 (74 pc× 70 pc) and showed that HCO+ J= 1–0 correlates with the thermal radio continuum in the intense central star- forming regions. Here, we present the first study of the molecular gas in the immediate vicinity of the central BH using high-resolution observations with ALMA. We describe the observations and data reduction in Section 2, the analysis and results in Section 3, and our conclusions in Section 4. 2. ALMA Observations and Data Reduction He 2–10 has been observed by ALMA in different bands for six different projects. The specifics of the ALMA observations are summarized in Table 1, encompassing essential information such as project codes, bands, frequency ranges, spectral lines, total integration time, bandpass, phase, and amplitude calibrators. The names of the PIs are presented in the caption located below the table. We downloaded all the data listed in Table 1 from the ALMA archive website.4 The calibrated data were processed using scriptForPI.py with the appropriate version of the Common Astronomy Software Application (CASA; Bean et al. 2022), as outlined in the README files. We also flagged some bad antennas, which improved the sensitivity. 2.1. Spectral Line Image Cubes Spectral line image cubes were generated using CASA 6.5.1 with the continuum-subtracted spectral line visibilities. The continuum subtraction was performed using the CASA task Figure 1. The four-color composite image of He 2–10 is generated from F330W (B), F606W (V ), F658N (Hα), and F814W (I) from HST. The white box shows the central 6″ × 4″ region, and the inset shows the Hα+continuum image presenting the outflow from the BH. (Credit to NASA, ESA, Zachary Schutte, and Amy Reines.) 4 The ALMA archival website is at https://almascience.nrao.edu/aq/. 2 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines https://almascience.nrao.edu/aq/ uvcontsub, which uses line-free channels and a linear fit of order 1 to estimate the continuum. The continuum-subtracted visibilities were then imaged using the CASA task tclean, with the Högbom algorithm for deconvolution of the point-spread function, briggsbwtaper, and perchanweightdensity. These options provide a flat rms noise and beam size distribution over the channels, and a beam size comparable to that of the continuum image. The velocity frame is set to BARY, which normalizes the velocity to the barycentric system. The cleaning process was repeated until the maximum of the residual map reached 2σ, where σ is the rms noise of the image cube per channel. The line-free channels for the continuum subtraction, image size, cell size, velocity spacing, robust value for the briggsbwtaper, and rest frequency are summarized in Table 2. The resulting image cubes were smoothed to the largest beam size within the selected velocity range including lines with the CASA task imsmooth. The resulting beam sizes and rms noises of the molecular line image cubes are summarized in Table 2. 2.2. Continuum Imaging We created four different continuum images: low-frequency Band 3 (low-Band 3), high-frequency Band 3 (high-Band 3), Band 6, and Band 7 continuum images. The Band 3 continuum data are from ALMA project codes 2011.0.00348.S, 2015.1.01569.S, 2016.1.00027.S, and 2019.1.01641.S, which were split into low- and high-frequency data at 110 GHz. This frequency was chosen to maximize the angular resolution of the high-Band 3 continuum image. The Band 6 and 7 continuum data are from ALMA project codes 2012.1.00413.S and 2016.1.00492.S, respectively. Table 1 ALMA Observations of He 2–10 (1) (2) (3) (4) (5) (6) (7) (8) Project Code Band Frequency Lines Integration Bandpassa Complex Gainb Flux Densityc (GHz) 2011.0.00348.S (1) 3 86.95–102.47 HCN (1–0) 204m 37s J053851-440507 J082601-223027 Titan HCO+ (1–0) Ceres 2012.1.00413.S (2) 6 248.79–267.73 HCN (3–2) 286m 47s J0750+1231 J0747-3310 J1037-295 HCO+ (3–2) J0538-4405 J0826-2230 J0538-4405 J1058+0133 2015.1.01569.S (3) 3 111.99–115.18 CO (1–0) 187m 50s J1037-2934 J0826-2230 J1107-4449 2016.1.00027.S (3) 3 99.72–115.17 J0538-4405 J0846-2607 2016.1.00492.S (4) 7 332.06–348.00 CO (3–2) 50m 29s J0538-4405 J0846-2607 J0538-4405 2019.1.01641.S (5) 3 96.91–112.51 L 103m 43s J0725-0054 J0846-2607 J0725-0054 J1037-2934 J1037-2934 Notes. PI names: (1) Kelsey Johnson; (2) Amy Reines; (3) Nia Imara; (4) Sara Beck; and (5) Zhiyu Zhang. a Bandpass calibrators. b Complex gain calibrators. c Flux density scale calibrators. Table 2 Parameters for Imaging and Image Statistics (1) (2) (3) (4) (5) (6) (7) (8) (9) Line Line-free Image Size Cell Size Velocity Spacing Robust Rest Frequency Syntdesized rms Channels Beam Size (pixel × pixel) (arcsec pixel−1) (km s−1) (GHz) (μJy beam−1) HCN (1–0) 301–750 1150–3539 512 × 512 0.25 10.0 0.5 88.6316 1 78 × 1 68 350a HCO+ (1–0) 300–3300 3500–3539 512 × 512 0.25 10.0 0.5 89.1885 1 84 × 1 68 350a HCN (3–2) 200–1600 2000–3640 1024 × 1024 0.025 20.0 1.5 265.886 0 19 × 0 16 <470a HCO+ (3–2) 200–1600 2500–3640 1024 × 1024 0.025 7.5 1.5 267.558 0 18 × 0 16 325a CO (1–0) 100–1700 2150–3740 1024 × 1024 0.05 1.693 0.0 115.208201 0 45 × 0 28 1030a CO (3–2) 10–90 161–230 1024 × 1024 0.05 1.693 0.5 345.79599 0 31 × 0 28 1000a 99 GHz L 1024 × 1024 0.05 L 0.0 L 0 55 × 0 32 13.0 113 GHz L 1024 × 1024 0.05 L 0.0 L 0 45 × 0 27 18.0 251 GHz L 1024 × 1024 0.03 L 1.5 L 0 18 × 0 15 12.4 340 GHz L 1024 × 1024 0.04 L 0.5 L 0 32 × 0 27 27.7 Note. a The unit is (μJy beam−1 (channel width)−1) for spectral image cubes. 3 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines The continuum images were created by combining the visibilities from the continuum spectral windows (SPWs) with wide channel widths of 15.625 MHz and the line-free channels from the SPWs containing spectral lines with narrow channel widths (as explained in Section 2.1). The latter were smoothed to a channel width of 15.625 MHz in the SPWs for the continuum using the CASA task mstransform. The continuum images were generated using the CASA 6.5.1 task tclean with the Briggs weighting function and the Högbom deconvolution algorithm. The standard multifre- quency synthesis (MFS) was applied for imaging the continuum data, except for low-Band 3. For the low-Band 3 data, we applied the multiterm MFS (Rau & Cornwell 2011), as its fractional bandwidth is about 23.7% (with a bandwidth of 23.4 GHz, or 87.1–110.5 GHz, relative to the effective frequency of 98.8 GHz). We used nterms= 3, which yields a more accurate intensity map by calculating the mean intensity, spectral index, and curvature over the frequency. The cleaning was stopped when it reached 2σ, where σ is the rms noise of the image. We summarize the parameters used in the CASA task tclean and the resulting synthesized beam sizes and rms noise in Table 2. The synthesized beam sizes and sensitivities of the continuum images in this paper are slightly different from those obtained in Costa et al. (2021) for 113 and 340 GHz, while the characteristics of the 251 GHz image are the same. Our 113 GHz continuum image was obtained by concatenating the data of 2015.1.01569.S, 2016.1.00027.S, and 2019.1.01641.S, which is different from Costa et al. (2021) who used 2016.1.00027.S for the 113 GHz continuum image. The synthesized beam size and rms noise were 0 38× 0 21 and 17 μJy beam−1, respectively (Costa et al. 2021), where the beam size is 18% better than our 113 GHz continuum image. This difference is attributed to adding two other observations to our data, 2015.1.01569.S and 2019.1.01641.S, since the robust value used in the CASA task tclean was the same as ours. The project codes 2015.1.01569.S and 2019.1.01641.S have maximum baselines of 1.4 km and 181.9 m, respectively, while the project 2016.1.00027.S has a maximum baseline of 2.6 km. Adding these two data sets with shorter baselines degraded the angular resolution when the data of 2016.1.00027.S are used solely. While the synthesized beam sizes are almost the same for the 340 GHz images, the rms noises are 27.7 and 52.0 μJy beam−1 for our image and that of Costa et al. (2021). The reason for this difference is unclear. 2.3. Astrometry It is important to facilitate a comparison of the multi- wavelength observations of He 2–10 taken with different telescopes. The location of the BH in He 2–10 was originally determined in the reference frame of the Two Micron All Sky Survey (2MASS), due to its absolute positional accuracy of 0 1. Images from the HST, VLA, and Chandra were shifted to match the locations with the astrometric frame of 2MASS (Reines et al. 2011). Here, we identified the offsets in our ALMA observations with respect to the 2MASS reference frame and applied them to the ALMA observations. We used the highest-angular-resolution ALMA continuum image at 251 GHz (with a synthesized beam size of 0 18× 0 16) and the 1.4 GHz continuum image observed at the LBA with a synthesized beam size of 0 11× 0 03 to estimate the offset. The astrometry of the LBA 1.4 GHz continuum image was corrected to that of 2MASS (Reines & Deller 2012). The shift was determined by matching the local maximum in the 251GHz continuum image to the center of the BH in the 1.4 GHz continuum image. We estimated a shift of ΔR.A.= 0 0735 west and Δdecl.= 0 095 north and applied this shift to all ALMA spectral line image cubes and continuum images. 3. Analysis and Results 3.1. Continuum Emission and Radio Spectral Energy Distribution of the BH Figure 2 presents the submillimeter continuum maps at 99 (first row), 113 (second row), 251 (third row), and 340 GHz (fourth row), respectively. The white contours show the 99 GHz continuum at 100, 150, 180, 200, and 250 μJy beam−1. The local maxima of the contours in the center Figure 2. Continuum images at 99 GHz (A), 113 GHz (B), 251 GHz (C), and 340 GHz (D). The central BH is marked with the black ellipse and red arrow in the 99 GHz continuum image (A). The white contours represent the 99 GHz continuum at the levels of 50, 100, 150, and 200 μJy beam−1. The beam size in each image is shown as the gray ellipse at the bottom right. 4 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines indicate the position of the BH, which coincides with the positions determined by the VLA and LBA (Reines et al. 2011; Reines & Deller 2012). The continuum images show two prominent blobs and one central source. The two strong emission blobs correspond to the starburst regions, and the central source indicates the massive BH (Reines et al. 2011; Reines & Deller 2012). The position of the BH is defined in the continuum images as the local maxima in the central region between the two prominent starburst regions, as indicated in the top panel of Figure 2. The two starburst regions are clearly visible in all the continuum maps; however, the central BH is not seen in the 340 GHz continuum map. In the 340 GHz continuum image, the BH is located between two small continuum sources located 0 217 (10 pc) and 0 526 (23 pc) away. An earlier study by Allen et al. (1976) showed that He 2–10 as a whole is dominated by nonthermal emission. However, high-angular-resolution VLA and ALMA observations have shown that the central few hundred parsecs are dominated by free–free emission from bursty star formation (Johnson & Kobulnicky 2003; Johnson et al. 2018; Costa et al. 2021). The central BH, on the other hand, is governed by nonthermal synchrotron emission (Johnson & Kobulnicky 2003; Reines et al. 2011; Reines & Deller 2012; Costa et al. 2021). Although the VLA + ALMA spectral energy distribution (SED) analysis for the BH has been done by Costa et al. (2021), we revisited this analysis with revised flux density measurements of the submillimeter continuum and the inclusion of the 1.4 GHz flux density observed at the LBA (Reines & Deller 2012). The radio SED of the BH from 1.4 to 340 GHz is shown in Figure 3. In addition to our new ALMA measurements in the submillimeter regime (see below), we include measurements from the literature, including 1.4 GHz data from the LBA (Reines & Deller 2012), 5 and 8 GHz data from the VLA (Reines et al. 2011), as well as 15, 22, and 33 GHz data from the VLA (Costa et al. 2021). Here, the flux density at 33 GHz is not included in this analysis, because the data suffered from decorrelation due to weather conditions and the low decl. of He 2–10 at this high frequency (Costa et al. 2021). Therefore, the flux density at 33 GHz is presented with a lower limit (upper arrow) in Figure 3. The flux densities of the BH in the submillimeter regime were estimated with our reimaged ALMA continuum data. First, to help mitigate the impact of the different beam sizes, we convolved all the continuum images to have the same synthesized beam size as the lowest one (0 55× 0 32 for the 99 GHz image) using the CASA task imsmooth. While the inherent limitations of interferometry may still result in different spatial sensitivities among the observations, as well as the omission of large-scale emission, this would not be a major concern if the emission associated with the BH is point- like. However, if the emission is marginally resolved, we may be missing flux density that could, for example, lead to the low measurement at 251 GHz seen in Figure 3. We performed aperture photometry on the submillimeter continuum using an elliptical aperture equivalent to one beam size centered on the BH position (i.e., 0 55× 0 32 with a position angle of −64°.0, marked with the blue ellipse in Figure 3). The arrow points to the BH, which is shown as the black ellipse. The measured flux densities of the BH are 140± 14, 134± 21, 29± 46, and 57± 48 μJy at 99, 113, 251, and 340 GHz, respectively. Here, it is worth noting that the uncertainty of the flux density may be underestimated, because we measured it as a mean of the rms noise within an aperture. The rms noise maps were produced with the AIPS task RMSD (Greisen 2003), which calculates the rms of signals after rejecting signals greater than 3σ within 150 pixels of each pixel. Our measurements are consistent with those of Costa et al. (2021) at 113 GHz, but not at 251 and 340 GHz. Costa et al. (2021) reported flux densities of 70± 30 and 130± 30 μJy at 258 and 340 GHz, respectively, while our measurements are 29± 46 and 57± 48 μJy at these frequencies. At 251 GHz, our photometry shows 2.4 times less flux density than that of Costa et al. (2021). Costa et al. (2021) likely measured the flux density of a point source shown within the contour in the center of Figure 2(C). However, this point source is smeared after convolution to the beam size of 0 55× 0 32 in our measurement, as shown in Figure 3(C). Since we applied Figure 3. Continuum images at 99 GHz (A), 113 GHz (B), 251 GHz (C), 340 GHz (D), and the SED analysis (E). Continuum images were convolved to have the same synthesized beam size of 0 55 × 0 32. The central BH is marked with the small black ellipse indicated by the arrow. The blue ellipse is the aperture of our flux density measurements. The beam size in each image is shown as the gray ellipse at the bottom right. The sources A and B in the 340 GHz continuum image (D) are adjacent sources of the BH. In (E), the flux densities are depicted as the black data points, with a black arrow indicating the lower limit at 33 GHz. The best-fit model is represented by a red line, and the 1σ range of the best-fit model is delineated by the red shaded region, providing a measure of the model’s uncertainty. 5 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines aperture photometry on the convolved image, we measured a smaller flux density. At 340 GHz, our measurement is 2.3 times fainter than the flux density obtained by Costa et al. (2021). This difference is due to the difference in measurements. There are two continuum sources near the BH position, with distances of 0 22 (source A) and 0 53 (source B) from the BH, as shown in Figure 3(D). Their flux densities were measured as 125± 80 and 342± 130 μJy for the sources A and B, respectively, using PyBDSF (Mohan & Rafferty 2015). Costa et al. (2021) likely reported the flux density of source A as that of the BH. Our aperture for the BH does not include source A and source B entirely, but includes a part of source A, as seen in Figure 3(D). Therefore, our reported values are more robust than the ones reported by Costa et al. (2021). Our measurements have larger uncertainties than those of Costa et al. (2021) at both 251 and 340 GHz. These larger uncertainties are due to the convolution to the beam size of 0 55× 0 32. The differences in beam area between the original and convolved images are 6.5 and 2.0 at 251 and 340 GHz, respectively. Convolution to the larger beam size introduces a larger uncertainty. Quantitatively, the 251 GHz image with the original synthesized beam size of 0 18× 0 14 has an rms noise of 12.4 μJy beam−1, but its rms noise increases to 24.7 μJy beam−1 after convolution to a synthesized beam size of 0 55× 0 32. Similarly, the 340 GHz image has an rms noise of 27.7 μJy beam−1 originally, but increases to 38.7 μJy beam−1 after convolution. Our measurements at 251 and 340 GHz had larger uncertainties, but it is important to match the beam sizes across the frequencies in order to accurately measure the spectral index. Additionally, it is important to measure the flux density of the BH at the correct position, rather than the flux density of a nearby yet different source. As a result, the flux densities we report are more accurately measured than previous values. In order to estimate the best-fit parameters and corresp- onding uncertainties for an SED ranging from 1.4 to 340 GHz, we employed the R function nls for nonlinear least-squares fitting and utilized confint2 to estimate the uncertainties of the best-fit parameters. Four models were tested, including: (i) a three-component model with nonthermal, free–free, and dust emission, n n n~ + +a -S S S S0,nth 0,ff 0.1 0,dust 4nth ; (ii) a two- component model with nonthermal and free–free emission, n n~ +a -S S S0,nth 0,ff 0.1nth ; (iii) a two-component model with nonthermal and dust emission, n n~ +aS S S0,nth 0,dust 4nth ; and (iv) a one-component model with nonthermal emission, n~ aS S0,nth nth. The models that included dust emission did not yield physical results. The two-component model with nonthermal and free–free emission had large uncertainties and poorly constrained solutions. The one-component model with nonthermal emission provided a meaningful solution, with a best-fit model of S = (1349± 217)ν−0.52±0.06 μJy, represented by the red line with the light red shaded region in Figure 3(E). 3.2. Molecular Line Emission 3.2.1. Moment Maps Moment maps of molecular gas can provide insight into the morphology and kinematics of the gas. Moment maps were generated using the 3D-Based Analysis of Rotating Objects via Line Observations (3D-Barolo; Di Teodoro & Fraternali 2015), which employs the Duchamp algorithm to identify pixels above a primary signal-to-noise ratio (S/N) and agglomerates neighboring pixels above a secondary S/N. The initial moment maps were created by running 3D-Barolo on the nonprimary beam attenuation-corrected image cubes, and the final moment maps were derived by multiplying the primary beam attenuation- corrected image cubes by the mask image cubes created during the generation of the initial moment maps. For this analysis, the primary and secondary S/N were set at 5 and 3, respectively. In Figure 4, the moment maps of the CO (1–0), CO (3–2), HCN (1–0), HCO+ (1–0), and HCO+ (3–2) lines are shown from top to bottom. The total intensity (moment 0), velocity field (moment 1), and velocity dispersion (moment 2) maps are displayed in the left, middle, and right panels, respectively. The black contours indicate the 99 GHz continuum, with the BH located at the local maximum between two major blobs of star formation. The moment maps of the molecular gas lines indicate that there is a significant amount of molecular gas beyond the continuum regions traced by the black contours. These extended features are consistent with previous observa- tions (Baas et al. 1994; Beck et al. 2018; Imara & Faesi 2019). Total intensity maps of the molecular gas (the left panels of Figure 4) show two prominent regions corresponding to the strong continuum blobs of Figure 2. The maximum intensities of the molecular gas in the left panels of Figure 4 coincide with the maximum positions in the continuum. The BH is located between the two prominent starburst regions in the molecular gas distributions, forming a spur-like structure visible in the high-resolution CO data. The detections of HCO+ (1–0), HCO+ (3–2), and HCN (1–0) above S/N > 5 indicate the existence of significant dense gas (with critical density 106 cm−3) in the vicinity of the BH. 3.2.2. Position–Velocity Diagram The position–velocity diagram (PVD) serves as a valuable tool for elucidating the kinematics of the gas. The ionized gas in the central region of He 2–10 exhibits a sinusoidal motion in the PVD, which is well fit by the bipolar outflow model presented in Schutte & Reines (2022). Schutte & Reines (2022) used HST/Space Telescope Imaging Spectrograph (STIS) spectroscopy, employing a long slit centered on the BH with a width of 0 2 and a position angle of 100° aligned with the ionized filament seen in the Hα imaging. In our investigation, we look for a similar pattern in the velocities of the molecular gas by generating PVDs along a line with a position angle of 100°, centered on the BH, from our CO (3–2) and CO (1–0) image cubes. We extracted PVDs along the line in bins corresponding to the synthesized beam sizes of the CO (3–2) and CO (1–0) image cubes, i.e., 0 3 and 0 45, respectively. However, these values are more than three times larger than the spatial resolution of the HST/STIS observations (∼0 1). Figure 5 illustrates the resulting PVDs for CO (3–2) in panel (A) and CO (1–0) in panel (B). Unlike the ionized gas observed by Schutte & Reines (2022), no sinusoidal motions were discerned in either PVD. This may be attributed to the coarser angular resolution of our observations, leading to a beam smearing effect that obscures outflow signals amid adjacent background emissions. Alternatively, the absence of sinusoidal motion in the molecular gas could signify a genuine lack of such motion, suggesting a potential insensitivity of the molecular gas to the BH outflow, particularly if the BH outflow lacks sufficient strength (Koudmani et al. 2022). We calculated the velocity gradient from the PVD of the CO (3–2) line in Figure 5(A), leveraging the smaller beam size of the 6 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines CO (3–2) image cube compared to that of CO (1–0). We restrict this analysis to the gas in the central region near the BH between −0 4 (−18 pc) and 0 8 (35 pc) with velocities in the range ∼870–900 km s−1. The resulting weighted mean velocity (the red squares in Figure 5) determined above 3σ with a bin width of 0 3 —where σ= 1.0 mJy beam−1 (1.693 km s−1)−1—yielded a velo- city gradient of −5.8 km s−1 arcsec−1. We repeated this analysis for CO (1–0), constructing a PVD along a line with a position Figure 4. Moment maps of CO (1–0), CO (3–2), HCN (1–0), HCO+ (1–0), and HCO+ (3–2), from top to bottom. Total intensity, velocity field, and velocity dispersion maps are shown in the left, middle, and right panels, respectively. The black contours indicate the flux density of the 99 GHz continuum, corresponding to 50 (5σ), 100 (10σ), 150 (15σ), and 200 (20σ) μJy per beam size of 0 55 × 0 3. The BH is located at the center of the innermost black contour of the continuum. The beam sizes are marked with solid ellipses in the bottom right corners. 7 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines angle of 100° and a width of 0 45 (matching the beam size of CO (1–0)), resulting in a similar velocity gradient of −4.6 km s−1 arcsec−1). Our derived velocity gradient is shallower than the −7.8 km s−1 arcsec−1 reported by Imara & Faesi (2019). The discrepancy can be attributed to differences in: (i) the length of the position axis used to estimate the velocity gradient; (ii) the position angle; and (iii) the bin size used to trace the gas velocity. Imara & Faesi (2019) estimated the velocity gradient over a range of −1 7 (∼70 pc) to 1 7 along a line with a position angle of 143° and a width of 4″. In our PVD (Figure 5), the molecular gas at approximately 1 5 is associated with the western star-forming region, contributing to a steeper velocity gradient due to its high velocity. The differences in position angle and bin size introduce additional inconsistencies, particularly given the asymmetrical distribu- tion of the molecular gas around the BH. 3.2.3. Molecular Gas in the Vicinity of the BH In this work, we are primarily interested in examining the molecular gas in the vicinity of the BH.5 We define this region (i.e., the BH vicinity) as follows. There is a distinct velocity gas blob near the BH in HCO+ (3–2), seen in the bottom middle panel of Figure 4. The velocity range of this blob is 852.1–904.6 km s−1 in HCO+(3–2). To help isolate the gas in the vicinity of the BH, both spatially and in velocity, we generated moment maps of the molecular gas tracers with higher angular resolution (CO (1–0), CO (3–2), and HCO+ (3–2)) within this restricted velocity range (Figure 6). The resulting moment maps reveal an elongated CO gas structure with a distinct velocity from the surrounding regions, Figure 5. PVDs of CO J = 3–2 (A) and CO J = 1–0 (B) are presented, with the weighted mean velocity in each position bin above 3σ (red squares) and the velocity gradient (black lines). The light blue contours indicate 3σ, 10σ, and 50σ in CO (3–2) and 3σ, 5σ, and 10σ in CO (1–0), respectively. The dotted vertical lines present the regions estimating the velocity gradient, e.g., −0 4 < x < 0 8 for CO (3–2) and −0 5 < x < 0 85 for CO (1–0). 5 The ALMA observations do not have sufficient angular resolution to probe the BH sphere of influence. Given the BH mass of ∼106 Me (Reines et al. 2016; Riffel 2020), the sphere of influence is ∼4 pc (0 1), which is at least three times smaller than the angular resolution of our observations. 8 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines as seen in the middle panels of Figure 6. This distinct velocity region is defined as the BH vicinity, represented by the black ellipse in Figure 6. The ellipse is centered at R.A. 08h 36m 15 108 and decl. −26° 24′ 34 343, with a size of 2 97× 0 79 (≈128 pc× 34 pc) and a position angle of 80°.34, as presented in Table 3. The moment maps in Figure 6 clearly demonstrate that the BH (small black contour) is located north of the BH vicinity (black ellipse) in molecular gases such as CO and HCO+ (3–2). In other words, the BH is not at the center of the molecular gas structure. The lack of molecular gas to the north of the BH is likely a result of star formation feedback occurring in the young massive star clusters with ages of ∼4–5Myr (Beck et al. 2018), as seen in the Hα image in Figure 7. This could help explain why the BH is accreting at such a low Eddington ratio (∼10−6; Reines et al. 2016). Furthermore, the asymmetric distribution of molecular gas around the massive BH resembles the central molecular zone of the Milky Way (Bally et al. 1988) and the central region of NGC 253 (Meier et al. 2015; Krieger et al. 2019; Levy et al. 2022). In the cases of the Milky Way and NGC 253, the origin of this asymmetry is attributed to feedback from star formation in their centers (Krieger et al. 2019; Henshaw et al. 2022). We have measured the velocity ranges within the BH vicinity using three different molecular lines: CO J= 1–0, CO J= 3–2, and HCO+ J= 3–2. The velocity ranges were found to be 865–899 km s−1, 871–904 km s−1, and 871–904 km s−1, respectively. Excluding the borders, the velocity ranges for CO J= 1–0 and J= 3–2 were reduced to 871–886 km s−1 and 871–885 km s−1, respectively. The mean velocities within the Figure 6. Moment maps of the molecular gas species observed at high angular resolution including CO (1–0), CO (3–2), and HCO+ (3–2) (top to bottom), where the moment maps are generated with a restrictive velocity range of 852.1–904.6 km s−1 (see Section 3.2). The total intensity, velocity field, and velocity dispersion are presented in the left, middle, and right panels, respectively. The small black contour indicates the position of the BH traced by the high-brightness-temperature radio core detected with the LBA (Reines & Deller 2012). The outer black ellipse defines the region we refer to as the BH vicinity, which has a near uniform velocity that is consistent with the systemic velocity of the galaxy, Vsys = 873 km s−1 (Kobulnicky et al. 1995). Table 3 Properties of Molecular Gas in the Vicinity of the BH (Black Ellipse in Figure 6) Parameter Value Center coordinate R.A. 08h 36m 15 108, decl. −26° 24′ 34 343 Major and minor axes 2 97 × 0 79 (≈128 pc × 34 pc) Position angle 83°. 34 Mean velocity 880.2 ± 4.7 km s−1 in CO (1–0) 880.3 ± 5.9 km s−1 in CO (3–2) 875.3 ± 6.2 km s−1 in HCO+ (3–2) H2 mass (0.21–1.53) × 106 Me 9 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines BH vicinity are 880± 5 (CO J= 1–0), 880± 6 (CO J= 3–2), and 875± 6 km s−1 (HCO+ J= 3–2), and are therefore consistent with one another within the uncertainties. For reference, the systemic velocity of the galaxy is 873 km s−1 (Kobulnicky et al. 1995). The velocity field maps do not show any significant rotation in the region defined as the BH vicinity. The velocity dispersion map of CO J= 3–2 shows some high- velocity dispersion along the SE and NW directions with Δσ≈ 7 km s−1. The high-velocity dispersion (Δσ≈ 50 km s−1) of stellar motions near the BH observed by Riffel (2020) is not seen in the molecular gas. The molecular gas mass within the BH vicinity is estimated using the scaling relation between the molecular gas mass and CO luminosity, a= ¢M LH CO CO2 . The conversion factor αCO depends on the metallicity (Bolatto et al. 2013 and references therein). Cresci et al. (2017) conducted spectro- scopic observations with MUSE to determine the metallicity of He 2–10, with a spatial sampling size of 0 2× 0 2. They created a map of the metallicity distribution, which revealed significant variation across the galaxy, largely due to changes in the ionization parameters. Through simultaneous line fitting on multiple lines, they found that the metallicity in the central region was supersolar, with 12+ [O/H]≈ 9.0, corresponding to Z= 0.044. Based on the relationship between αCO and metallicity outlined in Bolatto et al. (2013), we choose a range of αCO values between 0.6 and 4.3Me (K km s−1 pc2)−1, where αCO= 4.3 is the value for the Milky Way. Our range of αCO also includes the widely accepted value for the centers of starburst galaxies with solar metallicities, αCO≈ 1.1Me (K km s−1 pc2)−1 (Leroy et al. 2015; Krieger et al. 2019). These considerations ensure that our analysis is consistent with established standards in the field. Our estimated CO luminosity is ¢ =LCO ´ -3.56 10 K km s pc5 1 2, which, when combined with the αCO range, results in an estimated molecular gas mass of = ´M 0.21 1.53 10H 6 2 ( – ) Me. 3.2.4. CO Line Ratios The line ratios between different molecular gas tracers are good indicators of gas excitation (Privon et al. 2015). The HCN (1–0) and HCO+ (1–0) line image cubes have large beam sizes, and the HCO+ (3–2) line image cube has a relatively low S/N, so we restricted our analysis to the line ratio between CO J= 1–0 and J= 3–2, which have higher angular resolutions and better sensitivities. Given that the CO (3–2) line emission exhibits a more extended distribution compared to the CO (1–0) line emission, the distribution of the CO line ratio is constrained by the CO (1–0) line distribution. The line ratio is defined as n n= =- - - - -R L L I I ,31 CO 3 2 CO 1 0 32 10 2 CO 3 2 CO 1 0( ) ( )( ) ( ) ( ) ( ) where ICO is the velocity-integrated line intensity in units of Jy km s−1. Both image cubes were convolved to have the same synthesized beam size of 0 45× 0 30 using the CASA task imsmooth. We present the pixel-by-pixel CO line ratio map in the left panel and an HST Hα image (e.g., Schutte & Reines 2022) in the right panel of Figure 7. Our line ratio map of the central region of the galaxy indicates highly excited CO molecular gas, as evidenced by R31 1, exceeding the typical range observed in the Milky Way, 0.3< R31< 0.5 (Leroy et al. 2009; Carilli & Walter 2013). This result aligns with the findings of Beck et al. (2018), who reported a high flux ratio of SCO(3−2)/SCO(1−0)= 6.1 for this galaxy. 3.2.5. Shocked Molecular Clouds from the BH Outflow The discovery of a low-velocity bipolar outflow originating from the BH was reported in ionized gas (Schutte & Reines 2022). The precession angle and frequency of this outflow were estimated to be in the range of 2°.4–6°.1 and 3.0–7.5 revolutions per Myr, respectively. Notably, Schutte & Reines (2022) identified several impacts of the BH outflow, including the stimulation of star formation in the eastern star- forming region, as well as in the dark cloud (region 5) and the dense clouds within the western star-forming region (region 6). The left panel of Figure 7 reveals that the highest R31 line ratios Figure 7. Left: pixel-by-pixel line ratio map between CO (3–2) and CO (1–0), R31. Right: HST Hα plus continuum image with the same field of view. The BH position is marked with black contours from the LBA 1.4 GHz observations by Reines & Deller (2012), and the region defined as the BH vicinity is indicated with the orange dashed ellipse. Regions 5 and 6 from Schutte & Reines (2022) are presented as rectangles with the HST slit position for the ionized gas outflow. The highest line ratios are observed at the interface between the BH outflow and the eastern star-forming region, as well as region 5 to the west of the BH. This implies the cold molecular gas in these regions is shocked by the bipolar outflow, supporting the case for positive BH feedback. 10 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines are observed at the interface between the known ionized gas outflow from the BH and the eastern star-forming region (Schutte & Reines 2022), as depicted in the Hα image (the right panel of Figure 7). The line ratio at this interface is R31= 2.7± 0.2, while the region defined as the BH vicinity excluding this interface has R31= 1.3± 0.1. Additionally, we detect a high line ratio (R31= 2.5± 0.7) to the west of the BH at the boundary of a dark cloud of gas and dust (region 5 in Figure 7). This region was previously identified by Schutte & Reines (2022) and exhibits double-peaked optical emission lines, likely due to the other side of the bipolar BH outflow intercepting dense clouds and pushing them laterally (as is the case with the eastern interface region). The CO line ratio at region 5 is comparable to that of the interface between the outflow and the eastern star-forming region, consistent with this scenario. These findings suggest that the molecular gas impacted by the BH outflow to the east and west has both high density and high kinetic temperature (Carilli & Walter 2013). It is noteworthy that the elevated R31 may be an edge effect limited by low flux densities of CO (1–0). Nevertheless, certain regions at the peripheries, particularly in the northeast and south, exhibit low R31 values. The presence of these regions implies that the pixel-by-pixel R31 values are influenced by the relative distributions of CO (1–0) and CO (3–2) due to some physical origins, rather than being solely dictated by low CO (1–0) flux densities. Moreover, our derived R31 values are robust, exhibiting similar sensitivities, e.g., σ= 1.03 mJy beam−1 (1.693 km s−1)−1 for CO (1–0) and 1.00 mJy beam−1 (1.693 km s−1)−1 for CO (3–2), respectively. Therefore, the elevated R31 values observed in this study are not solely constrained by low CO (1–0) flux densities. Multiple studies have investigated the occurrence of elevated CO line ratios in different galaxies. García-Burillo et al. (2014) and Viti et al. (2014) examined the CO line ratio in NGC 1068 and demonstrated that the ratio between J= 3–2 and J= 1–0 varies significantly between the circumnuclear disk around the AGN (R31= 2.7) and the starburst ring (R31= 1.2). Further- more, they observed higher CO line ratios in more active star- forming regions compared to less active ones (R31= 0.7–1.0). Their analysis suggested that the elevated CO line ratio in the circumnuclear disk was a result of the AGN-driven outflow. The CO line ratio in IC 5063 was extensively studied by Dasyra et al. (2016) and Oosterloo et al. (2017). They identified excited molecular gas along the radio jet and noted a higher CO line ratio in the central region compared to the outer region. For instance, Dasyra et al. (2016) conducted a calculation of the line ratio between CO (4–3) and CO (2–1), revealing a value of 12 in the robust jet–cloud interacting regions compared to an average of 5. Oosterloo et al. (2017) demonstrated that R31> 1 along the jets, as opposed to 0.17 in the disk. Their study further indicated that the most elevated CO line ratio is proximate to the two radio lobes, implicating the shock generated by the outflow as the causative factor (Morganti et al. 2021). More recently, Audibert et al. (2023) identified a high CO line ratio between J = 3–2 and J = 2–1, specifically with TCO(3−2)/TCO(2−1)= 0.8 in the region perpendicular to the radio jet, in contrast to the value of 0.4 in other regions within the type 2 QSO Teacup galaxy. Their findings underscored that this augmented CO line ratio is a consequence of the jet-driven shock interacting with the molecular gas. Our observations align with the results reported by García- Burillo et al. (2014), where the R31 values stand at 2.7 in the circumnuclear disk, 1.2 in starburst regions, and 0.4−0.5 in the halo. Consequently, the regions exhibiting high R31 values in our study can be attributed to shock heating induced by the BH outflow. 3.2.6. The Western Star-forming Region While prominent at radio wavelengths (e.g., Figure 2), the star- forming region to the west of the BH appears to suffer from significant extinction at optical wavelengths (e.g., Figures 8 and 7). Schutte & Reines (2022) estimate a stellar age of <3Myr for region 6 (e.g., see Figure 7), based on the equivalent width of the Hα emission line. This is consistent with a scenario in which the star clusters have not had enough time to clear away their birth material, leading to high levels of extinction. We quantify their assertion with the molecular gas distribution. Figure 8 presents channel maps of the CO (3–2) line within velocities of 864–888 km s−1 in panels (a) to (o), a total intensity map in panel (p), an HST I-band image in panel (q), and an HST Hα plus continuum image in panel (r). The total intensity map was produced using 3D-Barolo (Di Teodoro & Fraternali 2015), with a primary S/N of 5 and a secondary S/N of 3 on the image cube shown in panels (a)–(o). In the channel maps and total intensity map, the bulk of the gas was masked to enhance the contrast and highlight fainter emission. The red boxes in each panel of Figure 8 (the white boxes in panel (r)) indicate regions 5 and 6 defined by Schutte & Reines (2022), respectively. Figure 8 demonstrates that there is a significant amount of molecular gas between regions 5 and 6, implying that the molecular gas is abundant enough to absorb the emission from young stellar populations embedded in the gas. We estimate the optical extinction in this region as follows. The total intensity of the CO (1–0) line within the span of regions 5 and 6 is measured as WCO(1−0) = 788.2± 57.4 K km s−1, resulting in an estimated molecular gas column density of N(H2) = (2.21± 0.16)× 1022− (1.58± 0.12)× 1023 cm−2. This estimation is derived from the mean velocity-integrated flux density ratio of CO (3–2) to CO (1–0), represented as ICO(3−2)/ICO(1−0)= 12.54 in this region. The conversion factor XCO= (0.28− 2)× 1020 cm−2 (Kkm s−1)−1, corresponding to αCO= 0.6–4.3Me (Kkm s−1 pc2)−1, has been employed for this calculation. Kobulnicky et al. (1995) conducted observations of atomic hydrogen (H I) using the D-, C-, and B-arrays of the VLA. In this region, the total intensity of H I was 1.52± 0.05 Jy beam−1 km s−1, leading to a corresponding column density N(H I) at (1.51± 0.05)× 1018 cm−2, considering a synthesized beam size of 31″× 30″. Assuming a uniform distribution of H I, the optical extinction can be computed using the relationship between N(H) and AV (Zhu et al. 2017), AV= 4.81× 10−22(N(H I)+ 2N(H2)). This calculation yields a range of AV∼ 21–152mag. If the H Iis more concentrated in this region, these values would represent lower limits. On the other hand, the optical extinction is estimated to be AV∼ 7–51 for the eastern star-forming region, supporting the optically dark region being predominantly observed between regions 5 and 6. The 251 and 340 GHz images presented in Figure 2 showed that the compact features of region 6 align with the peaks in CO (3–2) and HCO+ (3–2) shown in Figure 4. This observation is in line with the characteristics of bright and compact embedded 11 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines clusters seen in other young clusters, e.g., NGC 253 (Ando et al. 2017; Leroy et al. 2018; Levy et al. 2021; Mills et al. 2021) and NGC 4945 (Emig et al. 2020). To quantify our findings, we focused on the radio continuum blob located in the western star-forming region, corresponding to region 6 as defined in Schutte & Reines (2022). Dust and gas mass Figure 8. We present the channel maps of the CO (3–2) line within the velocity range of 864.1–887.8 km s−1. In order to highlight the faint emission beyond the BH vicinity (indicated by the dashed orange line), we have masked the majority of CO gas. The 99 GHz continuum emission is depicted by the black contours at levels of 50, 100, 200, and 300 μJy beam−1. The BH is represented by a solid black ellipse, corresponding to the central maximum of the 1.4 GHz continuum observed by the LBA (Reines & Deller 2012). The diagonal dotted lines illustrate the spectroscopic slit used in the HST observations, with a width of 0 2. Additionally, regions 5 (left) and 6 (right) previously discussed by Schutte & Reines (2022) are depicted by red boxes. The ALMA beam size is provided in the bottom right corner of each panel. The bottom row exhibits the total intensity map of the CO (3–2) line on the left, followed by the HST I-band image in the middle, and the HST Hα+continuum image on the right. In the Hα+continuum image (r), the white contours represent the 99 GHz continuum map, the cyan dotted lines depict the slit, and the green boxes highlight regions 5 and 6 for improved visibility. 12 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines estimations were carried out within this region. For the calculation of the dust mass, we utilized the 340 GHz flux density and employed the equation k = n n n M S D B T 1 d dust 2 ( ) ( ) (Launhardt & Henning 1997). We adopted a dust absorption coefficient of κν= 1.3 cm2 g−1 at 340 GHz (Costa et al. 2021) and the dust temperature of 29K (Santangelo et al. 2009). The measured flux density at 340GHz is S340 GHz= 610.6± 30.0μJy, resulting in a calculated dust mass of (2.37± 0.12)× 103Me. The associated molecular mass is =  ´M 7.79 0.15H2 ( ) -  ´10 5.65 0.11 104 5( ) Me, depending on the αCO. Conse- quently, the dust-to-molecular-gas-mass ratio falls within the range of 1:33–238, which is consistent with the typical value of the Milky Way and other starburst galaxies, 1:120 (Wilson et al. 2008). The existence of abundant molecular gas and a high optical extinction value in the optically dark regions (regions 5 and 6) indicates that the molecular gas in these regions is compressed. In He 2–10, the shock generated by the BH outflow emerges as the most plausible source for the compression of molecular gas, contributing to the formation of embedded star clusters. This is an example of the “positive feedback” of the BH outflow. The asymmetric distributions of molecular gas and embedded star formation provide a key to understanding the evolution of He 2–10. Given that He 2–10 has experienced a tidal interaction (Kobulnicky et al. 1995), the asymmetric distributions of molecular gas result from the tidal interaction by disturbing the galaxy potential. This asymmetric distribution of molecular gas might introduce the asymmetric embedded star formations via the asymmetric gas inflow, similar to the Central Molecular Zone of the Milky Way (Battersby et al. 2020; Henshaw et al. 2022). 4. Conclusions We reimaged the ALMA archival data of He 2–10 using newly developed algorithms and studied properties of the continuum and molecular gas, including CO (1–0), CO (3–2), HCN (1–0), HCO+ (1–0), and HCO+ (3–2), focusing on the massive BH and its vicinity. We reanalyzed the radio– submillimeter SED of the BH with updated flux densities and investigated the kinematics and morphology of molecular gas around the BH. Our main results are summarized below. 1. The radio–submillimeter SED of the BH from 1.4 to 340 GHz is dominated by synchrotron emission, with a best-fit model of Sν= (1349± 217)ν−0.52±0.06 μJy (Figure 3). 2. We find an elongated molecular gas structure (≈128 pc× 34 pc) in the vicinity of the BH, with a distinct velocity component in CO (1–0), CO (3–2), and HCO+ (3–2) from the surrounding regions (Figure 6). The mass of the molecular gas structure in the vicinity of the BH is = ´M M0.21 1.53 10H 6 2 ( – ) for a range of metallicity-dependent values of αCO. While the position of the BH is within this region, it is significantly offset from the peak intensity, which could explain why the BH is accreting at a low rate. 3. The ratio of CO (3–2) to CO (1–0) shows that the CO gas is highly excited in the BH vicinity (Figure 7). The highest excitation values are located at the interface between the BH outflow and the eastern region of triggered star formation identified by Schutte & Reines (2022). Our results suggest that the molecular gas is being compressed by the shock produced by the outflow from the BH, providing additional evidence for positive BH feedback. 4. The high excitation of the CO line ratio is also detected to the west of the BH, suggesting the other side of the bipolar flow is shocking dense molecular clouds in the western star-forming region as well. We find high molecular gas column densities (N(H2)∼ 1021−22 cm−2) in this region, as well as high levels of extinction (AV 2–11 mag). This is consistent with a scenario in which the infant star clusters with stellar ages 3Myr (Schutte & Reines 2022) have not yet had enough time to destroy or disperse the dense molecular clouds from which they formed. The study presented in this paper constitutes the first investigation of the molecular gas properties in the immediate vicinity of the massive BH in He 2–10. While our results provide additional support for positive BH feedback, future observations with higher angular resolution are needed to resolve the sphere of influence of the BH and trace any potential molecular gas outflow directly. Furthermore, analyz- ing the higher transitions of the CO molecular gas and applying models for the photodissociation region and/or X-ray-domi- nated region will provide additional information on physical parameters such as density, temperature, and outflow power, ultimately yielding more insight into the impact of BH feedback on nearby molecular gas. Acknowledgments The authors express their gratitude to the anonymous reviewer for valuable contributions enhancing the quality of this paper. We are grateful to Drs. John Hibbard, Loreto Barcos-Muñoz, Mark Lacy, Allison Costa, George Privon at NRAO, Dr. Seth Kimbrell at Carleton College, and Sheyda Salehirad at Montana State University for valuable discussions. A.E.R. acknowledges support for this work provided by NASA through EPSCoR grant No. 80NSSC20M0231. This paper makes use of the following ALMA data: ADS/JAO.ALMA #2011.0.00348.S, ADS/JAO. ALMA #2012.1.00413.S, ADS/JAO.ALMA #2015.1.01569. S, ADS/JAO.ALMA #2016.1.00027.S, ADS/JAO.ALMA #2016.1.00492.S, and ADS/JAO.ALMA #2019.1.01641.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. This research has made use of NASA’s Astrophysics Data System Bibliographic Services. Facilities: ALMA, VLA, HST. Software: CASA (Bean et al. 2022), AIPS (Greisen 2003), CARTA (Comrie et al. 2021), R (R Core Team 2022), 3D- Barolo (Di Teodoro & Fraternali 2015). ORCID iDs Hansung B. Gim https://orcid.org/0000-0003-1436-7658 Amy E. Reines https://orcid.org/0000-0001-7158-614X 13 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0003-1436-7658 https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X https://orcid.org/0000-0001-7158-614X References Allen, D. A., Wright, A. E., & Goss, W. M. 1976, MNRAS, 177, 91 Ando, R., Nakanishi, K., Kohno, K., et al. 2017, ApJ, 849, 81 Audibert, A., Ramos Almeida, C., García-Burillo, S., et al. 2023, A&A, 671, L12 Baas, F., Israel, F. P., & Koornneef, J. 1994, A&A, 284, 403 Bally, J., Stark, A. A., Wilson, R. W., & Henkel, C. 1988, ApJ, 324, 223 Battersby, C., Keto, E., Walker, D., et al. 2020, ApJS, 249, 35 Bayet, E., Gerin, M., Phillips, T. G., & Contursi, A. 2004, A&A, 427, 45 Bean, B., Bhatnagar, S., Castro, S., et al. 2022, PASP, 134, 114501 Beck, S. C., Turner, J. L., & Consiglio, S. M. 2018, ApJ, 867, 165 Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 Carilli, C. L., & Walter, F. 2013, ARA&A, 51, 105 Comrie, A., Wang, K.-S., Hsu, S.-C., et al. 2021, CARTA: The Cube Analysis and Rendering Tool for Astronomy, v2.0.0, Zenodo, doi:10.5281/zenodo. 4905459 Costa, A. H., Johnson, K. E., Indebetouw, R., et al. 2021, ApJ, 918, 76 Cresci, G., Vanzi, L., Telles, E., et al. 2017, A&A, 604, A101 Dasyra, K. M., Combes, F., Oosterloo, T., et al. 2016, A&A, 595, L7 Di Teodoro, E. M., & Fraternali, F. 2015, MNRAS, 451, 3021 Emig, K. L., Bolatto, A. D., Leroy, A. K., et al. 2020, ApJ, 903, 50 García-Burillo, S., Combes, F., Usero, A., et al. 2014, A&A, 567, A125 Greene, J. E., Strader, J., & Ho, L. C. 2020, ARA&A, 58, 257 Greisen, E. W. 2003, in Information Handling in Astronomy—Historical Vistas, ed. A. Heck, Vol. 285 (Dordrecht: Kluwer), 109 Henshaw, J. D., Krumholz, M. R., Butterfield, N. O., et al. 2022, MNRAS, 509, 4758 Imanishi, M., Nakanishi, K., Tamura, Y., Oi, N., & Kohno, K. 2007, AJ, 134, 2366 Imara, N., & Faesi, C. M. 2019, ApJ, 876, 141 Johnson, K. E., Brogan, C. L., Indebetouw, R., et al. 2018, ApJ, 853, 125 Johnson, K. E., & Kobulnicky, H. A. 2003, ApJ, 597, 923 Kobulnicky, H. A., Dickey, J. M., Sargent, A. I., Hogg, D. E., & Conti, P. S. 1995, AJ, 110, 116 Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511 Koudmani, S., Sijacki, D., & Smith, M. C. 2022, MNRAS, 516, 2112 Krieger, N., Bolatto, A. D., Walter, F., et al. 2019, ApJ, 881, 43 Launhardt, R., & Henning, T. 1997, A&A, 326, 329 Leroy, A. K., Bolatto, A. D., Ostriker, E. C., et al. 2015, ApJ, 801, 25 Leroy, A. K., Bolatto, A. D., Ostriker, E. C., et al. 2018, ApJ, 869, 126 Leroy, A. K., Walter, F., Bigiel, F., et al. 2009, AJ, 137, 4670 Levy, R. C., Bolatto, A. D., Leroy, A. K., et al. 2021, ApJ, 912, 4 Levy, R. C., Bolatto, A. D., Leroy, A. K., et al. 2022, ApJ, 935, 19 Meier, D. S., Turner, J. L., Crosthwaite, L. P., & Beck, S. C. 2001, AJ, 121, 740 Meier, D. S., Walter, F., Bolatto, A. D., et al. 2015, ApJ, 801, 63 Mills, E. A. C., Gorski, M., Emig, K. L., et al. 2021, ApJ, 919, 105 Mohan, N., & Rafferty, D. 2015, PyBDSF: Python Blob Detection and Source Finder, Astrophysics Source Code Library, ascl:1502.007 Morganti, R., Oosterloo, T., Murthy, S., & Tadhunter, C. 2021, AN, 342, 1135 Nguyen, D. D., Seth, A. C., Reines, A. E., et al. 2014, ApJ, 794, 34 Oosterloo, T., Raymond Oonk, J. B., Morganti, R., et al. 2017, A&A, 608, A38 Privon, G. C., Herrero-Illana, R., Evans, A. S., et al. 2015, ApJ, 814, 39 Rau, U., & Cornwell, T. J. 2011, A&A, 532, A71 R Core Team 2022, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, https://www.R-project.org/ Reines, A. E. 2022, NatAs, 6, 26 Reines, A. E., & Deller, A. T. 2012, ApJL, 750, L24 Reines, A. E., Reynolds, M. T., Miller, J. M., et al. 2016, ApJL, 830, L35 Reines, A. E., Sivakoff, G. R., Johnson, K. E., & Brogan, C. L. 2011, Natur, 470, 66 Reines, A. E., & Volonteri, M. 2015, ApJ, 813, 82 Riffel, R. A. 2020, MNRAS, 494, 2004 Santangelo, G., Testi, L., Gregorini, L., et al. 2009, A&A, 501, 495 Schutte, Z., & Reines, A. E. 2022, Natur, 601, 329 Somerville, R. S., & Davé, R. 2015, ARA&A, 53, 51 Viti, S., García-Burillo, S., Fuente, A., et al. 2014, A&A, 570, A28 Volonteri, M., Habouzit, M., & Colpi, M. 2021, NatRP, 3, 732 Wilson, C. D., Petitpas, G. R., Iono, D., et al. 2008, ApJS, 178, 189 Zhu, H., Tian, W., Li, A., & Zhang, M. 2017, MNRAS, 471, 3494 14 The Astrophysical Journal, 963:103 (14pp), 2024 March 10 Gim & Reines https://doi.org/10.1093/mnras/177.1.91 https://ui.adsabs.harvard.edu/abs/1976MNRAS.177...91A/abstract https://doi.org/10.3847/1538-4357/aa8fd4 https://ui.adsabs.harvard.edu/abs/2017ApJ...849...81A/abstract https://doi.org/10.1051/0004-6361/202345964 https://ui.adsabs.harvard.edu/abs/2023A&A...671L..12A/abstract https://ui.adsabs.harvard.edu/abs/2023A&A...671L..12A/abstract https://ui.adsabs.harvard.edu/abs/1994A&A...284..403B/abstract https://doi.org/10.1086/165891 https://ui.adsabs.harvard.edu/abs/1988ApJ...324..223B/abstract https://doi.org/10.3847/1538-4365/aba18e https://ui.adsabs.harvard.edu/abs/2020ApJS..249...35B/abstract https://doi.org/10.1051/0004-6361:20035614 https://ui.adsabs.harvard.edu/abs/2004A&A...427...45B/abstract https://doi.org/10.1088/1538-3873/ac9642 https://ui.adsabs.harvard.edu/abs/2022PASP..134k4501C/abstract https://doi.org/10.3847/1538-4357/aae481 https://ui.adsabs.harvard.edu/abs/2018ApJ...867..165B/abstract https://doi.org/10.1146/annurev-astro-082812-140944 https://ui.adsabs.harvard.edu/abs/2013ARA&A..51..207B/abstract https://doi.org/10.1146/annurev-astro-082812-140953 https://ui.adsabs.harvard.edu/abs/2013ARA&A..51..105C/abstract https://doi.org/10.5281/zenodo.4905459 https://doi.org/10.5281/zenodo.4905459 https://doi.org/10.3847/1538-4357/ac0e93 https://ui.adsabs.harvard.edu/abs/2021ApJ...918...76C/abstract https://doi.org/10.1051/0004-6361/201730876 https://ui.adsabs.harvard.edu/abs/2017A&A...604A.101C/abstract https://doi.org/10.1051/0004-6361/201629689 https://ui.adsabs.harvard.edu/abs/2016A&A...595L...7D/abstract https://doi.org/10.1093/mnras/stv1213 https://ui.adsabs.harvard.edu/abs/2015MNRAS.451.3021D/abstract https://doi.org/10.3847/1538-4357/abb67d https://ui.adsabs.harvard.edu/abs/2020ApJ...903...50E/abstract https://doi.org/10.1051/0004-6361/201423843 https://ui.adsabs.harvard.edu/abs/2014A&A...567A.125G/abstract https://doi.org/10.1146/annurev-astro-032620-021835 https://ui.adsabs.harvard.edu/abs/2020ARA&A..58..257G/abstract https://ui.adsabs.harvard.edu/abs/2003ASSL..285..109G/abstract https://doi.org/10.1093/mnras/stab3039 https://ui.adsabs.harvard.edu/abs/2022MNRAS.509.4758H/abstract https://ui.adsabs.harvard.edu/abs/2022MNRAS.509.4758H/abstract https://doi.org/10.1086/523598 https://ui.adsabs.harvard.edu/abs/2007AJ....134.2366I/abstract https://ui.adsabs.harvard.edu/abs/2007AJ....134.2366I/abstract https://doi.org/10.3847/1538-4357/ab16cc https://ui.adsabs.harvard.edu/abs/2019ApJ...876..141I/abstract https://doi.org/10.3847/1538-4357/aa9ff8 https://ui.adsabs.harvard.edu/abs/2018ApJ...853..125J/abstract https://doi.org/10.1086/378585 https://ui.adsabs.harvard.edu/abs/2003ApJ...597..923J/abstract https://doi.org/10.1086/117500 https://ui.adsabs.harvard.edu/abs/1995AJ....110..116K/abstract https://doi.org/10.1146/annurev-astro-082708-101811 https://ui.adsabs.harvard.edu/abs/2013ARA&A..51..511K/abstract https://doi.org/10.1093/mnras/stac2252 https://ui.adsabs.harvard.edu/abs/2022MNRAS.516.2112K/abstract https://doi.org/10.3847/1538-4357/ab2d9c https://ui.adsabs.harvard.edu/abs/2019ApJ...881...43K/abstract https://ui.adsabs.harvard.edu/abs/1997A&A...326..329L/abstract https://doi.org/10.1088/0004-637X/801/1/25 https://ui.adsabs.harvard.edu/abs/2015ApJ...801...25L/abstract https://doi.org/10.3847/1538-4357/aaecd1 https://ui.adsabs.harvard.edu/abs/2018ApJ...869..126L/abstract https://doi.org/10.1088/0004-6256/137/6/4670 https://ui.adsabs.harvard.edu/abs/2009AJ....137.4670L/abstract https://doi.org/10.3847/1538-4357/abec84 https://ui.adsabs.harvard.edu/abs/2021ApJ...912....4L/abstract https://doi.org/10.3847/1538-4357/ac7b7a https://ui.adsabs.harvard.edu/abs/2022ApJ...935...19L/abstract https://doi.org/10.1086/318782 https://ui.adsabs.harvard.edu/abs/2001AJ....121..740M/abstract https://ui.adsabs.harvard.edu/abs/2001AJ....121..740M/abstract https://doi.org/10.1088/0004-637X/801/1/63 https://ui.adsabs.harvard.edu/abs/2015ApJ...801...63M/abstract https://doi.org/10.3847/1538-4357/ac0fe8 https://ui.adsabs.harvard.edu/abs/2021ApJ...919..105M/abstract http://www.ascl.net/1502.007 https://doi.org/10.1002/asna.20210037 https://ui.adsabs.harvard.edu/abs/2021AN....342.1135M/abstract https://doi.org/10.1088/0004-637X/794/1/34 https://ui.adsabs.harvard.edu/abs/2014ApJ...794...34N/abstract https://doi.org/10.1051/0004-6361/201731781 https://ui.adsabs.harvard.edu/abs/2017A&A...608A..38O/abstract https://doi.org/10.1088/0004-637X/814/1/39 https://ui.adsabs.harvard.edu/abs/2015ApJ...814...39P/abstract https://doi.org/10.1051/0004-6361/201117104 https://ui.adsabs.harvard.edu/abs/2011A&A...532A..71R/abstract https://www.R-project.org/ https://doi.org/10.1038/s41550-021-01556-0 https://ui.adsabs.harvard.edu/abs/2022NatAs...6...26R/abstract https://doi.org/10.1088/2041-8205/750/1/L24 https://ui.adsabs.harvard.edu/abs/2012ApJ...750L..24R/abstract https://doi.org/10.3847/2041-8205/830/2/L35 https://ui.adsabs.harvard.edu/abs/2016ApJ...830L..35R/abstract https://doi.org/10.1038/nature09724 https://ui.adsabs.harvard.edu/abs/2011Natur.470...66R/abstract https://ui.adsabs.harvard.edu/abs/2011Natur.470...66R/abstract https://doi.org/10.1088/0004-637X/813/2/82 https://ui.adsabs.harvard.edu/abs/2015ApJ...813...82R/abstract https://doi.org/10.1093/mnras/staa903 https://ui.adsabs.harvard.edu/abs/2020MNRAS.494.2004R/abstract https://doi.org/10.1051/0004-6361/200911688 https://ui.adsabs.harvard.edu/abs/2009A&A...501..495S/abstract https://doi.org/10.1038/s41586-021-04215-6 https://ui.adsabs.harvard.edu/abs/2022Natur.601..329S/abstract https://doi.org/10.1146/annurev-astro-082812-140951 https://ui.adsabs.harvard.edu/abs/2015ARA&A..53...51S/abstract https://doi.org/10.1051/0004-6361/201424116 https://ui.adsabs.harvard.edu/abs/2014A&A...570A..28V/abstract https://doi.org/10.1038/s42254-021-00364-9 https://ui.adsabs.harvard.edu/abs/2021NatRP...3..732V/abstract https://doi.org/10.1086/590910 https://ui.adsabs.harvard.edu/abs/2008ApJS..178..189W/abstract https://doi.org/10.1093/mnras/stx1580 https://ui.adsabs.harvard.edu/abs/2017MNRAS.471.3494Z/abstract 1. Introduction 2. ALMA Observations and Data Reduction 2.1. Spectral Line Image Cubes 2.2. Continuum Imaging 2.3. Astrometry 3. Analysis and Results 3.1. Continuum Emission and Radio Spectral Energy Distribution of the BH 3.2. Molecular Line Emission 3.2.1. Moment Maps 3.2.2. Position–Velocity Diagram 3.2.3. Molecular Gas in the Vicinity of the BH 3.2.4. CO Line Ratios 3.2.5. Shocked Molecular Clouds from the BH Outflow 3.2.6. The Western Star-forming Region 4. Conclusions References