4. 光波导与量子模式工程

自发四波混频（SFWM）是集成量子光子芯片中产生纠缠光子对与压缩态的核心三阶非线性过程。它要求泵浦光在波导中同时满足相位匹配与模式重叠，而这两者又直接受波导几何、材料色散与模式分布支配。本节从三种典型波导的模式特性出发，逐步建立群速度色散、非线性系数与量子模式匹配之间的联系，并给出一个面向 SFWM 的色散工程波导设计流程。

一、三种波导的模式特性与适用边界

矩形波导（strip waveguide）由高折射率芯层被低折射率包层包围，横截面通常为数百纳米宽、二百纳米左右高。其模式被强限制在芯内，模场面积可小至 0.1 μm² 量级，因此非线性系数高，但侧壁粗糙度引起的散射损耗也大，且色散对宽度极为敏感。脊波导（rib waveguide）在芯层两侧保留一层薄平板，模式部分渗入平板区，模场面积略大、损耗更低，同时可通过脊高与平板厚度调节色散，工艺容差优于矩形波导。光子晶体波导（photonic crystal waveguide）通过周期孔阵列形成缺陷通道，利用慢光效应在带边附近显著增强非线性，但群速度强烈下降会同时放大损耗与色散，带宽窄，对结构无序极为敏感。

三种波导并非互相替代，而是对应不同的量子光源需求：矩形波导适合宽带、高非线性、短器件；脊波导适合低损耗、可调色散、中等非线性；光子晶体波导适合强慢光增强但窄带的场景。设计 SFWM 波导时，首先要根据目标带宽、损耗预算与工艺平台选择波导类型。

二、群速度色散与相位匹配条件

SFWM 的相位匹配要求泵浦、信号、闲频三波满足动量守恒。在连续波或准连续波近似下，相位失配量可写为 Δβ = 2β_p − β_s − β_i + 2γP_p，其中 β 为传播常数，γ 为非线性系数，P_p 为泵浦功率。将 β 在泵浦频率处作泰勒展开，可得 Δβ ≈ β₂Ω² + (β₄/12)Ω⁴ + 2γP_p，其中 Ω 为信号与泵浦的频率差，β₂ 为群速度色散参数，β₄ 为四阶色散。当 β₂ 与 β₄ 同号且泵浦位于反常色散区时，可在较宽带宽内实现相位匹配；若 β₂ 与 β₄ 异号，则相位匹配带宽被压缩，甚至出现多组解。

群速度色散 D 与 β₂ 的关系为 D = −(2πc/λ²)β₂。在硅基波导中，通过调节波导宽度，可将零色散波长（ZDW）移动至 1550 nm 附近。例如，对于 220 nm 高、500 nm 宽的硅矩形波导，其 ZDW 约在 1550 nm 附近；若将宽度减至 450 nm，ZDW 蓝移至约 1480 nm；若增至 550 nm，ZDW 红移至约 1620 nm。这些数值为典型仿真趋势，实际值需结合具体工艺与包层条件标定。

三、非线性系数与量子模式匹配

非线性系数 γ = 2πn₂/(λA_eff)，其中 n₂ 为非线性折射率，A_eff 为有效模场面积。硅在 1550 nm 处 n₂ ≈ 4.5×10⁻¹⁸ m²/W，若 A_eff = 0.12 μm²，则 γ 可达约 200 W⁻¹m⁻¹ 量级。高 γ 有利于降低泵浦功率阈值，但过强的非线性也会引入自相位调制与双光子吸收，需在设计中权衡。

量子模式匹配要求泵浦、信号、闲频三模在波导横截面上具有足够重叠，且模式的有效折射率满足相位匹配。对于 SFWM，模式重叠积分可写为 Γ = ∫∫ F_p² F_s F_i dA / (√(∫F_p⁴dA ∫F_s²F_i²dA))，其中 F 为归一化模场。若三模均为基模，Γ 接近 1；若信号或闲频落入高阶模，Γ 显著下降，导致转换效率降低。因此，色散工程不仅要调 β₂，还要确保在目标波段内模式保持单模传输。

