Search

The He II reionization epoch is expected to take place at $z\sim3-5$. In this stage, the helium and metals in the inter-galactic medium (IGM) are further ionized with additional contributions from harder non-stellar sources, and some large-scale gravitationally bound systems approach virialization. The "Probing the He II re-Ionization ERa via Absorbing C IV Historical Yield (HIERACHY)" program utilizes high- and medium-resolution spectra of bright background quasars at $z\approx3.9-5.2$ to investigate Ly$\alpha$, C IV, and other metal absorption lines during this epoch. Additionally, we employ narrow-band imaging to search for Ly$\alpha$ emitters associated with C IV absorbers, alongside multi-wavelength observations to identify and study particularly intriguing cases. In this paper, we present the design of the HIERACHY program, its current status, major scientific goals, and examples of initial data products from completed Magellan/MIKE, MagE spectroscopy, and MDM imaging observations. We also provide a brief outlook on future multi-wavelength observations that may significantly impact the related science., Comment: 25 pages, 10 figures, 2 tables, accepted for publication by ApJ

The kinematic information of the extraplanar diffuse ionized gas (eDIG) around galaxies provides clues to the origin of the gas. The eDIG-CHANGES project studies the physical and kinematic properties of the eDIG around the CHANG-ES sample of nearby edge-on disk galaxies. We use a novel multi-slit narrow-band spectroscopy technique to obtain the spatial distribution of spectral properties of the ionized gas around NGC 891, which is often regarded as an analog of the Milky Way. We developed specific data reduction procedures for the multi-slit narrow-band spectroscopy data taken with the MDM 2.4m telescope. The data presented in this paper covers the H$\alpha$ and [N II]$\lambda\lambda$6548,6583\AA emission lines. The eDIG traced by the H$\alpha$ and [N II] lines shows an obvious asymmetric morphology, being brighter in the northeastern part of the galactic disk and extending a few kpc above and below the disk. Global variations in the [N II]/H$\alpha$ line ratio suggest additional heating mechanisms for the eDIG at large heights beyond photoionization. We also construct position-velocity (PV) diagrams of the eDIG based on our optical multi-slit spectroscopy data and compare them to similar PV diagrams constructed with the H I data. The dynamics of the two gas phases are generally consistent with each other. Modeling the rotation curves at different heights from the galactic mid-plane suggests a vertical negative gradient in turnover radius and maximum rotation velocity, with magnitudes of approximately $3$ kpc kpc$^{-1}$ and $22-25$ km s$^{-1}$ kpc$^{-1}$, respectively. Measured vertical gradients of the rotation curve parameters suggest significant differential rotation of the ionized gas in the halo, or often referred to as the lagging eDIG. Systematic study of the lagging eDIG in our eDIG-CHANGES project, will help us to better understand the dynamics of the ionized gas in the halo., Comment: 11 pages, 8 figures, accepted for publication in A&A

An Early Enrichment Population (EEP) has been theorized to produce the observed metallicity in the intracluster medium (ICM) of galaxy clusters. This population likely existed at high redshifts (z$\sim$10), relics of which we posit exist today as dwarf galaxies. Previous work argues that the initial mass function (IMF) of the EEP must be flatter than those found at lower redshifts, but with considerable uncertainties. In this work, we present a more quantitative model for the EEP and demonstrate how observational constraints can be applied to the IMF using supernova Type Ia (SNIa) rates, delay time distribution (DTD), and the luminosity function (LF) of galaxy clusters. We determine best-fit values for the slope and mass break of the IMF by comparing IMFs from literature with observed DTDs and the low-luminosity component ($M(R)< -12$) of the Coma LF. We derive two best-fit IMFs: (1) $\alpha_{lo} = -0.13 \pm 0.24$ for $0.07 < M/M_\odot < 1.75$ and $\alpha_{hi} = 0.53 \pm 0.01$ for $1.75 < M/M_\odot < 150$, and (2) $\alpha_{lo} = 1.06 \pm 0.11$ for $0.07 < M/M_\odot < 6$ and $\alpha_{hi} = 0.53 \pm 0.01$ for $6 < M/M_\odot < 150$. We also compare these with sl-5 from Loewenstein (2013) with $\alpha=0.5$ for $0.07 < M/M_\odot < 8$ and $\alpha=0.3$ for $8 < M/M_\odot < 150$. This EEP model, along with stars formed at later times, can produce the observed metallicity, is consistent with other observations, and predicts a significant rise in the SNIa rate at increasing redshift.

No progenitor of a Type Ia supernova is known, but in old population early-type galaxies, one may find SN Ia associated with globular clusters, yielding a population age and metallicity. It also provides insight into the formation path and the SN enhancement rate in globular clusters. We sought to find such associations and identified SN 2019ein to be within the ground-based optical positional uncertainty of a globular cluster candidate within the early-type galaxy NGC 5353 at about 30 Mpc distance. We reduced the positional uncertainties by obtaining Hubble Space Telescope images with the Advanced Camera for Surveys, using filters F475W and F814W and obtained in June 2020. We find that the globular cluster candidate has a magnitude, color, and angular extent that are consistent with it being a typical globular cluster. The separation between the globular cluster and SN 2019ein is 0.43'', or 59 pc in projection. The chance occurrence with a random globular cluster is about 3%, favoring but not proving an association. If the SN progenitor originated in the globular cluster, one scenario is that SN 2019ein was previously a double degenerate white dwarf binary that was dynamically ejected from the globular cluster and exploded within 10 Myr; models do not predict this to be common. Another, but less likely scenario is where the progenitor remained bound to the globular cluster, allowing the double degenerate binary to inspiral on a much longer timescale before producing a SN., Comment: 7 pages with 1 figure

The extraplanar diffuse ionized gas (eDIG) represents ionized gases traced by optical/UV lines beyond the stellar extent of galaxies. We herein introduce a novel multi-slit narrow-band spectroscopy method to conduct spatially resolved spectroscopy of the eDIG around a sample of nearby edge-on disk galaxies (eDIG-CHANGES). In this paper, we introduce the project design and major scientific goals, as well as a pilot study of NGC 3556 (M108). The eDIG is detected to a vertical extent of a few kpc above the disk, comparable to the X-ray and radio images. We do not see significant vertical variation of the [N II]/H$\alpha$ line ratio. A rough examination of the pressure balance between different circum-galactic medium (CGM) phases indicates the magnetic field is in a rough pressure balance with the X-ray emitting hot gas, and may play an important role in the global motion of both the eDIG and the hot gas in the lower halo. At the location of an HST/COS observed UV bright background AGN $\sim29\rm~kpc$ from the center of NGC 3556, the magnetic pressure is much lower than that of the hot gas and the ionized gas traced by UV absorption lines, although the extrapolation of the pressure profiles may cause some biases in this comparison. By comparing the position-velocity diagrams of the optical and CO lines, we also find the dynamics of the two gas phases are consistent with each other, with no evidence of a global inflow/outflow and a maximum rotation velocity of $\sim150\rm~km~s^{-1}$., Comment: 19 pages, 9 figures, 1 table, accepted for publication in ApJ

This paper presents the multi-scale temperature structures in the Milky Way (MW) hot gas, as part of the XMM-Newton Line Emission Analysis Program (X-LEAP), surveying the O VII, O VIII, and Fe-L band emission features in the XMM-Newton archive. In particular, we define two temperature tracers, $I_{\rm OVIII}/I_{\rm OVII}$ (O87) and $I_{\rm FeL}/(I_{\rm OVII}+I_{\rm OVIII})$ (FeO). These two ratios cannot be explained simultaneously using single-temperature collisional ionization models, which indicates the need for multi-temperature structures in hot gas. In addition, we show three large-scale features in the hot gas: the eROSITA bubbles around the Galactic center (GC); the disk; and the halo. In the eROSITA bubbles, the observed line ratios can be explained by a log-normal temperature distribution with a median of $\log T/{\rm K} \approx 6.4$ and a scatter of $\sigma_T \approx 0.2$ dex. Beyond the bubbles, the line ratio dependence on the Galactic latitude suggests higher temperatures around the midplane of the MW disk. The scale height of the temperature variation is estimated to be $\approx$2 kpc assuming an average distance of $5$ kpc for the hot gas. The halo component is characterized by the dependence on the distance to the GC, showing a temperature decline from $\log\,T/{\rm K}\,\approx\, 6.3$ to $5.8$. Furthermore, we extract the auto-correlation and cross-correlation functions to investigate the small-scale structures. O87 and FeO ratios show a consistent auto-correlation scale of $\approx$$ 5^\circ$ (i.e., $\approx$$ 400$ pc at 5 kpc), which is consistent with expected physical sizes of X-ray bubbles associated with star-forming regions or supernova remnants. Finally, we examine the cross-correlation between the hot and UV-detected warm gas, and show an intriguing anti-correlation., Comment: 12 pages, 7 figures. ApJ in press

