Medições CELIV em células solares de Perovskita: Parte a
1. Visão geral
Neste tutorial começaremos introduzindo os conceitos básicos da técnica Charge Extraction by Linearly Increasing Voltage (CELIV), descrevendo como ela é usada para investigar mobilidade e densidade de portadores em semicondutores de filme fino. Após estabelecer a base teórica, passaremos rapidamente aos passos práticos, mostrando como configurar e executar uma simulação CELIV no OghmaNano. Essa combinação de teoria e simulação prática fornecerá tanto uma compreensão clara do método quanto as ferramentas para aplicá-lo aos seus próprios dispositivos.
2. Contexto de CELIV
A técnica Charge Extraction by Linearly Increasing Voltage (CELIV) é um método experimental usado para estudar a mobilidade e a densidade de portadores de carga em células solares orgânicas e de perovskita. Em uma configuração típica, uma célula solar é conectada a um gerador de sinais que aplica uma polarização reversa linearmente crescente, enquanto o transiente de corrente resultante é registrado em um osciloscópio. Iluminação contínua opcional ou um curto pulso de laser podem ser adicionados para gerar portadores de carga. Um esquema dessa configuração é mostrado em ??.
Quando a rampa de tensão é aplicada, a corrente medida consiste em duas contribuições distintas: uma linha de base retangular da capacitância geométrica do dispositivo e um pico adicional da extração de portadores de carga móveis no semicondutor. O termo de capacitância geométrica está sempre presente, enquanto o pico relacionado aos portadores fornece informação direta sobre mobilidade e densidade de portadores. Um exemplo desses transientes de tensão e corrente é mostrado em ??.
A mobilidade dos portadores pode ser determinada diretamente a partir do transiente CELIV medindo o tempo em que a corrente de extração atinge seu máximo. No caso mais simples de distribuição uniforme de portadores e recombinação desprezível, a mobilidade é dada pela equação CELIV: \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] onde \(d\) é a espessura da camada ativa, \(A\) é a taxa da rampa de tensão (\(A = \mathrm{d}V/\mathrm{d}t\)), e \(t_{\text{max}}\) é o tempo em que ocorre o pico de corrente. Essa relação conecta o transiente observado experimentalmente ao parâmetro fundamental de transporte de interesse.
3. Limitações de CELIV
A derivação da equação de mobilidade CELIV baseia-se em várias hipóteses simplificadoras. Primeiro, assume-se que apenas um tipo de portador de carga domina o transporte, de modo que apenas uma única espécie de portador se move quando a rampa de tensão é aplicada. Segundo, assume-se que os portadores estão distribuídos uniformemente, deslizando suavemente para fora do dispositivo de forma semelhante ao caso idealizado mostrado à esquerda de ??. Terceiro, os processos de recombinação são considerados desprezíveis durante a extração, de modo que não ocorre perda significativa de portadores. Na realidade, essas condições raramente são satisfeitas. Como ilustrado no lado direito de ??, os portadores de carga frequentemente ocupam armadilhas rasas, médias ou profundas, e sua liberação retardada produz transientes de extração alargados ou distorcidos. Como em muitas técnicas experimentais, a mobilidade extraída usando CELIV deve, portanto, ser considerada uma aproximação e não o verdadeiro valor microscópico, e na verdade a mobilidade aparente pode evoluir durante o próprio transiente (doi:10.1063/1.4818267). Essa limitação é especialmente importante em semicondutores orgânicos desordenados, onde desordem energética e aprisionamento afetam fortemente o transporte. Em contraste, materiais de perovskita frequentemente exibem mobilidades mais altas e menores densidades de armadilhas, tornando a análise CELIV padrão mais robusta e mais fácil de interpretar.
💡 Teste seu entendimento: Tente aplicar a teoria CELIV a um caso de uso simples.
Suponha que você tenha um dispositivo com espessura d = 200 nm,
uma taxa de rampa de tensão de A = 2 × 106 V/s,
e um pico de corrente observado em tmax = 5 µs.
Usando a equação de mobilidade CELIV, estime a mobilidade dos portadores.
Mostrar resposta
A mobilidade CELIV é dada por
\[
\mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}}
\]
Substituindo os valores:
d = 200 × 10-9 m,
A = 2 × 106 V/s,
tmax = 5 × 10-6 s.
\[ \mu = \frac{2 (200 × 10^{-9})^{2}}{3 (2 × 10^{6}) (5 × 10^{-6})^{2}} \approx 1.1 × 10^{-4} \; \text{cm}^2/\text{Vs} \]
Este exemplo simples mostra como o tempo de pico em um transiente CELIV pode ser convertido em uma estimativa de mobilidade.
4: Executando uma simulação CELIV no OghmaNano
Clique em Nova simulação. Isso abre a biblioteca de tipos de dispositivos disponíveis, mostrada em ??. Clique duas vezes em Células de Perovskita (destacado em vermelho) para abrir a pasta de exemplos de perovskita. Você verá uma lista de simulações predefinidas, incluindo Dispositivo MAPbI₃, Célula solar de Perovskita, e demonstrações CELIV dedicadas, como mostrado em ??. Para este tutorial, selecione o Exemplo de célula solar de Perovskita – CELIV. Quando solicitado, salve a simulação em uma pasta na qual você tenha permissão de escrita.
Após selecionar o Exemplo de célula solar de Perovskita – CELIV, a janela principal de simulação é aberta (veja ??). Para iniciar o cálculo, clique em Executar simulação (ícone azul de play) ou pressione F9. O OghmaNano então resolverá as equações de drift–diffusion dependentes do tempo e gerará o transiente CELIV.
jv.csv (dados da curva JV), time_j.csv (corrente de extração vs. tempo),
time_v.csv (tensão aplicada vs. tempo), e optical_output (distribuições de campo).
Clique duas vezes em qualquer arquivo para visualizar os resultados correspondentes nas ferramentas de plotagem integradas.
Etapa 5: Visualizar os resultados CELIV
Abra a aba Saída
(??)
e clique duas vezes em time_v.csv. Isso plota o programa de tensão aplicado —
a rampa CELIV — mostrada em
??.
Em seguida, clique duas vezes em time_j.csv para exibir o transiente de corrente de extração,
como em ??.
Você notará que o transiente tem a forma familiar de CELIV, mas aparece **invertido**:
por convenção de sinal no OghmaNano, a corrente que sai do dispositivo é negativa. Quando a rampa em
?? é aplicada,
a carga é puxada para fora do dispositivo, então o dispositivo fornece corrente negativa (pico de extração).
O pico posterior de sinal oposto (positivo em
??) ocorre quando a polarização
retorna e a carga corre de volta para dentro do dispositivo depois que o transiente retorna ao seu estado inicial.
O gráfico em ?? mostra o transiente CELIV após multiplicar a corrente simulada por –1 para corresponder à apresentação convencional usada em experimentos. Nessa forma, a linha de base corresponde à capacitância geométrica do dispositivo, enquanto o pico positivo acentuado indica a extração de portadores de carga móveis durante a rampa de tensão aplicada. A magnitude e o tempo desse pico são as quantidades-chave usadas para determinar a mobilidade e a densidade de portadores. Apresentar o transiente nessa orientação convencional facilita comparar diretamente simulações do OghmaNano com medições CELIV publicadas. Para extrair a mobilidade, identifique o tempo em que a corrente atinge seu máximo (\(t_\text{max}\)) e substitua esse valor na equação de mobilidade CELIV \(\mu = \tfrac{2d^2}{3At_\text{max}^2}\), onde \(d\) é a espessura da camada ativa e \(A\) é a taxa da rampa de tensão aplicada.
div style="clear: both;">🧪 Exemplo resolvido: Leia o tempo de pico CELIV em
?? e estime a mobilidade.
Assuma uma espessura de camada ativa de d = 600 nm, uma taxa de rampa de tensão de
A = 5 × 105 V/s, e um tempo de pico medido de
tmax = 3 µs. Que mobilidade você obtém?
Mostrar resposta
Use a relação CELIV \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] com \(d = 600\times10^{-9}\,\mathrm{m}\), \(A = 5\times10^{5}\,\mathrm{V\,s^{-1}}\), \(t_{\text{max}} = 3\times10^{-6}\,\mathrm{s}\).
\[ \mu = \frac{2(600\times10^{-9})^{2}}{3\,(5\times10^{5})\,(3\times10^{-6})^{2}} \approx 5.3\times10^{-8}\ \mathrm{m^{2}\,V^{-1}\,s^{-1}} \] Convertendo para \(\mathrm{cm^{2}\,V^{-1}\,s^{-1}}\) (\(1\ \mathrm{m^{2}} = 10^{4}\ \mathrm{cm^{2}}\)): \[ \mu \approx 5.3\times10^{-4}\ \mathrm{cm^{2}\,V^{-1}\,s^{-1}}. \]
Esse valor está na faixa típica relatada para semicondutores orgânicos desordenados. (Se o seu \({t_\text{max}}\) mudar com a taxa de rampa ou a iluminação, você verá \(\mu\) mudar de acordo.)
👉 Próximo passo: Vá para a próxima seção da simulação CELIV (C-LIV) — Parte B — onde iremos: ajustar a taxa da rampa de tensão, experimentar a intensidade de iluminação, e explorar como alterar as mobilidades de elétrons/lacunas remodela o transiente CELIV.