Galactic Center Research Articles

Multistructured Accretion Flow of Sgr A*. I. Examination of a Radiatively Inefficient Accretion Flow Model

The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test models of radiatively inefficient accretion flows (RIAFs). Previous fits to the quiescent Chandra ACIS-S spectrum found that a RIAF model with an equal inflow–outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution grating spectrum indicate an inflow balanced by an outflow (s ∼ 1) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following 1/r. The data require that the abundance of iron relative to solar is Z Fe < 0.32 Z ⊙ (90% credible interval), much lower than the 2 Z ⊙ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra’s unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic center in X-rays.

Multiwavelength identification of millisecond pulsar candidates in the Galactic bulge

Context. The existence of a population of millisecond pulsars in the Galactic bulge is supported, along with other evidence, by the Fermi GeV excess, an anomalous γ-ray emission detected almost 15 years ago in the direction of the Galactic center. However, radio surveys searching for pulsations have not yet revealed bulge millisecond pulsars. Aims. Identifying promising bulge millisecond pulsar candidates is key to motivating pointed radio pulsation searches. Candidates are often selected among steep-spectrum or polarized radio sources, but multiwavelength information can also be exploited: The aim of this work is to pinpoint strong candidates among the yet unidentified X-ray sources. Methods. We investigated the multiwavelength counterparts of sources detected by the Chandra X-ray observatory that have spectral properties expected for millisecond pulsars in the Galactic bulge. We considered that ultraviolet, optical, and strong infrared counterparts indicate that an X-ray source is not a bulge pulsar, while a radio or a faint infrared counterpart makes it a promising candidate. Results. We identify a large population of more than a thousand X-ray sources without optical, ultraviolet, or strong infrared counterparts. Among them, five are seen for the first time in unpublished radio imaging data from the Very Large Array. We provide the list of promising candidates, for most of which follow-up pulsation searches are ongoing.

The Three Hundred: The existence of massive dark matter-deficient satellite galaxies in cosmological simulations

The observation of a massive galaxy with an extremely low dark matter content (i.e. NGC 1277) has posed questions about how such objects form and evolve in a hierarchical universe. We here report on the finding of several massive, dark matter-deficient galaxies in a set of 324 galaxy clusters theoretically modelled by means of full-physics hydrodynamical simulations. We first focus on two example galaxies selected amongst the most massive and dark matter-deficient ones. By tracing the evolution of these galaxies, we find that their lack of dark matter is a result of multiple pericentre passages. While orbiting their host halo, tidal interactions gradually strip away dark matter while preserving the stellar component. A statistical analysis of all massive satellite galaxies in the simulated clusters shows that the stellar-to-total mass ratio today is strongly influenced by the number of orbits and the distance at pericentres. Galaxies with more orbits and closer pericentres are more dark matter-deficient. Additionally, we find that massive, dark matter-deficient galaxies at the present day are either the remnants of very massive galaxies at infall or former central galaxies of infalling groups. We conclude that such massive yet dark matter-deficient galaxies exist and are natural by-products of typical cluster galaxy evolution, with no specific requirement for an exotic formation scenario.

MICONIC: JWST/MIRI MRS observations of the nuclear and circumnuclear regions of Mrk 231

We present JWST/MIRI MRS spatially resolved ∼5 − 28 μm observations of the central ∼4 − 8 kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk 231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [Mg V] at 5.61 μm or [Ne V] at 14.32 μm) typically associated with the presence of an active galactic nucleus (AGN) are detected in the nuclear region of Mrk 231. This is likely due to the intrinsically X-ray weak nature of its quasar. Some intermediate ionization potential lines, for instance, [Ar III] at 8.99 μm and [S IV] at 10.51 μm, are not detected either, even though they are clearly observed in a star-forming region ∼920 pc south-east of the AGN. Thus, the strong nuclear mid-infrared (mid-IR) continuum is also in part hampering the detection of faint lines in the nuclear region. The nuclear [Ne III]/[Ne II] line ratio is consistent with values observed in star-forming galaxies. Moreover, we resolve for the first time the nuclear starburst in the mid-IR low-excitation line emission (size of ∼400 pc, FWHM). Several pieces of evidence also indicate that it is partly obscured even at these wavelengths. At the AGN position, the ionized and warm molecular gas emission lines have modest widths (W80 ∼ 300 km s−1). There are, however, weak blueshifted wings reaching velocities v02 ≃ − 400 km s−1 in [Ne II]. The nuclear starburst is at the center of a large (∼8 kpc), massive rotating disk with widely-spread, low velocity outflows. Given the high star formation rate of Mrk 231, we speculate that part of the nuclear outflows and the large-scale non-circular motions observed in the mid-IR are driven by its powerful nuclear starburst.

