Rates of convergence in long time asymptotics of an alignment model with symmetry breaking

1. Il Problema

Il paper studia il comportamento asintotico a lungo termine di un'equazione di Fokker-Planck non lineare derivata da un modello di Cucker-Smale per il "flocking" (greggiamento) in presenza di rumore browniano e una forza di auto-propulsione.

Il modello considera la densità di probabilità $f(t, v)$ delle velocità degli individui in uno spazio omogeneo (senza dipendenza spaziale), governata dall'equazione:
$$\frac{\partial f}{\partial t} = D \Delta f + \nabla_v \cdot \left[ (\nabla \psi_\alpha(v) + v - u_f) f \right]$$
dove $D$ è l'intensità del rumore, $\psi_\alpha$ è un potenziale esterno, e $u_f$ è la velocità media della popolazione.

Il sistema presenta una transizione di fase dipendente dal rumore $D$:

• Regime simmetrico ($D \ge D^*$): Esiste un'unica soluzione stazionaria isotropa ($G_0$, con velocità media nulla).
• Regime di rottura della simmetria ($D < D^*$): La soluzione isotropa diventa instabile. Esiste un continuo di soluzioni stazionarie polarizzate $G_u$, dove la velocità media $u$ appartiene a una sfera $S$ di raggio $r(D)$.

La questione aperta: Mentre era noto che per $D < D^*$ e condizioni iniziali con energia libera sufficientemente bassa, la soluzione converge verso l'insieme delle soluzioni stazionarie polarizzate, non era stato dimostrato se:

1. La velocità media $u_f(t)$ convergesse a un unico limite $u_\infty$ (selezione di una singola soluzione stazionaria).
2. La velocità di questa convergenza fosse esponenziale.

2. Metodologia

L'autore impiega un approccio analitico basato sull'analisi dell'energia libera (funzionale di Lyapunov) e tecniche di disuguaglianze funzionali.

• Energia Libera e Entropia Relativa: Viene utilizzata l'energia libera $F[f]$ come funzione di Lyapunov. La convergenza è studiata attraverso la discesa dell'energia libera verso il suo minimo $F^*$ e la convergenza dell'entropia relativa $H[f|G_u]$.
• Stime di tipo Polyak-Łojasiewicz: Un punto cruciale è l'analisi dettagliata del potenziale $V(u)$ (legato alla costante di normalizzazione delle distribuzioni di Gibbs). L'autore dimostra stime di tipo Polyak-Łojasiewicz che legano il valore dell'energia libera alla distanza dal set di minimi $S$, permettendo di controllare la velocità di decadimento.
• Disuguaglianze Funzionali: Vengono utilizzate disuguaglianze di Log-Sobolev e Poincaré generalizzate rispetto alla misura di Gibbs $G_u dv$. Queste sono essenziali per legare il tasso di decadimento dell'entropia (Fisher information) alla distanza dalla soluzione stazionaria.
• Linearizzazione e Proiezione: Per ottenere la convergenza esponenziale, l'autore evita la linearizzazione attorno alla soluzione stazionaria limite $G_{u_\infty}$ (che potrebbe non avere un gap spettrale se la convergenza è tangenziale alla sfera $S$). Invece, considera la perturbazione attorno alla proiezione dinamica $u^*(t)$ della velocità media sulla sfera $S$. Questo permette di mantenere una coercività uniforme dell'operatore linearizzato.

3. Contributi Chiave e Risultati

A. Identificazione del Limite Unico (Teorema 1.1)

Il primo risultato principale risolve un problema aperto da tempo. Sotto l'ipotesi che l'energia libera iniziale sia inferiore a quella dello stato isotropo ($F[f_{in}] < F[G_0]$) e che $D < D^*$:

• La velocità media $u_f(t)$ converge a un unico limite $u_\infty \in S$.
• Di conseguenza, la soluzione $f(t)$ converge a una singola distribuzione di Gibbs polarizzata $G_{u_\infty}$.
• La convergenza è garantita in norma $L^2$ pesata e in entropia relativa, con un tasso algebrico (per dati iniziali generici in $L^1$).

B. Tasso di Convergenza Esponenziale (Teorema 1.2)

Il risultato più significativo è la prova della convergenza esponenziale sotto condizioni aggiuntive sui dati iniziali (appartenenza a spazi $L^2$ pesati per tutte le direzioni possibili su $S$):
$$H[f(t)|G_{u_\infty}] + |u_f(t) - u_\infty|^2 = O(e^{-\beta t})$$
Questo dimostra che il sistema non solo si stabilizza su una direzione specifica, ma lo fa a un tasso esponenziale, superando le limitazioni dei risultati precedenti che fornivano solo tassi algebrici o non specificavano la selezione del limite.

C. Analisi della Transizione di Fase Critica

Il paper fornisce anche una prova alternativa e più generale per il caso isotropo ($D \ge D^*$), includendo il caso critico $D = D^*$, dimostrando la convergenza esponenziale (o algebrica a seconda delle condizioni) verso lo stato isotropo $G_0$.

4. Significato e Impatto

Questo lavoro è significativo per la teoria dei sistemi di particelle interagenti e la dinamica dei fluidi biologici per diversi motivi:

1. Selezione della Simmetria: Dimostra rigorosamente come, in presenza di rumore basso, un sistema simmetrico "rompa" la simmetria e selezioni spontaneamente una direzione di movimento specifica, risolvendo l'ambiguità sulla non unicità del limite.
2. Velocità di Convergenza: Stabilisce che la dinamica di allineamento è estremamente rapida (esponenziale) una volta che il sistema è vicino allo stato stazionario, il che ha implicazioni per la stabilità dei modelli di flocking in applicazioni reali.
3. Metodologia Robusta: L'uso combinato di stime di Polyak-Łojasiewicz e della proiezione dinamica sulla varietà dei minimi offre un nuovo strumento analitico per studiare la stabilità di sistemi non lineari con simmetrie continue, superando le difficoltà legate alla mancanza di un gap spettrale negli operatori linearizzati standard.
4. Generalità: I risultati si applicano a una classe di dati iniziali più ampia rispetto alla letteratura precedente, rafforzando la validità del modello di Cucker-Smale con rumore e auto-propulsione come descrizione affidabile del comportamento collettivo.

In sintesi, Surin fornisce una comprensione completa della dinamica a lungo termine di questo modello di allineamento, chiudendo il cerchio sulla convergenza verso stati stazionari unici e quantificando la rapidità con cui tale stato viene raggiunto.