4. 原子制备阶段仿真
原子制备是冷原子干涉仪的“源头工序”，它决定了进入干涉序列的原子数、温度与空间密度，也直接决定了后续拉曼脉冲的跃迁效率、干涉条纹对比度以及最终的测量灵敏度。本节把磁光阱（MOT）的装载、冷却与囚禁过程抽象为一组可数值积分的速率方程与动力学方程，通过仿真观察原子数 N(t)、温度 T(t) 与密度 n(t) 的演化，并据此评估制备参数对干涉阶段的影响。
一、把物理过程写成可积分的模型
MOT 的装载本质上是“俘获速率减去损失速率”的竞争。设背景气体碰撞损失率为 Γbg，两体碰撞损失系数为 β，则原子数演化可写为：
dN/dt = R_load − Γ_bg·N − β·n·N
n = N / V_eff, V_eff ≈ (2π)^(3/2)·σ_r²·σ_z其中 R_load 为从背景蒸汽或二维 MOT 注入的装载速率，V_eff 由冷却光束腰斑与磁场梯度决定的椭球体积近似。温度方面，采用简化的多普勒冷却平衡温度模型：T 由冷却光失谐 Δ 与饱和参数 s 决定，T ≈ T_D·(1+s₀/2+2Δ²/Γ²)/(2|Δ|/Γ)，并在装载初期因光压加热而偏高，随失谐优化逐步趋近多普勒极限 T_D = ℏΓ/(2k_B)。密度则由 N 与体积共同决定，压缩阶段（CMOT）通过提高磁场梯度与降低光强使 σ 收缩，密度随之上升。
需要强调的是，上述方程是教学用的简化模型，参数取值属于合理假设的示例，并非实测结果。真实仿真中还应加入亚多普勒冷却、再泵浦光比例、磁场零点漂移等修正项。
二、仿真流程与关键参数
下面给出一个可直接参考的仿真配置示例（参数为示例假设值，用于演示演化趋势）：
| 参数 | 符号 | 示例取值 | 作用 |
|---|---|---|---|
| 装载速率 | R_load | 2×10⁸ s⁻¹ | 决定 N 上升斜率 |
| 背景碰撞损失率 | Γ_bg | 1.0 s⁻¹ | 决定饱和原子数上限 |
| 两体损失系数 | β | 1×10⁻¹⁰ cm³/s | 高密度时限制 N |
| 冷却光失谐 | Δ | −2Γ → −6Γ | 影响温度与俘获率 |
| 饱和参数 | s₀ | 1.5 → 0.5 | 压缩阶段降低光强 |
| 磁场梯度 | ∂B/∂z | 10 → 30 G/cm | 压缩阶段提高密度 |
数值实现上，采用显式欧拉或四阶龙格–库塔法，时间步长取 10 μs 量级，总仿真时长约 2 s（装载 1.5 s + 压缩 0.3 s + 光学黏胶 0.2 s）。每一步先由当前 N 计算密度 n，再更新 N、T，最后根据磁场梯度与光强更新体积，形成闭环迭代。
三、演化过程与结果判读
下图给出仿真输出的典型演化曲线结构：装载段 N 快速上升并趋于饱和，温度在失谐优化后下降；压缩段密度出现明显跃升，温度略有回升；光学黏胶段温度进一步降低但原子数因损失而小幅下降。判读时应重点关注三个量：饱和原子数 Nsat、最终温度 Tf 与峰值密度 npeak。
若仿真显示 Nsat 明显低于预期，应优先检查 Γbg 是否过大（真空度不足）或 R_load 是否受限（冷却光功率、失谐偏离）；若 Tf 偏高，则需检查压缩段光强是否下降过快、失谐是否过近共振；若 npeak 不足，往往是磁场梯度提升幅度不够或体积模型高估了束腰。
四、制备参数如何传导到干涉阶段
制备结果并非孤立指标，而是通过三条路径影响干涉：其一，原子数 N 决定荧光收集信号的信噪比，N 越大条纹对比度越高，但过高密度会引入两体碰撞频移与柱状密度效应；其二，温度 T 决定原子云在飞行时间内的横向膨胀，T 过高会使拉曼光束覆盖不全、有效作用原子比例下降；其三，密度 n 与温度共同决定碰撞频移 δcoll ∝ n·T−1/2，这是干涉相位的主要系统误差来源之一。因此仿真评估时应把 (N, T, n) 作为一组联合指标，而不是单独追求最大原子数。
综上，原子制备阶段仿真的价值在于把“装载快慢、冷却好坏、密度高低”转化为可比较的数值轨迹，并通过参数扫描建立制备条件与干涉性能之间的映射关系，为后续拉曼脉冲时序与探测窗口的设计提供依据。