The high energy X-ray probe (HEX-P): Resolving the nature of Sgr A* flares, compact object binaries and diffuse X-ray emission in the Galactic center and beyond

Abstract

HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ( FWHM) and broad spectral coverage (0.2–80 keV) with an effective area far superior to current facilities’ (including XMM-Newton and NuSTAR). These capabilities will enable revolutionary new insights into a variety of important astrophysical problems. We present scientific objectives and simulations of HEX-P observations of the Galactic Center (GC) and Bulge. We demonstrate the unique and powerful capabilities of the HEX-P observatory for studying both X-ray point sources and diffuse X-ray emission. HEX-P will be uniquely equipped to explore a variety of major topics in Galactic astrophysics, allowing us to 1) investigate broad-band properties of X-ray flares emitted from the supermassive black hole (BH) at Sgr A* and probe the associated particle acceleration and emission mechanisms; 2) identify hard X-ray sources detected by NuSTAR and determine X-ray point source populations in different regions and luminosity ranges; 3) determine the distribution of compact object binaries in the nuclear star cluster and the composition of the Galactic Ridge X-ray emission; 4) identify X-ray transients and measure fundamental parameters such as black hole spin; 5) find hidden pulsars in the Galactic Center; 6) search for BH–OB binaries and hard X-ray flares from young stellar objects in young massive clusters; 7) measure white dwarf (WD) masses of magnetic CVs to deepen our understanding of CV evolution and the origin of white dwarf magnetic fields; 8) explore primary particle accelerators in the GC in synergy with future TeV and neutrino observatories; 9) map out cosmic-ray distributions by observing non-thermal X-ray filaments; 10) explore past X-ray outbursts from Sgr A* through X-ray reflection components from giant molecular clouds.

1 Introduction

The Galactic Center (GC) harbors an extremely dense and diverse population of stars, compact objects, X-ray binaries (XRBs), molecular clouds, magnetic filaments, and energetic cosmic-ray accelerators within its central few degrees. Sgr A*, the radiative counterpart of the supermassive black hole (SMBH), plays a crucial role in the gas dynamics and binary formation processes in the central few parsecs, where stars and XRBs cluster (Muno et al., 2003; Mori et al., 2021). Despite its faint quiescent X-ray emission, Sgr A* intermittently accelerates particles and emits daily X-ray flares. In the past, Sgr A* is believed to have been far more active, possibly fueling the Fermi GeV bubbles (Cheng et al., 2011), X-ray chimneys (Ponti et al., 2021), and, as recently as a few hundred years ago, X-ray emission from the GC molecular clouds (Ponti et al., 2010; Clavel et al., 2013). The population of relativistic cosmic rays, including leptons and hadrons, is evident in radio/X-ray filaments and TeV emission from molecular clouds (Zhang et al., 2014; Zhang et al., 2020).

X-ray observations provide valuable insights into compact objects, their binary systems, and particle accelerators in the GC region, since those sources are frequently observed to emit X-rays. However, conducting X-ray studies in this densely populated region, which is located at a large distance of 8 kpc, poses challenges. Furthermore, the GC is filled with numerous point and diffuse X-ray sources, requiring high-angular-resolution X-ray telescopes to resolve them. Over the past two decades, extensive surveys of the GC and bulge regions in the soft X-ray band ( keV) have been carried out by Chandra, XMM-Newton and Suzaku. In a 2° × 0.8° region, Chandra detected X-ray point sources below 8 keV (Wang et al., 2002; Muno et al., 2009; Zhu et al., 2018). Studies of diffuse X-ray emission in the GC and Bulge regions have been conducted by Suzaku and XMM-Newton, revealing extensive hot plasma emission (Koyama et al., 2009; Ponti et al., 2015b; Anastasopoulou et al., 2023). Frequent X-ray monitoring of the GC with Swift-XRT resulted in the detection of X-ray transients in the GC (Degenaar et al., 2015). Follow-up X-ray observations with Chandra, XMM-Newton and NuSTAR enabled the identification of some of those transients (Mori et al., 2013; Mori et al., 2019).

In the GC and bulge regions, soft X-ray data is limited to E > 2 keV due to significant X-ray absorption and scattering. Consequently, the narrow-band (2–10 keV) X-ray data obtained by Chandra and XMM-Newton are often insufficient for identifying source types and accurately measuring spectral parameters (e.g., plasma temperatures and photon indices) due to the parameter degeneracy with neutral hydrogen column density (N_H). In the GC, N_H is typically high (N_H ∼ (5–20) × 10²² cm⁻²) and variable across different regions (Johnson et al., 2009), complicating our ability to constrain other spectral parameters in the soft X-ray band. In the higher energy band above 10 keV, INTEGRAL/IBIS and Suzaku-HXD unveiled diffuse hard X-ray emission and a handful of bright XRBs (Revnivtsev et al., 2004; Yuasa et al., 2008). Subsequently, NuSTAR surveyed of the central 2° × 0.8° region, resolving the hard X-ray emission better with its sub-arcminute angular resolution. However, the NuSTAR GC observations are significantly affected by contamination from stray-light and ghost-ray background photons. As a result, several GC and bulge regions, including the Limiting Window (one of the low extinction windows where Chandra detected hundreds of soft X-ray point sources), Sgr C (one of the GC molecular clouds), and the southern part of the central 2° × 0.8° region, remain unobserved by NuSTAR due to the elevated background levels. While Chandra, XMM-Newton and Suzaku have provided a comprehensive soft X-ray view of the GC region, our broad-band knowledge of the X-ray point and diffuse sources in the GC remains relatively limited.

In this paper, we highlight the unique capabilities of HEX-P to significantly advance our understanding of X-ray-emitting objects within the GC, Bulge and other regions in our Galaxy. Specifically, HEX-P′s exceptional angular resolution and broad-band X-ray coverage extending up to 80 keV are crucial for resolving X-ray sources, characterizing their X-ray emission, and identifying their nature. Given its energy band—extending to 25 keV with 3.5″ angular resolution, LET will be highly advantageous for investigating X-ray sources in crowded fields like the GC. On the other hand, HET will allow us to investigate hard X-ray sources and non-thermal X-ray emission with far better sensitivity than NuSTAR. X-ray spectral and timing properties of transients will be fully characterized by HEX-P thanks to its fast readout time. Both LET and HET will provide excellent Fe line diagnostics at E = 6–7 keV; detecting and resolving highly-ionized Fe lines at 6.7 and 6.9 keV will enable us to discern between thermal and non-thermal X-ray sources. Moreover, a neutral Fe fluorescence line at 6.4 keV can be utilized to constrain X-ray reflection models for molecular clouds and magnetic CVs (mCVs). Additionally, both LET and HET possess high-resolution timing capabilities, making them valuable tools for identifying pulsars and millisecond pulsars (Alford et al., 2023; Jaodand et al., 2023). During the 2030s, HEX-P is poised to make unique contributions to multi-messenger astrophysics of the GC, Bulge and Plane, in synergy with existing and upcoming telescopes across various wavelengths. These include the CTAO for TeV gamma-ray observations, the IceCube gen2 neutrino experiment, radio telescopes such as EHT (Event Horizon telescope), SKA (Square Kilometer Array) and next-generation VLA, as well as infrared telescopes like JWST and GRAVITY.

The paper is organized as follows. §2 describes the current design of the HEX-P telescope and simulation tools. §3 briefly introduces key science cases in the primary GC observation program. More detailed descriptions can be found about Sgr A* flares (§4), X-ray source populations in the GC and Bulge (§5 and diffuse X-ray sources in the GC (§6). §7.1 and §7.2 demonstrate some examples for potential HEX-P GO program observations of star-forming regions and CVs. §8 concludes the paper. Additional figures and tables, including a table defining all acronyms used in this paper, can be found in the Supplementary Material.

2 HEX-P mission design and simulation

The High Energy X-ray Probe (HEX-P; Madsen et al., 2023) is a probe-class mission concept that offers sensitive broad-band X-ray coverage (0.2–80 keV) with exceptional spectral, timing and angular capabilities. It features two high-energy telescopes (HETs) that focus hard X-rays and one low-energy telescope (LET) that focuses lower-energy X-rays.

The LET consists of a segmented mirror assembly coated with Ir on monocrystalline silicon that achieves an angular resolution of 3.5″, and a low-energy DEPFET detector, of the same type as the Wide Field Imager (WFI; Meidinger et al., 2020) onboard Athena (Nandra et al., 2013). It has 512 × 512 pixels that cover a field of view of 11.3′ × 11.3′. The LET has an effective passband of 0.2–25 keV, and a full frame readout time of 2 ms, which can be operated in a 128 and 64 channel window mode for higher count rates to mitigate pile-up with faster readout. Pile-up effects remain below an acceptable limit of ∼1% for fluxes up to ∼100 mCrab in the smallest window configuration. Excising the core of the PSF, a common practice in X-ray astronomy, will allow for observations of brighter sources, with a typical loss of up to ∼60% of the total photon counts.

The HET consists of two co-aligned telescopes and detector modules. The optics are made of Ni-electroformed full shell mirror substrates, leveraging the heritage of XMM-Newton, and coated with Pt/C and W/Si multilayers for an effective passband of 2–80 keV. The high-energy detectors are of the same type as those flown on NuSTAR, and they consist of 16 CZT sensors per focal plane, tiled 4 × 4, for a total of 128 × 128 pixel spanning a field of view of 13.4′ × 13.4′. The HET utilizes the same optics technology as XMM-Newton and the PSF is energy-dependent with an HPD of 10″ at 3 keV, at 20 keV, and increases at higher energies. For the purpose of simulations in this paper, an average HPD of 17″ was used across the entire bandpass.

All simulations presented here were produced with a set of response files that represent the observatory performance based on current best estimates (see Madsen et al., 2023). The effective area is derived from a ray-trace of the mirror design including obscuration by all known structures. The detector responses are based on simulations performed by the respective hardware groups, with an optical blocking filter for the LET and a Be window and thermal insulation for the HET. The LET background was derived from a GEANT4 simulation (Eraerds et al., 2021) of the WFI instrument, and the one for the HET from a GEANT4 simulation of the NuSTAR instrument, both positioned at L1. Throughout the paper, we present our simulation results for HEX-P using the SIXTE (Dauser et al., 2019) and XSPEC toolkits (version 12.13.0 (Arnaud, 1996)). To ensure the most realistic simulation results, we incorporated recent high-resolution X-ray images (mostly from Chandra or other wavelength observations), the best-known spectral information, and theoretical model predictions. Various exposure times have been considered for the feasibility studies presented in the following sections.

3 HEX-P’s Galactic center observation program

Throughout the paper, we present how HEX-P will have a high impact on exploring and understanding a variety of high-energy sources in the GC. With its versatile and well-balanced capabilities, including 1) large effective area in 0.2–80 keV, 2) high energy resolution for resolving atomic lines, 3) low background levels, 4) (LET) and (HET) angular resolutions, and 5) ms timing resolution, HEX-P is an ideal X-ray probe mission for investigating a diverse class of X-ray point and diffuse sources. We have outlined the targets for HEX-P′s primary science program in Table 1 and briefly described them below. Figure 1 displays the survey regions covered by the primary science program, along with Chandra X-ray and MeerKAT radio images of the GC region.

FIGURE 1

1) Broad-band X-ray investigation of Sgr A* X-ray flares (§4): HEX-P observations of Sgr A* X-ray flares will play a unique and impactful role in the multi-messenger astrophysics of studying how the supermassive BH accelerates particles and emit flares in various ways. HEX-P will fully characterize the spectral and timing properties of bright flares thanks to the broadband and continuous coverage without earth occultation. In the 2030s, HEX-P in synergy with other telescopes in the radio (EHT), IR (GRAVITY) and TeV (CTA) bands will provide a wealth of multi-wavelength SED and lightcurve data on Sgr A* flares. These legacy data sets will make a significant impact on understanding the particle acceleration and emission mechanisms of Sgr A* flares.

2) Revealing X-ray point source populations in the GC (§5.1, 5.2, 5.3): The broadband X-ray spectral and timing capabilities of HEX-P will provide a variety of useful diagnostic tools for identifying X-ray point sources and transients in the GC, Bulge, and beyond. These diagnostic tools encompass Fe line analysis, broad-band X-ray spectroscopy, pulsation detection, and X-ray variability study. In the primary science program, HEX-P is expected to detect X-ray sources above 8 keV and identify/classify them with existing or future IR and X-ray observation data. In particular, HEX-P will probe the X-ray luminosity range of L_X ∼ 10³²–10³³ erg s⁻¹, poorly explored by NuSTAR due to the limited sensitivity. This is the regime where a large number of the Chandra sources, mostly composed of compact object binaries, remain unclassified. As demonstrated by NuSTAR, HEX-P will obtain broadband X-ray spectral and timing data of X-ray transients, will be instrumental in detecting pulsations (in the case of transient magnetars) and BH spins (in the case of BH transients) as well as characterizing X-ray properties of very faint X-ray transients (VFXTs).

3) Resolving the central 10 parsec region above 10 keV (§5.4): The nuclear star cluster (NSC), spanning over r ≲ 10 pc around Sgr A*, contains an extremely high concentration of stars and X-ray sources where the central hard X-ray emission (CHXE) was discovered by NuSTAR. HEX-P will be able to resolve the CHXE in the 20–40 keV band with better spatial resolution than NuSTAR, and determine whether the core of the CHXE is cuspy or not. This refined hard X-ray view of the NSC/CHXE will enable testing theoretical models on the formation of X-ray binaries through interactions with Sgr A* BH and other stars. Moreover, HEX-P will obtain the most pristine images of non-thermal X-ray sources in the vicinity of Sgr A* over 40 keV, including the central PWN, X-ray filaments, and unidentified hard X-ray sources.

4) Probing the composition of the Galactic Ridge X-ray emission (§5.5): The Limiting Window, known as one of the low extinction regions where Chandra resolved Fe line emission to a few hundred X-ray point sources, holds the key to unravel the composition of the Galactic Ridge X-ray emission (GRXE). Due to the narrow-band Chandra observations below 8 keV, their true plasma temperatures, possibly varying between different X-ray luminosities (which may reflect different types of CVs in the LW), remain unknown. HEX-P is expected to classify a majority of the Chandra sources by performing X-ray spectral analysis on bright sources individually and analyzing stacked spectral data. Consequently, we will be able to accurately measure the plasma temperatures of the X-ray sources across different X-ray luminosity ranges. Such measurements will lead to identifying which types of X-ray sources predominantly constitute the GRXE.

5) Surveying the central molecular zone (§6.1): Within the central molecular zone (CMZ), a handful of molecular clouds have exhibited X-ray emission. This emission consists of variable X-ray reflections of Sgr A* outbursts in the past and/or steady X-ray emission originating from cosmic-ray interactions with the clouds. In the primary observation program, HEX-P will conduct a survey of the Sgr A complex, Sgr B1, B2 and C, aiming to map their neutral Fe fluorescent line emission and X-ray continuum emission. The initial survey data will play a pivotal role in guiding subsequent HEX-P observations, dissecting the X-ray emission mechanisms, and gaining valuable insights into the past activities of Sgr A* and cosmic-ray distributions around the clouds.

