Simulação de Célula Solar Heteroestrutural AlGaAs/GaAs: Estrutura do Dispositivo, Alinhamento das Bandas de Energia, Curva J–V, Eficiência Quântica Externa e Relação entre a Irradiância e a Tensão de Circuito Aberto.

1. Introdução: por que usar uma heterojunção AlGaAs/GaAs?

O arsenieto de gálio (GaAs) é um semicondutor do grupo III–V com banda de energia direta, amplamente utilizado em dispositivos fotovoltaicos de alto desempenho. A sua capacidade de absorção relativamente elevada acima do limite da banda permite que uma camada fina absorva uma fração considerável da luz solar incidente em comprimentos de onda adequados. No entanto, uma junção homojunção convencional de GaAs pode sofrer absorção ótica indesejada e recombinação de portadores minoritários perto da superfície iluminada. Um semicondutor com uma banda de energia mais ampla, colocado no lado iluminado, pode reduzir essas perdas quando as suas propriedades óticas e eletrónicas são escolhidas de forma adequada.

Este tutorial apresenta um modelo de uma heterojunção p-Al0,3Ga0,7As / n-GaAs. As duas camadas superiores de AlGaAs formam uma janela/emissor de tipo p com uma banda proibida mais ampla, enquanto a camada inferior de GaAs de tipo n fornece a principal região de absorção ótica. A diferença na banda proibida e na afinidade eletrónica na sua interface introduz descontinuidades nas bandas de condução e de valência, para além da curvatura eletrostática produzida pelos dopantes ionizados. A heterojunção resultante deve, portanto, ser analisada utilizando tanto os deslocamentos de banda do material quanto os níveis de quasi-Fermi dos portadores; tratá-la como uma simples junção PN de GaAs deixaria de considerar aspetos físicos importantes.

Utilizaremos o solucionador de difusão-deriva unidimensional da OghmaNano, acoplado a um cálculo de geração óptica. O modelo fornecido inclui uma varredura J–V sob iluminação, opções para cálculos dependentes de EQE e intensidade, e perfis internos salvos. O fluxo de trabalho é intencionalmente baseado num exemplo já configurado, de modo que o objetivo é compreender, executar e analisar o dispositivo, em vez de inserir novamente todos os parâmetros desde o início.

Num dispositivo fotovoltaico, os fotões absorvidos criam pares eletrão-lacuna; os contactos seletivos extraem os portadores de carga e a divisão resultante do nível quase-Fermi interno produz uma tensão de circuito aberto. Um modelo útil deve estabelecer uma ligação entre a absorção ótica e a geração de carga e a eletrostática da junção, a extração e a recombinação dos portadores de carga. As secções que se seguem examinam diretamente cada uma dessas ligações.

2. Criar a simulação a partir dos exemplos de demonstração do GaAs.

Abra o OghmaNano e selecione «Nova simulação». No navegador de simulações, abra as demonstrações de GaAs e, em seguida, selecione a célula solar de heterojunção AlGaAs/GaAs. A caixa de diálogo de seleção e a categoria da demonstração são apresentadas em ?? e ??. O exemplo é aberto na janela principal da estrutura do dispositivo (??).

A janela principal contém uma representação gráfica tridimensional do dispositivo. A simulação física descrita aqui utiliza a variação espacial ao longo do eixo y, perpendicular às camadas. A placa renderizada é, portanto, uma forma de visualizar a estrutura e os contactos, e não uma evidência de uma solução elétrica tridimensional completa. Selecione a faixa de opções «Estrutura do dispositivo» para inspecionar ou ajustar as camadas; os separadores «Saída» e «Terminais» serão utilizados após a execução.

É útil guardar uma cópia separada do exemplo antes de realizar os testes. Mantenha a versão original intacta para que as alterações no processo de dopagem, na proporção da liga ou nos parâmetros da interface possam ser comparadas com uma referência de trabalho.

3. Inspecionar a estrutura epitaxial, a dopagem e a composição da liga.

Abra o editor de camadas para examinar as quatro regiões epitaxiais (??). A geometria inicial e a dopagem exatas deste exemplo estão resumidas abaixo. Todas as concentrações são expressas como densidades volumétricas em m−3, conforme armazenado no OghmaNano; para converter para cm−3, divida por 106.

A heterointerface encontra-se em y = 0,40 µm, onde o emissor de AlGaAs do tipo p se encontra com o absorvedor de GaAs do tipo n. Isto representa, simultaneamente, uma mudança na polaridade do dopante majoritário e uma mudança na composição do semicondutor. Em contraste, a interface de 0,10 µm no interior do AlGaAs e a interface de 3,20 µm no interior do GaAs representam, principalmente, mudanças na concentração de dopagem; essas fronteiras não introduzem um novo sistema de materiais.

