Gás de Bose

Em mecânica estatística, um gás de bósons é um sistema ideal de partículas que obedece à estatística de Bose-Einstein.

Potencial termodinâmico

Devido aos efeitos de troca, a maneira mais simples de trabalhar com gases quânticos é com o Ensemble Grand-Canônico:

<math>\mathcal{Z} = \sum_{N=0}^{\infty}\sum_{i}z^{N}e^{-\beta \varepsilon_i}</math>

que para um gás fica:

<math>\mathcal{Z} = \sum_{N=0}^{\infty}\sum_{\{n_i\}}\,^{\prime}\prod_{i}z^{n_i}e^{-\beta n_iE_i}</math>

A segunda soma é restrita ao número total de partículas ser <math>N</math>. Uma maneira de fazer tal soma é somar primeiro sobre todos os <math>N</math> possíveis e depois multiplicar todos os níveis. Para um sistema de bósons, qualquer valor de <math>N</math> é permitido, logo:

<math>\mathcal{Z} = \prod_{i}(1+ze^{-\beta E_i})^{-1}</math>

O potencial termodinâmico é então:

<math>-PV(z,\beta) = \Omega(z,\beta)= \frac{1}{\beta}\sum_i\ln (1-ze^{-\beta E_i})</math>

Se o gás possuir apenas graus de liberdade translacionais em <math>d</math> dimensões (os demais casos podem ser tratados de forma análoga):

<math>\Omega(z,\beta,V^{(d)}) = \frac{V^{(d)}}{\beta}\frac{2\pi^{d/2}}{h^d\Gamma(d/2)}\int_0^{\infty}p^{d-1}\ln (1-ze^{-\beta p^2/2m})-\frac{1}{\beta}\text{Li}_1(z)</math>

onde <math>\Gamma</math> é a função gama, <math>\text{Li}_s(z)</math> é a função polilogarítma e <math>V^{(d)}</math> é o volume d-dimensional que o gás ocupa.

<math>\Omega(z,\beta,V^{(d)}) = -\frac{V^{(d)}(d-1)\pi^{d/2}}{2mh^d}\left(\frac{2m}{\beta}\right)^{d/2+1}\text{Li}_{d/2+1}(z)-\frac{1}{\beta}\text{Li}_1(z)</math>

Note que a função polilogarítma só está definida para <math>z</math> reais menores ou iguais a 1. O segundo termo que já estava presente na expressão anterior é a contribuição de momento zero, ou seja, do estado de menor energia.

Condensação de Bose-Einstein

O gás de bósons é o sistema mais simples que apresenta o fenômeno de condensação de Bose-Einstein. Para ver esse efeito, escrevemos o número médio de partículas:

<math>N(z,\beta,V^{(d)}) = -\beta z \frac{\partial \Omega}{\partial z} = \frac{V^{(d)}(d-1)\pi^{d/2}}{h^d}\left(\frac{2m}{\beta}\right)^{d/2}\text{Li}_{d/2}(z)+\text{Li}_0(z)</math>

O maior valor da função polilogarítma acontece em <math>z=1</math> quando o número de partículas em estados excitados é:

<math>N^{(p>0)}(z,\beta,V^{(d)}) = \frac{V^{(d)}(d-1)\pi^{d/2}}{h^d}\left(\frac{2m}{\beta}\right)^{d/2}\zeta(d)</math>

Perceba que para <math>d>2</math> isso é um número finito que é atingido numa certa temperatura <math>T_0</math>. Todas as demais

<math>N^{(p=0)}=N\left[1-\left(\frac{T}{T_0}\right)^{d/2}\right]</math>

partículas deverão estar no estado fundamental, não importando quantas sejam (contanto que a aproximação de gás continue valendo).