6) Investigating energetic particle accelerators and mapping cosmic-ray distributions (§6.2 and Renolds et al. 2023): As an integral part of the major observation program aimed at exploring Galactic particle accelerators, HEX-P will observe prominent cosmic-ray acceleration sites, including the youngest SNR in our Galaxy (G1.9 + 0.3), the most luminous TeV sources (e.g., PWN G0.9 + 0.1), and star clusters (e.g., Arches) in synergy with future CTAO TeV observatory (Reynolds et al., 2023; Mori et al., 2021). Furthermore, HEX-P will characterize the synchrotron radiation emitted by X-ray filaments in conjunction with high-resolution radio data. Through broadband X-ray spectral data of these filaments, we will create a detailed map of the distribution of TeV–PeV electrons in the GC. Ultimately, HEX-P will explore how and where the relativistic cosmic rays are accelerated and then propagate through the GC region.

4 Sgr A^⋆ flares

Sgr A^⋆ is the radiative counterpart of the supermassive BH at the center of the Milky Way. Sgr A^⋆ has been proposed as one of the energetic particle accelerators powered by accretion flow among other PeVatrons candidates in the GC, such as PWNe and star clusters (HESS Collaboration et al., 2016). However, the quiescent emission from Sgr A^⋆ remains so faint with L_Bol ∼ 10³⁶ erg s⁻¹ and L_X ∼ 2 × 10³³ erg s⁻¹, that it casts doubts about its role as a persistent particle accelerator (Genzel et al., 2010). On the other hand, Sgr A^⋆ emits flares frequently, during which its X-ray luminosity can increase by up to two orders of magnitude (Baganoff et al., 2001; Neilsen et al., 2013; Barrière et al., 2014; Ponti et al., 2015a; Haggard et al., 2019). Similarly, IR flares have been observed daily (Genzel et al., 2003; Ghez et al., 2004). Interestingly, whenever an X-ray flare is observed, an IR counterpart is present, whereas the opposite is not always true. These flares are considered to originate from synchrotron or synchrotron self-Compton scattering radiation emitted from accelerated electrons to high Lorentz factors (γ_e ∼ 10^5–6 in the synchrotron case). Typically, X-ray flares last for a few thousand seconds before they decay via radiative cooling or particle escaping. While simultaneous, multi-wavelength observations of Sgr A^⋆ flares in the IR and X-ray band are important, HEX-P will uniquely explore the particle acceleration and emission mechanisms of the highest energy electrons produced by the supermassive BH at Sgr A^⋆.

Previous Chandra observations of Sgr A^⋆ flares suggest an average rate of 1-2 flares per day (Nowak et al., 2012; Neilsen et al., 2013). With fainter X-ray flares occurring more frequently, the expected number of detectable flares depends on the telescope’s sensitivity. For instance, Ponti et al. (2015a) detected 11 flares with XMM-Newton in 1.6 Ms until 2014, yielding a lower rate of 0.6 flares per day. Note that HEX-P will be more sensitive to detecting fainter X-ray flares than XMM-Newton. It was reported that, since 2014, bright X-ray flares with 2–10 keV fluence above erg cm⁻² have seemed to occur more frequently, at a rate of 0.8 flares per day (Mossoux et al., 2020). Assuming these flare detection rates, 3.5–5 flares will be detected by HEX-P during the 500 ks exposure allocated to HEX-P observations of Sgr A^⋆ in the primary science program. We anticipate more HEX-P observations of Sgr A^⋆ will be carried out through subsequent GO or PI-led GTO programs, jointly with other telescopes. Presently, NuSTAR, XMM-Newton and Chandra observe Sgr A^⋆ each year, in conjunction with EHT and GRAVITY, with a total exposure of ≳ 500 ks. Assuming that similar multi-wavelength observation campaigns will be undertaken in the 2030s, it is projected that HEX-P will detect 18–25 X-ray flares during the first 5 years of operation. The orbit of HEX-P at L1 enables uninterrupted coverage of Sgr A^⋆. This is a tremendous advantage for characterizing the flare lightcurve profiles, whereas previous NuSTAR observations were hampered as parts of the flares were missed due to earth occultation.

The primary objective of investigating synchrotron X-ray emission during Sgr A^⋆ flares is to determine their photon indices and potential spectral cutoffs. X-ray photon indices are directly linked to the intrinsic energy distributions of electrons, while X-ray spectral cutoffs reflect the maximum energy at which electrons are accelerated or a synchrotron cooling break. These X-ray spectral characteristics are likely to vary between different flares, such as faint vs. bright or short vs. long flares, and HEX-P is the most suitable for these measurements due to its broader energy coverage (compared to Chandra and XMM-Newton) and higher sensitivity (compared to NuSTAR). For example, Figure 2 shows simulated HEX-P spectra of a bright X-ray flare from Sgr A^⋆ assuming different spectral shapes. These simulation results represent a significant advance as HEX-P will be able to detect spectral cutoffs from Sgr A* flares, which was not achievable by NuSTAR due to the limited sensitivity. Overall, HEX-P will provide much greater insight into the particle acceleration mechanism near the event horizon.

FIGURE 2

Furthermore, HEX-P will eventually collect a large sample of X-ray flares from Sgr A^⋆  revealing their statistical properties. Based on seven bright flares observed with NuSTAR in ≃ 1 Ms, Zhang et al. (2017) reported a potential but unconfirmed correlation between the photon index and flux data. Figure 2C presents a scattered plot of the photon indices and the flare strengths measured by HEX-P, assuming the same flux distribution observed with NuSTAR and that the potential correlation reported in Zhang et al. (2017) exists. With about five-seven bright flares, HEX-P will be able to establish the correlation at ≳ 3σ, and more significantly from a larger number of flares detectable by HEX-P in the first few years.

4.1 Synergy with radio, IR, and TeV gamma-ray telescopes

Simultaneous X-ray and IR observations have proven to provide a more in-depth investigation of the particle acceleration and cooling mechanisms of Sgr A^⋆, as recently realized by a few such successful observations. For instance, EHT directly images the SMBH and its environment in the sub-millimeter band, potentially capturing hot spots swirling around the event horizon following X-ray flares (Wielgus et al., 2022). In the IR band, GRAVITY has pinpointed Sgr A^⋆ and detected flares with its exceptional sensitivity of the adaptive optics. Since the X-ray and IR emission represent different components of the entire synchrotron radiation emitted from Sgr A^⋆ flares, multi-wavelength SED data will enable not only distinguishing between various emission models but also determining the magnetic field strengths through cooling breaks between the IR and X-ray bands (Dodds-Eden et al., 2009; Ponti et al., 2017).

Previous IR and sub-millimeter observations have revealed that flares originate from an area that extends over one gravitational radius, and they orbit around Schwarzschild radii from the SMBH at a non-Keplerian speed (GRAVITY Collaboration et al., 2018; GRAVITY Collaboration et al., 2020). However, these observations are less sensitive in determining the acceleration mechanism since the electron cooling times in these bands are orders of magnitude longer than those for X-ray-emitting electrons. Therefore, only X-ray observations can be used to identify the central engine of accelerating particles to the highest energies and whether it is powered by magnetic reconnection, shock acceleration, etc. As mentioned earlier, by studying bright flares through HEX-P with simultaneous IR observations, we can directly measure the energy distribution and maximum energy of the accelerated electrons (γ ∼ 10⁴–10⁶ for synchrotron, combined with the magnetic field strength determined by IR and X-ray SED data). Tracking the evolution of the electron spectra during a very bright flare may be feasible as HEX-P will yield ample photon statistics of X-ray flares with continuous coverage. Sensitive X-ray observations with HEX-P can also reveal the formation and evolution of the acceleration source on a timescale comparable with the light crossing time, which is of the order of the Schwarzschild radius. These are unique opportunities that can be pursued by HEX-P but are not feasible with the limited bandpass of Chandra and XMM-Newton.

Recent NuSTAR observations suggested somehow softer photon indices of Sgr A* flares compared with those measured by XMM-Newton and Chandra. This could indicate the presence of a high energy cut-off above 10 keV (GRAVITY Collaboration et al., 2021) or an artifact caused by the high background level in NuSTAR data. In contrast, HEX-P observations of Sgr A* flares over a broader X-ray band and with significantly reduced background contamination provide more direct insights into how particles are accelerated at a few Schwarzschild radii from the SMBH and will help to constrain the acceleration models and their physical parameters. If Sgr A* is indeed accelerating electrons to high Lorentz factors above γ_e ∼ 10⁶, which will be evident from detecting hard X-rays with HEX-P, these electrons will also emit TeV gamma-rays through the inverse Compton scattering process. Thus, hard X-ray detections by HEX-P will motivate simultaneous TeV observations with the upcoming CTAO mission and ultimately determine whether the supermassive BH at Sgr A* is a PeVatron accelerator at the present time.

5 X-ray source population in the GC and bulge

Extensive Chandra surveys have detected over 10,000 X-ray sources in the 2° × 0.8° GC field (Wang et al., 2002; Muno et al., 2009; Zhu et al., 2018). Even excluding foreground sources, which constitute about 25% of the total, these Chandra sources represent the densest X-ray source population in our Galaxy. This is illustrated in Figure 3, which shows the Chandra X-ray flux distributions in different regions of the GC and Bulge. In the GC and bulge regions, where optical, UV, and soft X-ray (E ≲ 2 keV) observations are hindered due to significant extinction and dust scattering, the hard X-ray band provides a singular opportunity to explore the populations of compact objects and their binaries. Thus far, multi-decade X-ray observations of the GC, Bulge, and Ridge have revealed that 1) the X-ray point sources are predominantly mCVs, with a modest fraction of XRBs (Xu et al., 2019a); 2) a small population of non-thermal X-ray sources exist, and they are likely LMXBs or runaway pulsars (Hailey et al., 2018); 3) X-ray transients have been detected, a majority of them in the central few pc (Mori et al., 2021); 4) diffuse hard X-ray emission is largely composed of unresolved CVs (Hailey et al., 2016). As described briefly below, these results indicate distinct spatial distributions and compositions of X-ray sources in the GC and Bulge. Identifying these X-ray sources and the underlying X-ray populations is important for testing the fundamental theories and recent N-body simulation results of XRB/CV formation in the vicinity of Sgr A*, the NSC and beyond (Generozov et al., 2018; Szölgyén and Kocsis, 2018; Panamarev et al., 2019), with implications ranging from the rates of gravitational wave events to the nature of dark matter.

FIGURE 3

1) Magnetic CV population in the GC and Bulge: Stacked X-ray spectra of the Chandra sources were well fit with an optically thin thermal plasma model with kT ≳ 8 keV, indicating that they are predominantly mCVs (Hong et al., 2016). MCVs possess WD magnetic fields strong enough to truncate or quench accretion disks, and they are generally classified into two classes: intermediate polars (IPs) and polars. An IP is a type of mCV that has non-synchronized orbits and WD magnetic fields (B ∼ 0.1–10 MG) strong enough to truncate the inner accretion disk; they are copious emitters of hard X-rays (kT ∼ 20–40 keV). On the other hand, polars are magnetically synchronized mCVs with higher B-field strength (typically B ∼ 10–240 MG) and thus possess softer X-ray spectra (kT ∼ 10–20 keV) than IPs due to faster cyclotron cooling. By contrast, non-magnetic CVs (nmCVs) have lower B-fields (B ≲ 0.1 MG), and their accretion disks extend to the WD surface. NmCVs exhibit softer and fainter X-ray emission (kT ≲ 10 keV) than those of mCVs (Byckling et al., 2010). Since only a few % of the Chandra sources in the GC have been detected by NuSTAR above 10 keV (Hong et al., 2016), their source types remain elusive. To distinguish between nmCVs, IPs, and polars, it is important to characterize both broad-band X-ray continuum and Fe emission lines at 6–7 keV (Mukai, 2017; Xu et al., 2019b).

2) X-ray binaries and pulsars: A fraction of the Chandra sources are characterized by non-thermal X-ray spectra with no or weak Fe lines. These non-thermal X-ray sources are likely quiescent LMXBs (containing NSs or BHs) or pulsars (Mori et al., 2021). Subsequently, multi-decade X-ray observations of the Sgr A* region with Chandra, XMM-Newton, Swift and NuSTAR revealed a cusp of LMXBs in the central pc (Hailey et al., 2018; Mori et al., 2021). Finding a concentration of pulsars in the Galactic Bulge will have significant implications for the origin of the GeV gamma-ray excess emission detected by Fermi-LAT (Dexter and O’Leary, 2014).

3) X-ray transients: Daily Swift/XRT monitoring since 2006, combined with other frequent X-ray observations of the GC, resulted in the detection of 20 X-ray transients in the central pc (Muno et al., 2005; Degenaar et al., 2015). Six of the GC transients were identified as NS-LMXBs through the detection of type I X-ray bursts or pulsations (Degenaar et al., 2012), while NuSTAR observations identified two recent Swift transients as outbursting BH-LMXBs (Mori et al., 2019). Otherwise, VFXTs have been detected with peak luminosities of L_X ∼ 10³⁴–10³⁶ erg s⁻¹ (Degenaar et al., 2015) but their identity remains elusive despite extensive investigations by Swift-XRT (Degenaar et al., 2015).

4) Diffuse hard X-ray emission: In the Bulge, Ridge and central 10 pc regions, diffuse hard X-ray emission with different plasma temperatures ranging from kT ∼ 8 to 35 keV have been detected by Suzaku, INTEGRAL/IBIS and NuSTAR (Türler et al., 2010; Yuasa et al., 2012; Perez et al., 2019; 2015). This diffuse X-ray emission is believed to originate from unresolved populations of X-ray point sources. Only a few low-extinction regions in the bulge and ridge, such as the Limiting Window (LW), have been observed by Chandra and resolved into hundreds of X-ray point sources, most of which remain unidentified due to the lack of hard X-ray observations.

As we present below, HEX-P will greatly enhance our understanding of X-ray source populations in the GC and Bulge, covering a wide range of X-ray luminosity (L_X). HEX-P will conduct deep observations of prominent regions such as the Sgr A complex, the Arches cluster, Sgr B2, Sgr C, and the LW. These investigations will fully utilize the broadband spectral and timing capabilities of HEX-P to reveal the X-ray source populations within those regions. Our primary focus will be on classifying X-ray source populations in the range of L_X ∼ 10³²–10³³ erg s⁻¹ at the GC distance (8 kpc), and examining how they vary with L_X and distance from the GC. This unexplored parameter space of X-ray sources holds considerable interest, as we expect to unveil a number of mCVs, XRBs, and pulsars that have not yet been discovered in the GC. In contrast, a significant portion of fainter X-ray sources is likely background AGNs and coronally active dwarfs (Ebisawa et al., 2005; DeWitt et al., 2010; Morihana et al., 2022), which are unrelated to the compact object populations in the GC and Bulge.