Dopagem aberta/Íons provenientes dos controlos elétricos (??). O gráfico deve mostrar a elevada densidade de aceitadores perto do contacto iluminado, a menor densidade de aceitadores no restante emissor de AlGaAs, o absorvedor de GaAs dopado com doadores e o contacto traseiro de GaAs com maior concentração de dopantes. A mudança de sinal na dopagem líquida identifica a junção PN. A dopagem não é, em si, um parâmetro de absorção ótica; ela entra na equação de Poisson através da carga espacial ionizada e afeta a condutividade, a recombinação e a recolha de portadores.

Em seguida, verifique a composição da liga (??). Ambas as camadas de AlGaAs apresentam g_x0 = g_x1 = 0.30: uma fração constante de alumínio em toda a sua espessura. Ambas as camadas de GaAs apresentam g_x0 = g_x1 = 0 neste exemplo específico. Não existe uma variação gradual na composição na estrutura fornecida. Uma região de AlGaAs com uma variação gradual linear exigiria pontos finais de composição diferentes e um modelo de material Lua capaz de avaliar a composição local em cada ponto da malha.

O efeito da janela de banda larga depende do comprimento de onda. Os fotões com energias inferiores ao limiar de absorção do AlGaAs conseguem atravessar mais facilmente o semicondutor superior, enquanto os fotões com energias superiores a esse limiar também podem ser absorvidos no emissor. A magnitude deste benefício depende das constantes óticas reais, da espessura da camada e das propriedades de coleta interna; a fração de Al, por si só, não estabelece uma redução quantitativa das perdas na superfície frontal.

4. Modelos de materiais elétricos e os desfasamentos de banda na heterojunção.

Abra o editor de parâmetros elétricos e examine os scripts Lua dos materiais associados às regiões (?? e ??). Este exemplo avalia os parâmetros do material utilizando os scripts registados, em vez de definir cada parâmetro de banda individualmente com um valor fixo.

Os caminhos Lua para AlGaAs e GaAs são distintos: Inorganic/AlGaAs/AlGaAs.lua e Inorganic/GaAs/GaAs.lua. O valor da composição da camada superior, 0,30, está disponível para o modelo AlGaAs, permitindo o cálculo de grandezas elétricas dependentes da composição. Os conjuntos de dados óticos são configurados separadamente: Inorganic/AlGaAs/Adachi para as duas camadas de AlGaAs e Inorganic/GaAs/Adachi para as camadas de GaAs. Os parâmetros da banda de condução elétrica e os conjuntos de dados óticos de índice de refração/absorção não devem ser confundidos entre si.

Na heterointerface, a extremidade local da banda de condução está relacionada com a afinidade eletrónica através de uma convenção da forma \(E_C=-\chi-q\phi\), e a extremidade da banda de valência através de \(E_V=E_C-E_g\), desde que a mesma energia de referência seja utilizada em ambos os lados. As descontinuidades resultantes são convencionalmente escritas como \(\Delta E_C\) e \(\Delta E_V\). Elas dependem tanto da diferença de energia entre as bandas, que varia com a composição, como do modelo de afinidade eletrónica. Assim, a maior diferença de energia entre as bandas do AlGaAs não determina, por si só, como o deslocamento total é dividido entre as bandas de condução e de valência.

Fisicamente, um pico na banda de condução pode dificultar a extração de eletrões quando suficientemente grande, enquanto uma descontinuidade adequada na banda de valência pode dificultar o movimento indesejado de lacunas na direção oposta. Estes efeitos competem com a deriva, a difusão, a dopagem, a ativação térmica e a recombinação. É por isso que a estrutura de bandas simulada fornece mais informações do que simplesmente indicar uma diferença no intervalo de energia.

4.1. Parâmetros de transporte e recombinação.

O OghmaNano resolve as equações de Poisson e de continuidade de elétrons e lacunas, que estão acopladas. Na aproximação de deriva e difusão, as densidades de corrente de elétrons e lacunas dependem da mobilidade dos portadores de carga, dos gradientes de densidade dos portadores e do campo eletrostático. Na convenção de sinais usual, estes podem ser expressos como:

O campo elétrico e as concentrações de portadores de carga são calculados de forma autoconsistente. Sob iluminação, as equações de continuidade devem também considerar a fotogeração \(G\) e a recombinação \(R\). Por exemplo, utilizando a convenção usual de corrente unidimensional, \(dJ_n/dy=q(R-G)\) e \(dJ_p/dy=q(G-R)\) em estado estacionário. A sua soma deve, consequentemente, permanecer espacialmente constante numa solução ideal de estado estacionário, sem a necessidade de um termo adicional de armazenamento de carga distribuída.

O formato JSON permite definir o estado estacionário da recombinação SRH (ss_srh_enabled = "True") e fornece uma energia de armadilha nominal no meio da banda proibida, a densidade de armadilhas \(N_t=5\times10^{19}\,\mathrm{m}^{-3}\) e as seções de captura \(\sigma_n=\sigma_p=10^{-22}\,\mathrm{m}^2\) em cada camada. Estes campos descrevem o modelo de defeitos configurado; deve-se verificar se um campo é obtido diretamente do JSON ou recalculado pelo script de material escolhido, comparando-o com os valores efetivos reportados pelo OghmaNano. A taxa de um processo SRH padrão de um único nível é:

Num semicondutor com banda proibida direta, a recombinação radiativa de eletrões livres em lacunas livres também pode ser importante, sendo convencionalmente representada por \(R_{\mathrm{rad}}=B(np-n_i^2)\). O programa também pode fornecer coeficientes de Auger e outros canais através dos modelos de materiais. Não se deve inferir qual mecanismo é dominante apenas pela presença das suas configurações: utilize os ficheiros de recombinação gerados, as opções ativadas e os valores do modelo de materiais no exemplo em execução. Para este dispositivo, a interpretação adequada é obtida comparando as taxas no mesmo ponto de tensão e utilizando escalas verticais consistentes.

A vista do editor elétrico pode ser utilizada para analisar a mobilidade e outras propriedades locais (??). Ao experimentar com diferentes composições, mantenha as opções de substituição da base de dados de banda proibida e afinidade eletrónica desativadas, a menos que exista uma razão específica para substituir os valores do Lua: valores fixos manualmente \(E_g\) ou \(\chi\) podem, inadvertidamente, eliminar a dependência da composição pretendida.

5. Contatos elétricos e a malha de simulação

O modelo possui dois terminais elétricos externos. O ânodo p superior, em y0, é o contato de condução; o cátodo n inferior, em y1, é o terminal de referência (terra). Ambos utilizam uma barreira de maioria de portadores com comportamento ôhmico e uma barreira de minoria de portadores com comportamento de bloqueio. Esta é uma escolha importante no modelo: a estrutura não é representada por contatos simétricos genéricos que extraem eletrões e lacunas de forma indiscriminada. Abra o editor de contatos para visualizar as linhas de contato (??).

No ânodo, a maioria das cargas são lacunas e a minoria são eletrões. No cátodo, os papéis são invertidos. Assim, a aplicação de uma tensão variável altera o potencial no terminal p em relação ao terminal n aterrado. O modelo de contato é distinto das camadas de contato epitaxiais altamente dopadas, que fazem parte da geometria do semicondutor. As condições explícitas dos terminais determinam como as cargas abandonam o domínio computacional.

O solucionador elétrico está configurado para uma malha unidimensional no eixo y. Os seus quatro segmentos de malha acompanham as quatro espessuras epitaxiais, com um número nominal de pontos de 10, 17, 156 e 11, respetivamente. Isto corresponde a aproximadamente 194 pontos especificados, embora o número de nós únicos efetivamente utilizados possa depender da forma como o solucionador trata as fronteiras partilhadas. Ambas as direções da malha elétrica lateral estão desativadas. As configurações da malha podem ser visualizadas em ??.

A resolução da malha é especialmente importante perto da heterointerface, a 0,40 µm, onde a deflexão eletrostática e as diferenças de banda podem variar numa curta distância. Se a descontinuidade de banda ou a densidade de portadores calculadas apresentarem um pico estreito, uma verificação da convergência da malha é mais convincente do que simplesmente aceitar a configuração padrão. Repita o cálculo utilizando uma malha mais fina para o emissor/junção e compare os resultados de corrente-tensão, bem como os perfis espaciais relevantes. Se estes apresentarem alterações significativas, a malha de referência ainda não estava suficientemente refinada.

6. Executar a simulação J–V com iluminação.

O modo selecionado por padrão na configuração fornecida é J–V. A variação da tensão é configurada de -0,10 V a +1,00 V, em incrementos de 0,02 V. A iluminação utiliza o espectro AM1.5G e a escala de intensidade da fonte de luz de referência corresponde a uma unidade solar (Psun = 1). A geração óptica é acoplada à simulação elétrica através do modelo de geração da matriz de transferência configurado. O valor definido para a temperatura é de 300 K.

