Início Exemplos Capturas de ecrã Manual do utilizador Logótipo Bluesky YouTube ☰
OghmaNano Platforma symulacji wielofizycznych dla urządzeń optoelektronicznych i systemów fotonicznych DESCARREGAR Quick Start guide

Estrutura de bandas de perovskita (8 bandas k·p) — CsPbX3 cúbico, MAPbX3 e ligas de haletos mistos.

1. Introdução

Calculated MAPbI3 perovskite electronic band structure using an eight-band k-p Hamiltonian
Estrutura eletrónica calculada de MAPbI3 (MAPI) utilizando o Hamiltoniano de oito bandas k·p, que mostra a banda de valência e as bandas de condução com separação devido ao efeito spin-órbita.

Calcular a estrutura de bandas e a mobilidade dos portadores de carga das perovskites de haletos metálicos é importante para compreender o seu desempenho em células solares, LEDs e lasers. Materiais como CsPbI3, CsPbBr3 e MAPbI3 (MA = metilamónio) apresentam estruturas eletrónicas invulgares, ilustradas em ??. Ao contrário dos semicondutores convencionais do tipo III–V, como o GaAs, a sua banda de valência é predominantemente do tipo s e as suas bandas de condução são do tipo p. O forte acoplamento spin-órbita divide as bandas de condução em aproximadamente 1,4–1,5 eV, sendo que a banda de condução mais baixa resultante forma a banda proibida fundamental no ponto R da zona de Brillouin na fase cúbica.

Estas bandas eletrónicas determinam propriedades importantes, incluindo a massa efetiva, as estatísticas dos portadores de carga, o ganho ótico e a mobilidade. Embora a teoria funcional da densidade (DFT) baseada em princípios primeiros possa calcular a estrutura eletrónica, a avaliação de redes densas de vetores de onda para cálculos de transporte pode ser dispendiosa. O método k·p de oito bandas oferece uma alternativa eficiente, utilizando um Hamiltoniano compacto parametrizado com base na DFT. Aqui, utilizamos o modelo de Gawarecki et al., Physical Review Applied 22, 014058 (2024) (DOI; arXiv), que descreve seis perovskitas de haleto de chumbo cúbicas e pseudocúbicas.

Este tutorial explica como calcular a estrutura de bandas eletrónicas de CsPbCl3, CsPbBr3, CsPbI3, MAPbCl3, MAPbBr3, MAPbI3 e ligas de halo-elementos mistos, utilizando o OghmaNano. Começando com os parâmetros k·p publicados, calculamos a dispersão das bandas e, em seguida, utilizamos os estados eletrónicos resultantes para avaliar o espalhamento eletrão-fónon de Fröhlich, os tempos de relaxamento e a mobilidade dos portadores limitada pelos fónons. Também analisamos as estatísticas dos portadores e o nível de Fermi em função da densidade de portadores.

2. Criar uma nova simulação

Para iniciar o cálculo da estrutura de bandas e da mobilidade dos portadores de carga de uma perovskita, clique em «Nova simulação» na barra de ferramentas do OghmaNano para abrir a caixa de diálogo de tipo de simulação (veja ??). Clique duas vezes em «Espalhamento em massa» para exibir os oito exemplos de perovskitas disponíveis no modelo k·p (veja ??).

São fornecidos exemplos para CsPbCl3, CsPbBr3, CsPbI3, MAPbCl3, MAPbBr3, MAPbI3 e várias ligas de haletos mistos. Qualquer um destes pode ser usado com este tutorial, mas começaremos com o CsPbCl3. Clique duas vezes no exemplo correspondente de dispersão de 8 bandas k·p, escolha um diretório e salve a nova simulação.

OghmaNano new simulation dialog showing the Bulk scattering category
Nova janela de diálogo de simulação. Clique duas vezes em «Espalhamento em massa».
List of eight-band k-p perovskite scattering examples for CsPbX3 and MAPbX3
Exemplos de espalhamento de 8 bandas disponíveis. Selecione CsPbCl3 para este tutorial.
Main OghmaNano simulation window showing a single CsPbCl3 bulk material layer
Janela principal do OghmaNano após a criação da simulação de CsPbCl3.

