Tutorial de Domínio de Éxcitons (Parte A): Executar uma Simulação 3D de Domínio de Éxcitons

1. Introdução

Compreender como os éxcitons são gerados, transportados e dissociados é central para o funcionamento de células solares orgânicas de heterojunção em volume (BHJ). Nesses sistemas, a luz é absorvida predominantemente na fase doadora, gerando pares elétron–lacuna ligados conhecidos como éxcitons. Esses éxcitons devem difundir através do material doador e alcançar uma interface doador–aceitador antes de decair. Na interface eles podem se dissociar em elétrons e lacunas livres, com uma probabilidade determinada pela morfologia local, dimensionalidade e pelas constantes cinéticas relevantes.

O equilíbrio entre difusão de éxcitons, dissociação interfacial e processos de perda concorrentes desempenha portanto um papel decisivo na determinação da eficiência efetiva de fotogeração. Entretanto, a modelagem explícita de éxcitons foi historicamente evitada não por dificuldade conceitual, mas por incerteza nos parâmetros. Comprimentos de difusão de éxcitons, tempos de vida, taxas de dissociação interfacial, canais de decaimento radiativos e não radiativos e processos de aniquilação eram frequentemente pouco conhecidos para um determinado sistema material. Em muitos fluxos práticos de modelagem esses efeitos eram portanto absorvidos em um único fator escalar fator de eficiência de fóton, \(\eta_{\mathrm{photon}}\), representando perdas geminadas líquidas sem resolver explicitamente o transporte e a cinética subjacentes.

Essa situação está mudando. Avanços na caracterização experimental estão começando a fornecer medições diretas de tempos de vida de éxcitons, comprimentos de difusão e canais de perda em sistemas do tipo BHJ. À medida que esses parâmetros se tornam melhor determinados, a modelagem explícita de domínios de éxcitons torna-se cada vez mais informativa — não como meio de prever eficiências absolutas de dispositivos, mas como forma de explorar como tamanho de domínio, dimensionalidade e cinética de materiais determinam conjuntamente um rendimento efetivo de geração de carga. O modelo de domínio de éxcitons utilizado neste tutorial foi desenvolvido e publicado como parte desse trabalho recente de experimentação e modelagem (ver, por exemplo, Nature Materials 21, 55–61 (2022)).

Neste tutorial adotamos portanto uma geometria de célula unitária idealizada, mas totalmente tridimensional: um domínio doador embutido em uma matriz aceitadora, representado inicialmente como uma esfera doadora dentro de uma caixa aceitadora. Embora a geometria seja simplificada, o modelo trata explicitamente a geração de éxcitons, difusão, dissociação interfacial e processos de perda concorrentes. Isso o torna um quadro prático e fisicamente transparente para explorar como parâmetros medidos experimentalmente interagem, para verificar sua consistência interna e para desenvolver intuição sobre como morfologia e cinética controlam conjuntamente a eficiência efetiva de fotogeração em três dimensões.

2. Equações governantes

Neste tutorial, o transporte de éxcitons é tratado explicitamente resolvendo a equação tridimensional de difusão de éxcitons em todo o domínio de simulação. Quando o modelo de éxcitons está habilitado, a absorção óptica alimenta diretamente a população de éxcitons, que então evolui de acordo com a seguinte equação:

\[ \frac{\partial X}{\partial t} = \nabla \cdot \left( D \nabla X \right) + 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 (\(\mathrm{m^{-3}}\)) e \(D\) é o coeficiente de difusão de éxcitons (\(\mathrm{m^2\,s^{-1}}\)). O termo fonte \(G_{\mathrm{optical}}\) (\(\mathrm{m^{-3}\,s^{-1}}\)) representa a geração de éxcitons devido à absorção óptica. Os termos restantes descrevem os processos concorrentes que removem éxcitons da população:

• Dissociação (\(k_{\mathrm{dis}} X\)): conversão de éxcitons em elétrons e lacunas livres nas interfaces doador–aceitador. Este é o canal produtivo que contribui para a geração de carga.
• Transferência do tipo Förster (\(k_{\mathrm{FRET}} X\)): transferência não radiativa de energia de éxcitons para sítios ou domínios vizinhos, redistribuindo energia sem gerar diretamente portadores livres.
• Perda radiativa (\(k_{\mathrm{PL}} X\)): recombinação de éxcitons via emissão de fótons (fotoluminescência), representando um canal intrínseco de perda.
• Aniquilação (\(\alpha X^2\)): aniquilação bimolecular éxciton–éxciton em altas densidades, importante sob excitação forte ou em geometrias confinadas.

O coeficiente de difusão \(D\) é parametrizado usando um comprimento de difusão de éxcitons \(L\) (m) e um tempo de vida de éxcitons \(\tau\) (s) por meio de

\[ D = \frac{L^2}{\tau}. \]

Neste tutorial executamos uma simulação apenas de éxcitons: o solver calcula a densidade tridimensional de éxcitons \(X(\mathbf{r})\) e os termos de reação associados, mas não resolve o transporte de elétrons e lacunas. O termo de dissociação \(k_{\mathrm{dis}} X\) é portanto interpretado como uma taxa espacialmente resolvida potencial de geração de carga, indicando onde e com que eficiência os éxcitons se converteriam em portadores livres nas interfaces doador–aceitador.

3. Criar a simulação de domínio de éxcitons

Abra a janela New simulation (Figura ??). Ela pode ser acessada clicando em New simulation na janela principal. Se você clicar duas vezes em Exciton simulations, verá a biblioteca de éxcitons (Figura ??). Clique duas vezes em Exciton domain para abrir o projeto de exemplo.

4. Inspecionar a geometria e os parâmetros

Após abrir o exemplo, a janela principal exibe uma cena tridimensional simples: uma esfera doadora embutida dentro de uma caixa aceitadora (Figura ??). Os raios verdes indicam iluminação incidente de cima. Para inspecionar os parâmetros de éxcitons, clique em Electrical parameters no painel à esquerda (na aba Device structure). Isso abre o editor de parâmetros elétricos. Se você rolar até o final dessa janela, encontrará uma seção chamada Excitons, onde os parâmetros específicos de éxcitons são listados separadamente para cada objeto na cena (a camada circundante e a esfera embutida), conforme mostrado nas Figuras ?? e ??.

A seção de éxcitons define os parâmetros que aparecem na equação de difusão–reação. Estes incluem o comprimento de espalhamento de éxcitons \(L\), que define o coeficiente de difusão via \(D = L^{2}/\tau\), e o tempo de vida de éxcitons \(\tau\), que controla a escala temporal global para transporte na ausência de canais adicionais de perda. A recombinação radiativa é descrita pela taxa de fotoluminescência \(k_{\mathrm{PL}}\), enquanto processos de transferência de energia não radiativos são representados pela taxa do tipo Förster \(k_{\mathrm{FRET}}\). Em densidades mais altas de éxcitons, a aniquilação bimolecular é incluída através do coeficiente \(\alpha\), correspondente ao termo \(-\alpha X^{2}\) na equação governante. Finalmente, a taxa de dissociação \(k_{\mathrm{dis}}\) especifica quão eficientemente os éxcitons se convertem em estados de transferência de carga nas interfaces doador–aceitador e é o principal parâmetro que controla a eficiência efetiva de geração de carga neste modelo.

A principal característica conceitual deste exemplo é que a esfera (doador) e a camada (matriz aceitadora) utilizam física de éxcitons diferente. Nas capturas de tela acima, a esfera doadora inclui um \(k_{\mathrm{PL}}\) não nulo (um canal de perda). A região circundante é configurada com um \(k_{\mathrm{dis}}\) muito grande, de modo que quando os éxcitons alcançam (ou entram) nessa região eles se dissociam rapidamente. Isso faz com que a fronteira doador/aceitador se comporte como um sumidouro eficiente, produzindo o perfil característico “alto no centro, depletado perto da fronteira”.

5. Executar a simulação

Clique no triângulo azul Run simulation para iniciar o solver. A saída do terminal aparecerá na janela principal, como mostrado na Figura ??. As linhas iniciais relatam informações gerais sobre a configuração da simulação e a geometria. O diagnóstico principal aparece quando o solver de éxcitons começa a executar e tem a forma Exciton f() = … , steps = …. Essa linha relata o progresso do solver tridimensional de éxcitons. A quantidade f() é um residual que mede o quão distante o campo atual de densidade de éxcitons está da auto-consistência; à medida que o solver itera, esse valor deve diminuir. O contador steps indica o número de iteração. No exemplo mostrado, o residual cai de um valor inicial para 2.97 × 104 no passo 2, após o qual o solver relata que a equação foi resolvida e termina.

Em um laptop moderno típico, este exemplo deve completar em cerca de 5–10 segundos. Se o tempo de execução se estender para minutos, isso geralmente indica que a malha está excessivamente fina ou que a geometria ou parâmetros foram modificados de forma que aumentou muito o número de pontos de malha. Nesse caso, vale a pena revisar as configurações de malha ou refazer as etapas de configuração anteriores antes de prosseguir.

6. Plotando a saída

Quando a execução terminar, mude para a aba Output (Figura ??). Os itens mais relevantes para este tutorial são o diretório exciton_output/, que contém os resultados espaciais de éxcitons, e o arquivo exciton_sim_info.json, que resume estatísticas globais de geração e perda e será utilizado posteriormente na Parte C. Clicar duas vezes em exciton_output/ revela o conteúdo mostrado na Figura ??.

Abra G.csv para visualizar a taxa de geração de éxcitons produzida pelo modelo óptico (Figura ??). Na janela de gráfico 3D você pode rotacionar a cena com o mouse e usar os controles deslizantes Z e Y na parte inferior para cortar o volume. Esta é frequentemente a forma mais rápida de verificar se a geometria, a direção de iluminação e a malha estão se comportando como esperado.

Em seguida, abra exciton.csv (Figura ??), que mostra a densidade de éxcitons em regime estacionário. Você deve observar a assinatura característica de difusão em direção a um sumidouro interfacial: alta densidade de éxcitons no interior da esfera doadora e depleção em direção à fronteira doador–aceitador onde a dissociação é forte. Use os controles deslizantes para confirmar que essa depleção segue a geometria da interface e não é um artefato do ângulo de visualização.

Finalmente, abra Gn.csv (Figura ??), que mostra a taxa de geração de elétrons espacialmente resolvida resultante da dissociação de éxcitons. Neste exemplo você deve observar um anel (ou casca) característico de geração localizado próximo à interface doador–aceitador. Use os controles Z e Y para mover o plano de corte através do volume e explorar como esse perfil de dissociação evolui com a posição.

7. Quais são os arquivos em exciton_output/?

O diretório exciton_output/ contém dois tipos distintos de arquivos. O primeiro grupo consiste em campos calculados produzidos pelo solver de éxcitons, como a densidade de éxcitons e taxas de geração derivadas da dissociação. O segundo grupo contém mapas de parâmetros: cópias tridimensionais dos parâmetros de éxcitons definidos na interface gráfica. Eles são escritos principalmente como auxílio de consistência e depuração, permitindo verificar se os parâmetros de material pretendidos foram aplicados corretamente no espaço. As Tabelas 1 e 2 resumem essas duas categorias. Todas as quantidades são reportadas em unidades SI. (Alguns arquivos podem não ser escritos se os submodelos correspondentes estiverem desativados.)

Os arquivos de mapa de parâmetros na Tabela 2 não são saídas do solver no sentido usual. Eles simplesmente refletem os parâmetros de éxcitons definidos na seção Electrical → Excitons da interface gráfica, escritos na mesma malha tridimensional que os campos de solução. Seu objetivo é fornecer uma maneira transparente de verificar se os parâmetros foram atribuídos corretamente aos objetos e regiões do domínio, particularmente ao trabalhar com geometrias complexas.

8.1 Interpretando exciton_sim_info.json (eficiência de dissociação)

Enquanto os gráficos tridimensionais são inestimáveis para visualizar onde os éxcitons são gerados, perdidos e dissociados, eles não fornecem por si só resumos numéricos concisos. O arquivo exciton_sim_info.json cumpre esse papel ao relatar taxas de geração e perda integradas espacialmente para cada objeto na simulação (a layer circundante, a sphere embutida e quaisquer outras formas adicionais que você possa adicionar).

Para visualizar essas informações, abra a aba principal Output e clique duas vezes em exciton_sim_info.json. O arquivo contém uma estrutura JSON aninhada que lista, para cada objeto, a taxa total de geração de éxcitons e as taxas totais associadas a cada canal de perda ou dissociação concorrente. Essas quantidades permitem quantificar como mudanças em parâmetros de material, geometria ou morfologia se traduzem em mudanças na eficiência de geração de carga.

A chave de nível superior segment0 simplesmente marca o início do relatório. Dentro dela, o campo name identifica o objeto envolvente (aqui, a layer), enquanto o bloco aninhado shape0 corresponde ao objeto embutido (a sphere). Cada bloco relata a contribuição total desse objeto para processos de geração, dissociação e perda.

A Tabela 2 resume os campos que você encontrará neste arquivo e como interpretá-los.

Em conjunto, os gráficos tridimensionais e os resumos numéricos em exciton_sim_info.json fornecem uma visão complementar da simulação: os gráficos mostram onde os processos ocorrem, enquanto o arquivo JSON quantifica quanto cada processo contribui. Isso o torna uma ferramenta útil para analisar sistematicamente como mudanças em parâmetros ou geometria afetam a dissociação de éxcitons e a geração efetiva de carga.

9. Integração dos resultados

Em conjunto, os resultados tridimensionais contam uma história física coerente. A absorção óptica gera éxcitons principalmente no interior do domínio doador. Esses éxcitons então difundem através do doador, são depletados perto da interface doador–aceitador e se dissociam eficientemente quando atingem a região aceitadora circundante. Os gráficos espaciais mostram onde esses processos ocorrem; o relatório resumido exciton_sim_info.json quantifica quanto cada processo contribui.

Neste exemplo, o relatório separa claramente a física subjacente por região espacial, permitindo que processos de geração, perda e dissociação sejam interpretados de maneira transparente.

• Esfera (doador): A maior parte da absorção óptica — e portanto da geração de éxcitons — ocorre dentro da esfera doadora, refletida por um valor grande de tot_G no bloco shape0/sphere. Dentro do doador, o principal canal de perda concorrente neste exemplo é o decaimento radiativo, quantificado por tot_Rk_pl. A razão entre tot_Rk_pl e tot_G fornece uma medida direta de quão fortemente as perdas radiativas suprimem a população de éxcitons antes que ela alcance a interface doador–aceitador.
• Camada (matriz aceitadora): A camada circundante absorve apenas fracamente neste exemplo, portanto sua contribuição para geração de éxcitons (tot_G) é pequena. Seu papel principal é a conversão de éxcitons: a camada é configurada com um forte canal de dissociação, tornando tot_Rk_dis o termo dominante no bloco layer. Isto corresponde diretamente à dissociação interfacial observada nos gráficos espaciais, onde Gn.csv mostra uma região em anel ou casca de geração de elétrons concentrada na fronteira doador–aceitador.

A métrica final dis_eff reporta a eficiência global de dissociação do sistema (aqui aproximadamente \(89\%\)). Fisicamente, isso indica que a maioria dos éxcitons gerados no doador difunde até a interface e se dissocia, com apenas uma pequena fração perdida por decaimento radiativo. Outros canais, como transferência de Förster ou aniquilação éxciton–éxciton, estão inativos neste exemplo específico.

Em outras palavras, a simulação opera em um regime em que a dissociação interfacial supera o decaimento de éxcitons, resultando em um alto rendimento efetivo de conversão de fótons em portadores livres. Essa combinação de visualização tridimensional e relatórios quantitativos fornece uma descrição compacta mas fisicamente completa de como geometria e parâmetros cinéticos determinam conjuntamente a eficiência de geração de carga em um modelo de domínio de éxcitons.