🔬 Curso de Fotônica Integrada Módulo 04 • Teoria de Modos Acoplados (CMT) • Yariv • Lorentz • Bragg

Teoria de Modos Acoplados em Fotônica Integrada

Das origens em micro-ondas de Pierce à unificação de Yariv em óptica guiada: dedução rigorosa via reciprocidade de Lorentz, dinâmicas co e contra-direcionais, formalismo matricial e aplicações em PICs.

Autor: Prof. Marcelo Eduardo Vieira Segatto • IN-Foton & UFES
📥 Baixar Slides Oficiais (PDF Beamer) 🧪 Laboratório Virtual de Simulação CMT

1. Evolução Histórica da Teoria de Modos Acoplados

A Teoria de Modos Acoplados (Coupled-Mode Theory — CMT) constitui uma das ferramentas mais universais e elegantes da física aplicada e engenharia eletromagnética. Em vez de resolver as equações completas de Maxwell para geometrias complexas e heterogêneas, a CMT modela o sistema como uma coleção de subsistemas ideais desacoplados, cuja interação contínua transfere energia mediante perturbações espaciais ou temporais.

1954 • O Ponto de Partida com John R. Pierce:
Pierce publicou no Journal of Applied Physics o artigo pioneiro "Coupling of Modes of Propagation", analisando a troca de energia entre feixes eletrônicos modulados e ondas eletromagnéticas lentas em tubos de ondas progressivas (TWTs). Pierce demonstrou que, ao igualar as velocidades de fase (condição de sincronismo ou phase-matching), a potência flui continuamente do feixe para o campo guiado, amplificando sinais de micro-ondas.

Nos anos seguintes, a teoria expandiu-se rapidamente:

  • Stewart E. Miller (1954): Transpôs os conceitos para linhas de transmissão de micro-ondas multiguiadas nos Laboratórios Bell, derivando a física fundamental do acoplador direcional passivo e divisores de potência.
  • Dietrich Marcuse (1974): Sistematizou o formalismo estatístico de modos acoplados em fibras ópticas e guias dielétricos com perturbações contínuas, demonstrando como a rugosidade superficial (sidewall roughness) acopla modos confinados ao contínuo de radiação (perdas por espalhamento).
  • Amnon Yariv (1973): No célebre artigo "Coupled-Mode Theory for Guided-Wave Optics" (IEEE J. Quantum Electron.), Yariv unificou o formalismo da CMT para fotônica integrada, lasers de realimentação distribuída (DFB), grades de Bragg e interações ópticas não lineares ($\chi^{(2)}$ e $\chi^{(3)}$).
  • Hermann Haus e Wei-Ping Huang (1985–1991): Desenvolveram a CMT Vetorial Não-Ortogonal para plataformas de silício (SOI), demonstrando que em guias de alto contraste de índice ($\Delta n \approx 2{,}0$), a sobreposição intensa das caudas evanescentes viola a ortogonalidade estrita dos modos individuais.

2. Fundamentos Físicos e Dedução via Reciprocidade de Lorentz

Considere dois guias dielétricos paralelos (Guia 1 e Guia 2) alinhados ao longo do eixo $z$. O campo eletromagnético transversal total da estrutura é expandido como superposição dos modos não perturbados de cada guia:

$$\mathbf{E}(x,y,z,t) = \left[ A_1(z)\,\mathbf{\mathcal{E}}_1(x,y)\,e^{-j\beta_1 z} + A_2(z)\,\mathbf{\mathcal{E}}_2(x,y)\,e^{-j\beta_2 z} \right] e^{j\omega t}$$

A aplicação do Teorema da Reciprocidade de Lorentz entre o campo total perturbado e o modo não perturbado sobre a seção transversal infinita $S_\infty$ estabelece:

$$\frac{d}{dz} \iint_{S_\infty} \left( \mathbf{E} \times \mathbf{H}_1^* + \mathbf{E}_1^* \times \mathbf{H} \right) \cdot \hat{\mathbf{z}}\,dx\,dy = -j\omega \iint_{S_\infty} \Delta\epsilon_1(x,y)\, \mathbf{E} \cdot \mathbf{\mathcal{E}}_1^*\,dx\,dy$$

Sob a hipótese de variação lenta de amplitude (SVEA: $|d^2 A_m / dz^2| \ll \beta_m |dA_m / dz|$), obtemos as equações canônicas de acoplamento co-direcional:

$$\frac{dA_1}{dz} = -j\kappa_{12}\, A_2(z)\, e^{j\Delta\beta z}, \qquad \frac{dA_2}{dz} = -j\kappa_{21}\, A_1(z)\, e^{-j\Delta\beta z}$$

onde $\Delta\beta = \beta_1 - \beta_2$ é o descasamento de fase entre os modos. Pelo princípio da conservação de energia em guias sem perdas, $\kappa_{21} = \kappa_{12}^* \equiv \kappa$.

3. O Coeficiente de Acoplamento $\kappa$ e o Efeito do Gap

