Célula Solar de Silício Amorfo (a-Si:H) (1D) — p/i/n de Filme Fino (Limitada por defeitos)

1. Introdução

Células solares de silício amorfo hidrogenado (a-Si:H) são uma tecnologia fotovoltaica clássica de filme fino. Diferentemente do silício cristalino, o a-Si:H é estruturalmente desordenado e contém uma alta densidade de estados de defeito e de cauda. Como resultado, o comportamento do dispositivo é tipicamente limitado por recombinação de defeitos (processos do tipo Shockley–Read–Hall) e limitado por transporte devido a mobilidades menores. Isso faz do a-Si:H um excelente sistema para aprender como recombinação e densidade de portadores determinam V_oc, e como a óptica de filme fino molda J_sc.

Este tutorial fornece um fluxo de trabalho prático, baseado em física, para simular uma célula solar de a-Si:H em 1D usando OghmaNano. O modelo resolve o transporte de portadores por deriva–difusão acoplado à eletrostática de Poisson, com geração óptica dependente da profundidade e recombinação fisicamente motivada dominada por recombinação do tipo SRH por defeitos. O objetivo é conectar saídas padrão — curvas JV, Voc e comportamento Suns–Voc — a variáveis internas (bandas, níveis quase-Fermi, carga, e tempo efetivo de recombinação).

Você irá construir e simular uma célula solar unidimensional de a-Si:H de filme fino com uma arquitetura padrão p/i/n (veja ??). O dispositivo é definido como: Estrutura (frente → trás): p a-Si:H / i a-Si:H / n a-Si:H / contato metálico. Usando esse dispositivo de base, você gerará uma curva JV iluminada, extrairá métricas-chave de desempenho e executará varreduras de Voc e Suns–Voc para identificar perda de tensão limitada por recombinação e sua dependência da densidade de portadores.

2. Criando uma nova simulação

Para começar, crie uma nova simulação a partir da janela principal do OghmaNano. Clique no botão New simulation na barra de ferramentas. Isso abre a caixa de diálogo de seleção do tipo de simulação (veja ??).

Na caixa de diálogo do tipo de simulação, dê duplo clique em Si demos, depois dê duplo clique em Amorphous silicon (a-Si:H) (veja ??). O OghmaNano carregará uma simulação predefinida de célula solar de filme fino em a-Si:H.

A janela principal (veja ??) fornece uma visualização tridimensional da estrutura do dispositivo. Você pode usar o mouse para rotacionar, deslocar e ampliar a cena para inspecionar a geometria. Embora o presente tutorial use um modelo elétrico unidimensional, a visualização 3D oferece uma forma conveniente de visualizar a pilha de filme fino e os contatos.

Clique na aba Layer editor para abrir a tabela de camadas (veja ??). Aqui você pode inspecionar a estrutura vertical do dispositivo, incluindo a camada tipo p, o absorvedor intrínseco (camada i), a camada tipo n e o contato traseiro, juntamente com suas espessuras e materiais atribuídos.

3. Examinando o perfil de dopagem

O perfil de dopagem define a estrutura eletrostática do dispositivo antes que iluminação e polarização sejam aplicadas. Em uma célula solar de silício amorfo hidrogenado (a-Si:H) de filme fino, a dopagem é usada principalmente para formar contatos seletivos de portadores em vez de criar uma ampla região de depleção dentro de um absorvedor dopado. O princípio central de operação é o campo elétrico interno ao longo de uma camada intrínseca, que separa e transporta os portadores fotogerados.

Para visualizar a dopagem no OghmaNano, abra a faixa Electrical na janela principal e clique em Doping / Ions (veja ??). Isso abre o editor de perfil (veja ??), que plota a densidade de carga de doadores e aceitadores versus profundidade e lista os valores numéricos atribuídos a cada camada.