5.1 X-ray source population

The NuSTAR GC survey covered of the central 2° × 0.8° region, resulting in detecting a total of 77 hard X-ray sources above 10 keV (Hong et al., 2016). These sources exhibit either thermal X-ray spectra with kT > 20 keV or non-thermal X-ray spectra with Γ = 1.5–2, indicating that they are mCVs or XRBs, respectively (Hong et al., 2016). The NuSTAR data analysis of the GC sources has been largely limited below keV due to high stray-light background levels (Mori et al., 2015). Consequently, despite extensive studies combining NuSTAR observations with archived Chandra and XMM-Newton data, the majority of individual NuSTAR sources remain unidentified because their existing X-ray spectra do not allow us to distinguish between thermal (mCVs) or non-thermal sources (XRBs), and lack accurate measurements of plasma temperatures or photon indices. Furthermore, of the 2° × 0.8° GC region, as well as the LW, were not surveyed by NuSTAR due to substantial background contamination. In addition to identifying the known hard X-ray sources detected by NuSTAR, HEX-P will survey the unexplored regions, allowing us to detect new hard X-ray sources and determine their source types.

The goals and capabilities of HEX-P in classifying and identifying X-ray sources in the GC are specified in Supplementary Table S2 in the Supplementary Material section and in the numbered paragraphs below. The primary GC survey program involves surveying the central 1° × 0.7° region with multiple tiled observations with a typical exposure time of ks per field. Each tile will have a 50 ks exposure, likely divided into multiple, shorter exposures for obtaining X-ray variability data. These tiles will overlap with each other to achieve ks exposure at a given location in the survey area. Additionally, HEX-P plans to conduct deeper observations of Sgr A*, Sgr A, and the LW (Table 1).

We conducted extensive simulations based on existing data to demonstrate the remarkable capabilities of our proposed surveys. In order to determine the detectability of X-ray sources with HEX-P, we adopted the X-ray flux data of 6,760 Chandra sources from Muno et al. (2009), excluding foreground sources. We assumed that these sources have kT = 15–30 keV, which are typical plasma temperatures for mCVs. Assuming a 100 ks exposure per field and kT = 15 keV, HET will be able to detect sources at significance; with a temperature of kT = 30 keV, that number rises to sources within the 2° × 0.8° GC field (Figure 4). HET will be able to extract X-ray photometry data cleanly from these sources since they are separated from each other by at least 8.5″ (half the HET HPD), allowing us to obtain broadband logN-logS distribution, hardness ratio, and variability data in the 2–40 keV band. This is a significant improvement compared to the 77 hard X-ray sources detected by NuSTAR. As described below, HEX-P will be able to identify and classify X-ray sources in the GC using broadband X-ray spectral and timing data, combined with the archived Chandra and IR data. Table 2 (in Section 10) lists various types of X-ray sources that can be identified by extensive HEX-P observations in the GC and other regions in our Galaxy.

FIGURE 4

1) Identification of Chandra and IR counterparts: 18 out of the 77 NuSTAR sources do not have clear associations with Chandra sources because the limited angular resolution of NuSTAR (18″ FWHM) hinders their correlation with Chandra positions, as well as their IR counterparts. HEX-P′s angular resolution, superior to NuSTAR′s, will allow us to identify Chandra counterparts to hard X-ray sources with high confidence. Consequently, Chandra′s sub-arcsecond angular resolution will enable us to find the NIR counterparts to those X-ray sources and determine their nature using the combined X-ray and NIR spectral and temporal properties of those systems (Table 2). Accreting binaries are expected to display ellipsoidal modulations in their NIR light curves, which result from the distortion of Roche-lobe-filling donor stars due to the tidal forces exerted by the compact object. Currently ongoing—and upcoming—JWST surveys will be able to detect such periodic variability even from low-mass companions, thanks to their unprecedented depth (Schoedel et al., 2023). Combining deep NIR and X-ray data will open a new window into the nature and evolution of accreting compact objects. Furthermore, unambiguous Chandra counterpart identification will enable joint X-ray spectral and variability analysis with archived Chandra and XMM-Newton data, allowing us to understand the temporal properties of those sources over a timescale of decades.

2) X-ray spectral analysis: Hard X-ray detections above 10 keV suggest the NuSTAR sources could be either mCVs or XRBs (Hong et al., 2016). To accurately identify these sources, Fe line diagnostics are crucial. Typically, mCVs exhibit thermal X-ray spectra with strong Fe emission lines at 6–7 keV (Xu et al., 2016), while quiescent XRBs show non-thermal X-ray spectra with weak or no Fe lines. However, the majority of the NuSTAR sources lack clear Fe line determinations because their NuSTAR spectra are largely contaminated by diffuse X-ray emission, which also exhibits strong Fe lines. With both the LET and HET covering the 6–7 keV band, HEX-P will greatly improve Fe line diagnostics thanks to its superior angular resolution and significantly reduced background levels compared to NuSTAR. Figure 5 (top) illustrates how HEX-P can discern between thermal and non-thermal X-ray models by resolving the Fe line complex region at 6–7 keV. For faint X-ray sources that lack sufficient source counts for individual spectral analysis, HEX-P will be highly effective in classifying source types via hardness ratios and stacked X-ray spectra. Given its broad energy band, HEX-P is less susceptible to X-ray absorption below a few keV, which can cause parameter degeneracy with intrinsic spectral parameters like kT and Γ. Unlike the limited bandpass of Chandra and XMM-Newton data, HEX-P can investigate the broad-band hardness or quantile ratios, which are more intrinsic to source types and properties. For instance, one can define hardness ratios between 3–10 and 10–40 keV bands, which are nearly unaffected by X-ray absorption and dust scattering effects.

FIGURE 5

3) Detection of X-ray variability and periodic signals: Multi-epoch HEX-P observations of the GC regions can provide long-term X-ray variability data over the course of days and weeks. Year-scale X-ray variabilities can be investigated by combining HEX-P with archived X-ray data obtained by NuSTAR, Chandra  and XMM-Newtonover the past decades. These long-term X-ray variability data can be useful in identifying mCV types. For example: compared to IPs, polars often display significant flux variations (by more than an order of magnitude) in the optical and X-ray band due to higher variability in mass accretion rates. In the GC and bulge regions, optical observations are hindered by significant extinction and dust scattering, whereas IR emission of mCVs is usually dominated by companion stars. Chandra observations have shown promise in detecting orbital periods from a handful of mCVs in the LW (Hong et al., 2012; Bao and Li, 2020), which HEX-P could also achieve by using several sophisticated timing methods with only hundreds of source counts. Based on the simulation-based detection rate from Bao and Li (2020), in the case of a sinusoidal light curve, the detection completeness with significance over 99% by GL algorithm (Gregory and Loredo, 1992) was generally ≳ 20% and 50% for source counts ≳ 100 and 300, respectively. Moreover, at high amplitude variation (F_max/F_min ≳ 4), which is typical for polars, the detection completeness will be %.

5.2 X-ray transients

The GC region hosts a dense population of XRBs and X-ray transients, as revealed by deep X-ray surveys and long-term monitoring over the last two decades (Muno et al., 2005; Degenaar and Wijnands, 2010; Hailey et al., 2018; Mori et al., 2021). In the crowded GC region, HEX-P is uniquely suited to make these requisite hard X-ray ( keV) observations due to its superior angular resolution. HEX-P will be sensitive in obtaining broad-band X-ray spectral and timing data of VFXTs, which belong to a unique class of X-ray transients with peak luminosities of L_X ∼ 10³⁴–10³⁶ erg s⁻¹ (Degenaar et al., 2015). Given its combination of broad energy band and fast readout time, HEX-P is ideal for characterizing bright X-ray transients through spectral and timing analysis without significant deadtime effects. Some distinct features often observed in BH and NS transients, such as thermal disk emission, Comptonization, and relativistic Fe lines, will be fully characterized in HEX-P broad-band X-ray spectral data. Fitting the broad-band HEX-P spectra with relativistic X-ray reflection models will yield the spin measurements of a BH transient within 100 pc of the GC as demonstrated by NuSTAR observations of two BH transients in 2016 (Mori et al., 2019). Given the high concentration of XRBs in the GC, a combination of frequent X-ray monitoring (e.g., Star-X or Swift-XRT) and HEX-P ToO observations will provide an efficient way to measure BH spins from new stellar-mass BH transients (Connors et al. submitted to FrASS). In addition, HEX-P has superb capabilities of detecting pulsations and QPO signals in a wide frequency range from mHz to kHz, thanks to the drastic reduction of background noise which results in an increase of fractional r.m.s. Variability (e.g., Bachetti et al., this volume; Alford et al., this volume).

5.3 Pulsars

In 2009, Fermi-LAT detected an excess of GeV gamma-ray emission extending over the central of the Milky Way Galaxy (Goodenough and Hooper, 2009). The nature of the GeV excess still remains one of the great puzzles in astrophysics, with often-conflicting studies that suggest either a dark matter origin or a large population of unresolved pulsars (Hooper and Linden, 2011; Brandt and Kocsis, 2015; Gautam et al., 2022). In addition, recent N-body simulations predict run-away NSs from the NSC due to their natal kicks or interactions with the supermassive BH and surrounding stars (Bortolas et al., 2017). However, despite extensive searches, no such population of pulsars or milli-second pulsars has ever been observed, resulting in the “missing pulsar problem”. As demonstrated by the discovery of a transient magnetar in the vicinity of Sgr A* (Mori et al., 2013), finding even a few pulsars in the GC and bulge regions will have profound implications for the underlying NS population and test the hypothesis of NS ejection from the NSC (Dexter and O’Leary, 2014; Bortolas et al., 2017). However, it is extremely challenging to detect radio pulsars in the GC due to large dispersion measures until SKA becomes fully operational (Keane et al., 2015). With the sub-millisecond timing resolution, HEX-P is uniquely suited for finding pulsars in the X-ray band and complementary to future radio pulsar searches in the GC and Bulge.

There are a handful of pulsar candidates such as non-thermal, hard X-ray sources detected by Chandra and NuSTAR. For instance, a runaway NS called the “Cannonball”, detected in the radio and X-ray bands (Zhao et al., 2013), exhibited non-thermal X-ray emission extending up to 30 keV (Nynka et al., 2013). Furthermore, two PWNe (G359.95-0.04 and G0.13-0.11) are associated with TeV sources (Mori et al., 2015; Zhang et al., 2020). While Chandra has resolved PWNe X-ray point sources within both PWNe, no pulsations have been detected yet. Since the search for radio pulsations is severely hampered by significant dispersion measures, X-ray telescopes equipped with high angular and timing resolution, such as HEX-P, would be an ideal pulsar finder in the GC and Bulge. Figure 6 illustrates the sensitivity of HEX-P in detecting a  Hz pulsation from G359.95-0.04, located away from Sgr A*. Higher or lower pulse frequencies would change this plot only slightly, due to a weak dependence of significance on the number of trials used in a blind search^1,2 Compared to XMM-Newton EPIC-PN and NuSTAR, HEX-P will be far more sensitive in detecting pulsars due to its msec temporal resolution, high angular resolution, and hard X-ray coverage ( 10 keV), where pulsar emission is more dominant over the PWN and thermal diffuse X-ray emission. During the planned Sgr A* observations (500 ks total), both LET and HET will be able to detect a pulsation signal at level if the pulsed fraction is higher than %. For other pulsar candidates, such as the Cannonball and G0.13-0.11, we anticipate that pulsations will be more easily detectable due to lower background rates from the surrounding regions. We also note that HEX-P, being in a Lagrangian point, will also guarantee much longer stretches of uninterrupted observations with little or no periods of occultation or flaring, which will improve stability for the study of slow variability and the detection of the slowest pulsars (e.g. RCW 103 with a 6.7-h pulsation (De Luca et al., 2006)).

FIGURE 6

5.4 Probing the distribution of compact object binaries in the nuclear star cluster

The central 10 pc region of our Galaxy contains a dense concentration of stars and compact objects around the SMBH at Sgr A*. Within the NSC, which boasts the densest stellar environment in our Galaxy, Chandra detected a few hundred X-ray point sources (Zhu et al., 2018). Above 20 keV, NuSTAR discovered diffuse hard X-ray emission spanning an area of approximately 8 pc × 4 pc, which coincides with the NSC (Perez et al., 2015). This emission, known as the central hard X-ray emission (CHXE), is predominantly attributed to an unresolved population of mCVs with a mean WD mass of M_WD ∼ 0.9M_⊙ (Hailey et al., 2016). However, the core of the CHXE/NSC within the central few pc region could not be resolved by NuSTAR due to the limited angular resolution and presence of hard X-ray filaments.

The central pc region corresponds to the influence radius of Sgr A*, where the immense gravity of the SMBH affects the motion of stars and compact objects. Since NSCs in other galactic nuclei cannot be spatially resolved in the X-ray band due to Mpc distances, the NSC in the GC offers a unique opportunity to study its X-ray source compositions with HEX-P. The gravitational interactions between Sgr A* and surrounding stars significantly boost the rates of binary formation; evaporation is also significantly enhanced within r ≲ 3 pc. As a result, this region should be largely inhabited by binaries with more massive compact objects, such as stellar-mass BHs, and with tighter orbits, as suggested by recent Chandra studies of the non-thermal X-ray sources detected in r ≲ 1 pc (Hailey et al., 2018; Mori et al., 2021). If the CHXE/NSC is indeed composed of hundreds of mCVs, an extremely intriguing question is whether its core shows a lower or higher concentration of X-ray point sources (predominantly mCVs), compared to the cuspy stellar distribution revealed by IR observations (Schödel et al., 2018). MCVs with heavier WD masses and shorter binary separations are more likely to survive in the NSC core. While heavier WD masses make X-ray spectra harder due to higher shock temperatures in their accretion columns, smaller binary orbits should correspond to reduced mass accretion rates (Hillman et al., 2020), and thus lower X-ray luminosities. A hard X-ray view of the NSC and its core within r < 3 pc with higher angular resolution than NuSTAR may reveal energy-dependent spatial profiles above 10 keV, where diffuse soft X-ray emission from the SNR Sgr A East becomes negligible. Depending on the (unknown) WD mass and binary separation distributions within r < 3 pc, HEX-P may unveil a sharply lower or higher hard X-ray flux compared to the CHXE, which extends out to r ∼ 8 pc. In the NSC region, a unique population of mCVs, whose thermal X-ray emission is typically characterized by plasma temperatures kT ∼ 10–40 keV, can be probed most sensitively in the 20–40 keV band. Figure 7 shows simulated radial profiles of 20–40 keV hard X-ray emission resolved by HEX-P for two possible scenarios representing a lower (“deficit”) and higher (“cusp”) concentration of mCVs within r < 3 pc. The two radial profiles are clearly distinct from each other, thanks to the higher angular resolution of HEX-P, which can resolve the CHXE and its core from other non-thermal X-ray sources.

FIGURE 7

