Rotating two-dimensional Bose--Einstein condensates in power-law traps: collapse, Thomas--Fermi scaling, and the inverse-square endpoint

Resumo Técnico: Condensados de Bose–Einstein Bidimensionais Rotativos em Armadilhas de Lei de Potência

Definição do Problema
Este trabalho investiga as propriedades do estado fundamental e a estabilidade de condensados de Bose–Einstein (BECs) bidimensionais descritos pela equação de Gross–Pitaevskii (GP) sob potenciais de armadilha de lei de potência $V_s(r) = r^s/s$, onde o expoente $s$ varia de $-2$a$4$. O estudo aborda a interdependência entre energia cinética, interações de contato (atrativas ou repulsivas) e rotação externa através desta família de armadilhas. Um foco central é o comportamento de escala do tamanho do condensado próximo ao limiar de colapso para interações atrativas e o limite de Thomas–Fermi (TF) para interações repulsivas, particularmente na presença de rotação e potenciais singulares (especificamente o potencial inverso ao quadrado para $s=-2$).

Metodologia
Os autores utilizam uma combinação de argumentos analíticos de escala e simulações numéricas rigorosas:

1. Análise de Escala: O estudo utiliza argumentos de dilatação para derivar leis de escala para a largura $\ell$ e a energia do condensado. Estabelece que, sob dilatação preservadora de norma, as energias cinética e de contato escalam identicamente ($\lambda^2$), enquanto a energia da armadilha escala como $\lambda^{-s}$.
2. Ferramentas Numéricas:
   • Propagação no Tempo Imaginário 2D: Fluxo de gradiente normalizado usando divisão de Strang para computar estados estacionários tanto para o regime atrativo quanto para o repulsivo, com e sem rotação.
   • Solucionadores Radiais: Soluções radiais de alta precisão (volume finito, iteração de Newton fronteiriçada, extrapolação de Richardson) servem como referências para estados não rotativos.
   • Espectros de Bogoliubov: Cálculo das frequências de excitação para determinar limiares de estabilidade (especificamente a frequência de dipolo) para estados centrados.
   • Limites Variacionais: Minimização da energia sobre uma função de tentativa do tipo Townes com reforço de fase (gauge-boosted) para limitar as energias do estado fundamental por cima.
3. Regularização: Armadilhas singulares ($s \le 0$) são regularizadas com um corte $\epsilon$ (por exemplo, $V_{s,\epsilon} = (r+\epsilon)^s/s$) para lidar com divergências numéricas, com atenção específica ao limite $\epsilon \to 0$ para o caso do inverso ao quadrado.

Principais Contribuições e Resultados

1. Leis de Escala e Dinâmica de Colapso

• Atração: O condensado colapsa no acoplamento de Townes $g^*$ independentemente da forma da armadilha ou da velocidade de rotação (desde que a armadilha seja confinante). A aproximação ao colapso segue a lei $\ell \propto \delta^{1/(s+2)}$, onde $\delta = 1 - |g|/g^*$. Esta lei é válida para sete expoentes diferentes ($s \in \{1, 2, 3, 4, -0.5, -1, -1.5\}$).
• Núcleos Regularizados: Para armadilhas singulares com um corte finito $\epsilon$, se a largura do condensado tornar-se comparável a $\epsilon$, a lei de escala transita para uma "lei de cone" $\ell \propto \delta^{1/3}$, determinada pela forma linear local do núcleo regularizado, em vez da lei de potência de campo distante.
• Repulsão: Para gases repulsivos, a nuvem expande-se para um raio de Thomas–Fermi $R_{TF} \propto g^{1/(s+2)}$. Mostra-se que o erro de energia TF divide-se exatamente em uma parte de gradiente de densidade e uma parte de fluxo (desajuste entre o fluxo de vórtice e a rotação rígida), sendo a parte de fluxo comparável ou maior que a parte de gradiente em armadilhas singulares.

2. Rotação e Quebra de Simetria

