Adjoint-Based Gradient Evaluation for Metasurface Inverse Design via Affine Geometric Transformations

1. Problema e Contexto

O projeto de metasuperfícies (estruturas compostas por um grande número de espalhadores interagentes, como meta-átomos) para controle preciso de campos eletromagnéticos enfrenta desafios significativos devido à sua natureza multiescala. Embora métodos de otimização "inversa" sejam necessários para superar as limitações das abordagens tradicionais de "célula unitária" (que assumem variações lentas e têm baixa aceitação angular), a otimização direta de grandes metasuperfícies é computacionalmente proibitiva.

• Desafio Principal: O cálculo do gradiente de uma função de custo em relação a milhares de parâmetros de design (um para cada meta-átomo) usando métodos de diferenças finitas exigiria a solução de um número imenso de problemas diretos, tornando o processo inviável.
• Limitações de Métodos Existentes: A otimização de forma livre (freeform), embora poderosa, exige processos de binarização complexos e repetidos para viabilidade de fabricação, aumentando a complexidade computacional ao escalar para grandes áreas.

2. Metodologia Proposta

O artigo propõe uma metodologia sistemática que combina a Método da Variável Adjoint (Adjoint Variable Method - AVM) com uma parametrização geométrica baseada em transformações afins.

• Parametrização Geométrica: Em vez de otimizar a forma livre de cada meta-átomo (o que geraria um número enorme de variáveis), o método define o design através de um conjunto limitado de parâmetros geométricos para cada meta-átomo:

  • Posição (translação rígida).
  • Orientação (rotação).
  • Tamanho/Forma (expansão anisotrópica).
  • Isso reduz drasticamente a dimensionalidade do problema de otimização.

• Transformações Afins: O método permite que cada meta-átomo sofra uma transformação afim completa a cada iteração de otimização, combinando rotação, escala anisotrópica e translação.

• Cálculo do Gradiente via Método Adjoint:

  • O núcleo da contribuição é a derivação de uma expressão analítica para o gradiente da função de custo em relação a esses parâmetros geométricos.
  • Utilizando o Teorema da Reciprocidade de Lorentz, o problema é reformulado para calcular o gradiente resolvendo apenas um problema direto adicional (o problema adjunto), independentemente do número de meta-átomos ou parâmetros.
  • O gradiente é expresso como uma integral sobre a fronteira do meta-átomo, envolvendo o campo direto ($E_f$), o campo adjunto ($E_a$) e um termo de ponderação geométrica ($w$) que depende da transformação específica (rotação, expansão ou translação).
  • Desacoplamento: O método desacopla o cálculo em duas partes: a definição da fonte adjunta (dependente apenas da função de custo) e a função de ponderação geométrica (dependente apenas da transformação). Isso permite reutilizar a mesma solução adjunta para diferentes tipos de transformações geométricas.

3. Principais Contribuições

1. Formulação Geral para Transformações Afins: Diferente de trabalhos anteriores que focavam apenas em expansão isotrópica, este trabalho fornece expressões analíticas para o gradiente sob rotação, expansão anisotrópica e translação rígida.
2. Eficiência Computacional: Demonstra que o gradiente para todos os parâmetros de design de uma grande metasuperfície pode ser computado resolvendo apenas dois problemas diretos (um direto padrão + um adjunto), independentemente do número de meta-átomos.
3. Viabilidade de Fabricação: Ao restringir a otimização a parâmetros geométricos bem definidos, o método naturalmente evita designs irracionais ou de difícil fabricação, eliminando a necessidade de etapas de binarização pós-otimização típicas da otimização de forma livre.
4. Versatilidade de Funções de Custo: O artigo deriva fontes adjuntas específicas para diversas funções de custo comuns, incluindo:
   • Magnitude do produto escalar.
   • Norma da diferença entre campos.
   • Ângulo entre campos.
   • Diferença de quadrados das normas.
   • Intensidade em pontos específicos.

4. Resultados e Validação Numérica

Os autores validaram a metodologia através de simulações numéricas utilizando o Método dos Elementos de Fronteira (BEM) para problemas de propagação TEz com luz vermelha ($\lambda_0 = 660$ nm) e meta-átomos de dióxido de titânio ($TiO_2$).

• Precisão do Gradiente: Comparou-se o gradiente calculado pelo método adjunto com o obtido por diferenças finitas. Os resultados mostraram uma concordância quase perfeita para diferentes parâmetros (rotação, largura e posição), validando a precisão matemática da derivação.
• Exemplos de Design Inverso:
  • Focalização de Luz: Otimização de uma metasuperfície com 128 meta-átomos para focalizar luz em um ponto específico ($37\lambda_0$ de distância). O design resultante apresentou um padrão de rotação altamente não simétrico e variável, alcançando um foco com largura de meio pico (FWHM) de $0.84\lambda_0$.
  • Focalização Múltipla: O método foi aplicado com sucesso para focalizar luz simultaneamente em dois e três pontos distintos, demonstrando a capacidade de lidar com tarefas de design complexas e de múltiplos objetivos.

5. Significado e Conclusão

Este trabalho representa um avanço significativo na otimização inversa de grandes metasuperfícies. Ao integrar o método adjunto com uma parametrização geométrica baseada em transformações afins, os autores conseguiram:

• Reduzir drasticamente o custo computacional, tornando viável o design de estruturas de grande escala.
• Manter a flexibilidade para explorar um amplo espaço de designs (rotação, deformação, translação) sem sacrificar a eficiência.
• Fornecer um framework robusto que garante que os designs otimizados sejam geometricamente realizáveis, superando as limitações de escalabilidade dos métodos de forma livre.

Em suma, a abordagem proposta oferece uma ferramenta poderosa e eficiente para o desenvolvimento de metasuperfícies e metalentes de alto desempenho para aplicações como controle de feixe, imageamento e computação analógica.