Após ser guardado, o OghmaNano abre a janela principal de simulação (ver ??), exibindo uma única camada ativa de CsPbCl3. Uma vez que se trata de um cálculo de um material em massa, a camada é simplesmente utilizada para definir o material que está a ser estudado; a sua espessura física não afeta a estrutura de bandas em massa calculada. A simulação está agora pronta para que possamos analisar os parâmetros do material e o modelo de oito bandas k·p.

3. Análise da simulação e dos parâmetros dos materiais.

Antes de executar o cálculo da estrutura de bandas, vale a pena dedicar um momento para compreender como a simulação está configurada. Abra o editor de camadas (veja ??). Verá uma única camada de 40 nm de CsPbCl3, com o material ótico definido como Generic/n/2.0. Isto é apenas um espaço reservado: uma vez que estamos a calcular as propriedades eletrónicas de um material em massa, nem a espessura da camada nem as suas constantes óticas são utilizadas nos cálculos da estrutura de bandas ou da mobilidade. Portanto, não é necessário alterar nada no editor de camadas.

Em vez disso, abra os parâmetros elétricos na guia «Estrutura do dispositivo» (veja ??). Repará que os parâmetros de material habituais, incluindo a mobilidade, a largura de banda, a densidade efetiva de estados e a permissividade, estão desativados e exibem a mensagem «Ver script para o valor». Isto é intencional. Em vez de inserir valores individuais através da interface gráfica, esta simulação utiliza um script Lua para definir as propriedades físicas do material.

Clique na aba «Editor de código» para visualizar o script do material, CsPbCl3.lua (consulte ??). O início do ficheiro contém informações sobre o material, a sua estrutura cristalina, as fontes dos parâmetros e as convenções de unidades. As propriedades do material são definidas mais adiante no ficheiro e, quando apropriado, podem depender da temperatura ou da composição.

Desloque-se até ao final do código (veja ??). Aqui encontrará os parâmetros numéricos utilizados para calcular a estrutura de bandas eletrónicas, incluindo a diferença de energia entre as bandas, o efeito de acoplamento spin-órbita e vários parâmetros que controlam a forma das bandas. Cada valor é acompanhado por uma referência bibliográfica, o que permite rastrear o cálculo até à investigação original.

OghmaNano layer editor showing the single CsPbCl3 bulk material layer
Editor de camadas. A espessura e o material ótico não são relevantes para este cálculo global.
Electrical parameter editor showing greyed-out fields with See script for value
Parâmetros elétricos. Os valores são fornecidos pelo script do material, e não pelo editor gráfico.
Lua code editor showing the CsPbCl3 material file header and literature references
Editor de código que exibe o cabeçalho de CsPbCl3.lua e as suas fontes de parâmetros.
Band structure parameters in the CsPbCl3 material script showing gamma3 prime and its reference
Os parâmetros utilizados para calcular a estrutura de bandas podem ser encontrados no final do ficheiro do material.

Compreendendo os parâmetros utilizados para calcular a estrutura de bandas.

O cálculo da estrutura de bandas utiliza um conjunto de parâmetros de materiais retirados de Gawarecki et al. (2024). Estes descrevem a energia e a forma das bandas eletrónicas próximas à banda proibida. Os valores para as seis perovskites suportadas estão apresentados em ??. A parte superior da tabela contém os parâmetros utilizados diretamente no cálculo, enquanto a parte inferior apresenta valores adicionais derivados destes.

Parâmetro CsPbCl3 CsPbBr3 CsPbI3 MAPbCl3 MAPbBr3 MAPbI3
Eg (eV) 2.7442.0731.416 3.0072.2771.549
Δc (eV) 1.4441.4761.494 1.5061.5081.500
δ (eV) 000 0.0450.0570.075
ζ (eV) 000 0.0160.0300.045
P paralelo (eV·Å) 9.2338.9488.601 8.8788.6238.097
Pz (eV·Å) 9.2338.9488.601 9.8969.3029.058
γ′1 1.6432.1832.997 1.5812.1442.967
γ′2 0.1900.3940.683 0.1400.1900.286
γ′3 −0.691−1.081−1.169 −0.599−1.060−1.150
Parâmetros numéricos utilizados para construir o Hamiltoniano de oito bandas k·p para o composto cúbico CsPbX.3 e pseudo-cúbico MAPbX3 perovskitas. Os valores numéricos foram retirados da Tabela I, que está disponível gratuitamente. Pré-publicação no arXiv de Gawarecki. e outros. (arXiv:2306.08643).