5.5 Resolving the composition of the GRXE through the Limiting Window

While a majority of the thousands of sources detected by Chandra in the GC are believed to be mCVs (Muno et al., 2009), their specific source types, average WD masses, and their dependence on the X-ray luminosity (L_X) and radial distance from Sgr A* are still largely unknown. It is the study of diffuse X-ray emission that has revealed the presence of distinct populations of CVs between the GC and ridge. In contrast to the CHXE, which has kT ∼ 35 keV, the Galactic ridge exhibits softer diffuse X-ray emission with kT ∼ 8–15 keV (Türler et al., 2010; Yuasa et al., 2012; Perez et al., 2019). Several studies of this diffuse X-ray emission suggest that the softer ridge X-ray emission is due to an unresolved population of polars and/or nmCVs (Hailey et al., 2016; Xu et al., 2016). However, it is challenging to separate the diffuse X-ray emission into distinct CV populations (with varying L_X and kT distributions) solely through analysis of diffuse hard X-ray continuum or Fe line emission.

The Limiting Window (LW), located 1.4° south of Sgr A*, is a region with low extinction that allows for the study of X-ray sources in the GC and Bulge regions without significant obscuration. A deep exposure of the LW with Chandra has resolved a large portion of the Fe line emission (≳80%) into point sources (Revnivtsev et al., 2009). This discovery has been considered a crucial indication of the point source nature of the GRXE, as opposed to being truly diffuse in origin. However, it is still unclear what types of point sources constitute the GRXE. Revnivtsev et al. (2009) suggest that relatively faint coronal sources, such as active binaries (ABs), may constitute a major portion of the GRXE. On the other hand, it was suggested that relatively bright accretion sources, particularly mCVs, dominate the GRXE (Hong et al., 2012; Schmitt et al., 2022). These conflicting findings highlight the need for further investigation to determine the contribution of different source types in the LW and shed light on the understanding of the X-ray emission in the Bulge and Ridge regions. Indeed, the identification of the bright X-ray sources in the LW is crucial, and their individual spectral and timing properties can provide valuable insights. For instance, it is difficult to identify soft coronal sources like ABs in other typical GC and Bulge regions due to heavy obscuration, while the bright point source population in other low extinction windows in the Bulge such as Baade’s window is dominated by AGNs (Figure 3). HEX-P′s hard X-ray coverage above 10 keV can decisively determine the coronal vs. accretion nature of X-ray emission from these sources.

For our simulations, we considered different X-ray source compositions in the LW. For example, we assigned kT = 8, 15 or 35 keV to the X-ray sources in the LW. Figure 8 shows simulated LET (0.2–25 keV) and HET (3–40 keV) images of a section of the LW. Known mCVs from Hong et al. (2012) were given their best-fit model parameters specific to each source, while all other Chandra sources within this region were assumed to have spectra characterized by an absorbed APEC model with kT = 15 keV. HEX-P is expected to detect 53 and 29 sources above 4σ significance among 235 and 261 Chandra sources used as an input to SIXTE simulation for LET and HET, respectively. For different plasma temperatures (kT = 8 and 35 keV), we found that HEX-P can detect a similar number of Chandra sources in broad energy bands beyond 8 keV. Eventually, broad-band X-ray spectral data obtained by HEX-P will allow us to determine the X-ray source composition in the LW. More specifically, we aim to measure 1) plasma temperatures of the known mCVs, 2) the mean plasma temperature of other bright X-ray sources (L_X ≳ 10³² erg s⁻¹), and 3) the mean plasma temperature of faint X-ray sources (L_X ≲ 10³² erg s⁻¹). For group 1), we will analyze HEX-P and archived Chandra spectra of the known mCVs individually. For groups 2) and 3), stacking their LET and HET spectra, along with archived Chandra data, will allow us to measure their mean plasma temperatures. In all cases, the hard X-ray coverage by HEX-P is crucial to determine their plasma temperatures accurately (with % error), which can be compared with those measured in the Bulge (8 keV), Ridge (15 keV) and CHXE (35 keV). It is possible that the three groups of X-ray sources may reveal different plasma temperatures. For instance, group 1) and 2) could have higher temperatures kT ∼ 35 (predominantly IPs) and 15 keV (predominantly polars), respectively, while group 3) may have lower temperatures of kT ≲ 8 keV (predominantly nmCVs and ABs) for 3). This is how HEX-P is capable of dissecting the composition of X-ray source populations in different luminosity ranges. Moreover, HEX-P may be able to detect orbital periods from more X-ray sources in the LW as an alternate way of classifying them as polars or IPs (Hong et al., 2012; Bao and Li, 2020).

FIGURE 8

6 Diffuse X-ray sources in the Galactic center

The GC hosts a zoo of various, unique types of diffuse X-ray sources, containing both powerful particle accelerators (e.g., Sgr A* BH, SNRs, PWNe, star clusters) and targets bombarded by relativistic cosmic rays and X-rays (e.g., molecular clouds and filaments). HEX-P will survey all primary particle accelerators, molecular clouds (MCs), and X-ray filaments to obtain a whole picture of the CR spatial and energy distributions in the GC. Broad-band X-ray spectroscopy and excellent spatial resolution of HEX-P are crucial for distinguishing thermal and non-thermal diffuse X-ray emission. Combined with the CTAO GC survey, HEX-P will obtain the most accurate multi-wavelength SED data over the X-ray and TeV bands.

6.1 Molecular clouds

While the SMBH at Sgr A* is currently dormant, only occasionally emitting X-ray flares (see §4), a series of X-ray observations of the CMZ that Sgr A* used to be more active in the recent past. Following the first detection of hard X-ray emission from Sgr B2, it was suggested that the molecular cloud (MC) reflected bright X-ray flares, presumably emitted from Sgr A* (Sunyaev et al., 1993). The X-ray reflection exhibits Compton scattering, fluorescent Fe emission lines at 6.4 and 7.06 keV, and photo-absorption. In addition, both the X-ray continuum and Fe fluorescent lines have been observed to decrease over the last two decades as the X-ray flare propagated through the Sgr B2 cloud (Koyama et al., 1996; Sunyaev and Churazov, 1998; Terrier et al., 2010; 2018; Nobukawa et al., 2011; Zhang et al., 2015). Propagation of X-ray echoes is also detected from other MCs in the CMZ (e.g. Muno et al., 2007; Ponti et al., 2010; Clavel et al., 2013; Ryu et al., 2013; Chuard et al., 2018; Terrier et al., 2018). High energy observations with INTEGRAL (Sgr B2) and NuSTAR (Sgr A and B2) have shown that the reflection signal extends above 10 keV and has been detected up to 100 keV (Revnivtsev et al., 2004; Terrier et al., 2010; Mori et al., 2015; Zhang et al., 2015; Kuznetsova et al., 2022). NuSTAR played a crucial role in resolving hard X-ray continuum emission from Sgr B2 and Sgr A complex clouds, while Sgr C has not been observed with NuSTAR due to the high background contamination (Mori et al., 2015; Zhang et al., 2015; Rogers et al., 2022).

While the variable X-ray emission is primarily caused by single Compton scatterings, another non-variable X-ray component may emerge after the X-ray light front passes the clouds. This “constant” X-ray emission component may be caused by 1) multiply-scattered X-ray photons confined within dense cloud cores (Sunyaev and Churazov, 1998; Molaro et al., 2016; Khabibullin et al., 2020; Sazonov et al., 2020) and/or 2) non-thermal X-ray emission and collisional ionization by cosmic-rays (Tatischeff et al., 2012). Recent X-ray observations indicated that clouds in Sgr B2 (Kuznetsova et al., 2022; Rogers et al., 2022) and the Arches cluster (Kuznetsova et al., 2019) may have entered the constant X-ray emission phase. Diffuse TeV gamma-ray emission was observed to coincide with the MCs in the GC region, indicating that the clouds are bombarded by energetic CRs (Aharonian et al., 2006; van Eldik and H. E. S. S. Collaboration, 2008; Beilicke and VERITAS Collaboration, 2012; Sinha et al., 2022). Overall, extensive X-ray observations over the last two decades revealed that the MCs exhibit different time evolution and X-ray properties, likely due to multiple X-ray outbursts from Sgr A* and different cosmic-ray bombardment mechanisms.

HEX-P will survey all prominent MCs, including Sgr A, B1, B2 and C, in the primary GC observation program (Table 1). The HEX-P survey will capture the X-ray landscape of the MCs through both the Fe fluorescent line at 6.4 keV and non-thermal X-ray continuum emission in the early 2030s. HEX-P will likely conduct follow-up observations of these clouds later to monitor X-ray variabilities. All MC complexes with X-ray variability will be covered with at least 100 ks exposure. Detection of changes in the morphology and spectral parameters of the MC emission will determine whether X-ray echoes are still propagating through the clouds or the cosmic-ray bombardment component is dominating. The significantly lower stray-light background will allow us to investigate the hard X-ray morphology of the CMZ with HEX-P far more precisely than NuSTAR. To demonstrate the capability of HEX-P on the MC science, we present simulation results for Sgr B2, whose emission has been decreasing since the 1990s and has now likely reached a faint constant level. We adopted the latest NuSTAR (Zhang et al., 2015; Rogers et al., 2022) and INTEGRAL (Kuznetsova et al., 2022) observation results for XSPEC simulations (Zhang et al., 2015). Figure 9 shows the simulated HEX-P spectra of the Sgr B2 core. Our simulation demonstrates that the non-thermal X-ray component will be detected at significance, and the photon index and the equivalent width of the 6.4 keV emission line will be measured with error.

FIGURE 9

Furthermore, we investigated whether HEX-P can discriminate between different emission models based on broad-band X-ray spectral data. We considered two main scenarios: 1) reflection of past X-ray flares and 2) non-thermal bremsstrahlung emission and collisional excitation by low-energy cosmic-ray protons (LECRp). In the reflection scenario, the spectrum can be dominated by single or multiple scatterings, depending on the location of the X-ray light front. To clarify whether these scenarios can be discerned by HEX-P, we simulated Sgr B2 spectra for three models using 2013 NuSTAR flux measurements; the same low-energy component (absorption and astrophysical background) and either LECRp parameters (model described in Tatischeff et al., 2012) from NuSTAR (Zhang et al., 2015) or reflection parameters (CREFL16 model from Churazov et al., 2017) from INTEGRAL observations to account for the emission of Sgr B2 at high energies (Revnivtsev et al., 2004; Kuznetsova et al., 2022). Figure 10 illustrates the three different models and corresponding spectra simulated for the 90″-region of Sgr B2.

FIGURE 10

The Compton shoulder is clearly present below the Fe fluorescent lines (both for Kα and Kβ) for the reflection models (Khabibullin et al., 2020), while the LECR spectra exhibit only the narrow emission lines (Tatischeff et al., 2012). The multiple-scattering case in the middle panel can be identified well by the large line ratio between the Compton shoulder and 6.4 keV line, which should be as a result of multiple Compton scatterings (Khabibullin et al., 2020).

Overall, HEX-P′s Fe line spectroscopy and broad-band X-ray data can distinguish between the X-ray reflection and cosmic-ray models unambiguously. HEX-P′s spatial and spectral resolutions will allow us to trace the X-ray morphology evolution and characterize the Fe line profiles in detail, respectively. HEX-P′s broadband X-ray spectral data will lead to measuring the spectral index of past X-ray flares from Sgr A* and bombarding CRs. In the former case, this is a unique opportunity to investigate the mechanism of Sgr A* X-ray flares that occurred a few hundred years ago. For the latter, we can map out the population of CRs in the GC region in conjunction with future TeV observations by CTAO.

6.2 X-ray filaments

A unique phenomenon in the GC region is the existence of numerous radio and X-ray filaments (Yusef-Zadeh et al., 1984; Morris and Serabyn, 1996). The origin and formation of these filaments have been a long-standing question for decades. Within the central region, MeerKAT radio surveys have revealed over 150 filaments, which seem to be associated with bi-polar radio bubbles (Heywood et al., 2019; Yusef-Zadeh et al., 2023). The emission mechanism of non-thermal radio filaments has been pinned down to synchrotron emission thanks to the detection of linear polarization. Within the filaments, the magnetic field strengths are found to be 1–2 orders of magnitude higher than the surroundings at B ∼ 1 mG and aligned with the major axis of a filament (e.g. (Morris and Serabyn, 1996)). High-resolution JVLA observations further revealed complex sub-filaments entangled with each other within many of the radio filaments (Morris et al., 2014). It suggests that radio filaments are magnetic structures, where strong and highly-organized magnetic fields trap GeV electrons and produces synchrotron emission (Zhang et al., 2014).

At higher energies, filaments at smaller spatial scales have also been detected in the X-ray band. About 30 parsec-long X-ray filaments have been detected so far, some of which have radio counterparts (Lu et al., 2008; Muno et al., 2008; Johnson et al., 2009; Zhang et al., 2014; Ponti et al., 2015a; Zhang et al., 2020). Out of known X-ray filaments, NuSTAR was able to detect only the four brightest ones (Zhang et al., 2014; Mori et al., 2015; Zhang et al., 2020). The brightest X-ray filament Sgr A-E exhibits a featureless power-law spectrum extending up to 50 keV, which requires electron energies of multi-hundred TeV electrons that are most likely produced by hadronic interactions between primary PeV protons and nearby MCs (Zhang et al., 2020). Therefore, mapping out even a subset of the X-ray filaments could serve as a powerful way of indirectly probing PeVatron particle distribution in the GC.

HEX-P data above 40 keV will provide an unprecedentedly pristine view of non-thermal X-ray emission in the central 10 pc region. While the central PWN G359.95-0.04 and several X-ray filaments have been detected by NuSTAR above 40 keV (Mori et al., 2015; Nynka et al., 2015), other (potentially hidden) non-thermal X-ray sources remain largely unresolved. These non-thermal X-ray sources in the GC represent energetic particle accelerators (e.g., PWNe) or tracers of TeV–PeV electrons emitting synchrotron X-rays. It is still unclear just how high the maximum energy of X-rays emitted from the filaments around Sgr A* can reach (Li et al., 2013). This information is crucial for measuring the maximum energies of relativistic electrons in the region. HEX-P is poised to resolve high-energy X-ray emission from a much larger number of X-ray filaments than NuSTAR and potentially discover new filaments. For example, a 500 ks HEX-P simulation of the central 13′ × 13′field results in detecting 8 out of the 12 known X-ray filaments above 40 keV (Figure 11). A full list of all filaments included in this simulation—along with their spectral parameters—is provided in Supplementary Table S3 in the Supplementary Material. Figure 11B shows simulated HEX-P spectra for three filaments with 100 ks exposure: Sgr A-E (G359.889−0.081), G0.007−0.014 and G0.029−0.06, respectively. The simulated spectra can be fitted well with an absorbed power-law model extending beyond 50 keV. For brighter filaments like Sgr A-E and G0.029-0.06, spectral indices can be constrained to less than 5%. A spectral break or cutoff below 50 keV can be ruled out. Determination of the maximum X-ray photon energy is essential to infer the energy of parent electrons. Unlike the current and proposed soft X-ray telescopes operating below 8–10 keV, broadband X-ray spectroscopy with HEX-P can measure the X-ray spectral shapes robustly with minimal parameter degeneracy issues with N_H. Multi-wavelength SED and morphology studies of these filaments with MeerKAT/SKA and HEX-P will allow us to determine ambient B-field and broadband particle energy distribution. In summary, HEX-P observation of X-ray and radio filaments will help to reveal the origin of relativistic particles in the GC, whether they are driven by common PeVatron accelerators in the GC (e.g., Sgr A* BH) or associated with local particle accelerators such as PWNe, SNRs, stellar winds or magnetic reconnections within the filaments.