The XMM-Newton Line Emission Analysis Program (X-LEAP) is designed to study diffuse X-ray emissions from the Milky Way (MW) hot gas, as well as emissions from the foreground solar wind charge exchange (SWCX). This paper reports an all-sky survey of spectral feature intensities corresponding to the O VII, O VIII, and iron L-shell (Fe-L) emissions. These intensities are derived from 5418 selected XMM-Newton observations with long exposure times and minimal contamination from point or extended sources. For 90% of the measured intensities, the values are within $\approx$ 2-18 photons cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (line unit; L.U.), $\approx$ 0-8 L.U., and $\approx$ 0-9 L.U., respectively. We report long-term variations in O VII and O VIII intensities over 22 years, closely correlating with the solar cycle and attributed to SWCX emissions. These variations contribute $\sim30\%$ and $\sim20\%$ to the observed intensities on average and peak at $\approx$ 4 L.U. and $\approx$ 1 L.U. during solar maxima. We also find evidence of short-term and spatial variations in SWCX, indicating the need for a more refined SWCX model in future studies. In addition, we present SWCX- and absorption-corrected all-sky maps for a better view of the MW hot gas emission. These maps show a gradual decrease in oxygen intensity moving away from the Galactic center and a concentration of Fe-L intensity in the Galactic bubbles and disk., Comment: 23 pages, 13 figures, Accepted by ApJS

All-sky maps of the thermal Sunyaev-Zel'dovich effect (SZ) tend to suffer from systematic features arising from the component separation techniques used to extract the signal. In this work, we investigate one of these methods known as needlet internal linear combination (NILC) and test its performance on simulated data. We show that NILC estimates are strongly affected by the choice of the spatial localization parameter ($\Gamma$), which controls a bias-variance trade-off. Typically, NILC extractions assume a fixed value of $\Gamma$ over the entire sky, but we show there exists an optimal $\Gamma$ that depends on the SZ signal strength and local contamination properties. Then we calculate the NILC solutions for multiple values of $\Gamma$ and feed the results into a neural network to predict the SZ signal. This extraction method, which we call Deep-NILC, is tested against a set of validation data, including recovered radial profiles of resolved systems. Our main result is that Deep-NILC offers significant improvements over choosing fixed values of $\Gamma$., Comment: 18 pages, 9 figures, 2 Tables, Accepted by ApJ

Synthesized in the cores of stars and supernovae, most metals disperse over cosmic scales and are ultimately deposited well outside the gravitational potential of their host galaxies. Since their presence is well visible through their X-ray emission lines in the hot gas pervading galaxy clusters, measuring metal abundances in the intracluster medium (ICM) offers us a unique view of chemical enrichment of the Universe as a whole. Despite extraordinary progress in the field thanks to four decades of X-ray spectroscopy using CCD (and gratings) instruments, understanding the precise stellar origins of the bulk of metals, and when the latter were mixed on Mpc scales, requires an X-ray mission capable of spatial, non-dispersive high resolution spectroscopy covering at least the soft X-ray band over a large field of view. In this White Paper, we demonstrate how the Line Emission Mapper (LEM) probe mission concept will revolutionize our current picture of the ICM enrichment. Specifically, we show that LEM will be able to (i) spatially map the distribution of ten key chemical elements out to the virial radius of a nearby relaxed cluster and (ii) measure metal abundances in serendipitously discovered high-redshift protoclusters. Altogether, these key observables will allow us to constrain the chemical history of the largest gravitationally bound structures of the Universe. They will also solve key questions such as the universality of the initial mass function (IMF) and the initial metallicity of the stellar populations producing these metals, as well as the relative contribution of asymptotic giant branch (AGB) stars, core-collapse, and Type Ia supernovae to enrich the cosmic web over Mpc scales. Concrete observing strategies are also briefly discussed., Comment: 19 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity

Molecular gas plays a critical role in explaining the quiescence of star formation (SF) in massive isolated spiral galaxies, which could be a result of either the low molecular gas content and/or the low SF efficiency. We present IRAM 30m observations of the CO lines in the Sombrero galaxy (NGC~4594), the most massive spiral at $d\lesssim30\rm~Mpc$. We detect at least one of the three CO lines covered by our observations in all 13 observed positions located at the galactic nucleus and along a $\sim25\rm~kpc$-diameter dusty ring. The total extrapolated molecular gas mass of the galaxy is $M_{\rm H_2}\approx4\times10^{8}\rm~M_\odot$. The measured maximum CO gas rotation velocity of $\approx379\rm~km~s^{-1}$ suggests that NGC~4594 locates in a dark matter halo with a mass $M_{\rm200}\gtrsim10^{13}\rm~M_\odot$. Comparing to other galaxy samples, NGC~4594 is extremely gas poor and SF inactive, but the SF efficiency is apparently not inconsistent with that predicted by the Kennicutt-Schmidt law, so there is no evidence of enhanced SF quenching in this extremely massive spiral with a huge bulge. We also calculate the predicted gas supply rate from various sources to replenish the cold gas consumed in SF, and find that the galaxy must experienced a starburst stage at high redshift, then the leftover or recycled gas provides SF fuels to maintain the gradual growth of the galactic disk at a gentle rate., Comment: 21 pages, 13 figures, accepted for publication in MNRAS

We introduce the New-ANGELS program, an XMM-Newton survey of $\sim7.2\rm~deg^2$ area around M 31, which aims to study the X-ray populations in M 31 disk and the X-ray emitting hot gas in the inner halo of M 31 up to 30 kpc. In this first paper, we report the catalogue of 4506 detected X-ray sources, and attempt to cross-identify or roughly classify them. We identify 352 single stars in the foreground, 35 globular clusters and 27 supernova remnants associated with M 31, as well as 62 AGNs, 59 galaxies, and 1 galaxy clusters in the background. We uniquely classify 236 foreground stars and 17 supersoft sources based on their X-ray colors. X-ray binaries (83 LMXBs, 1 HMXBs) are classified based on their X-ray colors and X-ray variabilities. The remaining X-ray sources either have too low S/N to calculate their X-ray colors or do not have a unique classification, so are regarded as unclassified. The X-ray source catalogue is published online. Study of the X-ray source populations and the contribution of X-ray sources in the unresolved X-ray emissions based on this catalogue will be published in companion papers., Comment: 30 pages, 12 figures. Accepted for publication in APJS

The extraplanar diffuse ionized gas (eDIG) represents the cool/warm ionized gas reservoir around galaxies. We present a spatial analysis of H$\alpha$ images of 22 nearby edge-on spiral galaxies from the CHANG-ES sample (the eDIG-CHANGES project), taken with the APO 3.5m telescope, in order to study their eDIG. We conduct an exponential fit to the vertical intensity profiles of the sample galaxies, of which 16 can be decomposed into a thin disk plus an extended thick disk component. The median value of the scale height (h) of the extended component is $1.13\pm 0.14$ kpc. We find a tight sublinear correlation between h and the SFR. Moreover, the offset of individual galaxies from the best-fit SFR-h relation shows significant anti-correlation with SFR_SD. This indicates that galaxies with more intense star formation tend to have disproportionately extended eDIG. Combined with data from the literature, we find that the correlations between the eDIG properties and the galaxies' properties extend to broader ranges. We further compare the vertical extension of the eDIG to multi-wavelength measurements of other CGM phases. We find the eDIG to be slightly more extended than the neutral gas (HI 21-cm line), indicating the existence of some extended ionizing sources. Most galaxies have an X-ray scale height smaller than the h, suggesting that the majority of the X-ray emission detected in shallow observations are actually from the thick disk. The h is comparable to the L-band radio continuum scale height, both slightly larger than that at higher frequencies (C-band), where the cooling is stronger and the thermal contribution may be larger. The comparable H$\alpha$ and L-band scale height indicates that the thermal and non-thermal electrons have similar spatial distributions. This further indicates that the thermal gas, the cosmics rays, and the magnetic field may be close to energy equipartition.

