4. 拉曼激光系统强度与频率噪声预算

在拉曼型原子干涉重力仪中,拉曼激光同时承担三件事:选态、分束与合束、以及把原子内态相位差“印”到干涉条纹上。因此激光的强度起伏和频率抖动都会沿一条明确的传递链进入干涉相位。本节的目标不是泛泛讨论激光噪声,而是把强度噪声与频率噪声分别折算成等效相位噪声,给出可代入典型参数的估算方法,并说明如何判断哪一项在预算中占主导。

先约定符号与参照:设拉曼作用时间(脉冲间隔)为 T,有效波矢 keff,拉曼脉冲为 π/2–π–π/2 序列,总干涉时间 2T。相位噪声统一折算到干涉输出相位 φ,功率谱密度记为 Sφ(f),单边。灵敏度以每 √Hz 的相位噪声表示,最终换算为重力加速度噪声时除以 (keffT2)。

一、强度噪声的传递函数

强度噪声进入相位的机制主要有两条。第一条是光频移与有效拉比频率抖动:拉曼脉冲的翻转角 θ = Ωeffτ 正比于两束拉曼光场振幅乘积,强度相对起伏 ε(t)=δI/I 直接造成脉冲面积误差。对 π/2 脉冲,面积误差 δθ 引起的相位误差近似为 δφ ≈ δθ·sin(φ0),在条纹斜率最大处(φ0=π/2)退化为 δφ ≈ δθ。第二条是差分光频移:两束拉曼光强度不平衡时,AC 斯塔克频移差随强度起伏,等效为频率噪声,再经 2T 积分放大。

把脉冲序列展开,强度噪声到相位的传递可写成对三个脉冲时刻强度起伏的加权和。对低频(f ≪ 1/T)分量,权重近似为:首末 π/2 脉冲贡献 +1 与 −1,中间 π 脉冲贡献 −2(因 π 脉冲对相位的作用是两倍)。于是低频段传递函数模平方近似为

|H_I(f)|² ≈ (1² + 2² + 1²) = 6   (f ≪ 1/T,白噪声近似)
S_φ,I(f) ≈ 6 · (δθ/θ)² · S_ε(f) / (2T)   (单边,ε 为相对强度噪声)

这里 Sε 是相对强度噪声的单边功率谱密度。若激光相对强度噪声在相关带宽内为白噪声,其与 RIN 的关系为 Sε(f) = RIN(单位 1/Hz,RIN 以 dBc/Hz 给出时先换算)。对高频分量,三个脉冲的加权不再同相,传递函数随 f 出现零点与旁瓣,因此强度噪声预算必须说明所取频段,通常取 0.1 Hz 到 1/(2T) 的积分带宽。

二、频率噪声的传递函数

拉曼光频率抖动 δν 通过两种方式进入相位。其一是直接相位累积:原子在两次脉冲之间积累的相位差正比于拉曼频率失谐与时间的乘积,频率噪声在干涉时间内积分,等效相位噪声为 δφ = 2π·δν·(有效积分时间)。对 π/2–π–π/2 序列,低频频率噪声的传递函数模平方近似为

|H_ν(f)|² ≈ (2π)² · [ (2T)² · sinc²(πf·2T) ]   (低频极限 → (4πT)²)
S_φ,ν(f) ≈ (2π)² · S_ν(f) · (2T)² · sinc²(πf·2T)

其二是差分光频移与光频移耦合:单束光频率变化会改变 AC 斯塔克频移,若两束拉曼光频率噪声不共模,残余差分频移会叠加到有效拉曼频率上。工程上常把两束光锁定到同一参考(如同一台超稳激光或同一光频梳),使共模频率噪声在差分中抵消,剩余的是参考锁定带宽内的残余噪声。

需要注意,频率噪声的传递在低频是“积分型”的,即 Sφ ∝ Sν·T2,这意味着拉曼作用时间越长,频率噪声越被放大;而强度噪声的传递对 T 的依赖较弱(主要经脉冲面积误差)。这一差异是判断主导噪声源的关键。

