Algoritmo para simulação de zonas proibidas em estruturas eletrônicas de bandas

Algorithm for simulation of forbidden zones in electronic band structure

ARTIGOS GERAIS

Algoritmo para simulação de zonas proibidas em estruturas eletrônicas de bandas

Algorithm for simulation of forbidden zones in electronic band structure

Wagner Leite Araujo¹

Departamento de Engenharia Nuclear, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brasil

RESUMO

Este artigo visa apresentar um algoritmo descritivo de gaps energéticos em semicondutores utilizando o modelo de Kronig-Penney. O código desenvolvido proporciona uma ferramenta didática na descrição de zonas proibidas em redes cristalinas, tema comumente tratado em cursos de estrutura da materia e física de estado sólido. Os cálculos presentes no algoritmo seguem detalhados. O principal objetivo e descrever uma metodologia para se obter curvas da energia vs. número de onda, tratando o coeficiente de transmissão do material como constante, caso bem conhecido na literatura, e como função da energia.

Palavras-chave: modelo de Kronig-Penney, funções de bloch, gaps de energia, estruturas eletrônicas em um cristal.

ABSTRACT

This paper presents a descriptive algorithm of energy gaps in semiconductors applying the Kronig-Penney model. The code provides a teaching tool in the description of forbidden zones in crystalline lattices, which is an issue usually treated in structure of matter and solid state physic courses. All calculations contained in the algorithm is detailed. The main goal is to describe a methodology to obtain curves of energy vs. wave number, addressing the transmission coefficient of the material as constant, well known in the literature, and as a function of energy.

Keywords: Kronig-Penney model, bloch functions, energy gaps, electronic structures in a crystal.

1. Introdução

Elétrons de um único átomo ocupam orbitais que definem um ajuste discreto de energia. Quando átomos são colocados juntos em uma molecula, seus orbitais se separam como em uma oscilação acoplada. Isso produz um número de orbitais moleculares proporcional ao número de átomos [1]. Quando milhares de átomos interagem uns com os outros formando um sólido, o número de orbitais se torna excessivamente grande e a diferença de energia entre eles se torna muito pequena; assim, os níveis podem ser considerados, praticamente, bandas contínuas de energia e não níveis discretos como em átomos isolados. Contudo, alguns intervalos de energia não contém orbitais, não importando quantos átomos estejam agregados; isso da origem aos gaps de energia.

O objetivo desse artigo é descrever um algoritmo, baseado no modelo de Kronig-Penney, que conduz a curva da energia em função do número de onda para o coeficiente de transmissão constante, sendo essa, bem conhecida na literatura [2-5] e, também, dependente da energia, cujo esboço não foi encontrado em textos relacionados. Esse algoritmo pode ser facilmente modificado, possibilitando novos testes e diagnósticos. O modelo de Kronig-Penney serúa previamente abordado, após considerações sobre o teorema de Bloch [6]. O desígnio final e averiguar a presença das zonas proibidas. Os cálculos serão detalhados no Apêndice A. No Apêndice B, encontra-se um algoritmo projetado para o Scilab [7], que e um programa livre. Esse algoritmo pode ser facilmente modificado possibilitando novos testes e diagnósticos.

2. Considerações iniciais

O comportamento de elétrons transitando no interior de um cristal semicondutor pode ser comparado a partículas em uma caixa tridimensional cujo interior é demasiadamente complicado. Os átomos dessa estrutura vibram em posições específicas na rede. Em cristais reais existem defeitos estruturais que se originam, por exemplo, de átomos perdidos e impuros. Uma maneira de simplificar o modelamento matemático do problema é assumir que tais defeitos e as vibrações do núcleo do átomo conduzirão a uma perturbação de segunda ordem. De fato, essa hipótese expressa a estrutura da rede como perfeita e os átomos como fixados em posições bem determinadas.

O modelo de Kronig-Penney [8] será empregado para as análises tratadas nesse artigo. Este, admite soluções analíticas, permitindo simplificações algebricas. A estrutura de bandas eletrônicas está diretamente relacionada a várias propriedades macroscópicas de materiais e, por conseguinte, tem desempenhado um interesse vasto na área de semicondutores. Materiais hipoteticos são investigados por cálculos envolvendo tais estruturas [9-12]; quando esses materiais mostram propriedades atrativas, despertam interesse e, quando possível, são preparados em ensaios experimentais.