Dispersion and rotation measures from fast radio burst (FRB) host galaxies based on the TNG50 simulation

Context. Fast radio bursts (FRBs) are poised to become important cosmological tools in the near future, as the number of observed FRBs is increasing rapidly with multiple surveys underway. A large sample of FRBs will soon have available dispersion measures (DMs) and rotation measures (RMs), which can be used to study the cosmic baryon density and the intergalactic magnetic field. However, the observed DM and RM of FRBs consists of multiple contributions that must be quantified to estimate the DM and RM of the intergalactic medium (IGM). Aims. In this paper, we estimate one such contribution to DM and RM, namely, of FRB host galaxies. We show how this contribution changes with redshift, galaxy type, and the stellar mass of the galaxies. We also investigate its dependence on galaxy inclination and on an FRB’s offset from the center of the galaxy. Methods. Using the TNG50 simulation of the IllustrisTNG project, we selected 16 500 galaxies at redshifts of 0≤ ɀ ≤2, with stellar masses in the range of 9 ≤ log(M*/M⊙) ≤ 12. In each galaxy, we calculated the DM and RM contributions of 1000 sightlines; from these, we constructed the DM and RM probability density functions (PDFs). Results. We find that the rest frame DM distributions of all galaxies at a given redshift can be fitted by a log normal function and its median and width increase as a function of redshift. The rest-frame RM distribution is symmetric, with a median RMhost,rf=0 rad m–2 and it can be fitted by a combination of a Lorentzian and two Gaussian functions. The redshift evolution of the distribution width can be fitted by a curved power law. The parameters of these functions change for different subsets of galaxies with different stellar mass, inclination, and FRB offset. These changes are due to an increasing ne with redshift, SFR, and stellar mass. We do find a more ordered B field at lower ɀ compared to higher ɀ, as suggested by the presence of more galaxies with B field reversals and B fields dominated by random B field at higher ɀ. Conclusions. We estimated the FRB host DM and RM contributions, which can be used in the future to isolate the IGM contribution from the observed DM and RM of FRBs. We predict that to constrain a σRM,IGM of 2 rad m–2 to the 95% confidence level, we would need to observe 95 000 FRBs at ɀ = 0.5, but only 9 500 FRBs at ɀ = 2.

MeerKAT reveals a ghostly thermal radio ring towards the Galactic Centre

We present the serendipitous discovery of a new radio-continuum ring-like object nicknamed Kýklos (J1802–3353), with MeerKAT UHF and L-band observations. The radio ring, which resembles the recently discovered odd radio circles (ORCs), has a diameter of ∼80″ and is located just ∼6° from the Galactic plane. However, Kýklos exhibits an atypical thermal radio-continuum spectrum (α = −0.1 ± 0.3), which led us to explore different possible formation scenarios. We concluded that a circumstellar shell around an evolved massive star, possibly a Wolf-Rayet, is the most convincing explanation with the present data.

NOEMA formIng Cluster survEy (NICE): Characterizing eight massive galaxy groups at 1.5 &lt; z &lt; 4 in the COSMOS field

