Estatística de Bose-Einstein

Termodinâmica
Processo de convecção.
Glossário de termodinâmica
Ramos Clássica Estatística Química
Leis Zero Primeira Segunda Terceira
Sistema termodinâmico Sistema termodinâmico Estado Equação de estado Gás ideal Gás real Estados físicos da matéria Equilíbrio Volume de controle Transformações Transformação termodinâmica Isobárica Isocórica Isotérmica Adiabática Isentrópica Isentálpica Politrópica Expansão livre Reversibilidade Irreversibilidade Ciclos Ciclo termodinâmico Máquinas térmicas Ciclo de refrigeração Eficiência termodinâmica
Grandezas físicas Funções de estado Nota: Variáveis conjugadas em itálico Propriedades extensivas e intensivas Calor Calor específico Condutividade térmica Densidade Energia Entropia Fluxo Potência Potencial químico Pressão Temperatura Título Trabalho Vazão Volume
Propriedades Capacidade térmica específica  $c=$ $T$ $\partial S$ $N$ $\partial T$ Compressibilidade  $\beta =-$ $1$ $\partial V$ $V$ $\partial p$ Dilatação térmica  $\alpha =$ $1$ $\partial V$ $V$ $\partial T$
Equações termodinâmicas Teorema de Carnot Lei dos gases ideais Relações de Maxwell Teorema da equipartição Teorema do virial Lei de Henry Lei de Boyle-Mariotte Equação de Van der Waals Raio de Van der Waals Lei de Charles
Potencial termodinâmico Potencial termodinâmico Energia livre Energia interna $U(S,V)$ Entalpia $H(S,p)=U+pV$ Energia livre de Helmholtz $A(T,V)=U-TS$ Energia livre de Gibbs $G(T,p)=H-TS$
História Cultura História Geral Leis dos gases Filosofia Demônio de Maxwell Teorias Teoria calórica Teoria do calor Vis viva ("força viva") Equivalente mecânico do calor Potência motriz
Cientistas Bernoulli Boltzmann Carnot Clapeyron Clausius Carathéodory Duhem Gibbs von Helmholtz Joule Maxwell von Mayer Onsager Rankine Smeaton Stahl Thompson Thomson van der Waals Waterston

Em um sistema quântico constituído de muitas partículas idênticas com spin inteiro, a estatística de Bose-Einstein, ou estatística BE, é utilizada para descrever o sistema e calcular os valores médios das grandezas físicas.

Em um sistema de $N$ bósons idênticos de massa $m$, que possuem interação mútua desprezível, contidos em um recipiente de volume $V$, a uma temperatura $T$, em equilíbrio, o número médio de partículas ${\bar {n}}_{s}$ num estado de energia $s$ é dado por

${\bar {n}}_{s}={\frac {g_{s}}{e^{\beta (\epsilon _{s}-\mu )}-1}}$ ,

em que $g_{s}$ é a degenerescência quântica do estado $s$, $\epsilon _{s}$ é a energia do estado $s$, $\mu$ é o potencial químico, e $\beta =1/k_{B}T$, em que $k_{B}$ é a constante de Boltzmann.

Formulação matemática

Seja um gás de $N$ partículas idênticas confinadas em um volume $V$. Sendo $Q_{i}$ o conjunto das coordenadas generalizadas da i-ésima partícula e $s_{i}$ o índice rotulador dos possíveis estados quânticos desta única partícula, o estado do gás inteiro é então descrito pelo conjunto de números quânticos

$\{s_{1},s_{2},\dots ,s_{N}\}$

os quais definem a função de onda $\Psi$ do gás nesse estado.

$\Psi =\Psi _{\{s_{1},s_{2},\dots ,s_{N}\}}(Q_{1},Q_{2},\dots ,Q_{N})$

A estatística BE trata de partículas em que o spin total é medido em unidade inteiras de $\hbar$ — os chamados bósons — como, por exemplo, o átomo de Hélio-4 e o fóton. Partículas quânticas são indistinguíveis, ou seja, a troca de duas partículas não altera o estado do sistema. Além disto, a função de onda total para um sistema de bósons é simétrica sobre a troca de duas partículas

$\Psi (\dots ,Q_{i},\dots ,Q_{j})=\Psi (\dots ,Q_{j},\dots ,Q_{i})$