O que significam os parâmetros?

Os principais parâmetros utilizados no cálculo da estrutura de bandas são explicados abaixo. Não é necessário alterar estes valores para seguir o tutorial, mas compreender o que representam será útil ao comparar diferentes materiais perovskitas.

Todos estes parâmetros já estão incluídos nos ficheiros de materiais do OghmaNano, pelo que não é necessário introduzi-los manualmente. Quando é selecionado um material diferente, os respetivos parâmetros são carregados automaticamente. Na secção seguinte, veremos como alternar entre perovskitas e explorar diferentes composições de materiais.

4. Alteração dos parâmetros dos materiais

Nestes cálculos da estrutura de bandas em grande escala, o material é definido integralmente pelo seu script Lua. A simulação é, na prática, um cálculo elétrico de dimensão zero (0D), o que significa que as dimensões físicas apresentadas no editor de camadas não são utilizadas. Em vez disso, todas as propriedades do material necessárias para calcular a estrutura de bandas, a dispersão e a mobilidade são obtidas diretamente do script.

Existem duas formas de alterar os parâmetros do material. Pode editar o script existente diretamente, por exemplo, se encontrar valores alternativos na literatura científica ou se quiser investigar como a alteração de um determinado parâmetro afeta os resultados. Em alternativa, pode carregar um perovskite completamente diferente a partir da base de dados de materiais da OghmaNano.

Para importar um novo material, clique na pequena seta ao lado do botão «Script» no editor de parâmetros elétricos e selecione «Importar script de material» (veja ??). Isso abrirá a biblioteca de materiais em materials/Perovskites/ (veja ??). A biblioteca contém definições de materiais para CsPbCl3, CsPbBr3, CsPbI3, MAPbCl3, MAPbBr3, MAPbI3 e perovskitas de haletos mistos.

Electrical parameter editor showing the Script menu and Import material script option
Abra o menu suspenso «Script» e selecione «Importar script de materiais».
OghmaNano perovskite materials database showing available Lua material definitions
Base de dados de materiais perovskitas. Clique duas vezes num ficheiro de material para importar os seus parâmetros.

Ao clicar duas vezes num script de material, todos os seus parâmetros são importados para a simulação em curso. O script importado é copiado para a pasta da simulação, onde pode ser visualizado ou modificado através do editor de código. Por exemplo, ao importar CsPbBr3.lua, os parâmetros originais de CsPbCl3 são substituídos pelos de CsPbBr3. Os novos valores serão utilizados na próxima vez que a simulação for executada.

Se editar manualmente algum dos valores, é recomendável atualizar as referências bibliográficas correspondentes, de modo que a origem de cada parâmetro permaneça clara. Lembre-se de que, para estes cálculos em lote, a alteração do script Lua é o que modifica o material; não é necessário modificar o editor de camadas.

5. Executar a simulação e calcular a estrutura de bandas.

Estamos agora prontos para calcular a estrutura de bandas eletrónicas do nosso perovskite. Clique em «Executar simulação» na janela principal do OghmaNano. O separador «Terminal» apresenta os parâmetros do material que estão a ser utilizados, seguidos do progresso dos cálculos da estrutura de bandas e da dispersão (ver ??). É importante verificar estes valores, especialmente se tiver alterado o script do material ou selecionado uma composição mista de haletos.

Neste exemplo, estamos a calcular as propriedades de CsPb(I0,5Br0,25Cl0,25)3. Os parâmetros apresentados no Terminal podem ser comparados com os valores em ?? utilizando o método de interpolação linear descrito na secção anterior. Por exemplo, o efeito de divisão spin-órbita é calculado da seguinte forma:

\[ \Delta_c = 0.5(1.494) + 0.25(1.476) + 0.25(1.444) = 1.477\ \mathrm{eV}. \]