O coeficiente de acoplamento mútuo $\kappa$ quantifica a sobreposição ponderada da cauda evanescente do campo do Guia 2 com a perturbação de índice do núcleo do Guia 1:

$$\kappa = \frac{\omega \epsilon_0}{4 P_0} \iint_{\text{núcleo 1}} \left( n_{\text{core}}^2 - n_{\text{clad}}^2 \right) \mathbf{\mathcal{E}}_1^*(x,y) \cdot \mathbf{\mathcal{E}}_2(x,y)\,dx\,dy$$

Como o campo no cladding decai com taxa evanescente $\gamma = k_0 \sqrt{n_{\text{eff}}^2 - n_{\text{clad}}^2}$, o coeficiente $\kappa$ decai exponencialmente com a largura do gap $g$:

$$\kappa(g) \approx \kappa_0\, \exp(-\gamma\, g)$$
Sensibilidade Nanométrica de Fabricação em SOI:
Em silício ($n_{\text{eff}} \approx 2{,}45$, $n_{\text{clad}} = 1{,}444$ em $1550\text{ nm}$), a taxa de decaimento evanescente é $\gamma \approx 8{,}03\,\mu\text{m}^{-1}$, o que implica em um comprimento característico de decaimento de apenas $1/\gamma \approx 124\text{ nm}$. Uma variação litográfica de apenas $\pm 10\text{ nm}$ no processo DUV acarreta uma variação de mais de $8\%$ no valor de $\kappa$, alterando substancialmente a razão de divisão de potência.

4. Soluções Analíticas: Sincronismo vs. Descasamento de Fase

Supondo injeção unitária de luz exclusivamente no Guia 1 [$A_1(0) = 1, A_2(0) = 0$]:

4.1. Caso Síncrono ($\Delta\beta = 0$ — Guias Idênticos)

As potências ópticas guiadas nas portas Bar (direta) e Cross (cruzada) são dadas por:

$$P_{\text{bar}}(z) = |A_1(z)|^2 = \cos^2(\kappa z), \qquad P_{\text{cross}}(z) = |A_2(z)|^2 = \sin^2(\kappa z)$$
  • Comprimento de Acoplamento Total ($100\%$ Cross): $L_c = \frac{\pi}{2\kappa}$.
  • Divisor $3\text{ dB}$ ($50:50$ Cross/Bar): $L_{3\text{dB}} = \frac{L_c}{2} = \frac{\pi}{4\kappa}$.
  • Defasagem Intrínseca de $90^\circ$ ($-\pi/2$): A onda que acopla para a porta cross adquire defasagem de quadratura exata em relação à porta direta.

4.2. Caso Assíncrono ($\Delta\beta \neq 0$ — Guias Descasados)

Quando as larguras ou os índices são desiguais ($\Delta\beta = \beta_1 - \beta_2 \neq 0$), definindo o parâmetro de batimento $q = \sqrt{\kappa^2 + (\Delta\beta/2)^2}$:

$$P_{\text{cross}}(z) = \frac{\kappa^2}{q^2} \sin^2(q z) = \frac{\kappa^2}{\kappa^2 + (\Delta\beta/2)^2} \sin^2\left( \sqrt{\kappa^2 + (\Delta\beta/2)^2}\,z \right)$$
Filtragem e Multiplexação WDM:
A transferência máxima de potência é limitada por $F_{\text{max}} = \frac{\kappa^2}{\kappa^2 + (\Delta\beta/2)^2} < 1$. Se $|\Delta\beta| \gg \kappa$, a troca de potência é virtualmente nula. Esse efeito permite projetar multiplexadores WDM assimétricos onde apenas o comprimento de onda com $\Delta\beta(\lambda_0) = 0$ é transferido.

5. Acoplamento Contra-Direcional e Grades de Bragg

No acoplamento contra-direcional, a onda progressiva $A^+(z)$ transfere potência para a onda refletida regressiva $A^-(z)$ mediante uma perturbação periódica de passo $\Lambda$:

$$\frac{dA^+}{dz} = -j\kappa\, A^-(z)\, e^{j\Delta\beta z}, \qquad \frac{dA^-}{dz} = +j\kappa^*\, A^+(z)\, e^{-j\Delta\beta z}$$

onde o descasamento de Bragg é $\Delta\beta = 2\beta - \frac{2\pi m}{\Lambda}$. A ressonância de 1ª ordem ($m=1$) impõe a célebre condição de Bragg:

$$\lambda_{\text{Bragg}} = 2\, n_{\text{eff}}\, \Lambda$$

As soluções hiperbólicas fornecem a refletividade de pico na ressonância ($\Delta\beta = 0$):

$$R = \left| \frac{A^-(0)}{A^+(0)} \right|^2 = \tanh^2(|\kappa|\, L)$$

com largura espectral da banda proibida (photonic stopband):