O dispositivo usa uma arquitetura padrão de a-Si:H p+/i/n+. Uma camada fina e fortemente dopada p+ é colocada na superfície frontal iluminada (densidade de aceitadores da ordem de \(10^{24}\,\mathrm{m^{-3}}\)), fornecendo um contato seletivo para lacunas e definindo a curvatura de bandas na frente. Em seguida vem uma camada absorvedora de a-Si:H intrínseca muito mais espessa, que é apenas muito levemente dopada (aqui \(\sim 10^{20}\,\mathrm{m^{-3}}\)) e contém a maior parte da geração óptica. Na parte traseira do dispositivo, uma camada fina e fortemente dopada n+ (densidade de doadores novamente \(\sim 10^{24}\,\mathrm{m^{-3}}\)) forma um contato traseiro seletivo para elétrons adjacente ao metal.

A característica-chave a reconhecer no editor de perfil é que o forte campo elétrico interno se estende pela camada intrínseca, e não pelas regiões dopadas. Os portadores fotogerados são criados principalmente no a-Si:H intrínseco e são separados por esse campo interno: elétrons são conduzidos ao contato traseiro n+, enquanto lacunas são conduzidas ao contato frontal p+. As camadas dopadas são mantidas finas para minimizar recombinação assistida por defeitos, ao mesmo tempo que fornecem boa seletividade elétrica e condutividade de contato.

4. Parâmetros elétricos: transporte, eletrostática e recombinação em a-Si:H

Abra o editor de parâmetros elétricos a partir da janela principal: Device structure → Electrical parameters. Os parâmetros são definidos por região usando as abas p+, i e n+, correspondendo à camada frontal de contato, ao absorvedor intrínseco e à camada traseira de contato da pilha de a-Si:H. As configurações neste editor determinam os coeficientes de transporte por deriva–difusão, o alinhamento energético, e os mecanismos de recombinação usados durante as varreduras JV e Suns–Voc.

O transporte em silício amorfo é caracterizado por baixas mobilidades devido à desordem e ao espalhamento forte. A demonstração, portanto, usa mobilidades muito menores do que no silício cristalino. Na camada intrínseca (veja ??), a mobilidade de elétrons é definida como μ_n = 7×10⁻⁴ m²V⁻¹s⁻¹ e a mobilidade de lacunas como μ_p = 7×10⁻⁶ m²V⁻¹s⁻¹. As camadas adjacentes aos contatos usam mobilidades diferentes (tipicamente menores) (veja ?? e ??), refletindo o fato de que as camadas dopadas/seletivas atuam principalmente como regiões finas de extração em vez de absorvedores de alta mobilidade. Em deriva–difusão, mobilidade reduzida aumenta o tempo de residência dos portadores e torna a recombinação mais consequente próximo de circuito aberto.

A base eletrostática usada em toda a pilha é visível em cada aba: afinidade eletrônica χ = 4.05 eV, gap de banda E_g = 1.75 eV, e permissividade relativa ε_r = 11. O gap mais largo reflete o gap óptico do a-Si:H e define o limite superior da fototensão atingível, enquanto a afinidade eletrônica e a permissividade definem o alinhamento de bandas e a penetração do campo. Juntamente com o perfil de dopagem p+/i/n+ da Seção ??, esses parâmetros estabelecem um forte campo interno ao longo da camada intrínseca, que é o principal mecanismo de separação de portadores em dispositivos p–i–n de filme fino.

A recombinação em a-Si:H é dominada por processos mediados por defeitos, então o modelo central de perda é a recombinação Shockley–Read–Hall (SRH) através de estados de armadilha localizados. Isso é habilitado em todas as regiões por meio do bloco Equilibrium SRH traps (veja os campos de parâmetros SRH em ??, ??, e ??). Nesta formulação, a densidade de armadilhas N_t, as seções de captura σ_n e σ_p, e o nível de energia de armadilha E_t definem as taxas de captura de portadores e, portanto, o tempo de vida efetivo: \[ \tau_n = \frac{1}{\sigma_n v_{\mathrm{th}} N_t}, \qquad \tau_p = \frac{1}{\sigma_p v_{\mathrm{th}} N_t}, \] onde \(v_{\mathrm{th}}\) é a velocidade térmica. A taxa de recombinação SRH correspondente assume a forma padrão \[ R_{\mathrm{SRH}} = \frac{np - n_i^2}{\tau_p (n+n_1) + \tau_n (p+p_1)}, \] com \(n_1\) e \(p_1\) determinados pelo nível de energia da armadilha \(E_t\). Como a maior parte da fotogeração ocorre no absorvedor intrínseco, a recombinação SRH na camada i é o mecanismo primário que limita a divisão dos níveis quase-Fermi e, portanto, V_oc neste tutorial.

