The current design当前设计
This design is taken from E032–E045 and the scripts they cite. It separates three physical origins of M31's 0.6–0.8 keV disk emission: old stars (traced by WISE \(W1\)), star formation/ISM (traced by WISE \(W4\) dust), and diffuse hot halo/CGM gas.
此设计取自 E032–E045 及其引用的脚本。它将 M31 0.6–0.8 keV 圆盘发射的三个物理起源分开:旧恒星(由 WISE \(W1\) 追踪)、恒星形成/ISM(由 WISE \(W4\) 尘埃追踪)、以及弥散热晕/CGM 气体。
Variable-W1 elliptical annuli变 W1 椭圆环带
The radial coordinate is an isophotal ellipse built from WISE \(W1\): \(q(a)\) and \(PA(a)\) vary with radius (45 nested boundaries, \(a=4'\)–\(180'\)); beyond \(a=40'\) the geometry is frozen at \(q=0.280\), \(PA=37.7^\circ\) east of north. The accepted extraction annuli are 9 non-overlapping 8' bins from \(a=16'\) to \(a=88'\) (\(a016\_024\) … \(a080\_088\)); \(a088\_096\) (\(a=92'\)) is excluded as unreliable (PI direction 2026-09-10), and \(a056\_064\) (\(a=60'\)) is included with an elevated-rstat caveat (PI direction 2026-09-11).
半径坐标是由 WISE \(W1\) 构建的等光度椭圆:\(q(a)\) 和 \(PA(a)\) 随半径变化(45 个嵌套边界,\(a=4'\)–\(180'\));\(a=40'\) 之后几何冻结于 \(q=0.280\)、\(PA=37.7^\circ\)(北偏东)。接受的提取环带是 9 个不重叠的 8' 箱,从 \(a=16'\) 到 \(a=88'\)(\(a016\_024\) … \(a080\_088\));\(a088\_096\)(\(a=92'\))作为不可靠被排除(PI 指示 2026-09-10),\(a056\_064\)(\(a=60'\))带升高的 rstat 说明被纳入(PI 指示 2026-09-11)。
Animation: the M31 ellipse with annuli and four quadrant sectors动画:带环带和四个象限扇区的 M31 椭圆
MOS-only spectra仅 MOS 谱
Spectra are MOS1/MOS2 only (no pn), extracted with ESAS/SAS on Virgo from the run6 campaign: 33 ObsIDs, 490 complete MOS task bundles (3,920 FITS) across 252 ObsID-region combinations. Each bundle supplies source, background (QPB-subtracted particle background), RMF and ARF. Point-source exclusion uses edetect \(DET\_ML>6\) point-like rows (user decision 2026-07-19). The MOS 1.4–2.0 keV interval (Al K\(\alpha\) + Si K\(\alpha\) + SP) is included in the current campaign (user decision 2026-07-19).
谱仅为 MOS1/MOS2(无 pn),由 Virgo 上的 ESAS/SAS 从 run6 运动提取:33 个 ObsID,490 个完整 MOS 任务包(3,920 个 FITS),横跨 252 个 ObsID-区域组合。每个包提供源谱、背景(QPB 减去的粒子背景)、RMF 和 ARF。点源排除使用 edetect \(DET\_ML>6\) 点状行(用户决定 2026-07-19)。MOS 1.4–2.0 keV 区间(Al K\(\alpha\) + Si K\(\alpha\) + 软质子)包含在当前运动中(用户决定 2026-07-19)。
XSherpa fit template and stagesXSherpa 拟合模板与阶段
The no-background SrcApec fit contract (validated): LB kT=0.1 keV fixed; MWhalo kT free 0.10–0.35 keV, Z=0.3; SrcApec kT free 0.5–1.0 keV (must exceed MWhalo.kT), Z=1; CXB photon index \(\Gamma=1.4\) fixed; fit band 0.4–7.2 keV; abundances angr throughout.
无背景 SrcApec 拟合契约(已验证):LB kT=0.1 keV 固定;MWhalo kT 自由 0.10–0.35 keV,Z=0.3;SrcApec kT 自由 0.5–1.0 keV(必须超过 MWhalo.kT),Z=1;CXB 光子指数 \(\Gamma=1.4\) 固定;拟合带 0.4–7.2 keV;丰度表全程 angr。
scale_skyarea; never Const1D sky area (it adds an energy-bin-width factor and corrupts absolute normalisations).XSherpa 契约(AGENTS.md §8):sky area 是通过 scale_skyarea 的 Python 浮点乘数;绝不要用 Const1D sky area(它会加上一个能量-bin 宽度因子并破坏绝对归一化)。Stages: (1) source/QPB/RMF/ARF extraction QA; (2) XSherpa inspect + spectrum overlay; (3) joint MOS1+MOS2 fit with shared sky and nuisance parameters; (4) conditional \(\Delta\chi^2=1\) errors with the foreground frozen; (5) profile-likelihood errors with the foreground free (required after R008). Vertical errors are conditional \(\Delta\chi^2=1\) with saved nuisance parameters fixed — not profile-likelihood intervals.
阶段:(1)源/QPB/RMF/ARF 提取 QA;(2)XSherpa inspect + 谱叠加;(3)共享 sky 与 nuisance 参数的 MOS1+MOS2 联合拟合;(4)前景冻结时的条件性 \(\Delta\chi^2=1\) 误差;(5)前景自由后的 profile-likelihood 误差(R008 后需要)。垂直误差是条件性 \(\Delta\chi^2=1\),nuisance 参数固定于已保存值——不是 profile-likelihood 区间。
Foreground model前景模型
Three foreground terms per spectrum: Local Bubble (kT=0.1 keV), Milky Way halo (kT=0.25 keV, Z=0.3), and MWhotHalo (kT=0.7 keV). The disk fits freeze foreground absorption at nH=0.03e22 cm\(^{-2}\), but Lit L02/R016 found the MW-only column toward M31 is \(\approx 0.054\)–\(0.067\times10^{22}\) cm\(^{-2}\), about twice the frozen value (08 revision 14). V1's nH bracket 0.048 / 0.0585 / 0.076 stands.
A second, more consequential freeze: the disk template sets MWhotHalo to zero. So any Galactic 0.7 keV emission is not fitted — it lands inside the SrcApec term. That is exactly why the disk soft floor \(s_0 = 3.11\times10^{-16}\) equals the halo-field 0.7 keV term (median \(2.93\times10^{-16}\)) to within 6% (R024), and why the three numbers of I021 (halo-field level, Milky Way star prediction, disk floor) can all be the same foreground.
每谱三个前景项:局部泡(kT=0.1 keV)、银河晕(kT=0.25 keV,Z=0.3)和 MWhotHalo(kT=0.7 keV)。圆盘拟合将前景吸收冻结于 nH=0.03e22 cm\(^{-2}\),但 Lit L02/R016 发现指向 M31 的仅 MW 柱密度为 \(\approx 0.054\)–\(0.067\times10^{22}\) cm\(^{-2}\),约为冻结值两倍(08 修订版 14)。V1 的 nH 括号 0.048 / 0.0585 / 0.076 保持。
还有第二个更要紧的冻结:圆盘模板把 MWhotHalo 设为零。因此任何银河系 0.7 keV 辐射都没有被拟合——它落进了 SrcApec 项内。这正是为什么圆盘软底板 \(s_0 = 3.11\times10^{-16}\) 与晕场 0.7 keV 项(中位数 \(2.93\times10^{-16}\))在 6% 以内相等(R024),也正是为什么 I021 的三个数字(晕场水平、银河系恒星预测、圆盘底板)可以是同一个前景。
Template decomposition模板分解
The spatial model decomposes the 0.7 keV surface brightness into old stellar (\(A_\mathrm{old}\cdot W1\)), young (\(A_\mathrm{young}\cdot W4\)dust), and a constant \(C\). E034 completed old_only (\(\chi^2=11771.817\), dof 7550) and young_only (12541.776); mixed never completed (killed 2026-09-17). E038's identifiability forecast: old+C and old+FUV+C are identifiable (\(\sigma\)(A_young,FUV) 3.3% of A_old as a lower bound); W4dust+mixed is not identifiable (corr(A_young,C)=−1.00) because both WISE templates carry the sky pedestal (B026: W1 off-disk 6.87 DN = 53–97% of the template; W4dust \(\approx\)99% pedestal). V6 rebuilds with net (sky-subtracted) templates.
空间模型将 0.7 keV 表面亮度分解为旧恒星(\(A_\mathrm{old}\cdot W1\))、年轻(\(A_\mathrm{young}\cdot W4\)dust)和常数 \(C\)。E034 完成了 old_only(\(\chi^2=11771.817\),dof 7550)和 young_only(12541.776);mixed 未完成(2026-09-17 被终止)。E038 的可辨识性预测:old+C 和 old+FUV+C 可辨识(\(\sigma\)(A_young,FUV) 为 A_old 的 3.3%,下界);W4dust+mixed 不可辨识(corr(A_young,C)=−1.00),因为两个 WISE 模板都带有天空底座(B026:W1 场外 6.87 DN = 模板的 53–97%;W4dust \(\approx\)99% 底座)。V6 用净(减天空)模板重建。
Annulus \(\times\) quadrant sectors环带 \(\times\) 象限扇区
E041 (V8, PI D21) extracts 5 annuli (\(a032\_048\) … \(a096\_112\)) \(\times\) 4 symmetric 90° quadrants about the M31 centre (PA 37.7°, q=0.28): NE [−7.3°, 82.7°), SE [82.7°, 172.7°), SW [172.7°, 262.7°), NW [262.7°, 352.7°). Expected area shares: NE+SW = 41.3% each, SE+NW = 8.7% each; the four must sum to the annulus within 0.5%. This replaces the earlier sector fits (R011: they were ObsID subsets of full-annulus spectra, 23–72% of area in the other quadrants).
E041(V8,PI D21)提取 5 个环带(\(a032\_048\) … \(a096\_112\))\(\times\) 4 个围绕 M31 中心的对称 90° 象限(PA 37.7°,q=0.28):NE [−7.3°, 82.7°),SE [82.7°, 172.7°),SW [172.7°, 262.7°),NW [262.7°, 352.7°)。预期面积份额:NE+SW 各 41.3%,SE+NW 各 8.7%;四者之和与环带相差在 0.5% 以内。这取代了早先的扇区拟合(R011:它们是完整环带谱的 ObsID 子集,23–72% 面积在其他象限)。
Newest source read: cognition-os/08_RUNNING.md revision 49, updated 2026-11-05 03:56 EST; 07_REVISION_LOG.md R054–R131 and 04_DECISION_LOG.md D26–D30 read in full.
最新读取的来源:cognition-os/08_RUNNING.md 修订版 49,更新于 2026-11-05 03:56 EST;07_REVISION_LOG.md R054–R131 与 04_DECISION_LOG.md D26–D30 已全文读取。
Dual-path extraction (D23/D24)双路径提取(D23/D24)
Every early extraction is made twice (PI D23, 2026-11-01): the native ESAS mos-spectra/
mos_back products (the delivered science product, built with our campaign engine) and DXC's fast_regions
candidate on the same ObsID \(\times\) camera \(\times\) DET region. Each pair is compared channel by channel with DXC's read-only
dualpath comparator (object/QPB spectra, EXPOSURE, BACKSCAL, CCD support, ARF, RMF; verdicts PASS/FAIL/HOLD/
NOT_COMPARABLE/NOT_RUN, failures never dropped). Receipts go to DXC's cumulative table (DXC D099). fast_regions
is developed only in esas_data_reduction, never patched here. V8b (the full sector extraction) starts only after
the pilot pair has been compared and reported.
每次早期提取都做两遍(PI D23,2026-11-01):原生 ESAS mos-spectra/mos_back
产品(交付的科学产品,由我们的运动引擎构建)和 DXC 的 fast_regions 候选对象,作用于相同的 ObsID \(\times\) 相机 \(\times\) DET 区域。
每对用 DXC 的只读 dualpath 比较器逐道比较(对象/QPB 谱、EXPOSURE、BACKSCAL、CCD 支持、ARF、RMF;
裁定 PASS/FAIL/HOLD/NOT_COMPARABLE/NOT_RUN,失败绝不丢弃)。收据转交 DXC 的累积表(DXC D099)。fast_regions
仅在 esas_data_reduction 中开发,从不在此补丁。V8b(完整扇区提取)只有在试点对子被比较并报告后才开始。
D24 (2026-11-02) amends D23: the native side of every dual-path pair, and every new extraction of this project,
uses the current production release esas_data_reduction v1.7.11 (production since DXC D089, 2026-10-29; engine
SHA-256 7899484c…; SAS 20, SAS_RAND_SEED 314159 by default), run end to end from the ODF (calibration, filter, source
detection, background/spectra). The Fornax v1.6.1 fork and the run6 checkpoints it made are no longer the native reference;
run6 products stay valid for the fits already made and serve as a cross-version control. fast_regions candidate
v1.7.12-rc4 is the fast side (built on the v1.7.11 line).
D24(2026-11-02)修订 D23:每对双路径的原生方,以及本项目的每次新提取,都使用当前生产版本
esas_data_reduction v1.7.11(自 DXC D089,2026-10-29 起生产;引擎 SHA-256 7899484c…;SAS 20,默认
SAS_RAND_SEED 314159),从 ODF 到端运行(校准、滤波、源检测、背景/谱)。Fornax v1.6.1 分支及其制作的 run6
检查点不再是原生参考;run6 产品对已完成的拟合保持有效,并作为交叉版本控制。fast_regions 候选
v1.7.12-rc4 是快方(构建于 v1.7.11 线上)。
The pre-registered decision rule预注册决策规则
V7/E040 (DONE, R115): the outer annuli (\(a064\_072\), \(a072\_080\), \(a080\_088\)) were fit under two foreground variants — anchored (MWhotHalo at the v19 level) and released-MWhotHalo — and came back undecided (E040 rule v2): the intercept from SrcApec alone is \(-1.32\sigma\) for HF1 and 0 for HF2/HF3, so a SrcApec boundary at zero with an adequate fit is not evidence that the foreground exceeds the total. The outcome classes are:
V7/E040(已完成,R115):外环带(\(a064\_072\)、\(a072\_080\)、\(a080\_088\))在两个前景变体下拟合—— 锚定(MWhotHalo 在 v19 水平)和释放 MWhotHalo——裁定为未决(E040 规则 v2):仅用 SrcApec 的截距在 HF1 为 \(-1.32\sigma\)、在 HF2/HF3 为 0,因此“SrcApec 在零处触边且拟合充分”并不构成前景超过总量的证据。结果 类别是:
| Outcome结果 | Meaning含义 | Next step下一步 |
|---|---|---|
| AA | Foreground explains the outer rise (obs/pred returns to \(\approx 1\))前景解释外层上升(obs/pred 回到 \(\approx 1\)) | I021 reading stands; the rise is foreground leakageI021 解读成立;上升是前景泄漏 |
| BB | M31 excess survives (\(≥2\sigma\) in two of three outer annuli)M31 超额存留(外层三个环带中两个以上 \(\ge2\sigma\)) | I021 breaks; a real M31 component or a hot halo remainsI021 破裂;一个真实 M31 成分或热晕仍存留 |
| undecided未决 | Neither A nor B at \(2\sigma\)既非 A 也非 B(\(2\sigma\) 处) | V2 decides whether the halo-field 0.7 keV term survives its nullV2 决定晕场 0.7 keV 项能否通过其零假设 |