A teoria de Kronig-Penney é um arqétipo unidimensional da mecânica quântica na análise de redes cristalinas. Apesar das simplificações, a estrutura de bandas eletrônicas oriunda desse modelo proporciona o estudo de várias características referentes a composições de bandas que resultam de modelos mais complexos envolvendo cálculos numéricos extensos.

A análise reportada por este trabalho terá como base o estudo de uma partícula presente em uma estrutura unidimensional dentro de uma rede cristalina periódica. Este problema pode ser simplificado tratando uma barreira de potencial em três dimensões (uma partícula em uma caixa [13]) de forma unidimensional. O potencial que define tal barreira e causado pelos íons na estrutura periódica do cristal. Essa estrutura cria um campo elétromagnético de forma que os elétrons ficam sujeitos a um potencial regular. Isso e uma extensão do modelo do elétron livre [14], que assume um potencial nulo dentro da estrutura.

3. Formação de bandas no modelo de elétron livre

3.1. O teorema de Bloch

Se o potencial U (x) e uma função periódica tal que U (x + a) = U (x), o teorema de Bloch pode ser escrito como (x + a) = e ^ika(x), onde (x) representa a função de onda unidimensional, k o número de onda ²π /λ), α a distância interatômica, x a distância na rede cristalina e i = . De forma equivalente, o teorema ainda pode ser escrito como (x) = e^ikxu (x) , sendo que a função de onda da célula tem a mesma periodicidade do potêncial, isso siginifica que u (x + a) = u (x). Então, (x + a) = e^ik (x + a) u (x + a) = e^ika e^ikx u (x) = e ^ika (x). A função de onda (x) claramente não é periódica. Ela tem a forma de uma onda plana ( emodulada por uma função u (x) que reflete a periodicidade da célula cristalina e a periodicidade do potencial.

O uso de condições de contorno periódicas em um arranjo de N átomos implica que (x) = (x+ Na) = e^ikNa (x), então e^ikNa= 1 =e²ηπ assim k=²ηπ/Na, para η = 0, ±1, ±2,± ^N/2. Para um cristal com dimensões finitas, o número de onda k pode assumir apenas valores inteiros limitados por ± ^π/a. Em um cristal real, a grande quantidade de átomos implica em valores de k muito próximos e, assim, definem praticamente uma continuidade na energia.

Vale ressaltar que razao do número de onda k ser proporcional a^ηπ/α é devido a formação de ondas estacionárias para os estados nas extremidades das barreiras. As funções de onda interferem construtivamente, com diferença de fase de 2π com a função de onda do vizinho mais próximo como mostrado na Fig. 1. A onda ₁ (x) interfere construtivamente com ₂ (x). Então, pode-se dizer que λ= 2α. Com essa consideração, a interferencia será construtiva sempre que 2α = λ, 2λ, 3λ, ... = πη, o que implica que λ = ^2α/η. Porém λ = ^2π/k, logo, k = ηπ/α.

3.2. O modelo de Kronig-Penney

O modelo de Kronig-Penney é fundamentado em uma rede cristalina unidimensional e infinita, onde se utiliza uma solução analítica conhecida da equação de Schrödinger para uma barreira de potencial [15]. Esta equação, originada do trabalho de Erwin Schrödinger [16], descreve a evolução espacial e/ou temporal de um sistema quântico e, quando se remove a dependência temporal, é dada pela Eq. (1).

O modelo de Kronig-Penney se basea em resolver a Eq. (1) na célula ilustrada na Fig. 2, região A periodicidade da rede cristalina e por conseguinte, a periodicidade dos potenciais influenciam nas funções de onda. Tal comportamento é considerado por intermédio do teorema de Bloch, onde os potenciais periódicos reais são aproximados por potenciais periódicos em forma de barreira, como descrito na Fig. 2.

3.2.1 . Solução da equação de Schrödinger em uma barreira de potencial

Sejam _e (x) e _d (x) as soluções da Eq. (1), para uma onda propagante da esquerda para direita _e (x) (Fig. 3A) e da direita para esquerda _d (x) (Fig. 3B). Definem-se os números complexos r e t como o coeficiente de reflexão e de transmissão, respectiavamente. Na Fig. 3, i representa a onda incidente, r a porção refletida e t a parte transmitida da onda. Com isso, as soluções da Eq. (1) para os potenciais representados pelas Figs. 3A e 3B serão:

Através do princípio da superposição, a combinação linear (x) = α _e (x)+B_d (x) é solução da Eq. (1), já que _e (x) e _e(x) também são soluções.

O teorema de Bloch postula que

e que

Este teorema, em

, conduz a

Substituindo _e (x) e _d (x) das Eqs. (2) e (3) nas Eqs. (4) e (5), chega-se em um sistema de duas equações com duas varióveis que será resolvido para as incógnitas α e B (Eq. (6)). A solução deste sistema resulta na Eq. (7).

Ao considerar o potencial periódico igual a zero em todo o cristal, tem-se |t|² = 1 e |r|² =0. Ou seja, a transmitância será 1 e não haverá nenhuma reflexão. De fato, a Eq. (7) é consistente com essa consideração. No limite onde não há nenhum potencial cos(α) = 1/2 e^ika + 1/2 e^-ika. Pontua-se que a relação: cos (ax) =e^iax+e^-iax/2, permite escrever cos (α) = cos(Ka), isto é: = k. Assim, o perfil da energia será o mesmo do modelo de elétrons livres, E=; o que acontece quando se toma U_p (x + Na) = 0.

Desenvolvendo a Eq. (7), como detalhado no Apêndice 1, obtóem-se:

4. Bandas de energia e zonas de Brillouin

Através da Eq. (8), a razao matemática da formação linear de gaps de energia se torna clara. Nota-se que a função cos (α) só pode assumir valores no intervalo [-1,1], como |t| < 1, a função pode ser um número maior que 1. A energia em unidades de é dada pela equação E = k². As zonas proibidas (os gaps de energia) resultam, naturalmente, deste evento, como mostra a Fig. 4.

A Fig. 5 mostra uma representação das bandas na zona reduzida. E possível esboçar as curvas de energia em um diagrama extendido como mostrado na Fig. 6. Nesse diagrama a energia E pode crescer indefinidamente ao longo do eixo do número de onda k. O esquema da zona extendida possibilita uma comparação mais prática com a curva característica do modelo do elétron livre que e definida pela parábola E = .

A Fig. 6 permite definir uma terminologia relevante em física do estado sólido, a zona de Brillouin. Tais zonas são numeradas em função de intervalos de energia ou de número de onda. Um exemplo e mostrado na Fig. 4, onde as zonas de Brillouin são definidas. Na Fig. 6, a primeira zona de Brillouin e a região entre

5. O algoritmo

A construção do algoritmo se baseia na Eq. (8). Atraves desta, o algoritmo fornece a curva da energia vs. ka. A edificação da curva se dara em duas conjunturas:

• A primeira com o coeficiente de transmissão t = 0, 8.

• A segunda com a transmitancia variando com o k; para tal será escolhido um valor para a altura da barreira, U₀.

O algoritmo empregara um processo interativo para gerar pares de coordenadas que esboçarão a curva pretendida. O eixo da energia será dado em unidades , possibilitando certa arbitrariedade na escolha de unidades de medida. Por exemplo, ao se escolher 2 Å para α, que e um valor plausível, diretamente se infere nas unidades do eixo da energia; isso significa que esse acoplamento de unidades permite que se escolha α = 1 m, que apesar de absurdo, não traz prejuízo qualitativo na representação gráfica.

O perfil da curva se dara no esquema da zona reduzida e na região positiva do produto ka. Para generalizarmos toda a primeira zona de Brillouin, bastaria espelhar o resultado para valores negativos de ka.

Primeiramente e necessário definir linhas de comando para determinar os valores para a largura do poco de potencial a, coeficiente de transmissão t, fase δ e a acuidade do plot definida como dk (espacamento entre pares de coordenadas). Na linha seguinte, define-se o k₀ para representar o valor inicial de k. Em seguida e necessário gerar um comando que impõe um processo iterativo. Tal processo faz o i variar entre zero e o número de iterações determinada pelo usuário. Essa operação se inicia pela Eq. (9)

Na Eq. (10), f é definido para representar o lado esquerdo da Eq. (8)

e dado por uma função arco-cosseno. Logo, torna-se necessário uma linha de comando condicional que limita o resultado de f entre- 1 e 1, antes da Eq. (10)

Na Eq. (11), o k e isolado. A Eq. (12) fornece o valor da energia

Em seguida, gera-se um comando para determinar a energia em função de ka e/ou as coordenadas E e k. Temos assim a curva pretendida na primeira conjuntura.

A segunda conjuntura segue a mesma feição da primeira; contudo, deve-se definir a transmitôncia em relação à energia. Isso e realizado por meio da Eq. (13) para a energia maior que a barreira de potencial e pela Eq. (14) para a barreira maior que a energia. As Eqs. (13) e (14) são descritas por Eisberg e Resnick [15]

onde T representa o inverso do coeficiente de transmissão.

Como na primeira conjuntura, e necessário gerar uma linha solicitando entradas para os parâmetros a, dk, k_o, U_o, δ. O mesmo comando condicional aplicado na primeira conjuntura atua na Eq. (9), depois na Eq. (12) e em seguida um novo comando condicional e requerido. Este comando seleciona as energias que serão aplicadas na Eq. (13) ou na Eq. (14), isto é , se E > U₀, a energia será aplicada na Eq. (13); se E < U₀, aplica-se E na Eq. (14). Deve ser garantido que os denominadores das Eqs. (13) e (14) não assumam valores nulos. Para resolver esse empecilho, basta coloca uma linha de comando que selecione valores diferentes de zero para o denominador dessas equações. A proxima linha de comando faz com que f receba o valor da transmitância, como na Eq. (10). Em seguida, calcula-se por meio da Eq. (11) e se determina o comando para gerar a curva e/ou as coordenadas provedoras da curva.

6. Simulações e resultados

A Fig. 7 mostra a curva obtida no caso da transmitância constante. Os dados de entrada foram α = 1, dk = 0, 01, k_o = 0, δ = 0, 2, t = 0, 8 e i variando de 0 a 500. A Fig. 8 mostra a simulação com i variando de 0 a 1500. Os gaps de energia ficaram bem evidentes.

Na primeira simulação para a transmitancia variavel, a curva apresentou o perfil ilustrado na Fig. 9. Nesta, observa-se dois gaps de energia com um comportamento peculiar.

A Fig. 10 apresenta um domínio maior na extensão da curva, para isso foi usado o dobro de iterações da simulação apresentada na Fig. 9. Ao aumentar o número de iterações não se observa mais zonas proibidas, isto e , a curva apresentou apenas dois gaps visíveis, ficando evidenciado após determinar vários intervalos iniciando de valores distintos de k_o e averiguando minuciosamente os diagnosticos.

7. Conclusão

Nesse trabalho abordamos o teorema de Bloch e o modelo de Kronig-Penney seguido de um algoritmo que fornece curvas da energia vs. número de onda para a transmitância constante e variando com a energia, onde o escopo foi esboçar a existência dos gaps de energia em semicondutores. As zonas proibidas foram observadas em ambas as análises. A curva para o coeficiente de transmissão constante atingiu os mesmos resultados das obtidas na literatura. Para a conjuntura que utiliza o coeficiente de transmissão variavel, não foram encontrados dados na literatura para comparação. E deixado um algoritmo completo que segue a mesma id eia do descrito nesse artigo, onde o leitor pode experimentar novos paraômetros em suas proprias simulações.

Agradecimentos

Agradeço a Andre Tanus Ces ario pelo auxílo a esse trabalho, principalmente pelos c alculos apresentados no Apêndice 8.1, e aos pesquisadores Mario Sergio de Carvalho Mazzoni e Luiz Alberto Cury devido ao conhecimento e àas id eias compartilhadas.

Recebido em 2/3/2010; Aceito em 2/8/2011; Publicado em 2/3/2011

8. Apêndices

8.1 . Apêndice: desenvolvimento analítico da Eq. (7)

Solução do sistema descrito pela Eq. (6) fornece

Substituindo as equações para _e (x) e _d (x) (Eqs. (2) e (3)) na Eq. (15), chega-se ao sistema de duas equações com duas variáveis descrito abaixo

Somando e subtraindo as equações acima, tem-se

Multiplicando a primeira equação do sistema (17) por e a segunda equação por tem-se

Pode-se escrever a Eq. (18) da seguinte forma

Definindo o sistema acima na forma matricial como

Assim, o sistema linear homogêneo

pode ser escrito como

Para que um sistema de equações lineares homogêneo tenha solução não -trivial, deve-se ter det (A) = 0

Rearranjando os termos da Eq. (23), obtem-se

Considerando t e r números complexos tais que com sendo e então . Assim, , como , logo . Finalmente, aplicando esse resultado na Eq. (24), tem-se

A Eq. (25) equivale a Eq. (8).

8.2. Apêndice: algoritmo no Scilab

8.2.1. Coeficiente de transmissão constante

• Notação para o algoritmo: k=k e Kappa=.

a=1; dk=0.01; k0=0; delta=0.2; for i=0:1000 do k=k0+i*dk; t=0.8; f=T*cos(k*a+delta); if (abs(f)<=1) then Kappa=acos(f )/a; E=(k^2); printf("%f %f \n",Kappa*a,E); plot(Kappa*a,E,'.'); end end

8.2.2. Coeficiente de transmissão variável

• Notação para o algoritmo: k =k e Kappa =

a=1; dk=0.01; k0=0; U0=3; delta=0.2; for i=0:350 do k=k0+i*dk; E=(k^2); if E>U0 then w=(sin(k*a))^2; x=(4*(E/U0)*((E/U0)-1)); end if E<U0 then T=abs(1+(w/x)); f=T*cos(k*a+delta); if (abs(f)<= 1) then Kappa=(acos(f))/a; printf("%f %f \n",Kappa*a,E) plot(Kappa*a,E,'.') end, end, end,

• [1] J.M. Bass and C.C. Matthai, Phys. Rev. B 52, 4712 (1995).
• [2] C. Kittel, Introduction to Solid State Physics (Wiley, New York, 1996).
• [3] R.F. Pierret, Advanced Semiconductor Fundamentals (Prentice Hall, Upper Saddle River, 2002).
• [4] M.N. Rudden and J. Wilson, Elements of Solid State Physics (Wiley, New York, 1993).
• [5] J. Singleton, Band Theory and Electronic Properties of Solids (Oxford University Press, New York, 2001).
• [6] F. Bloch, Zeitschrift fur Physik A Hadrons and Nuclei, 52, 555 (1929).
• [7] Scilab, The Free Software for Numerical Computation, http://www.scilab.org, accessado em 23/8/2011.
• [8] R. Kronig and W.G. Penney, Proceedings of the Royal Society of London. Series A 130, 499 (1931).
• [9] A. Benzair and H. Aourag, Physica Status Solidi (B) 231, 411 (2002).
• [10] S.H. Liu, Surface Science 231, 633 (1982).
• [11] S.B. Fagan, R. Mota, R.J. Baierle, G. Paiva, A.J.R. da Silva and A. Fazzio, Journal of Molecular Structure: THEOCHEM 539, 101 (2001).
• [12] M.J. Bucknum and R. Hoffmann, Journal of the American Chemical Society 116, 11456 (1994).
• [13] L. Pauling and E.B. Wilson, Introduction to Quantum Mechanics with Applications to Chemistry (Dover Publications, New York, 1985).
• [14] S. Fujita and K. Ito, Quantum Theory of Conducting Matter: Newtonian Equations of Motion for a Bloch Electron (Springer, New York, 2007).
• [15] R. Eisberg and R. Resnick, Quantum Physics ofAtoms, Molecules, Solids, Nuclei, and Particles (Wiley, New York, 1985).
• [16] E. Schrödinger, Journal of the American Chemical Society 116, 11456 (1926).

Datas de Publicação

• Publicação nesta coleção
  16 Maio 2012
• Data do Fascículo
  Mar 2012

Histórico

• Recebido
  02 Mar 2010
• Aceito
  02 Ago 2011