The NOrthern Extended Millimeter Array (NOEMA) formIng Cluster survEy (NICE) is a NOEMA large programme targeting 69 massive galaxy group candidates at z &gt; 2 over six deep fields with a total area of 46 deg2. Here we report the spectroscopic confirmation of eight massive galaxy groups at redshifts 1.65 ≤ z ≤ 3.61 in the Cosmic Evolution Survey (COSMOS) field. Homogeneously selected as significant overdensities of red IRAC sources that have red Herschel colours, four groups in this sample are confirmed by CO and [CI] line detections of multiple sources with NOEMA 3 mm observations, three are confirmed with Atacama Large Millimeter Array (ALMA) observations, and one is confirmed by Hα emission from Subaru/FMOS spectroscopy. Using rich ancillary data in the far-infrared and sub-millimetre, we constructed the integrated far-infrared spectral energy distributions for the eight groups, obtaining a total infrared star formation rate (SFR) of 260–1300 M⊙ yr−1. We adopted six methods for estimating the dark matter masses of the eight groups, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass densities. We find that the radial stellar mass densities of the eight groups are consistent with a NFW profile, supporting the idea that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are log(Mh/M⊙) = 12.8 − 13.7 with a general uncertainty of 0.3 dex. Based on the halo mass estimates, we derived baryonic accretion rates (BARs) of (1 − 8)×103 M⊙/yr for this sample. Together with massive groups in the literature, we find a quasi-linear correlation between the integrated SFR/BAR ratio and the theoretical halo mass limit for cold streams, Mstream/Mh, with SFR/BAR = 10−0.46 ± 0.22(Mstream/Mh)0.71 ± 0.16 with a scatter of 0.40 dex. Furthermore, we compared the halo masses and the stellar masses with simulations, and find that the halo masses of all structures are consistent with those of progenitors of Mh(z = 0) &gt; 1014 M⊙ galaxy clusters, and that the most massive central galaxies have stellar masses consistent with those of the brightest cluster galaxy progenitors in the TNG300 simulation. Above all, the results strongly suggest that these massive structures are in the process of forming massive galaxy clusters via baryonic and dark matter accretion.

Different behaviour of the gas-phase and stellar metallicity in the central part of MaNGA galaxies

We quantified the disparity between gas-phase and stellar metallicity in a large galaxy sample obtained from the MaNGA DR17 survey. We found that the gas metallicity is on average closely aligned with the stellar metallicity in the centers of intermediate-mass galaxies. Conversely, the difference is notably larger within the center of massive galaxies. It reaches about −0.18 dex on average for the most massive galaxies, while for low-mass galaxies, the gas metallicity exhibits a slightly lower value than the stellar metallicity. Moreover, the most prominent instances of a reduced gas-phase metallicity in relation to stellar metallicity were observed within the centers of massive red galaxies with low specific star formation rates. Because of the absence of a correlation between the integral mass fraction of neutral gas and the disparity between gas and stellar metallicity, we suggest that the diminished gas-phase metallicity in the centers of massive galaxies might be attributed to the replenishment of gas-depleted central regions through processes such as radial gas flows or accretion from the circumgalactic medium rather than gas infall from the intergalactic medium.

SOFIA/FORCAST Galactic Center Source Catalog

The central regions of the Milky Way constitute a unique laboratory for a wide swath of astrophysical studies; consequently, the inner ∼400 pc have been the target of numerous large surveys at all accessible wavelengths. In this paper, we present a catalog of sources at 25 and 37 μm located within all of the regions observed with the SOFIA/FORCAST instrument in the inner ∼200 pc of the Galaxy. The majority of the observations were obtained as part of the SOFIA Cycle 7 Galactic Center Legacy program survey, which was designed to complement the Spitzer/MIPS 24 μm catalog in regions saturated in the MIPS observations. Due to the wide variety of source types captured by our observations at 25 and 37 μm, we do not limit the FORCAST source catalog to unresolved point sources, or treat all sources as if they are pointlike sources. The catalog includes all detectable sources in the regions, resulting in a catalog of 950 sources, including point sources, compact sources, and extended sources. We also provide the user with metrics to discriminate between the source types.

The Galactic Center in Color: Measuring Extinction with High-proper-motion Stars

The Milky Way’s central parsec is a highly extinguished region with a population of high-proper-motion stars. We have tracked 145 stars for ∼10 yr at wavelengths between 1 and 4 μm to analyze extinction effects in color–magnitude space. Approximately 30% of this sample dims and reddens over the course of years, likely from the motion of sources relative to an inhomogeneous screen of dust. We correct previous measurements of the intrinsic variability fraction for differential extinction effects, resulting in a reduced stellar variability fraction of 34%. The extinction variability subsample shows that the extinguishing material has subarcsecond scales, much smaller variations than previously reported. The observed extinction events imply a typical cross section of 500 au and a density of around 3 × 104 atoms cm−3 for the extinguishing material, which are consistent with measurements of filamentary dust and gas at the Galactic center. Furthermore, given that the stars showing extinction variability tend to be more highly reddened than the rest of the sample, the extinction changes are likely due to material localized to the Galactic center region. We estimate the relative extinction between 1 and 4 μm as AH:AK′:AL′=1.67±0.05:1:0.69±0.03 . Our measurement of extinction at longer wavelengths—L’ (3.8 μm)—is inconsistent with recent estimations of the integrated extinction toward the central parsec. One interpretation of this difference is that the dust variations this experiment are sensitive to—which are local to the Galactic center—are dominated by grains of larger radius than the foreground.

