Simulação de célula solar orgânica de grande área (PM6:Y6): de redes de contato 3D a dispositivos iluminados
1. Introdução
Neste tutorial, continuamos a série sobre simulação de dispositivos de grande área, com foco em dispositivos que já ultrapassaram pequenas células de laboratório e passaram para geometrias em escala de aplicação. Ele se situa entre dois fluxos de trabalho já existentes no manual: projetar contatos de grande área (sem dispositivos ativos) e módulos perovskita 3D complexos. Este exemplo faz a ponte entre os dois, introduzindo diodos e absorção de luz em uma estrutura de contato previamente projetada, mas sem chegar a um módulo interconectado completo.
Aqui simulamos uma única célula solar orgânica de grande área baseada em uma arquitetura do tipo PM6:Y6 com contatos metálicos hexagonais impressos. O objetivo é entender como um dispositivo que apresenta bom desempenho em pequena área se comporta quando ampliado, onde o fluxo lateral de corrente, a queda de tensão e a resistência de contato passam a dominar. Embora este tutorial use materiais OPV, o fluxo de trabalho é geral: ao substituir os dados ópticos n/k e ajustar os valores dos componentes, o mesmo modelo pode representar o seu dispositivo de pesquisa.
Para aplicar este fluxo de trabalho ao seu próprio dispositivo, basta fornecer dois conjuntos de entradas: resistividades das camadas (tipicamente obtidas da literatura ou de estruturas de teste simples) e um pequeno número de parâmetros de diodo, como o fator de idealidade, a fotocorrente e a corrente de saturação. Depois que esses valores são definidos, o modelo permite escalar virtualmente um dispositivo de laboratório e identificar quais mecanismos de perda limitarão sua eficiência ao ser transferido para geometrias reais de grande área.
Aqui não usamos o solver completo de drift-diffusion 3D, embora o OghmaNano o ofereça. Para dispositivos de grande área, o drift-diffusion 3D totalmente resolvido rapidamente se torna impraticável do ponto de vista computacional e muitas vezes é desnecessário quando as limitações dominantes surgem do transporte lateral e da resistência de contato, em vez do movimento de portadores através de cada camada. Em vez disso, o OghmaNano constrói um modelo de circuito 3D sensível à geometria, no qual o transporte lateral e os contatos são representados por elementos resistivos e o comportamento fotovoltaico é capturado usando diodos embutidos.
Enquanto o problema elétrico é resolvido usando uma formulação de circuito, a óptica é tratada de forma rigorosa. A fotogeração é calculada usando um solver óptico por transfer-matrix avaliado em fatias unidimensionais ao longo da geometria tridimensional, mantendo efeitos de interferência e absorção em filmes finos, enquanto o solver de circuito trata a coleta de corrente em grande área.
2. Criar a simulação
Na janela principal, clique em New simulation. Você verá uma lista de categorias de simulação (veja ??). Em seguida, clique duas vezes em Large area PM6:Y6 solar cell (veja ??).
Após criar a simulação, a janela principal é aberta e exibe a geometria do dispositivo em três dimensões (veja ??). O contato superior é uma malha metálica hexagonal (favo de mel) rotulada como Ag. Essa camada de prata atua como uma grade de coleta de corrente altamente condutiva. Como ela é implementada como malha, a maior parte da superfície permanece aberta, permitindo que a luz atravesse o contato. Abaixo da malha de prata está o PEDOT:PSS, que serve como camada de espalhamento lateral de corrente. Como no tutorial de contatos de grande área, a carga é primeiro coletada lateralmente na camada polimérica e depois transferida para a grade de prata de baixa resistência para extração eficiente.
A camada fotoativa é PM6:Y6. Eletricamente, essa camada é representada por um elemento diodo com geração óptica habilitada, enquanto componentes resistivos adicionais garantem que o solver leve em conta sua condutividade finita. Na base da pilha está um contato contínuo de Ag que cobre toda a área do dispositivo e completa a conexão elétrica.
A barra verde visível na janela principal pode parecer estar flutuando, mas é apenas um indicador visual da condição de contorno de extração. Ela marca a região onde o solver aplica o contato elétrico usado para extrair corrente da rede de circuito.
3. Inspecionar a pilha epitaxial de camadas e contatos
Esta simulação (e o exemplo anterior de contato de grande área) é definida usando uma estrutura epitaxial em camadas. O dispositivo é construído como uma pilha de camadas de filmes finos usando o editor de camadas, o que garante que cada camada fique corretamente sobre a seguinte, sem lacunas e com tratamento consistente da espessura. Essa abordagem é adequada para dispositivos planares de grande área, onde a pilha vertical é fixa e o comportamento lateral é introduzido por meio do modelo de circuito. Ela difere de simulações posteriores no estilo de módulo (veja o tutorial de módulo interconectado), onde dedos, barramentos e layouts mais complexos são definidos explicitamente usando geometria 3D completa.
Para inspecionar a estrutura, abra o Layer editor para visualizar a pilha de camadas (veja ??), e em seguida abra o painel Contacts no editor de estrutura do dispositivo (veja ??). Geometricamente, a configuração é simples, mas importante: o contato inferior cobre toda a área do dispositivo, enquanto o contato superior é muito menor e define a região onde a corrente é extraída, consistente com o marcador verde de extração mostrado na visualização 3D. Entender essa geometria de contato é essencial, pois ela define os caminhos laterais de corrente que dominam o comportamento de dispositivos de grande área.
4. Inspecionando a malha
Mude para a aba Circuit diagram e clique no ícone refresh no canto inferior esquerdo (o botão no estilo reciclagem). Essa ação constrói a representação em circuito do dispositivo diretamente a partir da geometria definida e da pilha de camadas. A malha de circuito resultante é mostrada em ??. Cada ligação na malha corresponde a um elemento de circuito, como um resistor ou diodo, com a cor indicando a camada ou região de material da qual ele se origina.
Ao girar a visualização, você pode identificar claramente as regiões de contato superior e inferior. Em ?? e ??, os pontos azuis indicam os nós onde a corrente é extraída da malha de circuito. Esses nós definem as condições de contorno dos contatos e representam os locais nos quais o solver de circuito coleta corrente do dispositivo.
5. Parâmetros elétricos e geração óptica
Abra Electrical parameters no painel de estrutura do dispositivo. Esse editor controla as propriedades de resistor e diodo usadas pelo modelo de circuito de grande área. Exemplos de configurações são mostrados em ?? e ??. Para a camada PM6:Y6 você verá parâmetros específicos de um elemento diodo, incluindo o fator de idealidade n e a corrente de polarização reversa (saturação) I0. Se você substituir esses valores por valores das suas próprias curvas experimentais, verá como o seu dispositivo será ampliado. Você também pode habilitar a geração óptica de portadores de carga. Quando a geração óptica está habilitada, o solver óptico por transfer-matrix fornece um termo de fotogeração que entra no modelo elétrico como uma fotocorrente Iph na equação do diodo iluminado:
I(V) = I0 [ exp( qV / (n k T) ) − 1 ] − Iph
Neste tutorial, o seletor controla efetivamente o termo Iph. Com a geração óptica desabilitada, o circuito se comporta como uma rede de diodos no escuro. Com ela habilitada, a iluminação injeta corrente no circuito por meio do elemento diodo, com Iph determinado pelo fluxo de fótons absorvidos a partir da óptica por transfer-matrix.
💡 Detalhe de implementação: a camada ativa pode ser discretizada em múltiplos elementos resistivos ao longo de sua espessura para capturar transporte em volume e resistência em série. Entretanto, o modelo de circuito sempre contém uma única camada de diodo, normalmente colocada na base da região ativa. Isso reflete o fato de que o diodo representa a resposta JV global e a fotogeração do dispositivo, enquanto os elementos resistivos contabilizam a queda de tensão e o espalhamento de corrente dentro do material.
6. Executando a simulação
Para iniciar a simulação, clique em Run simulation (o triângulo azul na barra de ferramentas principal), ou pressione F9. O solver começará imediatamente a execução, e o terminal de saída será atualizado em tempo real (veja ??).
As linhas coloridas que aparecem no início da saída correspondem ao cálculo óptico. Nesse estágio, o solver por transfer-matrix está avaliando a fotogeração ao longo do dispositivo, fatia por fatia, para construir o perfil espacialmente resolvido de geração. Quando isso é concluído, o solver passa ao problema elétrico e começa a resolver a rede de circuito.
Vale a pena prestar atenção à saída do terminal enquanto a simulação é executada. Essa informação não é ruído diagnóstico: ela fornece uma visão útil sobre se a simulação está se comportando fisicamente. No instantâneo mostrado, podemos ver que aproximadamente −9.48 × 102 A m−2 está fluindo para fora de um contato, enquanto −9.48 × 101 A m−2 está fluindo para fora do outro. Ao mesmo tempo, o residual reportado pelo solver F = 1.7 × 10−14 indica que a continuidade de corrente foi satisfeita com precisão extremamente alta, com o primeiro ponto JV resolvido em cerca de 17 ms.
A diferença de densidade de corrente entre os dois contatos é totalmente esperada. O contato inferior cobre toda a área do dispositivo, enquanto o contato superior ocupa apenas uma pequena região. Como resultado, a mesma corrente total é extraída através de uma área muito menor no contato superior, levando a uma densidade de corrente significativamente maior ali. Ler a saída do terminal dessa maneira fornece uma verificação rápida e eficaz de que a geometria e as condições de contorno da simulação são fisicamente consistentes. Se o erro residual aumentasse para ordem de unidade, ou se o tempo de execução se tornasse incomumente longo, isso geralmente indicaria um problema, como uma malha quebrada ou um contato desconectado. Quando a simulação terminar, abra a aba Output (veja ??). As saídas disponíveis são descritas em mais detalhes no tutorial anterior e exploradas mais a fundo no próximo tutorial. Aqui, focamos apenas nas curvas JV resolvidas por contato.
Cada contato elétrico tem seu próprio arquivo JV, jv_contact0.csv e jv_contact1.csv, que registram a corrente fluindo para fora daquele contato específico. Embora as densidades de corrente sejam diferentes, a corrente total extraída é a mesma quando as áreas de contato são levadas em conta. Abrir jv_contact0.csv produz a curva JV mostrada abaixo.
7. Editando parâmetros do material
As propriedades do material podem ser modificadas diretamente por meio do editor de objeto. Ao clicar com o botão direito em qualquer camada na visualização 3D e selecionar Edit object, o editor de objeto é aberto (veja ?? e ??).
A partir daqui, você pode alterar o material óptico atribuído ao objeto, modificar sua forma geométrica (por exemplo, ao usar contatos padronizados como malhas tipo favo de mel) ou ajustar atributos visuais, como cor. Nesta simulação em estilo epitaxial, a posição absoluta do objeto é em grande parte definida pela pilha de camadas, de modo que parâmetros como translação e rotação normalmente têm pouco efeito. O que mais importa são a atribuição do material óptico e a geometria do objeto.
💡 Experimente: clique com o botão direito na camada PM6:Y6 e substitua seu material óptico por uma perovskita. Em seguida, execute novamente a simulação e observe como a curva JV muda. Se desejar, você também pode atualizar os parâmetros do diodo — como o fator de idealidade, a corrente de saturação e a fotocorrente — usando valores reportados na literatura de perovskitas.