• Critério de Dipolo: Para nuvens atrativas em armadilhas com $s > 2$, o estado centrado torna-se instável quando a frequência de rotação $\Omega$ excede a frequência de dipolo $\omega_1$ da nuvem.
  • Para $s=2$ (harmônica), $\omega_1$ é constante (modo de Kohn), de modo que a rotação não induz quebra de simetria abaixo do limite centrífugo.
  • Para $s > 2$, a armadilha é "suave" (soft) perto do centro; $\omega_1$ diminui conforme o condensado se concentra ($\omega_1 \propto \ell^{(s-2)/2}$). Consequentemente, a rotação pode desestabilizar o estado centrado mesmo para uma atração fraca.
• Armadilha Quartica ($s=4$): O limiar para a quebra de simetria é encontrado como $\delta_c \simeq 0.324 \Omega^6$. Este resultado, derivado de uma tentativa de Townes, é confirmado por cálculos de Bogoliubov dentro de 3%.
• Estado Pós-Instabilidade: Uma vez instável, a nuvem atrativa desloca-se do eixo e colapsa sobre o mínimo do anel do potencial no referencial rotativo $V_{eff} = V - \frac{1}{2}\Omega^2 r^2$. A largura deste fragmento localizado no anel é definida pela curvatura da armadilha pura, não pelo potencial efetivo.
• Rotação Repulsiva: Em gases repulsivos, a rotação nucleia vórtices. Em armadilhas quárticas, a densidade TF esvazia-se no centro para formar um anel em uma rotação reduzida crítica $\zeta_c = (3/\pi)^{1/6}$, independente da força de interação $g$.

3. O Ponto Final do Inverso ao Quadrado ($s = -2$)

• Covariância: Em $s=-2$, a energia regularizada é exatamente covariante com o corte $\epsilon$ ($E \propto \epsilon^{-2}$). Não existe escala de comprimento intrínseca; o sistema cai para o centro tanto para interações atrativas quanto repulsivas.
• Estrutura de Halo: Para interações repulsivas, o condensado forma um núcleo denso rodeado por um halo universal com perfil de densidade $|\psi|^2 \propto 1/(gr^2)$. A extensão deste halo, $L_h$, cresce exponencialmente com o acoplamento: $L_h \propto \exp(g / [\pi(1+\kappa)])$, onde $\kappa$ caracteriza a força da cauda.
• Universalidade: O perfil do halo e sua taxa de crescimento exponencial são universais (independentes da forma específica de regularização), enquanto as propriedades do núcleo dependem dos detalhes da regularização.

4. Artefatos Numéricos Próximo ao Colapso

• O artigo identifica um viés estacionário de primeira ordem na propagação de tempo imaginário por divisão de passos (split-step) próximo ao limiar de colapso. Este viés faz com que a largura computada seja ligeiramente maior que a verdadeira função de estado fundamental.
• Crucialmente, este erro é invisível para testes de resíduo padrão (que medem o resíduo da equação GP) porque o erro escala com o modo de respiração (breathing mode). O viés cresce como $1/\delta$ próximo ao limiar. Um esquema de passo simétrico ou a extrapolação para o passo de tempo zero são necessários para eliminar este erro.

Significado e Alegações
Os autores alegam três novos resultados principais:

1. Colapso Induzido por Rotação: O mecanismo preciso e o limiar ($\delta \propto \Omega^6$) pelo qual a rotação desloca um condensado atrativo do eixo em armadilhas suaves ($s>2$), levando ao colapso em anel.
2. Ponto Final do Inverso ao Quadrado: A covariância exata da energia com o corte e a derivação do halo exponencial universal para gases repulsivos no limite $s=-2$.
3. Viés Numérico: A identificação e caracterização de um viés estacionário de primeira ordem em métodos de tempo imaginário próximo ao colapso, que é indetectável pelos critérios de convergência padrão, mas afeta significativamente a largura.

O artigo não alega novidade para resultados estabelecidos, como o limiar de Townes, a lei de escala geral $\ell \propto \delta^{1/(s+2)}$ ou a origem do momento angular no centro de massa em nuvens atrativas; em vez disso, testa estas leis conhecidas quantitativamente através da família de leis de potência e as integra num quadro unificado. O trabalho fornece um relato abrangente de como a geometria da armadilha, o sinal da interação e a rotação determinam conjuntamente a estabilidade e a estrutura de BECs 2D, particularmente em regimes envolvendo potenciais singulares e colapso próximo ao crítico.