强度噪声与频率噪声到干涉相位的传递链 激光相对强度噪声 RIN (dBc/Hz) 脉冲面积抖动 δθ/θ = δI/I 三脉冲加权 +1, −2, +1 相位噪声 Sφ,I ≈ 6·(δθ/θ)²·Sε/(2T) 拉曼光频率噪声 Sν (Hz²/Hz) 干涉时间积分 δφ = 2π·δν·2T 差分光频移残余 共模锁定抵消 相位噪声 Sφ,ν ≈ (4πT)²·Sν 典型参数下的量级估算(示例假设,非实测) 强度噪声:RIN = −120 dBc/Hz,T = 50 ms Sφ,I ≈ 1.2×10⁻⁴ rad²/Hz 频率噪声:Sν = 1 Hz²/Hz,T = 50 ms Sφ,ν ≈ 3.9×10⁻² rad²/Hz 结论:该示例中频率噪声经 2T 积分放大后主导,需优先压缩拉曼光频率噪声或缩短 T。

三、典型参数下的量级估算

下面给出一组教学示例参数(属合理假设,非实测):T = 50 ms,keff ≈ 1.6×107 rad/m,拉曼光相对强度噪声 RIN = −120 dBc/Hz(即 Sε ≈ 1×10−12 /Hz),拉曼光频率噪声 Sν = 1 Hz²/Hz(单边,白噪声近似)。

强度噪声:取 δθ/θ ≈ 0.5·√Sε(脉冲面积对强度起伏的响应系数约 0.5),则 Sφ,I ≈ 6·(0.5)²·1×10−12/(2×0.05) ≈ 1.5×10−11 rad²/Hz,对应相位噪声约 3.9×10−6 rad/√Hz。这个量级相当小。

频率噪声:Sφ,ν ≈ (4π×0.05)²·1 ≈ (0.628)² ≈ 0.39 rad²/Hz,对应约 0.63 rad/√Hz。两者相差约五个数量级。换算成重力噪声,δg/√Hz = √Sφ/(keffT2),频率噪声项约 0.63/(1.6×107×2.5×10−3) ≈ 1.6×10−5 m/s²/√Hz,即约 1.6 mGal/√Hz 量级——这显然偏大,说明示例中的频率噪声参数偏保守,实际系统通常需要把 Sν 压到 10−4 Hz²/Hz 以下,或采用共模锁定与差分测量。

估算时务必区分单边/双边谱密度、dBc/Hz 与线性 RIN 的换算,以及传递函数所适用的频段。低频“积分型”频率噪声传递与高频“零点型”强度噪声传递不能混用同一套权重。

四、预算表的组织与判读

把上述两项与探测噪声、振动噪声并列,可形成一张相位噪声预算表。表中每一行应给出:噪声源、谱密度参数、传递函数、折算相位噪声、以及占总噪声的百分比。判读规则很简单:若某一项折算后比其余项之和高出一个量级以上,它就是当前灵敏度瓶颈;若两项接近,则需同时优化,且要注意它们在积分带宽内的形状不同,不能只比较单频点数值。

噪声源输入谱密度(示例)传递关系折算 Sφ (rad²/Hz)占比
拉曼强度噪声RIN = −120 dBc/Hz6·(δθ/θ)²/(2T)≈1.5×10⁻¹¹<0.1%
拉曼频率噪声Sν = 1 Hz²/Hz(4πT)²≈0.39>99%
探测噪声(参考)——待本节外补充—

实际工作中,频率噪声往往不是白噪声,而是有 1/f 或特定线谱(如电源谐波、参考腔振动耦合)。此时不能直接用 (4πT)² 乘一个常数,而应把 Sν(f) 与 |Hν(f)|² 逐频点相乘后积分。一个实用的做法是:先在 0.1–10 Hz 内取若干频点,分别计算强度项与频率项的贡献,画出两条折算相位噪声曲线,观察交叉频率。若交叉点落在探测带宽内,说明两项都不可忽略,需要分别压制。

最后提醒一个常见误区:把激光频率噪声直接当作“频率不稳定度”代入,而忽略拉曼两束光的共模抑制。若两束光来自同一激光经声光调制器频移产生,则源激光的频率噪声在差分中大部分抵消,剩余的是调制器驱动噪声与路径非共模部分。预算时应按“共模抑制比”折减,而不是把源激光的 Sν 全额计入。这一折减系数通常需要单独测量或由锁定环路带宽估算,属于预算中不确定性最大的一项,应明确标注为假设值。