Antes de executar a simulação, examine a guia «Resultados» (??). Clique em «Executar simulação». A guia «Terminal» exibe o progresso e os resultados do solucionador; as capturas de tela do terminal correspondentes (?? e ??) mostram um estado inicial e um estado posterior do cálculo. A alteração da exibição no terminal indica que o solucionador está em funcionamento, mas, por si só, não comprova que todos os pontos convergiram. Verifique se a simulação foi concluída, se ocorreram erros e quais são os arquivos de resultados gerados.

Após a execução, volte para a secção «Resultados» e abra jv.csv (??). A curva J–V iluminada difere, em termos conceituais, da curva da díada em polarização direta e no escuro. Uma expressão simples de primeira ordem é \(J(V)=J_{\mathrm{dark}}(V)-J_{\mathrm{ph}}(V)\), para uma convenção de sinais em que a fotocorrente extraída é negativa. O sinal representado pelo simulador segue a convenção de contacto configurada, pelo que as quantidades fisicamente relevantes devem ser identificadas a partir dos valores de tensão e corrente, em vez de se assumir que todos os gráficos fotovoltaicos são desenhados com os mesmos sinais.

A densidade de corrente de curto-circuito, \(J_{SC}\), é a densidade de corrente no ponto \(V=0\), enquanto a tensão de circuito aberto, \(V_{OC}\), é a tensão na qual a densidade de corrente no terminal é zero. O ponto de máxima potência ocorre na tensão de polarização aplicada que maximiza a magnitude da potência fornecida a uma carga externa, e o fator de preenchimento é definido por:

Com uma densidade de potência ótica incidente \(P_{\mathrm{in}}\) consistentemente especificada, a eficiência de conversão segue o padrão \(\eta=|J_{\mathrm{mpp}}V_{\mathrm{mpp}}|/P_{\mathrm{in}}\). Não se deve indicar uma eficiência apenas com base na forma da curva: é também necessário considerar a convenção da potência de iluminação e um ponto de operação máximo de potência fiável e devidamente identificado.

A curva J–V é a forma mais simples de resumir o desempenho, mas não explica por que razão se verifica uma perda de corrente. Uma diminuição na fotocorrente pode resultar de uma absorção inadequada, de uma recolha deficiente de portadores ou de um aumento na recombinação; uma tensão de circuito aberto mais baixa pode resultar de uma recombinação mais intensa ou de alterações na eletrostática da junção. Por esse motivo, as secções seguintes utilizam dados espacialmente resolvidos como evidência diagnóstica, em vez de se basearem apenas no gráfico J–V.

A visualização device.csv (??) oferece uma ferramenta de verificação independente e útil para garantir que as regiões físicas simuladas e os seus limites correspondem à geometria pretendida. O painel de informações da simulação (??) é útil para comparar diferentes execuções e verificar as configurações associadas a um determinado resultado.

7. Diagramas de bandas, densidades de portadores e transporte através da heterojunção.

Abra as capturas de tela salvas no visualizador de resultados e mova o controle deslizante da captura de tela pelos pontos de ajuste. As quatro capturas de tela de cada quantidade mostradas abaixo representam diferentes estados salvos da mesma simulação; as suas escalas numéricas e os rótulos de ajuste devem ser lidos diretamente no visualizador correspondente. Ao comparar dois estados, determine primeiro se eles partilham os mesmos limites de eixo e qual a quantidade representada que está a ser comparada.

7.1. Bandas de condução e de valência, níveis de Fermi.

O diagrama de bandas é o local mais adequado para identificar a heterojunção AlGaAs/GaAs (??–??). A interface entre os materiais está localizada em 0,40 µm. Uma mudança abrupta nas extremidades da banda local nesse ponto reflete o desalinhamento dos materiais, enquanto a variação suave da curvatura em outros pontos reflete a eletrostática. Essa distinção é fundamental: uma heterojunção não pode ser interpretada apenas em termos de uma banda de energia do GaAs que se curva continuamente.

No equilíbrio térmico, os potenciais eletroquímicos dos eletrões e das lacunas coincidem, resultando num nível de Fermi comum. Fora do equilíbrio, os níveis de quasi-Fermi dos eletrões e das lacunas podem separar-se. A diferença entre eles é uma medida do desequilíbrio local dos portadores de carga e está relacionada com o produto local das populações de eletrões e lacunas. Em circuito aberto, a tensão nos terminais reflete a diferença entre os potenciais eletroquímicos estabelecidos nos contactos seletivos; não se deve confundir uma única separação local dos níveis de quasi-Fermi numa profundidade arbitrária com a tensão medida em circuito aberto.