Parâmetros de recombinação Auger também estão presentes nas abas p+ e n+ (veja ?? e ??), com coeficientes C_n = C_p = 1×10⁻⁴³ m⁶s⁻¹. A recombinação Auger representa processos de perda de alta densidade envolvendo três partículas da forma \[ R_{\mathrm{Auger}} = C_n\,n^2 p + C_p\,p^2 n, \] que se tornam relevantes quando as densidades de portadores são altas, como em regiões de contato fortemente dopadas ou sob iluminação extrema em simulações Suns–Voc. Em operação padrão de 1 sol em silício amorfo, a recombinação Auger é tipicamente um efeito secundário; a perda de tensão é dominada pela recombinação Shockley–Read–Hall (SRH) mediada por defeitos. Incluir termos Auger fornece, ainda assim, um canal fisicamente consistente em alta injeção e evita acumulação não física de portadores quando o dispositivo é levado a densidades muito altas (por exemplo, em centenas a milhares de sóis). Uma discussão detalhada da recombinação Auger e de sua implementação é dada em Teoria da recombinação Auger.

Em conjunto, as baixas mobilidades de portadores, o gap efetivo relativamente largo e a alta densidade de estados localizados de defeito colocam o a-Si:H firmemente em um regime de operação limitado por recombinação. Como resultado, tanto a curva JV quanto a resposta Suns–Voc são governadas principalmente por como o tempo de vida SRH efetivo colapsa com o aumento da densidade de portadores. Este tutorial, portanto, interpreta o comportamento de tensão diretamente em termos de divisão dos níveis quase-Fermi e recombinação controlada por defeitos, em vez de assumir hipóteses típicas do silício cristalino de alta mobilidade e operação limitada por transporte. Uma discussão detalhada da física de recombinação SRH e de seu papel na modelagem de dispositivos semicondutores é fornecida na seção teórica Teoria da recombinação Shockley–Read–Hall.

Para clareza e transparência numérica, este tutorial modela a recombinação por defeitos usando um único nível efetivo de armadilha SRH com seções de captura fixas. Esta é uma escolha de modelagem deliberada para o presente tutorial: ela captura a física dominante de recombinação relevante para células solares de a-Si:H, ao mesmo tempo que mantém compacto o espaço de parâmetros. O OghmaNano também suporta tratamentos mais avançados nos quais uma distribuição de armadilhas armazena carga ativamente e evolui dinamicamente sob polarização e iluminação. Tais tratamentos de não equilíbrio tornam-se essenciais quando se considera comportamento no domínio do tempo, quando o carregamento do dispositivo é importante ou quando é necessário acordo quantitativo com dados experimentais transientes ou com histerese. Essas extensões e sua motivação física são descritas em Por que armadilhas são necessárias em semicondutores desordenados e Recombinação Shockley–Read–Hall de não equilíbrio.

5. Perfil de geração óptica

Esta demonstração usa um cálculo óptico unidimensional para gerar um termo-fonte dependente da profundidade para o solucionador de deriva–difusão. Na prática, o modelo óptico calcula como o espectro AM1.5 incidente se propaga pela pilha de filme fino, quanto é absorvido em função do comprimento de onda e da profundidade, e converte essa potência absorvida em uma taxa de geração de pares elétron–lacuna. Em a-Si:H de filme fino, o absorvedor normalmente tem apenas algumas centenas de nanômetros de espessura, de modo que o perfil óptico é fortemente acoplado ao projeto da pilha: grande parte da absorção útil ocorre próxima ao lado iluminado, e estratégias de retroreflexão/aprisionamento de luz afetam fortemente J_sc.