FIGURE 11

7 HEX-P science in other regions of the Milky Way

Besides the primary observation program of the GC and Bulge regions described above, below we list several unique HEX-P science cases that can be explored through GO observations. These observation ideas have not been feasible with NuSTAR or stem from expanding the previous NuSTAR observations.

7.1 Investigation of X-ray sources in star-forming regions

Outside the GC and Bulge, a HEX-P survey of several young massive clusters (YMCs) and star-formation regions in our Galaxy can be conducted to search for hidden BH–OB binaries, detect hard X-ray flares from young stellar objects (affecting the dynamics of protostellar disks) and identify hard X-ray sources previously detected by NuSTAR. YMCs, ranging in age from approximately 1 to 10 million years, are characterized by the presence of massive star-forming regions. These clusters offer valuable insights into ongoing star formation and often host compact objects such as NSs and stellar-mass BHs. This is due to the relatively short lifetimes of massive stars, which eventually undergo core-collapse supernovae. The Orion Nebula is particularly well-suited for studying a wide range of X-ray activities, including hard X-ray emission from BH binaries and flares from massive young stars. The Norma spiral arm, where NuSTAR detected hard X-ray sources, could be surveyed with HEX-P searching for HMXBs (Fornasini et al., 2017).

7.1.1 Search for black hole binaries with OB stars

It is estimated that the Milky Way harbors a substantial number of stellar mass BHs, ranging from 10⁸ to 10⁹. However, only around 20 BHs have been dynamically confirmed, and approximately 50 candidates have been identified mainly based on their X-ray characteristics. This discrepancy highlights the scarcity of confirmed BHs compared to the estimated population. In addition, these confirmed BH systems are predominantly LMXBs, in contrast to the dearth of BH-HMXBs. Currently, there are only six known HMXBs that contain a BH companion. Among them, Cygnus X-1 is the only system with an undisputed BH companion. Thus, the discovery of X-ray emission from just a few new BH–OB systems would significantly impact our understanding of these systems and their evolution.

The majority of massive stars are formed as binary or high multiple systems, as supported by studies such as Moe and Di Stefano (2017). Langer et al. (2020) estimated that % of massive OB stars have companions in the form of BHs, resulting in BH–OB systems existing within our Galaxy. The discovery of these systems would provide valuable insights into various aspects of astrophysics, including star formation processes, stellar and binary evolution, galaxy evolution, supernova rates, and gravitational wave events. Through precise proper motion measurements conducted by the Gaia mission, two BH binary systems named Gaia BH1 and BH2 have been recently discovered using their astrometric orbital solutions (El-Badry et al., 2023a; b; Chakrabarti et al., 2023). These groundbreaking findings have opened new avenues of research in the field. It is anticipated that Gaia Data Release 3 will unveil BH–OB binaries (Janssens et al., 2022), which may revolutionize our understanding of these systems. It is worth noting that Gaia would be sensitive to detecting the BH binary systems with long orbital periods (≳ 10 days). These binaries are non-interacting in nature, similar to the initial two detections (Gaia BH1 with an orbital period of 185.6 days and Gaia BH2 with 1,277 days; see Figure 12 in El-Badry et al., 2023b).

FIGURE 12

While Gaia′s sensitivity does extend to systems with shorter orbital periods, the detection of BH–OB systems with periods ≲ 10 days can be facilitated by observing X-ray emission resulting from wind accretion onto the BH (Gomez and Grindlay, 2021). This is because BH–OB systems in close binaries with short orbital periods are more likely to have accretion onto the BH from the powerful stellar wind of the OB star, leading to the generation of X-ray emission. In particular, the presence of hard X-ray emission above 10 keV serves as a definitive signature of the BH’s existence, especially when combined with knowledge of the orbital mass (which allows for differentiation between BH–OB and NS–OB systems). This is due to the fact that single OB stars or (massive) active binaries rarely emit hard X-rays. Eta Carinae stands as the sole known massive binary system exhibiting hard X-rays above 10 keV. Thus, observing such hard X-ray emission can provide crucial evidence for the presence of a BH in an BH–OB system.

In Gomez and Grindlay (2021), a novel method for identifying BH–OB systems is presented. The approach involves initially compiling a list of single-line spectroscopic binaries that exhibit no observable stellar spectrum of the secondary (with absorption lines) companion (hence, the binary is a SB1). Such systems are strong candidates to have secondary companions that are evolved stellar remnants: either a NS or a stellar-mass BH. These selected systems are then subjected to follow-up observations using NuSTAR to investigate the presence of hard X-ray emission originating from the binary. The critical signature of a NS or BH secondary is that the X-ray spectrum is detectable above 10 keV since isolated O stars have soft thermal spectra with kT ∼ 1–3 keV and L_X ∼ 10^–7L_bol ∼ 10³¹ erg s⁻¹. According to El-Badry et al. (2023a,b), BH binary systems with P_orb ≲ 10 days are expected to exhibit X-ray fluxes greater than ∼ 10^–13 erg s⁻¹ cm⁻². These estimations suggest that both soft and hard X-ray emission from close BH–OB systems in nearby massive clusters are detectable with HEX-P, for example, by surveying the central 23′ × 23′field of the Orion nebula, containing 1,600 Chandra X-ray sources (Getman et al., 2005), with 2 × 2 pointings (60 ks exposure each) and an additional 60 ksec observation at the center. HEX-P′s hard X-ray search for BH–OB binary candidates in several YMCs (alternatively, hard X-ray sources identified by HEX-P in these fields can be further investigated using follow-up optical/infrared spectroscopic studies to measure orbital mass) will complement the BH systems discovered through astrometric orbital solutions with Gaia, which are primarily sensitive to non-interacting, long-period systems.

7.1.2 Hard X-ray flares from young stellar objects

Flares from stars are the most energetic events in their star systems, but our understanding of them, including flare distributions and their long-term influence on orbiting planets, is limited, primarily to our own sun. The situation becomes even more uncertain for flares from young stellar objects (YSOs), such as pre-main sequence stars (PMS). Deep X-ray surveys, like COUP, offer valuable opportunities to gain knowledge about flares originating from YSOs. According to the standard picture, a flare occurs due to magnetic reconnection, which accelerates non-thermal electrons along magnetic field lines, resulting in the heating of the chromosphere and corona. Based on the largest flares observed in COUP, Favata et al. (2005) suggested the existence of additional flare channels in YSOs, where reconnection events in magnetic field lines connect the star and the protoplanetary disk. Getman et al. (2008a,b) also found that the brightest flare distribution, in terms of peak luminosity or temperatures, significantly deviates from expectations based on extrapolations from main sequence stars.

The evolution of a protoplanetary disk relies on the strength of the interaction between the magnetic field and ionized matter within the disk. High-temperature flares emit energetic photons that play a significant role in ionizing the disk, creating what is known as an active zone since they penetrate deep into the disk. The active zone, in turn, facilitates the transport of matter throughout the disk (Glassgold et al., 2004) through mechanisms such as the magnetorotational instability (Glassgold et al., 1997). The influence of hard X-rays is particularly important as they penetrate deeper into the disk, making the hard X-ray luminosity a key factor in determining the accretion rate of a protostar.

The simultaneous observation of the Orion Nebula using NuSTAR and Chandra aimed to study the properties of flares from YSOs in the hard X-ray range (10 keV) (Guenther 2023 in preparation). The NuSTAR data revealed a wide range of hard X-ray activity in the region, as depicted in Figure 12. However, the angular resolution of NuSTAR is relatively poor, leading to contamination from nearby sources in this densely populated area, which limits the hard X-ray investigation. Although flares and variables were detected, and some exhibited a similar pattern in their lightcurves between soft and hard X-rays, the NuSTAR lightcurves exhibited significant fluctuations due to low statistics and contamination from neighboring sources (Figure 12). In contrast, due to its significantly improved sensitivities compared to NuSTAR, HEX-P will enable a comprehensive study of five times more flares in both the soft and hard X-ray bands. Such a study will provide valuable insights into the evolution of protoplanetary disks and the formation of planets.

7.2 Magnetic cataclysmic variables

The knowledge of WD masses, metallicities, and source types (e.g., magnetic vs. nmCVs) offers valuable insights into the formation and evolution of CVs in various environments within our Galaxy (Mukai, 2017). The high sensitivity and wide spectral coverage of HEX-P present a unique opportunity to investigate the spatial evolution of the seemingly distinct CV populations in the GC, Bulge, Ridge and solar neighborhood. Also, as described in the next section, accurately measuring the WD masses of mCVs—a task that HEX-P is uniquely suited for—allows us to probe the fundamental properties of these systems, such as 1) the long-term mass evolution of mCVs toward the Chandrasekhar mass limit with IPs, and 2) the origin of WD magnetic fields with polars. In the near future, eROSITA, ZTF, and Rubin all-sky surveys are expected to discover a large number of mCVs in the Galactic Disk, including rare classes of CVs such as AM CVn stars, novae, and fast-spinning CVs (Rodriguez et al., 2023). HEX-P will provide broad-band X-ray spectral and timing data of these exotic CVs, allowing us to study and understand those new populations with physically meaningful interpretations.

7.2.1 White dwarf mass measurements from magnetic CVs

Whether mCVs gain or lose mass is of great significance for probing the two competing channels for type Ia supernovae (SNIa) (Pala and Gänsicke, 2017). In the double-degenerate case, two WDs in a binary system merge to form a WD with M > M_Ch (the Chandrasekhar mass limit, M_WD ∼ 1.4M_⊙); in the alternative single-degenerate model, a CO WD increases its mass at it accretes matter from its companion, eventually exceeding the Chandrasekhar limit. It is, however, still unknown whether CVs increase in mass over time due to accretion or lose mass through episodic nova eruptions. Interestingly, some of the newly identified IPs originally detected by INTEGRAL were found to have WD masses close to the Chandrasekhar limit (Tomsick et al., 2016; Tomsick et al., 2023). However, the INTEGRAL sources are intrinsically biased toward a population of more massive IPs, which generally possess harder X-ray spectra. In order to robustly test the CV mass evolution models such as the eCAML model (Schreiber et al., 2016), it is necessary to collate WD masses and characterize their distribution from a larger, unbiased sample of IPs. Polars, on the other hand, are more suited for studying the WD mass and B-field correlation, in order to test the common envelope (CE) model. The CE model is a leading but unconfirmed candidate for WD B-field formation where differential rotation between the decaying binary orbit and the CE leads to a magnetic dynamo, enhancing the B-field in the WD. One of the most robust predictions of CE models is an anti-correlation between B and M_WD (Briggs, 2019). Measuring the WD masses of polars will allow us to explore the B vs. M correlation because their B-fields (B = 7–240 MG) are well determined in the optical band. Since polars are typically fainter and more variable than IPs, HEX-P′s high sensitivity in the broad X-ray band will be crucial for determining their WD masses accurately.

X-ray spectral modeling of the accretion flow over a broad energy range (extending beyond keV) is a particularly effective method for measuring WD masses (Shaw et al., 2020; Vermette et al., 2023). In mCVs, infalling material is funneled onto the WD poles along B-field lines and heated to temperatures (kT_shock ∼ 10–80 keV) that scale with the WD mass. Below a standoff shock where the plasma temperature is highest, a column of infalling material cools via thermal bremsstrahlung and cyclotron radiation as it approaches the WD surface with varying temperatures and densities. Since hard X-rays are emitted from the post-shock region, hard X-ray data above E ∼ 10 keV provide the most accurate measurements of WD masses (Hailey et al., 2016). With the advent of more sophisticated X-ray spectral models for mCVs (e.g., Suleimanov et al., 2019; Hayashi et al., 2021), NuSTAR′s broadband X-ray spectral data led to measuring WD masses of IPs (Shaw et al., 2020), including new IPs discovered in the Bulge accurately (≲ 20%) (Mondal et al., 2022; Mondal et al., 2023). HEX-P will expand this program to measure WD masses from over 100 mCVs in the solar neighborhood (d ≲few × 100 pc), allowing us to parse WD mass data into different types, regions, and accretion rates. Figure 13A shows simulated HEX-P spectra of an IP yielding WD mass measurements with 10% error.

FIGURE 13

7.2.2 Exotic classes of magnetic CVs

Among the thousands of CVs discovered by extensive optical and X-ray surveys, several rare classes of CVs have been identified, such as AM CVn stars (compact CVs with short orbital periods), recurrent novae, and fast-spinning CVs. While rare, the populations and high-energy emission mechanisms of these exotic CVs are considered extremely important in contemporary astrophysics. For example, given their tight orbits, AM CVn stars are among the most likely gravitational wave sources detectable by the future LISA observatory. In another example, recurrent novae such as RS Oph, which undergo episodic thermonuclear explosions on the WD surface that produce shocks from the ejecta slamming into the surrounding circumstellar winds, are of great interest for their strong X-ray and gamma-ray emission (Cheung et al., 2022). Of particular and increasing interest is a rare class of fast-spinning CVs (FSCVs) with P_spin ≲ a few minutes. These FSCVs discovered so far are unique and distinct from regular accretion-powered CVs, including the WD “pulsar” systems (AR Sco and J191213.72-441045.1) and propeller CV (AE Aqr) (Pelisoli et al., 2023). In the near future, a larger population of undetected FSCVs may be revealed by the Rubin and eROSITA all-sky surveys (Pretorius and Mukai, 2014). HEX-P′s higher sensitivity and broad-band X-ray coverage will allow us to study a large number of FSCVs better than NuSTAR. For example, HEX-P will be able to detect the putative non-thermal X-ray component from AE Aqr above 10 keV (Figure 13B). As demonstrated for AR Sco (Takata et al., 2018), multi-wavelength observations of new WD pulsars will be invaluable for understanding their particle acceleration and non-thermal emission mechanisms.

8 Conclusion