À medida que a polarização aplicada se altera, o potencial eletrostático e a distribuição de portadores de carga respondem em conjunto. Sob polarização elétrica direta, a barreira da junção é reduzida em relação ao equilíbrio; durante o funcionamento fotovoltaico, os fotoportadores estabelecem populações fora do equilíbrio e uma fototensão. Qualquer variação aparente no traçado do nível de quasi-Fermi deve ser analisada em conjunto com o fluxo de corrente e a resolução da malha antes de ser interpretada como um potencial eletroquímico abrupto.

7.2. Densidades de eletrões e lacunas.

Os gráficos da densidade de portadores (??–??) permitem distinguir as regiões com a maioria dos portadores dos portadores minoritários injetados ou fotogerados. Através do emissor do tipo p, a espécie majoritária em equilíbrio são as lacunas; em todo o absorvedor do tipo n, são os eletrões. Os perfis dos portadores minoritários contêm informações importantes sobre a coleta e a recombinação, particularmente no lado de GaAs da junção.

As concentrações de eletrões e lacunas são sensíveis à densidade de estados efetiva, ao deslocamento das bandas, ao potencial eletrostático e ao equilíbrio entre geração e recombinação. As alterações na concentração não podem, portanto, ser atribuídas exclusivamente ao perfil de dopagem. Uma mudança abrupta perto da heterointerface pode resultar da combinação dos efeitos do alinhamento local das bandas e do transporte por deriva e difusão.

7.3. Densidades de corrente de eletrões e lacunas.

Os gráficos das correntes de eletrões e lacunas (??–??) mostram como a corrente total nos terminais é distribuída entre os dois tipos de portadores em função da profundidade. A fotogeração cria tanto eletrões quanto lacunas, mas a forma como cada um desses elementos contribui para a corrente varia ao longo do dispositivo, de acordo com a sua geração, recombinação, mobilidade e campo elétrico locais.

Num cálculo unidimensional em estado estacionário, a soma das correntes de eletrões e lacunas, considerando os seus sinais, deve ser constante em qualquer ponto que não esteja próximo de uma fonte de corrente distribuída ou de uma perda modelada explicitamente. No entanto, cada componente não precisa de ser constante: à medida que os portadores de carga são gerados e recombinam, os componentes da corrente de eletrões e lacunas trocam contribuições, mantendo a corrente total. Este é um método eficaz para verificar a consistência de uma solução completa.

Em alguns casos, determinados padrões podem apresentar variações acentuadas em condições de baixa corrente. Evite diagnosticar uma barreira física de transporte com base em um único traçado ruidoso: analise as extremidades da banda, as densidades de portadores e os componentes de corrente com sinal em conjunto e verifique se o padrão se altera quando a malha é refinada. O gráfico de mobilidade separado (??) oferece uma forma prática de visualizar os coeficientes de transporte local modelados e verificar se eles variam em função do material ou da composição.

8. Recombinação, taxas de portadores e fotogeração.

A corrente gerada pelas células solares é o resultado de uma competição entre a geração e a perda de carga. Os fotões absorvidos no dispositivo criam portadores de carga, mas apenas uma fração desses portadores consegue atingir os contactos seletivos. As equações locais de continuidade de eletrões e lacunas refletem este equilíbrio. A OghmaNano apresenta uma variedade de taxas de produção, que podem ser comparadas espacialmente, e não apenas através de uma corrente medida no próprio dispositivo.

8.1. Recombinação versus polarização aplicada

Abra os perfis de recombinação (??–??). A taxa de recombinação é uma taxa volumétrica, geralmente expressa em portadores por metro cúbico por segundo; não deve ser comparada numericamente com uma densidade de corrente sem integrar ao longo de uma distância e multiplicar pela carga elementar. Para um dispositivo unidimensional, uma densidade de corrente de recombinação ilustrativa é \(J_{\mathrm{rec}}=q\int R(y)\,dy\), sujeita à definição escolhida da taxa líquida e dos limites.

Compare não apenas as alturas máximas, mas também a largura e a localização das taxas. Um pico elevado e estreito pode contribuir menos para a perda total do que uma taxa moderada que se estende por todo o absorvedor. Utilize a mesma voltagem e região ao comparar os mecanismos. Numa heterojunção, o acúmulo de portadores próximo a uma descontinuidade de banda pode aumentar localmente a recombinação, mas o fator determinante é se esse acúmulo é realmente visível nos gráficos de densidade e de banda.

