铌酸锂晶体与BICs结合实现高效二次谐波产生

1. 项目背景与核心价值

在非线性光学领域，铌酸锂晶体因其优异的二阶非线性系数和稳定的物理特性，一直是实现高效频率转换的首选材料。而近年来，连续域束缚态（Bound States in the Continuum, BICs）作为一种特殊的电磁模式，因其理论上无限大的品质因子（Q因子）和极强的场局域能力，为非线性光学增强提供了全新思路。

这个项目最吸引我的地方在于将BICs的场增强特性与铌酸锂的二阶非线性效应相结合，通过光子晶体超表面结构实现高效的二次谐波产生（SHG）。传统非线性光学器件往往需要毫米级长度才能获得可观的转换效率，而这种基于超表面的方案有望在亚波长厚度下实现可比拟甚至更高的转换效率——这对集成光子学器件的小型化具有革命性意义。

2. 物理原理与技术路线

2.1 BICs的物理本质与分类

连续域束缚态本质上是一种存在于辐射连续谱中却保持局域化的特殊电磁模式。其形成机制主要分为三类：

1. 对称性保护型BIC：依靠结构对称性阻止能量泄漏
2. 参数空间型BIC：通过调节波矢参数实现模式解耦
3. 拓扑型BIC：源于能带结构的拓扑特性

在光子晶体超表面中，我们主要利用对称性保护型BIC。当结构具有特定的对称性（如C4对称）时，某些模式无法与辐射连续谱耦合，从而形成完美的束缚态。实际加工中不可避免的结构不对称会使BIC退化为准BIC（quasi-BIC），但仍能保持极高的Q因子。

2.2 铌酸锂的非线性耦合机制

铌酸锂（LiNbO₃）作为负单轴晶体，其二阶非线性极化率张量χ⁽²⁾具有以下非零分量：

• χ₁₅ = χ₂₄
• χ₃₁ = χ₃₂
• χ₃₃
• χ₂₂ = -χ₂₁ = -χ₁₆

在COMSOL中模拟时，我们需要通过以下本构关系描述非线性极化：

matlab复制P_i^(2ω) = ε_0 χ_ijk^(2) E_j^ω E_k^ω

其中ω和2ω分别代表基频光和倍频光的角频率。

2.3 超表面结构设计要点

为实现高效SHG，我们的光子晶体超表面需要同时满足：

1. 相位匹配条件：Δk = k₂ω - 2kω ≈ 0
2. 模式重叠优化：最大化基频模与倍频模的空间重叠积分
3. Q因子平衡：基频和倍频模式的Q因子需要满足Q₂ω