四、面向 SFWM 的色散工程波导设计示例

以下给出一个基于硅脊波导的 SFWM 设计流程，目标为在 1550 nm 附近产生宽带纠缠光子对。设计参数为示例假设，实际需根据工艺库调整。

步骤 1：确定目标波段与泵浦波长。 设泵浦 λ_p = 1550 nm，信号与闲频分别位于 1530 nm 与 1570 nm，对应 Ω/2π ≈ 2.5 THz。要求相位匹配带宽覆盖 ±5 THz。

步骤 2：选择波导类型与初始尺寸。 选用脊波导，芯层 Si 高度 220 nm，脊宽 500 nm，平板厚度 50 nm。该结构在 1550 nm 附近支持单模，且色散可通过脊宽与平板厚度调节。

步骤 3：扫描脊宽与平板厚度，提取 β₂ 与 β₄。 使用模式求解器计算不同几何下的有效折射率 n_eff(λ)，拟合得到 β₂ 与 β₄。示例扫描结果如下表所示（数值为合理假设，用于说明趋势）。

步骤 4：计算相位失配与匹配带宽。 取脊宽 500 nm、平板 50 nm，β₂ = −0.2 ps²/m，β₄ = 0.01 ps⁴/m，γ = 200 W⁻¹m⁻¹，P_p = 100 mW，A_eff = 0.12 μm²。则 2γP_p ≈ 2×200×0.1 = 40 m⁻¹。在 Ω/2π = 2.5 THz 处，β₂Ω² ≈ −0.2×(2π×2.5×10¹²)² ≈ −4.9×10⁴ m⁻¹，β₄Ω⁴/12 ≈ 0.01×(2π×2.5×10¹²)⁴/12 ≈ 1.6×10⁴ m⁻¹，总 Δβ ≈ −4.9×10⁴ + 1.6×10⁴ + 40 ≈ −3.3×10⁴ m⁻¹。该值较大，说明需要进一步调整几何使 β₂ 更接近零，或利用 β₄ 补偿。

步骤 5：优化几何使 Δβ 在目标带宽内接近零。 将脊宽微调至 520 nm，平板厚度保持 50 nm，重新提取 β₂ ≈ −0.05 ps²/m，β₄ ≈ 0.008 ps⁴/m。此时在 Ω/2π = 2.5 THz 处，β₂Ω² ≈ −1.2×10⁴ m⁻¹，β₄Ω⁴/12 ≈ 1.3×10⁴ m⁻¹，Δβ ≈ 0.1×10⁴ + 40 ≈ 1.04×10³ m⁻¹，已接近匹配。进一步在 ±5 THz 范围内扫描，Δβ 保持在 ±2×10³ m⁻¹ 以内，对应相干长度 L_coh = π/Δβ ≈ 1.5 mm，满足典型器件长度需求。

步骤 6：验证模式匹配与单模条件。 在 1530–1570 nm 范围内检查模式分布，确认泵浦、信号、闲频均为基模，Γ > 0.95。若出现高阶模，需减小脊宽或调整平板厚度。

五、结果判读与迭代方向

完成上述设计后，可通过以下指标判断是否达到 SFWM 目标：相位匹配带宽是否覆盖目标信号/闲频范围；Δβ 在带宽内是否保持在 ±2×10³ m⁻¹ 以内；模式重叠 Γ 是否大于 0.95；有效模场面积与 γ 是否满足泵浦功率预算。若带宽不足，优先微调脊宽使 β₂ 更接近零；若 Γ 下降，检查是否出现高阶模；若损耗过高，可适当增加平板厚度或改用脊波导。整个设计过程以模式求解与色散提取为核心，迭代 3–5 轮即可收敛到合理几何。