Gás de Elétrons Livres PDF
Gás de Elétrons Livres PDF
Gás de Elétrons Livres PDF
1 Estados eletrônicos
Metais
Os metais possuem propriedades singulares quando comparados com outros
tipos de sólidos. Possuem ductibilidade, maleabilidade e apresentam o ca-
racterı́stico brilho metálico. Entretanto, o que nos interessa particularmente
é que eles são excelentes condutores de eletricidade e calor. A condução
elétrica em metais se deve aos elétrons de condução, os quais também são
responsáveis por grande parte da condução térmica. Se imaginarmos átomos
sendo agrupados para a formar um sólido metálico, vemos que os elétrons
de valência se desligam dos átomos e adquirem a mobilidade para percorrer
todo o sólido e por isso são denominados elétrons de condução. Surpreenden-
temente, a grande maioria dos sólidos simples, compostos por um único tipo
de átomo, apresentem caráter metálico. Alguns sólidos metálicos simples
contribuem com um elétron de condução por átomo, como acontece com os
alcalinos (Li, Na, Rb, Cs), o cobre, a prata e o ouro. Outros contribuem com
dois, como os alcalinos-terrosos (Be, Mg, Ca, Sr, Ba), o ferro, o manganês,
o zinco, o cádmio e o mercúrio no estado sólido. O alumı́nio e o gálio con-
tribuem com três, o estanho e o chumbo com quatro e o bismuto e estanho
com cinco.
Um sólido metálico pode ser visto como constituı́do por ı́ons positivos
localizados nos sı́tios de uma estrutura espacial e por elétrons de condução
não localizados. Os ı́ons não se encontram em repouso absoluto mas vibram
1
em torno de suas posições de equilı́brio. A quantização dos modos normais
de vibração dos átomos dão origem aos fônons. Dessa forma podemos ima-
ginar um sólido metálico como um sistema formado por fônons e elétrons
de condução. Quatro tipos de interações devem ser considerados: (a) entre
elétrons e fônons, (b) entre fônons e fônons, (c) entre elétrons e elétrons, e
(d) entre elétrons e os ı́ons estáticos.
Um modelo simples para um sólido metálico consiste em desprezar a in-
teração entre elétrons e fônons. De acordo com esse modelo, o sólido metálico
se torna equivalente a dois sistemas fracamente interagentes, um deles cons-
tituı́do apenas por fônons e outro apenas por elétrons. Dessa forma, a ca-
pacidade térmica, por exemplo, se torna a soma da capacidade térmica dos
fônons, ou da rede cristalina, e a dos elétrons. A termodinâmica dos fônons
já foi objeto de estudo em capı́tulo anterior. Aqui vamos considerar apenas
a termodinâmica dos elétrons de condução e em particular determinar sua
capacidade térmica.
A parte eletrônica é tratada considerando a como um gás de elétrons não
interagentes, isto é, desprezando as colisões entre elétrons. Isso se justifica
pois as colisões entre elétrons num metal ocorrem com pouca frequência em-
bora a distância entre eles seja da mesma ordem de grandeza da distância
interatômica. Essa incrı́vel propriedade pode ser explicada basicamente a
partir do princı́pio de exclusão de Pauli. Após a colisão os elétrons devem
ocupar orbitais que, pelo princı́pio de exclusão, não estejam ocupados. Essa
restrição em conjunto com a conservação de energia reduzem drasticamente
as possibilidades da ocorrência de orbitais não ocupados que permitam a
colisão.
Densidade de orbitais
O modelo que adotamos para o estudo da parte eletrônica consiste de uma
coleção de elétrons não interagentes mas sujeitos a um potencial devido aos
ı́ons, considerados imóveis em suas posições de equilı́brio. Tendo em vista a
ausência de interação, os orbitais, isto é, os possı́veis estados eletrônicos, são
determinados considerando apenas um elétron sujeito ao potencial devido aos
ı́ons. Num metal cristalino esse potencial é periódico e os estados eletrônicos
correspondem a ondas planas moduladas por funções periódicas e descritas
pelo vetor de onda ~k. Inicialmente estudamos um modelo em que o potencial
devido aos ı́ons é constante. Dessa forma a única energia a ser considerada
é a energia cinética dos elétrons livres.
2
Para determinar os orbitais dos elétrons livres devemos resolver a equação
de Schrödinger independente do tempo relativa a um único életron confinado
numa determinada região que consideramos ser uma região cúbica de lado
L e volume V = L3 . O elétron possui apenas energia cinética de modo que
equação de Schrödinger independente do tempo para a função de onda ψ do
elétron é dada por
h̄2 2
− ∇ ψ = ǫψ, (1)
2m
em que m é a massa do elétron, h̄ = h/2π sendo h a constante de Planck e ǫ
é o autovalor da energia do elétron. As autofunções são dadas por
1
ψ(x, y, z) = √ ei(kx x+ky y+kz z) (2)
V
e os autovalores por
h̄2 k 2
ǫ~k = , (3)
2m
que é também chamada relação de dispersão.
Os possı́veis valores do vetor de onda ~k = (kx , ky , kz ) são determina-
dos pelas condições de contorno que adotamos como sendo periódicas. Para
satisfazer essas condições devemos impor que ψ(x + L, y, z) = ψ(x, y, z),
ψ(x, y + L, z) = ψ(x, y, z) e ψ(x, y, z + L) = ψ(x, y, z), que resulta em
h̄2 2
ǫ~k = (k + ky2 + kz2 ) ≤ ǫ (6)
2m x
3
e multiplicando o resultado por 2 pois a cada vetor de onda corresponde dois
orbitais. Se associarmos a cada possı́vel vetor um ponto no espaço (kx , ky , kz ),
então o número de vetores de onda que satisfazem (6) é igual ao número de
cubos elementares de volume (2π/L)3 = 8π 3 /V dentro da esfera de raio
q
2mǫ/h̄2 e portanto
3/2
4 2mε
3
π h̄2 V
2m
3/2
N (ǫ) = 2 = 2 ǫ3/2 . (7)
8π 3
V
3π h̄2
2 Propriedades termodinâmicas
Distribuição de Fermi-Dirac
Para determinar as propriedades termodinâmicas de um gás de elétrons livres
consideramos que essas partı́culas estejam confinadas numa região de volume
V à temperatura T em contato com um reservatório de partı́culas que fixa o
potencial quı́mico µ dos elétrons. Tendo em vista que os elétrons são férmions
o número de ocupação médio hns i de um orbital s é dado pela distribuição
de Fermi-Dirac
1
hns i = β(ǫs −µ) = f (ǫs ), (9)
e +1
em que ǫs é a energia do orbital s. O grande potencial termodinâmico
Φ(T, V, µ) é dado por
X
Φ = −kB T ln(1 + f (ǫs )). (10)
s
e X
N= f (ǫ). (12)
s
4
f
1
1/2
0
0 µ ε
5
Gás de elétrons a temperatura nula
À temperatura zero, a função de Fermi vale
1, ǫ < µ,
f (ǫ) = 1/2, ǫ = µ, (19)
0, ǫ > µ,
que equivale a dizer que os orbitais com energia menores do que µ estão
ocupados e equeles com energia maiores do que µ estão vazios. Substituindo
esse resultado em (15), vemos que o número de elétrons N é dado por
Z µ
N= D(ǫ)dǫ. = N (µ), (20)
0
6
e depende apenas do número de elétrons por unidade de volume N/V . A
temperatura de Fermi TF é definida por TF = ǫF /kB .
A razão entre U, dada por (23), e N, dada por (21), fornece o resultado
U/N = (3/5)µ ou, tendo em vista que o potencial quı́mico se identifica com
a energia de Fermi,
U 3
= ǫF . (26)
N 5
À temperatura zero, a energia média por elétron é igual a três quintos da
energia de Fermi.
Resumindo, à temperatura zero, o potencial quı́mico, a energia e a pressão
de um gás de elétrons livres, como funções de N e V , são dados por
3 2
µ = ǫF , U = NǫF , pV = NǫF . (27)
5 5
O último resultado foi obtido usando pV = (2/3)U e nos revela que a pressão
de um gás de elétrons livres a temperatura zero é não nula.
Expansão de Sommerfeld
Para determinar as propriedades termodinâmicas tais como a energia e o
número médio de partı́culas, dadas pelas expressões (14) e (15), devemos
efetuar integrais do tipo
∞ ∞ g(ǫ)
Z Z
I(T ) = g(ǫ)f (ǫ)dǫ = dǫ, (28)
0 0 eβ(ǫ−µ) + 1
em que g(ǫ) não depende de T nem de µ. Para baixas temperaturas desejamos
mostrar a seguinte expansão assintótica
π2 ′ µ
Z
I(T ) = g (µ)(kB T )2 ,
g(ǫ)dǫ + (29)
0 6
denominada expansão de Sommerfeld, válida até termos quadráticos em
kB T /µ.
Inicialmente vemos que
Z µ
I(0) = g(ǫ)dǫ, (30)
0
que subtraı́do de (28) resulta numa expressão que pode ser escrita como
µ g(ǫ) ∞ g(ǫ) µ
Z Z Z
I(T ) − I(0) = dǫ + dǫ − g(ǫ)dǫ, (31)
0 eβ(ǫ−µ) + 1 µ eβ(ǫ−µ) + 1 0
7
gf
B
A
0
0 µ ε
ou
I(T ) − I(0) = IA − IB , (32)
em que
g(ǫ)
Z ∞
IA = dǫ, (33)
µ eβ(ǫ−µ) + 1
que é igual à a área da região A da figura 2, e
µ µg(ǫ) µ g(ǫ)
Z Z Z
IB = g(ǫ)dǫ − β(ǫ−µ)
dǫ = dǫ (34)
0 0 e +1 0 1 + e−β(ǫ−µ)
que é igual à área da região B da figura 2. Vale notar que as áreas das regiões
A e B não são iguais. A área de A é ligeiramente maior do que a área de B.
Em seguida efetuamos as seguintes mudanças de variáveis. Na integral
IA fazemos a mudança ξ = β(ǫ − µ), de modo que ǫ = µ + kB T ξ, para obter
Z ∞ g(ǫ) Z ∞
g(µ + kB T ξ)
IA = β(ǫ−µ)
dǫ = k B T dξ. (35)
µ e +1 0 eξ + 1
Na integral IB fazemos a mudança ξ = −β(ǫ−µ), de modo que ǫ = µ−kB T ξ,
para obter
µ g(ǫ) βµ g(µ − kB T ξ)
Z Z
IB = dǫ = kB T dξ. (36)
0 1 + e−β(ǫ−µ) 0 1 + eξ
8
Como estamos interessados numa expansão de baixas temperaturas po-
demos estender o limite superior da integral IB até o infinito e assim obter
a seguinte expressão para a diferença entre as áreas das regiões A e B,
∞ g(µ + kB T ξ) − g(µ − kB T ξ)
Z
IA − IB = kB T dξ. (37)
0 eξ + 1
A expressão no numerador pode ser expandida em série de potências de ξ.
A expansão até termos lineares em ξ vale
g(µ + kB T ξ) − g(µ + kB T ξ) = 2g ′(µ)kB T ξ, (38)
que nos leva ao seguinte resultado
∞ ξ π2 ′
Z
2
′
IA − IB = 2g (µ)(kB T ) ξ
dξ = g (µ)(kB T )2 , (39)
0 e +1 6
em que levamos em conta que a integral vale π 2 /12. Como I(T ) − I(0) =
IA − IB , obtemos finalmente
π2 ′
I(T ) = I(0) + g (µ)(kB T )2 , (40)
6
que é equivalente à expressão desejada (29) tendo em vista o resultado (30)
para I(0).
9
válidas para kB T << µ.
Experimentalmente as propriedades termodinâmicas são medidas a partir
de uma amostra de sólido que obviamente possui um número de elétrons que
é mantido inalterado. Para obter as grandezas termodinâmicas como funções
de N, em particular a capacidade térmica Cv , eliminamos o potencial quı́mico
das expressões obtidas acima para obter U como função de T, V e N. Com
essa finalidade começamos por obter µ como função de T, V e N. A partir
de (42) e usando o resultado (24), obtemos
!2
3π 2 kB T
µ = ǫF {1 + }−2/3 . (43)
24 µ
Nessa equação, substituı́mos o potencial quı́mico que está fora dos parênteses
pela expressão dada por (44) e aqueles que estão entre parênteses por ǫF =
kB TF . Esse procedimento resulta na expressão
2 2 2
U π2 T 15π 2 T 3π 2 T
= ǫF {1 − }{1 + }{1 + }−1 . (46)
N 12 TF 24 TF 24 TF
10
Tabela 1: Energia de Fermi ǫF calculada por (25), temperatura de Fermi
TF = ǫ/kB e razão entre o coeficiente do calor especı́fico eletrônico observado
experimentalmente γ ∗ e o coeficiente γ dado por (52), para vários metais.
Kittel, 1971.
que é a expressão desejada para U. Vale a pena notar que a pressão dos
elétrons é dada por pV = 2U/3, isto é,
2 5π 2 T 2
pV = ǫF {1 + }. (48)
5 12 TF
Para determinar a capacidade térmica Cv basta efetuar a derivada Cv =
(∂U/∂T )V,N . Derivando U relativamente à temperatura,
π2 T
Cv = NkB . (49)
2 TF
Esse resultado nos diz que a capacidade térmica de elétrons livres degenerados
possui comportamento linear com a temperatura. A medida da capacidade
térmica eletrônica nos permite obter a temperatura de Fermi.
11
entre elétrons e fônons sejam desprezı́veis. Somando a capacidade térmica
eletrônica, dada por (49), com a capacidade térmica da rede, calculado de
acordo com o modelo de Debye, vista no capı́tulo 4, obtemos o seguinte
resultado para o calor especı́fico cv de metais a baixas temperaturas
cv = γT + αT 3 , (50)
em que
12π 4 kB
α= (51)
5 TD3
e
π 2 kB
γ= z , (52)
2 TF
sendo z é o número de elétrons livres por átomo. O calor especı́fico é definido
como a capacidade térmica a volume constante dividido pelo número de
átomos.
Um gráfico dos dados experimentais de cv /T como função de T 2 deve
ser linear. O coeficiente linear determina α e a intersecção com a ordenada
determina γ. O resultados experimentais para gamma, que denotamos por
γ ∗ , diferem de γ dado por (52), embora sejam da mesma ordem de grandeza.
Como a temperatura de Fermi é inversamente proporcional à massa m do
elétron, o coeficiente γ é proporcional m. A discrepância pode ser explicada
admitindo que os elétrons de condução possuem massa efetiva m∗ , que pode
ser distinta de m e definida por
γ∗
m∗ = m. (53)
γ
Exercı́cios
1. A relação de dispersão de um gás de elétrons livres ultra-relativı́ticos é
dada por q
ǫ~k = ch̄ kx2 + ky2 + kz2 .
12
Determine o número e a densidade de orbitais. Mostre que a pressão p se
relaciona com a energia U por meio de pV = U/3. Determine a energia U no
regime de elétrons livres degenerados e mostre que a a capacidade térmica
vale
T
Cv = π 2 NkB .
TF
2. Mostre que o potencial quı́mido de um gás de elétrons livres degenerados
é dado por
π 2 D ′ (ǫF )
µ = ǫF − (kB T )2 .
6 D(ǫF )
A partir dessa expressão obtenha o resultado (44).
3. Use a identidade termodinâmica
! ! !
∂U ∂U ∂U (∂N/∂T )µ
= −
∂T N
∂T µ
∂µ T
(∂N/∂µ)T
π2 2
Cv = D(ǫF )kB T.
3
A partir dessa expressão obtenha o resultado (49).
4. Determine a compressibilidade isotérmica κT = −(1/V )(∂V /∂p)T e o
coeficiente de expansão térmica α = (1/V )(∂V /∂T )p de um gás de elétrons
livres degenerados. A partir desses resultados ache a capacidade térmica a
pressão constante Cp = Cv + T V α2 /κT e determine a razão Cp /Cv .
5. Determine a energia livre de Helholtz F de um gás de elétrons degenerados.
A partir dela ache a entropia S = −∂F/∂T .
6. Mostre que para um gás de elétrons livres a razão S/N entre a entropia
e o número de elétrons só depende da razão µ/T entre o potencial quı́mico
e a temperatura. A partir desse resultado mostre que pV 5/3 é constante ao
longo de uma adiabática.
13