The soft X-ray emission toward the inner halo of M31 contains contributions from both M31 and the Milky Way, whose thermal spectra are nearly degenerate at CCD resolution. We analyze 22 XMM-Newton pointings at projected radii of approximately 10–31 kpc and treat the directly measured cool component as the full line-of-sight sum rather than assigning it to either galaxy. Our homogeneous baseline excludes the MOS 1.4–2.0 keV instrumental-line interval, uses a pointing-dependent 15 arcmin HI4PI atomic hydrogen column, and retains 14 quality-selected dual-MOS fields in the primary sample. The primary spectra give \(kT=0.175\pm0.004\,\mathrm{keV}\) and an inverse-variance all-field absorbed 0.4–1.25 keV surface flux of \((1.092\pm0.052)\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). The unequal radial coverage of the 10 North/NW and four South/SE fields prevents a clean side comparison; the measured contrast is not statistically compelling. Under the named reference Combined (\(\beta=0.5\)) Milky Way model of Locatelli et al., a response-matched projection through the actual footprint gives an estimator-matched all-field \(x=1.602\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), with individual-field values of \(1.431\)–\(1.668\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). This exceeds the measured total and therefore leaves no non-negative central M31 residual under that model. We interpret this as a conditional tension among the Milky Way model, line-to-band conversion, and measured broadband sum, not as evidence for negative M31 emission. The direct result remains a spatially widespread line-of-sight component; any distance decomposition requires an explicit spatial or foreground prior.
Keywords: Andromeda Galaxy; circumgalactic medium; diffuse X-ray background; X-ray astronomy; Milky Way Galaxy
Introduction
Hot circumgalactic gas is expected to contain a substantial fraction of the baryons associated with massive spiral galaxies, but its low density makes direct X-ray measurements difficult. Deep observations of individual massive spirals and stacked samples demonstrate extended soft X-ray halos, while also showing that inferred profiles depend on galaxy mass, unresolved-source subtraction, and background treatment (Anderson et al. 2016; Li et al. 2017, 2018; Zhang, Comparat, Ponti, Merloni, Nandra, Haberl, Locatelli, et al. 2024; Zhang, Comparat, Ponti, Merloni, Nandra, Haberl, Truong, et al. 2024). M31 is close enough that its inner halo can be sampled with multiple XMM-Newton fields, but its large angular extent creates a different limitation: emission from the M31 halo is superposed on spatially structured soft emission from the Milky Way and potentially the Local Group.
The central few kiloparsecs of M31 contain multiphase diffuse X-ray emission associated with the bulge and disk (Shirey et al. 2001; Takahashi et al. 2004; Li and Wang 2007). Those measurements do not determine the amplitude of cooler gas at projected radii of 10–31 kpc. Off-disc XMM-Newton observations provide a closer spectral precedent. Kavanagh et al. (2020) measured a component near \(0.17\)–\(0.18\,\mathrm{keV}\) and explicitly modeled it as the combined emission of the Galaxy, M31, and the intergalactic medium. A possible Local Hot Bridge in the direction of M31 has a similar inferred temperature (Qu et al. 2021). Milky Way surveys also find characteristic temperatures near \(0.15\)–\(0.22\,\mathrm{keV}\) but substantial spatial variation (Henley and Shelton 2013a, 2015; Ponti et al. 2023; Kaaret et al. 2020; Ueda et al. 2022; Locatelli et al. 2024).
These overlapping temperature scales make distance attribution from an EPIC spectrum alone unreliable. We therefore organize the analysis around the directly observed line-of-sight component, \[\begin{equation} S_{\rm obs}(E)=S_{\rm MW}(E)+S_{\rm M31}(E), \label{eq:direct-sum} \end{equation}\] with an optional third Local-Group term in explicitly labeled conditional scenarios. We ask two assumption-light questions. First, what cool thermal surface flux is measured toward the North/NW and South/SE inner-halo fields? Second, how strongly do the data constrain a difference between the two sampled sides? The M31-versus-Milky-Way decomposition is discussed only after these direct measurements.
Observations and Data Reduction
XMM-Newton program and field selection
We use observations from XMM-Newton proposal 080073, “XMM-Newton Legacy Survey of M31 Halo: Searching for the Missing Accreted Hot CGM” (PI J.-T. Li), obtained with the EPIC cameras (Jansen et al. 2001; Turner et al. 2001). The selected observations were taken between 2017 June 28 and 2018 January 17. All 43 detector products used here have the Thin1 optical blocking filter. The v19 analysis is MOS-only: 21 fields have one MOS1 and one MOS2 spectrum, while 0800730101 contributes MOS2 only.
EPIC-pn products are not used in this homogeneous baseline because
the retained preliminary pn branch does not provide an equivalently
validated particle-background set. In the audited
PATTERN<=4 pn branch of the local
SAS 20.0.0 reduction, all 22 fields have a source spectrum,
but only 21 have a matched QPB spectrum in the same processing branch
(the 22nd, 0800730101, has a paired source and QPB in a neighboring
PATTERN<=0 branch but not in this one, so “matched” is
defined within the same event-selection and processing branch), and 17
of the 22 source PHA files record a cleaned exposure shorter than 10 ks.
This is the regime in which the current XMM-ESAS Cookbook identifies a
specific pn limitation, known as “The Corner Problem”: the small
out-of-field-of-view corner sample has poor counting statistics, is
affected by out-of-time events, and is more sensitive than previously
recognized to residual Soft Proton flares, so the inferred pn QPB can be
under- or over-estimated, including possible oversubtraction (XMM-Newton Science
Analysis Centre 2025). We therefore defer pn inclusion until
those spectra and their corner-normalized QPB products are reprocessed
and validated independently across the full sample. This is a
product-level background-systematics choice for the local
SAS 20.0.0 pn branch, not a claim that EPIC-pn data are
generally unusable or that all current SAS 22.0 pn products
share this limitation. The per-field audit is provided in
m31_cgmsum_v19_pn_qpb_readiness.csv.
The 22-field parent set is the proposal-080073 subset for which the
upstream project screening retained the required diffuse-field products.
Observation 0800731701 is absent from the local reduction tree,
0800732401 lacks the required 4.background products, and
0800731201, 0800731401, 0800731801, 0800732501, 0800732901, and
0800733001 were rejected upstream because of extreme 2.0–3.2 keV
observed-to-QPB ratios. Observation 0800733101 was rejected because of
extreme residual Soft Proton contamination and incomplete MOS coverage.
These exclusions define the parent set before the v19 spectral-quality
cuts described in Section 3.4; they are
not selected on fitted \(\mathrm{CGM}_{\mathrm{sum}}\) temperature
or surface flux.
Pointing centers are taken from the final PHA headers. For an adopted M31 center and distance of 785 kpc (McConnachie et al. 2005), their circular projected radii span 10.5–30.8 kpc. The pointings form two approximately minor-axis chains, labeled North/NW and South/SE, rather than a complete azimuthal annulus. In the final primary sample, the median radii are 19.28 and 14.97 kpc on the two sides, respectively, and no South/SE primary field lies at 20–25 kpc. Table 1 gives the observation-level inventory.
Table 1. XMM-Newton Observation and v19 Product Inventory
| ObsID | R.A. (deg) | Decl. (deg) | Start date (UTC) | Side | Rproj (kpc) | MOS1 (ks) | MOS2 (ks) | NHI (10²⁰ cm⁻²) | v19 use |
|---|---|---|---|---|---|---|---|---|---|
| 0800730101 | 9.2301 | 43.2408 | 2017-06-28 | North/NW | 30.78 | – | 10.01 | 5.58 | Numerical only |
| 0800730201 | 8.8793 | 42.8879 | 2017-07-25 | North/NW | 28.79 | 13.92 | 11.49 | 5.78 | Primary |
| 0800730301 | 8.4932 | 42.5826 | 2017-07-14 | North/NW | 28.69 | 16.11 | 18.08 | 5.76 | Primary |
| 0800730401 | 8.7369 | 42.2086 | 2017-07-27 | North/NW | 23.71 | 14.51 | 14.86 | 5.46 | Excluded |
| 0800730501 | 9.0777 | 42.4873 | 2017-07-28 | North/NW | 23.40 | 11.83 | 14.28 | 5.53 | Primary |
| 0800730601 | 9.3863 | 42.8204 | 2017-08-12 | North/NW | 25.03 | 11.57 | 11.10 | 5.79 | Primary |
| 0800730701 | 9.5865 | 42.4375 | 2017-08-12 | North/NW | 19.54 | 13.61 | 14.35 | 5.73 | Primary |
| 0800730801 | 9.2710 | 42.1599 | 2017-08-10 | North/NW | 18.92 | 14.49 | 14.75 | 5.60 | Primary |
| 0800730901 | 9.0152 | 41.8738 | 2018-01-07 | North/NW | 19.01 | 18.41 | 18.80 | 5.59 | Primary |
| 0800731001 | 9.1456 | 41.5646 | 2017-08-10 | North/NW | 16.32 | 14.82 | 15.09 | 5.79 | Excluded |
| 0800731101 | 9.8256 | 42.0115 | 2018-01-07 | North/NW | 13.45 | 16.48 | 18.05 | 5.67 | Primary |
| 0800731301 | 10.6467 | 42.2062 | 2018-01-11 | North/NW | 12.85 | 6.14 | 7.06 | 7.64 | Excluded |
| 0800731501 | 9.7074 | 41.4971 | 2017-08-12 | North/NW | 10.52 | 14.73 | 14.98 | 5.85 | Primary |
| 0800731601 | 9.4171 | 41.2226 | 2017-08-02 | North/NW | 13.07 | 20.57 | 20.27 | 6.30 | Primary |
| 0800731901 | 11.6666 | 40.9259 | 2018-01-13 | South/SE | 11.17 | 17.08 | 16.83 | 6.48 | Primary |
| 0800732001 | 11.9666 | 41.3051 | 2018-01-13 | South/SE | 13.21 | 16.51 | 16.85 | 6.92 | Primary |
| 0800732101 | 12.2198 | 40.9845 | 2018-01-13 | South/SE | 16.31 | 7.99 | 10.33 | 5.90 | Excluded |
| 0800732201 | 12.0572 | 40.7280 | 2018-01-14 | South/SE | 16.01 | 5.92 | 5.88 | 5.58 | Excluded |
| 0800732301 | 11.3785 | 40.1663 | 2018-01-15 | South/SE | 16.73 | 8.39 | 9.74 | 5.29 | Primary |
| 0800732601 | 12.6676 | 40.3879 | 2018-01-16 | South/SE | 23.84 | 3.04 | 2.80 | 5.17 | Excluded |
| 0800732701 | 12.2969 | 39.9949 | 2018-01-17 | South/SE | 24.20 | 5.33 | 2.20 | 5.27 | Excluded |
| 0800732801 | 12.1075 | 39.7012 | 2018-01-17 | South/SE | 26.10 | 6.36 | 7.98 | 4.81 | Primary |
Exposures are cleaned PHA detector exposures, not independent spacecraft time. The 43 MOS products sum to 533.571 ks. “Numerical only” is the MOS2-only field 0800730101; excluded fields fail technical or numerical-quality requirements.
SAS processing provenance
The spectra were derived from the existing project reduction with the
XMM-Newton Science Analysis System (SAS; Gabriel et al. (2004)). We
audited the headers of every v19 input rather than assigning generic
“standard SAS” settings. All 43 source PHA files report
xmmsas_20211130_0941-20.0.0, evselect-3.71.1;
their ARFs report arfgen-1.102.1, and their RMFs report
rmfgen-2.8.5. The event selections stored in every PHA
require PATTERN<=12 and FLAG==0. The
spectra therefore include valid MOS single–quadruple events and reject
flagged events.
The final PHAs contain the applied good-time intervals as embedded
GTI extensions, but the local tree does not retain the upstream
flare-light-curve thresholds or a command log that uniquely reconstructs
those GTIs for the full sample. One partial command record survives for
0800732701, but it cannot define a homogeneous 22-field prescription.
Likewise, the local ccf.cif files are not bound to the
PHA/ARF/RMF products by a recorded identifier or hash, and their file
dates do not uniquely establish which calibration index was active
during each extraction. We therefore report the verified SAS task
versions and final cleaned exposures, but do not invent sample-wide
flare thresholds or a CCF release.
The retained tree is sufficient to reproduce the paper analysis from
the 43 final spectral product groups. It is not sufficient for a bitwise
ODF-to-PHA rebuild: complete ODF material is retained locally for only
three selected observations, and original mos-spectra/QPB
command files are absent for most fields. Recovery of the homogeneous
flare-screening, calibration-index, and extraction-command provenance
remains a submission requirement.
Source masks and extraction footprint
The analysis uses the products under
rpc/4.background/bkg_mos[12]_00500-02000_exclude_extent_source_mask.
The extraction is the usable MOS field of view after detector-coordinate
chip/gap masks and point/extended-source exclusions; it is not a simple
15 arcmin circle. The authoritative mask geometry is embedded in each
PHA as REGION extensions and in the SLCTEXPR selection. The
embedded regions contain negative detector-coordinate ellipses for
excluded sources together with the valid chip and field-of-view
geometry. The reduction tree also retains the associated
edetect_stack source-position and extended-source region
products. The exact source-detection likelihood threshold and the rule
that maps detections to exclusion-ellipse size are not recoverable from
the current command history and must be frozen before submission.
We convert the PHA BACKSCAL value with the production
relation \(\Omega_i=\texttt{BACKSCAL}_i(1/20/60)^2\)
arcmin\(^2\), corresponding to the 0.05
arcsec DETX/DETY sampling. The resulting detector-specific areas span
238.3–565.8 arcmin\(^2\). MOS1 and MOS2
areas are kept separate, so failed CCDs, chip gaps, bad detector
regions, and source holes propagate directly into each model through its
own area factor. Relative to the 706.9 arcmin\(^2\) area of a nominal 15 arcmin circle,
the individual detector footprints cover 33.7–80.0%. The
one-PHA-per-field convention used in the earlier absorption-footprint
audit covers 33.7–52.6% (median 45.2%). The 15 arcmin aperture is used
only for the baseline HI4PI column estimate, not for spectral
extraction.
Observed, QPB, and response spectra
For each retained detector we use four matched products:
*-obj.pi for the observed diffuse-field spectrum,
*-back.pi for the ESAS quiescent-particle background (QPB),
and the corresponding *.arf and *.rmf. The QPB
is generated in the same detector footprint as the observed spectrum. It
represents the particle background and is not a blank-sky or
astrophysical foreground spectrum. In Sherpa, obj.pi is
loaded as the source dataset, back.pi is loaded as its
background with the stored errors, and QPB subtraction is applied before
grouping and fitting. The remaining spectrum therefore contains
celestial emission plus residual detector background and Soft Protons,
which are modeled explicitly below.
The cleaned detector exposures range from 2.20 to 20.57 ks and sum to
533.571 ks over the 43 detector products. The full product paths,
task-version strings, exposures, areas, embedded-mask checks, and
quality classifications are written to the machine-readable file
m31_cgmsum_v19_observation_products.csv. Its generator
reads the frozen production manifest and the actual FITS headers, and
fails unless it finds exactly 22 pointings, 43 detector products, Thin1
in every PHA, and the verified PATTERN/FLAG selections.
Energy filtering and response-aware grouping
The final MOS fitting domain is \[\begin{equation} 0.4\text{--}1.4~\mathrm{keV}\ \cup\ 2.0\text{--}8.0~\mathrm{keV}. \label{eq:fitband} \end{equation}\] The 1.4–2.0 keV interval is excluded because it contains the strong MOS Al K\(\alpha\) and Si K\(\alpha\) instrumental features and is sensitive to residual Soft Protons. Every fit-stage energy selection first resets the broad band and then reapplies this gap; the final Sherpa sessions preserve the two disjoint intervals.
Spectra are grouped after QPB subtraction. Starting from individual channels, a bin is closed only after both (1) the QPB-subtracted signal-to-noise ratio \[\begin{equation} {\rm S/N}=\frac{C_{\rm obj}-C_{\rm QPB}} {\sqrt{C_{\rm obj}+C_{\rm QPB}}} \end{equation}\] reaches 6 and (2) its channel width samples the local RMF FWHM by no more than a factor of 6. Group accumulation is reset at ignored channels, and group starts are forced at 0.4, 1.4, 2.0, and 8.0 keV. Consequently, no fitted bin bridges the excluded MOS interval.
Spectral Analysis
Response construction and model topology
We fit each pointing independently in Sherpa 4.18 (Freeman et al. 2001) using XSPEC model components (Arnaud 1996). Within one pointing, the MOS1 and MOS2 celestial parameters are linked, whereas the residual instrumental-line and Soft Proton parameters remain detector specific. No parameter is linked between different pointings. If \(R_i\) denotes the detector RMF folded with its ARF, \(D_i\) denotes the same redistribution with a unit effective area, and \(\Omega_i\) is the PHA solid angle, the count model for detector \(i\) is \[\begin{equation} \begin{split} M_i(E)={}&R_i\Big\{\Omega_i\big[G_i(E)+S_{\rm LB}(E)+S_{\rm SWCX}(E)\\ &+A(N_{{\rm HI},i})\{S_{\rm CXB}(E)+S_{\mathrm{CGM}_{\mathrm{sum}}}(E)\\ &\qquad\qquad\quad+S_{\rm hot}(E)\}\big]\Big\}\\ &+D_i\{\Omega_i P_i(E)\}. \end{split} \label{eq:model} \end{equation}\] Here \(G_i\) is the residual detector-line model, \(P_i\) is the Soft Proton continuum, and the remaining terms are sky surface-brightness components. Multiplication by \(\Omega_i\) converts every surface-brightness normalization to the actual detector footprint before response folding. The multiplicative sky-background scale present in the implementation is fixed to unity.
The component called MWhalo in the fitting code is
denoted \(\mathrm{CGM}_{\mathrm{sum}}\)
here. It is the cool absorbed APEC component required by the line of
sight, not a spectroscopic assignment to the Milky Way. The historical
target-specific SrcApec, SrcPo, and
Srcabs components remain in the implementation for
compatibility but are set to zero and frozen in this independent
diffuse-field analysis. Thus no additional M31-only spectral component
is fitted on top of \(\mathrm{CGM}_{\mathrm{sum}}\).
Table 2. Baseline Spectral Components and Parameter Treatment
| Component | XSPEC/Sherpa form | Absorption/response | MOS linking | Baseline parameter treatment |
|---|---|---|---|---|
| Local Bubble | apec |
unabsorbed, \(R_i\) | linked | \(kT=0.08\)–0.12 keV free; \(Z=1\) and \(z=0\) fixed; normalization free. |
| \(\mathrm{CGM}_{\mathrm{sum}}\) | phabs*apec |
absorbed, \(R_i\) | linked | \(kT=0.14\)–0.40 keV and normalization free; \(Z=0.3\) and \(z=0\) fixed. |
| Hot thermal | phabs*apec |
absorbed, \(R_i\) | linked | \(kT=0.7\) keV, \(Z=1\), and \(z=0\) fixed; non-negative normalization free. |
| CXB | phabs*pegpwrlw |
absorbed, \(R_i\) | linked | Photon index 1.46 fixed; 0.5–2.0 keV pegged normalization free. |
| SWCX lines | three gaussians |
unabsorbed, \(R_i\) | linked | O VII, O VIII, and Fe-L normalizations fixed to zero. |
| Residual 1.3 keV line | gaussian |
\(R_i\) | independent | Centroid free in 1.2–1.4 keV, \(\sigma=0.01\) keV fixed, normalization free. |
| Al/Si detector lines | two gaussians |
\(R_i\) | independent | Normalizations fixed to zero because 1.4–2.0 keV is excluded. |
| Soft Proton | bknpower |
\(D_i\) (unit ARF) | independent | Break fixed at 3.2 keV; normalization and low-energy slope free; final high-energy slope fixed after an intermediate fit. |
Note. “Linked” means shared between MOS1 and MOS2 within one pointing only. All celestial and detector-continuum normalizations are expressed per unit solid angle before multiplication by the detector-specific \(\Omega_i\).
Ri denotes the full ARF/RMF response; Di denotes the RMF plus unit-ARF response used for residual Soft Protons.
Absorption and atomic configuration
For each pointing, the absorption column is fixed to the uniform-area mean of the HI4PI atomic \(N_{\rm HI}\) map within a 15 arcmin aperture (HI4PI Collaboration 2016). The adopted values span \(4.81\times10^{20}\) to \(7.65\times10^{20}\,\mathrm{cm^{-2}}\). The map has a native angular resolution of 16.2 arcmin; a nominal aperture therefore contains only about 2.38 independent beams. The aperture is centered on the pointing coordinate, has complete valid-map coverage, and is not weighted by the MOS source mask or exposure map.
These inputs are atomic H I columns, not total equivalent X-ray absorbing columns. Molecular gas, dust-to-gas conversion, 21-cm optical-depth corrections, and the mismatch between the circular HI4PI aperture and the irregular MOS footprints are not included in the formal errors. They are retained as controlled absorption-systematics branches.
The XSPEC configuration is fixed to the angr
solar-abundance table, vern photoelectric cross sections,
and APEC/AtomDB 3.0.9 (Foster et al. 2012). The exact
APEC line and continuum files are present in the CIAO installation and
their SHA-256 hashes are stored in the production manifest.
Staged fitting procedure
All fits use the Gehrels chi-square statistic
(chi2gehrels) and the Sherpa Levenberg–Marquardt optimizer.
The staged procedure is identical for every pointing:
In the initial \(0.2\)–\(1.4\cup2.0\)–\(8.0\) keV fit, set the Local Bubble to \(kT=0.1\) keV and \(Z=1\), set \(\mathrm{CGM}_{\mathrm{sum}}\) to \(kT=0.2\) keV and \(Z=0.3\), and fix the hot component to \(kT=0.7\) keV and \(Z=1\). The Local Bubble, \(\mathrm{CGM}_{\mathrm{sum}}\), hot-component, and CXB normalizations are free. The per-detector Soft Proton normalization and both slopes are initially frozen with zero normalization. The three SWCX line normalizations remain zero.
Thaw the Local Bubble and \(\mathrm{CGM}_{\mathrm{sum}}\) temperatures within the bounds in Table 2, keep both abundances fixed, and refit the linked MOS spectra.
Apply the final domain in Equation 2 and refit the thermal, CXB, and residual-line parameters.
For each detector, thaw the broken-power-law normalization and both photon indices while keeping the 3.2 keV break fixed, and refit all active parameters jointly.
Freeze the fitted high-energy Soft Proton index. Refit the low-energy index and normalization with the physical ordering \(0\leq\Gamma_{\rm low}\leq\Gamma_{\rm high}\); if the unconstrained intermediate fit violates the ordering, initialize \(\Gamma_{\rm low}=\Gamma_{\rm high}/2\) before the final fit.
The production worker publishes a session only when Sherpa reports
succeeded=True. Every published sidecar records the fit
statistic, method, number of bins, degrees of freedom, all-stage
filters, input/output hashes, and final parameter states. The current
baseline uses one deterministic starting sequence and local
optimization; multi-start and profile-likelihood checks are
intentionally listed as submission work rather than implied here.
Quality selection and validation
All 22 fields were attempted before the quality mask was applied. The structural validator requires the frozen input detector inventory, unique file paths, exact code/input/output/APEC hashes, final fit success and message, finite fit dimensions, all-stage energy filters, grouping parameters, the adopted frozen \(N_{\rm HI}\), and the expected frozen Al/Si line states. A second validator restores each saved Sherpa session in CIAO and rechecks the actual dataset paths, final fit result, filters, grouping starts and gap, absorption, and forbidden-line parameters. Thus a valid JSON sidecar alone cannot promote a field into the scientific sample.
Twenty fields produced structurally valid Sherpa sessions, and all 20
passed the independent restore test. Eighteen passed the complete
technical gate. Fields 0800731301 and 0800732201 have restorable
sessions but non-positive final degrees of freedom. For 0800732601 and
0800732701, the final Sherpa result has succeeded=False
with an “improper input parameters” message; the worker retained
machine-readable failure records rather than publishing sessions.
A field enters the numerical sample if it passes technical validation, has positive degrees of freedom, has \(Q\geq0.01\), and has finite \(\mathrm{CGM}_{\mathrm{sum}}\) temperature and normalization with positive local errors. The primary sample additionally requires exactly one MOS1 and one MOS2 input. The strict sensitivity sample additionally requires at least ten degrees of freedom. The resulting counts are listed in Table 3.
Among the technically valid fields, 0800730401, 0800731001, and 0800732101 fail the \(Q\geq0.01\) numerical-quality criterion. Field 0800730101 passes the numerical gate but remains outside the primary sample because it has MOS2 only.
Table 3. Sample Accounting
| Category | Number of fields |
|---|---|
| Selected pointings | 22 |
| Restorable production sessions | 20 |
| Complete technical validation | 18 |
| Numerically valid | 15 |
| Dual-MOS primary | 14 |
| Strict sensitivity sample | 12 |
The numerical sample includes MOS2-only 0800730101. The primary sample contains 10 North/NW and 4 South/SE fields.
Derived fluxes and local uncertainties
After restoring each valid session, we extract the fitted \(\mathrm{CGM}_{\mathrm{sum}}\) temperature,
APEC normalization, and their block of the saved Levenberg–Marquardt
covariance matrix. For a band \(b\), a
unit-normalization absorbed phabs*apec model with the
pointing-specific \(N_{\rm HI}\)
defines the conversion \(C_b(kT)\).
Because the fitted APEC normalization is already per arcmin\(^2\), the component surface flux is \[\begin{equation}
F_b=K_{\mathrm{CGM}_{\mathrm{sum}}}C_b(kT_{\mathrm{CGM}_{\mathrm{sum}}}).
\end{equation}\] The conversion is evaluated on a model energy
grid with 0.001 keV spacing from 0.05 to 10 keV, using the same
angr, vern, \(Z=0.3\), \(z=0\), and APEC 3.0.9 configuration as the
fit.
We propagate the saved \(kT\)–normalization covariance with a first-order delta method, \[\begin{equation} \begin{split} \sigma_F^2={}&\left(\frac{\partial F}{\partial kT}\right)^2\sigma_{kT}^2+ \left(\frac{\partial F}{\partial K}\right)^2\sigma_K^2\\ &+2\frac{\partial F}{\partial kT}\frac{\partial F}{\partial K}\, \mathrm{Cov}(kT,K), \end{split} \end{equation}\] where the temperature derivative is evaluated by a centered finite difference with a step no smaller than \(10^{-5}\) keV. These are local Gaussian statistical approximations. They do not include fixed-\(N_{\rm HI}\), foreground, detector-background, model-topology, or calibration uncertainty and are not profile-likelihood confidence intervals.
The primary reported quantity is the absorbed component-only 0.4–1.25 keV surface flux. This band lies wholly inside the directly fitted soft interval. We also calculate 0.5–2.0 keV for comparison with the literature, but label it secondary because 1.4–2.0 keV is supplied by model extrapolation across the excluded MOS interval.
Population summaries and side comparison
For temperature and flux separately, we summarize each side and the full primary sample with a normal random-effects model. For measurements \(x_i\) with local errors \(\sigma_i\), the likelihood uses variances \(\sigma_i^2+\tau^2\); at fixed \(\tau\), the mean is weighted by \((\sigma_i^2+\tau^2)^{-1}\). We determine the non-negative intrinsic scatter \(\tau\) by maximizing the normal likelihood and quote \([\sum_i(\sigma_i^2+\tau^2)^{-1}]^{-1/2}\) as the formal error on the mean. The reported heterogeneity diagnostic is the fixed-effect statistic \[\begin{equation} Q=\sum_i \frac{(x_i-\bar{x})^2}{\sigma_i^2}. \end{equation}\] North-minus-South differences are formed from the two fitted side means, with their formal errors added in quadrature; the ratio and fractional asymmetry are propagated to first order from the same side-mean errors. For the current temperature and flux measurements, the maximum-likelihood extra scatter is zero. This means only that the available local errors do not require additional field-to-field scatter; it does not establish a physically uniform plasma.
The baseline side summaries do not include radius as a covariate. Because the two sides have different radial coverage, they are descriptive diagnostics rather than the final publication likelihood. The required shared radial term plus side offset, boundary-aware intervals, and controlled systematic branches are listed in Section 6.
Results
A common soft spectral component
The 14 primary fields give a common temperature of \[\begin{equation} kT_{\mathrm{CGM}_{\mathrm{sum}}}=0.17531\pm0.00436\,\mathrm{keV}. \end{equation}\] The heterogeneity statistic is \(Q=11.631\) for 13 degrees of freedom, corresponding to \(p=0.558\), and the fitted extra scatter is zero. The separate means are \(0.17469\pm0.00444\,\mathrm{keV}\) in the North/NW and \(0.19314\pm0.02382\,\mathrm{keV}\) in the South/SE. Their difference is \(-0.01845\pm0.02423\,\mathrm{keV}\). The data therefore support a common phenomenological temperature near 0.18 keV but do not establish that all fields contain one spatially uniform plasma. With only 14 fields, four on the South/SE side, and comparatively large South/SE errors, this test has limited power to detect modest intrinsic scatter; the zero maximum-likelihood value is not affirmative evidence for physical uniformity.
North/NW and South/SE surface fluxes
The primary absorbed 0.4–1.25 keV surface fluxes are summarized in Table 4. The North/NW best fit is higher, but the difference is only \(1.01\sigma\). The fractional contrast is also consistent with zero. Figure 1 shows that the two sides have different radial sampling and that several individual uncertainties are large.
Table 4. Primary v19 Measurements with Local-Covariance Uncertainties
| Quantity | North/NW | South/SE | North minus South |
|---|---|---|---|
| kTCGMsum (keV) | 0.17469 ± 0.00444 | 0.19314 ± 0.02382 | −0.01845 ± 0.02423 |
| Absorbed F0.4–1.25 (10⁻¹⁵ erg cm⁻² s⁻¹ arcmin⁻²) | 1.106 ± 0.054 | 0.894 ± 0.203 | 0.212 ± 0.210 |
The flux ratio is 1.238 ± 0.288 and the fractional asymmetry is 0.212 ± 0.230. Errors are propagated local-covariance errors, not profile-likelihood intervals.
The numerically valid sample contains 15 fields after adding the single-MOS sensitivity case. Its North-minus-South flux difference is \(0.198\pm0.210\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), nearly identical to the primary result. The strict 12-field sample gives \(0.147\pm0.274\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). These checks do not turn the mild positive best fit into a significant asymmetry.
Effect of the corrected analysis
For the same 14 fields, the legacy v18 analysis gave a North-minus-South difference of \(0.173\pm0.094\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). The v19 revision changes the effective band, resets grouping across the MOS gap, removes unconstrained forbidden-band line components, and uses pointing-dependent HI4PI columns. The central contrast changes by only \(0.039\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), whereas its formal uncertainty more than doubles. Thus the corrections primarily alter how strongly the contrast is constrained rather than reversing its best-fitting sign, although their individual contributions cannot be separated from this combined regression. The legacy result is retained only as a regression comparison.
Discussion
The directly measured quantity is a line-of-sight sum
The approximately 0.18 keV component is widespread across the sampled footprint and has no detected excess temperature scatter. This consistency does not determine its distance. The Milky Way halo temperature distribution overlaps the current measurement (Henley and Shelton 2013a, 2015; Ponti et al. 2023; Locatelli et al. 2024), and target-specific work toward M31 has explicitly identified a similar combined foreground and target component (Kavanagh et al. 2020). At EPIC resolution, fitting two freely varying APEC components at nearly the same temperature would allow the optimizer to divide one spectral shape arbitrarily. Such a decomposition would not be a physical distance measurement.
Equation 1 is therefore the central interpretation rule. If Local-Group or bridge emission contributes, the observable becomes \[\begin{equation} S_{\rm obs}=S_{\rm MW}+S_{\rm M31}+S_{\rm LG}. \end{equation}\] The additional term increases the physical ambiguity unless it is constrained by an external spatial prior. It is therefore carried only as the explicitly named Local-Group branch of the conditional hierarchy in Section 5.3.
A mild best-fitting contrast, not an asymmetry detection
The North/NW best-fitting flux exceeds the South/SE value, but the \(1.01\sigma\) difference is not compelling. A small observed contrast constrains only \[\begin{equation} \Delta S_{\rm obs}=\Delta S_{\rm MW}+\Delta S_{\rm M31}, \label{eq:contrast} \end{equation}\] with an analogous Local-Group term if required. It constrains M31 asymmetry only if the Milky Way contribution is assumed uniform over the footprint. Conversely, it constrains Milky Way variation only if the sampled M31 halo is assumed symmetric. Neither assumption follows from the current data.
The geometric imbalance is equally important. The primary sample contains 10 North/NW but only four South/SE fields, and their median radii differ by 4.31 kpc. A side-only mean can mix radial structure with azimuthal structure. We therefore treat the values in Table 4 as a descriptive empirical baseline. The next analysis must fit a common radial function and a side offset in one likelihood, with the same nuisance conventions on both sides.
Figure 2 makes the geometric mixing explicit. We define the side-symmetrized total and contrast as \[\begin{equation} T=\frac{F_{\rm N}+F_{\rm S}}{2},\qquad D=F_{\rm N}-F_{\rm S}. \end{equation}\] The primary measurements give \(T=1.000\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\) and \(D=0.212\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). Because the same two side means enter both quantities, their errors are strongly covariant; the confidence ellipses in Figure 2 use the full linear transformation of the two local-covariance errors rather than independent horizontal and vertical error bars.
Projecting the Zhang, Comparat, Ponti, Merloni, Nandra, Haberl, Locatelli, et al. (2024) symmetric M31-mass beta profile through the actual 14-field estimator gives \(D/T=-0.088\): the inverse-variance-weighted North/NW fields lie at larger M31 radii and are predicted to be slightly fainter. In contrast, the disk-like plus spherical Milky Way models of Locatelli et al. (2024), including the HI4PI absorption variation across the same pointings, give \(D/T=+0.058\) to \(+0.063\) because the North/NW fields lie closer to the Galactic plane. The observed positive sign is therefore more compatible with a Milky Way latitude gradient than with dominance by the symmetric M31 template. Quantitatively, however, the residuals from a uniform component, the symmetric-M31 template, and the central Milky Way template are only \(1.01\sigma\), \(1.49\sigma\), and \(0.69\sigma\), respectively, after propagating the \(T\)–\(D\) covariance. The sign is at most a weak discriminator: the present uncertainty excludes neither template nor any non-negative mixture between them. Figure 2 is consequently an exploratory geometry diagnostic, not evidence against a Milky Way or M31 spatial model. Field-to-field cosmic variance in the unresolved CXB can also perturb both the soft-band fluxes and their contrast, and it remains part of the required controlled systematic grid.
Foreground and M31 priors in a two-component constraint plane
Figure 3 temporarily sets \(S_{\rm LG}=0\) and represents the Milky Way and M31 terms by characteristic footprint-averaged values \(x\) and \(y\). Under this explicit approximation, the North/NW and South/SE measurements become two parallel constraints, \(x+y=1.106\pm0.054\) and \(x+y=0.894\pm0.203\), respectively. Their overlap constrains a sum, not a distance decomposition. Horizontal or vertical literature bands become conditional external priors, and an intersection with a diagonal band is an inference under that named prior rather than an independent measurement of either component.
The Henley and Shelton (2013a) continuum survey is retained only as population context and is no longer drawn as a Milky Way interval in Figure 3. Its intrinsic 0.5–2.0 keV range samples many high-latitude directions rather than the M31 footprint: the catalog requires \(|b|>30^\circ\), M31 lies at \(b=-21.57^\circ\), and the nearest detected catalog field is about \(29.7^\circ\) away (Henley and Shelton 2013b). The M31-region oxygen-line catalog entries are observed totals from M31 pointings and are likewise not reused as a pure Milky Way prior (Henley and Shelton 2012a, 2012b).
The beta model summarized by Zhang et al. (2026) does not solve this foreground problem directly. Their illustrative \(\beta=0.4\) and \(r_c=5\) kpc assumptions convert stacked surface-brightness profiles of external Milky-Way-mass galaxies into three-dimensional density profiles as seen by an outside observer. They are informative for the horizontal M31-mass population comparison, but they are not an internal-observer Milky Way geometry and cannot be inserted as the vertical foreground band toward M31.
A Milky Way reference model leaves no non-negative central M31 residual
For the vertical Milky Way prediction, we instead adopt the reference Combined (\(\beta=0.5\)) realization of Locatelli et al. (2024), rather than their combined-plus-SWCX systematic test. The reference model is a spherical \(n_{\rm h}=C r^{-3\beta}\) halo plus an exponential \(n_{\rm d}=n_0\exp(-R/R_h)\exp(-|z|/z_h)\) disk, with \(C=(4.6\pm0.1)\times10^{-2}\), \(n_0=(3.2\pm0.4)\times10^{-2}\,\mathrm{cm^{-3}}\), \(R_h=6.2\pm0.4\) kpc, \(z_h=1.1\pm0.1\) kpc, \(kT=0.15\) keV, and \(Z=0.1\,Z_\odot\). We transform every primary pointing from \((l,b,s)\) to Galactic \((R,z,r)\) and integrate the separately emitting halo and disk terms to 350 kpc.
The absolute conversion requires an extra check because Locatelli et al. fit a 0.614–0.694 keV O VIII-band map after smoothing the APEC spectrum with an approximately 80 eV FWHM kernel, rather than an unsmoothed absorbed broadband energy flux. We therefore apply the same smoothing to the target APEC model before comparing narrow-band emissivities. For \(Z=0.1\,Z_\odot\), the published emissivity is \(1.649\times10^{-16}\,\mathrm{ph\,cm^3\,s^{-1}}\) and the response-matched target value is \(1.736\times10^{-16}\,\mathrm{ph\,cm^3\,s^{-1}}\), giving a closure factor of 0.950. We then apply each pointing’s HI4PI absorption and integrate 0.4–1.25 keV. The resulting reference geometry predicts \(x=1.431\)–\(1.668\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\) across the 14 actual field positions, with \(x_{\rm N}=1.608\), \(x_{\rm S}=1.511\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), and side-balanced \(x=1.560\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). This field-to-field range is the hatched orange band in Figure 3. A dotted bar shows the \(1.507\)–\(1.756\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\) direct-emission-measure/target-APEC mapping before the 5 percent emissivity-closure correction; it is retained as a conversion check. The pale orange envelope, \(1.162\)–\(1.958\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), is only a diagonal first-order propagation of the published Table 1 marginal errors. Because the posterior covariance was not published and Figure D.2 shows correlations, it is not a formal credible interval. Both conversion branches exceed the measured all-field sum across the footprint, so taking the model literally leaves no non-negative central M31 residual; this is a tension among the external O VIII geometry, emissivity convention, and the v19 broadband sum, not evidence for negative M31 emission.
Two transfer assumptions make this comparison conservative but incomplete. First, the \(\mathrm{CGM}_{\mathrm{sum}}\) measurement is obtained with \(Z=0.3\,Z_\odot\), whereas the Locatelli reference model and its O VIII-to-broadband conversion use \(Z=0.1\,Z_\odot\). Both results are expressed as observed-band energy fluxes, so they are not dimensionally incommensurate, but the count-to-flux and line-to-band mappings remain abundance- and temperature-dependent. Holding the fitted O VIII line normalization fixed while changing only the target APEC conversion to \(Z=0.3\,Z_\odot\) gives \(x=1.358\)–\(1.583\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), an all-field estimator value of \(1.520\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), and a side-balanced value of \(1.480\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), about 5 percent below the corresponding reference conversion and still above the measured total. This is a conversion-only sensitivity, not a refit of the v19 spectra; a matched abundance/temperature fitting grid remains required for a complete systematic comparison. Second, the calculation places the entire pointing-dependent HI4PI column in front of the Milky Way component. Real Milky Way emission is distributed through the absorbing material, whereas M31 emission lies behind the Galactic column and may also suffer M31-internal absorption. Replacing the full screen with a distributed absorber would transmit more Milky Way flux for a fixed intrinsic model and would therefore increase, not remove, the \(x>S_{\rm obs}\) tension.
Three interpretations remain open. The all-sky O VIII-fitted Milky Way normalization may overpredict this low-latitude footprint, where the exponential disk term has substantial leverage: at \(b=-21.57^\circ\), one reference vertical scale height (\(z_h=1.1\) kpc) is reached after only about 3.0 kpc along the sightline. The spectral, abundance, or absorption conversion may overpredict the absorbed broadband flux; or the v19 sum may be biased low by residual detector-background, Soft Proton, CXB, or model-topology effects. The present analysis does not adjudicate among them. Its stronger statement is conditional: the reference Locatelli realization, under the documented conversion, leaves no non-negative central M31 residual.
We also audited the X-LEAP products of Pan et al. (2024; Qu et al. 2024). X-LEAP I contains 52 catalog rows within \(4.4^\circ\) of M31, but all are flagged as M31-contaminated and none belongs to the clean sample; the authors further note an M31-associated enhancement extending roughly \(10^\circ\)–\(20^\circ\). Clean 10–20 and 20–30 degree annuli both have a raw detected O VIII median near 1.41 L.U., but those public values are observed line intensities rather than component-separated absorbed broadband fluxes and retain direction-dependent SWCX/absorption issues. At M31’s \(118.78^\circ\) angular distance from the Galactic center, the X-LEAP II empirical radial formula gives a nominal emission measure of \(-0.0011\pm0.0129\,\mathrm{pc\,cm^{-6}}\), demonstrating that it is unconstraining there; that paper explicitly leaves the three-dimensional density decomposition to future work. X-LEAP is therefore retained as a local line/temperature validation data set, not plotted as an independent Figure 3 Milky Way broadband prior.
The empirical horizontal reference comes from the M31-stellar-mass eROSITA stack of Zhang, Comparat, Ponti, Merloni, Nandra, Haberl, Locatelli, et al. (2024). Evaluating its nominal beta profile at 20 kpc and converting its intrinsic 0.5–2.0 keV luminosity surface density with the same APEC template gives \(y=0.942\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). Its intersection with the all-field sum corresponds to a conditional Milky Way residual of \(x\simeq0.150\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). This line is not a formal prior interval: the published marginalized beta-profile errors lack the covariance required for a reliable 10–30 kpc envelope, and the population stack remains sensitive to mass matching and unresolved-source subtraction. The associated scaling and star-forming/quiescent analyses demonstrate the same normalization dependence (Zhang, Comparat, Ponti, Merloni, Nandra, Haberl, Truong, et al. 2024, 2025).
The simulation predictions in Grayson et al. (2025) span a much larger range. Digitizing their 10–30 kpc M31-mass bin and applying the same band conversion gives \(y=0.257\) for EAGLE-AGNdT9, \(0.512\) for fiducial EAGLE, \(1.076\) for EAGLE-NoAGN, and \(1.859\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\) for fiducial SIMBA; the bootstrap bands are \(0.228\)–\(0.281\), \(0.450\)–\(0.568\), \(0.920\)–\(1.194\), and \(1.549\)–\(2.174\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), respectively. SIMBA-NoAGN gives \(y=12.13\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\) and is off scale. Under \(x\geq0\), EAGLE-AGNdT9 and fiducial EAGLE require substantial Milky Way residuals, EAGLE-NoAGN lies close to the full observed sum, and both SIMBA predictions exceed it before any Milky Way term is added. This comparison is potentially discriminating, but remains conditional because one v19 APEC spectrum was used to convert intrinsically multiphase simulation emission.
Conditional complements under named external models
Table 5 collects the central-value bookkeeping implied by the all-field sum \(x+y=1.092\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\). It is not a table of confidence intervals: the literature rows have heterogeneous uncertainty definitions and the v19 value still has local-covariance and systematic limitations. A negative complement means only that the named central external model is incompatible with a non-negative two-component decomposition under the adopted conversion.
Table 5. Conditional Milky Way–M31 Complements at the Central All-Field Sum
| Named scenario | Externally fixed component | Fixed value | Implied complement | Central x, y ≥ 0? |
|---|---|---|---|---|
| Non-negative bookkeeping (U0) | none | – | x + y = 1.092 | continuous feasible set |
| Locatelli reference Combined | MW x | 1.602 | M31 y = −0.510 | no |
| Zhang M31-mass stack | M31 y | 0.942 | MW x = 0.150 | yes |
| EAGLE-AGNdT9 | M31 y | 0.257 | MW x = 0.835 | yes |
| EAGLE fiducial | M31 y | 0.512 | MW x = 0.580 | yes |
| EAGLE-NoAGN | M31 y | 1.076 | MW x = 0.016 | yes, near boundary |
| SIMBA fiducial | M31 y | 1.859 | MW x = −0.767 | no |
| SIMBA-NoAGN | M31 y | 12.13 | MW x = −11.04 | no |
Fluxes are absorbed 0.4–1.25 keV surface brightnesses in 10⁻¹⁵ erg cm⁻² s⁻¹ arcmin⁻². The Locatelli value uses the same all-field inverse-variance estimator as the measured total; its 14-field range is 1.431–1.668. These central complements are conditional model checks, not detections, confidence intervals, or calibrated upper limits.
More massive spiral halos provide useful shape precedents but not a direct M31 normalization (Anderson et al. 2016; Li et al. 2017, 2018). Any inferred M31 residual must retain the name, band conversion, and uncertainty of the foreground and target-profile priors used to obtain it.
A useful conditional hierarchy is:
no external prior: report only the full line-of-sight sum and its spatial contrast;
uniform Milky Way prior: interpret \(\Delta S_{\rm obs}\) as a diluted M31 side contrast;
symmetric M31 prior: interpret \(\Delta S_{\rm obs}\) as Milky Way variation;
external Milky Way distribution: profile a non-negative M31 component and report a prior-labeled upper limit;
external galaxy profile: fit a prior-labeled M31 radial template plus a Milky Way residual.
Local-Group/bridge template: add a named non-negative \(S_{\rm LG}\) term and report how it changes the conditional M31 and Milky Way ranges.
No conditional case should be presented as a direct decomposition.
Analysis Required before Submission
The v19 baseline is sufficiently mature to support the manuscript structure, but five items remain submission blockers.
Matched spatial likelihood. Fit all accepted fields with a common radial dependence plus a North/NW–South/SE offset. Report the joint covariance, likelihood interval for the side offset, and sensitivity to radial parameterization. The raw side means cannot be the final headline statistic.
Boundary-aware inference and optimizer checks. Use multi-start optimization, profile the non-negative thermal normalization and side contrast, and calibrate one-sided coverage with simulations under representative nuisance conditions. Replace every local-covariance interval used for a final claim.
Controlled systematic branches. At minimum vary the mask/exposure-weighted absorption treatment, full-screen versus distributed Milky Way absorption, HI/dust or total-column prescription, Soft Proton model, CXB normalization/cosmic variance, SWCX treatment, abundance and O VIII-to-broadband conversion, hot-component/model topology, and fitting band. Report branch shifts separately before defining any combined uncertainty model.
Quantitative constraint product. Report a calibrated interval on the side offset or fractional asymmetry and at least the conservative U0 upper envelope obtained by assigning the full cool line-of-sight component to M31. A named-prior U1 or smooth-foreground U2 M31 limit would substantially strengthen a full ApJ article, but must remain explicitly conditional.
Submission documentation and figures. Recover the original flare-screening thresholds, bind each extraction to its calibration index, and freeze the source-detection likelihood and mask-radius rule. Quantify any leakage from the bright inner disk/bulge, PSF wings, and unresolved compact sources into the innermost retained fields rather than assuming it is negligible. Add a sky-footprint figure, representative rebinned spectra and residuals, the radial-plus-side fit, and a systematic-branch summary. Archive the current observation table, code, and machine-readable product inventory with persistent versioning.
An absolute M31-versus-Milky-Way decomposition is not a required blocker if the paper remains explicitly about the line-of-sight sum. It becomes a blocker only for claims of a detected M31 halo component or a foreground-subtracted M31 luminosity. A matched off-M31 control field or an external foreground prior would strengthen that second, conditional layer.
Conclusions
We have constructed a homogeneous XMM-Newton baseline for 22 sightlines through the 10–31 kpc projected inner halo of M31. The direct conclusions are:
Fourteen quality-selected dual-MOS primary fields show a cool component with a common phenomenological temperature of \(kT=0.175\pm0.004\,\mathrm{keV}\).
The absorbed 0.4–1.25 keV surface fluxes are \(F_{\rm N}=1.106\pm0.054\) and \(F_{\rm S}=0.894\pm0.203\), in units of \(10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\).
The North-minus-South difference is \(0.212\pm0.210\,10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), so the positive best fit is not a statistically compelling asymmetry detection.
The component is the unresolved line-of-sight sum of Milky Way and M31 emission, with a possible additional Local-Group term. EPIC spectra do not provide a unique distance decomposition.
Under the named Locatelli reference Combined model and documented response-matched conversion, the predicted Milky Way flux exceeds the measured total and leaves no non-negative central M31 residual. This is a model–conversion–data tension, not a direct M31 measurement.
Publication-level spatial constraints require a shared radial-plus-side likelihood, boundary-aware intervals, and a homogeneous systematic grid.
The stable result is therefore not a detection of M31’s hot halo by itself. It is a reproducible measurement of the foreground-degenerate soft component toward M31 and a weak constraint on its large-scale spatial contrast.
Data and Code Availability
The XMM-Newton observations are available from the XMM-Newton Science Archive. The machine-readable measurement table, production manifest, fitting code, validation scripts, and figure-generation script will be archived with a persistent identifier before submission.
This work made use of XMM-Newton data and the HI4PI atomic hydrogen survey. Full acknowledgments, proposal information, and funding statements will be added after the author list is finalized.
