Draft version August 2, 2021 Typeset using LATEX preprint2 style in AASTeX63 Hot Extended Galaxy Halos Around Local L* Galaxies From Sunyaev-Zeldovich Measurements Joel N. Bregman ,1 Edmund Hodges-Kluck,2 Zhijie Qu,1 Cameron Pratt,1 Jiang-Tao Li,1 and Yansong Yun1 1Department of Astronomy University of Michigan Ann Arbor, MI 48109 USA 2NASA/GSFC Mail Code: 662 Greenbelt , MD 20771 USA ABSTRACT Most of the baryons in L∗ galaxies are unaccounted for and are predicted to lie in 6:5 hot gaseous halos (T ∼ 10 K) that may extend beyond R200. A hot gaseous halo will produce a thermal Sunyaev-Zeldovich signal that is proportional to the product of the gas mass and the mass-weighted temperature. To best detect this signal, we used a Needlet Independent Linear Combination all-sky Planck map that we produced from the most recent Planck data release, also incorporating WMAP data. The sample is 12 L* spiral galaxies with distances of 3 − 10 Mpc, which are spatially resolved so that contamination from the optical galaxy can be excluded. One galaxy, NGC 891, has a particularly strong SZ signal, and when excluding it, the stack of 11 galaxies is detected at about 4σ (declining with radius) and is extended to at least 250 kpc (≈ R200) at > 99% confidence. The gas mass within a spherical volume to a radius of 250 kpc is 10 6 9:8 ± 2:8 × 10 M , for Tavg = 3 × 10 K. This is about 30% of the predicted baryon 11 content of the average galaxy (3:1×10 M ), and about equal to the mass of stars, disk 11 gas, and warm halo gas. The remaining missing baryons (≈ 1:4 × 10 M , 40-50% of the total baryon content) are likely to be hot and extend to the 400 − 500 kpc volume, if not beyond. The result is higher than predictions, but within the uncertainties. Keywords: Spiral galaxies, Sunyaev-Zeldovich effect, Circumgalactic medium arXiv:2107.14281v1 [astro-ph.GA] 29 Jul 2021 1. INTRODUCTION mass of the dark matter halo from the rotation One of the striking realizations of the past curves of galaxies. The rotation curves do not 30 years is that galaxies are missing most of extend to the virial radius (we use R200 hence- their baryons (McGaugh et al. 2010; Dai et al. forth, as it is similar to the virial radius), but a 2010). This comes about because we know the number of rotation curves extend to radii where dark matter dominates the potential. The dark matter distribution is understood from models Corresponding author: Joel N. Bregman [email protected] (Salucci 2019), so we can determine the halo 2 Bregman et al. mass, Mh(R200). The cosmic baryon fraction is remains bound as a hot Circumgalactic Medium 15.7% (Planck Collaboration et al. 2018), cor- (CGM; for spiral galaxies AGN heating is likely responding to a dark matter to baryon ratio of less important). Models are in agreement that 5.37, so the associated baryonic mass is directly the CGM is turbulent and multi-temperature predicted. due to this feedback process. However, mod- An accounting of the mass within the optical els differ significantly as to the mass of the galaxy requires one measure the stellar mass of CGM, its temperature distribution, the frac- stars, plus corrections for no longer visible stel- tion of baryons that lie beyond R200, and the lar remnants (e.g., Loewenstein 2013). To this rate at which gas is cooling onto the disk of the is added the mass of the various ISM phases, galaxy (e.g., Pandya et al. 2020; Davies et al. such as the HI, HII and H2 disks of spirals, or 2020). The reason for these broad uncertainties the hot atmosphere for early-type galaxies. The is that the observational constraints on the hot ratio of the sum of these baryonic components gas properties of the CGM are limited. to the predicted baryonic mass is ∼ 30% for an Gas at Tvirial for a L* spiral has a temperature L* galaxy and can be as little as a few percent typically of 2 − 3 × 106 K, such as in the Milky for less massive galaxies (McGaugh et al. 2010). Way (Henley & Shelton 2012, 2013) or NGC This constitutes the missing baryon problem for 891 (Hodges-Kluck et al. 2018). Some less mas- galaxies, which has motivated a number of ob- sive galaxies, such as NGC 4631 have a lower servational and theoretical studies. temperature halo (Wang et al. 2001; Yamasaki Theoretical models provide context for this et al. 2009), as expected. At these tempera- problem and show that it is not just an account- tures, a significant fraction of the thermal emis- ing curiosity, but a crucial clue to understanding sion falls below the Galactic absorption thresh- galaxy formation and evolution. The simplest old (0.2-0.3 keV), so these halos can be X-ray model for the accretion of baryons were spher- faint. Also, the emission from external galaxies ical in nature, which led to an accretion shock occurs in the energy band where the Galactic at large radius, producing a hot halo at ∼Tvirial X-ray background is highest (mostly due to lo- that might contain significant mass (Mo et al. cal hot gas processes), adding to the challenge 2010). This picture is modified by radiative of observing such objects. losses, leading to significant cooling in the inner Despite these challenges, diffuse X-ray emis- region (50-100 kpc) while the outer part of the sion is commonly detected around nearby edge- halo has a cooling time greater than the Hubble on L* spirals (T¨ullmann et al. 2006; Li & Wang time (Oppenheimer et al. 2018; Hopkins et al. 2013). This X-ray gas is found to extend per- 2018; Nelson et al. 2018; Soko lowska et al. 2018). pendicularly to the disk, and can be traced to These numerical simulations show that this about 5-10 kpc above the plane (e.g., Hodges- model is modified in several major ways for Kluck et al. 2018). Around the Milky Way, such galaxies near L* (e.g., Soko lowska et al. 2018; hot gas is detected in both emission and absorp- Lochhaas et al. 2020; Pandya et al. 2021). The tion, where it is dominated by material within accretion of gas is generally highly asymmetric, 50 kpc, about 0:2R200 (Miller & Bregman 2015). flowing inward along low entropy streams that Another approach to detecting the extended terminate within 0.1R200. The ensuing star for- hot gas halo is through the thermal Sunyaev- mation leads to supernovae that heat their sur- Zel'dovich (SZ) effect(Sunyaev & Zeldovich roundings, driving winds from the galaxies, so 1980), whereby the CMB photons are scattered some of the gas is rendered unbound and some by the hot electrons, thus distorting the spec- Hot Halos in L* Galaxies 3 trum. Previous works have used ensembles of gether, reveal a clear signal to distances of about galaxies to identify a signal, using Planck data. R200 (250 kpc; we note that for our sample, DA The effort by the Planck Collaboration et al. equals the comoving radial distance, to high ac- (2013) and Greco et al.(2015) stacked ∼ 104 curacy). This nearby galaxy sample has a ma- galaxies per bin and detected a signal down to jor advantage in that the SZ signal can be re- about logM∗ = 11.3. These are massive galaxies solved at the resolution of Planck, allowing one (generally early-type) and the stacking proce- to avoid contamination from the optical galaxy dure sometimes includes a signal from the larger and providing spatial information about the ex- scale environment in which they are found (e.g., tent of the hot medium. galaxy clusters). Another stack of slightly less 12:6 13 2. SAMPLE SELECTION massive galaxies (10 . M500 M . 10 M ) by Singh et al.(2018) at typical redshifts of 0.1- We based our target selection on the pre- 0.14 combine X-ray and SZ data and infer that dicted strength of the SZ signal, as given by the 5=3 −2 ∼ 20 − 30% of baryons in these galaxies are hot Mh DA scaling relationship. We chose galax- and produce the observed signal. Their inferred ies close enough that, at Planck resolution, the Y value is less than 2σ significance in most optical galaxy can easily be distinguished from M∗ bins. The study of Ma et al.(2015) cross- the extended SZ halo. The warm dust in the correlated a galaxy sample with weak gravita- optical galaxy can cause a signal that contami- tional lensing (Heymans et al. 2012) with the nates the SZ signal and must be separated. One Planck data and obtain a signal in broad mass can determine Mh with a sufficient amount of 12 14 bins (i.e., 10 − 10 M ) at about 3σ, with a dynamical data, but that does not exist uni- signal extending beyond the virial radius. These formly for each of the target galaxies or can be studies do not provide information about the too uncertain. Instead, we adopted the stellar more ordinary L* spirals, nor do they address mass in the K band as a proxy for Mh (Ktot from the emission region within R200, issues that are 2MASS; Skrutskie et al. 2006). The range in the addressed in this work. K band luminosity, or equivalently, M∗, is mod- The total SZ signal for a single object, Y, est, with a median of logM∗ = 10:83 and where is roughly proportional to MgasTvir, and when the quartile values are at 10.76 and 10.90 (0.07 Mgas / Mh, simple scaling relationships show difference with the median; Table 1).