3. 测量流程分解与仿真需求定义
冷原子干涉仪的测量流程看似是一条连续的时序链,但在仿真建模时,必须把它拆成若干职责清晰、接口明确的模块。拆分的意义不在于把步骤写得更细,而在于让每个模块都有可独立验证的输入、输出和精度指标——只有这样,仿真结果才能与真实装置逐段比对,误差才能定位到具体环节,而不是笼统地说"仿真和实验不一致"。
本节以一套典型的自由落体式冷原子干涉重力仪为对象,把测量流程分解为原子制备、干涉序列、探测与信号处理三大模块,并给出各模块的输入输出定义与仿真精度要求。示例中的参数均为教学假设值,用于说明建模方法,不代表任何实际装置的性能。
一、模块划分与接口定义
分解的第一原则是按物理过程分段、按可观测信号定接口。原子制备模块负责把热原子云冷却、囚禁并释放为具有确定初态和初速度的冷原子团;干涉序列模块负责三脉冲拉曼光作用,使原子波包分束、偏转、合束并积累相位;探测与信号处理模块负责收集荧光或吸收信号,提取跃迁概率,反演干涉相位并换算为重力加速度。三个模块之间传递的不是抽象数据,而是具有明确物理意义的量:原子团的位置、速度、温度、态布居数,以及最终的光强时间序列。
接口定义要具体到"仿真器之间交换什么"。例如原子制备模块的输出应包含:原子团质心位置与速度(三维分量)、位置与速度分布的标准差(等效温度)、基态布居比例、原子数,以及释放时刻。干涉序列模块以此为初值,输出合束后两态的布居差随拉曼失谐或干涉相位的映射关系。探测模块则把布居差转换为带噪声的光电信号。接口一旦确定,各模块就可以独立开发、独立测试。
二、原子制备模块:从热云到确定初态
制备模块的仿真目标,是给出一个"可信的初始条件"。真实装置中,原子先经二维磁光阱预冷却,再进入三维磁光阱装载,随后用偏振梯度冷却把温度降到微开尔文量级,最后关闭囚禁场让原子团自由下落。仿真时不必逐帧模拟每一个光子散射事件,而应抓住对干涉结果影响最大的量:原子团质心位置、质心速度、速度分布宽度(对应等效温度)、原子数,以及两态布居。
一个可操作的建模方式是:把释放时刻的原子团描述为相空间中的高斯分布,位置标准差 σr 与速度标准差 σv 满足 σv = √(kBT/m)。以铷-87 为例,若假设冷却后温度 T = 3 μK,则 σv ≈ 1.7 cm/s;若释放后自由演化 100 ms,位置展宽约 1.7 mm。这个量级直接决定了拉曼光腰斑和探测光斑的尺寸设计,也决定了仿真网格的空间分辨率。精度要求上,制备模块对温度的仿真误差应控制在 5% 以内,初速度误差控制在 1 mm/s 量级,否则后续干涉相位会出现系统性偏移。
三、干涉序列模块:相位积累的核心
干涉序列模块是整个仿真的核心。三脉冲时序为:第一束 π/2 脉冲把原子分束到两个动量态,间隔 T 后 π 脉冲交换两态动量,再间隔 T 后第二束 π/2 脉冲合束。合束后两态布居差随干涉相位振荡,而相位中携带重力信息:φ = keff·g·T²,其中 keff 为有效波矢,T 为脉冲间隔。
仿真这一模块时,需要显式建模三类效应:一是拉曼光与原子团的有限作用时间,它决定脉冲是否为理想 π/2 或 π 脉冲,可用脉冲面积 Θ = Ωeff·τ 描述,Ωeff 为有效拉比频率,τ 为脉宽;二是原子团在脉冲间的自由演化,包括重力下落与波包扩散;三是光场相位噪声与脉冲时序抖动,它们直接转化为相位误差。以 T = 100 ms、keff ≈ 1.6×10⁷ rad/m 估算,g 变化 1 μGal 对应相位变化约 1.6×10⁻⁴ rad,因此仿真中相位误差需压到 10 mrad 以下,才能分辨毫伽量级的重力变化。脉冲时序抖动应控制在 1 μs 以内,否则会引入可观的相位噪声。
四、探测与信号处理模块:从光强到重力值
探测模块把"原子处在哪个态"转换为可读的电信号。常用做法是归一化探测:先测某一态(如 F=2)的荧光强度,再测另一态,用两者之差除以之和得到跃迁概率 P。仿真中需要建模探测光强、光电探测器增益、散粒噪声与背景噪声。若假设原子数 N = 10⁵,单次探测的散粒噪声极限约为 1/√N ≈ 0.3%,对应信噪比约 300;实际装置受背景光与探测器噪声限制,信噪比常在 100 左右。仿真精度要求信噪比不低于 100,采样率不低于 1 MSa/s,以完整记录荧光脉冲波形。
信号处理环节则完成相位反演:对 P 随扫描相位的变化做正弦拟合,提取条纹对比度与相位偏置,再结合 keff、T 换算为 g。这里的关键是条纹对比度——它反映干涉质量,若仿真中对比度明显低于预期,往往说明制备模块的温度或初速设置不合理,或干涉模块的脉冲面积偏离了 π/2。因此对比度不仅是输出指标,也是模块间联合调试的诊断量。
五、仿真需求定义表
把上述分析整理为可执行的需求定义,是本节最直接的交付物。下表给出各模块的输入、输出、关键参数与精度要求,可直接作为仿真器开发与验收的依据。表中数值均为教学示例假设值,实际项目应按装置指标重新标定。
| 模块 | 主要输入 | 主要输出 | 关键参数(示例) | 仿真精度要求 |
|---|---|---|---|---|
| 原子制备 | 光场/磁场时序、装载率 | 位置、速度、温度、态布居、原子数 | T ≈ 3 μK,N ≈ 10⁵ | 温度误差 ≤ 5%,初速误差 ≤ 1 mm/s |
| 干涉序列 | 初态原子团、拉曼光参数 | 干涉相位 φ、跃迁概率 P | T = 100 ms,keff ≈ 1.6×10⁷ rad/m | 相位误差 ≤ 10 mrad,时序抖动 ≤ 1 μs |
| 探测与信号处理 | 态布居、探测光、增益 | 光强时序、条纹对比度、g 值 | 信噪比 ≈ 100,采样率 1 MSa/s | 信噪比 ≥ 100,对比度偏差 ≤ 3% |
需要强调的是,模块划分不是越细越好。把拉曼光路再拆成激光频率稳定、功率稳定、光束准直等子模块,对光学系统设计有价值,但对测量流程仿真而言会引入大量与相位积累无直接关系的细节。合理的粒度应满足:模块边界与可观测信号对齐,模块内部的物理过程可以用少量关键参数刻画,模块之间的耦合通过明确接口传递。按这一原则,三模块划分既能覆盖从冷原子到重力值的完整链路,又保留了足够的可验证性,是后续建立仿真器、分配计算资源和设定验收指标的基础。