The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole feedback and flows of baryonic matter into and out of galaxies. These processes are best studied in X-rays, and emission-line mapping is the pressing need in this area. LEM will use a large microcalorimeter array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft X-ray line emission from objects that constitute galactic ecosystems. These include supernova remnants, star-forming regions, superbubbles, galactic outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs in other galaxies), the Circumgalactic Medium in the Milky Way and other galaxies, and the Intergalactic Medium at the outskirts and beyond the confines of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV band will make it possible to disentangle the faintest emission lines in those objects from the bright Milky Way foreground, providing groundbreaking measurements of the physics of these plasmas, from temperatures, densities, chemical composition to gas dynamics. While LEM's main focus is on galaxy formation, it will provide transformative capability for all classes of astrophysical objects, from the Earth's magnetosphere, planets and comets to the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling gas in galaxy clusters. In addition to pointed observations, LEM will perform a shallow all-sky survey that will dramatically expand the discovery space., Comment: 18 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. v2: All-sky survey figure expanded, references fixed. v3: Added energy resolution measurements for prototype detector array. v4: Author list and reference fixes

The variability of blazars in the X-ray and optical regions informs both the physics of their emitting region and places demands on the observer if a program requires that the object be bright or faint. The extensive simultaneous X-ray and optical observation by Swift provides the best insight into the variable nature of these objects. This program uses \textit{Swift} data for 19 X-ray-bright blazars, generally at $z > 0.1$, to determine their variability properties. The analysis is based on structure functions and provides insight into the nature of the variability and how it depends on time, luminosity, and redshift. We also consider strategies for observing blazars at or above average brightness, given a time delay between planning an observation and obtaining the data. This is critical to observations with orbiting X-ray telescopes, current or future. The variability in the soft X-ray band is typically three to eight times larger than at UV-optical wavelengths, at fixed time differences (i.e., 30 or 100 days). There is almost no difference in the amplitude of variation (X-ray or UV-optical) as a function of redshift (time delay of 30 days) and a modest positive correlation with luminosity. In the X-ray band, blazars that become brighter than normal typically remain bright for at least 2-3 months, although with significant flickering. One can avoid observing objects that are significantly below the average X-ray flux by scheduling the observation when the $F_X > 0.9 F_{X,avg}$, which requires monitoring observations near the time of the scheduling activity., Comment: 24 pages, 10 figures

AGN bubbles could play an important role in accelerating high-energy CRs and galactic feedback. Only in nearby galaxies could we have high enough angular resolution in multi-wavelengths to study the sub-kpc environment of the AGN, where the bubbles are produced and strongly interact with the surrounding ISM. In this paper, we present the latest Chandra observations of the Virgo cluster galaxy NGC 4438, which hosts multi-scale bubbles detected in various bands. The galaxy also has low current star formation activity, so these bubbles are evidently produced by the AGN rather than a starburst. We present spatially resolved spectral analysis of the Chandra data of the $\sim3^{\prime\prime}\times5^{\prime\prime}$ ($\sim200{\rm~pc}\times350\rm~pc$) nuclear bubble of NGC 4438. The power law tail in the X-ray spectra can be most naturally explained as synchrotron emission from high-energy CR leptons. The hot gas temperature increases, while the overall contribution of the non-thermal X-ray emission decreases with the vertical distance from the galactic plane. We calculate the synchrotron cooling timescale of the CR leptons responsible for the non-thermal hard X-ray emission to be only a few tens to a few hundreds of years. The thermal pressure of the hot gas is about three times the magnetic pressure, but the current data cannot rule out the possibility that they are still in pressure balance. The spatially resolved spectroscopy presented in this paper may have important constraints on how the AGN accelerates CRs and drives outflows. We also discover a transient X-ray source only $\sim5^{\prime\prime}$ from the nucleus of NGC 4438. The source was not detected in 2002 and 2008, but became quite X-ray bright in March 2020, with an average 0.5-7 keV luminosity of $\sim10^{39}\rm~ergs~s^{-1}$., Comment: 20 pages, 7 figures, 2 tables, MNRAS in press

Solar wind charge exchange (SWCX) is the primary contamination to soft X-ray emission lines from the Milky Way (MW) hot gas. We report a solar-cycle ($\approx 10$ yr) temporal variation of observed \ion{O}{7} and \ion{O}{8} emission line measurements in the {\it XMM-Newton} archive, which is tightly correlated with the solar cycle traced by the sunspot number (SSN). This temporal variation is expected to be associated with the heliospheric SWCX. Another observed correlation is that higher solar wind (SW) fluxes lead to higher O VII or O VIII fluxes, which is due to the magnetospheric SWCX. We construct an empirical model to reproduce the observed correlation between the line measurements and the solar activity (i.e., the SW flux and the SSN). With this model we discovered a lag of $0.91_{-0.22}^{+0.20}$ yr between the O VII flux and the SSN. This time lag is a combination of the SW transit time within the heliosphere, the lag of the neutral gas distribution responding to solar activity, and the intrinsic lag between the SSN and the launch of a high-energy SW (i.e., $\rm O^{7+}$ and $\rm O^{8+}$). MW O VII and O VIII fluxes have mean values of 5.4 L.U. and 1.7 L.U., which are reduced by $50\%$ and $30\%$, compared to studies where the SWCX contamination is not removed. This correction also changes the determination of the density distribution and the temperature profile of the MW hot gas., Comment: 10 pages, 5 figures. ApJ accepted

Dwarf galaxies are missing nearly all of their baryons and metals from the stellar disk, presumed to be in a bound halo or expelled beyond the virial radius. The virial temperature for galaxies with $M_{\rm h} \sim 10^9 - 10^{10}$ $M_{\odot}$ is similar to the collisional ionization equilibrium temperature for the C IV ion. We searched for UV absorption from C IV in six sightlines toward three dwarf galaxies in the anti-M31 direction and at the periphery of the Local Group ($D \approx$ 1.3 Mpc; Sextans A, Sextans B, and NGC 3109). The C IV doublet is detected in only one of six sightlines, toward Sextans A, with $\log N({\rm C IV})$ = $13.06 \pm 0.08$. This is consistent with our gaseous halo models, where the halo gas mass is determined by the cooling rate, feedback, and the star formation rate; the inclusion of photoionization is an essential ingredient. This model can also reproduce the higher detection rate of O VI absorption in other dwarf samples (beyond the Local Group), and with C IV only detectable within $\sim 0.5R_{\rm vir}$., Comment: Accepted for publication by ApJ, in press

To understand structure formation in the universe and impose stronger constraints on the cluster mass function and cosmological models, it is important to have large galaxy cluster catalogs. The Swift AGN and Cluster Survey is a serendipitous X-ray survey aimed at building a large statistically selected X-ray cluster catalog with 442 cluster candidates in its first release. Our initial SDSS follow-up study confirmed $50\%$ of clusters in the SDSS footprint as z $<$ 0.5 clusters. Here, we present further optical follow-up analysis of 248 (out of 442) cluster candidates from the Swift cluster catalog using multi-band imaging from the MDM $2.4m$ telescope and the Pan-STARRS survey. We report the optical confirmation of 55 clusters with $> 3\sigma$ galaxy overdensities and detectable red sequences in the color-magnitude space. The majority of these confirmed clusters have redshifts z $<$ 0.6. The remaining candidates are potentially higher redshift clusters that are excellent targets for infrared observations. We report the X-ray luminosity and the optical richness for these confirmed clusters. We also discuss the distinction between X-ray and optical observables for the detected and non-detected cluster candidates., Comment: 18 pages; 6 Figures; 1 Table; accepted for publication in ApJ Supplement series

Most baryonic matter in the universe exists in gaseous form and can be found in structures such as galactic halos and the low-density intergalactic medium. proposed-ray spectroscopy missions such as Athena, Arcus, and Lynx will have the capability to identify absorption lines in spectra toward bright active galactic nuclei (AGNs), which can be used as a tool to probe this missing matter. In this study, we examine the optical fields surrounding 15 primary observing targets and identify the foreground galaxies and galaxy groups that are potential hosts of absorption. We record the basic properties of the potential host and their angular and physical separation from the AGN line of sight. This process is done by marking the location of various galaxies and groups in optical images of the field surrounding the target and plotting their angular separation vs. redshift to gauge physical proximity to the background source. We identify the surrounding objects according to those which have measured redshifts and those that require them., Comment: 38 pages, 29 figures

