Simulação de Filtros Ópticos Utilizando o Método da Matriz de Transferência
Neste tutorial, utilizamos o solucionador de matriz de transferência ótica da OghmaNano para simular a reflexão, a transmissão e a interferência ótica em filtros óticos multicamadas de película fina. Essas estruturas são amplamente utilizadas em refletores de Bragg, espelhos dielétricos, revestimentos antirreflexo, OLEDs e células solares de película fina, onde a resposta ótica é controlada pelo índice de refração e pela espessura de cada camada na estrutura.
1. Contexto
O método da matriz de transferência (MMT) é amplamente utilizado para simular estruturas ópticas multicamadas finas, como microcavidades OLED, células solares de filme fino, espelhos dielétricos, refletores de Bragg distribuídos (DBR), revestimentos antirreflexo e filtros ópticos seletivos em comprimento de onda. Nestas estruturas, em incidência normal, a luz é parcialmente refletida e transmitida em cada interface entre materiais com diferentes índices de refração, produzindo interferência óptica entre ondas eletromagnéticas que se propagam para frente e para trás. O método da matriz de transferência modela esta interferência de forma rigorosa, permitindo que a reflexão, a transmissão e a absorção óptica de toda a pilha multicamadas sejam calculadas com precisão. A resposta óptica de estruturas multicamadas finas é governada pela interferência de filme fino entre ondas eletromagnéticas refletidas e transmitidas no interior da pilha. Uma representação esquemática da estrutura multicamadas e dos campos ópticos que se propagam é mostrada na figura ??.
No âmbito do formalismo da matriz de transferência, o campo elétrico em cada camada é expresso como a soma de uma onda que se propaga para a frente e de uma onda que se propaga para trás.
\[ E_i(z)=E_i^{+}e^{-\zeta_i z}+E_i^{-}e^{\zeta_i z} \]onde \(E_i^{+}\) e \(E_i^{-}\) representam, respetivamente, as amplitudes do campo que se propaga para a frente e para trás, e
\[ \zeta_i=\frac{2\pi}{\lambda}\bar{n}_i \]onde \(\bar{n}=n+j\kappa\) representa o índice de refração complexo do material. Ao aplicar as equações de Maxwell e as condições de contorno eletromagnéticas em cada interface, obtêm-se os coeficientes de reflexão e transmissão de Fresnel.
\[ r_{ij}=\frac{n_i-n_j}{n_i+n_j}, \qquad t_{ij}=\frac{2n_i}{n_i+n_j} \]permitindo que os campos óticos de ambos os lados de uma interface sejam relacionados por meio de uma matriz de transferência.
\[ \begin{pmatrix} E_i^{+} \\ E_i^{-} \end{pmatrix} = \frac{1}{t_{ij}} \begin{pmatrix} 1 & r_{ij} \\ r_{ij} & 1 \end{pmatrix} \begin{pmatrix} E_j^{+} \\ E_j^{-} \end{pmatrix} \]enquanto a propagação através de uma camada de espessura \(d_i\) pode ser expressa da seguinte forma:
\[ \begin{pmatrix} E_i^{+}(z+d_i) \\ E_i^{-}(z+d_i) \end{pmatrix} = \begin{pmatrix} e^{\zeta_i d_i} & 0 \\ 0 & e^{-\zeta_i d_i} \end{pmatrix} \begin{pmatrix} E_i^{+}(z) \\ E_i^{-}(z) \end{pmatrix} \]Dessa forma, a estrutura ótica multicamadas completa é obtida multiplicando as matrizes de interface e de propagação de cada camada.
\[ M=M_1M_2M_3 \cdots M_n \]a partir do qual a potência ótica refletida, transmitida e absorvida pode ser calculada diretamente. Esta abordagem é computacionalmente eficiente, numericamente estável e requer apenas operações simples de multiplicação de matrizes, o que a torna adequada para simular filtros óticos multicamadas e estruturas de interferência de filmes finos.
Ao ajustar as espessuras e os índices de refração das camadas, é possível conceber filtros óticos com propriedades espectrais personalizadas. Por exemplo, uma única camada dielétrica de um quarto de onda pode reduzir a reflexão da superfície num comprimento de onda específico, enquanto camadas alternadas de alto e baixo índice de refração formam um refletor de Bragg com uma forte banda de bloqueio ótico.
2. Primeiros passos
Para iniciar o seu primeiro cálculo de filtro ótico, abra a janela «Nova simulação» na faixa de opções «Ficheiro» do menu principal. Clique duas vezes no exemplo de «Filtro ótico» (veja ??) e guarde a simulação numa pasta no seu disco. Em seguida, verá a janela principal (veja ??) com uma pilha de múltiplas camadas, composta por cerca de dez camadas alternadas. Clique no botão «Executar simulação» (reproduzir) para calcular os espectros; quando terminar, os resultados aparecerão no separador «Resultados» (veja ??).
3. Análise dos resultados.
Após executar a simulação, clique duas vezes em «Saída Óptica» em ??. Isso abrirá o Editor de Simulação Óptica. O editor contém várias abas. A primeira aba, «Distribuição de Fótons», é exibida automaticamente (veja ??). Aqui, é mostrada a densidade de fótons dentro da cavidade, e a estrutura em camadas do filtro é claramente visível como faixas verticais. A segunda aba, «Distribuição de Fótons Absorvida» (veja ??), mostra onde os fótons são absorvidos. Neste exemplo, a absorção é fraca, mas não nula, porque um dos materiais foi configurado com um pequeno coeficiente de absorção \(\alpha\). Finalmente, a aba «Luz Refletida» (veja ??) exibe o espectro de refletância. O resultado demonstra uma forte reflexão entre aproximadamente 500 nm e 800 nm, enquanto a luz fora desta faixa é transmitida em graus variáveis, característica de um filtro do tipo Bragg.
O espectro da luz transmitida confirma o comportamento de rejeição de banda do filtro. A luz com comprimentos de onda entre aproximadamente 300 nm e 500 nm é transmitida de forma eficaz, enquanto os comprimentos de onda na faixa de 500 a 800 nm são fortemente bloqueados. Em comprimentos de onda mais longos, acima de 800 nm, observa-se uma nova transmissão, o que demonstra a natureza multibanda da resposta do filtro.
4. Editar o filtro
Para inspecionar ou modificar a pilha, abra a guia «Estrutura do dispositivo» na janela principal e clique em «Editor de camadas» (consulte ??). O editor lista cada camada do dispositivo, juntamente com a sua espessura, material ótico e configurações. Pode editar as espessuras das camadas diretamente na tabela, alterar os materiais, adicionar ou remover camadas e reorganizá-las conforme necessário para o seu projeto de filtro.
5. Medição da espessura do filtro ótico por digitalização.
Em seguida, analisaremos a possibilidade de realizar digitalizações automáticas das larguras dos filtros. Isso permite que você crie filtros, realize digitalizações automáticas e estude como a estrutura do filtro afeta tanto a transmissão quanto a reflexão. Ao realizar digitalizações sistemáticas, você pode gerar rapidamente uma visão geral de vários projetos de filtros diferentes e escolher aquele que melhor se adapta à sua aplicação.
Acesse a faixa de opções «Automação» na janela principal e clique em «Análise de parâmetros» (consulte ??). A janela «Análise de parâmetros» será aberta e deverá já conter uma nova entrada de análise chamada «novo». Clique duas vezes nesta entrada e a janela de configuração da análise será exibida (consulte ??). Por fim, clique em «Iniciar análise» para começar o cálculo. A análise pode demorar um pouco para ser concluída, dependendo do número de valores de parâmetros que selecionou.
Ao clicar em «Executar análise», executou o pequeno programa apresentado em ??. Explicaremos este programa mais tarde, mas, em essência, ele varia a espessura da camada com o índice de refração mais elevado dentro de um conjunto de valores. Abra o separador «Resultados» na janela «Análise de parâmetros» (veja ??): verá quatro diretórios, cada um dos quais corresponde a uma das espessuras analisadas. Cada diretório contém uma simulação completa com os ficheiros habituais, diferenciando-se apenas na espessura dessa camada.
Na pasta principal da análise, também encontrará ícones especiais de «curvas múltiplas» (ficheiros CSV com um símbolo de várias linhas). Estes agrupam as curvas correspondentes de todas as subsimulações. Clique duas vezes em «optical_output» para o abrir (veja ??); verá reflect.csv e transmit.csv. Ao abrir reflect.csv, será apresentado um gráfico com a refletância de todas as espessuras analisadas, conforme mostrado em ??.
reflect.csv e transmit.csv.
6. Compreensão do programa de análise da janela de parâmetros.
Se observar novamente ??, verá que estão listadas cinco linhas. Cada linha especifica um parâmetro a ser analisado. Neste exemplo, o parâmetro é a espessura (dy) de camadas específicas na estrutura de epitaxia. As entradas são epitaxy.layer0.dy, epitaxy.layer2.dy, epitaxy.layer4.dy, epitaxy.layer6.dy e epitaxy.layer8.dy. Estas linhas correspondem às camadas com alto índice de refração na estrutura do dispositivo.
Na primeira linha, pode ver uma lista de valores de espessura: 1.0e-7, 2e-7, 4e-7, 8e-7. Estes são os valores que o programa irá atribuir a layer0. A operação para esta linha está definida como «varredura», o que significa que o programa irá variar sistematicamente a espessura de layer0 ao longo destes valores especificados.
Para as outras camadas (layer2, layer4, layer6 e layer8), a operação está associada a epitaxy.layer0.dy. Na coluna de valores, isto aparece como duplicado. Isto indica ao programa para copiar o valor atualmente definido para layer0 e aplicá-lo às outras camadas de alto índice. Na prática, sempre que layer0 for alterado, as outras camadas de alto índice são atualizadas automaticamente para a mesma espessura.
Em resumo, este programa analisa um conjunto de valores de espessura para uma camada de alto índice de refração e, em seguida, replica essa espessura em todas as outras camadas de alto índice de refração na estrutura. O solucionador é executado para cada caso, permitindo que você explore como o desempenho do filtro varia em função da espessura da camada.
7. Resumo
Nesta parte do tutorial, aprendeu como automatizar o processo de conceção de filtros óticos através da utilização de varreduras de parâmetros. Variou sistematicamente as espessuras das camadas, duplicou parâmetros em várias camadas e explorou como estas alterações afetam a transmissão e a reflexão. Com estas ferramentas, pode agora avaliar rapidamente vários projetos de filtro possíveis e identificar as estruturas que melhor se adequam à sua aplicação.