Arbitrary harmonic functions as Bose--Einstein condensates

Resumo Técnico: Funções Harmônicas Arbitrárias como Condensados de Bose–Einstein

Enunciado do Problema
O artigo investiga a natureza dos condensados de Bose–Einstein (BEC) que podem surgir no limite termodinâmico de um gás de Bose ideal quando as condições de contorno para o Laplaciano na fronteira do domínio são escolhidas livremente. Embora observações anteriores (por exemplo, em [5]) tenham indicado que gases de Bose ideais de volume infinito admitem muitos estados KMS diferentes, caracterizados por formas sesquilineares positivas arbitrárias na função de dois pontos — correspondendo a funções harmônicas que descrevem condensados — permanecia uma questão em aberto se tais estados poderiam ser obtidos rigorosamente via o procedimento padrão do limite termodinâmico (tomando o limite dos estados de Gibbs de sistemas de volume finito). Além disso, embora estados KMS não invariantes por translação com condensados lineares fossem conhecidos [6], a possibilidade de realizar condensados descritos por funções harmônicas arbitrárias, ou múltiplos condensados simultaneamente, através de condições de contorno específicas, ainda não havia sido estabelecida.

Metodologia
Os autores abordam este problema construindo uma sequência específica de realizações autoadjuntas positivas do Laplaciano, denotadas por $-\Delta_{\Omega_j}^\Phi$, em uma sequência de domínios limitados $\Omega_j \subset \mathbb{R}^d$ que esgotam $\mathbb{R}^d$.

1. Construção do Operador:
   O núcleo da metodologia envolve definir o inverso do Laplaciano em um domínio $\Omega$ via uma perturbação do inverso do Laplaciano de Dirichlet. Dada uma sequência de funções harmônicas $\Phi = \{\phi_k\}_k$ em $\mathbb{R}^d$ satisfazendo $\sum_k |\phi_k(x)|^2 < \infty$, o operador é definido como:
   $$(-\Delta_{\Omega}^\Phi)^{-1} := \left( (-\Delta_{\Omega}^0)^{-1} + \sum_k |\chi_\Omega \phi_k\rangle \langle \phi_k \chi_\Omega| \right)^{-1}$$
   onde $-\Delta_{\Omega}^0$ é o Laplaciano de Dirichlet e $\chi_\Omega$ é a função característica de $\Omega$. Esta construção garante que o novo operador seja uma extensão autoadjunta positiva do Laplaciano em $C_0^\infty(\Omega)$.

2. Condições de Contorno:
   Os autores derivam as condições de contorno explícitas associadas a este operador. Eles mostram que funções no domínio de $-\Delta_{\Omega}^\Phi$ satisfazem uma condição de contorno do tipo Robin generalizada envolvendo o operador de Dirichlet-para-Neumann $H_\Omega$ e um operador específico $r_\Omega^\Phi$ determinado pelas funções harmônicas $\phi_k$. Especificamente, a derivada normal $\partial_n f$ na fronteira está relacionada aos valores de fronteira de $f$ e à projeção de $f$ sobre o subespaço gerado pelas funções harmônicas.

3. Análise de Convergência:
   Para estabelecer o limite termodinâmico, os autores provam um resultado de convergência de resolvente forte (Teorema 2). Eles demonstram que, para qualquer sequência de extensões autoadjuntas não negativas do Laplaciano em domínios $\Omega_j$ convergindo para $\mathbb{R}^d$, os operadores $F(-\Delta_{\Omega_j})$ convergem fortemente para $F(-\Delta)$ em $L^2(\mathbb{R}^d)$ para qualquer função contínua limitada $F$. Este resultado geral é então aplicado à função específica $F(x) = (e^{\beta x} - 1)^{-1} - (\beta x)^{-1}$, que se relaciona com a função de dois pontos do gás de Bose.

Resultados Principais
O principal resultado é apresentado no Teorema 1. Os autores provam que, para qualquer sequência de funções harmônicas $\Phi = \{\phi_k\}_k$ (satisfazendo a condição de somabilidade quadrática), existe uma sequência de condições de contorno tal que a função de dois pontos do gás de Bose ideal no limite termodinâmico converge para:
$$\lim_{j\to\infty} \left\langle f \left| \frac{1}{e^{-\beta \Delta_{\Omega_j}^\Phi} - 1} \right| g \right\rangle = \int_{\mathbb{R}^d} \frac{\hat{f}(p)\overline{\hat{g}(p)}}{e^{\beta|p|^2} - 1} dp + \frac{1}{\beta} \sum_k \left( \int_{\mathbb{R}^d} f \phi_k \right) \left( \int_{\mathbb{R}^d} \overline{g \phi_k} \right)$$
Características fundamentais deste resultado incluem:

• Condensados Arbitrários: O sistema suporta um número arbitrário (finito ou infinito) de condensados de Bose–Einstein.
• Perfis Harmônicos: Cada condensado é descrito por uma função harmônica arbitrária $\phi_k$. A densidade do $k$-ésimo condensado é proporcional a $|\phi_k(x)|^2$.
• Potencial Químico Zero: O resultado é válido para potencial químico zero. Os autores observam que o resultado se estende a potenciais químicos levemente negativos que desaparecem suficientemente rápido conforme o volume aumenta.
• Restrições Dimensionais: Em dimensões $d \in \{1, 2\}$, a convergência requer a suposição adicional de que as transformadas de Fourier das funções de teste se anulem em zero ($\hat{f}(0) = \hat{g}(0) = 0$) para garantir a integrabilidade do termo do gás livre em baixos momentos.

Significância e Alegações
O artigo afirma resolver a questão de se os diversos estados KMS observados em volume infinito podem ser realizados como limites de estados de Gibbs de volume finito. Ao construir explicitamente as condições de contorno, os autores mostram que:

1. É possível projetar um sistema onde as funções de onda dos condensados não estão restritas a constantes ou funções lineares, mas podem ser quaisquer funções harmônicas.
2. Múltiplos condensados podem coexistir no limite termodinâmico, um fenômeno distinto dos "múltiplos condensados" observados em volumes finitos com condições de contorno de Dirichlet (onde o efeito desaparece no limite da função de dois pontos).
3. Os autovalores dos Laplacianos construídos devem diminuir como $(\int_\Omega |\phi_k|^2)^{-1}$ conforme o domínio se expande, garantindo que as projeções espectrais associadas a esses pequenos autovalores produzam contribuições finitas e não nulas para a densidade do condensado no limite.

Os autores mantêm um escopo modesto, focando estritamente na existência matemática de tais estados dentro do arcabouço de gases de Bose ideais e do limite termodinâmico, sem propor realizações experimentais ou estender seus resultados para sistemas interagentes. O trabalho baseia-se na teoria espectral de extensões autoadjuntas e na convergência de resolventes para unir problemas de valor de contorno de volume finito e mecânica estatística de volume infinito.