Comment on "Angle insensitive filters based on Fabry-Perot resonance structures" [J. Appl. Phys. 136, 193102 (2024)]

Resumo Técnico do Comentário sobre "Filtros insensíveis ao ângulo baseados em estruturas de ressonância Fabry–Pérot"

Declaração do Problema
Este Comentário aborda uma discrepância significativa em relação aos resultados numéricos apresentados em um artigo recente de Cao et al. (J. Appl. Phys. 136, 193102, 2024), que alegou demonstrar filtros de banda estreita insensíveis ao ângulo utilizando uma estrutura de metassuperfície em cascata. O manuscrito original relatou que um pico de transmissão com uma largura de linha de ~5 nm permanecia em um comprimento de onda central constante (aproximadamente 1150 nm) para ângulos de incidência variando de 0° a 70°, com apenas uma leve diminuição na transmissão de pico. O autor deste Comentário, Michele Cotrufo, afirma a incapacidade de replicar esses resultados apesar de utilizar os parâmetros exatos fornecidos pelos autores originais, levantando preocupações sobre a validade do comportamento independente do ângulo relatado.

Metodologia
O Comentário emprega uma abordagem de simulação comparativa rigorosa para investigar as alegações do artigo original:

1. Tentativas de Replicação: Utilizando os parâmetros geométricos divulgados por Cao et al. ($d_g=240$ nm, $d_w=240$ nm, $d_m=310$ nm, $f=0.6$, $\Lambda=640$ nm), o autor realizou simulações usando três solvers eletromagnéticos distintos: Lumerical FDTD (correspondendo ao software original), Comsol (Método dos Elementos Finitos) e Reticolo (RCWA).
2. Análise de Dados: O autor analisou os dados brutos fornecidos por Cao et al. para examinar a sobreposição espectral e a dependência da intensidade em relação ao ângulo.
3. Varreduras de Parâmetros e Condições de Contorno: Para descartar artefatos numéricos, o autor conduziu campanhas sistemáticas variando a resolução da malha, harmônicos de Fourier e parâmetros geométricos.
4. Teste de Hipótese de Condição de Contorno (CC): Especificamente abordando a sugestão dos autores originais de que a discrepância poderia decorrer de "condições de contorno", o autor simulou intencionalmente a estrutura usando CCs incorretas (definindo periodicidade em vez das condições Bloch/Floquet necessárias para injeção de onda inclinada) para testar se esse erro poderia reproduzir os resultados originais.

Principais Contribuições e Resultados

• Falha em Replicar a Independência do Ângulo: Simulações usando condições de contorno físicas corretas (Bloch/Floquet) em três pacotes de software diferentes mostraram consistentemente que o pico de transmissão se desloca substancialmente à medida que o ângulo de incidência excede 5°. Aos 20°, o pico havia se deslocado para fora da faixa de 1150 nm, e picos estreitos adicionais apareceram. Este comportamento é consistente com a física padrão de grades de guias de onda ressonantes e contradiz a resposta plana e independente do ângulo relatada por Cao et al.
• Análise dos Dados Originais: Os dados brutos fornecidos por Cao et al. mostraram uma sobreposição perfeita "extremamente incomum" das formas de linha de transmissão em diferentes ângulos, incluindo oscilações e ruídos fracos idênticos. Além disso, a dependência da intensidade de pico em relação ao ângulo nos dados originais ajustou-se a uma função cosseno reescalada com um $R^2$ de 0,99996, uma perfeição matemática não típica de estruturas ressonantes físicas.
• Identificação de uma Fonte Específica de Erro: O Comentário demonstra que, quando o autor utilizou intencionalmente condições de contorno incorretas (periodicidade em vez de Bloch/Floquet) no Lumerical, os resultados da simulação corresponderam estreitamente aos dados apresentados por Cao et al. Esta configuração incorreta efetivamente injeta a onda em incidência normal, independentemente da configuração do ângulo da porta, reduzindo artificialmente a amplitude (explicando a queda de transmissão tipo cosseno) enquanto mantém a posição espectral.
• Inconsistência nos Parâmetros Relatados: O Comentário destaca uma contradição entre os parâmetros fornecidos para os espectros de transmissão (Fig. 7) e os perfis de campo (Fig. 8) no artigo original. A espessura da camada de sílica divulgada ($d_m=310$ nm) é mais fina que a espessura combinada da grade, contudo, os diagramas de campo na Fig. 8 retratam uma camada de sílica significativamente mais espessa que as metassuperfícies, sugerindo que as figuras podem não corresponder à mesma estrutura.
• Falta de Explicação Teórica: O artigo original não oferece nenhum mecanismo físico causal para a independência do ângulo alegada, afirmando meramente que os padrões de ondas estacionárias permanecem inalterados, sem explicar por que tal fenômeno deveria ocorrer nesta geometria específica.

Significado e Conclusão
O artigo conclui que o comportamento independente do ângulo alegado por Cao et al. é provavelmente não físico e resulta de uma grande discrepância entre os parâmetros de simulação relatados e os reais, ou, mais especificamente, do uso de condições de contorno incorretas nas simulações numéricas. O autor observa que os autores originais recusaram-se a compartilhar os arquivos de simulação, citando preocupações com propriedade intelectual, apesar de já terem publicado os parâmetros e resultados.

O significado deste Comentário reside em sua demonstração de que o desempenho "notável" do dispositivo relatado é um artefato de erro de simulação e não um fenômeno físico novo. Ao reproduzir os resultados errôneos através da aplicação deliberada de condições de contorno incorretas, o autor fornece uma explicação plausível para a discrepância e sublinha a necessidade de protocolos de simulação rigorosos e transparência na pesquisa de fotônica computacional. O Comentário foi publicado no arXiv após a retratação do artigo original e a falta de uma explicação clara por parte do conselho editorial sobre o papel deste Comentário no processo de retratação.