Como consequência, os bósons não obedecem ao princípio de exclusão de Pauli, não havendo limitação no número de partículas que podem ocupar um dado estado quântico. Tomando este gás a uma temperatura $T$ e que cada partícula está em um estado $r$ de energia $\epsilon _{r}$ onde há $n_{r}$ partículas neste estado, o gás inteiro possui $R$ estados possíveis. Desprezando as interações mútuas entre as partículas, a energia do gás no estado $R$ será dada por

$E_{R}=n_{1}\epsilon _{1}+n_{2}\epsilon _{2}+\dots =\sum _{r}n_{r}\epsilon _{r}$

em que a soma se estende sobre todos os estados $r$ possíveis de uma partícula. Como o número de partículas é fixo, também temos que

$\sum _{r}n_{r}=N$

A função partição $Z$ do gás será dada por

$Z=\sum _{R}e^{-\beta E_{R}}=\sum _{R}e^{-\beta n_{1}\epsilon _{1}+\dots }$

Essa soma é sobre todos os estados $R$ possíveis do gas inteiro, isto é, sobre todos os possíveis números $n_{1},n_{2},n_{3},\dots$. Como exp$[-\beta (n_{1}\epsilon _{1}+n_{2}\epsilon _{2}+...)]$ é a probabilidade relativa de encontrar o gás em um estado partículas onde há $n_{1}$ partículas em 1, $n_{2}$ em partículas em 2 e assim sucessivamente, pode-se escrever o número médio de partículas em um estado $s$ como

${\bar {n_{s}}}={\dfrac {\sum _{R}n_{s}e^{-\beta (n_{1}\epsilon _{1}+\dots )}}{\sum _{R}e^{-\beta (n_{1}\epsilon _{1}+\dots )}}}.$

Somando todos os possíveis valores de $n_{s}$, usando a propriedade multiplicativa da exponencial e rearranjando, pode-se escrever ${\bar {n_{s}}}$ como

${\bar {n_{s}}}={\dfrac {\displaystyle \sum _{n_{s}}n_{s}e^{-\beta n_{s}\epsilon _{s}}\cdot \sum _{n_{1},n_{2},\dots }^{(s)}e^{-\beta (n_{1}\epsilon _{1}+\dots )}}{\displaystyle \sum _{n_{s}}e^{-\beta n_{s}\epsilon _{s}}\cdot \sum _{n_{1},n_{2},\dots }^{(s)}e^{-\beta (n_{1}\epsilon _{1}+\dots )}}}$

em que o sobrescrito $(s)$ no somatório indica a soma com exceção do estado $s$ em particular. Adotando a notação:

$\sum _{n_{1},n_{2},\dots }^{(s)}e^{-\beta (n_{1}\epsilon _{1}+\dots )}=Z_{s}(N)$

A restrição de $N$ fixo implica que se uma partícula está no estado $s$, a soma $\sum ^{(s)}$ se estende pelas $N-1,N-2,\dots$ partículas restantes que podem ser colocadas nos estados $n_{r}=1,2,3,\dots$. Ao executar explicitamente a soma sobre ${\bar {n_{s}}}$ ter-se-á

${\bar {n_{s}}}={\dfrac {0+e^{-\beta \epsilon _{s}}Z_{s}(N-1)+2e^{-2\beta \epsilon _{s}}Z_{s}(N-2)+\dots }{Z_{s}(N)+e^{-\beta \epsilon _{s}}Z_{s}(N-1)+2e^{-2\beta \epsilon _{s}}Z_{s}(N-2)+\dots }}$

Estabelecendo a seguinte relação entre $Z_{s}(N-\Delta N)$ e $Z_{s}(N)$ onde $N>>\Delta N$

$\ln {Z_{s}(N-\Delta N)}=\ln {Z_{s}(N)}-{\dfrac {\partial \ln {Z_{s}}}{\partial N}}\cdot \Delta N=\ln {Z_{s}}-\alpha _{s}\Delta N$

em que

$\alpha _{s}={\dfrac {\partial \ln {Z_{s}}}{\partial N}}$

e portanto,

$Z_{s}(N-\Delta N)=Z_{s}(N)e^{-\alpha _{s}\Delta N}$

Mas como $Z_{s}(N)$ é uma soma sobre muitos estados, espera-se que a variação de seu logaritmo natural com o número total de partícula seja imperceptível para o qual um estado particular $s$ seja omitido da soma. Vamos introduzir a aproximação de que $\alpha _{s}$ é independente de $s$, então podemos só escrever $\alpha _{s}=\alpha$ para todo $s$. Logo,