8.2. Taxas de recombinação eletrão-lacuna.

Os dados relacionados com as taxas de geração e aniquilação de pares elétron-lacuna (??–??) ajudam a comparar as taxas apresentadas no visualizador de imagens. É importante distinguir a recombinação líquida de qualquer processo individual de captura, emissão, geração ou transição de portadores livres. Nem todas as taxas representadas correspondem necessariamente a um mecanismo de perda independente; somar as curvas exibidas sem verificar as suas definições pode levar a uma dupla contagem dos caminhos de transferência.

No caso de um canal de recombinação radiativa, a lei bem conhecida \(R_{\mathrm{rad}}=B(np-n_i^2)\) ilustra por que a recombinação aumenta à medida que ambas as populações de portadores de carga aumentam. A recombinação de Shockley-Read-Hall (SRH), por outro lado, depende da ocupação das armadilhas e da cinética de captura. A importância relativa desses processos depende dos parâmetros efetivos do material utilizados na simulação. Este tutorial não assume que a recombinação radiativa seja o processo dominante em todas as camadas ou em todos os pontos de polarização, apenas porque o GaAs é um semicondutor de banda direta.

8.3. Perfil da geração.

A taxa de geração de portadoras ópticas é fundamentalmente diferente de uma taxa de recombinação: é uma fonte, e não uma perda. Para um fluxo de fotões dependente do comprimento de onda \(\Phi(\lambda,y)\) e um coeficiente de absorção \(\alpha(\lambda,y)\), uma relação local simples é \(G(y)=\int \alpha(\lambda,y)\Phi(\lambda,y)\,d\lambda\), sendo que o cálculo preciso depende da interferência, da reflexão e das definições de campo do modelo ótico escolhido.

Os perfis de geração armazenados (?? e ??) permitem visualizar onde os fotoportadores são produzidos. Uma geração mais intensa perto do lado iluminado pode resultar da absorção exponencial; a interferência e a absorção dependente do comprimento de onda podem produzir uma dependência da profundidade mais complexa. Um portador gerado opticamente não é automaticamente um portador coletado, razão pela qual a geração deve ser analisada em conjunto com os gráficos de corrente e recombinação.

9. Eficiência quântica externa: coleta de portadores dependente do comprimento de onda.

A eficiência quântica externa (EQE) mede o número de eletrões coletados através do circuito externo por cada fotão incidente no dispositivo, num comprimento de onda especificado, tendo em conta o sinal e a normalização utilizados na saída. Num cenário ideal de recolha de um eletrão por fotão, \(EQE(\lambda)=J_{\mathrm{photo}}(\lambda)/(q\Phi_{\mathrm{incident}}(\lambda))\), com a corrente espectral e o fluxo de fotões expressos de forma consistente.

Para abrir o cálculo do EQE, selecione a opção de simulação de EQE no OghmaNano e execute-a. Em seguida, examine o resultado guardado em «Resultados» (??). O ficheiro de simulação fornecido também contém um comprimento de onda de teste de EQE configurado de 710 nm; esta configuração não deve ser confundida com um espetro de EQE de banda larga completo, medido ou calculado, a menos que tenha sido realmente realizada uma varredura de comprimentos de onda e o resultado contenha esse espetro.

No caso do GaAs, as energias dos fotões próximas e superiores à sua banda proibida correspondem à extremidade do espectro infravermelho próximo/visível da resposta fotovoltaica. O emissor de AlGaAs, com uma banda proibida mais ampla, pode ser mais transparente em parte da faixa de absorção útil do GaAs, aumentando a fração de luz disponível para o absorvedor. Em comprimentos de onda muito curtos, a absorção no emissor de AlGaAs frontal e a recombinação de portadores gerados perto da superfície podem reduzir a eficiência de coleta. Em comprimentos de onda mais longos, próximos à extremidade de absorção do GaAs, a absorção incompleta num absorvedor de espessura finita de 2,8 µm pode tornar-se um fator importante. As conclusões espectrais quantitativas devem ser obtidas a partir da curva de eficiência quântica externa (EQE) calculada e de conjuntos de dados óticos verificados, e não apenas inferidas a partir da banda proibida.

A eficiência quântica externa (EQE) não é equivalente à absorvância ótica: a absorvância contabiliza os fotões absorvidos; a EQE contabiliza os portadores coletados. A comparação entre ambas revela perdas na coleta elétrica. Se a absorção for alta, mas a EQE for baixa em um determinado comprimento de onda, é necessário analisar o perfil de geração espacial e o comportamento de recombinação/transporte da camada na qual esses fotões estão sendo absorvidos.

