[论文解读] The Nonlinear Analytical Envelope Equation in quadratic nonlinear crystals

Morten Bache|arXiv (Cornell University)|Mar 1, 2016

Optical Network Technologies参考文献 1被引用 5

一句话总结

本文推导了超快脉冲在二次非线性晶体中传播的非线性解析包络方程（NAEE），在子载波分辨率下整合了χ⁽²⁾、χ⁽³⁾以及延迟拉曼响应。主要贡献在于提出了一种物理解释清晰的模型，其中非线性项明确源自包络形式，揭示了先前模型因拉曼项归一化错误，导致拉曼效应被低估了2/3倍。

ABSTRACT

We here derive the so-called Nonlinear Analytical Envelope Equation (NAEE) inspired by the work of Conforti et al. [M. Conforti, A. Marini, T. X. Tran, D. Faccio, and F. Biancalana, "Interaction between optical fields and their conjugates in nonlinear media," Opt. Express 21, 31239-31252 (2013)], whose notation we follow. We present a complete model that includes $χ^{(2)}$ terms [M. Conforti, F. Baronio, and C. De Angelis, "Nonlinear envelope equation for broadband optical pulses in quadratic media," Phys. Rev. A 81, 053841 (2010)], $χ^{(3)}$ terms, and then extend the model to delayed Raman effects in the $χ^{(3)}$ term. We therefore get a complete model for ultrafast pulse propagation in quadratic nonlinear crystals similar to the Nonlinear Wave Equation in Frequency domain [H. Guo, X. Zeng, B. Zhou, and M. Bache, "Nonlinear wave equation in frequency domain: accurate modeling of ultrafast interaction in anisotropic nonlinear media," J. Opt. Soc. Am. B 30, 494-504 (2013)], but where the envelope is modelled rather than the electrical field while still keeping a sub-carrier level resolution. The advantage of the envelope formation is that the physical origin of the additional terms that are included to model the physics at the carrier level becomes more clear, in contrast to the electric field equations that are more "black box" expansions of the electrical field. We also point out that by comparing our results to a very similar model and widely used model [G. Genty, P. Kinsler, B. Kibler, and J. M. Dudley, "Nonlinear envelope equation modeling of sub-cycle dynamics and harmonic generation in nonlinear waveguides," Opt. Express 15, 5382-5387 (2007)], the Raman terms presented there will most likely lead to an artificially lower Raman effect.

研究动机与目标

• 开发一种用于超快脉冲在二次非线性晶体中传播的物理解释清晰的模型，同时保持子载波分辨率。
• 在单一解析包络框架内统一χ⁽²⁾、χ⁽³⁾和延迟拉曼响应。
• 阐明包络形式中非线性项的物理起源，避免电场基方程的“黑箱”特性。
• 识别并纠正广泛使用的模型（如Genty等，2007年）中拉曼效应的系统性低估问题。

提出的方法

• 通过使用解析信号形式将非线性波动方程转换到包络域，推导出非线性解析包络方程（NAEE）。
• 引入包络 $ A(z,t) = \tilde{\boldsymbol{\rho}}(z,t) e^{-i\beta_0 z + i\tilde{\boldsymbol{\rho}}_0 t} $，其中 $ \tilde{\boldsymbol{\rho}}_0 $ 为中心频率，$ \beta_0 $ 为传播常数。
• 将非线性极化建模为 $ P_{\text{NL}} = P_{\text{NL}}^{(2)} + P_{\text{NL}}^{(3)} $，其中 $ P_{\text{NL}}^{(2)} $ 为瞬时项，$ P_{\text{NL}}^{(3)} $ 包含电子和拉曼贡献。
• 通过归一化响应函数 $ g(t) = \alpha\delta(t) + (1-\alpha)h_R(t) $ 表达拉曼响应，其中 $ h_R(t) $ 已归一化，且 $ f_R = 1 - \alpha $。
• 应用自陡峭算子 $ \hat{S}_\tau = 1 + i\omega_0^{-1}\partial_\tau $ 和色散算子 $ \hat{D}_\tau $，在共动参考系中推导出最终的NAEE。
• 使用辅助场 $ a(\zeta,\tau) = A(\zeta,\tau)e^{-i\Delta k_{\text{pg}}\zeta} $ 简化相位失配项，实现高效的数值实现。

实验结果

研究问题

• RQ1 如何为包含χ⁽²⁾、χ⁽³⁾和延迟拉曼响应的超快脉冲传播，推导出一个统一的解析包络模型，且具备子载波分辨率？
• RQ2 包络形式中的非线性项的物理起源是什么？与电场基模型中的非线性项有何不同？
• RQ3 为何广泛使用的Genty等（2007年）模型导致拉曼效应被人为降低？
• RQ4 在非线性包络方程中，拉曼项的正确归一化方式是什么，以确保拉曼响应建模的准确性？
• RQ5 在包络方程中，$ e^{i\omega_0\tau} $ 相位因子如何影响拉曼项的结构？

主要发现

• NAEE模型成功在一个统一且物理解释清晰的框架中，实现了χ⁽²⁾、χ⁽³⁾和延迟拉曼响应的整合，且具备子载波分辨率。
• 该模型揭示了Genty等（2007年）模型中的拉曼项因拉曼响应函数归一化错误，系统性地被低估了2/3倍。
• Genty等（2007年）模型中的拉曼贡献使用系数 $ \frac{2}{3}f_R $，这意味着有效拉曼强度低于真实物理值。
• 本文表明，Genty等（2007年）中用于熔融石英的 $ f_R = 0.18 $ 值导致拉曼效应被人为降低，而非物理真实值。
• 正确的拉曼项应使用 $ f_R $ 作为拉曼响应的分数，而非 $ \frac{2}{3}f_R $，这意味着该模型中拉曼效应被低估了33%。
• 推导出的NAEE方程（31）提供了一种数值高效且物理解释准确的超快非线性动力学模拟形式，适用于二次晶体。