$\alpha ={\frac {\partial \ln {Z}}{\partial N}}$

Do ensemble canônico, sabe-se que essa derivada parcial resulta em $-\beta \mu$, em que $\mu$ é o potencial químico do gás. Portanto,

$\alpha =-\beta \mu$

Substituindo estes resultados em ${\bar {n_{s}}}$ teremos

${\bar {n_{s}}}={\dfrac {e^{-\beta (\epsilon _{s}-\mu )}Z_{s}(N)+2e^{-2\beta (\epsilon _{s}-\mu )}Z_{s}(N)+\dots }{Z_{s}(N)+e^{-\beta (\epsilon _{s}-\mu )}Z_{s}(N)+2e^{-2\beta (\epsilon _{s}-\mu )}Z_{s}(N)+\dots }}$

Cancelando $Z_{s}(N)$ ter-se-á

${\bar {n_{s}}}={\dfrac {e^{-\beta (\epsilon _{s}-\mu )}+2e^{-2\beta (\epsilon _{s}-\mu )}+\dots }{1+e^{-\beta (\epsilon _{s}-\mu )}+2e^{-2\beta (\epsilon _{s}-\mu )}+\dots }}$

Que é uma série geométrica. Usando a função partição de acordo com isso resultará em

${\bar {n_{s}}}={\frac {1}{e^{\beta (\epsilon _{s}-\mu )}-1}}$

Se o estado $s$ é degenerado com degenerescência $g_{s}$ então o número médio de partículas com energia $\epsilon _{s}$ é obtido multiplicando a expressão anterior por $g_{s}$. Finalmente,

${\bar {n_{s}}}={\frac {g_{s}}{e^{\beta (\epsilon _{s}-\mu )}-1}}$

A estatística de Bose-Einstein reduz-se à estatística de Maxwell-Boltzmann para energias: $(\epsilon _{s}-\mu )>>kT$.

Gás de fótons

Um caso especial da estatística de Bose-Einstein é o gás de fótons. Fótons possuem spin inteiro igual a 1, então desta forma são considerados bósons. O caso é especial devido ao fato de que se considerarmos vários fótons dentro de um recipiente com volume V, o número destes fótons não será constante, pois conforme estes fótons interagem com as paredes do recipiente estes são absorvidos ou emitidos. Desta forma, não podemos impor um vínculo ao número total de fótons no sistema. Neste caso, precisaremos realizar as somas sobre todos os possíveis números de partículas em cada estado, da forma:

$n_{r}=0,1,2,3,...$ para todo r.

A função partição para o gás de fótons é dada por:

$Z=\sum _{R}e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+...)}$

No qual R são todos os estados possíveis do gás. Como não há vínculos agora para o número de partículas por estado, podemos utilizar as propriedades das funções exponenciais e reescrever a soma acima como:

$Z=\sum _{n_{1},n_{2},...}e^{-\beta n_{1}\varepsilon _{1}}e^{-\beta n_{2}\varepsilon _{2}}...$

Ou ainda:

$Z=\left(\sum _{n_{1}=0}^{\infty }e^{-\beta n_{1}\varepsilon _{1}}\right)\left(\sum _{n_{2}=0}^{^{\infty }}e^{-\beta n_{2}\varepsilon _{2}}\right)...$

Como não há restrição para o número de fótons, as somas acima são consideradas até um número muito grande de partículas por estado de energia, no qual matematicamente isto se traduz a realizarmos a soma até o infinito, embora fisicamente, estejamos carregando a soma até um número muito grande de partículas. Como estamos tratando de uma exponencial com argumento negativo, após um certo valor os termos serão desprezíveis, não tendo problemas com divergências. Se olharmos com cuidado para as somas dentro dos colchetes acima, percebe-se que podemos escrevê-las como abaixo:

$\sum _{n_{1}=0}^{\infty }e^{-\beta n_{1}\varepsilon _{1}}=1+e^{-\beta \varepsilon _{1}}+e^{-2\beta \varepsilon _{1}}+...=\left({\frac {1}{1-e^{-\beta \varepsilon _{1}}}}\right)$

Ou seja, cada termo é uma soma geométrica. Assim, podemos reescrever a função partição da seguinte forma:

$Z=\left({\frac {1}{1-e^{-\beta \varepsilon _{1}}}}\right)\left({\frac {1}{1-e^{-\beta \varepsilon _{2}}}}\right)...$

E o logaritmo natural da função partição, que é o qual estamos interessados é dado por:

$\ln {Z}=-\sum _{r}\ln {(1-e^{-\beta \varepsilon _{r}})}$

Nosso próximo objetivo agora é encontrar o número médio de partículas em um estado de energia $\varepsilon _{s}$. Tal resultado pode ser obtido através da expressão abaixo:

${\bar {n_{s}}}=-{\frac {1}{\beta }}{\frac {\partial \ln {Z}}{\partial \varepsilon _{s}}}$

Realizando a derivação do logaritmo natural da função partição, temos o resultado abaixo:

${\bar {n_{s}}}={\frac {e^{-\beta \varepsilon _{s}}}{1-e^{-\beta \varepsilon _{s}}}}$

O resultado acima é este pois a única derivada que não é zero é o termo da soma no qual $r=s$. Podemos simplificar o resultado acima, multiplicando e dividindo a expressão acima pela exponencial com o mesmo argumento, porém positivo, e assim obtemos o importante resultado dado por:

${\bar {n_{s}}}={\frac {1}{e^{\beta \varepsilon _{s}-1}}}$

Este resultado é conhecido como a distribuição de Planck, e fornece o número médio de fótons em um determinado estado s. Uma das aplicações mais famosas do resultado acima é no problema da radiação de corpo negro.

Radiação de corpo negro

Ver artigo principal: Radiação de corpo negro

Radiância espectral em função da frequência

Todo corpo a uma temperatura $T$ emite radiação eletromagnética. A distribuição de Planck fornece o espectro de emissão para uma classe especifica de corpos, os chamados corpos negros, definidos como os corpos que absorvem toda a radiação incidente. Pode-se modelar um corpo negro como uma cavidade metálica com volume $V$, tal que haja apenas um pequeno orifício em uma de suas paredes. Logo, esta cavidade absorve toda a radiação que entra por ali, e radiação emitida pelo orifício que é oriunda das emissões a partir das superfícies internas da cavidade se comporta como se fosse um corpo negro.

Busca-se a chamada radiância espectral $R_{T}(f)df$, que fornece a potência irradiada por unidade de área com frequência entre $f$ e $f+df$ pelo corpo estudado a uma data temperatura. Através de uma análise física do problema, pode-se mostrar que a radiância espectral está diretamente ligada com a densidade de energia dentro da cavidade. A relação entre as duas grandezas é dada por:

$R_{T}(f)={\frac {c}{4}}u_{T}(f)$

Pode-se então obter a densidade de energia, e assim resolve-se o problema da mesma maneira. A densidade de energia pode ser obtida a partir da probabilidade de um nível com energia $\varepsilon _{s}$ estar ocupado por $n_{s}$ fótons, sendo assim se multiplicarmos este valor pelo número de médio fótons por unidade de volume naquele estado, teremos a densidade de energia dentro desta cavidade na forma:

$u_{T}(\varepsilon )={g(\varepsilon ){\bar {\varepsilon _{s}}}}$

Onde $g(\varepsilon )$ é a densidade de estados ou degenerescência. Como temos o número médio de fótons em um estado $s$, basta multiplicar este número pela energia do estado $s$, assim:

${\bar {\varepsilon _{s}}}={\frac {\varepsilon _{s}}{e^{\beta \varepsilon _{s}}-1}}$

Como os estados do sistema estão muito próximos um dos outros, por conta da cavidade ocupar um volume V macroscópico, podemos tratar as variáveis como sendo contínuas. O número de estados por unidade de volume dentro da cavidade com frequência entre $f$ e $f+df$ é dada por:

$g(f)df={\frac {8\pi }{c^{3}}}f^{2}df$

Sabemos que a energia e a frequência de um fóton estão ligado pela expressão de Planck $E=hf$, desta forma se fizermos a substituição sugerida, obtemos assim a expressão para a densidade de energia dentro de uma cavidade.

$u_{T}(f)={\frac {8\pi h}{c^{3}}}{\frac {f^{3}}{e^{\frac {hf}{kT}}-1}}$

E a radiância espectral é dada diretamente por:

$R_{T}(f)={\frac {2\pi h}{c^{2}}}{\frac {f^{3}}{e^{\frac {hf}{kT}}-1}}$

