2、光波导基础理论：麦克斯韦方程组与波动方程、平板波导模式分析、条形波导与脊形波导

光波导，说白了就是光走的“高速公路”。没有它，集成光学芯片就是一堆散落的零件。这一章我们得把地基打牢——从麦克斯韦方程组一路聊到三种最常见的波导结构。我个人习惯是先理解物理图像，再动手仿真，这样遇到问题才不会慌。

2.1 麦克斯韦方程组与波动方程

所有光波导行为的根源，都在这四个方程里。你想想看，光也是一种电磁波，那它怎么传播、怎么被约束，自然得听麦克斯韦的指挥。

在无源、各向同性介质中，麦克斯韦方程组可以写成：

∇ × E = -∂B/∂t
∇ × H = ∂D/∂t
∇ · D = 0
∇ · B = 0

其中 D = εE，B = μH。对第一个方程两边取旋度，再代入第二个，就能得到波动方程：

∇²E - (με) ∂²E/∂t² = 0

嗯，这里要注意——这个方程假设介质是均匀的。实际波导有折射率变化，处理起来会更复杂一些。

2.2 平板波导模式分析

平板波导是最简单的结构——三层介质叠在一起，中间是芯层，上下是包层。光在芯层里全反射传播。为什么它能约束光？因为芯层折射率比包层高。

模式分析的核心，就是求解波动方程在分层介质中的本征解。对于TE模式（电场平行于界面），电场分布可以写成：

E_y(x) = A cos(κx + φ)   (芯层内)
E_y(x) = B exp(-γ|x|)    (包层内)

其中 κ 是横向波数，γ 是衰减系数。它们的关系由特征方程决定：

κd = mπ + 2 arctan(γ/κ)

m = 0, 1, 2, ... 就是模式阶数。m=0是基模，能量最集中。

V参数（归一化频率）是个关键指标：

V = (2π/λ) · d · √(n₁² - n₂²)

V < π/2 时，波导只支持基模。我设计单模波导时，一定会先算这个值。

2.3 条形波导与脊形波导

平板波导只在垂直方向约束光，水平方向是自由的。实际芯片需要二维约束——条形波导和脊形波导就是干这个的。

条形波导：芯层被完全刻蚀，四周都是包层。约束最强，但侧壁粗糙度会导致散射损耗。

脊形波导：只刻蚀一部分芯层，留下一个脊。约束稍弱，但损耗更低，工艺容差更大。

这两种波导没有解析解，得用数值方法。我常用的方法是有效折射率法（EIM）——把二维问题拆成两个一维问题。

# 有效折射率法计算条形波导模式
import numpy as np

def effective_index_method(n_core, n_clad, w, h, lam):
    """
    条形波导有效折射率计算
    n_core: 芯层折射率
    n_clad: 包层折射率
    w: 波导宽度
    h: 波导高度
    lam: 波长
    """
    # 第一步：计算垂直方向有效折射率
    # 假设平板波导厚度为h
    V_vert = (2*np.pi/lam) * h * np.sqrt(n_core**2 - n_clad**2)
    # 基模有效折射率（近似）
    n_eff_vert = np.sqrt(n_core**2 - (lam/(2*h))**2)
    
    # 第二步：计算水平方向模式
    # 将垂直有效折射率作为水平方向的"芯层折射率"
    n_eff_horiz = np.sqrt(n_eff_vert**2 - (lam/(2*w))**2)
    
    return n_eff_horiz

# 示例：硅波导，1550nm
n_eff = effective_index_method(3.48, 1.44, 0.5e-6, 0.22e-6, 1.55e-6)
print(f"有效折射率: {n_eff:.4f}")

我个人做硅光芯片时，90%的情况用条形波导。只有在需要低损耗或大模场时才换脊形。你想想看，硅和二氧化硅折射率差接近2，条形波导能把光约束在亚微米尺度——这是硅光能实现高集成度的根本原因。

下面这张图展示了三种波导的结构和模式分布：

最后说一句——波导设计没有银弹。我见过有人死磕条形波导的侧壁损耗，折腾了三个月，最后换成脊形波导一周就搞定了。关键是要根据你的应用场景选对结构，然后针对性地优化。