Fermi and eROSITA Bubbles as Persistent Structures of the Milky Way

The Fermi and eROSITA bubbles (FBs and eRBs), large diffuse structures in our Galaxy, can be the by-products of steady star formation activity. To simultaneously explain the star formation history of the Milky Way (MW) and the metallicity of ∼Z ⊙ at the Galactic disk, a steady Galactic wind driven by cosmic rays (CRs) is required. For tenuous gases with a density of ≲10−3 cm−3, CR heating dominates over radiative cooling, and the gas can maintain the virial temperature of ∼0.3 keV, ideal for escape from the Galactic system as the wind. A part of the wind falls back onto the disk like a Galactic fountain flow. We model the wind dynamics according to the Galactic evolution scenario and find that the scale height and surface brightness of the X-ray and the hadronic gamma-ray emissions from such fountain flow region can be consistent with the observed properties of the FBs and eRBs. This implies that the bubbles are persistent structures of the MW existing over (at least) the last ∼1 Gyr rather than evanescent structures formed by nontrivial, ∼10 Myr past Galactic center transient activities.

Do neutrinos bend? Consequences of an ultralight gauge field as dark matter

An ultralight gauge boson could address the missing cosmic dark matter, with its transverse modes contributing to a relevant component of the galactic halo today. We show that, in the presence of a coupling between the gauge boson and neutrinos, these transverse modes affect the propagation of neutrinos in the galactic core. Neutrinos emitted from galactic or extra-galactic supernovae could be delayed by δt=10−8–101s for the gauge boson masses mA′=10−23–10−19eV and the coupling with the neutrino g=10−27–10−20. While we do not focus on a specific formation mechanism for the gauge boson as the dark matter in the early Universe, we comment on some possible realizations. We discuss model-dependent current bounds on the gauge coupling from fifth-force experiments, as well as future explorations involving supernovae neutrinos. We consider the concrete case of the DUNE facility, where the coupling can be tested down to g≃10−27 for neutrinos coming from a supernova event at a distance d=10kpc from Earth.

Imaging fermionic dark matter cores at the centre of galaxies

ABSTRACT Current images of the supermassive black hole (SMBH) candidates at the centre of our Galaxy and M87 have opened an unprecedented era for studying strong gravity and the nature of relativistic sources. Very-long-baseline interferometry data show images consistent with a central SMBH within General Relativity (GR). However, it is essential to consider whether other well-motivated dark compact objects within GR could produce similar images. Recent studies have shown that dark matter (DM) haloes modelled as self-gravitating systems of neutral fermions can harbour very dense fermionic cores at their centres, which can mimic the space–time features of a black hole (BH). Such dense, horizonless DM cores can satisfy the observational constraints: they can be supermassive and compact and lack a hard surface. We investigate whether such cores can produce similar observational signatures to those of BHs when illuminated by an accretion disc. We compute images and spectra of the fermion cores with a general-relativistic ray tracing technique, assuming the radiation originates from standard $\alpha$ discs, which are self-consistently solved within the current DM framework. Our simulated images possess a central brightness depression surrounded by a ring-like feature, resembling what is expected in the BH scenario. For Milky Way-like haloes, the central brightness depressions have diameters down to ${\sim} 35\, \mu \text{as}$ as measured from a distance of approximately $8\,$ kpc. Finally, we show that the DM cores do not possess photon rings, a key difference from the BH paradigm, which could help discriminate between the models.

Testing the Origins of Neutrino Mass with Supernova-Neutrino Time Delay.

The origin of neutrino masses remains unknown. Both the vacuum mass and the dark mass generated by the neutrino interaction with dark matter (DM) particles or fields can fit the current oscillation data. The dark mass squared is proportional to the DM number density and, therefore, varies on the galactic scale with much larger values around the Galactic Center. This affects the group velocity and the arrival time delay of core-collapse supernovae (SN) neutrinos. This time delay, especially for the ν_{e} neutronization peak with a sharp time structure, can be used to distinguish the vacuum and dark neutrino masses. For illustration, we explore the potential of the Deep Underground Neutrino Experiment (DUNE), which is sensitive to ν_{e}. Our simulations show that DUNE can distinguish the two neutrino mass origins at more than 5σ C.L., depending on the observed local value of neutrino mass, the neutrino mass ordering, the DM density profile, and the SN location.