Todo o caminho feito para o obtenção do resultado foi a partir da análise das propriedades quânticas e estatísticas de fótons dentro de um volume V, e em equilíbrio térmico a uma temperatura T. Pode-se chegar no mesmo resultado analisando a interação da radiação eletromagnética dentro do volume V com as paredes do recipiente.

Comparação entre as estatísticas

Comparação das distribuições de Maxwell-Boltzmann (MB), Bose-Einstein (BE) e Fermi-Dirac (FD).

Influência do potencial químico na distribuição BE.

Influência da temperatura na distribuição de BE.

Nota-se que um sistema de bósons possui uma maior ocupação dos estados, devido aos efeitos quânticos de simetria da função de onda do sistema . Nesta distribuição identifica-se que para o regime $(\epsilon -\mu )/kT>>1$ as distribuições quânticas de Bose-Einstein (BE) e Fermi-Dirac (FD) se aproximam da distribuição clássica de Maxwell-Boltzmann (MB), que representa o regime de baixas densidades e altas temperaturas, ou seja, o limite clássico .

Em diferentes valores do potencial químico, nota-se que os estados de menor energia (próximos ao estado fundamental) são os mais populados. Para baixas temperaturas as partículas se concentram nos estados de menor energia. No limite de temperatura tendendo a zero, todas as partículas vão para o estado fundamental ($s=1$), logo, ${\bar {n}}_{1}\approx N$ enquanto ${\bar {n}}_{s}\approx 0$ para os demais estados . Já para altas temperaturas, devido à energia térmica do sistema, as partículas tem maior probabilidade de atingir estados mais energéticos.

Condensação de Bose-Einstein

O condensado de Bose-Einstein é uma fase da matéria formada por bósons a uma temperatura muito próxima do zero absoluto. Nestas condições, uma grande fração de átomos atinge o mais baixo estado quântico, e nestas condições os efeitos quânticos podem ser observados em escala macroscópica. Sistemas em baixa temperatura ou com densidade relativamente alta de partículas são mais prováveis de apresentarem comportamentos quânticos, mesmo em sistemas onde a interação intramolecular é desprezível.

Temperatura crítica para um gás ideal de bósons

Um gás ideal de bósons não está sujeito ao princípio de exclusão de Pauli. Logo, os bósons podem se condensar no seu estado de menor energia. A densidade de estados $g(\epsilon )$ é dada por,

$g(\epsilon )={\frac {V}{4\pi ^{2}}}{\frac {(2m)^{\frac {3}{2}}}{\hbar ^{3}}}\epsilon ^{1/2}$

Logo, o número de partículas $N$ pode ser reescrito como

$N=\int _{0}^{\infty }{\frac {g(\epsilon )}{e^{\beta (\epsilon -\mu )}-1}}d\epsilon$

A distribuição continua pode ser utilizada pois os níveis do sistema quântico discreto são numerosos e estão muito próximos. Portanto,

$N={\frac {V}{4\pi ^{2}}}{\frac {(2m)^{\frac {3}{2}}}{\hbar ^{3}}}\int _{0}^{\infty }{\frac {\epsilon ^{1/2}}{e^{\beta (\epsilon -\mu )}-1}}d\epsilon$

Porém, como a integração dar-se-á de $\epsilon =0$, a informação do estado fundamental é perdida. Todavia, quando a temperatura do sistema diminui, o potencial químico aumenta, e o número de partículas no estado fundamental é dada por

$N_{0}={\frac {1}{e^{-\beta \mu }-1}}$

Seja o número de partículas em estados excitados dado por $N_{e}$, temos

$N=N_{0}+N_{e}$

$N={\frac {1}{e^{-\beta \mu }-1}}+{\frac {V}{4\pi ^{2}}}{\frac {(2m)^{\frac {3}{2}}}{\hbar ^{3}}}\int _{0}^{\infty }{\frac {\epsilon ^{1/2}}{e^{\beta (\epsilon -\mu )}-1}}d\epsilon$

Porém, para temperaturas muito próximas de zero $N\approx N_{0}$, ou seja,

$N\approx {\frac {1}{e^{-{\frac {\mu }{kT}}}-1}}$

Variação da temperatura de No/N (preto) e de Ne/N (vermelho) para uma gás ideal de bósons.

E, considerando $N$ grande,

