NGC 2613 热 halo：现有数据与 XMM-Newton 曝光可行性

AO-26 proposal 技术论证草稿 · 更新于 2026-09-25 · 数据： ObsID 0149160201, v1.7.1/v1.7.2 ESAS MOS

结论（初版）：现有 23 ks MOS 数据在盘两侧 1–3′（6.6–20 kpc）只以 Δχ²≈8.5 看到弱的盘外热成分，温度未约束。按现有亮度，110 ks 申请（约 79 ks clean）能在 N+S 1–3′ 或 1–2′ 得到一个温度（约 ±15–25%），但 2–3′ 即使 300 ks 也难以可靠检出，NGC 2613 做不出温度的垂直剖面。它更适合回答“盘外亮度是否低于预测”。

23.0 ks现有 MOS clean 曝光

9.1×10^-7N+S 1–3′ SrcApec norm/arcmin²（kT=0.3，现有数据）

Δχ² = 8.5现有数据中 SrcApec 显著性

±14%110 ks clean 下 kT 5–95% 分布半宽（N+S 1–3′）

1. 现有 XMM 数据与有效时间

本页只用 ObsID 0149160201 的 MOS 谱。0149160101 受严重背景耀斑影响（proposal 草稿称不可用），未纳入。0149160201 的筛选后 GTI 尚未像 NGC 4565 那样逐段核对，下表给出最终谱 EXPOSURE。

ObsID	仪器	筛选后 GTI (ks)	谱 EXPOSURE (ks)	备注
0149160201	MOS1	—	22.97	本页使用
0149160201	MOS2	—	23.01	本页使用
0149160201	pn	—	17.85	未用于本页
0149160101	MOS1 / MOS2 / pn	—	8.10 / 2.06 / 9.61	严重耀斑，未使用

2. 光学与 X 射线图像

30′×30′ 视场，DSS2 red（CDS hips2fits）与 ObsID 0149160201 的 0.5–2 keV MOS1+MOS2+pn 计数图（未去点源、未做曝光与背景改正）。第二张为 10′ 视场并标出文献中的结构。

NGC 2613：DSS2 red 与 XMM 0.5–2 keV（ObsID 0149160201）。

3. 抽谱区域

盘两侧（N、S）各一组沿盘长 6′ 的垂直 box，点源以 cheese scale=0.3 挖除。N 与 S 同高度的谱联合拟合（各自响应，共享 SrcApec）。BKG 区为 r<14′ 且在 3×D25 椭圆外。

NGC 2613 盘两侧垂直抽谱区域。所有 box 都进入下面的剖面；模拟使用 1–2′、2–3′、1–3′（N+S 联合）。

区域组	区域	离盘高度 (kpc)	有效面积 (arcmin²，去点源后，MOS 平均，多侧相加)	用途
N+S 1–3′	N_1_3 + S_1_3	6.6–19.9	22.0	主方案（宽 halo 区）
N+S 1–2′	N_1_2 + S_1_2	6.6–13.3	10.9	高度分 bin
N+S 2–3′	N_2_3 + S_2_3	13.3–19.9	11.1	高度分 bin
N+S 3–5′	N_3_5 + S_3_5	19.9–33.2	22.7	敏感性对比
N+S 4–6′	N_4_6 + S_4_6	26.5–39.8	20.5	敏感性对比
N+S 5–7′	N_5_7 + S_5_7	33.2–46.4	18.3	敏感性对比
BKG	BKG_r14_out3D25	—	406 / 398 (MOS1/MOS2)	r<14′ 且在 3×D25 椭圆外

4. 垂直面亮度剖面（现有数据）

每个抽谱 box 单独拟合（kT 固定 0.3 keV；用带符号的幅度让 norm 可以为负，接近零的 box 得到可靠误差）。只有 S 1–2′ 有明确的正信号（Δχ²≈11），N 侧与零相容。3′ 以外 S 侧显著为负（S 4–6′：−2.2±0.4 ×10^-6，Δχ²≈27）：当前 BKG 区（r<14′、3×D25 外、两侧）比南侧外围天空更亮。背景在约 2×10^-6 的水平上不均匀，比内侧信号（≈0.9×10^-6）还大，因此 NGC 2613 内侧信号目前不可靠。正在抽谱的最外侧 S 7–10′ 将作为背景（方案 a）来检验这一点。