We compose a 265-sight line MW C IV line shape sample using the HST/COS archive, which is complementary to the existing Si IV samples. C IV has a higher ionization potential ($47 - 64$ eV) than Si IV ($33 - 45$ eV), so it also traces warm gas, which is roughly cospatial with Si IV. The spatial density distribution and kinematics of C IV is identical with Si IV within $\approx 2 \sigma$. C IV is more sensitive to the warm gas density distribution at large radii with a higher element abundance. Applying the kinematical model to the C IV sample, we find two possible solutions of the density distribution, which are distinguished by the relative extension along the disk mid-plane and the normal-line direction. Both two solutions can reproduce the existing sample, and suggest a warm gas disk mass of $\log M(M_\odot) \approx 8$ and an upper limit of $\log M(M_\odot) < 9.3$ within 250 kpc, which is consistent with Si IV. There is a decrease of the C IV/Si IV column density ratio from the Galactic center to outskirts by $0.2-0.3$ dex, which may suggest a phase transition or different ionization mechanisms for C IV and Si IV. Also, we find that the difference between C IV and Si IV is an excellent tracer of small-scale features, and we find a typical size of $5^\circ-10^\circ$ for possible turbulence within individual clouds ($\approx 1\rm~kpc$)., Comment: 14 pages, 9 figures, ApJ Accepted