Para visualizar (e regenerar) a solução óptica, vá para a faixa Optical na janela principal e clique em Transfer matrix. Isso abre o editor da simulação óptica. Pressione o botão azul Run optical simulation (ícone de play) para calcular o campo óptico e atualizar os gráficos (veja ??).

O gráfico em ?? é um mapa comprimento de onda–profundidade. O eixo horizontal é a profundidade no dispositivo (rotulado como posição y), e o eixo vertical é o comprimento de onda. Cores brilhantes próximas à superfície iluminada indicam alta população de fótons; o enfraquecimento com a profundidade indica absorção dentro da pilha.

Para a-Si:H, a absorção no visível é forte, então uma grande fração dos fótons úteis é absorvida em uma distância relativamente curta a partir da frente. O limite de comprimento de onda mais longo é controlado pelo gap efetivo do a-Si:H: fótons abaixo do gap não geram portadores de forma eficiente, de modo que o mapa de fótons mostrará um limite espectral além do qual absorção e geração caem fortemente. Essa é uma das razões pelas quais dispositivos de a-Si:H de filme fino normalmente trocam J_sc absoluto por gap mais alto e melhor seletividade espectral em comparação com o silício cristalino.

O editor óptico fornece várias abas porque “densidade de fótons” não é a mesma coisa que “taxa de geração”. A visualização Photon distribution diz onde o campo óptico existe. As abas Photon distribution absorbed e Generation rate são as que importam para o modelo elétrico: o mapa de fótons absorvidos é convertido em uma taxa de geração de pares elétron–lacuna dependente da profundidade, que é o termo-fonte injetado nas equações de deriva–difusão. Para esta célula de a-Si:H, o perfil de geração resultante é fortemente concentrado na frente e confinado ao absorvedor fino, que é exatamente o regime em que recombinação e transporte próximos de circuito aberto se tornam decisivos para V_oc.

6. Executando a simulação, curvas JV e extração de parâmetros

Quando a estrutura do dispositivo, os parâmetros elétricos e a geração óptica estiverem definidos, a simulação pode ser executada diretamente da janela principal. Clique no botão Run simulation na barra de ferramentas para iniciar o solucionador. Durante a execução, o progresso do solucionador e as informações de convergência são gravados na janela do terminal (veja ??).

Para cada ponto de polarização, a tensão aplicada no contato superior é listada primeiro, seguida pela densidade de corrente resultante. Sob iluminação, a corrente é inicialmente negativa (geração de potência). À medida que a tensão aplicada aumenta, a magnitude da corrente diminui até cruzar zero na tensão de circuito aberto. Além desse ponto, a corrente se torna positiva, correspondendo à operação de diodo em polarização direta. A saída do terminal também informa o residual do solucionador e o tempo de convergência para cada passo de polarização. Residuais pequenos indicam que as equações acopladas de deriva–difusão e Poisson estão sendo resolvidas com precisão.

Para inspecionar as características elétricas, abra a aba Output e dê duplo clique em jv.csv. Isso exibe a curva corrente-densidade versus tensão (JV) (veja ??). A curva JV é o principal diagnóstico do comportamento do dispositivo: ela deve passar pela corrente de curto-circuito em polarização zero e cruzar corrente zero na tensão de circuito aberto.

Dar duplo clique em siminfo.dat abre a janela de informações da simulação (veja ??), que informa métricas extraídas de desempenho, incluindo fator de forma, eficiência de conversão de potência, ponto de máxima potência, V_oc e J_sc. Quantidades diagnósticas adicionais, como densidades de portadores livres em circuito aberto, também são úteis para interpretar comportamento limitado por defeitos em a-Si:H.

