[论文解读] Generalized Nonlinear Wave Equation in Frequency Domain

Hairun Guo, Xianglong Zeng+1 more|arXiv (Cornell University)|Jan 8, 2013

Optical Network Technologies参考文献 9被引用 6

一句话总结

本文提出了一种广义的时谐域非线性波方程（NWEF），用于在非线性介质中建模电场动力学，涵盖至三阶极化率的情况，采用慢变谱振幅近似。推导出两种形式：一种用于块体材料中的平面波传播，另一种用于波导，通过频率依赖的等效极化率项纳入空间模场分布及非线性过程（如三波混频、四波混频和拉曼效应）。

ABSTRACT

We interpret the forward Maxwell equation with up to third order induced polarizations and get so called nonlinear wave equation in frequency domain (NWEF), which is based on Maxwell wave equation and using slowly varying spectral amplitude approximation. The NWEF is generalized in concept as it directly describes the electric field dynamics rather than the envelope dynamics and because it concludes most current-interested nonlinear processes such as three-wave mixing, four-wave-mixing and material Raman effects. We give two sets of NWEF, one is a 1+1D equation describing the (approximated) planar wave propagation in nonlinear bulk material and the other corresponds to the propagation in a waveguide structure.

研究动机与目标

• 开发一个完整的时谐域非线性光学波传播建模框架，超越包络近似。
• 在一个统一的波方程中统一描述关键非线性过程——三波混频、四波混频和拉曼效应。
• 通过推导模依赖的非线性波方程，考虑波导中的空间模场分布。
• 利用慢变谱振幅近似，直接引入至三阶的感应极化。
• 通过直接建模电场动力学而非包络动力学，将标准非线性薛定谔方程推广为更一般的形式。

提出的方法

• 通过傅里叶变换从时域麦克斯韦方程推导出时谐域麦克斯韦波方程。
• 使用时域极化率函数将非线性感应极化表示为至三阶项的和，并将其转换为时谐域表达式。
• 应用慢变谱振幅近似（SVSA），将场包络的空间与谱变化视为相对于载波频率缓慢变化。
• 引入空间积分因子（如 $\tilde{\Theta}^{(2)}$、$\tilde{\Theta}^{(3)}$），作为依赖于波导中模场分布的有效非线性极化率。
• 在 SVSA 和平面波假设下，通过简化完整波方程，推导出用于块体材料中平面波传播的 1+1D NWEF。
• 通过将电场投影到导模场分布上，并引入模依赖的非线性耦合项，推导出波导专用的 NWEF。

实验结果

研究问题

• RQ1 如何在不依赖包络近似的情况下，描述非线性介质中电场动力学的完整时谐域行为？
• RQ2 包含至三阶非线性效应并考虑波导中空间模场分布的非线性波方程的形式是什么？
• RQ3 三波混频、四波混频和拉曼效应如何自然地统一于一个时谐域波方程中？
• RQ4 空间模场分布在波导结构中对有效非线性极化率的形成起何作用？
• RQ5 慢变谱振幅近似能否一致地应用于推导广义的时谐域非线性波方程？

主要发现

• 推导出的 1+1D NWEF 能够准确建模非线性块体材料中的平面波传播，直接在时谐域中包含至三阶感应极化。
• 对于波导结构，NWEF 以模特定的非线性耦合系数（$\bar{\Theta}^{(2)}$、$\bar{\Theta}^{(3)}$）表示，其由模场分布上的空间积分导出。
• 二阶非线性项表示为 $\varepsilon_0 \sum_{\alpha_1\alpha_2} \bar{\Theta}_{j;\alpha_1\alpha_2}^{(2)} F[A_{\alpha_1}^\varphi A_{\alpha_2}^\varphi]$，可捕捉和频与差频生成。
• 三阶非线性项同时包含瞬时与拉曼贡献：$\bar{\Theta}_{j;\alpha_1\alpha_2\alpha_3}^{(3)} \left[ (1-f_R) F[A_{\alpha_1}^\varphi A_{\alpha_2}^\varphi A_{\alpha_3}^\varphi] + f_R F[A_{\alpha_3}^\varphi F^{-1}[\tilde{h}_R F[A_{\alpha_1}^\varphi A_{\alpha_2}^\varphi]]] \right]$，支持拉曼效应建模。
• 在慢变模分布近似（SVMDA）下，有效非线性极化率项（$\tilde{\Theta}^{(2)}$、$\tilde{\Theta}^{(3)}$）可近似为常数，从而简化数值实现。
• 最终的波导 NWEF（公式 39）明确将沿传播方向的场演化与非线性耦合项关联，支持在集成光子器件中模拟复杂非线性光学现象。