Resumos
Os condensados de Bose-Einstein foram previstos teoricamente nos anos de 1924–25. Contudo, sua observação experimental só foi obtida em 1995, laureando seus descobridores com o prêmio Nobel de 2001. No ano de 2015 comemora-se o aniversário de 90 anos de sua previsão teórica e 20 anos de sua realização experimental. Os condensados de Bose-Einstein têm sido usados nos mais variados contextos, que vão da matéria condensada à astrofísica. Este trabalho será dedicado a uma interpretação intuitiva para a condensação. Discutiremos também como modelar armadilhas modernas que aprisionam partículas de um gás de bósons fracamente interagentes para formar os condensados.
condensados de Bose-Einstein; mecânica estatística; função partição
The Bose-Einstein condensates were theoretically predicted in 1924–25. However, the experimental observation was reported in 1995, laureating their discoverers with the 2001 Nobel prize. The year of 2015 celebrates the ninetieth anniversary of the theoretical prediction and the twentieth anniversary of its experimental realization. The Bose-Einstein condensates have been used in different contexts, from condensed matter to astrophysics. This work is dedicated to provide a naïve interpretation of the condensation process. We will also model modern traps which imprison the particles of a weak interacting gas of bosons to form the condensates.
Bose-Einstein condensates; statistical mechanics; partition function
Introdução
- a)
Eles ocorrem a temperaturas próximas ao zero absoluto, num regime propício para efeitos que rompem com a termodinâmica usual, próximo ao limite da terceira lei.
- b)
Seguindo esta linha, próximo ao zero absoluto os bósons estão na iminência de violar o princípio da incerteza: uma vez que não há mais agitação térmica conseguiríamos detectar posição e momento do condensado. Isto significa testar limites de validade para a mecânica quântica.
- c)
A baixas temperaturas, um gás de bósons aprisionado é bem descrito pela equação de Gross-Pitaeviskii [6[6] A.G. de Sousa e V.S. Bagnato, Revista Brasileira de Ensino de Física 26, 43 (2004).]. Numa certa aproximação, os condensados têm sua evolução temporal descrita por uma equação do tipo de Dirac. Dessa forma, os CBE foram utilizados para se observar experimentalmente o chamado Zitterbewegung [7[7] L.J. LeBlanc, M.C. Beeler, K. Jiménez-García, A.R. Perry, S. Sugawa, R.A. Williams and I.B. Spielman, New J. Phys. 15, 073011 (2013).–9[9] V. Achilleos, D.J. Frantzeskakis and P.G. Kevrekidis, Phys. Rev. A 89, 033636 (2014).].2 2 Zitterbewegung foi um movimento oscilatório dos elétrons previsto teoricamente por Schrödinger ao analisar a dinâmica dos operadores de posição oriundos da equação de Dirac. Além do movimento retilíneo e uniforme, os elétrons apresentariam também este movimento oscilatório com amplitude da ordem do comprimento de onda Compton do elétron (∼ 10−12 m) e frequência também da ordem de frequência Compton do elétron (∼ 1020 s−1). Estes valores poderiam impedir a detecção experimental do Zitterbewegung. Para mais detalhes, veja a Ref. [10].
- d)
O fenômeno da localização de Anderson [11[11] P.W. Anderson, Phys. Rev. 109, 1492 (1958).] tem sido estudo de intensa pesquisa, veja por exemplo [12[12] T. Micklitz, C.A. Müller, A. Altland, Phys. Rev. Lett. 112, 110602 (2014).]. Foram utilizados átomos ultra-frios de um CBE para observação deste efeito [13[13] Giacomo Roati, Chiara D'Errico, Leonardo Fallani, Marco Fattori, Chiara Fort, Matteo Zaccanti, Giovanni Modugno, Michele Modugno, Nature 453, 895 (2008).].
- e)
Recentemente foi sugerido uma possível condensação em objetos compactos da astrofísica, em particular, nas estrelas de nêutrons [14[14] Christine Gruber and Axel Pelster, arXiv: 1403.3812v1 [gr-qc].].
Há também razões pedagógicas para se estudar os condensados: estudantes ainda na graduação dos bacharelados em ciências exatas têm oportunidade de entrar em contato com assuntos na fronteira do conhecimento, como os itens (c), (d) e (e) acima. Os únicos prérequisitos necessários para se analisar a condensação de Bose-Einstein são fundamentos de estatística quântica. O estudo de tópicos modernos é motivante para alunos em geral, já que ao longo de um bacharelado em física, por exemplo, são fornecidos os pilares básicos para carreiras científicas incluindo mecânicas clássica e quântica, eletromagnetismo e física estatística, com pouca ênfase em física contemporânea.
Nosso objetivo neste trabalho será mostrar como ocorre a condensação. O intuito aqui não é substituir a literatura padrão do assunto [3[3] Luca Peliti, Statistical Mechanics in a Nutshell (Princeton University Press, Princeton, 2011)., 15[15] R.K. Pathria, Statistical Mechanics (Butterworth Heinemann, Oxford, 1996).–18[18] M. Kardar, Statistical Physics of Particles (Cambridge University Press, Cambridge, 2007).] e sim fornecer uma interpretação intuitiva e introdutória para os CBE. Apresentaremos também uma análise teórica mais realista para armadilhas modernas onde estuda-se o gás de bósons fracamente interagentes.
2.Descrição de um sistema com N partículas com energia discreta
Vamos assumir daqui para frente que o leitor tenha certa familiaridade com fundamentos de física estatística quântica. Do contrário, direcionamos o leitor para a Seção IV (ou Conclusão), onde discutimos a sentido físico dos CBE. Os cálculos desta Seção e das Seções III, IV e V serão apresentados por completude.
Consideremos um gás de N partículas não-interagentes de massa m em uma caixa cúbica com aresta L (seu volume é dado por V = L3). Da mecânica quântica [15[15] R.K. Pathria, Statistical Mechanics (Butterworth Heinemann, Oxford, 1996).], sabemos que a energia do sistema é discreta,
onde Ĥ é o hamiltoniano do sistema, εi corresponde aos valores de energia e n̂i é o operador número de partículas. n̂i atua num espaço de Hilbert ℋ “dividido” em sub-espaços de número de partículas, ℋ = ℋ1 ⊕ ℋ2 ⊕ ⋯ de forma que cada vetor de estado tem a seguinte estrutura, Por exemplo, |ψ〉 = |0,2,5,0,0...〉 representa duas partículas em ℋ2, cinco partículas em ℋ3. Assim, n^1|ψ〉 = 0, n^2|ψ〉 = 2|ψ〉, etc. Esta construção é válida para bósons. Para férmions, pelo princípio da exclusão de Pauli só podemos ter zero ou uma partícula em cada ℋi. Neste caso, o espaço de Hilbert ℋ, como construído acima, é conhecido na literatura como espaço de Fock [18[18] M. Kardar, Statistical Physics of Particles (Cambridge University Press, Cambridge, 2007).]. Esta é a construção que permite tratar número de partículas como um observável. 3.Fundamentos de estatística quântica aplicados ao gás de bósons
Em estatística quântica, o valor esperado de um observável A é dado por
onde  é o operador que representa o observável, Tr denota o traço do operador entre parêntesis e ρ corresponde ao operador densidade.3 3 Descrevemos o operador densidade com mais detalhes no Apêndice A. Como trataremos uma condensação, é natural utilizar o ensemble grande canônico, uma vez que há variação do número de partículas passando do estado gasoso para o condensado. Assim, β = (kBT)−1 com kB a constante de Boltzmann, T a temperatura do sistema e μ é o potencial químico. Por fim, denota o operador número total de partículas. Para o caso em questão, devido às Eqs. (1) e (5), temos Valores esperados podem ser calculados diretamente a partir da função partição,4 4 É costume também usar a letra Z para a função de partição pois no fundo estamos calculando a soma de estados do sistema, que em alemão traduz-se por Zustandsumme [3, 15]. Ao longo do texto, usaremos a letra Ξ por tratar-se do ensemble grande canônico. pois, por exemplo, onde usamos α ≡ −μβ.Usando a base |ni〉 de ℋ de autovetores de N̂ (e também de Ĥ), temos (os detalhes estão no Apêndice B Apêndice B Neste apêndice detalharemos os passos para a obtenção da função partição (9). Sabendo que o traço de um operador é independente da base do espaço que se utiliza para calculá-lo [23], para se obter a Eq. (9) usaremos {|n̂i〉}, a base de ℋ formada por autovetores comuns a N̂ e Ĥ. Temos então, (64)Ξ=Tr(e−β∑iεin^i−α∑in^i)=Tr(e−∑i(βεi+α)n^i)=〈…,n2,n1|e−∑i(βεi+α)n^i|n1,n2,…〉.A expressão acima é composto pelo produto de fatores da forma, (65)〈ni|e−∑ni(βεi+α)n^i|ni〉.Atuando com o operador n̂i, chegamos à expressão (66)Ξ=∑n1,n2,…∏ie−(βεi+α)ni=(∑n1e−(βε1+α)n1)(∑n2e−(βε2+α)n2)⋯Finalmente, como os índices ni, i = 0,1,2, ... assumem os mesmos valores em cada um das somas acima, podemos escrever a função de partição como o seguinte produto (67)Ξ=∏i∑n=0+∞e−(βεi+α)n. )
Tomando o logaritmo de Ξ (veja a Eq. (8)) e levando em conta que a série (9) é a série geométrica, temos, No estado fundamental, quando T → 0, os bósons têm maior probabilidade de serem encontrados em seu menor estado de energia. Assim, separamos do somatório (10) o termo com i = 0, quando ε0 = 0, Para entender a expressão acima, poderíamos fazer uma analogia com uma mistura que contenha duas fases da mesma substância, sendo cada uma das fases descritas pelos dois termos no lado direito da Eq. (11). Por exemplo, consideremos um recipiente contendo água nas formas líquida e gasosa. A fase líquida (condensada) corresponderia ao primeiro termo e o vapor ao segundo termo.Nos interessa analisar o comportamento do termo correspondente à fração com εi ≠ 0 do gás de bósons. Para estimar seu valor, vamos tomar o limite termodinâmico, que consiste em assumir que N, V → +∞ com a fração N/V constante. Deste modo pode-se aproximar o somatório em (11) por uma integral. Precisamos neste caso determinar a densidade de estados em função da energia do sistema. Consideramos o gás como partículas livres em uma caixa cúbica de aresta L. Usando condições de contorno periódicas, as funções de onda têm a forma
com nx, ny e nz inteiros não negativos. A energia por sua vez é dada por onde usamos p⃗ = ħk⃗. A expressão (13) fornece a estrutura geométrica do espaço de estados: fixado n⃗ = (nx, ny, nz), temos um esfera de raio . Logo, o volume Ω no espaço de estados com energia menor ou igual a E é dado por Por fim, a densidade de estados toma a forma Portanto, tende a Mas a série para o logaritmo é conhecida, Então, Assim, a integral em (16) é reescrita como Integrando e retornando à expressão (11), chegamos a onde as funções partições ficam divididas entre condensado (ΞC) e gás (ΞG) Quando T → 0, 〈n̂0〉 → N, isto é, todo o gás encontrase no estado de energia ε = 0. Assim, onde α0 = α(T → 0). Por outro lado, a partir da Eq. (20) vemos que deve existir uma temperatura crítica TC quando começa a condensação (as partículas do gás aglomeram-se e passam a ser descritas por ΞC no lugar de ΞG. Nas palavras do próprio Einstein [4[4] Sílvio R. Dahmen, Revista Brasileira de Ensino de Física 27, 271 (2005); 27, 283 (2005).]: “Uma separação surge; parte (do gás) se condensa e o resto permanece como um gás ideal saturado” [1[1] A. Einstein, Sitzungsberichte der Preussischen Akademie der Wissenschaften I, 3–14 (1925). A tradução deste trabalho foi feita na Revista Brasileira de Ensino de Física 27, 113 (2005) pelo Prof. Sílvio R. Dahmen.]). Como quando a temperatura diminui, suponhamos que uma pequena parcela p do gás se condensou quando T ≈ TC, digamos 〈n̂〉C ≈ pN, com , isto é : pois N >> 1. Usando que α(TC) ≈ 0 e calculando 〈n̂〉G para TC, A série que aparece na última expressão é a função ζ de Riemann , Assim, A expressão (27) foi obtida com algumas aproximações. Uma fórmula mais acurada que fornece TC, conhecida como temperatura de Bose-Einstein pode ser vista na Ref. [17[17] Sílvio R.A. Salinas, Introdução à Física Estatística (Editora da Universidade de São Paulo, São Paulo, 2005).]. 4.Interpretando o condensado
Vamos focar nossa atenção na expressão (27), obtida na Seção anterior. Primeiro, estimamos TC,
Para valores típicos, 〈n̂〉G ≈ N ∼ 1024, V ∼ 125, cm3 (cubo com 5 cm de aresta), m ∼ 1, 8 × 10−27 kg (como só queremos uma estimativa, tomamos a massa de um próton) e sabendo que ζ(3/2) ≈ 2,6, temos TC ∼ 6 K. Vamos reescrever a expressão que obtivemos para TC da seguinte maneira, onde Com esta fórmula chegamos ao objetivo central deste trabalho, a interpretação intuitiva do condensado. O fator tem unidades de momento linear. Logo, tem unidades de distância: interpretamos este fator como a aresta de um cubo com tamanho da ordem do comprimento de onda térmico de de Broglie de uma partícula, onde elas se agrupam. Quando T diminui (T < TC), λ aumenta. Para manter o lado esquerdo da expressão (29) constante, o valor esperado do número de partículas do gás 〈n̂〉G diminui, isto é, o gás condensa-se. 5.Formulação mais realista para os CBE
Anteriormente descrevemos um gás de bósons em uma caixa. A condição de contorno dada pela caixa não é uma boa modelagem para experiências modernas de condensação. Assim, nosso objetivo agora é fazer uma descrição mais realista dos CBE. Inicialmente, vamos discutir de maneira sucinta um método fortemente utilizado para aprisionamento de átomos [19[19] A.H. Iavaronni, E.A.L. Henn, E.R. F. Ramos, J.A. Seman, T. Amthor e V.S. Bagnato, Revista Brasileira de Ensino de Física 29, 209 (2007).], as armadilhas magnéticas. Maiores detalhes podem ser encontrados em [20[20] C.J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge, Cambridge University Press, 2008).].5 5 Os autores agradecem ao referee pela sugestão desta referência. Na presença de campo magnético B⃗, um átomo com momento magnético μ⃗ (proporcional ao seu spin) ganha energia potencial
Em teoria da perturbação, a presença de V no Hamiltoniano gera um split nos níveis de energia do átomo. Este é o famoso efeito Zeeman. De acordo com o ângulo θ entre μ⃗ e B⃗, podemos ter átomos com energia crescente ou decrescente com o módulo de B⃗ já que V = −|μ⃗‖B⃗| cosθ. O grupo de átomos cuja energia aumenta com o módulo de B⃗ são aqueles que serão aprisionados. De fato, pela tendência de um sistema físico de minimizar sua energia, os átomos serão atraídos para um ponto de mínimo de |B⃗|. O contrário também seria possível: aprisionar átomos cuja energia decresce com |B⃗| já que eles buscariam um máximo de |B⃗|. Contudo descartamos esta última possibilidade uma vez que só é possível gerar mínimos de campo na ausência de correntes [19[19] A.H. Iavaronni, E.A.L. Henn, E.R. F. Ramos, J.A. Seman, T. Amthor e V.S. Bagnato, Revista Brasileira de Ensino de Física 29, 209 (2007).,20[20] C.J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge, Cambridge University Press, 2008).]. O confinamento de átomos nestas armadilhas é bem aproximado por um oscilador harmônico (isotrópico) [20[20] C.J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge, Cambridge University Press, 2008).,21[21] Franco Dalfovo, Stefano Giorgini, Lev P. Pitaevskii and Sandro Stringari, Reviews of Modern Physics 71, 463 (1999).], cuja energia potencial de um átomo de massa m aprisionado é dada por, onde x, y e z são as coordenadas de posição do átomo e a subtração do fator atribui energia 0 para o estado fundamental. Seguiremos os mesmos passos da Seção anterior para mostrar que a temperatura crítica quando começa a condensação é da ordem de mili Kelvins. Como veremos, a estrutura da função de onda do oscilador no estado fundamental nos permite prever que a condensação ocorre numa região espacial bem delimitada. Um outro ponto que discutiremos nesta Seção, por razões pedagógicas, é a contagem de estados para um sistema quântico. Como temos uma estrutura discreta na mecânica quântica, este procedimento, central em Física estatística,6 6 A contagem de estados e sua relação com a entropia é tão fundamental que a expressão S = k log W está gravada no túmulo de Boltzmann no Zentralfriehof, em Viena. Veja a foto na Ref. [17]. é feito de maneira enumerável.Para obter a função partição grande canônica de um gás de bósons na armadilha, precisamos determinar a densidade de estados em função da energia. Os autovalores de energia do oscilador isotrópico são dados por
uma vez que a energia é zero no estado fundamental. nx, ny e nz são inteiros não negativos. No espaço dos n’s, a Eq. (33) representa um plano para cada valor de n, como pode ser visto na Fig. 1.Aumentando o valor de n, teremos o seguinte número de sítios,
e assim sucessivamente, de forma que7 7 Observemos que para n = 3, além dos 9 sítios nas arestas do tetraedro, temos ainda um deles pertencente ao interior do plano. Isto pode ser mostrado notando que a diagonal de um cubo com lado unitário tem o mesmo tamanho que a altura do tetraedro com base no plano definido por n = 3 e vértice na origem n = 0. Na física estatística clássica, a contagem de microestados acessíveis a determinado sistema é feita de forma “contínua”: efetivamente calculamos um volume no espaço de fase. Já na estatística quântica, estados acessíveis são literalmente enumerados uma vez que temos uma estrutura discreta dos estados, fixada a energia. A densidade de estados é então dada por onde usamos a Eq. (33) e ΔE = hν corresponde à variação de energia do sistema quando passamos de um estado n para n + 1. A função partição é dada por cujo logaritmo fornece, Assim como antes o primeiro termo corresponde ao condensado e o segundo, ao gás, que pode ser aproximado pela integral utilizando ω(E) (já expandimos o logaritmo em série) Supondo ainda que a energia das partículas do gás seja maior que a energia do estado fundamental E >> hν, fazemos . Assim, Como antes, o índice G refere-se ao gás. Desta maneira, Quando começa a condensação em T = TC, α(TC) ≈ 0, (veja a Eq. (24)) e 〈n̂〉G, = N isto é, ainda não há partículas condensadas ou, Para estimar TC, consideramos que os átomos estejam no estado fundamental, com função de onda em forma de gaussiana [22[22] David J. Griffiths, Introduction to Quantum Mechanics (Prentice Hall, New Jersey, 1995).], R corresponde a uma estimativa do raio de uma região esférica onde ocorre a condensação. Neste caso, Assim, . Por outro lado, A temperatura de condensação então é dada por, Supondo que a condensação ocorra em uma região espacial da ordem R ≈ 1 μm, e tomando dados fornecidos pelos experimentos de 1995 [5[5] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Science 269, 198 (1995).]: população da ordem de N ≈ 107 átomos de rubídio com massa m = 85 · 1, 66 · 10−27 kg, encontramos TC ∼ 3 μK. Para o Rubídio, TC = 170 nK. Nossa estimativa foi feita para um gás não-interagente e ainda assim o resultado é razoável. 6.Conclusão
Neste artigo apresentamos uma interpretação intuitiva para os condensados de Bose-Einstein. Usando fundamentos de estatística quântica, obtivemos a expressão (29): as partículas do gás de bósons se alocam dentro de um cubo com aresta da ordem do comprimento de onda térmico de de Broglie de uma partícula formadora do gás. Uma outra maneira para se interpretar a Eq. (29) é a seguinte: fornece a densidade do gás quando começa a condensação. é o volume que aproximadamente cada partícula vai ocupar já que dλ3 ≈ 1.
As experiências modernas que envolvem armadilhas que formam os CBE são bem descritas pelo potencial de um oscilador harmônico isotrópico. Assim, descrevemos também como são formados os condensados. Com dados da literatura, obtivemos, por exemplo, a temperatura em que ocorre o condensado, da ordem de μK. Esperamos que este trabalho possa ser utilizado tanto por alunos de graduação e pós-graduação ao cursar a disciplina de mecânica estatística quanto por professores do ensino médio que queriam ter uma visão global sobre os famosos condensados de Bose-Einstein.
Apêndice A
Vamos descrever neste apêndice a estrutura das expressões (3) e (4), fornecendo uma interpretação para o operador densidade. Consideremos um observável A representado por um operador auto-adjunto Â. Usando a decomposição espectral de  [24[24] R.C. de Cerqueira Leite e A.R. Britto de Castro, Física do Estado Sólido (Editora Edgard Blücher, São Paulo, 1978). ], temos
onde os ai’s são os autovalores associados aos autovetores ortonormais |ai〉. Cada 𝒫ai é o projetor ao autosubespaço gerado por |ai〉. Se um sistema está em um estado representado por um vetor |ψ〉, então o valor esperado 〈Â〉 de uma medida de A é dado pela média dos possíveis valores da medida, ponderado por pesos, que correspondem ao quadrado das componentes de |ψ〉 em cada autosubespaço Como 〈ai|ai〉 = 1, temos Em geral, o estado de um sistema é representado não só por um vetor, mas por uma sobreposição de vetores de estado com respectivas probabilidades de entrarem na mistura, Assim, o valor esperado de A quando o sistema encontra-se no estado dado pela Eq. (51) é Seja {|φi〉} uma base do espaço onde  atua. Inserindo a relação de completeza [23[23] C. Cohen-Tannoudji, B. Diu, F. Laloë, Quantum Mechanics (John Wiley and Sons, New York, 1977), v. I.] na Eq. (52), encontramos Reagrupando os termos Definindo ρ = ∑lpl|ψl〉 〈ψl| e lembrando que a representação matricial de um operador é feita pelo seu sanduíche com bras e kets, (Â)ij = 〈φi|Â|φj〉, a Eq. (55) toma a forma já que o traço de um operador é a soma dos elementos da diagonal principal da matriz que o representa. Deixemos o formalismo da estatística quântica de lado por um momento. No caso clássico, o valor esperado de uma variável dinâmica definida no espaço de fase B = B(qi, pi); i = 1, ..., 3N é dada por onde dΓ = dq1 ⋯ dq3N dp1 ⋯ dp3N e ϱ corresponde à densidade de probabilidade de encontrarmos o sistema em questão em determinado microestado definido por um ponto (qi, pi) no espaço de fase e tem a forma dos pesos de Boltzmann Para o caso particular do ensemble grande canônico, em que o sistema tem Hamiltoniana H e troca energia mediante a troca de partículas, Fazendo as identificações abaixo a partir das expressões (56) e (57), é sugestivo chamarmos ρ na Eq. (56) de operador densidade. Vamos construí-lo para o ensemble grande canônico, utilizado ao longo deste trabalho. Tomemos uma base {|φl〉} do espaço de Fock discutido na Seção II e utilizando os pesos pl = Ce−(ε l−μnl)) (C é uma constante determinada por normalização) Utilizando novamente o fato que ∑l|φl〉 〈φl| = 𝕀 e normalizando a constante C chegamos finalmente aApêndice B
Neste apêndice detalharemos os passos para a obtenção da função partição (9). Sabendo que o traço de um operador é independente da base do espaço que se utiliza para calculá-lo [23[23] C. Cohen-Tannoudji, B. Diu, F. Laloë, Quantum Mechanics (John Wiley and Sons, New York, 1977), v. I.], para se obter a Eq. (9) usaremos {|n̂i〉}, a base de ℋ formada por autovetores comuns a N̂ e Ĥ. Temos então,
A expressão acima é composto pelo produto de fatores da forma, Atuando com o operador n̂i, chegamos à expressão Finalmente, como os índices ni, i = 0,1,2, ... assumem os mesmos valores em cada um das somas acima, podemos escrever a função de partição como o seguinte produto-
2
Zitterbewegung foi um movimento oscilatório dos elétrons previsto teoricamente por Schrödinger ao analisar a dinâmica dos operadores de posição oriundos da equação de Dirac. Além do movimento retilíneo e uniforme, os elétrons apresentariam também este movimento oscilatório com amplitude da ordem do comprimento de onda Compton do elétron (∼ 10−12 m) e frequência também da ordem de frequência Compton do elétron (∼ 1020 s−1). Estes valores poderiam impedir a detecção experimental do Zitterbewegung. Para mais detalhes, veja a Ref. [10[10] A. Deriglazov, B. Rizzuti and G. Zamudio, Spinning Particles: Possibility of Space-Time Interpretation for the Inner Space of Spin (Lap Lambert Academic Publishing, Saarbrücken, 2012).].
-
3
Descrevemos o operador densidade com mais detalhes no Apêndice A Apêndice A Vamos descrever neste apêndice a estrutura das expressões (3) e (4), fornecendo uma interpretação para o operador densidade. Consideremos um observável A representado por um operador auto-adjunto Â. Usando a decomposição espectral de  [24], temos (48)A^=∑iai𝒫ai; 𝒫ai=|ai〉〈ai|,onde os ai’s são os autovalores associados aos autovetores ortonormais |ai〉. Cada 𝒫ai é o projetor ao autosubespaço gerado por |ai〉. Se um sistema está em um estado representado por um vetor |ψ〉, então o valor esperado 〈Â〉 de uma medida de A é dado pela média dos possíveis valores da medida, ponderado por pesos, que correspondem ao quadrado das componentes de |ψ〉 em cada autosubespaço (49)〈A^〉=∑iai‖𝒫ai|ψ‖2=∑iai〈ψ|ai〉〈ai|ai〉〈ai|ψ〉.Como 〈ai|ai〉 = 1, temos (50)〈A^〉=〈ψ|(∑iai|ai〉〈ai|)|ψ〉=〈ψ|A^|ψ〉.Em geral, o estado de um sistema é representado não só por um vetor, mas por uma sobreposição de vetores de estado com respectivas probabilidades de entrarem na mistura, (51)|ψ〉=∑lpl|ψl〉; pl≥0,∑lpl=1.Assim, o valor esperado de A quando o sistema encontra-se no estado dado pela Eq. (51) é (52)〈A^〉=∑lpl〈ψl|A^|ψl〉.Seja {|φi〉} uma base do espaço onde  atua. Inserindo a relação de completeza [23] (53)𝕀=∑i|φi〉 〈φi|na Eq. (52), encontramos (54)〈A^〉∑lpl〈ψl|(∑i|φi〉 〈φi|)A^|ψl〉.Reagrupando os termos (55)〈A^〉∑i〈φl|(A^∑lpl|ψl〉〈ψl|)|φi〉.Definindo ρ = ∑lpl|ψl〉 〈ψl| e lembrando que a representação matricial de um operador é feita pelo seu sanduíche com bras e kets, (Â)ij = 〈φi|Â|φj〉, a Eq. (55) toma a forma (56)〈A^〉=Tr(A^ρ),já que o traço de um operador é a soma dos elementos da diagonal principal da matriz que o representa. Deixemos o formalismo da estatística quântica de lado por um momento. No caso clássico, o valor esperado de uma variável dinâmica definida no espaço de fase B = B(qi, pi); i = 1, ..., 3N é dada por (57)〈B〉=1N𝒫h3N∫B ϱdΓ1N!h3N∫ ϱdΓ,onde dΓ = dq1 ⋯ dq3N dp1 ⋯ dp3N e ϱ corresponde à densidade de probabilidade de encontrarmos o sistema em questão em determinado microestado definido por um ponto (qi, pi) no espaço de fase e tem a forma dos pesos de Boltzmann (58)ϱ=e−βE.Para o caso particular do ensemble grande canônico, em que o sistema tem Hamiltoniana H e troca energia mediante a troca de partículas, (59)E=H−μN.Fazendo as identificações abaixo a partir das expressões (56) e (57), (60)Estatística qua^nticaEstatísica cla´ssicaA^↔B(qi,pi)Tr↔1N!h3N∫ dΓρ↔ϱé sugestivo chamarmos ρ na Eq. (56) de operador densidade. Vamos construí-lo para o ensemble grande canônico, utilizado ao longo deste trabalho. Tomemos uma base {|φl〉} do espaço de Fock discutido na Seção II e utilizando os pesos pl = Ce−(ε l−μnl)) (C é uma constante determinada por normalização) (61)ρ=∑lpl|φl〉 〈φl|=∑lCe−(εl−μnl)|φl〉 〈φl|=Ce−β(H^−μN^)∑l|φl〉 〈φl|.Utilizando novamente o fato que ∑l|φl〉 〈φl| = 𝕀 e normalizando a constante C (62)C=1Tr(e−β(H^−μN^)),chegamos finalmente a (63)ρ=e−β(H^−μN^)Tr(e−β(H^−μN^)). .
-
4
É costume também usar a letra Z para a função de partição pois no fundo estamos calculando a soma de estados do sistema, que em alemão traduz-se por Zustandsumme [3[3] Luca Peliti, Statistical Mechanics in a Nutshell (Princeton University Press, Princeton, 2011)., 15[15] R.K. Pathria, Statistical Mechanics (Butterworth Heinemann, Oxford, 1996).]. Ao longo do texto, usaremos a letra Ξ por tratar-se do ensemble grande canônico.
-
5
Os autores agradecem ao referee pela sugestão desta referência.
-
6
A contagem de estados e sua relação com a entropia é tão fundamental que a expressão S = k log W está gravada no túmulo de Boltzmann no Zentralfriehof, em Viena. Veja a foto na Ref. [17[17] Sílvio R.A. Salinas, Introdução à Física Estatística (Editora da Universidade de São Paulo, São Paulo, 2005).].
-
7
Observemos que para n = 3, além dos 9 sítios nas arestas do tetraedro, temos ainda um deles pertencente ao interior do plano. Isto pode ser mostrado notando que a diagonal de um cubo com lado unitário tem o mesmo tamanho que a altura do tetraedro com base no plano definido por n = 3 e vértice na origem n = 0.
References
-
[1]A. Einstein, Sitzungsberichte der Preussischen Akademie der Wissenschaften I, 3–14 (1925). A tradução deste trabalho foi feita na Revista Brasileira de Ensino de Física 27, 113 (2005) pelo Prof. Sílvio R. Dahmen.
-
[2]S.N. Bose, Zeit. Phys. 26, 178 (1924).
-
[3]Luca Peliti, Statistical Mechanics in a Nutshell (Princeton University Press, Princeton, 2011).
-
[4]Sílvio R. Dahmen, Revista Brasileira de Ensino de Física 27, 271 (2005); 27, 283 (2005).
-
[5]M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Science 269, 198 (1995).
-
[6]A.G. de Sousa e V.S. Bagnato, Revista Brasileira de Ensino de Física 26, 43 (2004).
-
[7]L.J. LeBlanc, M.C. Beeler, K. Jiménez-García, A.R. Perry, S. Sugawa, R.A. Williams and I.B. Spielman, New J. Phys. 15, 073011 (2013).
-
[8]Chunlei Qu, Chris Hamner, Ming Gong, Chuanwei Zhang and Peter Engels, Phys. Rev. A 88, 021604 (2013).
-
[9]V. Achilleos, D.J. Frantzeskakis and P.G. Kevrekidis, Phys. Rev. A 89, 033636 (2014).
-
[10]A. Deriglazov, B. Rizzuti and G. Zamudio, Spinning Particles: Possibility of Space-Time Interpretation for the Inner Space of Spin (Lap Lambert Academic Publishing, Saarbrücken, 2012).
-
[11]P.W. Anderson, Phys. Rev. 109, 1492 (1958).
-
[12]T. Micklitz, C.A. Müller, A. Altland, Phys. Rev. Lett. 112, 110602 (2014).
-
[13]Giacomo Roati, Chiara D'Errico, Leonardo Fallani, Marco Fattori, Chiara Fort, Matteo Zaccanti, Giovanni Modugno, Michele Modugno, Nature 453, 895 (2008).
-
[14]Christine Gruber and Axel Pelster, arXiv: 1403.3812v1 [gr-qc].
-
[15]R.K. Pathria, Statistical Mechanics (Butterworth Heinemann, Oxford, 1996).
-
[16]F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill Book Company, New York, 1965).
-
[17]Sílvio R.A. Salinas, Introdução à Física Estatística (Editora da Universidade de São Paulo, São Paulo, 2005).
-
[18]M. Kardar, Statistical Physics of Particles (Cambridge University Press, Cambridge, 2007).
-
[19]A.H. Iavaronni, E.A.L. Henn, E.R. F. Ramos, J.A. Seman, T. Amthor e V.S. Bagnato, Revista Brasileira de Ensino de Física 29, 209 (2007).
-
[20]C.J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge, Cambridge University Press, 2008).
-
[21]Franco Dalfovo, Stefano Giorgini, Lev P. Pitaevskii and Sandro Stringari, Reviews of Modern Physics 71, 463 (1999).
-
[22]David J. Griffiths, Introduction to Quantum Mechanics (Prentice Hall, New Jersey, 1995).
-
[23]C. Cohen-Tannoudji, B. Diu, F. Laloë, Quantum Mechanics (John Wiley and Sons, New York, 1977), v. I.
-
[24]R.C. de Cerqueira Leite e A.R. Britto de Castro, Física do Estado Sólido (Editora Edgard Blücher, São Paulo, 1978).
Figure
Datas de Publicação
-
Publicação nesta coleção
30 Mar 2015 -
Data do Fascículo
Jan-Mar 2015
Histórico
-
Recebido
05 Nov 2014 -
Aceito
19 Dez 2014