10. Cálculos da intensidade da iluminação para os modelos Suns–Voc e Suns–Jsc.

10.1. Suns – Voc

As sondas Suns–Voc investigam a dependência da tensão de circuito aberto em relação à intensidade da iluminação. O circuito aberto é definido pela ausência de corrente externa líquida nos terminais, pelo que os portadores de carga fotogerados acumulam uma tensão em vez de fornecerem energia a uma carga. O exemplo fornecido contém uma configuração Suns–Voc que abrange um multiplicador de intensidade de 0,01 a 1000 unidades de iluminação, utilizando incrementos multiplicativos de 1,2. Acesse este modo de simulação e examine a seleção de resultados (??) e a curva (??).

Num modelo idealizado de díodo, com uma física de recombinação aproximadamente independente da iluminação, a relação \(V_{OC}\approx(n_{\mathrm{id}}k_BT/q)\ln(J_{\mathrm{ph}}/J_0+1)\) prevê um aumento aproximadamente logarítmico da tensão de circuito aberto com a intensidade da luz. Aqui, \(n_{\mathrm{id}}\) é um fator de idealidade efetivo. Num dispositivo de heterojunção, a distribuição local de portadores, os desníveis de banda, os canais de recombinação e o comportamento em condições de alta injeção podem todos variar com a intensidade, pelo que a inclinação efetiva não precisa de ser constante em toda a faixa de variação.

O parâmetro Suns–Voc é útil para analisar a resposta implícita do diodo em condições de corrente nula. Por si só, não inclui todas as perdas de tensão causadas pela passagem de corrente através de um contacto resistivo ou de uma camada de transporte em condições de geração de energia. Assim, a comparação do comportamento de Suns–Voc com a curva J–V sob iluminação pode ajudar a distinguir as limitações relacionadas com a recombinação das perdas de transporte dependentes da corrente, embora a interpretação de tal diferença exija a consistência da iluminação, dos modelos de recombinação e das definições de tensão.

10.2. Suns–JSC

O Suns–Jsc analisa a corrente de curto-circuito em relação à intensidade da iluminação. O exemplo está configurado para uma varredura que começa em 0,01 unidades de iluminação solar e termina em 10 unidades, com um fator de intensidade multiplicativo de 1,05 e o incremento inicial armazenado sundp = 0.01. Utilize o seletor de simulação Suns–Jsc e verifique o resultado em ?? e ??.

No regime mais simples de eficiência constante na coleta, a corrente de curto-circuito varia aproximadamente de forma linear com o fluxo de fotões incidente, portanto, \(|J_{SC}|\propto P_{\mathrm{sun}}\). Um desvio sistemático da linearidade pode indicar que a coleta, a recombinação ou o transporte de portadores dependem da intensidade. No entanto, a interpretação também deve considerar a forma como o código altera a escala da iluminação ótica e a convergência numérica de cada ponto.

As duas curvas de variação da intensidade respondem a questões complementares: a curva Suns–Voc analisa a tensão estabelecida sem corrente externa, enquanto a curva Suns–Jsc analisa a corrente coletada com uma tensão aplicada próxima de zero. Ambas devem ser interpretadas em conjunto com a curva J–V padrão, obtida sob uma única fonte de luz, e não como substitutas desta.

11. Modelagem óptica e distribuição espacial da fotogeração

O solucionador elétrico requer uma taxa de fotogeração espacial, mas essa taxa é fornecida pelo cálculo ótico. No exemplo atual, a fonte de luz está configurada para o espectro solar AM1.5G, incidindo pela parte superior (lado y0), e o modelo de geração é uma matriz de transferência. A malha espectral abrange de 320 a 1100 nm, com 201 pontos de comprimento de onda configurados; uma malha espacial ótica separada abrange todo o dispositivo de 3,4 µm, com 200 pontos configurados. Estas são discretizações óticas e são distintas da malha elétrica descrita anteriormente.

Abra a interface ótica e selecione os separadores de saída disponíveis (??–??). As capturas de ecrã ilustram as diferentes vistas de gráficos óticos, incluindo gráficos que podem apresentar poucos dados ou estar vazios para determinadas quantidades selecionadas. Um painel em branco não deve ser interpretado automaticamente como absorção ótica nula: primeiro, verifique se o cálculo ótico relevante foi executado, se existe um conjunto de dados de saída adequado e se os seus intervalos de comprimento de onda ou espaciais estão abrangidos pelos eixos selecionados.