Este é o valor exibido pelo programa de resolução. Os restantes parâmetros da estrutura de bandas são calculados da mesma forma, o que permite verificar se a composição da liga selecionada foi carregada corretamente.

OghmaNano terminal displaying perovskite band structure parameters and progress of the Frohlich scattering calculation
Saída do terminal que mostra os parâmetros dos materiais utilizados no cálculo, seguida do progresso do cálculo do espalhamento de Fröhlich.
OghmaNano output browser showing the calculated band structure, mobility and scattering results
O separador «Resultados» apresenta os ficheiros que contêm a estrutura de bandas calculada, juntamente com os diretórios que armazenam os resultados relativos à mobilidade e à dispersão.
Calculated band structure of mixed-halide CsPbI0.5Br0.25Cl0.25 perovskite showing valence and conduction band energies against wavevector
Estrutura de bandas eletrónicas calculada para o composto CsPb(I0,5Br0,25Cl0,25)3. Vermelho: banda de valência; ciano: banda de condução mais baixa; curvas superiores: bandas de condução mais altas. O vetor de onda é medido em relação ao ponto R.

Após a conclusão da simulação, abra a guia «Resultados» (consulte ??). O OghmaNano guarda os resultados calculados como ficheiros CSV, que podem ser visualizados diretamente no software ou exportados para análise posterior. As bandas eletrónicas individuais são guardadas no diretório bands, enquanto três ficheiros adicionais oferecem formas práticas de visualizar os resultados:

Clique duas vezes em Ec_Ev.csv para exibir a estrutura de bandas de perovskita calculada, conforme mostrado em ??. Este gráfico mostra como a energia dos estados eletrónicos disponíveis varia com o vetor de onda (k). Em termos simples, mostra a forma das bandas eletrónicas, a partir da qual podem ser determinadas propriedades importantes, como a largura da banda e as massas efetivas.

Para compreender os resultados, comece por analisar o centro do gráfico, onde k = 0 corresponde ao ponto R da zona de Brillouin. O ponto mais alto da banda de valência é utilizado como referência de energia e, por isso, é definido como 0 eV. Ao analisar o gráfico de baixo para cima, podem ser identificadas três características principais:

Uma característica importante da estrutura de bandas é a curvatura das bandas. As bandas com maior curvatura correspondem geralmente a portadores de carga com massas efetivas menores, enquanto as bandas mais planas correspondem a portadores mais pesados. Isso torna-se particularmente importante ao calcular a mobilidade dos portadores de carga, pois a estrutura de bandas determina como os eletrões e as lacunas respondem a um campo elétrico aplicado.

Por fim, é importante lembrar que este modelo descreve a estrutura eletrónica próxima à banda proibida, e não em toda a zona de Brillouin. O intervalo do vetor de onda apresentado aqui estende-se até aproximadamente 109 m−1 (0,1 Å−1), valor comparável ao intervalo no qual Gawarecki e colaboradores demonstram uma boa concordância com os cálculos de DFT (a sua Figura 3). As previsões que se afastam muito do ponto R devem, portanto, ser consideradas com cautela.

6. Cálculo da estrutura de bandas de perovskitas mistas de haletos

OghmaNano Electrical ribbon showing the Gradients and Strain button used to define mixed-halide perovskite compositions
A fita elétrica. Clique em «Gradientes/Deformação» para definir a composição mista de haletos.

As perovskitas de haletos mistos oferecem uma forma prática de ajustar as propriedades eletrónicas de um material, alterando as proporções relativas de iodo, bromo e cloro. Por exemplo, substituir parte do iodo em CsPbI3 por bromo geralmente aumenta a largura da banda proibida, permitindo que as propriedades óticas e eletrónicas sejam ajustadas para diferentes aplicações. Na prática, os investigadores geralmente trabalham com composições que se situam entre os três compostos de haletos puros.

O OghmaNano permite estimar a estrutura de bandas de composições mistas de haletos arbitrárias, sem exigir um conjunto separado de parâmetros de estrutura de bandas publicados para cada mistura possível. Isto é particularmente útil ao estudar materiais preparados em laboratório: ao inserir a composição de uma amostra, é possível calcular a sua estrutura de bandas prevista, analisar como a diferença de energia entre as bandas varia com a composição e, posteriormente, explorar a sua dispersão e a mobilidade dos portadores de carga. Também permite investigar como a substituição gradual de um haleto por outro afeta as propriedades eletrónicas do material.

