ZTF 类星体光变阻尼随机游走参数的高斯过程重估与基线长度偏置的仿真定量检验

1 · 研究问题

用 celerite 类高斯过程（Gaussian process）对数百条 ZTF 类星体光变拟合阻尼随机游走（damped random walk, DRW）参数后，以 ZTF 真实采样模式做的注入仿真能否定量重现 Kozłowski (2017) 的基线偏置判据（特征时标 τ ≲ 10% 基线才可靠）？在该判据筛选后的"合格"子样本上，τ 与光度的标度关系与 MacLeod et al. (2010) 是否一致？

2 · 研究背景与空白

类星体紫外/光学光变可用 DRW（阻尼指数协方差的高斯过程）近似描述，特征时标 τ 被认为联系吸积盘热时标，τ–光度–黑洞质量标度关系（MacLeod et al. 2010，SDSS Stripe 82）是该领域基准。已知的方法学陷阱：Kozłowski (2017) 证明当 τ 接近光变基线长度时拟合系统偏置；后续工作（如 2024 年 arXiv:2409.09637）进一步指出即使 τ < 10% 基线仍可能有可观偏置；Suberlak et al. (2021) 用 SDSS+PS1 延长基线改进参数；2025 年 arXiv:2508.12076 用 SDSS+PS1+ZTF 共 22 年基线研究超出 DRW 的行为。

空白在问题方向：已有工作要么用超长联合基线回避偏置，要么在附录里讨论偏置。把"单一巡天（ZTF，约 6–7 年基线）究竟能可靠回收什么"作为主问题、用真实采样窗口的注入仿真给出定量答案、再在合格子样本上独立检验 MacLeod 标度——这是评价维度转换 + 换样本的组合。风险高于前面各题：数据清洗重、且"合格子样本"可能偏小。

3 · 可检验假设

• H1：ZTF 真实采样注入仿真中，τ 恢复偏置的绝对值随 τ/基线单调上升，在 τ/T = 0.1 处偏置绝对值 > 20%（方向与幅度以仿真为准并与 Kozłowski 2017 的解析预期对照）。
• H2：仅保留仿真判定"可回收"的子样本后，τ–绝对星等关系斜率与 MacLeod et al. (2010) 在 95% 置信区间内一致；若不一致，偏离量与波段/红移分布差异的关系可被仿真解释。

4 · 量化验收标准

1. 方法学校验（硬门槛）：对 200 条合成 DRW 光变（已知 τ、σ，先规则采样后 ZTF 真实采样窗口）做 GP 拟合，规则采样情形恢复参数中位偏差 ≤ 10%。此步不过关，后续全部结论无效。
2. 光变质控口径预先写死：观测点数 ≥ 100、时基 ≥ 4 年、剔除坏点标志；样本选择判据（红移、星等范围）预注册。
3. τ 不确定度双口径：MCMC 后验宽度与 bootstrap 重采样都报告；样本级标度关系拟合用含测量误差的回归（如 LINMIX 类方法）而非普通最小二乘。
4. 全部注入网格（τ/T ∈ [0.01, 1]，10 点对数均匀 × 每点 ≥ 100 实现）预注册。
5. 真实样本 ≥ 300 条（SDSS DR16Q 类星体表仅作目标选择标签，不计入本项目的数据贡献）。
6. 全部脚本与拟合表开源。

5 · 数据与工具

用途	来源 / 工具
光变曲线	ZTF 公开数据释放（IRSA，irsa.ipac.caltech.edu：光变曲线 API 按坐标/编号查询，单条 KB–MB 级；覆盖 2018 至今、g/r 双波段、北天，节奏约 1–3 天）。API 批量查询限速需第 8 周核实
目标选择	SDSS DR16Q 类星体星表（SDSS 官网/CDS）——仅作标签，不计入贡献
GP 拟合	celerite2（内置 O(N) DRW/SHO 核与似然；MCMC 用 emcee；含误差回归需自行实现或用 linmix 移植，需核实）
仿真	自写 DRW 生成器（精确离散递推,约 20 行）+ ZTF 真实时间戳复用
对照关系	MacLeod et al. (2010) 标度系数——仅作对照，不计入贡献

算力：celerite2 单条拟合（含 MCMC）秒到分钟级；300 条真实 + 10³ 条仿真 ≈ 1–2 天后台，可行。降规模版本：真实样本 300 → 100、仿真每点 100 → 50 实现 → H1（仿真偏置曲线）完整保留，H2 统计力下降。

6 · 方法路径

1. 写 DRW 生成器并用解析自协方差验证；装 celerite2/emcee 跑通文档示例。
2. 完成硬门槛校验（规则采样恢复 ≤ 10%）。
3. 用 ZTF 真实时间戳做全网格注入仿真，产出偏置曲面与"可回收"判据（H1）。
4. 从 DR16Q 选样、经 IRSA 拉取光变、按预注册质控过滤。
5. 对真实样本做 GP 拟合，应用可回收判据分割样本。
6. 在合格子样本上做含误差回归，检验 H2 并与 MacLeod 2010 对照。
7. 独立交叉校验：对 30 条光变用结构函数（structure function）法独立估 τ，与 GP 结果秩相关 ρ > 0.7 方视为一致。

7 · 新颖性边界

• 本课题不声称发现 DRW 行为或标度关系（MacLeod 2010 已发表），不声称提出新的偏置理论（Kozłowski 2017 已给出）。
• 已有工作：Kozłowski 2017（偏置判据）、Suberlak 2021（SDSS+PS1 改进）、arXiv:2409.09637（阻尼检验）、arXiv:2508.12076（22 年联合基线、超越 DRW）已覆盖长基线路线。
• 本项目贡献（主结论）：以 ZTF 单巡天为对象、真实采样窗口下偏置的系统仿真定标 + 合格子样本上 MacLeod 关系的独立检验。"单巡天能回收什么"的定量答案对即将到来的大巡天时代（LSST 早期数据同样是短基线）有直接参考价值；H2 一致与否都有效。
• 历届获奖图谱无类星体光变先例。

8 · 决策门槛（go / no-go）

• 第 8 周末：硬门槛通过 + IRSA 批量获取速率核实（≥ 500 条/周可达）。获取受阻 → 具名降级：改用 ZTF 按天区打包的公开光变文件本地切片（数据量升至数十 GB 但一次性），或样本收窄到 100 条。
• 第 20 周末：统计合格子样本规模。若 < 100 条 → 具名降级：课题重心正式移到 H1（纯仿真偏置定标为主结论,真实数据仅作个例演示），仿真框架与统计口径完整保留。
• 第 30 周末：冻结样本；预算裁剪顺序：波段数（g+r → r）→ 结构函数交叉校验样本数 → 仿真网格不可裁（它是主结论底座）。
• 适合学生：统计成熟度较高、能读英文方法论文的队伍；风险在数据获取与子样本规模，两条退路都已具名。

---