Uma regra prática é sempre inspecionar a curva JV antes de interpretar as métricas numéricas. Se a curva JV não passar de forma limpa por curto-circuito e circuito aberto, ou se a corrente tiver um sinal ou forma inesperados, as quantidades derivadas em siminfo.dat também serão pouco confiáveis. A inspeção visual da curva JV é a forma mais rápida de diagnosticar problemas de configuração ou modelagem.

7. Análise Suns–Voc: tensão limitada por recombinação em a-Si:H

Uma medida de Suns–Voc investiga como a tensão de circuito aberto evolui em função da intensidade de iluminação sob condições estritas de corrente zero. Como a corrente no terminal é forçada a se anular, Suns–Voc isola a física de recombinação de efeitos de transporte e resistência de contato. A curva resultante é, portanto, um dos diagnósticos mais diretos dos mecanismos de perda de tensão em uma célula solar.

No OghmaNano, Suns–Voc é implementado como um modo dedicado de simulação. Em vez de varrer a tensão aplicada, o solucionador varre a intensidade de iluminação e, em cada intensidade, ajusta a tensão terminal até que a corrente líquida seja zero. Isso calcula diretamente o ponto de operação em circuito aberto em cada nível de iluminação.

Para habilitar Suns–Voc, abra a faixa Simulation type na janela principal (veja ??) e selecione Suns–Voc. Em seguida, clique em Run simulation. O solucionador grava automaticamente os resultados em disco; o arquivo suns_voc.csv contém V_oc em função da intensidade de iluminação.

Em baixa iluminação, o aumento de V_oc é governado pelo acúmulo de densidade excessiva de portadores no absorvedor intrínseco. No limite não degenerado, a divisão dos níveis quase-Fermi é \[ qV_{\mathrm{oc}} = E_{Fn} - E_{Fp} = kT \ln\!\left(\frac{np}{n_i^2}\right), \] de modo que V_oc aumenta apenas logaritmicamente com a densidade de portadores.

Em silício amorfo, esse regime de baixa intensidade é frequentemente modificado por efeitos de ocupação de defeitos. Em taxas de geração muito baixas, uma fração substancial dos portadores fotogerados é capturada por estados profundos e de cauda antes de contribuir para as populações de portadores livres. Isso pode produzir uma resposta inicial fraca ou ligeiramente achatada na curva Suns–Voc, como observado em ??. Esse comportamento é uma consequência física do preenchimento de armadilhas e não um artefato numérico.

Para explorar o regime de alta injeção, estenda o intervalo de iluminação. Abra o Suns–Voc editor a partir da faixa Editors e aumente a stop intensity para 1000 suns. Execute novamente a simulação.

Abrir suns_voc.csv após executar novamente a simulação revela o comportamento em toda a faixa de iluminação (veja ??). Em alta iluminação, V_oc continua a aumentar, mas com inclinação progressivamente menor, aproximando-se eventualmente da saturação.

A origem da saturação de tensão se torna clara quando a curva Suns–Voc é visualizada juntamente com a carga armazenada e o tempo de vida efetivo. À medida que a iluminação aumenta, a densidade total de portadores em excesso cresce (veja ??), mas o tempo efetivo de recombinação simultaneamente diminui (veja ??).

Em a-Si:H, a recombinação é dominada por processos do tipo SRH mediados por defeitos. À medida que a densidade de portadores aumenta, a recombinação assistida por armadilhas acelera e o tempo de vida efetivo \[ \tau_{\mathrm{eff}} = \frac{\Delta n}{R} \] cai. A partir desse ponto, fotogeração adicional aumenta principalmente a recombinação em vez de aumentar a separação dos níveis quase-Fermi.

Esse equilíbrio entre geração e recombinação impõe um limite intrínseco superior a V_oc. A tensão máxima permanece abaixo do gap efetivo do a-Si:H porque a recombinação impede que os níveis quase-Fermi de elétrons e lacunas se aproximem simultaneamente de suas respectivas bordas de banda. A diferença E_g/q − V_oc representa, portanto, a perda fundamental de tensão induzida por recombinação neste dispositivo.