São fornecidas duas formulações de materiais com haletos mistos: CsPb(I1−x−yBrxCly)3 e a sua versão equivalente com metilamónio (MA). Aqui, x representa a fração de brometo, y representa a fração de cloreto e o valor restante, 1 − x − y, representa a fração de iodeto. As frações devem satisfazer as condições x ≥ 0, y ≥ 0 e x + y ≤ 1, de modo que a composição total de haletos seja sempre igual a 100%.

Para calcular as propriedades de uma liga, primeiro importe o script do material misto de haletos apropriado, utilizando o método descrito na seção anterior. Em seguida, abra a faixa de ferramentas «Elétrica» e clique em «Gradientes/Deformação» (veja ??). Isso abrirá o editor de composição do material (veja ??), onde as frações de brometo e cloreto podem ser inseridas.

O editor permite especificar a composição no início (x0, y0, z0) e no fim (x1, y1, z1) da camada. Uma vez que estamos a realizar um cálculo em lote, a composição deve ser uniforme, pelo que deve inserir os mesmos valores em ambos os extremos. Por exemplo, ao definir x0 = x1 = 0,25 e y0 = y1 = 0,25, obtemos CsPb(I0,5Br0,25Cl0,25)3, que contém 50% de iodeto, 25% de brometo e 25% de cloreto. Isto é ilustrado em ??. Os valores de z não são utilizados nestes scripts de materiais.

Mixed-halide perovskite composition editor with bromide and chloride fractions set to 0.25
Editor de composição que mostra x = y = 0,25 para CsPb(I0,5Br0,25Cl0,25)3.
Mixed-halide perovskite material script showing the composition rules and parameter interpolation method
Cabeçalho do script do material de haletos mistos, que apresenta as regras de composição e o método utilizado para calcular os parâmetros da liga.
Lua material script showing linear interpolation of relative permittivity between the pure iodide bromide and chloride materials
Exemplo de interpolação linear no script do material, mostrando como a permissividade relativa é calculada.

Para estimar as propriedades eletrónicas de uma liga, a OghmaNano utiliza a interpolação linear entre os parâmetros dos três compostos puros. Estes são os mesmos parâmetros dos materiais que foram publicados em ??. Em vez de interpolar a estrutura de bandas final, o software calcula primeiro um novo conjunto de parâmetros de materiais para a composição escolhida e, em seguida, utiliza esses parâmetros para realizar o cálculo da estrutura de bandas.

Por exemplo, a largura da banda proibida, o efeito de divisão spin-órbita e os parâmetros que controlam a forma das bandas eletrónicas são calculados, em cada caso, como uma média ponderada dos valores do iodeto, brometo e cloreto. A expressão geral é:

\[ p(x,y) = (1-x-y)\,p_{\mathrm{I}} + x\,p_{\mathrm{Br}} + y\,p_{\mathrm{Cl}}, \]

onde p representa um parâmetro específico do material, e pI, pBr e pCl são os valores correspondentes para os três compostos puros. No exemplo acima, o cálculo utiliza 50% do valor do iodeto, 25% do valor do brometo e 25% do valor do cloreto.

Como exemplo, utilizando os valores das diferenças de energia da tabela de parâmetros, a diferença de energia estimada para o composto CsPb(I0.5Br0.25Cl0.25)3 é a seguinte:

\[ E_g = 0.5(1.416) + 0.25(2.073) + 0.25(2.744) = 1.912\ \mathrm{eV}. \]

A mesma abordagem de interpolação é aplicada às outras propriedades numéricas no script do material, incluindo a permissividade relativa mostrada em ??. Quando uma propriedade depende da temperatura, os seus coeficientes correspondentes são interpolados. Isto permite gerar um conjunto consistente de parâmetros para qualquer composição mista de haletos escolhida e, em seguida, utilizá-los para calcular a estrutura de bandas eletrónicas e as propriedades de transporte subsequentes.