$-{\frac {\mu }{kT}}\approx \ln \left(1+{\frac {1}{N}}\right)\approx {\frac {1}{N}}$

Logo, $e^{-\mu /kT}\approx 1$. E $N_{e}$ se torna,

$N_{e}={\frac {V}{4\pi ^{2}}}{\frac {(2m)^{\frac {3}{2}}}{\hbar ^{3}}}\int _{0}^{\infty }{\frac {\epsilon ^{1/2}}{e^{\beta \epsilon }-1}}d\epsilon$

Introduzindo a mudança de variável $x=\beta \epsilon$, temos

$N_{e}={\frac {2}{\sqrt {\pi }}}\left({\frac {2\pi mkT}{h^{2}}}\right)^{3/2}V\int _{0}^{\infty }{\frac {x^{1/2}}{e^{x}-1}}d\epsilon$

A integral pode ser escrita em termos da função gama $\Gamma (x)$ e da função zeta de Riemann $\zeta (x)$. De fato, ela é igual a $\Gamma (3/2)\zeta (3/2)$, onde $\Gamma (3/2)={\sqrt {\pi }}/2$, logo,

$N_{e}=\zeta (3/2)V\left({\frac {2\pi mkT}{h^{2}}}\right)^{3/2}$

A temperatura crítica, ou de Bose-Einstein, $T_{c}$, é a temperatura em que acima dela todos os bósons estão em um estado excitado, e pode ser encontrada tomando $N=N_{e}$, onde $T=T_{c}$, assim, da última equação,

$T_{c}={\frac {h^{2}}{2\pi mk}}\left[{\frac {N}{\zeta (3/2)V}}\right]^{2/3}$

Se $0<T<T_{c}$, a razão de bósons no estado excitado em relação ao total, é

${\frac {N_{e}}{N}}=\left({\frac {T}{T_{c}}}\right)^{3/2}$

Consequentemente, a razão para os bósons no estado fundamental é,

${\frac {N_{0}}{N}}=1-\left({\frac {T}{T_{c}}}\right)^{3/2}$

Hélio líquido

No caso do hélio líquido, que possui uma massa específica de 0,124 g/cm$^{3}$ na temperatura de ebulição de 4,2K, sua densidade é 1,87.10$^{28}$ átomos /m$^{3}$. Tal densidade é muito alta para que a distribuição de Boltzmann descreva o sistema corretamente .

Considerando o isótopo mais comum do hélio, o He$^{4}$, temos um átomo formado por dois prótons, dois neutrons e dois elétrons, todos férmions de spin 1/2. Entretanto, a combinação dos spins dessas partículas resulta em um átomo com spin zero, ou seja, o He$^{4}$ é um átomo de spin inteiro, ou seja, um bóson .

Os átomos de hélio interagem apenas através da força de van der Waals, uma vez que é um gás nobre. Essa força é fraca, e como a distância entre os átomos é grande, o mesmo pode ser modelado de maneira simplificada como um gás ideal, como foi sugerido por F. London em 1938 .

A temperatura crítica para o hélio líquido, dado pela equação da seção anterior, é de aproximadamente 3,14 K, próxima da temperatura do ponto $\lambda$ do hélio superfluido (2,17K). Isto é uma consequência da semelhança que a fase superfluida do hélio possui com um gás ideal de Bose na fase condensada, mesmo sendo pouco realista desprezar as interações entre as partículas num líquido quântico .

Bibliografia

• Eisberg, R.M (1974). Física Quântica, 4º edição. [S.l.]: Editora Campus LTDA. ISBN 8570013094 
• Griffiths, D.J (2005). Introduction to quantum mechanics, 2º edição. [S.l.]: Pearson Education. ISBN 0131118927 
• Pathria, R. K (1996). Statistical mechanics, second edition. [S.l.]: Elsevier. ISBN 0750624698 
• Reif, F (1965). Fundamentals of statistical and thermal physics. [S.l.]: McGraw-Hill. ISBN 9780070856158 
• Salinas, S. R. A (1999). Física estatística. [S.l.]: Editora da Universidade de São Paulo. ISBN 8531403863 
• Tipler, P. A, Llewellyn, R. A. (2001). Física Moderna, 3º edição. [S.l.]: LTC - Livros Técnicos e Científicos S.A. ISBN 9788521626077  !CS1 manut: Nomes múltiplos: lista de autores (link)