Universal self-similar evolution of two-dimensional Bose-Einstein condensates in the acoustic regime

Planteamiento del Problema
El artículo investiga la evolución no estacionaria y autosemejante de condensados de Bose-Einstein (BEC) bidimensionales (2D) impulsados fuera del equilibrio en el régimen acústico. Específicamente, aborda cómo un condensado inicialmente de gran escala evoluciona hacia un cascada turbulenta, transfiriendo energía hacia escalas menores (ultravioleta, UV). Si bien la evolución autosemejante es un fenómeno común en la física (por ejemplo, la difusión, las explosiones nucleares), la naturaleza del escalamiento en sistemas de ondas no lineales se categoriza a menudo en dos clases: de "primer tipo", determinado por el análisis dimensional y las leyes de conservación, y de "segundo tipo", donde los exponentes de escala son determinados por problemas de autovalores no lineales. Los autores pretenden caracterizar la dinámica autosemejante de los BECs 2D, particularmente en presencia de ondas acústicas, y determinar si el sistema sigue una autosemejanza de primer o segundo tipo. Un desafío específico abordado es el comportamiento de los BECs de polaritones, que sufren de vidas finitas y una fuerte disipación, en comparación con los BECs atómicos.

Metodología
El estudio emplea un enfoque dual que combina simulaciones numéricas de la ecuación de Gross–Pitaevskii (GPE) y la Ecuación de Onda Cinética (WKE):

1. Ecuación de Gross–Pitaevskii (GPE): Los autores simulan la GPE 2D, que describe la dinámica de los BECs atómicos y de polaritones. La ecuación incluye términos para la velocidad del sonido ($c_s$), la longitud de curación ($\xi$), el potencial de trampa externo ($V_{trap}$) y un término de amortiguamiento lineal ($\gamma$) para modelar las pérdidas de los polaritones.

   • Las simulaciones se realizaron utilizando el código pseudoespectral FROST.
   • Las condiciones iniciales consistieron en un estado fundamental perturbado por ondas de Bogoliubov débiles y aleatorias en grandes escalas.
   • Se utilizaron dos configuraciones: una trampa cuadrada homogénea de tamaño $75\xi$ y una caja periódica más grande de tamaño $1024\xi$ sin trampa para observar mejor las características universales.

2. Ecuación de Onda Cinética (WKE): Para obtener una visión teórica, los autores utilizaron el marco de la Teoría de la Turbulencia de Ondas Débiles (WWTT). Derivaron y resolvieron numéricamente la WKE para ondas de Bogoliubov de gran escala en 2D.

   • La WKE describe la evolución temporal del espectro de energía ($e_k$) en el límite de ondas débilmente no lineales.
   • Las simulaciones se realizaron utilizando el código WavKinS sobre un rango espectral que abarca cuatro décadas.
   • La WKE incluye un término de amortiguamiento añadido ad-hoc para modelar la disipación, reconociendo que la derivación rigurosa usualmente asume la ausencia de forzamiento o disipación.

3. Análisis: Los autores rastrearon la evolución del espectro de energía y definieron un vector de onda característico $k^*(t)$ (el promedio ponderado por energía) para cuantificar la propagación del frente turbulento. Buscaron soluciones autosemejantes de la forma $e_k(t) = E_0 f(t) \Phi(k/k^*(t))$.

Contribuciones Clave y Resultados

• Identificación del Régimen Autosemejante: Tanto las simulaciones de GPE como las de WKE revelan una evolución autosemejante donde un frente se propaga hacia vectores de onda más altos, dejando atrás un espectro de ley de potencia $e_k \sim k^{-1}$. Este exponente corresponde a la solución de Kolmogorov-Zakharov (KZ) para una cascada directa de energía en ondas acústicas 2D.
• Autosemejanza de Segundo Tipo: Contrario a la expectativa para sistemas de "capacidad infinita" (que típicamente exhiben autosemejanza de primer tipo determinada por el análisis dimensional), el sistema acústico 2D exhibe autosemejanza de segundo tipo. Los exponentes de escala no pueden derivarse únicamente de argumentos dimensionales; en su lugar, emergen de un problema de autovalores no lineales.
• Constante Universal $\beta$: La dinámica del frente en propagación está gobernada por una única constante universal adimensional, $\beta$. Los autores derivaron una ecuación para la posición del frente $k^*(t)$:
  $$\frac{1}{k^*} \frac{dk^*}{dt} \ln\left(\frac{k^*}{k_0}\right) = \frac{\beta}{\tau} e^{-t/\tau_D}$$
  donde $\tau$ es una escala de tiempo característica y $\tau_D$ es el tiempo de disipación.
• Determinación Numérica de $\beta$: Al ajustar la evolución numérica de $k^*(t)$ de las simulaciones de GPE y WKE a la predicción teórica, los autores determinaron el valor universal:
  $$\beta \approx 8.9$$
• Independencia de la Disipación: Un resultado sorprendente es que la forma autosemejante $\Phi(\eta)$ y la constante $\beta$ son independientes del coeficiente de amortiguamiento $\gamma$. A pesar de la presencia de disipación (relevante para los polaritones), el sistema preserva la misma forma autosemejante universal. Esto se atribuye al hecho de que el grado de homogeneidad del término de disipación coincide con el del integral de colisión en la WKE.
• Estructura Espectral: La función autosemejante $\Phi(\eta)$ exhibe un comportamiento de ley de potencia $\eta^{-1}$ para $\eta$ pequeños (consistente con la solución KZ) y un decaimiento exponencial $\exp(-\beta \eta / D)$ para $\eta$ grandes (escalas pequeñas), donde $D$ es un coeficiente de difusión.

Significancia y Reivindicaciones
El artículo afirma establecer la dinámica de los BECs acústicos 2D como un ejemplo específico de un tipo degenerado de solución autosemejante de segundo tipo. La significancia radica en los siguientes puntos:

1. Comportamiento Universal: El sistema está completamente caracterizado por la pareja $(\Phi, \beta)$, donde $\beta \approx 8.9$ es una constante universal independiente de las condiciones iniciales o de la presencia de disipación.
2. Sistema Límite: Los autores sugieren que las ondas acústicas 2D representan un "sistema de capacidad infinita límite". Aunque la integral de energía para la escala KZ es solo logarítmicamente divergente (a diferencia de la divergencia UV en sistemas estándar de capacidad infinita), el sistema hereda propiedades tanto de los regímenes de autosemejanza de primer como de segundo tipo.
3. Validación Teórica: El trabajo demuestra que la Teoría de la Turbulencia de Ondas Débiles (WWTT) predice con éxito la evolución autosemejante de la GPE, incluso en presencia de amortiguamiento, cerrando la brecha entre la dinámica microscópica de la GPE y la descripción macroscópica de la WKE.
4. Relevancia Experimental: Los hallazgos son relevantes tanto para los BECs atómicos (caso sin viscosidad) como para los BECs de polaritones (caso disipativo), proporcionando un marco teórico para interpretar experimentos fuera del equilibrio en estos sistemas, particularmente respecto a la formación de estructuras de pequeña escala y la propagación de frentes turbulentos.

Los autores concluyen que los exponentes de escala en este sistema son seleccionados dinámicamente en lugar de ser impuestos por restricciones dimensionales, resaltando la complejidad de la evolución fuera del equilibrio en fluidos cuánticos 2D.