现有数据中每个抽谱 box 的 SrcApec norm（kT 固定 0.3 keV）。横向误差棒为该 box 的完整高度范围（下轴 arcmin，上轴 kpc）；实心为独立的 1′ bin，空心为重叠的更宽窗口；纵向误差 1σ。右轴为 0.5–2 keV 未吸收面亮度。

Box	高度 (kpc)	面积 (arcmin²，MOS 平均)	谱	每 arcmin² norm (1σ)	Δχ²（1 自由度）	SB 0.5–2 keV (erg/s/cm²/arcmin²)
N 1–2′	6.6–13.3	5.5	MOS1+MOS2	2.6×10^-7 ± 6.4×10^-7	0.2	1.6×10^-16
N 2–3′	13.3–19.9	5.4	MOS1+MOS2	9.9×10^-7 ± 7.1×10^-7	2.0	6.1×10^-16
N 3–4′	19.9–26.5	5.5	MOS1+MOS2	-8.2×10^-7 ± 6.4×10^-7	1.6	-5.1×10^-16
N 4–5′	26.5–33.2	5.4	MOS1+MOS2	5.8×10^-7 ± 7.5×10^-7	0.6	3.6×10^-16
N 5–6′	33.2–39.8	5.4	MOS1	-8.3×10^-7 ± 9.6×10^-7	0.7	-5.1×10^-16
N 1–3′	6.6–19.9	10.9	MOS1+MOS2	6.4×10^-7 ± 4.7×10^-7	1.9	4.0×10^-16
N 2–4′	13.3–26.5	10.9	MOS1+MOS2	2.9×10^-7 ± 4.6×10^-7	0.4	1.8×10^-16
N 3–5′	19.9–33.2	10.8	MOS1+MOS2	-2.7×10^-7 ± 4.7×10^-7	0.3	-1.7×10^-16
N 4–6′	26.5–39.8	9.8	MOS1+MOS2	1.0×10^-7 ± 5.1×10^-7	0.0	6.3×10^-17
N 5–7′	33.2–46.4	7.8	MOS1+MOS2	-3.7×10^-8 ± 6.0×10^-7	0.0	-2.3×10^-17
S 1–2′	6.6–13.3	5.4	MOS1+MOS2	2.4×10^-6 ± 7.1×10^-7	11.2	1.5×10^-15
S 2–3′	13.3–19.9	5.7	MOS1+MOS2	-4.0×10^-8 ± 6.1×10^-7	0.0	-2.5×10^-17
S 3–4′	19.9–26.5	5.9	MOS1+MOS2	2.1×10^-7 ± 6.2×10^-7	0.1	1.3×10^-16
S 4–5′	26.5–33.2	6.0	MOS1+MOS2	-2.0×10^-6 ± 5.4×10^-7	13.6	-1.2×10^-15
S 5–6′	33.2–39.8	5.4	MOS1	-2.2×10^-6 ± 9.5×10^-7	5.5	-1.4×10^-15
S 1–3′	6.6–19.9	11.1	MOS1+MOS2	1.2×10^-6 ± 4.5×10^-7	6.6	7.2×10^-16
S 2–4′	13.3–26.5	11.6	MOS1+MOS2	8.1×10^-8 ± 4.3×10^-7	0.0	5.0×10^-17
S 3–5′	19.9–33.2	11.9	MOS1+MOS2	-8.5×10^-7 ± 4.1×10^-7	4.3	-5.3×10^-16
S 4–6′	26.5–39.8	10.7	MOS1+MOS2	-2.2×10^-6 ± 4.2×10^-7	27.2	-1.4×10^-15
S 5–7′	33.2–46.4	10.4	MOS1+MOS2	-9.3×10^-7 ± 5.0×10^-7	3.5	-5.8×10^-16

5. 背景谱拟合（天空＋仪器）

