3. 原子源与初态制备的噪声贡献

在量子重力仪的灵敏度预算中，原子源与初态制备环节常被误认为"只要原子够多就行"，实际上它同时决定了干涉条纹的对比度上限和单次测量的信噪比基线。本节把这一环节拆成三条可量化的噪声通道：原子数目涨落、温度（速度）分布、态制备效率，并给出把它们折算进总预算的具体做法。

三条噪声通道的物理来源

原子数目涨落直接体现为探测信号的散粒噪声与经典涨落。若单次参与干涉的原子数为 N，理想情况下荧光或吸收探测的散粒噪声极限为 1/√N，对应相位灵敏度 σ_φ,SN ≈ 1/(C√N)，其中 C 为条纹对比度。但真实原子源的 N 本身在逐次抛射间波动，相对涨落 σ_N/N 会以经典噪声形式叠加，使有效信噪比降为 1/√(N·(1+(σ_N/N)²N)) 量级——当 σ_N/N 超过约 1/√N 时，经典涨落开始主导，继续增加平均原子数的收益迅速饱和。

温度分布通过两个途径进入预算。其一是纵向速度展宽导致干涉时间与有效拉曼脉冲间隔的离散，使不同速度群积累不同相位，条纹被"抹平"，对比度按速度分布的傅里叶分量衰减；其二是横向速度展宽造成原子云在飞行中膨胀、与探测光束模式失配，降低收集效率。对典型三维磁光阱后经亚多普勒冷却的样品，纵向温度从 10 μK 降到 2 μK，可使对比度损失从约 20% 降到 5% 以内。

态制备效率包括原子在 |F=2, m_F=0⟩ 等磁不敏感态的布居比例、态制备后残余的杂态原子比例，以及制备脉冲本身的保真度。杂态原子不参与目标干涉跃迁，却贡献背景荧光，等效于把对比度乘以制备纯度 η。若 η=0.85，则即便其他环节完美，对比度上限也只有 0.85。

把参数折算成预算条目

预算表的关键不是罗列参数，而是把每个参数换算成"相位噪声贡献"这一统一量纲。以一次合理的教学示例设定为例（以下数值为示例假设，非实测）：目标态纯度 η=0.85，纵向温度 T_z=3 μK 对应速度分布衰减因子 f(T_z)≈0.95，横向模式失配因子 0.98，则对比度上限 C_max=0.85×0.95×0.98≈0.79。若平均原子数 N=1×10⁶，散粒噪声项 σ_φ,SN=1/(0.79×1000)≈1.27 mrad。再设逐次原子数相对涨落 σ_N/N=0.02，则经典项 σ_φ,cl≈(σ_N/N)/√N×修正系数，量级约 0.6 mrad，两项平方和开根得 σ_φ,原子源≈1.4 mrad。这个数字就可以直接填入总预算表，与拉曼光相位噪声、振动噪声、探测噪声并列比较。

实际操作中建议按下面的顺序做，避免"先测总噪声再反推"的盲目性：

1. 固定态制备参数并测纯度。用微波或拉曼选态后，通过飞行时间荧光或态选择性推光测 |F=2,m_F=0⟩ 的布居比例，重复 20～30 次取均值与标准差，得到 η 及其稳定性。
2. 测温度而非"估计温度"。用飞行时间法（释放后两个时刻的云尺寸）或拉曼光谱测速度分布，得到 T_z 与 σ_v⊥，再代入对比度衰减模型。
3. 统计原子数涨落。连续记录 50～100 次抛射的探测总计数，计算均值与相对标准差，注意区分真实原子数涨落与探测增益漂移。
4. 合成预算条目。把 η、f(T_z)、σ_N/N 分别折算为对比度因子与相位噪声项，按平方和合成，标注各项占比。

参数敏感度与优化方向

三条通道的敏感度并不相同。对比度对 η 是线性关系，η 每降 1%，C 就降约 1%，因此态制备纯度的提升收益直接且可预测；对温度则呈非线性，当 T_z 从 5 μK 降到 2 μK 时对比度改善明显，但继续降到 1 μK 以下收益递减，此时应转向抑制横向膨胀。原子数涨落的影响则与 N 耦合：N 越大，同样的相对涨落造成的经典项绝对贡献越小，但散粒噪声项也越小，两者的相对占比需要每次重新核算。

一个常见的预算失衡是：把大量精力放在提高 N 上，却忽略了 σ_N/N 随 N 增大而同步恶化（例如 MOT 装载饱和、推光效率波动）。此时总信噪比曲线会出现"平台"，继续增加平均原子数不再改善灵敏度。识别这一平台的方法很简单：在预算表中把 N 作为横轴、σ_φ,原子源 作为纵轴画一条曲线，若曲线在目标 N 附近已趋平，说明瓶颈在涨落而非原子数。

完成上述折算后，原子源与初态制备这一条目就以"对比度因子 + 相位噪声项"的双重身份进入总预算：对比度因子影响所有下游噪声项的换算系数，相位噪声项则与其他独立噪声源平方和合成。把这两部分分开列账，才能在后续排查中快速判断——当总灵敏度不达标时，究竟是原子源本身贡献过大，还是它通过降低对比度放大了其他环节的噪声。