1. O que é Rigorosamente um Modo Óptico?
Em uma estrutura dielétrica planar ou tridimensional, a radiação eletromagnética confinada não se propaga com distribuições arbitrárias. Ela organiza-se em autofunções espaciais estáveis denominadas modos ópticos guiados.
Um modo guiado é uma solução particular das equações de onda de Helmholtz que satisfaz as condições de contorno de continuidade eletromagnética nas interfaces do guia e propaga-se ao longo da direção longitudinal $z$ mantendo seu perfil transversal espacial invariante:
- $\mathbf{e}(x, y)$ e $\mathbf{h}(x, y)$ são os perfis de campo vetoriais transversais;
- $\beta$ é a constante de propagação longitudinal $[\text{rad/m}]$;
- $n_{\text{eff}} = \frac{\beta}{k_0}$ é o índice de refração efetivo associado ao modo, onde $k_0 = \frac{2\pi}{\lambda_0}$.
A Condição Fundamental de Confinamento Óptico
Para que uma solução permaneça confinada e ligada (bound mode) sem perda contínua de energia por radiação no revestimento (cladding), seu índice efetivo deve necessariamente satisfazer: $$n_{\text{clad}} < n_{\text{eff}} < n_{\text{core}}$$ Se $n_{\text{eff}} \le n_{\text{clad}}$, o vetor de onda transversal no revestimento deixa de ser imaginário puro (evanescente) e torna-se real, caracterizando o corte do modo (cutoff) e o vazamento de potência óptica para o contínuo de radiação.
2. Relação com os Módulos de Eletromagnetismo
A transição do estudo de linhas de transmissão e guias metálicos para a fotônica integrada apoia-se em conceitos tratados nos módulos fundamentais:
Guias Metálicos ([MOD-EM06])
Confinamento delimitado por condutores elétricos perfeitos (PEC). A condição $E_\parallel = 0$ nas paredes metálicas anula estritamente o campo fora da cavidade. Todos os modos apresentam uma frequência de corte inferior estrita.
Guias Dielétricos (Fotônica)
Confinamento por reflexão interna total ([MOD-EM03]). O campo não é nulo fora do núcleo: existe uma cauda evanescente com decaimento exponencial no revestimento, crucial para o funcionamento de acopladores e sensores.
A Ausência de Corte do Modo Fundamental TE₀
Diferentemente de guias metálicos retangulares onde o modo $\text{TE}_{10}$ possui frequência de corte inferior finita ($f_c = \frac{c}{2a}$), em um guia dielétrico planar perfeitamente simétrico ($n_{\text{cover}} = n_{\text{substrate}} = n_{\text{clad}}$) o modo fundamental $\text{TE}_0$ propaga-se para qualquer espessura finita $d > 0$, possuindo frequência de corte nula ($V_c = 0$).
3. Modelo de Raios em Ziguezague e Deslocamento de Goos-Hänchen
A abordagem da óptica de raios representa a propagação guiada como um feixe de ondas planas que sofre sucessivas reflexões internas totais na interface núcleo-revestimento com ângulo $\theta > \theta_c = \arcsin(n_2/n_1)$.
Para que o campo interfira construtivamente após um ciclo transversal duplo (ida e volta em $x$), a defasagem acumulada deve ser múltiplo de $2\pi$:
O deslocamento de fase nas reflexões com polarização TE e TM é dado por:
Essa fase finita corresponde à penetração do campo óptico na camada de revestimento antes de ser refletido. No espaço real, isso se manifesta como um deslocamento espacial lateral ao longo da interface, conhecido como deslocamento de Goos-Hänchen ($\Delta z_{\text{GH}}$). Fisicamente, a onda comporta-se como se estivesse confinada em uma espessura aparente maior: $d_{\text{eff}} = d + 2/q$.
4. Dedução Vetorial das Equações de Dispersão (TE e TM)
Para uma estrutura planar 1D uniforme no plano $y$ ($\partial/\partial y = 0$), as componentes de campo desacoplam-se nas equações de Helmholtz:
Impondo a continuidade dos campos elétrico e magnético tangenciais nas interfaces $x = \pm d/2$, obtêm-se as relações transcendentais de dispersão:
Polarização TE ($E_y, H_x, H_z$)
Polarização TM ($H_y, E_x, E_z$)
O fator de escala $(n_1/n_2)^2$ nas equações de TM decorre da descontinuidade de $E_x$ nas bordas dielétricas. Em guias de silício sobre isolante ($\text{SOI}$), onde $n_1 \approx 3{,}48$ e $n_2 \approx 1{,}44$, tem-se $(n_1/n_2)^2 \approx 5{,}81$. Essa disparidade produz uma acentuada birrefringência de forma entre modos TE e TM na mesma estrutura física.
5. Parâmetros Universais: Frequência Normalizada ($V$) e Confinamento ($\Gamma$)
Para comparar guias de onda de diferentes tecnologias independentemente do comprimento de onda absoluto, utilizam-se as variáveis normalizadas:
O número total de modos suportados pelo guia simétrico é dado por $M = 1 + \lfloor V/\pi \rfloor$. Portanto, a condição de operação monomodo para guias planares simétricos é:
O fator de confinamento óptico ($\Gamma$) expressa a fração da potência total que viaja confinada no interior do núcleo:
6. Plataformas Materiais e Dimensões Físicas Típicas
A tabela a seguir apresenta as dimensões físicas típicas recomendadas para operação monomodo em comprimento de onda de telecomunicações ($\lambda_0 = 1550\text{ nm}$):
| Plataforma Material | Núcleo / Revestimento | Espessura ($d$) | Largura ($w$) | Índice Efetivo ($n_{\text{eff}}$) |
|---|---|---|---|---|
| Silício sobre Isolante (SOI) | $\text{Si}$ / $\text{SiO}_2$ | $220\text{ nm}$ | $450\text{ nm}$ | $\approx 2{,}44$ (TE$_0$) / $1{,}68$ (TM$_0$) |
| Silício de Camada Larga | $\text{Si}$ / $\text{SiO}_2$ | $3{,}0\ \mu\text{m}$ | $3{,}0\ \mu\text{m}$ | $\approx 3{,}45$ (Multimodal) |
| Nitreto de Silício ($\text{Si}_3\text{N}_4$) | $\text{Si}_3\text{N}_4$ / $\text{SiO}_2$ | $800\text{ nm}$ | $1000\text{ nm}$ | $\approx 1{,}75$ (TE$_0$) |
| Niobato de Lítio (TFLN) | $\text{LiNbO}_3$ / $\text{SiO}_2$ | $400\text{ nm}$ | $1000\text{ nm}$ | $\approx 1{,}85$ (TE$_0$) |
| Fosfeto de Índio (InP) | $\text{InGaAsP}$ / $\text{InP}$ | $1{,}8\ \mu\text{m}$ | $1{,}8\ \mu\text{m}$ | $\approx 3{,}24$ (TE$_0$) |
| Polímeros / Sílica | Polímero / $\text{SiO}_2$ | $3{,}0\ \mu\text{m}$ | $3{,}0\ \mu\text{m}$ | $\approx 1{,}48$ (TE$_0$) |
7. Laboratório Virtual 1D: Solucionador de Modos em Guia Planar (Slab)
Utilize o simulador analítico abaixo para calcular as raízes da equação transcendental de dispersão e renderizar o perfil transversal de campo $E(x)$ ou $H(x)$ e a densidade de potência $S_z(x)$ para diferentes materiais e dimensões:
8. Laboratório Virtual 2D: Visualizador Interativo da Distribuição Espacial de Modos ($x, y$)
Explore a distribuição contínua bidimensional de intensidade óptica $|E(x,y)|^2$ e campo transversal em guias dielétricos de seção retangular (Strip Waveguide). O simulador emprega o Método do Índice Efetivo (EIM) em tempo real. Mova o cursor sobre o mapa de calor 2D para inspecionar os perfis de corte 1D horizontal e vertical nas posições de mira.
9. Estudo Numérico Canônico: Guia SOI 500 nm × 220 nm (Ensaios FEM & FDTD)
Enquanto os métodos analíticos como o EIM decompõem o problema em fatias planares independentes, o projeto de circuitos integrados em fundição requer a solução vetorial exata das equações de Maxwell. A seguir, apresentamos os ensaios canônicos obtidos no Estudo 02: Análise Modal e Propagação FDTD em Guia de Silício, utilizando resolvedores modais por Elementos Finitos (Femwell / FEM) e propagação temporal por Diferenças Finitas (PyMeep / FDTD).
Os ensaios computacionais individuais (Estudos 01 a 10) desenvolvidos para a caracterização física de dispositivos fotônicos integrados estão integralmente liberados para consulta didática e análise de dados neste curso, fornecendo dados experimentais sintéticos de alta fidelidade para o aprendizado dos estudantes.
9.1. Distribuição de Campos Eletromagnéticos Transversais 2D
A simulação vetorial completa em guia de silício em fita ($500\text{ nm} \times 220\text{ nm}$) em $\lambda_0 = 1550\text{ nm}$ sobre óxido enterrado de $2\ \mu\text{m}$ confirma que a estrutura opera em regime estritamente monomodo para cada polarização:
Figura 9.1: Distribuição de intensidade $|E|$ e componentes transversais dominantes $E_x$ e $E_y$ para os modos fundamentais Quasi-TE₀ e Quasi-TM₀ calculados por FEM vetorial.
| Parâmetro (@ 1550 nm) | Modo Fundamental Quasi-TE₀ | Modo Fundamental Quasi-TM₀ | Impacto de Projeto |
|---|---|---|---|
| Índice Efetivo ($n_{\text{eff}}$) | 2,4149 | 2,0858 | Forte birrefringência de forma ($\Delta n_{\text{eff}} = 0,3291$) |
| Índice de Grupo ($n_g$) | 4,3353 | 4,0868 | Velocidade de propagação de pulsos ($v_g = c / n_g$) |
| Dispersão Cromática ($D$) | -0,0 ps/(nm·km) | -10.340,5 ps/(nm·km) | Modo TE opera no ponto de dispersão nula (ZDW) |
| Fator de Confinamento ($\Gamma$) | 76,4% no núcleo de Si | 44,5% no núcleo de Si | Modo TM possui cauda evanescente 2,3× mais intensa no cladding |
| Pureza de Polarização | 96,7% (Fração TE) | 97,1% (Fração TM) | Campos híbridos quase puros sem acoplamento espúrio |
9.2. Simulação Temporal FDTD: Propagação Estável e Ausência de Perdas
A injeção da autofunção modal calculada na malha temporal FDTD comprova a estabilidade física da solução guiada ao longo de $10{,}7\ \mu\text{m}$ de silício:
Figura 9.2: Mapa de campo $E_z$ em propagação FDTD contínua e perfil transversal de intensidade $|E|^2$ coletado no monitor de saída.
A transmitância medida no monitor de saída foi de 100,00%, confirmando atenuação nula por espalhamento numérico ou vazamento em substrato de óxido com fronteiras absorventes perfeitamente casadas (PML).
9.3. Curvas Espectrais de Dispersão e Índice de Grupo na Banda C
A varredura espectral vetorial entre $1530\text{ nm}$ e $1560\text{ nm}$ (resolução de $0{,}5\text{ nm}$, 61 comprimentos de onda) evidencia o papel dominante da dispersão de guia geométrica sobre a dispersão material do silício:
Figura 9.3: Curvas espectrais de índice efetivo $n_{\text{eff}}(\lambda)$, índice de grupo $n_g(\lambda)$ e parâmetro de dispersão cromática $D(\lambda)$ para as polarizações fundamentais.
Para o modo Quasi-TE₀, a inclinação de $n_{\text{eff}}(\lambda)$ gera uma curvatura que anula a dispersão de segunda ordem em $1550\text{ nm}$ ($D \approx 0$). Em contraste, o modo Quasi-TM₀ exibe forte dispersão normal ($D \approx -10.340\text{ ps/(nm}\cdot\text{km)}$), o que impõe restrições severas ao tráfego de dados de alta taxa sem compensação de dispersão.
10. Referências Bibliográficas
- Marcuse, D. (1974). Theory of Dielectric Optical Waveguides. Academic Press, New York (2ª ed., 1991).
- Yariv, A., & Yeh, P. (2007). Photonics: Optical Electronics in Modern Communications (6th ed.). Oxford University Press.
- Chrostowski, L., & Hochberg, M. (2015). Silicon Photonics Design: From Devices to Systems. Cambridge University Press.
- Osgood Jr., R. M., & Meng, X. (2021). Principles of Photonic Integrated Circuits. Springer Nature.
- Rumpf, R. C. (2022). Electromagnetic and Photonic Simulation for the Beginner: Finite-Difference Frequency-Domain in MATLAB. Artech House.
- Vivien, L., & Pavesi, L. (Eds.) (2013). Handbook of Silicon Photonics. CRC Press, Boca Raton.