In conclusion, HEX-P, as a probe-class mission, will offer an unprecedented opportunity to revolutionize our understanding of various astrophysical phenomena in the GC and broadly in Galactic astrophysics. With its high spatial resolution X-ray imaging, broad spectral coverage, and superior effective area, HEX-P is poised to provide groundbreaking insights into a variety of important questions in the GC. These include the past and current X-ray flares from the SMBH at Sgr A*, the populations of compact object binaries from the NSC to the GRXE, the primary particle accelerators, and cosmic-ray distributions in the GC. Outside the GC, HEX-P is expected to accurately measure WD masses of hundreds of mCVs, explore broadband X-ray spectral and timing properties of X-ray transients (including determining BH spins from Galactic BH transients), detect pulsars, and search for BH–OB binaries in young star clusters. These scientific objectives can be uniquely achieved by HEX-P or in synergy with other currently operating or future telescopes such as EHT, GRAVITY, MeerKAT, Roman, CTAO, and IceCube gen2. This paper, along with other Galactic science papers, highlights HEX-P′s enormous potential to uncover significant new insights into the most important astrophysical problems in the field of Galactic astrophysics in the 2030s.

References

• 1

  AharonianF.AkhperjanianA. G.Bazer-BachiA. R.BeilickeM.BenbowW.BergeD.et al (2006). Discovery of very-high-energy γ-rays from the galactic centre ridge. Nature439, 695–698. 10.1038/nature04467

  • CrossRef
  • Google Scholar

• 2

  AnastasopoulouK.PontiG.SormaniM. C.LocatelliN.HaberlF.MorrisM. R.et al (2023). Study of the excess Fe XXV line emission in the central degrees of the Galactic centre using XMM-Newton data. Astronomy Astrophysics671, A55. 10.1051/0004-6361/202245001

  • CrossRef
  • Google Scholar

• 3

  ArnaudK. A. (1996). “XSPEC: the first ten years,” in Astronomical data analysis software and systems V, A.S.P. Conference series, vol. 101. Editors JacobyG. H.BarnesJ., 17.

  • Google Scholar

• 4

  BaganoffF. K.BautzM. W.BrandtW. N.ChartasG.FeigelsonE. D.GarmireG. P.et al (2001). Rapid X-ray flaring from the direction of the supermassive black hole at the Galactic Centre. Nature413, 45–48. 10.1038/35092510

  • CrossRef
  • Google Scholar

• 5

  BaoT.LiZ. (2020). Periodic X-ray sources in the Galactic bulge: application of the Gregory-Loredo algorithm. MNRAS498, 3513–3534. 10.1093/mnras/staa2603

  • CrossRef
  • Google Scholar

• 6

  BarrièreN. M.TomsickJ. A.BaganoffF. K.BoggsS. E.ChristensenF. E.CraigW. W.et al (2014). NuSTAR detection of high-energy X-ray emission and rapid variability from Sagittarius A* flares. Astrophysical J.786, 46. 10.1088/0004-637X/786/1/46

  • CrossRef
  • Google Scholar

• 7

  BeilickeM.VERITAS Collaboration (2012). The Galactic Center region imaged by VERITAS. Nucl. Instrum. Methods Phys. Res. A692, 208–211. 10.1016/j.nima.2012.01.045

  • CrossRef
  • Google Scholar

• 8

  BortolasE.MapelliM.SperaM. (2017). Supernova kicks and dynamics of compact remnants in the Galactic Centre. MNRAS469, 1510–1520. 10.1093/mnras/stx930

  • CrossRef
  • Google Scholar

• 9

  BrandtT. D.KocsisB. (2015). Disrupted globular clusters can explain the galactic center gamma-ray excess. ApJ812, 15. 10.1088/0004-637X/812/1/15

  • CrossRef
  • Google Scholar

• 10

  BriggsG. P. (2019). The genesis of magnetic fields in white dwarfs. arXiv e-prints, arXiv:1907.04513.

  • Google Scholar

• 11

  BycklingK.MukaiK.ThorstensenJ. R.OsborneJ. P. (2010). Deriving an X-ray luminosity function of dwarf novae based on parallax measurements. MNRAS408, 2298–2311. 10.1111/j.1365-2966.2010.17276.x

  • CrossRef
  • Google Scholar

• 12

  ChakrabartiS.SimonJ. D.CraigP. A.ReggianiH.BrandtT. D.GuhathakurtaP.et al (2023). A noninteracting galactic black hole candidate in a binary system with a main-sequence star. Astrophysical J.166, 6. 10.3847/1538-3881/accf21

  • CrossRef
  • Google Scholar

• 13

  ChengK. S.ChernyshovD. O.DogielV. A.KoC. M.IpW. H. (2011). Origin of the Fermi bubble. Astrophysical J. Lett.731, L17. 10.1088/2041-8205/731/1/L17

  • CrossRef
  • Google Scholar

• 14

  CheungC. C.JohnsonT. J.JeanP.KerrM.PageK. L.OsborneJ. P.et al (2022). Fermi LAT gamma-ray detection of the recurrent nova RS ophiuchi during its 2021 outburst. Astrophysical J.935, 44. 10.3847/1538-4357/ac7eb7

  • CrossRef
  • Google Scholar

• 15

  ChuardD.TerrierR.GoldwurmA.ClavelM.SoldiS.MorrisM. R.et al (2018). Glimpses of the past activity of Sgr A* inferred from X-ray echoes in Sgr C. A&A610, A34. 10.1051/0004-6361/201731864

  • CrossRef
  • Google Scholar

• 16

  ChurazovE.KhabibullinI.PontiG.SunyaevR. (2017). Polarization and long-term variability of Sgr A* X-ray echo. MNRAS468, 165–179. 10.1093/mnras/stx443

  • CrossRef
  • Google Scholar

• 17

  ClavelM.TerrierR.GoldwurmA.MorrisM. R.PontiG.SoldiS.et al (2013). Echoes of multiple outbursts of Sagittarius A^⋆revealed byChandra. A&A558, A32. 10.1051/0004-6361/201321667

  • CrossRef
  • Google Scholar

• 18

  DauserT.FalknerS.LorenzM.KirschC.PeilleP.CucchettiE.et al (2019). SIXTE: a generic X-ray instrument simulation toolkit. A&A630, A66. 10.1051/0004-6361/201935978

  • CrossRef
  • Google Scholar

• 19

  DegenaarN.WijnandsR. (2010). A four-year baseline Swift study of enigmatic X-ray transients located near the Galactic center. A&A524, A69. 10.1051/0004-6361/201015322

  • CrossRef
  • Google Scholar

• 20

  DegenaarN.WijnandsR.CackettE. M.HomanJ.in’t ZandJ. J. M.KuulkersE.et al (2012). A four-year XMM-Newton/Chandra monitoring campaign of the Galactic centre: analysing the X-ray transients. A&A545, A49. 10.1051/0004-6361/201219470

  • CrossRef
  • Google Scholar

• 21

  DegenaarN.WijnandsR.MillerJ. M.ReynoldsM. T.KenneaJ.GehrelsN. (2015). The Swift X-ray monitoring campaign of the center of the Milky Way. J. High Energy Astrophysics7, 137–147. 10.1016/j.jheap.2015.03.005

  • CrossRef
  • Google Scholar

• 22

  De LucaA.CaraveoP. A.MereghettiS.TiengoA.BignamiG. F. (2006). A long-period, violently variable X-ray source in a young supernova remnant. Science313, 814–817. 10.1126/science.1129185

  • CrossRef
  • Google Scholar

• 23

  DeWittC.BandyopadhyayR. M.EikenberryS. S.BlumR.OlsenK.SellgrenK.et al (2010). Near-infrared counterparts of Chandra X-ray sources toward the galactic center. Astrophysical J.721, 1663–1679. 10.1088/0004-637X/721/2/1663

  • CrossRef
  • Google Scholar

• 24

  DexterJ.O’LearyR. M. (2014). The peculiar pulsar population of the central parsec. Astrophysical J. Lett.783, L7. 10.1088/2041-8205/783/1/L7

  • CrossRef
  • Google Scholar

• 25

  Dodds-EdenK.PorquetD.TrapG.QuataertE.HauboisX.GillessenS.et al (2009). Evidence for X-ray synchrotron emission from simultaneous mid-infrared to X-ray observations of a strong Sgr A* flare. Astrophysical J.698, 676–692. 10.1088/0004-637X/698/1/676

  • CrossRef
  • Google Scholar

• 26

  EbisawaK.TsujimotoM.PaizisA.HamaguchiK.BambaA.CutriR.et al (2005). Chandra deep X-ray observation of a typical galactic plane region and near-infrared identification. Astrophysical J.635, 214–242. 10.1086/497284

  • CrossRef
  • Google Scholar

• 27

  El-BadryK.RixH.-W.CendesY.RodriguezA. C.ConroyC.QuataertE.et al (2023a). A red giant orbiting a black hole. MNRAS521, 4323–4348. 10.1093/mnras/stad799

  • CrossRef
  • Google Scholar

• 28

  El-BadryK.RixH.-W.QuataertE.HowardA. W.IsaacsonH.FullerJ.et al (2023b). A Sun-like star orbiting a black hole. MNRAS518, 1057–1085. 10.1093/mnras/stac3140

  • CrossRef
  • Google Scholar

• 29

  EraerdsT.AntonelliV.DavisC.HallD.HetheringtonO.HollandA.et al (2021). Enhanced simulations on the athena/wide field imager instrumental background. J. Astronomical Telesc. Instrum. Syst.7, 034001. 10.1117/1.JATIS.7.3.034001

  • CrossRef
  • Google Scholar

• 30

  FavataF.FlaccomioE.RealeF.MicelaG.SciortinoS.ShangH.et al (2005). Bright X-ray flares in orion young stars from COUP: evidence for star-disk magnetic fields?Astrophysical J. Suppl. Ser.160, 469–502. 10.1086/432542

  • CrossRef
  • Google Scholar

• 31

  FornasiniF. M.TomsickJ. A.HongJ.GotthelfE. V.BauerF.RahouiF.et al (2017). The NuSTAR hard X-ray survey of the norma arm region. Astrophysical J. Suppl. Ser.229, 33. 10.3847/1538-4365/aa61fc

  • CrossRef
  • Google Scholar

• 32

  GautamA.CrockerR. M.FerrarioL.RuiterA. J.PloegH.GordonC.et al (2022). Millisecond pulsars from accretion-induced collapse as the origin of the Galactic Centre gamma-ray excess signal. Nat. Astron.6, 703–707. 10.1038/s41550-022-01658-3

  • CrossRef
  • Google Scholar

• 33

  GenerozovA.StoneN. C.MetzgerB. D.OstrikerJ. P. (2018). An overabundance of black hole X-ray binaries in the Galactic Centre from tidal captures. MNRAS478, 4030–4051. 10.1093/mnras/sty1262

  • CrossRef
  • Google Scholar

• 34

  GenzelR.EisenhauerF.GillessenS. (2010). The Galactic Center massive black hole and nuclear star cluster. Rev. Mod. Phys.82, 3121–3195. 10.1103/RevModPhys.82.3121

  • CrossRef
  • Google Scholar

• 35

  GenzelR.SchödelR.OttT.EckartA.AlexanderT.LacombeF.et al (2003). Near-infrared flares from accreting gas around the supermassive black hole at the Galactic Centre. Nature425, 934–937. 10.1038/nature02065

  • CrossRef
  • Google Scholar

• 36

  GetmanK. V.FeigelsonE. D.BroosP. S.MicelaG.GarmireG. P. (2008a). X-ray flares in orion young stars. I. Flare characteristics. Astrophysical J.688, 418–436. 10.1086/592033

  • CrossRef
  • Google Scholar

• 37

  GetmanK. V.FeigelsonE. D.MicelaG.JardineM. M.GregoryS. G.GarmireG. P. (2008b). X-ray flares in orion young stars. II. Flares, magnetospheres, and protoplanetary disks. ApJ688, 437–455. 10.1086/592034

  • CrossRef
  • Google Scholar

• 38

  GetmanK. V.FlaccomioE.BroosP. S.GrossoN.TsujimotoM.TownsleyL.et al (2005). Chandra orion ultradeep project: observations and source lists. Astrophysical J. Suppl. Ser.160, 319–352. 10.1086/432092

  • CrossRef
  • Google Scholar

• 39

  GhezA. M.WrightS. A.MatthewsK.ThompsonD.Le MignantD.TannerA.et al (2004). Variable infrared emission from the supermassive black hole at the center of the Milky way. ApJL601, L159–L162. 10.1086/382024

  • CrossRef
  • Google Scholar

• 40

  GlassgoldA. E.NajitaJ.IgeaJ. (1997). X-ray ionization of protoplanetary disks. ApJ480, 344–350. 10.1086/303952

  • CrossRef
  • Google Scholar

• 41

  GlassgoldA. E.NajitaJ.IgeaJ. (2004). Heating protoplanetary disk atmospheres. ApJ615, 972–990. 10.1086/424509

  • CrossRef
  • Google Scholar

• 42

  GomezS.GrindlayJ. E. (2021). Optical analysis and modeling of HD96670, a new black hole X-ray binary candidate. ApJ913, 48. 10.3847/1538-4357/abf24c

  • CrossRef
  • Google Scholar

• 43

  GoodenoughL.HooperD. (2009). Possible evidence for dark matter annihilation in the inner Milky way from the Fermi gamma ray space telescope. arXiv e-prints, arXiv:0910.2998. 10.48550/arXiv.0910.2998

  • CrossRef
  • Google Scholar

• 44

  GRAVITY CollaborationAbuterR.AmorimA.BauböckM.BergerJ. P.BonnetH.et al (2018). Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA. A&A618, L10. 10.1051/0004-6361/201834294

  • CrossRef
  • Google Scholar

• 45

  GRAVITY CollaborationBauböckM.DexterJ.AbuterR.AmorimA.BergerJ. P.et al (2020). Modeling the orbital motion of Sgr A*’s near-infrared flares. A&A635, A143. 10.1051/0004-6361/201937233

  • CrossRef
  • Google Scholar

• 46

  GRAVITY CollaborationAbuterR.AmorimA.BauböckM.BaganoffF.BergerJ. P.et al (2021). Constraining particle acceleration in Sgr A^⋆ with simultaneous GRAVITY, spitzer, NuSTAR, and Chandra observations. A&A654, A22. 10.1051/0004-6361/202140981

  • CrossRef
  • Google Scholar

• 47

  GregoryP. C.LoredoT. J. (1992). A new method for the detection of a periodic signal of unknown shape and period. ApJ398, 146. 10.1086/171844

  • CrossRef
  • Google Scholar

• 48

  HaggardD.NynkaM.MonB.de la Cruz HernandezN.NowakM.HeinkeC.et al (2019). Chandra spectral and timing analysis of Sgr A*’s brightest X-ray flares. ApJ886, 96. 10.3847/1538-4357/ab4a7f

  • CrossRef
  • Google Scholar

• 49

  HaileyC. J.MoriK.BauerF. E.BerkowitzM. E.HongJ.HordB. J. (2018). A density cusp of quiescent X-ray binaries in the central parsec of the Galaxy. Nature556, 70–73. 10.1038/nature25029

  • CrossRef
  • Google Scholar

• 50

  HaileyC. J.MoriK.PerezK.CanipeA. M.HongJ.TomsickJ. A.et al (2016). Evidence for intermediate polars as the origin of the galactic center hard X-ray emission. ApJ826, 160. 10.3847/0004-637X/826/2/160

  • CrossRef
  • Google Scholar

