Expansion of a free Fermi gas released from an isotropic trapping potential

Resumo Técnico: Expansão de um Gás de Fermi Livre Liberado de um Potencial de Armadilha Isotrópica

Enunciado do Problema
Este trabalho investiga a dinâmica fora do equilíbrio de um sistema de $N$ férmions não interagentes, sem spin (ou com degenerescência de spin), inicialmente preparados no estado fundamental de muitos corpos de um potencial de armadilha isotrópico $V(\mathbf{x})$. O estudo foca no regime que segue a liberação instantânea do potencial em $t=0$. Especificamente, os autores analisam a expansão balística da nuvem fermiônica em tempos tardios, visando derivar os perfis espaciais e de momento completos do gás no limite de grande número de férmions $N$. O objetivo principal é determinar como a forma do potencial de armadilha inicial (caracterizada por um expoente de lei de potência $a$) e a dimensão espacial $d$ influenciam o perfil de densidade assintótico e o comportamento da borda da nuvem em expansão.

Metodologia
Os autores empregam uma combinação de soluções analíticas exatas para casos específicos e aproximações semiclássicas (WKB) para o caso geral de grande $N$.

1. Soluções Exatas (Potenciais Harmônico e de Poço Quadrado):

   • Para um potencial harmônico ($V \propto x^2$), os autores utilizam a dependência temporal exata da função de onda de muitos corpos, que é governada pela equação de Ermakov-Pinney. Eles demonstram que as posições reescalonadas dos férmions seguem a distribuição de autovalores de um conjunto de matrizes unitárias gaussianas (GUE).
   • Para um potencial de poço quadrado unidimensional, funções de onda exatas (envolvendo funções de Bessel em 2D e funções trigonométricas em 1D) são utilizadas para calcular a distribuição de momento $f_N(p)$ para $N$ finito, permitindo uma análise de comportamentos de borda e correções de tamanho finito.

2. Aproximação Semiclássica (Potenciais Gerais):

   • Para potenciais arbitrários de lei de potência isotrópica $V(r) \propto r^a$ em qualquer dimensão $d$, os autores aplicam a aproximação de fase estacionária à função de onda evoluída no tempo. Isso estabelece que, em tempos tardios, a posição $\mathbf{x}$ de um férmion está relacionada ao seu momento inicial $\mathbf{p}$ por escalonamento balístico: $\mathbf{x} \approx \mathbf{p}t$.
   • A distribuição de momento $f_N(p)$ é derivada pela soma do módulo ao quadrado das transformadas de Fourier dos estados próprios ocupados de partícula única. No limite de grande $N$, essa soma é convertida em uma integral usando a regra de quantização semiclássica e a aproximação de Thomas-Fermi para a densidade de estados.
   • A distribuição de momento microcanônica em um $E$ fixo é calculada com base no volume de fase clássico, que é então integrado sobre os estados de energia ocupados até a energia de Fermi $E_F$.

Principais Contribuições e Resultados

1. Perfis de Densidade Generalizados:
   O artigo deriva expressões explícitas para o perfil de densidade espacial tardia $\rho_N(r; t)$ e a distribuição de momento $f_N(p)$ para potenciais arbitrários de lei de potência $V(r) \propto r^a$ em $d$ dimensões.

   • Distribuição de Momento: Para $p < p_F$, a distribuição escala como:
     $$f_N(p) \sim \frac{N}{p_F^d} \left( 1 - \frac{p^2}{p_F^2} \right)^{d/a}$$
   • Perfil Espacial: Consequentemente, o perfil espacial tardio da nuvem (para $r < p_F t$) é:
     $$\rho_N(r; t) \sim \frac{N}{(p_F t)^d} \left( 1 - \frac{r^2}{(p_F t)^2} \right)^{d/a}$$
     Aqui, $p_F$ é o momento de Fermi, que escala com $N$ como $p_F \sim N^{a/((a+2)d)}$.

2. Expoente de Borda Universal:
   Uma descoberta central é a identificação de um expoente de borda universal $d/a$. Os perfis de densidade exibem singularidades de lei de potência nas bordas móveis ($r = p_F t$), em vez das singularidades de raiz quadrada características da lei do semicírculo de Wigner (que corresponde ao caso específico $d=1, a=2$) ou da distribuição de Thomas-Fermi ($d=3, a=2$). Os resultados generalizam essas leis conhecidas para todos os potenciais de lei de potência e dimensões.

3. Comportamento de Borda e Funções de Escalonamento:

   • Armadilha Harmônica ($a=2$): A borda da lei do semicírculo de Wigner é mostrada como sendo alargada em um crossover suave com uma largura relativa escalando como $N^{-2/3}$. O perfil próximo à borda é descrito por uma função de escalonamento envolvendo a derivada da função de Airy, $A(\xi) = \text{Ai}'(\xi)^2 - \xi \text{Ai}(\xi)^2$. O número médio de férmions fora do intervalo de Wigner é uma constante pequena e finita ($\approx 0,061$).
   • Potencial de Poço Quadrado ($a \to \infty$): No limite de um poço infinito, a distribuição de momento torna-se uniforme dentro da bola de Fermi. No entanto, a queda para zero na borda é descrita por uma função de escalonamento diferente $F(\xi)$ envolvendo funções digama e trigama. Diferente do caso harmônico, o potencial abrupto leva a uma cauda algébrica longa na distribuição ($F(\xi) \sim 1/\xi$), fazendo com que o número de férmions fora do intervalo de Fermi cresça logaritmicamente com $N$ ($N_{out} \sim \ln N$).

4. Efeitos de Dimensões Superiores:
   Em dimensões $d \ge 2$, a presença de degenerescência orbital é abordada. Os autores mostram que, para números específicos de férmions correspondentes a camadas completas de estados degenerados, o sistema possui um estado fundamental isotrópico único. Para um poço circular infinito em 2D, o artigo analisa o componente oscilatório da densidade de estados integrada, demonstrando que o momento de Fermi exibe oscilações irregulares mesmo quando restrito a números de camadas completas, uma característica ligada às órbitas periódicas do bilhar clássico.

Significância e Alegações
Os autores alegam que esses resultados são inéditos, fornecendo um arcabouço analítico abrangente para a expansão de férmions livres que generaliza a lei do semicírculo de Wigner e a distribuição de Thomas-Fermi. O trabalho destaca que a expansão balística em tempos tardios mapeia diretamente a distribuição de momento inicial para o perfil espacial, com o comportamento da borda sendo universalmente determinado pela razão entre a dimensão espacial e o expoente do potencial ($d/a$).

O artigo enfatiza a utilidade dos sistemas de férmions livres em fornecer acesso analítico detalhado à dinâmica fora do equilíbrio, que é frequentemente indisponível em sistemas interagentes. Ele estabelece uma conexão clara entre a geometria do confinamento inicial (via o expoente $a$) e a classe de universalidade da borda da nuvem em expansão. As descobertas sobre o crescimento logarítmico de partículas fora da borda de Fermi em potenciais de poço quadrado versus o número constante em armadilhas harmônicas ilustram como a nitidez do potencial inicial dita a natureza das flutuações de borda. O trabalho serve como um acompanhamento de estudos anteriores sobre probabilidades de retorno, estendendo a análise para perfis de densidade completos e dimensões superiores.