A interpolação linear é um ponto de partida útil para explorar a composição de materiais, especialmente quando os dados experimentais são limitados. No entanto, as perovskites mistas de haletos reais nem sempre apresentam um comportamento linear. Por exemplo, a largura da banda proibida pode apresentar uma curvatura, em que o valor medido difere do valor previsto pela interpolação simples. O modelo também não considera explicitamente a desordem local da liga e mantém as premissas da fase cristalina dos parâmetros originais. Para comparações quantitativas com uma amostra específica de laboratório, as propriedades calculadas devem, portanto, ser verificadas em relação às medições experimentais, ajustando o script Lua, quando necessário.

6. Cálculo da mobilidade dos portadores de carga e da dispersão em perovskitas.

Após calcular a estrutura de bandas eletrónicas, a OghmaNano pode agora utilizar estes resultados para determinar a facilidade com que os eletrões e as lacunas se movem através da perovskite. Para tal, é necessário considerar não só a forma das bandas eletrónicas, mas também a forma como os portadores de carga interagem com as vibrações da rede cristalina, conhecidas como fonões. Estas interações dispersam os portadores, limitando a sua mobilidade.

Nos perovskites de haletos de chumbo, uma importante fonte de dispersão é a interação entre os portadores de carga e os fonões óticos longitudinais (FOL), conhecida como dispersão de Fröhlich entre eletrões e fonões. A intensidade desta interação depende da energia dos fonões e das propriedades dielétricas do material. A OghmaNano utiliza as bandas eletrónicas calculadas, juntamente com os parâmetros dos fonões definidos no script Lua do material, para avaliar estes processos de dispersão e estimar a mobilidade resultante dos portadores de carga. Uma explicação mais detalhada da física subjacente pode ser encontrada no nosso guia sobre a dispersão de Fröhlich entre eletrões e fonões.

Os parâmetros dos fónons são obtidos a partir de fontes bibliográficas distintas, citadas nos documentos do material; não fazem parte do conjunto original de parâmetros da estrutura de bandas de Gawarecki et al. As mobilidades calculadas aqui são, portanto, mobilidades limitadas por Fröhlich, descrevendo o transporte de carga quando este mecanismo de dispersão é considerado.

Após a conclusão da simulação, abra a guia «Resultados» (consulte ??). Além dos ficheiros da estrutura de bandas mencionados anteriormente, encontrará duas pastas importantes: scattering e mobility. Estas contêm as taxas de dispersão calculadas, os tempos de relaxamento, as estatísticas dos portadores de carga e os resultados da mobilidade.

Ao abrir o diretório scattering (ver ??), é possível observar uma série de ficheiros CSV que descrevem o espalhamento de Fröhlich para as diferentes bandas eletrónicas. Os ficheiros que contêm frohlich_rate indicam a taxa de espalhamento, enquanto os que contêm frohlich_tau fornecem o tempo de relaxamento correspondente. Uma taxa de espalhamento mais elevada significa, geralmente, que os portadores de carga sofrem o fenómeno de espalhamento com maior frequência, o que resulta num tempo de relaxamento mais curto.

Em seguida, abra o diretório mobility (consulte ??). Este contém as mobilidades calculadas de eletrões e lacunas, as taxas de espalhamento efetivas, os tempos de relaxamento e os níveis de Fermi, expressos como funções da densidade de portadores de carga. Estes ficheiros permitem-nos investigar como as propriedades de transporte se alteram à medida que são introduzidos mais portadores de carga no material.

OghmaNano output tab showing the calculated band structure, scattering and mobility results
Resultado da simulação que mostra os diretórios scattering e mobility.
OghmaNano scattering output directory containing Fröhlich scattering rates and relaxation times
O diretório de dispersão contém as taxas de dispersão de Fröhlich calculadas e os tempos de relaxamento para cada banda individual.
OghmaNano mobility output directory containing electron and hole mobility, scattering and Fermi level results
O diretório de mobilidade contém as propriedades calculadas de transporte de eletrões e lacunas, expressas em função da densidade de portadores de carga.

Para analisar os resultados, clique duas vezes em kp8_perovskite_electron_fermi_level.csv, kp8_perovskite_electron_mobility.csv e kp8_perovskite_hole_mobility.csv. Os gráficos resultantes são apresentados abaixo para CsPb(I0.5Br0.25Cl0.25)3.

