Tutorial de Simulação de Éxcitons

1. Introdução

Neste tutorial, investigaremos como simular a dinâmica de éxcitons em dispositivos fotovoltaicos orgânicos (OPV). Quando um fóton é absorvido em um semicondutor orgânico, ele gera um par elétron–lacuna ligado conhecido como éxciton. Em dispositivos OPV, esses éxcitons devem derivar ou se difundir até uma interface doador–aceitador onde possam se dissociar em portadores de carga livres. Esse processo pode ser modelado usando uma equação de dissociação de éxcitons, que descreve tanto o transporte de éxcitons quanto sua conversão final em elétrons e lacunas. No OghmaNano, essa dinâmica de éxcitons é introduzida sobre as equações padrão de drift–diffusion para capturar o quadro fotofísico completo. Uma discussão mais detalhada sobre quando a modelagem de éxcitons é necessária (e quando pode ser omitida) é fornecida em Éxcitons e recombinação geminada.

2. Primeiros passos

Na aba New simulation da faixa file, clique para abrir a janela New simulation (veja Figure 1a). Nessa janela, você encontrará muitas categorias de exemplos pré-configurados. Dê um duplo clique na pasta Exciton simulations para abrir o submenu de exemplos relacionados a éxcitons disponíveis (veja Figure 1b). Para este tutorial, usaremos o modelo Exciton device, que fornece um ponto de partida simples para simular geração, transporte e dissociação de éxcitons em estruturas OPV.

3. Execute a simulação e inspecione as saídas de éxcitons

Depois que o dispositivo for aberto, você poderá ver a pilha na janela principal. Antes de executar a simulação, vá para a aba Electrical e verifique se o botão Exciton solver está pressionado — geralmente ele já está por padrão, mas vale a pena confirmar (??). Isso garante que a dinâmica de éxcitons esteja ativada no modelo. Depois disso, volte ao menu File e pressione o botão Run simulation para executar o modelo (??).

Quando a simulação terminar, você poderá inspecionar os resultados na aba Output (??). Ela lista todos os arquivos gravados em disco, incluindo saídas relacionadas a éxcitons. Ao dar um duplo clique em jv.csv, você poderá plotar a curva JV (??). O gráfico se parece com uma curva JV padrão que você obteria de qualquer simulação de dispositivo, mas neste caso o Exciton solver estava ativado, o que significa que a dinâmica de éxcitons está incluída na física subjacente.

Se você der um duplo clique na pasta exciton_output dentro do diretório de saída, poderá acessar os resultados detalhados do solver de éxcitons (veja ??). Esse diretório contém todas as saídas do solver, incluindo constantes em função da posição e quantidades calculadas derivadas da dinâmica de éxcitons. Por exemplo, dar um duplo clique em exciton.csv produz um gráfico da distribuição de éxcitons ao longo da espessura do dispositivo (veja ??). De forma semelhante, Gn.csv fornece a taxa de geração de elétrons em função da posição, enquanto Gp.csv mostra a taxa de geração de lacunas (veja ??). No caso de um dispositivo 1D simples, esses dois arquivos de taxa de geração são efetivamente idênticos.

4. Equação de transporte de éxcitons e parâmetros

A distribuição de éxcitons dentro do dispositivo é governada pela equação de transporte de éxcitons:

\[ \frac{\partial X}{\partial t} = \nabla \!\cdot \!\big(D\,\nabla X\big) + G_{\mathrm{optical}} - k_{\mathrm{dis}}\,X - k_{\mathrm{FRET}}\,X - k_{\mathrm{PL}}\,X - \alpha\,X^{2} \]

Aqui \(X(\mathbf{r},t)\) é a densidade de éxcitons (m\(^{-3}\)); \(D\) é o coeficiente de difusão de éxcitons (m\(^2\)s\(^{-1}\)); \(G_{\mathrm{optical}}\) é a taxa local de geração de éxcitons fornecida pelo modelo óptico (proporcional aos fótons absorvidos); \(k_{\mathrm{dis}}\) é a taxa de dissociação em cargas livres; \(k_{\mathrm{FRET}}\) é a taxa de transferência de energia por ressonância de Förster; \(k_{\mathrm{PL}}\) é a taxa de decaimento radiativo; e \(\alpha\) é o coeficiente de aniquilação éxciton–éxciton (m\(^3\)s\(^{-1}\)). Quando esse modelo está ativado, as equações eletrônicas de drift–diffusion usam um termo de geração definido como \(G = k_{\mathrm{dis}}\,X\), opcionalmente restrito a regiões interfaciais se configurado.

Todos esses parâmetros podem ser visualizados e modificados no Electrical parameter editor, acessível a partir do menu principal (veja ??). Os campos específicos de éxcitons são agrupados sob o título Excitons. Eles estão resumidos na tabela abaixo.

5. Como o solver se encaixa no processo de simulação

Quando o solver de éxcitons está desligado (parte superior de ??), a óptica transfer-matrix calcula a taxa de geração de portadores e a passa diretamente ao solver drift–diffusion. Quando o solver de éxcitons está ligado (parte inferior de ??), a óptica fornece em vez disso uma taxa de geração de éxcitons ao solver de éxcitons. O solver de éxcitons então evolui essa população — difundindo, transferindo (FRET), dissociando, decaindo radiativamente e aniquilando — usando os parâmetros definidos no Electrical parameter editor, e produz a taxa final de geração de portadores para o solver drift–diffusion. Em resumo, o solver de éxcitons se insere entre a óptica e o transporte elétrico para que você possa modelar a física de éxcitons sem alterar a configuração óptica nem as equações drift–diffusion.

A sequência de operações na simulação com éxcitons ativados é ilustrada em ??. Primeiro, o solver óptico é executado, calculando o perfil de absorção de fótons fatia por fatia ao longo do dispositivo. Em seguida, o solver de éxcitons é executado, onde os éxcitons gerados são propagados, transferidos ou dissociados até que o solver convirja, tipicamente em algumas dezenas de passos. Finalmente, a taxa de geração de éxcitons é examinada em função tanto da profundidade do dispositivo quanto do comprimento de onda, como calculado pelo modelo de matriz de transferência. Juntas, essas saídas formam o pipeline que conecta o modelo óptico ao solver de éxcitons e, por fim, às equações drift–diffusion.