O método da matriz de transferência avalia como as ondas óticas que se propagam para a frente e para trás interagem com uma pilha de camadas que possuem índices de refração complexos \(\tilde n(\lambda)=n(\lambda)+ik(\lambda)\). A solução do campo pode considerar a reflexão e a interferência; a absorção local é então convertida numa taxa de geração, utilizando a energia dos fotões. No caso do absorvedor de GaAs, o coeficiente de absorção dependente do comprimento de onda é um parâmetro essencial, pois determina onde os fotões de diferentes energias são absorvidos no interior da camada de 2,8 µm.

Preste especial atenção ao interpretar os dados óticos de AlGaAs. O modelo de material elétrico utiliza uma fração de liga especificada explicitamente, enquanto a base de dados de materiais óticos selecionada é referenciada de forma independente. Antes de utilizar o cálculo para o projeto espectral quantitativo, verifique se o conjunto de dados Inorganic/AlGaAs/Adachi fornece, de facto, a dependência de composição adequada nos comprimentos de onda relevantes. Uma fração de liga elétrica de 0,30 não comprova que todas as constantes óticas variem corretamente com a composição.

De forma semelhante, a estrutura 1D simulada idealiza a geometria ótica. Não modela de forma independente uma grelha metálica frontal fabricada, a dispersão lateral da corrente, a textura da superfície, os revestimentos antirreflexo detalhados ou os defeitos introduzidos por um processo epitaxial real, a menos que essas características sejam especificamente incluídas nos modelos óticos e elétricos. A análise ótica é, portanto, especialmente útil para compreender as tendências e identificar as possíveis perdas espectrais, enquanto a comparação com células medidas requer uma estrutura cuidadosamente ajustada e dados de materiais calibrados.

12. Premissas, verificações e limitações

Este exemplo deve ser encarado como uma demonstração numérica prática de uma heterojunção AlGaAs/GaAs, e não como uma previsão validada da eficiência certificada de uma célula solar específica. A estrutura de quatro camadas, os contactos e a malha provêm da configuração de simulação fornecida, e as capturas de ecrã mostram que vários cálculos e visualizações de resultados foram preparados. Uma previsão quantitativa fiável exige, adicionalmente, materiais, ótica, coeficientes de recombinação, qualidade dos contactos e convergência numérica verificados de forma independente.

Geometria e dimensionalidade. O visualizador 3D é uma representação visual; a dependência da profundidade elétrica e ótica neste tutorial é unidimensional. Efeitos como sombreamento da grade frontal, perdas nas bordas, queda de tensão lateral e contactos não uniformes estão fora desta descrição básica unidimensional.

Consistência entre o material e o modelo. A largura da banda proibida e a afinidade eletrónica numa liga cuja composição varia devem ser obtidas a partir de modelos consistentes. Os scripts Lua para cálculos elétricos e os conjuntos de dados óticos \(n,k\) são entradas distintas. Qualquer alteração aplicada apenas a um dos lados pode criar, inadvertidamente, descontinuidades artificiais, especialmente na heterointerface.

Recombinação e interfaces. Estão presentes as configurações de recombinação em massa do modelo SRH, mas isso não comprova que uma interface real de AlGaAs/GaAs tenha sido representada com uma velocidade de recombinação na interface, validada experimentalmente. A qualidade real da interface pode influenciar a tensão medida e a coleta de portadores. As alterações nas armadilhas e nos processos radiativos/de Auger devem ser testadas uma a uma e analisadas através dos perfis de taxa.

Reprodutibilidade numérica. Verifique a convergência dos resultados para toda a curva J–V e para a variação da intensidade, realize um refinamento da malha próximo à heterointerface e teste se as principais métricas fotovoltaicas e os perfis internos são estáveis. A falha de um ponto específico de alta intensidade não é evidência de um efeito físico no dispositivo. Utilize uma linha de base inalterada para distinguir as tendências físicas da sensibilidade do solucionador.

Analisando as capturas de tela. Cada imagem é um estado selecionado a partir dos resultados salvos da simulação, e os números das figuras neste tutorial identificam o índice da captura de tela, e não uma tensão aplicada garantida. Consulte a anotação da tensão real no gráfico do OghmaNano antes de atribuir um valor numérico de tensão a um perfil específico. Evite inferir um valor numérico para \(V_{OC}\), \(J_{SC}\), fator de preenchimento ou eficiência a partir de uma pequena miniatura da imagem, quando o arquivo de dados de saída fornece os valores exatos.

13. Exames complementares sugeridos