A calibrated model for N-body dynamical friction acting on supermassive black holes

Abstract Black holes are believed to be crucial in regulating star formation in massive galaxies, which makes it essential to faithfully represent the physics of these objects in cosmological hydrodynamics simulations. Limited spatial and mass resolution and the associated discreteness noise make following the dynamics of black holes especially challenging. In particular, dynamical friction, which is responsible for driving massive black holes towards the centres of galaxies, cannot be accurately modelled with softened N-body interactions. A number of subgrid models have been proposed to mimic dynamical friction or directly include its full effects in simulations. Each of these methods has its individual benefits and shortcomings, while all suffer from a common issue of being unable to represent black holes with masses below a few times the simulated dark matter particle mass. In this paper, we propose a correction for unresolved dynamical friction, which has been calibrated on simulations run with the code KETJU, in which gravitational interactions of black holes are not softened. We demonstrate that our correction is able to sink black holes with masses greater than the dark matter particle mass at the correct rate. We show that the impact of stochasticity is significant for low-mass black holes (MBH ≤ 5MDM) and propose a correction for stochastic heating. Combined, this approach is applicable to next generation cosmological hydrodynamics simulations that jointly track galaxy and black hole growth with realistic black hole orbits.

Central molecular zones in galaxies: thirco (6-5) and molecular gas conditions in bright nearby galaxies