XSherpa 配方 xartatoms_0900170701_background（profile 状态 STRUCTURAL_PREFLIGHT_ONLY，不声称与 XARTATOMS 数值等价），唯一改动为 Galactic N_H=0.0574×10²² cm^-2（HI4PI）。天空成分 MOS1/MOS2 共享；仪器线与 broken-PL 软质子逐谱独立；net S/N≥6、RMF-FWHM oversampling≤6 分组；0.4–7.2 keV；chi2gehrels；CIAO 4.18 / Sherpa 4.18 / XSPEC 12.14.0k，angr。χ²/dof = 110.4/78；经验质量审计 WARN（见注）。norm 为每 arcmin² 值。

BKG 区 MOS1+MOS2 联合拟合（sb_flux 单位）。数据与总模型用拟合分组；各成分在原生响应网格上绘制；软质子为 detector-space。

成分	kT (keV)	norm（每 arcmin²）
LHB	0.108	2.7×10^-6
MW halo	0.255	2.0×10^-6
CXB (pegpwrlw, Γ=1.4)	—	1.4×10^-3

质量审计 WARN 来自与 XARTATOMS 单一天区经验区间的比较（非先验）：MW halo norm 低于该天区区间但落在 M31CGM 区间内。无参数撞边界。CIAO 4.14 与 4.18 拟合结果一致。

6. 现有数据中的 halo 信号（SrcApec）

与模拟完全相同的估计器：SRC 与 BKG 谱联合拟合，LHB/MW halo/CXB 的 norm 自由（kT 固定为背景拟合值），MOS 仪器线与软质子取背景拟合的每 arcmin² 值，与 QPB 一起作为已知背景扣除；SrcApec 丰度固定 0.3 solar。Δχ² 为去掉 SrcApec 后 χ² 的增加（kT、norm 两个自由度，3σ 对应 11.8）。误差为 90%（covariance）。面亮度为 kT=0.3 keV 时 0.5–2 keV 未吸收能流：norm 1×10^-6/arcmin² ≈ 6.2×10^-16 erg s^-1 cm^-2 arcmin^-2。

区域组	Δχ²	kT 自由拟合 (keV)	norm @kT=0.3（每 arcmin²，90%）	SB 0.5–2 keV (erg/s/cm²/arcmin²)	χ²/dof
N+S 1–3′	8.5	0.87	9.1×10^-7 ± 5.4×10^-7	5.6×10^-16	83.1/86
N+S 1–2′	8.1	未约束（撞边界 3.00）	1.2×10^-6 ± 7.9×10^-7	7.4×10^-16	86.1/81
N+S 2–3′	1.2	未约束（撞边界 0.08）	4.0×10^-7 ± 1.3×10^-6	2.5×10^-16	82.2/80
N+S 3–5′	-0.0	未约束（撞边界 0.08）	0 ± 1.8×10^-7	0	107.1/82
N+S 4–6′	-0.0	未约束（撞边界 0.08）	0 ± 1.3×10^-7	0	135.9/81
N+S 5–7′	-0.6	未约束（撞边界 0.08）	0 ± 2.4×10^-7	0	101.6/80

N+S 1–3′ 与 1–2′ 仅 Δχ²≈8（约 2.5σ），kT 自由拟合未约束；kT 固定 0.3 keV 时 1–3′ norm = (0.9±0.5)×10^-6 arcmin^-2（90%），约为 NGC 4565 同角高度区域的 1/3。2–3′ 与零相容。

N+S 1–3′：现有数据 SRC 谱（扣除已知 QPB＋仪器项），SrcApec kT 固定 0.3 keV；灰线为被扣除的 QPB＋仪器水平。

N+S 1–2′：现有数据 SRC 谱（扣除已知 QPB＋仪器项），SrcApec kT 固定 0.3 keV；灰线为被扣除的 QPB＋仪器水平。

N+S 2–3′：现有数据 SRC 谱（扣除已知 QPB＋仪器项），SrcApec kT 固定 0.3 keV；灰线为被扣除的 QPB＋仪器水平。

7. 模拟方法

