1. 非线性光学三次谐波产生的物理基础
三次谐波产生（Third Harmonic Generation, THG）是非线性光学中最具代表性的三阶非线性效应之一。当强激光通过非线性介质时，介质极化强度与电场强度的高次项不再可忽略，导致出射光中会出现三倍于入射光频率的新光谱成分。
1.1 三阶非线性极化过程
在传统线性光学中，介质的电极化强度P与电场强度E呈线性关系：
code复制P = ε₀χ⁽¹⁾E
其中χ⁽¹⁾为线性极化率。但当光强达到MW/cm²量级时，必须考虑非线性项：
code复制P = ε₀(χ⁽¹⁾E + χ⁽²⁾E² + χ⁽³⁾E³ + ...)
THG主要来源于χ⁽³⁾项的作用。设入射光电场为E = E₀cos(ωt)，则三阶极化项将产生：
code复制P⁽³⁾ ∝ E₀³cos³(ωt) = (3E₀³/4)cos(ωt) + (E₀³/4)cos(3ωt)
其中3ω项即为三次谐波的来源。
1.2 相位匹配条件
THG转换效率η可表示为：
code复制η ∝ |χ⁽³⁾|²I²L²sinc²(ΔkL/2)
其中Δk = k₃ω - 3kω为相位失配量。实现高效THG需要满足Δk=0的相位匹配条件。常见方法包括：
- 角度相位匹配：利用双折射晶体中o光和e光的折射率差异
- 准相位匹配：周期性极化晶体（如PPLN）通过倒格矢补偿Δk
- 波导模式匹配：通过波导色散调控模式的有效折射率
提示：在COMSOL中计算Δk时，需要准确获取基波与三次谐波的波矢量k值，这要求对材料色散关系有精确建模。
2. COMSOL中THG模拟的建模方法
2.1 多物理场耦合设置
在COMSOL中建立THG模型需要配置以下物理场接口：
- 电磁波，频域（Wave Optics模块）
- 非线性极化材料属性
- 热力学耦合（考虑热致非线性效应时）
关键步骤：
matlab复制% 材料属性设置示例
material1 = mphgetmaterial(model, 'nonlinear_crystal');
material1.property('n', '1.5 + 0.1*(E^2)'); % 非线性折射率
material1.property('chi3', '2e-22'); % 三阶非线性极化率
2.2 几何结构与边界条件
典型THG仿真几何包含：
- 非线性晶体（长度5-20mm）
- 入射端：端口边界（Port）设置基波光源
- 出射端：端口边界检测三次谐波
- 侧面：完美匹配层（PML）吸收散射光
网格划分建议：
- 基波波长范围内至少5个网格单元
- 三次谐波波长需更密网格（建议8-10单元/波长）
- 边界层网格增强近场解析度
2.3 求解器配置技巧
- 频域扫描：
matlab复制study1 = model.study.create('std1');
study1.feature.create('freq', 'Frequency');
study1.feature('freq').set('plist', '[omega, 3*omega]');
- 非线性迭代设置：
- 最大迭代次数：20-30次
- 相对容差：1e-6
- 启用阻尼因子（damping factor）0.7-0.9
- 内存优化：
- 使用对称性简化模型
- 开启几何多重网格（Geometric Multigrid）
- 限制求解频段范围
3. 高转换效率THG系统设计
3.1 晶体参数优化
通过参数化扫描寻找最优配置：
matlab复制model.param.set('L', '[5:1:20]mm'); % 晶体长度扫描
model.param.set('theta', '[0:1:90]deg'); % 相位匹配角扫描
典型优化结果特征：
- 转换效率随晶体长度呈振荡增长
- 存在明确的最优相位匹配角
- 温度变化敏感度约0.1%/℃（需热稳定性控制）
3.2 谐振腔增强技术
法布里-珀罗腔可显著提升有效相互作用长度：
- 输入镜反射率R₁≈99.9%（基波）
- 输出镜反射率R₂≈30%（三次谐波）
- 腔长匹配条件：L = nλ/2
COMSOL中腔体建模要点：
- 使用"电磁波，边界模式"分析
- 设置周期性边界条件模拟往返传播
- 添加损耗项（散射、吸收等）
3.3 波导结构设计
条形波导THG系统优势：
- 光场限制增强功率密度
- 波导色散辅助相位匹配
- 集成化潜力大
设计参数示例：
matlab复制wg_width = 3; % 波导宽度(μm)
wg_height = 1.5; % 波导高度(μm)
delta_n = 0.02; % 折射率对比度
4. 典型问题排查与验证
4.1 收敛性问题处理
常见收敛失败原因及对策：
| 现象 | 可能原因 | 解决方案 |
|---|---|---|
| 解振荡 | 非线性太强 | 降低输入功率分步求解 |
| 发散 | 网格太粗 | 加密谐波频段网格 |
| 伪解 | 初始值不当 | 先求解线性问题作为初始值 |
4.2 实验验证方法
仿真结果验证指标：
- 光谱分析：用OSA检测3ω成分
- 转换效率：功率计测量P₃ω/Pω
- 近场分布：CCD相机记录光斑
误差来源分析：
- 材料参数不确定性（χ⁽³⁾误差可达±20%）
- 表面粗糙度散射（Ra<λ/10）
- 温度梯度影响（ΔT<1℃）
4.3 计算资源优化
大型THG模型加速技巧：
- 使用对称性简化模型
- 频域分解：先低频后高频
- 分布式计算配置：
matlab复制model.sol('sol1').feature('fs1').set('numproc', '4');
内存占用对比：
- 2D模型：~8GB RAM
- 3D全模型：>32GB RAM（建议使用HPC）
5. 进阶应用与扩展
5.1 超表面THG设计
亚波长结构可增强非线性效应：
- 金属纳米结构局域场增强
- 介电超表面相位调控
- 共振模式匹配设计
COMSOL建模要点：
- 使用RF模块的周期性条件
- 添加表面电流密度项
- 启用边界元法（BEM）求解
5.2 脉冲THG模拟
对于飞秒脉冲情况：
- 时域建模：
matlab复制model.study.create('time');
model.study('time').feature.create('time', 'Transient');
model.study('time').feature('time').set('tlist', '0:1e-15:100e-15');
- 关键参数：
- 群速度匹配（GVM）
- 自相位调制（SPM）
- 交叉相位调制（XPM）
5.3 量子关联THG研究
通过二阶关联函数g²(τ)分析：
- 添加量子光学模块
- 定义算符期望值
- 计算强度关联：
matlab复制model.result.numerical.create('g2', 'Global');
model.result.numerical('g2').set('expr', 'emw.g2tau');