Solutions with one dimensional concentration for a two dimensional Gross-Pitaevskii model with general potential

1. Il Problema Studiato

Il lavoro si concentra sull'analisi delle soluzioni d'onda stazionarie per l'equazione di Gross-Pitaevskii bidimensionale in presenza di un potenziale di intrappolamento generale $W(y)$ e sotto il vincolo di massa unitaria. L'equazione, che descrive i condensati di Bose-Einstein con interazioni attrattive, è data da:
$$i\partial_t \Psi = -\Delta \Psi + W(y)\Psi - a|\Psi|^2\Psi, \quad \text{con} \quad \int_{\mathbb{R}^2} |\Psi|^2 = 1.$$
Cercando soluzioni della forma $\Psi(t, y) = \tilde{u}(y)e^{-i\lambda t}$, il problema si riduce a un'equazione ellittica stazionaria con un parametro di frequenza $\lambda$ da determinare. L'articolo studia il comportamento asintotico delle soluzioni eccitate (non fondamentali) quando il parametro di interazione $a$ tende a un valore critico, il che implica che la frequenza $\lambda \to -\infty$ (o equivalentemente, il parametro di perturbazione singolare $\epsilon = \sqrt{-1/\lambda} \to 0$).

L'obiettivo principale è costruire soluzioni che mostrano un fenomeno di concentrazione lungo curve chiuse lisce (dimensione 1) in $\mathbb{R}^2$, in un contesto di potenziale non simmetrico (non radiale).

2. Metodologia

Gli autori adottano un approccio sofisticato basato sulla riduzione di Lyapunov-Schmidt infinita-dimensionale, estendendo tecniche sviluppate in lavori precedenti (in particolare di Del Pino, Kowalczyk e Wei) per gestire la mancanza di simmetria radiale.

Le fasi principali della metodologia includono:

• Coordinate di Fermi e Trasformazione: Il problema viene riformulato utilizzando coordinate di Fermi $(t, \theta)$ attorno alla curva di concentrazione $\Gamma_\epsilon$. Viene introdotta una trasformazione di Liouville per mappare la lunghezza della curva (che dipende da $\epsilon$) su un dominio periodico fisso $[0, 2\pi)$, gestendo così la lunghezza non uniformemente limitata della curva.
• Costruzione di Soluzioni Approssimate: Vengono costruite soluzioni approssimate di ordine superiore attraverso un processo iterativo di correzione. Si parte da un profilo locale basato sulla soluzione fondamentale $U$ dell'equazione $U'' - U + U^3 = 0$, corretto da termini che dipendono dai parametri di deformazione della curva ($f$) e dall'ampiezza della perturbazione ($e$).
• Procedura di "Gluing" (Incollaggio): Viene utilizzato uno schema di incollaggio interno-esterno per unire la soluzione approssimata locale (vicino alla curva) con il comportamento asintotico globale (decadimento esponenziale all'infinito). Questo porta a un problema proiettato su una striscia infinita.
• Riduzione a un Sistema di Equazioni Differenziali: Attraverso la proiezione degli errori residui sugli autovettori del linearizzato (specificamente $U'$ e $Z$, dove $Z$ è l'autofunzione associata all'autovalore positivo), il problema infinito-dimensionale viene ridotto a un sistema accoppiato di equazioni differenziali ordinarie non lineari per i parametri $f(\tau)$ e $e(\tau)$.
• Condizioni di Non-Degenerazione e Gap: Per risolvere il sistema ridotto, vengono imposte condizioni di non degenerazione sulla curva stazionaria e una "condizione di gap" (gap condition) per evitare fenomeni di risonanza legati agli autovalori del problema lineare, che dipendono dalla lunghezza pesata della curva.

3. Contributi Chiave e Risultati Principali

A. Condizioni Necessarie (Teorema 2.2)

Gli autori derivano le condizioni necessarie affinché esista una soluzione che concentra lungo una curva chiusa liscia $\Gamma_\epsilon$:

1. La lunghezza della curva $\ell$ deve essere limitata inferiormente da una costante positiva indipendente da $\epsilon$.
2. Il parametro di massa $a$ deve tendere a $+\infty$ quando $\epsilon \to 0$.
3. Il potenziale $W$ e la curvatura $\kappa$ della curva devono soddisfare una relazione asintotica specifica:
   $$\frac{3\epsilon^2 W_t(0, \theta)}{2[1 + \epsilon^2 W(0, \theta)]} + \kappa(\theta) = o(1).$$
   Questo risultato risponde parzialmente a una congettura di Ambrosetti, Malchiodi e Ni riguardo alle condizioni necessarie per la concentrazione su sottovarietà.

B. Esistenza di Soluzioni (Teorema 2.5 e 2.6)

Il risultato principale è la dimostrazione dell'esistenza di soluzioni positive che concentrano lungo curve $\Gamma_\epsilon$ che soddisfano le condizioni necessarie sopra citate.

• Potenziale Generale: A differenza di lavori precedenti che richiedevano simmetria radiale, questo risultato vale per potenziali $W$ non simmetrici.
• Lunghezza Non Limitata: Una novità significativa è che la lunghezza della curva di concentrazione $\ell$ non è uniformemente limitata; può crescere al diminuire di $\epsilon$. Questo introduce sfide tecniche uniche rispetto ai casi studiati in letteratura (es. [21]), dove le curve avevano lunghezza limitata.
• Stima Asintotica: Viene fornita la forma esplicita della soluzione:
  $$\tilde{u}_\lambda(y) \approx \sqrt{\frac{-\lambda}{a}} \left[1 - \frac{1}{\lambda}W(0, \theta)\right]^{1/2} U\left( \frac{[1 - \frac{1}{\lambda}W(0, \theta)]^{1/2} t}{\sqrt{-1/\lambda}} \right).$$

4. Significato e Impatto

• Estensione della Teoria di Riduzione: Il lavoro estende la teoria della riduzione di Lyapunov-Schmidt infinita-dimensionale a problemi con vincoli di massa e potenziali non simmetrici in due dimensioni, un caso molto più complesso rispetto ai modelli con potenziale radiale.
• Gestione della Risonanza: Gli autori affrontano con successo fenomeni di risonanza più delicati rispetto ai lavori precedenti, causati dalla dipendenza degli autovalori dal parametro $\epsilon$ e dalla lunghezza variabile della curva. L'introduzione di nuove norme pesate e la condizione di gap (B2) sono cruciali per garantire la risolubilità del sistema ridotto.
• Modelli Fisici: I risultati offrono una comprensione più profonda della formazione di strutture vorticali o di concentrazione in condensati di Bose-Einstein con interazioni attrattive in potenziali di intrappolamento generici, dove la simmetria non è garantita.

In sintesi, questo articolo rappresenta un avanzamento teorico significativo nella teoria delle perturbazioni singolari per equazioni ellittiche non lineari, fornendo un quadro rigoroso per l'esistenza di soluzioni concentrate su curve in contesti fisici realistici e non simmetrici.