Calculated electron Fermi level as a function of carrier density for a mixed-halide perovskite
Nível de Fermi dos eletrões calculado em função da densidade de portadores de carga.
Calculated Fröhlich-limited electron mobility versus carrier density for a mixed-halide perovskite
Mobilidade eletrónica limitada pelo modelo de Fröhlich, calculada em função da densidade de portadores de carga.
Calculated Fröhlich-limited hole mobility versus carrier density for a mixed-halide perovskite
Mobilidade dos portadores de carga, limitada pelo modelo de Fröhlich, calculada em função da densidade de portadores.

O primeiro gráfico (ver ??) mostra como o nível de Fermi dos eletrões varia com a densidade de portadores de carga. O nível de Fermi descreve a forma como os eletrões se distribuem pelos estados eletrónicos disponíveis. À medida que mais eletrões são introduzidos no material, os estados de menor energia tornam-se cada vez mais ocupados e o nível de Fermi desloca-se para cima. Para esta liga, a banda de condução mais baixa começa aproximadamente a 1,91 eV. À medida que o nível de Fermi se aproxima e entra nesta banda, a distribuição dos eletrões altera-se, o que também pode afetar as propriedades de transporte calculadas.

Os dois gráficos restantes mostram as mobilidades correspondentes dos eletrões e das lacunas. É interessante notar que a mobilidade não se mantém constante à medida que a densidade de portadores aumenta. Neste exemplo, a mobilidade dos eletrões diminui ligeiramente no início, antes de aumentar em densidades mais elevadas (ver ??). A mobilidade das lacunas calculada aumenta na maior parte do mesmo intervalo de densidade (ver ??).

A mobilidade calculada pode aumentar ou diminuir com a densidade de portadores de carga, dependendo do material. Este comportamento resulta de dois fatores principais:

No caso da perovskite de haletos mistos considerada aqui, a mobilidade eletrónica calculada diminui ligeiramente no início, antes de aumentar, enquanto a mobilidade das lacunas aumenta na maior parte da faixa de densidade representada no gráfico. Outras composições de perovskite podem apresentar tendências diferentes, uma vez que as suas estruturas de banda e taxas de dispersão são diferentes. Estes cálculos permitem-nos, portanto, analisar como a alteração da composição do material afeta o transporte de portadores de carga.

7. Âmbito e Limitações do Modelo

A abordagem apresentada neste tutorial oferece um cálculo baseado na física da estrutura de bandas eletrónicas e da mobilidade dos portadores de carga em perovskites de haleto metálico tridimensionais. O espalhamento eletrão-fónon de Fröhlich é amplamente reconhecido como um dos principais mecanismos intrínsecos que limitam a mobilidade dos portadores de carga nestes materiais, especialmente à temperatura ambiente. A inclusão desta interação oferece, portanto, uma forma fisicamente relevante de investigar o transporte de carga e comparar diferentes composições de perovskites.

O modelo combina a estrutura de bandas eletrónicas calculada com a dispersão de fotões óticos polares para determinar as taxas de dispersão, os tempos de relaxamento e a mobilidade dos portadores de carga, em função da energia. As principais premissas são:

Ao combinar uma estrutura eletrónica de bandas detalhada com um dos principais mecanismos de dispersão intrínsecos em perovskites de haleto metálico, a OghmaNano oferece uma abordagem útil e computacionalmente eficiente para calcular a mobilidade limitada por fónons. É particularmente adequada para explorar tendências de materiais, comparar composições de haletos e compreender como as alterações na estrutura eletrónica influenciam o transporte de carga.

8. Exploração Adicional

Agora que consegue calcular as estruturas de banda de perovskitas, as taxas de dispersão e a mobilidade dos portadores de carga, experimente alterar os parâmetros do material para verificar como diferentes composições influenciam os resultados. Seguem-se três pontos de partida úteis:

Estes exemplos servem como ponto de partida. Também pode editar parâmetros individuais nos scripts de materiais Lua para analisar materiais descritos na literatura ou composições preparadas no seu próprio laboratório.