真值：天空取背景拟合 anchor；SrcApec（kT、每 arcmin² norm，Z=0.3）；MOS 仪器线＋软质子取背景拟合每 arcmin² 值；每个区域自己的 ESAS QPB 谱。均通过该区域真实 ARF/RMF 折叠。
每次实现：在新曝光 T（MOS1、MOS2 各自的 clean 时间）下对 SRC 与 BKG 区同时 Poisson 抽样；QPB＋仪器项作为已知无噪声背景扣除；xsherpa net S/N≥6 分组；chi2gehrels 联合拟合，SrcApec kT、norm 与三项天空 norm 自由；另拟合无 SrcApec 模型得 Δχ²。
检出：Δχ²≥11.8（2 自由度 3σ），并用 norm=0 模拟核对假阳性率。温度：拟合 kT 的 5–95% 分布半宽 / 真值，以及单次实现 90% 误差在 ±30% 以内的比例。
网格：T = 50、110、200、300 ks（clean）；kT = 0.3 keV（主），0.2 与 0.5 keV 仅在 110 ks；norm = 0.25–5 ×10^-6/arcmin² 以及“现有数据测得值”情形；每格 30–100 次实现。
只用 MOS；加入 pn 预计计数约 ×2.5（余量，尚未模拟）。仪器项与 QPB 视为完全已知，天空背景误差通过联合拟合传播。

8. 模拟结果

曝光 T 为 MOS 的 clean 时间。申请的 110 ks 扣除 SOC 建议的约 40% 背景膨胀后约为 79 ks clean。

已完成 3730 次实现（108 个网格单元）。

以现有数据测得的亮度为真值（kT=0.3 keV）：拟合 kT 的中位数与 5–95% 范围（上）和检出概率（下）随 clean MOS 曝光的变化。虚线为 110 ks 申请扣除 40% 背景膨胀后的约 79 ks clean。

上：检出概率（Δχ²≥11.8）；下：kT 5–95% 分布半宽 / 真值（kT=0.3 keV）。橙色竖线与阴影为现有数据测得的 norm 及其 90% 范围。颜色为 clean MOS 曝光。

按现有数据测得的亮度（真值 kT=0.3 keV）

区域组	clean T (ks)	真值 norm	检出概率	拟合 kT 中位数 [5%, 95%]	kT 5–95% 分布半宽	单次 90% 误差在 ±30% 内	N
N+S 1–3′	50	9.1×10^-7	85%	0.29 [0.23, 0.45]	±37%	33%	60
N+S 1–3′	110	9.1×10^-7	100%	0.29 [0.26, 0.35]	±14%	60%	60
N+S 1–3′	200	9.1×10^-7	100%	0.31 [0.25, 0.37]	±20%	83%	60
N+S 1–3′	300	9.1×10^-7	100%	0.31 [0.27, 0.37]	±17%	92%	60
N+S 1–2′	50	1.2×10^-6	65%	0.28 [0.21, 0.51]	±51%	22%	60
N+S 1–2′	110	1.2×10^-6	98%	0.29 [0.24, 0.36]	±20%	62%	60
N+S 1–2′	200	1.2×10^-6	100%	0.30 [0.27, 0.36]	±15%	77%	60
N+S 1–2′	300	1.2×10^-6	100%	0.30 [0.27, 0.37]	±17%	85%	60
N+S 2–3′	50	4.0×10^-7	0%	0.26 [0.08, 0.97]	±148%	5%	60
N+S 2–3′	110	4.0×10^-7	8%	0.30 [0.14, 0.71]	±96%	10%	60
N+S 2–3′	200	4.0×10^-7	18%	0.31 [0.16, 0.65]	±83%	12%	60
N+S 2–3′	300	4.0×10^-7	37%	0.33 [0.18, 0.61]	±71%	12%	60

所需最低亮度（kT=0.3 keV）

在 norm 网格上对数插值；“—”表示网格上限 5×10^-6 内未达到。