Most of the baryons in L* galaxies are unaccounted for and are predicted to lie in hot gaseous halos (T ~ 3E6 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 4sigma (declining with radius) and is extended to at least 250 kpc (~R_{200}) at > 99% confidence. The gas mass within a spherical volume to a radius of 250 kpc is 9.8 +/- 2.8 E10 Msun, for Tavg = 3E6 K. This is about 30% of the cosmic baryon content of the average galaxy (3.1E11 Msun), and about equal to the mass of stars, disk gas, and warm halo gas. The remaining missing baryons (~ 1.4E11 Msun, 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., Comment: 25 pages, 11 figures, accepted by ApJ (30 Jan. 2022)

Outflows from super-massive black holes (SMBHs) play an important role in the co-evolution of themselves, their host galaxies, and the larger scale environments. Such outflows are often characterized by emission and absorption lines in various bands and in a wide velocity range blueshifted from the systematic redshift of the host quasar. In this paper, we report a strong broad line region (BLR) outflow from the z~4.7 quasar BR 1202-0725 based on the high-resolution optical spectrum taken with the Magellan Inamori Kyocera Echelle (MIKE) spectrograph installed on the 6.5m Magellan/Clay telescope, obtained from the `Probing the He II re-Ionization ERa via Absorbing C IV Historical Yield' (HIERACHY) project. This rest-frame ultraviolet (UV) spectrum is characterized by a few significantly blueshifted broad emission lines from high ions; the most significant one is the C IV line at a velocity of -6500 km/s relative to the H{\alpha} emission line, which is among the highest velocity BLR outflows in observed quasars at z > 4. The measured properties of UV emission lines from different ions, except for O I and Ly{\alpha}, also follow a clear trend that higher ions tend to be broader and outflow at higher average velocities. There are multiple C IV and Si IV absorbing components identified on the blue wings of the corresponding emission lines, which may be produced by either the outflow or the intervening absorbers., Comment: Accepted for publication in the MNRAS. 15 pages, 7 figures, 3 tables

X-ray observations provide a unique probe of the accretion disk corona of supermassive black holes (SMBHs). In this paper, we present a uniform \emph{Chandra} X-ray data analysis of a sample of 152 $z\geq4.5$ quasars. We firmly detect 46 quasars of this sample in 0.5-2~keV above 3~$\sigma$ and calculate the upper limits of the X-ray flux of the remaining. We also estimate the power law photon index of the X-ray spectrum of 31 quasars. 24 of our sample quasars are detected in the FIRST or NVSS radio surveys; all of them are radio-loud. We statistically compare the X-ray properties of our $z\geq4.5$ quasars to other X-ray samples of AGN at different redshifts. The relation between the rest-frame X-ray luminosity and other quasar parameters, such as the bolometric luminosity, UV luminosity, or SMBH mass, show large scatters. These large scatters can be attributed to the narrow luminosity range at the highest redshift, the large measurement error based on relatively poor X-ray data, and the inclusion of radio-loud quasars in the sample. The $L_{\rm X}-L_{\rm UV}$ relationship is significantly sub-linear. We do not find a significant redshift evolution of the $L_{\rm X}-L_{\rm UV}$ relation, expressed either in the slope of this relation, or the departure of individual AGNs from the best-fit $\alpha_{\rm OX}-L_{\rm UV}$ relation ($\Delta\alpha_{\rm OX}$). The median value of the X-ray photon index is $\Gamma\approx1.79$, which does not show redshift evolution from $z=0$ to $z\sim7$. The X-ray and UV properties of the most distant quasars could potentially be used as a standard candle to constrain cosmological models. The large scatter of our sample on the Hubble diagram highlights the importance of future large unbiased deep X-ray and radio surveys in using quasars in cosmological studies., Comment: 26 pages, 13 figures, 3 tables. MNRAS in press. Online data temporarily available at the following link before the formal publication of the paper: https://sites.google.com/a/umich.edu/jiangtaoli2/research/online-data?authuser=0 Please contact the first author if you have troubles downloading the data

Rich and poor galaxy clusters have the same measured halo metallicity, 0.35-0.4 $Z_\odot$, even though they are an order of magnitude apart in stellar fraction, $M_*/M_{gas}$. The measured intracluster medium (ICM) metallicity in high-mass clusters cannot be explained by the visible stellar population as stars typically make up 3-20% of the total baryon mass. The independence of metallicity of $M_*/M_{gas}$ suggests an external and universal source of metals such as an early enrichment population (EEP). Galaxy cluster RX J1416.4+2315, classified as a fosil system, has a stellar fraction of $M_*/M_{gas}=0.054\pm0.018$, and here we improve the halo metallicity determination using archival Chandra and XMM Newton observations. We determine the ICM metallicity of RXJ1416 to be $0.303\pm0.053$ $Z_\odot$ within $0.3

We report a large-scale ($r\approx 20^\circ$) X-ray and Sunyaev-Zeldovich (SZ)-bright diffuse enhancement toward M31, which might be a Local Hot Bridge connecting the Milky Way (MW) with M31. We subtract the Galactic emission from the all-sky O VII and O VIII emission line measurement survey, and find that the emission of these two ions is enhanced within $r\approx20^\circ$ around M31. The mean emission enhancements are $5.6\pm 1.3$ L.U., and $2.8\pm0.6$ L.U. for O VII and O VIII, respectively ($>4\sigma$ for both ions). We also extract the SZ signal around M31, which suggests a surface brightness $y$ of $2-4\times10^{-7}$, an enhancement $>2.5\sigma$ (and a best fit of $5.9\sigma$). These three measurements trace the hot gas with a temperature $\log~T({\rm K})> 6$, showing similar plateau shapes (flat within $\approx15^\circ$, and zero beyond $\approx30^\circ$). A single-phase assumption leads to a temperature of $\log~T({\rm K})=6.34\pm0.03$, which is determined by the O VII/O VIII line ratio. Combining X-ray and SZ measurements, we suggest that this feature is unlikely to be the hot halo around M31 (too massive) or in the MW (too high pressure and X-ray bright). The plateau shape may be explained by a cylinder connecting the MW and M31 (the Local Hot Bridge). We constrain its length to be about 400 kpc, with a radius of 120 kpc, a density of $\approx 2\times10^{-4}-10^{-3} ~\rm cm^{-3}$, and a metallicity of $0.02-0.1~ Z_\odot$. The baryon mass is $\gtrsim10^{11}~M_\odot$, and the oxygen mass is about $\gtrsim10^8~M_\odot$, which contribute to the baryon or metal budget of the Local Group., Comment: Submitted to ApJ

We report Chandra detection of three UV bright radio quiet quasars at $z\gtrsim5$. We have collected a sufficient number of photons to extract an X-ray spectrum of each quasar to measure their basic X-ray properties, such as the X-ray flux, power law photon index ($\Gamma$), and optical-to-X-ray spectral slope ($\alpha_{\rm OX}$). J074749+115352 at $z=5.26$ is the X-ray brightest radio-quiet quasar at $z>5$. It may have a short timescale variation (on a timescale of $\sim3800\rm~s$ in the observer's frame, or $\sim600\rm~s$ in the rest frame) which is however largely embedded in the statistical noise. We extract phase folded spectra of this quasar. There are two distinguishable states: a "high soft" state with an average X-ray flux $\sim2.7$ times of the "low hard" state, and a significantly steeper X-ray spectral slope ($\Gamma=2.40_{-0.32}^{+0.33}$ vs $1.78_{-0.24}^{+0.25}$). We also compare the three quasars detected in this paper to other quasar samples. We find that J074749+115352, with a SMBH mass of $M_{\rm SMBH}\approx1.8\times10^9\rm~M_\odot$ and an Eddington ratio of $\lambda_{\rm Edd}\approx2.3$, is extraordinarily X-ray bright. It has an average $\alpha_{\rm OX}=-1.46\pm0.02$ and a 2-10 keV bolometric correction factor of $L_{\rm bol}/L_{\rm2-10keV}=42.4\pm5.8$, both significantly depart from some well defined scaling relations. We compare $\Gamma$ of the three quasars to other samples at different redshifts, and do not find any significant redshift evolution based on the limited sample of $z>5$ quasars with reliable measurements of the X-ray spectral properties., Comment: 12 pages, 5 figures, 2 tables, accepted for publication by ApJ

Galactic winds are associated with intense star formation and AGNs. Depending on their formation mechanism and velocity they may remove a significant fraction of gas from their host galaxies, thus suppressing star formation, enriching the intergalactic medium, and shaping the circumgalactic gas. However, the long-term evolution of these winds remains mostly unknown. We report the detection of a wind from NGC 3079 to at least 60 kpc from the galaxy. We detect the wind in FUV line emission to 60 kpc (as inferred from the broad FUV filter in GALEX) and in X-rays to at least 30~kpc. The morphology, luminosities, temperatures, and densities indicate that the emission comes from shocked material, and the O/Fe ratio implies that the X-ray emitting gas is enriched by Type II supernovae. If so, the speed inferred from simple shock models is about 500 km/s, which is sufficient to escape the galaxy. However, the inferred kinetic energy in the wind from visible components is substantially smaller than canonical hot superwind models., Comment: 15 pages, 8 figures, 3 tables. Accepted for publication in ApJ

We develop a kinematical model for the Milky Way Si IV-bearing gas to determine its density distribution and kinematics. This model is constrained by a column density line shape sample extracted from the {\it HST}/COS archival data, which contains 186 AGN sight lines. We find that the Si IV ion density distribution is dominated by an extended disk along the $z$-direction (above or below the midplane), i.e., $n(z)=n_0\exp(-(z/z_0)^{0.82})$, where $z_0$ is the scale height of $6.3_{-1.5}^{+1.6}$ kpc (northern hemisphere) and $3.6_{-0.9}^{+1.0}$ kpc (southern hemisphere). The density distribution of the disk in the radial direction shows a sharp edge at $15-20$ kpc given by, $n(r_{\rm XY})=n_0\exp(-(r_{\rm XY}/r_0)^{3.36})$, where $r_0 \approx 12.5\pm0.6$ kpc. The difference of density distributions over $r_{\rm XY}$ and $z$ directions indicates that the warm gas traced by \ion{Si}{4} is mainly associated with disk processes (e.g., feedback or cycling gas) rather than accretion. We estimate the mass of the warm gas (within 50 kpc) is $\log (M(50 {\rm kpc})/M_\odot)\approx8.1$ (assuming $Z\approx0.5Z_\odot$), and a $3\sigma$ upper limit of $\log (M(250 {\rm kpc})/M_\odot)\approx9.1$ (excluding the Magellanic system). Kinematically, the warm gas disk is nearly co-rotating with the stellar disk at $v_{\rm rot}=215\pm3\rm~km~s^{-1}$, which lags the midplane rotation by about $8\rm~km~s^{-1}~kpc^{-1}$ (within 5 kpc). Meanwhile, we note that the warm gas in the northern hemisphere has significant accretion with $v_{\rm acc}$ of $69\pm 7\rm ~km~s^{-1}$ at 10 kpc (an accretion rate of $-0.60_{-0.13}^{+0.11}~M_\odot\rm~yr^{-1}$), while in the southern hemisphere, there is no measurable accretion, with an upper limit of $0.4~M_\odot\rm~yr^{-1}$., Comment: Accepted by ApJ

SZ scaling relations have been used to test the self-similar prediction for massive galaxy clusters, but little attention has been given to individual galaxy groups. We investigate the scaling relations of galaxy groups and clusters near the North Ecliptic Pole using X-ray and SZ observations. This region of the sky is where both the ROSAT and Planck satellites achieved their deepest observations, permitting the investigation of lower mass systems. Our sample consists of 62 X-ray detected groups and clusters, spanning a mass range of $10^{13.4}M_{\odot}<~M_{500}<10^{15}M_{\odot}$ and redshifts of $0.03\lesssim z \lesssim 0.82$. We extract the total SZ flux from unresolved Planck data and estimate the fraction of the SZ flux within $R_{500}$ assuming two different pressure profiles. The SZ scaling relations were derived using a Bayesian technique that accounts for censored data. We find a power law slope of $1.73^{+0.19}_{-0.18}$ for the $Y_{SZ}-M_{500}$ relation which is consistent with the self-similar prediction of 5/3. The slope of $0.89^{+0.09}_{-0.08}$ for the $Y_{SZ}-L_{X,500}$ relation is in agreement with other observational studies but not the self-similar prediction of 5/4, and the $Y_{SZ}-Y_{X}$ relation lies below the 1:1 relation when the slope is fixed to unity. The determined scaling relations are dependent on the selected pressure profile, so resolved data are needed to determine the effects of AGN feedback. In addition, we find a number of potential cluster candidates in the Planck Compton maps that were not identified in our X-ray sample., Comment: 15 pages, 5 figures, 4 tables, Accepted for publication in the Astrophysical Journal

The warm ($\log~T\approx5$) gas is an important gaseous component in the galaxy baryonic cycle. We built a 2-dimension disk-CGM model to study the warm gas distribution of the Milky Way (MW) using the absorption line surveys of Si IV and O VI. In this model, the disk component of both ions has the same density profile ($n(r, z)=n_0\exp(-|z|/z_0)\exp(-r/r_0)$) with a scale height of $z_0=2.6\pm0.4\rm~kpc$ and a scale length of $r_0=6.1\pm1.2\rm~kpc$. For this disk component, we calculate the warm gas mass of $\log(M/M_\odot)=(7.6\pm0.2)-\log\frac{Z}{Z_\odot}$. The similar disk density profiles and total masses of Si IV and O VI-bearing gas set constraints on the ionization mechanisms. We suggest that the warm gas disk might be dominated by the Galactic fountain mechanism, which ejects and recycles gas to set both the scale height and the scale length of the warm gas disk. The CGM component in our model has a dependence on Galactic latitude with a higher column density along the direction perpendicular to the Galactic plane ($b=90^\circ$) than the column density along the radial direction ($b=0^\circ$). The column density difference between these two directions is $0.82\pm0.32\rm~dex$ at $6.3\sigma$ for both ions. This difference may be due to the enrichment of Galactic feedback to the entire CGM, or an additional interaction layer between the warm gas disk and the CGM; existing data cannot distinguish between these two scenarios. If this higher column density at $b=90^\circ$ is for the entire CGM, the total warm CGM mass is $\log(M/M_\odot)\approx(9.5-9.8)-\log\frac{Z}{0.5Z_\odot}$ within the MW virial radius of $250\rm~kpc$., Comment: ApJ in press. 16 pages, 5 figures, 4 tables

We present IRAM 30m observations of molecular lines of CO and its isotopologues from the massive spiral galaxy NGC 5908 selected from the CGM-MASS sample. $^{12}$CO $J=1-0$, $^{12}$CO $J=2-1$, and $^{13}$CO $J=1-0$ lines have been detected in most of the positions along the galactic disk. The total molecular gas mass of NGC 5908 is $\sim7\times10^9\rm~M_\odot$ and the total cool gas mass adding atomic hydrogen is $\sim1.3\times10^{10}\rm~M_\odot$, comparable to the upper limit of the mass of the X-ray emitting hot gas in the halo. Modeling the rotation curves constructed with all three CO lines indicates that NGC 5908 has a dark matter halo mass of $M_{\rm vir}\sim10^{13}\rm~M_{\rm \odot}$, putting it among the most massive isolated spiral galaxies. The $^{12}$CO/$^{13}$CO $J=1-0$, $^{12}$CO $J=2-1$/$J=1-0$ line ratios and the estimated molecular gas temperature all indicate normal but non-negligible star formation in this fairly gas-rich massive isolated spiral galaxy, consistent with the measured star formation intensity and surface densities. The galaxy is probably at an early evolutionary stage after a fast growth stage with mergers and/or starbursts, with plenty of leftover cool gas, relatively high SFR, low hot CGM cooling rate, and low X-ray emissivity., Comment: 13 pages, 8 figures, 2 tables, accepted for publication by AJ

The metallicity of galactic gaseous halos provides insights into accretion and feedback of galaxies. The nearby edge-on galaxy NGC 891 has a multi-component gaseous halo and a background AGN (LQAC 035+042 003) projected 5 kpc above the disk near the minor axis. Against the UV continuum of this AGN, we detect lines from 13 ions associated with NGC 891 in new {\it HST}/COS spectra. Most of the absorption is from the warm ionized gas with $\log T=4.22\pm0.04$, $\log n_{\rm H}=-1.26\pm0.51$, and $\log N_{\rm H}=20.81\pm 0.20$. The metallicity of volatile elements (i.e., C, N, and S) is about half solar ($\rm[X/H] \approx -0.3\pm 0.3$), while Mg, Fe, and Ni show lower metallicities of $\rm[X/H]\approx-0.9$. The absorption system shows the depletion pattern seen for warm Galactic diffuse clouds, which is consistent with a mixture of ejected solar metallicity disk gases and the hot X-ray emitting halo ($Z=0.1-0.2Z_\odot$). The warm ionized gases are about 5 times more massive than the cold \ion{H}{1} emitting gases around the galactic center, which might lead to accretion with a mean rate of $10^2~M_\odot\rm~yr^{-1}$ for a period of time. We also detect low metallicity ($\approx 0.1~Z_\odot$) gases toward LQAC 035+042 003 at $110\rm~km~s^{-1}$ (a high velocity cloud) and toward another sight line (3C 66A; 108 kpc projected from NGC 891) at $30\rm~km~s^{-1}$. This low metallicity material could be the cold mode accretion from IGM or the tidal disruption of satellites in the NGC 891 halo., Comment: 17 pages, 9 figures, accepted by ApJ

A large fraction of the baryons and most of the metals in the Universe are unaccounted for. They likely lie in extended galaxy halos, galaxy groups, and the cosmic web, and measuring their nature is essential to understanding galaxy formation. These environments have virial temperatures >10^5.5 K, so the gas should be visible in X-rays. Here we show the breakthrough capabilities of grating spectrometers to 1) detect these reservoirs of hidden metals and mass, and 2) quantify hot gas flows, turbulence, and rotation around the Milky Way and external galaxies. Grating spectrometers are essential instruments for future X-ray missions, and existing technologies provide 50-1500-fold higher throughput compared to current orbiting instruments., Comment: Submitted to the Astro2020 Decadal (6 pages, 5 figures)

Stellar feedback -- stars regulating further star formation through the injection of energy and momentum into the interstellar medium -- operates through a complex set of processes that originate in star clusters but shape entire galaxies. A mature theory of stellar feedback is essential to a complete theory of star and galaxy formation, but the energy and momentum injected by hot gas into its surroundings remains unclear. With a next-generation X-ray observatory, breakthrough progress can be made through precision measurements of the temperature, density, velocity, and abundances of hot gas on scales of star clusters to galactic superwinds., Comment: white paper submitted to the Astro2020 Decadal Survey

Whenever gas is blown out of a galaxy, chances are that it contains some cosmic dust. This dust is an important part of the metals budget for the circumgalactic and intergalactic media, and traces the outflow and stripping history of the galaxy. The dust is also interesting in its own right, as dust plays an essential role in many astrophysical processes. We have only begun to learn about circumgalactic dust, and in particular we do not know how (and when) it escapes its host galaxy. Here we describe the prospects for measuring the dust mass and properties around many individual galaxies, which will form the basis for understanding how the dust got there., Comment: white paper submitted to the Astro2020 Decadal Survey

We report new Chandra hard X-ray ($>2\rm~keV$) and JVLA C-band observations of the nuclear superbubble of NGC 3079, an analog of the "Fermi bubble" in our Milky Way. We detect extended hard X-ray emission on the SW side of the galactic nucleus with coherent multi-wavelength features in radio, H$\alpha$, and soft X-ray. The hard X-ray feature has a cone shape with possibly a weak cap, forming a bubble-like structure with a diameter of $\sim1.1\rm~kpc$. A similar extended feature, however, is not detected on the NE side, which is brighter in all other wavelengths such as radio, H$\alpha$, and soft X-ray. Scattered photons from the nuclear region or other nearby point-like X-ray bright sources, inverse Compton emission from cosmic ray electrons via interaction with the cosmic microwave background, or any individually faint stellar X-ray source populations, cannot explain the extended hard X-ray emission on the SW side and the strongly NE/SW asymmetry. A synchrotron emission model, plus a thermal component accounting for the excess at $\sim1\rm~keV$, can well characterize the broadband radio/hard X-ray spectra. The broadband synchrotron spectra do not show any significant cutoff, and even possibly slightly flatten at higher energy. This rules out a loss-limited scenario in the acceleration of the cosmic ray electrons in or around this superbubble. As the first detection of kpc-scale extended hard X-ray emission associated with a galactic nuclear superbubble, the spatial and spectral properties of the multi-wavelength emissions indicate that the cosmic ray leptons responsible for the broad-band synchrotron emission from the SW bubble must be accelerated in situ, instead of transported from the nuclear region of the galaxy., Comment: 7 pages, 3 figures, ApJ in press

Galaxies are surrounded by halos of hot gas whose mass and origin remain unknown. One of the most challenging properties to measure is the metallicity, which constrains both of these. We present a measurement of the metallicity around NGC 891, a nearby, edge-on, Milky Way analog. We find that the hot gas is dominated by low metallicity gas near the virial temperature at $kT=0.20\pm0.01$ keV and $Z/Z_{\odot} = 0.14\pm0.03$(stat)$^{+0.08}_{-0.02}$(sys), and that this gas co-exists with hotter ($kT=0.71\pm0.04$ keV) gas that is concentrated near the star-forming regions in the disk. Model choices lead to differences of $\Delta Z/Z_{\odot} \sim 0.05$, and higher $S/N$ observations would be limited by systematic error and plasma emission model or abundance ratio choices. The low metallicity gas is consistent with the inner part of an extended halo accreted from the intergalactic medium, which has been modulated by star formation. However, there is much more cold gas than hot gas around NGC 891, which is difficult to explain in either the accretion or supernova-driven outflow scenarios. We also find a diffuse nonthermal excess centered on the galactic center and extending to 5 kpc above the disk with a 0.3-10 keV $L_X = 3.1\times 10^{39}$ erg s$^{-1}$. This emission is inconsistent with inverse Compton scattering or single-population synchrotron emission, and its origin remains unclear., Comment: 23 pages, 13 figures, 6 tables. Accepted to ApJ

We mapped 12CO J=1-0, 12CO J=2-1, 13CO J=1-0, and 13CO J=2-1 lines toward supernova remnant (SNR) Cassiopeia A with the IRAM 30m telescope. The molecular clouds (MCs) along the line of sight of Cas A do not show optically thin, shock-broadened 12CO lines ($\Delta V \le 7$ km s$^{-1}$ toward Cas A), or high-temperature features from shock heating ($T_k \le 22$ K toward Cas A). Therefore, we suggest that there is no physical evidence to support that the SNR is impacting the molecular gas. All the detected MCs are likely in front of Cas A, as implied by the HCO+ absorption line detected in the same velocity ranges. These MCs contribute H$_2$ column densities of $5\times 10^{21}$ cm$^{-2}$, $5\times 10^{21}$ cm$^{-2}$, and $2\times 10^{21}$ cm$^{-2}$ in the west, south, and center of the SNR, respectively. The 20 K warm gas at $V_{LSR}\sim -47$ km s$^{-1}$ is distributed along a large-scale molecular ridge in the south of Cas A. Part of the gas is projected onto Cas A, providing a foreground H$_2$ mass of $\sim 200 (d/3 kpc)^2$ Msun, consistent with the mass of cold dust (15--20 K; 2--4 Msun) found in front of the SNR. We suggest that the 20 K warm gas is heated by background cosmic rays with an ionization rate of $\zeta({\rm H_2})\sim 2\times 10^{-16}$ s$^{-1}$. The cosmic rays or X-ray emission from Cas A are excluded as the heating sources of the clouds., Comment: 18 pages, 8 figures, 2 tables; Accepted to ApJ

We present ~400ks NuSTAR observations of the northeast (NE) and southwest (SW) non-thermal limbs of the Galactic SNR SN1006. We discovered three sources with X-ray emission detected at >50keV. Two of them are identified as background AGN. We extract the NuSTAR spectra from a few regions along the non-thermal limbs and jointly analyze them with the XMM-Newton spectra and the radio data. The broad-band radio/X-ray spectra can be well described with a synchrotron emission model from a single population of CR electrons with a power law energy distribution and an exponential cutoff. The power law index of the electron particle distribution function (PDF) is ~1.88-1.95 for both the NE and SW limbs, and we do not find significant evidence for a variation of this index at different energy (curvature). There are significant spatial variations of the synchrotron emission parameters. The highest energy electrons are accelerated in regions with the lowest expansion velocity, which is opposite to what has been found in the Tycho's SNR. In addition to a gradual steepening of synchrotron emission from the center of the non-thermal limbs to larger azimuthal angles, we also find that both the emission spectrum and the PDF are significantly flatter in three regions in the SW limb where the shock encounters higher density ambient medium. The NE limb also shows significantly higher cutoff energy in the PDF than the SW limb. By comparing with the roughly symmetric TeV emission and largely asymmetric GeV emission from the two non-thermal limbs, we conclude that the asymmetry in the ambient medium and magnetic fields may have largely modified the acceleration and emission of CR leptons., Comment: 21 pages, 14 figures, 7 tables, ApJ in press

The high ionization-state ions trace the hot gases in the universe, of which gaseous halos around galaxies are a major contributor. Following Qu & Bregman (2018), we calculate the gaseous halo contribution to the observed column density distributions for these ions by convolving the gaseous halo model with the observed stellar mass function. The predicted column density distribution reproduces the general shape of the observed column density distribution -- a broken power law with the break point at $\log N=14.0$ for {\OVI}. Our modeling suggests that the high column density systems originate from galaxies for which the virial temperature matches the temperature of the ionization fraction peak. Specifically, this mass range is $\log M_\star=8.5-10$ for {\OVI}, $\log M_\star=9.5-10.5$ for {\NeVIII}, and higher for higher ionization state ions (assuming $T_{\rm max}=2T_{\rm vir}$). A comparison with the observed {\OVI} column density distribution prefers a large radius model, where the maximum radius is twice the virial radius. This model may be in conflict with the more poorly defined {\NeVIII} column density distribution, suggesting further observations are warranted. The redshift evolution of the high column density systems is dominated by the change of the cosmic star formation rate, which decreases from $z=1.0$ to the local universe. Some differences at lower columns between our models and observations indicate that absorption by the intra-group (cluster) medium and intergalactic medium are also contributors to the total column density distributions., Comment: Match the accepted version

We summarize and reanalyze observations bearing upon missing galactic baryons, where we propose a consistent picture for halo gas in L >~ L* galaxies. The hot X-ray emitting halos are detected to 50-70 kpc, where typically, M_hot(<50 kpc) ~ 5E9 Msun, and with density n \propto r^-3/2. When extrapolated to R200, the gas mass is comparable to the stellar mass, but about half of the baryons are still missing from the hot phase. If extrapolated to 1.9-3 R200, the baryon to dark matter ratio approaches the cosmic value. Significantly flatter density profiles are unlikely for R < 50 kpc and they are disfavored but not ruled out for R > 50 kpc. For the Milky Way, the hot halo metallicity lies in the range 0.3-1 solar for R < 50 kpc. Planck measurements of the thermal Sunyaev-Zeldovich effect toward stacked luminous galaxies (primarily early-type) indicate that most of their baryons are hot, near the virial temperature, and extend beyond R200. This stacked SZ signal is nearly an order of magnitude larger than that inferred from the X-ray observations of individual (mostly spiral) galaxies with M_* > 10^11.3 Msun. This difference suggests that the hot halo properties are distinct for early and late type galaxies, possibly due to different evolutionary histories. For the cooler gas detected in UV absorption line studies, we argue that there are two absorption populations: extended halos; and disks extending to ~50 kpc, containing most of this gas, and with masses a few times lower than the stellar masses. Such extended disks are also seen in 21 cm HI observations and in simulations., Comment: 22 pages, 20 figures, 2 tables, submitted to ApJ

The baryon content around local galaxies is observed to be much less than is needed in Big Bang nucleosynthesis. Simulations indicate that a significant fraction of these "missing baryons" may be stored in a hot tenuous circum-galactic medium (CGM) around massive galaxies extending to or even beyond the virial radius of their dark matter halos. Previous observations in X-ray and Sunyaev-Zel'dovich (SZ) signal claimed that $\sim(1-50)\%$ of the expected baryons are stored in a hot CGM within the virial radius. The large scatter is mainly caused by the very uncertain extrapolation of the hot gas density profile based on the detection in a small radial range (typically within 10\%-20\% of the virial radius). Here we report stacking X-ray observations of six local isolated massive spiral galaxies from the CGM-MASS sample. We find that the mean density profile can be characterized by a single power law out to a galactocentric radius of $\approx 200\rm~kpc$ (or $\approx130\rm~kpc$ above the 1~$\sigma$ background uncertainty), about half the virial radius of the dark matter halo. We can now estimate that the hot CGM within the virial radius accounts for $(8\pm4)\%$ of the baryonic mass expected for the halos. Including the stars, the baryon fraction is $(27\pm16)\%$, or $(39\pm20)\%$ by assuming a flattened density profile at $r\gtrsim130\rm~kpc$. We conclude that the hot baryons within the virial radius of massive galaxy halos are insufficient to explain the "missing baryons"., Comment: 7 pages, 4 figures, 1 table, accepted for publication by ApJL

The galactic gaseous halo is a gas reservoir for the interstellar medium in the galaxy disk, supplying materials for star formation. We developed a gaseous halo model connecting the galaxy disk and the gaseous halo by assuming the star formation rate on the disk is balanced by the radiative cooling rate of the gaseous halo, including stellar feedback. In addition to a single-temperature gaseous halo in collisional ionization equilibrium, we also consider the photoionization effect and a steady-state cooling model. Photoionization is important for modifying the ion distribution in low-mass galaxies and outskirts of massive galaxies due to the low densities. The multi-phase cooling model dominates the region within the cooling radius, where t_cooling=t_Hubble. Our model reproduces most of the observed high ionization state ions for a wide range of galaxy masses (i.e., OVI, OVII, NeVIII, MgX, and OVIII). We find that the OVI column density has a narrow range around ~10^14 cm^-2 for halo masses from M_star ~ 3 * 10^10 Msun to 6*10^12 Msun, which is consistent with some but not all observational studies. For galaxies with halo masses <~ 3 * 10^11 Msun, photoionization produces most of the OVI, while for more massive galaxies, the OVI is from the medium that is cooling from higher temperatures. Fitting the Galactic (Milky-Way) OVII and OVIII suggests a gaseous halo model where the metallicity is ~0.55 Zsun and the gaseous halo has a maximum temperature of ~ 1.9 * 10^6 K. This gaseous halo model does not close the census of baryonic material within R200., Comment: matched with accepted version

Double-lobed radio galaxies (DLRGs) often have radio lobes which subtend an angle of less than 180 degrees, and these bent DLRGs have been shown to associate preferentially with galaxy clusters and groups. In this study, we utilize a catalog of DLRGs in SDSS quasars with radio lobes visible in VLA FIRST 20 cm radio data. We cross-match this catalog against three catalogs of galaxies over the redshift range $0 < z < 0.70$, obtaining 81 tentative matches. We visually examine each match and apply a number of selection criteria, eventually obtaining a sample of 44 securely detected DLRGs which are paired to a nearby massive galaxy, galaxy group, or galaxy cluster. Most of the DLRGs identified in this manner are not central galaxies in the systems to which they are matched. Using this sample, we quantify the projected density of these matches as a function of projected separation from the central galaxy, finding a very steep decrease in matches as the impact parameter increases (for $\Sigma \propto b^{-m}$ we find $m = 2.5^{+0.4}_{-0.3}$) out to $b \sim 2$ Mpc. In addition, we show that the fraction of DLRGs with bent lobes also decreases with radius, so that if we exclude DLRGs associated with the central galaxy in the system the bent fraction is 78\% within 1 Mpc and 56\% within 2 Mpc, compared to just 29\% in the field; these differences are significant at $3.6\sigma$ and $2.8\sigma$ respectively. This behavior is consistent with ram pressure being the mechanism that causes the lobes to bend., Comment: accepted to AJ

We present the analysis of the XMM-Newton data of the Circum-Galactic Medium of MASsive Spirals (CGM-MASS) sample of six extremely massive spiral galaxies in the local Universe. All the CGM-MASS galaxies have diffuse X-ray emission from hot gas detected above the background extending $\sim(30-100)\rm~kpc$ from the galactic center. This doubles the existing detection of such extended hot CGM around massive spiral galaxies. The radial soft X-ray intensity profile of hot gas can be fitted with a $\beta$-function with the slope typically in the range of $\beta=0.35-0.55$. This range, as well as those $\beta$ values measured for other massive spiral galaxies, including the Milky Way (MW), are in general consistent with X-ray luminous elliptical galaxies of similar hot gas luminosity and temperature, and with those predicted from a hydrostatic isothermal gaseous halo. Hot gas in such massive spiral galaxy tends to have temperature comparable to its virial value, indicating the importance of gravitational heating. This is in contrast to lower mass galaxies where hot gas temperature tends to be systematically higher than the virial one. The ratio of the radiative cooling to free fall timescales of hot gas is much larger than the critical value of $\sim10$ throughout the entire halos of all the CGM-MASS galaxies, indicating the inefficiency of gas cooling and precipitation in the CGM. The hot CGM in these massive spiral galaxies is thus most likely in a hydrostatic state, with the feedback material mixed with the CGM, instead of escaping out of the halo or falling back to the disk. We also homogenize and compare the halo X-ray luminosity measured for the CGM-MASS galaxies and other galaxy samples and discuss the "missing" galactic feedback detected in these massive spiral galaxies., Comment: 25 pages, 17 figures, 8 tables, accepted for publication by ApJS

Based on the recently completed {\it Chandra}/ACIS survey of X-ray point sources in nearby galaxies, we study the X-ray luminosity functions (XLFs) for X-ray point sources in different types of galaxies and the statistical properties of ultraluminous X-ray sources (ULXs). Uniform procedures are developed to compute the detection threshold, to estimate the foreground/background contamination, and to calculate the XLFs for individual galaxies and groups of galaxies, resulting in an XLF library for 343 galaxies of different types. With the large number of surveyed galaxies, we have studied the XLFs and ULX properties across different host galaxy types, and confirm with good statistics that the XLF slope flattens from lenticular ($\alpha\sim1.50\pm0.07$) to elliptical ($\sim1.21\pm0.02$), to spirals ($\sim0.80\pm0.02$), to peculiars ($\sim0.55\pm0.30$), and to irregulars ($\sim0.26\pm0.10$). The XLF break dividing the neutron star and black hole binaries is also confirmed, albeit at quite different break luminosities for different types of galaxies. A radial dependency is found for ellipticals, with a flatter XLF slope for sources located between $D_{25}$ and 2$D_{25}$, suggesting the XLF slopes in the outer region of early-type galaxies are dominated by low-mass X-ray binaries in globular clusters. This study shows that the ULX rate in early-type galaxies is $0.24\pm0.05$ ULXs per surveyed galaxy, on a $5\sigma$ confidence level. The XLF for ULXs in late-type galaxies extends smoothly until it drops abruptly around $4\times10^{40}$ erg s$^{-1}$, and this break may suggest a mild boundary between the stellar black hole population possibly including 30 $M_\odot$ black holes with super-Eddington radiation and intermediate mass black holes., Comment: 59 pages, 23 figures, 3 tables. Accepted for publication in ApJ

The hot gas in galaxy halos may account for a significant fraction of missing baryons in galaxies, and some of these gases can be traced by high ionization absorption systems in QSO UV spectra. Using high S/N ratio $ HST$/COS spectra, we discovered a high ionization state system at $z=1.1912$ in the sightline toward \object{LBQS 1435-0134}, and two-components absorption lines are matched for Mg X, Ne VIII, Ne VI, O VI, Ne V, O V, Ne IV, O IV, N IV, O III, and H I. Mg X, detected for the first time ($5.8 \sigma$), is a particularly direct tracer of hot galactic halos, as its peak ion fraction occurs near $10^{6.1}\rm~ K$, about the temperature of a virialized hot galaxy halo of mass $\sim 0.5 M^*$. With Mg X and NeVIII, a photoionization model cannot reproduce the observed column densities with path lengths of galaxy halos. For collisional ionization models, one or two temperature models do not produce acceptable fits, but a three temperature model or a power law model can produce the observed results. In the power law model, ${\rm d}N/{\rm d}T = 10^{4.4\pm 2.2-[Z/X]} T^{1.55\pm 0.41}$ with temperatures in the range $10^{4.39\pm0.13} {\rm~K} < T < 10^{6.04\pm0.05}~ {\rm K}$, the total hydrogen column density is $8.2 \times 10^{19} (0.3Z_{\odot}/Z) \rm~ cm^{-2}$ and the positive power law index indicates most of the mass is at the high temperature end. We suggest that this absorption system is a hot volume-filled galaxy halo rather than interaction layers between the hot halo and cool clouds. The temperature dependence of the column density is likely due to the local mixture of multiple phase gases., Comment: 13 pages, 6 figures, 2 tables, accepted to ApJ

The Circum-Galactic Medium of MASsive Spirals (CGM-MASS) is a project studying the overall content, physical and chemical properties, and spatial distributions of the multi-phase circum-galactic medium (CGM) around a small sample of the most massive isolated spiral galaxies in the local Universe. We introduce the sample and present a detailed case study of the XMM-Newton observation of the hot gas halo of NGC5908. After data calibration, point source removal, and background analysis, we find that the diffuse soft X-ray emission of NGC5908 is significantly more extended than the stellar light in the vertical direction. The 0.5-1.25keV radial intensity profile tracing hot gas emission can be detected above the background out to about $2^\prime$, or $30\rm~kpc$ from the nucleus. The unresolved soft X-ray emission can be characterized with a $\beta$-model with a slope of $\beta\approx0.68$. The unresolved 0.5-2keV luminosity within $1^\prime$ is $6.8\times10^{39}\rm~ergs~s^{-1}$, but young stellar sources may contribute significantly to it. Assuming a metallicity of 0.2solar, an upper limit (without subtracting the very uncertain young stellar contribution) to the mass of hot gas within this radius is $2.3\times10^9\rm~M_\odot$. The cooling radius is $r_{\rm cool}\approx25\rm~kpc$ or $\approx0.06r_{\rm 200}$, within which the hot gas could cool radiatively in less than 10Gyr, and the cooling of hot gas could significantly contribute in replenishing the gas consumed in star formation. The hot gas accounts for $\approx1.9\%$ of the baryon detected within the cooling radius. By comparing NGC5908 to other galaxies, we find that its X-ray luminosity per stellar mass is consistent with lower-mass non-starburst field spiral galaxies. However, a large scatter in hot gas soft X-ray emissivity is indicated for spiral galaxies with $M_*\gtrsim2\times10^{11}\rm~M_\odot$., Comment: 14 pages, 8 figures, 2 tables, accepted for publication in the The Astrophysical Journal