$$\Delta\lambda_{\text{gap}} = \frac{\lambda_{\text{Bragg}}^2}{\pi n_g} \sqrt{\kappa^2 + \left(\frac{\pi}{L}\right)^2}$$

6. Formalismo Matricial de Espalhamento [S] e Transferência [T]

Para simulação circuital fotônica (EPDA), o acoplador direcional de 4 portas é modelado por sua matriz de espalhamento unitária:

$$\begin{bmatrix} b_1 \ b_2 \ b_3 \ b_4 \end{bmatrix} = \begin{bmatrix} 0 & 0 & t & -j\kappa_{\text{m}} \ 0 & 0 & -j\kappa_{\text{m}} & t \ t & -j\kappa_{\text{m}} & 0 & 0 \ -j\kappa_{\text{m}} & t & 0 & 0 \end{bmatrix} \begin{bmatrix} a_1 \ a_2 \ a_3 \ a_4 \end{bmatrix}$$

onde $t = \cos(\kappa L)$ e $\kappa_{\text{m}} = \sin(\kappa L)$, satisfazendo $|t|^2 + |\kappa_{\text{m}}|^2 = 1$ para componentes sem perdas ópticas.

7. Aplicações Tecnológicas da Teoria de Modos Acoplados em PICs

1. Divisores 3 dB e Roteadores

Comprimento ajustado para $L = L_c/2$, fornecendo divisão equitativa de potência ($50:50$) com quadratura intrínseca de fase de $90^\circ$.

2. Chaves Eletro-Ópticas $\Delta\beta$

Indução eletro-óptica de $\Delta n$ resulta em $\Delta\beta L = \sqrt{3}\pi$, comutando a luz da porta Cross para a porta Bar em nanossegundos.

3. Microrressonadores em Anel

Acoplamento evanescente no ponto tangente guia-anel. O acoplamento crítico ocorre quando a perda interna por volta do anel $a$ iguala $r$ ($r=a$).

4. Grades de Bragg & DFB

Reflexão de banda estreita para cavidades ressonantes, filtros rejeita-faixa e estabilização de comprimento de onda em lasers semicondutores.

5. Conversores de Modo & Polarização

Acoplamento intermodal com perturbações periódicas assimétricas para conversão $\mathrm{TE}_0 \to \mathrm{TE}_1$ (MDM) e rotação de polarização $\mathrm{TE} \leftrightarrow \mathrm{TM}$.

6. Óptica Não Linear ($\chi^{(2)}, \chi^{(3)}$)

Casamento de fase ($\Delta\beta = 0$) para geração de segundo harmônico (SHG), mistura de quatro ondas (FWM) e geração de pentes de frequência de Kerr.

8. Tabela Síntese das Configurações da CMT

Regime de Acoplamento Sentido Soluções Matemáticas Dispositivos Típicos Métrica de Mérito
Síncrono Simétrico ($\Delta\beta = 0$) Co-direcional ($+z, +z$) $\sin(\kappa z), \cos(\kappa z)$ Acopladores 3 dB, Divisores 50:50 Comprimento $L_c = \pi/(2\kappa)$
Assíncrono Descasado ($\Delta\beta \neq 0$) Co-direcional ($+z, +z$) Batimento ponderado $q$ Filtros WDM, Chaves $\Delta\beta$ Fração máx. $F = \frac{\kappa^2}{\kappa^2 + (\Delta\beta/2)^2}$
Periódico Contra-Propagante Contra-direcional ($+z, -z$) Hiperbólicas $\sinh, \cosh$ Grades de Bragg, Lasers DFB Refletividade $R = \tanh^2(|\kappa| L)$
Cavidade / Anel Óptico Barramento $\leftrightarrow$ Anel Auto-consistência recirculante Microrressonadores add-drop Fator $Q$, Acoplamento Crítico ($r=a$)
Intermodal Assimétrico Co-direcional intermodal Matriz fora da diagonal Conversores $\mathrm{TE}_0 \to \mathrm{TE}_1$, TE-TM Eficiência de conversão $> 98\%$
Paramétrico Não Linear Co / Contra-direcional Amplitudes acopladas não lineares SHG, FWM, Pentes de Kerr Casamento de fase $\Delta\beta = 0$

9. Laboratório Virtual Interativo: Simulação da Dinâmica CMT

Explore interativamente a troca contínua de potência luminosa entre dois guias de silício acoplados em função do espaçamento (gap $g$), do comprimento de interação $L$ e do descasamento de fase ($\Delta\beta$):

Coef. Acoplamento $\kappa$
0.086 µm⁻¹
Comprimento $L_c$ (100% Cross)
18.2 µm
Potência Bar ($P_1$)
12.5 %
Potência Cross ($P_2$)
87.5 %
Entrada $z=0$ (Guia 1 excitado) Região Acoplada (Interação Evanescente) Saída $z=L$
← Módulo Anterior: [MOD-PH03] Métodos Numéricos Próximo Módulo: [MOD-PH05] Dispositivos Passivos →