区域组	clean T (ks)	Ω·T (arcmin²·ks)	90% 概率检出的 norm	kT ±30%：50% 实现	kT ±30%：90% 实现	90% 检出对应 SB (erg/s/cm²/arcmin²)
N+S 1–3′	50	1102	1.0×10^-6	1.2×10^-6	2.0×10^-6	6.2×10^-16
N+S 1–3′	110	2424	7.9×10^-7	7.8×10^-7	1.8×10^-6	4.9×10^-16
N+S 1–3′	200	4407	5.0×10^-7	6.5×10^-7	9.6×10^-7	3.1×10^-16
N+S 1–3′	300	6610	4.5×10^-7	5.4×10^-7	8.8×10^-7	2.8×10^-16
N+S 1–2′	50	545	1.8×10^-6	2.0×10^-6	3.4×10^-6	1.1×10^-15
N+S 1–2′	110	1199	9.3×10^-7	1.2×10^-6	1.9×10^-6	5.7×10^-16
N+S 1–2′	200	2180	8.6×10^-7	7.6×10^-7	1.6×10^-6	5.3×10^-16
N+S 1–2′	300	3270	5.9×10^-7	7.7×10^-7	1.6×10^-6	3.7×10^-16
N+S 2–3′	50	557	1.8×10^-6	2.0×10^-6	3.7×10^-6	1.1×10^-15
N+S 2–3′	110	1225	9.2×10^-7	1.1×10^-6	2.0×10^-6	5.7×10^-16
N+S 2–3′	200	2227	8.9×10^-7	9.0×10^-7	1.7×10^-6	5.5×10^-16
N+S 2–3′	300	3340	6.8×10^-7	6.6×10^-7	1.3×10^-6	4.2×10^-16

固定面亮度下，kT 精度随“源区面积 × clean 曝光”的变化。背景主导时精度主要由 Ω·T 决定。

3′ 以外：检出与测温所需亮度随高度的变化

正在用现有外侧 box 计算（已完成 460 / 约 1440 次实现），完成后显示图和表。

假阳性核对：norm=0 模拟中 Δχ²≥11.8 的比例为 0%–0%（期望约 0.3%）。

9. 给 Joel 的回答

现有有效时间：ObsID 0149160201 的 MOS 谱曝光 22.97/23.01 ks、pn 17.85 ks；0149160101 受强耀斑影响未用。
现有数据中的 halo 信号：N+S 1–3′（6.6–19.9 kpc）SrcApec 仅 Δχ²≈8.5；kT=0.3 keV 时 norm = (0.9±0.5)×10^-6 arcmin^-2（90%），约 5.6×10^-16 erg s^-1 cm^-2 arcmin^-2（0.5–2 keV 未吸收）。
温度精度与曝光：以现有亮度为真值，N+S 1–3′ 在 50 ks clean 时检出概率 85%、kT 分布半宽 ±37%；110 ks 时检出 100%、约 ±14%。1–2′ 单独在 110 ks 时约 ±20%。2–3′ 在 300 ks 时检出概率仍只有约 37%。
所需最低亮度（亮度网格，kT=0.3 keV）：110 ks clean 时，N+S 1–3′（≈22 arcmin²）90% 检出需 norm≈0.8×10^-6 arcmin^-2，N+S 1′ bin 需≈0.9×10^-6；一半实现中 kT 在 ±30% 内需≈0.8–1.2×10^-6。50 ks 时分别升至≈1.0–1.8×10^-6 和≈1.2–2.0×10^-6。现有 N+S 1–3′ 亮度（0.9×10^-6）正好在 110 ks 阈值附近，且低于 S 侧约 2×10^-6 的背景不均匀，必须先解决背景问题。网格阈值含约 20% 的蒙特卡罗/插值噪声。
对 proposal 的含义：NGC 2613 能给出一个内侧温度以及亮度剖面/上限，但不能给出温度随高度的变化；若保留它，科学目标应写成“亮度不足是否在更深数据中持续”（Joel 的 eROSITA 差异方向），曝光按检出或排除某一亮度来计算。
与 eROSITA 的比较：本页给出面亮度换算，可与 Zhang 等的堆叠剖面对比；尚未完成该对比。

10. 局限与下一步

只用 MOS；加入 pn 后计数约增加 2.5 倍。
QPB 与仪器项视为完全已知；天空背景不确定度通过联合拟合传播。
单温 APEC、丰度固定 0.3 solar；Li et al. (2006) 的整体谱需要双温（0.08、0.81 keV），真实温度结构可能与 0.3 keV 单温不同，0.2/0.5 keV 情形仍在补充。
距离 22.8±3.3 Mpc（Tully–Fisher），kpc 换算有约 15% 不确定度。
现有亮度只有约 2.5σ，真值可能更低；请以亮度网格结果为准。