This paper summarizes all presently available $J_ upp thirco and accompanying co measurements of galaxy centers including new $J$=6-5 thirco and co observations of eleven galaxies with the Atacama Pathfinder EXperiment (APEX) telescope and also $Herschel$ high-$J$ measurements of both species in five galaxies. The observed $J$=6-5/$J$=1-0 co integrated temperature ratios range from 0.10 to 0.45 in matching beams. Multi-aperture data indicate that the emission of thirco (6-5) is more centrally concentrated than that of co (6-5). The intensities of co (6-5) suggest a correlation with those of $ but not with those of HCN. The new data are essential in refining and constraining the parameters of the observed galaxy center molecular gas in a simple two-phase model to approximate its complex multi-phase structure. In all galaxies except the Seyfert galaxy NGC 1068, high-$J$ emission from the center is dominated by a dense ($n and relatively cool (20-60 K) high-pressure gas. In contrast the low-$J$ lines are dominated by low-pressure gas of a moderate density ($n and more elevated temperature (60-150 K) in most galaxies. The three exceptions with significant high-pressure gas contributions to the low-$J$ emission are all associated with active central star formation.

Mono- and oligochromatic extreme-mass-ratio inspirals

The gravitational capture of a stellar-mass object by a supermassive black hole represents a unique probe of warped spacetime. The small object, typically a stellar-mass black hole, describes a very large number of cycles before crossing the event horizon. Because of the mass difference, we call these captures extreme-mass-ratio inspirals (EMRIs). Their merger event rate at the Galactic Center is negligible, but the amount of time spent in the early inspiral is not. Early EMRIs (E-EMRIs) spend hundreds of thousands of years in band during this phase. At very early stages, the peak of the frequency will not change during an observational time. At later stages in the evolution, it will change a bit and finally the EMRI explores a wide range of them when it is close to merger. We distinguish between “monocromatic” E-EMRIs, which do not change their (peak) frequency, oligochromatic E-EMRIs, which explore a short range and polychromatic ones, and the EMRIs which have been discussed so far in the literature. We derive the number of E-EMRIs at the Galactic Center, and we also calculate their signal-to-noise ratios (SNR) and perform a study of parameter extraction. We show that parameters such as the spin and the mass can be extracted with an error which can be as small as 10−11 and 10−5M⊙. There are between hundreds and thousands of E-EMRIs in their monochromatic stage at the Galactic Centre (GC), and tens in their oligochromatic phase. The SNR ranges from a minimum of 10 (larger likelihood) to a maximum of 106 (smaller likelihood). Moreover, we derive the contribution signal corresponding to the incoherent sum of continuous sources of oligochromatic E-EMRIs with two representatives masses, 10M⊙ and 40M⊙, and show that their curves will cover a significant part of Laser Interferometer Space Antennas (LISAs) sensitivity curve. Depending on their level of circularization, they might be detected as individual sources or form a foreground population. Published by the American Physical Society 2024

Binary Mergers in the Centers of Galaxies: Synergy between Stellar Flybys and Tidal Fields

Galactic centers (GCs) are very dynamically active environments, often harboring a nuclear star cluster and supermassive black hole at their cores. Binaries in these environments are subject to strong tidal fields that can efficiently torque its orbit, exciting near-unity eccentricities that ultimately lead to their merger. In turn, frequent close interactions with passing stars impulsively perturb the orbit of the binary, generally softening their orbit until their evaporation, potentially hindering the role of tides to drive these mergers. In this work, we study the evolution of compact object binaries in the GC and their merger rates, focusing for the first time on the combined effect of the cluster’s tidal field and flyby interactions. We find a significant synergy between both processes, where merger rates increase by a factor of ∼10−30 compared to models in which only flybys or tides are taken into account. This synergy is a consequence of the persistent tides-driven eccentricity excitation that is enhanced by the gradual diffusion of j z driven by flybys. The merger efficiency peaks when the diffusion rate is ∼10−100 slower than the tides-driven torquing. Added to this synergy, we also find that the gradual softening of the binary can lift the relativistic quenching of initially tight binaries, otherwise unable to reach extreme eccentricities, and thus expanding the available phase space for mergers. Cumulatively, we conclude that despite the gradual softening of binaries due to flybys, these greatly enhance their merger rates in GCs by promoting the tidal-field-driven eccentricity excitation.

Observation of the Galactic Center PeVatron beyond 100 TeV with HAWC

We report an observation of ultrahigh-energy (UHE) gamma rays from the Galactic center (GC) region, using 7 yr of data collected by the High-Altitude Water Cherenkov (HAWC) Observatory. The HAWC data are best described as a point-like source (HAWC J1746-2856) with a power-law spectrum ( dN/dE=ϕE/26TeVγ ), where γ = −2.88 ± 0.15stat − 0.1sys and ϕ = 1.5 × 10−15 (TeV cm2 s)−1 ±0.3stat−0.13sys+0.08sys extending from 6 to 114 TeV. We find no evidence of a spectral cutoff up to 100 TeV using HAWC data. Two known point-like gamma-ray sources are spatially coincident with the HAWC gamma-ray excess: Sgr A* (HESS J1745-290) and the Arc (HESS J1746-285). We subtract the known flux contribution of these point sources from the measured flux of HAWC J1746-2856 to exclude their contamination and show that the excess observed by HAWC remains significant (>5σ), with the spectrum extending to >100 TeV. Our result supports that these detected UHE gamma rays can originate via hadronic interaction of PeV cosmic-ray protons with the dense ambient gas and confirms the presence of a proton PeVatron at the GC.

The H i Reservoir in Central Spiral Galaxies and the Implied Star Formation Process ∗ ∗ This is the fourth paper in the “From Halos to Galaxies” series.

The cold interstellar medium (ISM) as the raw material for star formation is critical to understanding galaxy evolution. It is generally understood that galaxies stop making stars when, in one way or another, they run out of gas. However, here we provide evidence that central spiral galaxies remain rich in atomic gas even if their star formation rate (SFR) and molecular gas fraction have dropped significantly compared to “normal” star-forming galaxies of the same mass. Since H i is sensitive to external processes, here we investigate central spiral galaxies using a combined sample from the Sloan Digital Sky Survey, Arecibo Legacy Fast ALFA survey, and the extended GALEX Arecibo SDSS Survey. After proper incompleteness corrections, we find that the key H i scaling relations for central spirals show significant but regular systematic dependence on stellar mass. At any given stellar mass, the H i gas mass fraction is about constant with changing specific star formation rate (sSFR), which suggests that H i reservoir is ubiquitous in central spirals with any star formation status down to M * ∼ 109 M ⊙. Together with the tight correlation between the molecular gas mass fraction and sSFR for galaxies across a wide range of different properties, it suggests that the decline of SFR of all central spirals in the local Universe is due to the halt of H2 supply, though there is plenty of H i gas around. These hence provide critical observations of the dramatically different behavior of the cold multiphase ISM, and a key to understand the star formation process and quenching mechanism.