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

Enunciato del Problema
Il documento investiga l'evoluzione non stazionaria e auto-similare di condensati di Bose-Einstein (BEC) bidimensionali (2D) guidati fuori dall'equilibrio nel regime acustico. Nello specifico, affronta come un condensato inizialmente su larga scala evolva in un cascata turbolenta, trasferendo energia verso scale più piccole (ultravioletto, UV). Sebbene l'evoluzione auto-simile sia un fenomeno comune in fisica (ad esempio, diffusione, esplosioni nucleari), la natura della scalatura nei sistemi d'onda non lineari è spesso categorizzata in due classi: "del primo tipo", determinata dall'analisi dimensionale e dalle leggi di conservazione, e "del secondo tipo", dove gli esponenti di scalatura sono determinati da problemi di autovalori non lineari. Gli autori mirano a caratterizzare la dinamica auto-simile dei BEC 2D, in particolare in presenza di onde acustiche, e a determinare se il sistema segua una auto-similarità del primo o del secondo tipo. Una sfida specifica affrontata riguarda il comportamento dei BEC di polaritoni, che soffrono di tempi di vita finiti e forte dissipazione, rispetto ai BEC atomici.

Metodologia
Lo studio impiega un approio duale che combina simulazioni numeriche dell'equazione di Gross–Pitaevskii (GPE) e dell'Equazione Cinetica delle Onde (WKE):

1. Equazione di Gross–Pitaevskii (GPE): Gli autori simulano la GPE 2D, che descrive la dinamica dei BEC atomici e di polaritoni. L'equazione include termini per la velocità del suono ($c_s$), la lunghezza di coerenza ($\xi$), il potenziale di intrappolamento esterno ($V_{trap}$) e un termine di smorzamento lineare ($\gamma$) per modellare le perdite dei polaritoni.

   • Le simulazioni sono state eseguite utilizzando il codice pseudospettrale FROST.
   • Le condizioni iniziali prevedevano uno stato fondamentale perturbato da deboli onde di Bogoliubov casuali a grandi scale.
   • Sono stati utilizzati due setup: una trappola quadrata omogenea di dimensione $75\xi$ e un box periodico più grande di dimensione $1024\xi$ senza trappola, per osservare meglio le caratteristiche universali.

2. Equazione Cinetica delle Onde (WKE): Per ottenere un approfondimento teorico, gli autori hanno utilizzato il framework della Teoria della Turbolenza d'Onda Debole (WWTT). Hanno derivato e risolto numericamente la WKE per onde di Bogoliubov a grande scala in 2D.

   • La WKE descrive l'evoluzione temporale dello spettro di energia ($e_k$) nel limite di debole non linearità.
   • Le simulazioni sono state eseguite utilizzando il codice WavKinS su un intervallo spettrale che copre quattro decadi.
   • La WKE include un termine di smorzamento aggiunto ad-hoc per modellare la dissipazione, riconoscendo che la derivazione rigorosa assume solitamente l'assenza di forzanti o dissipazione.

3. Analisi: Gli autori hanno tracciato l'evoluzione dello spettro di energia e definito un vettore d'onda caratteristico $k^*(t)$ (la media pesata sull'energia) per quantificare la propagazione del fronte turbolento. Cercavano soluzioni auto-similari della forma $e_k(t) = E_0 f(t) \Phi(k/k^*(t))$.

Contributi Chiave e Risultati

• Identificazione del Regime Auto-simile: Sia le simulazioni GPE che quelle WKE rivelano un'evoluzione auto-simile in cui un fronte si propaga verso vettori d'onda più elevati, lasciandosi alle spalle uno spettro a legge di potenza $e_k \sim k^{-1}$. Questo esponente corrisponde alla soluzione di Kolmogorov-Zakharov (KZ) per un cascata diretta di energia in onde acustiche 2D.
• Auto-similarità del Secondo Tipo: Contrariamente all'aspettativa per i sistemi a "capacità infinita" (che tipicamente esibiscono auto-similarità del primo tipo determinata dall'analisi dimensionale), il sistema acustico 2D esibisce una auto-similarità del secondo tipo. Gli esponenti di scalatura non possono essere derivati solo da argomenti dimensionali; essi emergono invece da un problema di autovalore non lineare.
• Costante Universale $\beta$: La dinamica del fronte in propagazione è governata da una singola costante adimensionale universale, $\beta$. Gli autori hanno derivato un'equazione per la posizione del fronte $k^*(t)$:
  $$\frac{1}{k^*} \frac{dk^*}{dt} \ln\left(\frac{k^*}{k_0}\right) = \frac{\beta}{\tau} e^{-t/\tau_D}$$
  dove $\tau$ è una scala temporale caratteristica e $\tau_D$ è il tempo di dissipazione.
• Determinazione Numerica di $\beta$: Effettuando il fitting dell'evoluzione numerica di $k^*(t)$ dalle simulazioni GPE e WKE alla previsione teorica, gli autori hanno determinato il valore universale:
  $$\beta \approx 8.9$$
• Indipendenza dalla Dissipazione: Un risultato sorprendente è che la forma auto-simile $\Phi(\eta)$ e la costante $\beta$ sono indipendenti dal coefficiente di smorzamento $\gamma$. Nonostante la presenza di dissipazione (rilevante per i polaritoni), il sistema preserva la stessa forma auto-simile universale. Ciò è attribuito al fatto che il grado di omogeneità del termine dissipativo coincide con quello dell'integrale collisionale nella WKE.
• Struttura Spettrale: La funzione auto-simile $\Phi(\eta)$ esibisce un comportamento a legge di potenza $\eta^{-1}$ per piccoli $\eta$ (coerente con la soluzione KZ) e un decadimento esponenziale $\exp(-\beta \eta / D)$ per grandi $\eta$ (piccole scale), dove $D$ è un coefficiente di diffusione.

Significato e Rivendicazioni
L'articolo sostiene di aver stabilito la dinamica dei BEC acustici 2D come un esempio specifico di una soluzione auto-simile di tipo degenerato del secondo tipo. Il significato risiede nei seguenti punti:

1. Comportamento Universale: Il sistema è completamente caratterizzato dalla coppia $(\Phi, \beta)$, dove $\beta \approx 8.9$ è una costante universale indipendente dalle condizioni iniziali o dalla presenza di dissipazione.
2. Sistema di Confine: Gli autori suggeriscono che le onde acustiche 2D rappresentino un "sistema a capacità infinita di confine". Sebbene l'integrale di energia per la scalatura KZ sia solo divergentemente logaritmico (a differenza della divergenza UV nei sistemi standard a capacità infinita), il sistema eredita proprietà sia dai regimi di auto-similarità del primo che del secondo tipo.
3. Validazione Teorica: Il lavoro dimostra che la Teoria della Turbolenza d'Onda Debole (WWTT) predice con successo l'evoluzione auto-simile della GPE, anche in presenza di smorzamento, colmando il divario tra la dinamica microscopica della GPE e la descrizione macroscopica della WKE.
4. Rilevanza Sperimentale: Le scoperte sono rilevanti sia per i BEC atomici (caso inviscido) che per i BEC di polaritoni (caso dissipativo), fornendo un quadro teorico per interpretare esperimenti fuori dall'equilibrio in questi sistemi, in particolare riguardo alla formazione di strutture a piccola scala e alla propagazione di fronti turbolenti.

Gli autori concludono che gli esponenti di scalatura in questo sistema sono dinamicamente selezionati piuttosto che imposti da vincoli dimensionali, evidenziando la complessità dell'evoluzione fuori dall'equilibrio nei fluidi quantistici 2D.