Full-Wave Harmonic Balance Framework for Dispersive Time-Varying Photonic Structures

Resumo Técnico: Estrutura de Equilíbrio Harmônico de Onda Completa para Estruturas Fotônicas Dispersivas e Variáveis no Tempo

Definição do Problema
Estruturas fotônicas variáveis no tempo oferecem capacidades únicas, como conversão de frequência, não reciprocidade e controle dinâmico de frente de onda, ao quebrar a simetria de translação temporal. No entanto, modelar precisamente esses sistemas permanece um desafio significativo quando os materiais envolvidos são fortemente dispersivos (por exemplo, óxidos condutores transparentes, grafeno, materiais com epsilon próximo de zero). Nesses sistemas, a modulação temporal não pode ser tratada como uma simples perturbação instantânea de um parâmetro independente da frequência. Em vez disso, a causalidade dita que a resposta do material depende do histórico de excitação, o que significa que cada harmônico de Floquet gerado experimenta uma resposta de material distinta, permanecendo simultaneamente acoplado a todos os outros harmônicos através da modulação. Os métodos computacionais existentes enfrentam compromissos (trade-offs): métodos no domínio do tempo (como FDTD) frequentemente exigem simulações longas para atingir o estado estacionário e têm dificuldade com atualizações causais e dispersivas, enquanto métodos analíticos ou de ondas acopladas especializados carecem de flexibilidade geométrica para lidar com estruturas 3D complexas, meios inhomogêneos ou espectros de modulação arbitrários.

Metodologia
Os autores propõem uma estrutura de equilíbrio harmônico de onda completa que resolve a resposta eletromagnética em regime estacionário diretamente no domínio da frequência, evitando a progressão temporal. O núcleo da metodologia envolve a reformulação do problema variante no tempo como um conjunto acoplado de subproblemas no domínio da frequência dentro de um espaço de Floquet estendido.

Elementos principais da formulação incluem:

• Representação de Fonte Secundária: Componentes variantes no tempo são representados não como condições de contorno dependentes do tempo, mas como fontes secundárias induzidas. Densidades de polarização volumétrica ($P$) são usadas para meios volumétricos modulados, e densidades de corrente superficial ($J_s$) são usadas para folhas condutoras moduladas.
• Expansão de Floquet: Todos os campos e fontes induzidas são expandidos em uma série de Floquet truncada ($F(r, t) = \sum F_n(r)e^{j\omega_n t}$), onde $\omega_n = \omega_0 + n\Omega$.
• Formulação de Sistema Acoplado: O problema é resolvido como um sistema linear global $(L - SC)X = Y_{ext}$.
  • $L$ representa o operador de Maxwell no domínio da frequência para o background estático, avaliado na frequência física de cada harmônico ($\omega_n$) para considerar a dispersão do material.
  • $S$ representa os operadores de injeção de fonte.
  • $C$ representa os operadores de geração de fonte derivados da expansão de Fourier temporal dos parâmetros do material. Esses operadores introduzem acoplamento fora da diagonal no espaço harmônico, representando o mecanismo de conversão de frequência.
• Implementação: A estrutura é implementada utilizando o Método de Elementos Finitos (FEM). O vetor de incógnitas padrão do FEM é elevado para um vetor de domínio harmônico, permitindo que cada grau de liberdade espacial carregue múltiplos coeficientes harmônicos. Isso preserva a flexibilidade geométrica da modelagem de onda completa enquanto captura o acoplamento inter-harmônico.

Principais Contribuições

1. Estrutura Unificada: O método fornece uma única estrutura algébrica para tratar tanto meios volumétricos modulados quanto folhas condutoras moduladas, suportando espectros de modulação temporal arbitrários (tom único, múltiplos tons ou não sinusoidais).
2. Tratamento de Dispersão: Diferente de métodos que assumem dispersão fraca ou usam aproximações de envelope, esta estrutura avalia a dispersão do material e os operadores de radiação na frequência específica de cada harmônico gerado.
3. Observáveis Resolvidos por Harmônico: A solução fornece diretamente amplitudes, fases e perfis de campo em regime estacionário para todos os harmônicos de Floquet retidos, permitindo o cálculo de espalhamento, absorção, fluxo de potência e distribuições de campo próximo específicos para cada harmônico sem necessidade de pós-processamento de dados do domínio do tempo.

Resultos
A estrutura foi validada e demonstrada através de dois exemplos representativos:

1. Cilindro de Grafeno Rolado com Bombeamento Paramétrico:

   • Configuração: Um cilindro dielétrico revestido com uma folha de grafeno com um nível de Fermi modulado no tempo ($\Omega = 2\omega_0$), operando em um regime de amplificação paramétrica degenerada.
   • Validação: Os resultados numéricos para espalhamento e absorção resolvidos por harmônico mostraram excelente concordância com uma solução analítica de espalhamento de Floquet derivada para a geometria cilíndrica.
   • Descobertas: A simulação capturou com sucesso a redistribuição de energia do portador para o harmônico descendente ($n=-1$), confirmando a capacidade da estrutura em modelar forte dispersão e modulação rápida em sistemas de plasmons de superfície ressonantes.

2. Metasuperfície de Espaço-Tempo Óptico (baseada em ITO):

   • Configuração: Uma metasuperfície reflexiva consistindo em uma grade de $TiO_2$ binária em uma camada de ITO (óxido de índio e estanho) com epsilon próximo de zero, modulada no tempo e apoiada por ouro. A frequência de plasma do ITO foi modulada para impulsionar mudanças ultrarrápidas na permissividade.
   • Aplicação: A estrutura foi utilizada para projetar e otimizar a metasuperfície para deflexão de feixe por conversão de frequência. O objetivo era converter uma onda incidente em bandas laterais e redirecioná-las para uma ordem de difração espacial específica ($m=1$).
   • Descobertas: O design otimizado redirecionou com sucesso as bandas laterais geradas ($\omega_{\pm 1}$) para a primeira ordem de difração espacial. Os resultados demonstraram a interação entre a dispersão ENZ, o acoplamento de Drude dependente do tempo e a difração espacial, validando a utilidade do método para o design de dispositivos ópticos complexos e realistas.

Significância
O artigo afirma que esta estrutura estabelece uma plataforma geral para modelar e projetar estruturas fotônicas ativas e dispersivas. Ao combinar a transparência física das abordagens de equilíbrio harmônico com a flexibilidade geométrica do FEM de onda completa, o método supera as limitações dos solvers existentes para sistemas variantes no tempo. Ele permite que pesquisadores quantifiquem diretamente como a modulação redistribui a energia eletromagnética entre canais deslocados em frequência, fornecendo observáveis essenciais (potência, absorção e difração resolvidas por harmônico) para o design de sistemas fotônicos de espaço-tempo de próxima geração baseados em grafeno, óxidos condutores transparentes e metasuperfícies ressonantes. O trabalho posiciona a estrutura não apenas como uma ferramenta de análise, mas também para o design inverso de dispositivos que requerem conversão de frequência, amplificação paramétrica e controle dinâmico de frente de onda.