• 51

  HayashiT.KitaguchiT.IshidaM. (2021). Application of a new X-ray reflection model to V1223 Sagittarii. MNRAS504, 3651–3661. 10.1093/mnras/stab809

  • CrossRef
  • Google Scholar

• 52

  HESS CollaborationAbramowskiA.AharonianF.BenkhaliF. A.AkhperjanianA. G.AngünerE. O.et al (2016). Acceleration of petaelectronvolt protons in the galactic centre. Nature531, 476–479. 10.1038/nature17147

  • CrossRef
  • Google Scholar

• 53

  HeywoodI.CamiloF.CottonW. D.Yusef-ZadehF.AbbottT. D.AdamR. M.et al (2019). Inflation of 430-parsec bipolar radio bubbles in the Galactic Centre by an energetic event. Nature573, 235–237. 10.1038/s41586-019-1532-5

  • CrossRef
  • Google Scholar

• 54

  HeywoodI.RammalaI.CamiloF.CottonW. D.Yusef-ZadehF.AbbottT. D.et al (2022). The 1.28 GHz MeerKAT galactic center mosaic. ApJ925, 165. 10.3847/1538-4357/ac449a

  • CrossRef
  • Google Scholar

• 55

  HillmanY.SharaM. M.PrialnikD.KovetzA. (2020). A unified theory of cataclysmic variable evolution from feedback-dominated numerical simulations. Nat. Astron.4, 886–892. 10.1038/s41550-020-1062-y

  • CrossRef
  • Google Scholar

• 56

  HongJ.MoriK.HaileyC. J.NynkaM.ZhangS.GotthelfE.et al (2016). NuSTAR hard X-ray survey of the galactic center region. II. X-ray point sources. ApJ825, 132. 10.3847/0004-637X/825/2/132

  • CrossRef
  • Google Scholar

• 57

  HongJ.van den BergM.GrindlayJ. E.ServillatM.ZhaoP. (2012). Discovery of a significant magnetic cataclysmic variable population in the limiting window. ApJ746, 165. 10.1088/0004-637X/746/2/165

  • CrossRef
  • Google Scholar

• 58

  HongJ. S.van den BergM.GrindlayJ. E.LaycockS. (2009). Radial distribution of X-ray point sources near the galactic center. ApJ706, 223–237. 10.1088/0004-637X/706/1/223

  • CrossRef
  • Google Scholar

• 59

  HooperD.LindenT. (2011). On the origin of the gamma rays from the galactic center. Phys. Rev. D.84, 1230005. 10.1103/PhysRevD.84.123005

  • CrossRef
  • Google Scholar

• 60

  JanssensS.ShenarT.SanaH.FaiglerS.LangerN.MarchantP.et al (2022). Uncovering astrometric black hole binaries with massive main-sequence companions with Gaia. A&A658, A129. 10.1051/0004-6361/202141866

  • CrossRef
  • Google Scholar

• 61

  JohnsonS. P.DongH.WangQ. D. (2009). A large-scale survey of X-ray filaments in the Galactic Centre. Mon. Notices R. Astronomical Soc.399, 1429–1440. 10.1111/j.1365-2966.2009.15362.x

  • CrossRef
  • Google Scholar

• 62

  KeaneE.BhattacharyyaB.KramerM.StappersB.KeaneE. F.BhattacharyyaB.et al (2015). “A cosmic census of radio pulsars with the SKA,” in Advancing Astrophysics with the Square Kilometre Array (AASKA14), Giardini Naxos, Italy, 8 – 13 June, 2014, 40. 10.22323/1.215.0040

  • CrossRef
  • Google Scholar

• 63

  KhabibullinI.ChurazovE.SunyaevR.FederrathC.SeifriedD.WalchS. (2020). X-raying molecular clouds with a short flare: probing statistics of gas density and velocity fields. MNRAS495, 1414–1432. 10.1093/mnras/staa1262

  • CrossRef
  • Google Scholar

• 64

  KitaguchiT.AnH.BeloborodovA. M.GotthelfE. V.HayashiT.KaspiV. M.et al (2014). NuSTAR and Swift observations of the fast rotating magnetized white dwarf AE aquarii. ApJ782, 3. 10.1088/0004-637X/782/1/3

  • CrossRef
  • Google Scholar

• 65

  KoyamaK.MaedaY.SonobeT.TakeshimaT.TanakaY.YamauchiS. (1996). ASCA view of our galactic center: remains of past activities in X-rays?PASJ48, 249–255. 10.1093/pasj/48.2.249

  • CrossRef
  • Google Scholar

• 66

  KoyamaK.TakikawaY.HyodoY.InuiT.NobukawaM.MatsumotoH.et al (2009). Spatial distribution of the galactic center diffuse X-rays and the spectra of the brightest 6.4keV clumps. PASJ61, S255–S262. 10.1093/pasj/61.sp1.S255

  • CrossRef
  • Google Scholar

• 67

  KuznetsovaE.KrivonosR.ClavelM.LutovinovA.ChernyshovD.HongJ.et al (2019). Investigating the origin of the faint non-thermal emission of the Arches cluster using the 2015-2016 NuSTAR and XMM-Newton X-ray observations. MNRAS484, 1627–1636. 10.1093/mnras/stz119

  • CrossRef
  • Google Scholar

• 68

  KuznetsovaE.KrivonosR.LutovinovA.ClavelM. (2022). Sgr B2 hard X-ray emission with INTEGRAL after 2009: still detectable?MNRAS509, 1605–1613. 10.1093/mnras/stab3004

  • CrossRef
  • Google Scholar

• 69

  LangerN.SchürmannC.StollK.MarchantP.LennonD. J.MahyL.et al (2020). Properties of OB star-black hole systems derived from detailed binary evolution models. A&A638, A39. 10.1051/0004-6361/201937375

  • CrossRef
  • Google Scholar

• 70

  LiZ.MorrisM. R.BaganoffF. K. (2013). Evidence for a parsec-scale jet from the galactic center black hole: interaction with local gas. ApJ779, 154. 10.1088/0004-637X/779/2/154

  • CrossRef
  • Google Scholar

• 71

  LuF. J.YuanT. T.LouY.-Q. (2008). An imaging and spectral study of 10 X-ray filaments around the galactic center. ApJ673, 915–927. 10.1086/524369

  • CrossRef
  • Google Scholar

• 72

  MadsenK. K.GarcíaJ. A.SternD.ArminiR.BassoS.CoutinhoD.et al (2023). The high energy X-ray probe (HEX-P): instrument and mission profile. arXiv e-printsarXiv:2312.04678

  • Google Scholar

• 73

  MeidingerN.AlbrechtS.BeitlerC.BonholzerM.EmbergerV.FrankJ.et al (2020). “Development status of the wide field imager instrument for Athena,” in Space telescopes and instrumentation 2020: ultraviolet to gamma ray (SPIE). 10.1117/12.2560507

  • CrossRef
  • Google Scholar

• 74

  MoeM.Di StefanoR. (2017). Mind your ps and qs: the interrelation between period (P) and mass-ratio (Q) distributions of binary stars. ApJS230, 15. 10.3847/1538-4365/aa6fb6

  • CrossRef
  • Google Scholar

• 75

  MolaroM.KhatriR.SunyaevR. A. (2016). Probing the clumping structure of giant molecular clouds through the spectrum, polarisation and morphology of X-ray reflection nebulae. A&A589, A88. 10.1051/0004-6361/201527760

  • CrossRef
  • Google Scholar

• 76

  MondalS.PontiG.HaberlF.AnastasopoulouK.CampanaS.MoriK.et al (2022). An intermediate polar candidate toward the Galactic plane. A&A666, A150. 10.1051/0004-6361/202244264

  • CrossRef
  • Google Scholar

• 77

  MondalS.PontiG.HaberlF.MoriK.ReaN.MorrisM. R.et al (2023). Discovery of periodicities in two highly variable intermediate polars towards the Galactic centre. A&A671, A120. 10.1051/0004-6361/202245553

  • CrossRef
  • Google Scholar

• 78

  MoriK.GotthelfE. V.ZhangS.AnH.BaganoffF. K.BarrièreN. M.et al (2013). NuSTAR discovery of a 3.76 s transient magnetar near Sagittarius A. ApJL770, L23. 10.1088/2041-8205/770/2/L23

  • CrossRef
  • Google Scholar

• 79

  MoriK.HaileyC. J.KrivonosR.HongJ.PontiG.BauerF.et al (2015). NuSTAR hard X-ray survey of the galactic center region I: hard X-ray morphology and spectroscopy of the diffuse emission. ApJ814, 94. 10.1088/0004-637X/814/2/94

  • CrossRef
  • Google Scholar

• 80

  MoriK.HaileyC. J.MandelS.SchuttY. E.BachettiM.CoerverA.et al (2019). NuSTAR and Chandra observations of new X-ray transients in the central parsec of the Galaxy. ApJ885, 142. 10.3847/1538-4357/ab4b47

  • CrossRef
  • Google Scholar

• 81

  MoriK.HaileyC. J.SchuttT. Y. E.MandelS.HeuerK.GrindlayJ. E.et al (2021). The X-ray binary population in the galactic center revealed through multi-decade observations. ApJ921, 148. 10.3847/1538-4357/ac1da5

  • CrossRef
  • Google Scholar

• 82

  MorihanaK.TsujimotoM.EbisawaK.GandhiP. (2022). Deep near-infrared imaging observation of the faint X-ray point sources constituting the Galactic bulge X-ray emission. PASJ74, 283–297. 10.1093/pasj/psab124

  • CrossRef
  • Google Scholar

• 83

  MorrisM.SerabynE. (1996). The galactic center environment. Annu. Rev. Astronomy Astrophysics34, 645–701. 10.1146/annurev.astro.34.1.645

  • CrossRef
  • Google Scholar

• 84

  MorrisM. R.ZhaoJ. H.GossW. M. (2014). “Nonthermal filamentary radio features within 20 pc of the Galactic center,” in The galactic center: feeding and feedback in a normal galactic nucleus. Editors SjouwermanL. O.LangC. C.OttJ., 369–373.

  • Google Scholar

• 85

  MossouxE.FinocietyB.BeckersJ. M.VincentF. H. (2020). Continuation of the X-ray monitoring of Sgr A*: the increase in bright flaring rate confirmed. A&A636, A25. 10.1051/0004-6361/201937136

  • CrossRef
  • Google Scholar

• 86

  MukaiK. (2017). X-ray emissions from accreting white dwarfs: a review. Publ. Astronomical Soc. Pac.129, 062001. 10.1088/1538-3873/aa6736

  • CrossRef
  • Google Scholar

• 87

  MunoM. P.BaganoffF. K.BautzM. W.BrandtW. N.BroosP. S.FeigelsonE. D.et al (2003). A deep Chandra catalog of X-ray point sources toward the galactic center. ApJ589, 225–241. 10.1086/374639

  • CrossRef
  • Google Scholar

• 88

  MunoM. P.BaganoffF. K.BrandtW. N.MorrisM. R.StarckJ.-L. (2008). A catalog of diffuse X-ray-emitting features within 20 pc of Sagittarius A*: twenty pulsar wind nebulae?ApJ673, 251–263. 10.1086/521641

  • CrossRef
  • Google Scholar

• 89

  MunoM. P.BaganoffF. K.BrandtW. N.ParkS.MorrisM. R. (2007). Discovery of variable iron fluorescence from reflection nebulae in the galactic center. ApJL656, L69–L72. 10.1086/512236

  • CrossRef
  • Google Scholar

• 90

  MunoM. P.BauerF. E.BaganoffF. K.BandyopadhyayR. M.BowerG. C.BrandtW. N.et al (2009). A catalog of X-ray point sources from two megaseconds of Chandra observations of the galactic center. ApJS181, 110–128. 10.1088/0067-0049/181/1/110

  • CrossRef
  • Google Scholar

• 91

  MunoM. P.PfahlE.BaganoffF. K.BrandtW. N.GhezA.LuJ.et al (2005). An overabundance of transient X-ray binaries within 1 parsec of the galactic center. ApJL622, L113–L116. 10.1086/429721

  • CrossRef
  • Google Scholar

• 92

  NandraK.BarretD.BarconsX.FabianA.den HerderJ.-W.PiroL.et al (2013). The hot and energetic universe: a white paper presenting the science theme motivating the Athena+ mission. arXiv e-prints, arXiv:1306.2307. 10.48550/arXiv.1306.2307

  • CrossRef
  • Google Scholar

• 93

  NeilsenJ.NowakM. A.GammieC.DexterJ.MarkoffS.HaggardD.et al (2013). A Chandra/HETGS census of X-ray variability from Sgr A* during 2012. ApJ774, 42. 10.1088/0004-637X/774/1/42

  • CrossRef
  • Google Scholar

• 94

  NobukawaM.RyuS. G.TsuruT. G.KoyamaK. (2011). New evidence for high activity of the supermassive black hole in our galaxy. ApJL739, L52. 10.1088/2041-8205/739/2/L52

  • CrossRef
  • Google Scholar

• 95

  NowakM. A.NeilsenJ.MarkoffS. B.BaganoffF. K.PorquetD.GrossoN.et al (2012). Chandra/HETGS observations of the brightest flare seen from Sgr A. ApJ759, 95. 10.1088/0004-637X/759/2/95

  • CrossRef
  • Google Scholar

• 96

  NynkaM.HaileyC. J.MoriK.BaganoffF. K.BauerF. E.BoggsS. E.et al (2013). High-energy X-rays from j174545.5-285829, the Cannonball: a candidate pulsar wind nebula associated with Sgr A East. ApJL778, L31. 10.1088/2041-8205/778/2/L31

  • CrossRef
  • Google Scholar

• 97

  NynkaM.HaileyC. J.ReynoldsS. P.AnH.BaganoffF. K.BoggsS. E.et al (2015). G359.97-0.038: a hard X-ray filament associated with a supernova shell-molecular cloud interaction. ApJ800, 119. 10.1088/0004-637X/800/2/119

  • CrossRef
  • Google Scholar

• 98

  PalaA. F.GänsickeB. T. (2017). “Dynamical masses of accreting white dwarfs,” in 20th European white dwarf workshop. Editors TremblayP. E.GaensickeB.MarshT. (San Francisco: Astronomical Society of the Pacific), 91.

  • Google Scholar

• 99

  PanamarevT.JustA.SpurzemR.BerczikP.WangL.Arca SeddaM. (2019). Direct N-body simulation of the Galactic centre. MNRAS484, 3279–3290. 10.1093/mnras/stz208

  • CrossRef
  • Google Scholar

• 100

  PelisoliI.MarshT. R.BuckleyD. A. H.HeywoodI.PotterS. B.SchwopeA.et al (2023). A 5.3-min-period pulsing white dwarf in a binary detected from radio to X-rays. Nat. Astron.7, 931–942. 10.1038/s41550-023-01995-x

  • CrossRef
  • Google Scholar

• 101

  PerezK.HaileyC. J.BauerF. E.KrivonosR. A.MoriK.BaganoffF. K.et al (2015). Extended hard-X-ray emission in the inner few parsecs of the Galaxy. Nature520, 646–649. 10.1038/nature14353

  • CrossRef
  • Google Scholar

• 102

  PerezK.KrivonosR.WikD. R. (2019). The galactic bulge diffuse emission in broadband X-rays with NuSTAR. ApJ884, 153. 10.3847/1538-4357/ab4590

  • CrossRef
  • Google Scholar

• 103

  PfeifleR. W.BoormanR. G.WeaverK. A.BuchnerJ.CivanoF.MadsenK.et al (2023). The high energy X-ray Probe (HEX-P): the future of hard X-ray dual AGN science. arXiv e-prints. 10.48550/arXiv.2311.05154

  • CrossRef
  • Google Scholar

• 104

  PontiG.De MarcoB.MorrisM. R.MerloniA.Muñoz-DariasT.ClavelM.et al (2015a). Fifteen years of XMM-Newton and Chandra monitoring of Sgr A*: evidence for a recent increase in the bright flaring rate. MNRAS454, 1525–1544. 10.1093/mnras/stv1537

  • CrossRef
  • Google Scholar

• 105

  PontiG.GeorgeE.ScaringiS.ZhangS.JinC.DexterJ.et al (2017). A powerful flare from Sgr A* confirms the synchrotron nature of the X-ray emission. MNRAS468, 2447–2468. 10.1093/mnras/stx596

  • CrossRef
  • Google Scholar

• 106

  PontiG.MorrisM. R.ChurazovE.HeywoodI.FenderR. P. (2021). The Galactic center chimneys: the base of the multiphase outflow of the Milky Way. A&A646, A66. 10.1051/0004-6361/202039636

  • CrossRef
  • Google Scholar

• 107

  PontiG.MorrisM. R.TerrierR.HaberlF.SturmR.ClavelM.et al (2015b). The XMM-Newton view of the central degrees of the Milky Way. MNRAS453, 172–213. 10.1093/mnras/stv1331

  • CrossRef
  • Google Scholar

• 108

  PontiG.TerrierR.GoldwurmA.BelangerG.TrapG. (2010). Discovery of a superluminal Fe K echo at the galactic center: the glorious past of Sgr A* preserved by molecular clouds. ApJ714, 732–747. 10.1088/0004-637X/714/1/732

  • CrossRef
  • Google Scholar

• 109

  PretoriusM. L.MukaiK. (2014). Constraints on the space density of intermediate polars from the Swift-BAT survey. MNRAS442, 2580–2585. 10.1093/mnras/stu990

  • CrossRef
  • Google Scholar

• 110

  RevnivtsevM.SazonovS.ChurazovE.FormanW.VikhlininA.SunyaevR. (2009). Discrete sources as the origin of the galactic x-ray ridge emission. Nature458, 1142–1144. 10.1038/nature07946

  • CrossRef
  • Google Scholar

• 111

  RevnivtsevM. G.ChurazovE. M.SazonovS. Y.SunyaevR. A.LutovinovA. A.GilfanovM. R.et al (2004). Hard X-ray view of the past activity of Sgr A* in a natural Compton mirror. A&A425, L49–L52. 10.1051/0004-6361:200400064

  • CrossRef
  • Google Scholar

• 112

  RodriguezA. C.GaliullinI.GilfanovM.KulkarniS. R.KhamitovI.BikmaevI.et al (2023). SRGeJ045359.9+622444: a 55-min period eclipsing AM CVn discovered from a joint SRG/eROSITA + ZTF search. arXiv e-prints, arXiv:2306.13133. 10.48550/arXiv.2306.13133

  • CrossRef
  • Google Scholar

• 113

  RogersF.ZhangS.PerezK.ClavelM.TaylorA. (2022). New constraints on cosmic particle populations at the galactic center using X-ray observations of the molecular cloud Sagittarius B2. ApJ934, 19. 10.3847/1538-4357/ac7717

  • CrossRef
  • Google Scholar

• 114

  RyuS. G.NobukawaM.NakashimaS.TsuruT. G.KoyamaK.UchiyamaH. (2013). X-ray echo from the Sagittarius C complex and 500-year activity history of Sagittarius A. PASJ65, 33. 10.1093/pasj/65.2.33

  • CrossRef
  • Google Scholar

• 115

  SazonovS.PaizisA.BazzanoA.ChelovekovI.KhabibullinI.PostnovK.et al (2020). The galactic LMXB population and the galactic centre region. New Astron. Rev.88, 101536. 10.1016/j.newar.2020.101536

  • CrossRef
  • Google Scholar

• 116

  SchmittJ. H. M. M.CzeslaS.SchneiderP. C.FreundS.RobradeJ. (2022). The nature of the X-ray sources constituting the 6.7 keV Galactic ridge emission. A&A661, A88. 10.1051/0004-6361/202142876

  • CrossRef
  • Google Scholar

• 117

  SchödelR.Gallego-CanoE.DongH.Nogueras-LaraF.Gallego-CalventeA. T.Amaro-SeoaneP.et al (2018). The distribution of stars around the Milky Way’s central black hole. II. Diffuse light from sub-giants and dwarfs. A&A609, A27. 10.1051/0004-6361/201730452

  • CrossRef
  • Google Scholar

• 118

  SchoedelR.LongmoreS.HenshawJ.GinsburgA.BallyJ.FeldmeierA.et al (2023). The JWST galactic center survey – a white paper. arXiv e-prints, arXiv:2310.11912. 10.48550/arXiv.2310.11912

  • CrossRef
  • Google Scholar

• 119

  SchreiberM. R.ZorotovicM.WijnenT. P. G. (2016). Three in one go: consequential angular momentum loss can solve major problems of CV evolution. MNRAS455, L16–L20. 10.1093/mnrasl/slv144

  • CrossRef
  • Google Scholar

• 120

  ShawA. W.HeinkeC. O.MukaiK.TomsickJ. A.DoroshenkoV.SuleimanovV. F.et al (2020). Measuring the masses of magnetic white dwarfs: a NuSTAR legacy survey. MNRAS498, 3457–3469. 10.1093/mnras/staa2592

  • CrossRef
  • Google Scholar

• 121

  SinhaA.BaghmanyanV.PeronG.GallantY.CasanovaS.HollerM.et al (2022). “Search for enhanced TeV gamma ray emission from Giant Molecular Clouds using H.E.S.S,” in 37th International Cosmic Ray Conference, Berlin, Germany, 12-23 July 2021. 10.22323/1.395.0790

  • CrossRef
  • Google Scholar

• 122

  SuleimanovV. F.DoroshenkoV.WernerK. (2019). Hard X-ray view on intermediate polars in the Gaia era. MNRAS482, 3622–3635. 10.1093/mnras/sty2952

  • CrossRef
  • Google Scholar

• 123

  SunyaevR.ChurazovE. (1998). Equivalent width, shape and proper motion of the iron fluorescent line emission from molecular clouds as an indicator of the illuminating source X-ray flux history. MNRAS297, 1279–1291. 10.1046/j.1365-8711.1998.01684.x

  • CrossRef
  • Google Scholar

• 124

  SunyaevR. A.MarkevitchM.PavlinskyM. (1993). The center of the Galaxy in the recent past - a view from GRANAT. ApJ407, 606–610. 10.1086/172542

  • CrossRef
  • Google Scholar

• 125

  SzölgyénÁ.KocsisB. (2018). Black hole disks in galactic nuclei. Phys. Rev. Lett.121, 101101. 10.1103/PhysRevLett.121.101101

  • CrossRef
  • Google Scholar

• 126

  TakataJ.HuC. P.LinL. C. C.TamP. H. T.PalP. S.HuiC. Y.et al (2018). A non-thermal pulsed X-ray emission of AR scorpii. ApJ853, 106. 10.3847/1538-4357/aaa23d

  • CrossRef
  • Google Scholar

• 127

  TatischeffV.DecourchelleA.MaurinG. (2012). Nonthermal X-rays from low-energy cosmic rays: application to the 6.4 keV line emission from the Arches cluster region. A&A546, A88. 10.1051/0004-6361/201219016

  • CrossRef
  • Google Scholar

• 128

  TerrierR.ClavelM.SoldiS.GoldwurmA.PontiG.MorrisM. R.et al (2018). An X-ray survey of the central molecular zone: variability of the Fe Kα emission line. A&A612, A102. 10.1051/0004-6361/201730837

  • CrossRef
  • Google Scholar

• 129

  TerrierR.PontiG.BélangerG.DecourchelleA.TatischeffV.GoldwurmA.et al (2010). Fading hard X-ray emission from the galactic center molecular cloud Sgr B2. ApJ719, 143–150. 10.1088/0004-637X/719/1/143

  • CrossRef
  • Google Scholar

• 130

  TomsickJ. A.KumarS. G.CoughenourB. M.ShawA. W.MukaiK.HareJ.et al (2023). Classifying IGR J15038-6021 as a magnetic CV with a massive white dwarf. arXiv e-prints, arXiv:2306.04691. 10.48550/arXiv.2306.04691

  • CrossRef
  • Google Scholar

• 131

  TomsickJ. A.RahouiF.KrivonosR.ClavelM.StraderJ.ChomiukL. (2016). Identifying IGR J14091-6108 as a magnetic CV with a massive white dwarf using X-ray and optical observations. MNRAS460, 513–523. 10.1093/mnras/stw871

  • CrossRef
  • Google Scholar

• 132

  TürlerM.ChernyakovaM.CourvoisierT. J.-L.LubińskiP.NeronovA.ProduitN.et al (2010). INTEGRAL hard X-ray spectra of the cosmic X-ray background and Galactic ridge emission. A&A512, A49. 10.1051/0004-6361/200913072

  • CrossRef
  • Google Scholar

• 133

  van EldikC.H. E. S. S. Collaboration (2008). The H.E.S.S. view of the Galactic Centre region. Nucl. Instrum. Methods Phys. Res. A588, 72–75. 10.1016/j.nima.2008.01.025

  • CrossRef
  • Google Scholar

• 134

  VermetteB.SalcedoC.MoriK.GerberJ.YoonK. D.BridgesG.et al (2023). Constraining white dwarf mass and magnetic field strength of a new intermediate polar through X-ray observations. arXiv e-prints, arXiv:2307.11910. 10.48550/arXiv.2307.11910

  • CrossRef
  • Google Scholar

• 135

  WangQ. D.GotthelfE. V.LangC. C. (2002). A faint discrete source origin for the highly ionized iron emission from the Galactic Centre region. Nature415, 148–150. 10.1038/415148a

  • CrossRef
  • Google Scholar

• 136

  WielgusM.MoscibrodzkaM.VosJ.GellesZ.Martí-VidalI.FarahJ.et al (2022). Orbital motion near Sagittarius A*. Constraints from polarimetric ALMA observations. A&A665, L6. 10.1051/0004-6361/202244493

  • CrossRef
  • Google Scholar

• 137

  XuX.-j.LiZ.ZhuZ.ChengZ.LiX.-d.YuZ.-l. (2019a). Massive white dwarfs in the galactic center: a Chandra X-ray spectroscopy of cataclysmic variables. ApJ882, 164. 10.3847/1538-4357/ab32df

  • CrossRef
  • Google Scholar

• 138

  XuX.-j.WangQ. D.LiX.-D. (2016). Fe line diagnostics of cataclysmic variables and Galactic Ridge X-Ray emission. ApJ818, 136. 10.3847/0004-637X/818/2/136

  • CrossRef
  • Google Scholar

• 139

  XuX.-j.YuZ.-l.LiX.-d. (2019b). The Fe line flux ratio as a diagnostic of the maximum temperature and the white dwarf mass of cataclysmic variables. ApJ878, 53. 10.3847/1538-4357/ab1fe1

  • CrossRef
  • Google Scholar

• 140

  YuasaT.MakishimaK.NakazawaK. (2012). Broadband spectral analysis of the Galactic Ridge X-Ray emission. ApJ753, 129. 10.1088/0004-637X/753/2/129

  • CrossRef
  • Google Scholar

• 141

  YuasaT.TamuraK.-I.NakazawaK.KokubunM.MakishimaK.BambaA.et al (2008). Suzaku detection of extended/diffuse hard X-ray emission from the galactic center. PASJ60, S207–S221. 10.1093/pasj/60.sp1.S207

  • CrossRef
  • Google Scholar

• 142

  Yusef-ZadehF.ArendtR. G.WardleM.HeywoodI. (2023). The population of the galactic center filaments: position angle distribution reveals a degree-scale collimated outflow from Sgr A* along the galactic plane. ApJL949, L31. 10.3847/2041-8213/acd54b

  • CrossRef
  • Google Scholar

• 143

  Yusef-ZadehF.MorrisM.ChanceD. (1984). Large, highly organized radio structures near the galactic centre. Nature310, 557–561. 10.1038/310557a0

  • CrossRef
  • Google Scholar

• 144

  ZhangS.BaganoffF. K.PontiG.NeilsenJ.TomsickJ. A.DexterJ.et al (2017). Sagittarius A* high-energy X-ray flare properties duringNuStarMonitoring of the galactic center from 2012 to 2015. ApJ843, 96. 10.3847/1538-4357/aa74e8

  • CrossRef
  • Google Scholar

• 145

  ZhangS.HaileyC. J.BaganoffF. K.BauerF. E.BoggsS. E.CraigW. W.et al (2014). High-energy X-ray detection of g359.89-0.08 (Sgr A-E): magnetic flux tube emission powered by cosmic rays?ApJ784, 6. 10.1088/0004-637X/784/1/6

  • CrossRef
  • Google Scholar

• 146

  ZhangS.HaileyC. J.BaganoffF. K.BauerF. E.BoggsS. E.CraigW. W.et al (2015). Hard X-ray morphological and spectral studies of the galactic center molecular cloud Sgr B2: constraining past Sgr A* flaring activity. Astrophysical J.815, 132. 10.1088/0004-637X/815/2/132

  • CrossRef
  • Google Scholar

• 147

  ZhangS.ZhuZ.LiH.PashamD.LiZ.ClavelM.et al (2020). NuSTAR and Chandra observations of the galactic center nonthermal X-ray filament g0.13-0.11: a pulsar-wind-nebula-driven magnetic filament. ApJ893, 3. 10.3847/1538-4357/ab7dc1

  • CrossRef
  • Google Scholar

• 148

  ZhaoJ.-H.MorrisM. R.GossW. M. (2013). Radio detection of a candidate neutron star associated with galactic center supernova remnant Sagittarius A East. ApJ777, 146. 10.1088/0004-637X/777/2/146

  • CrossRef
  • Google Scholar

• 149

  ZhuZ.LiZ.MorrisM. R. (2018). An ultradeep Chandra catalog of X-ray point sources in the galactic center star cluster. ApJS235, 26. 10.3847/1538-4365/aab14f

  • CrossRef
  • Google Scholar