Sharp asymptotic behaviour of symmetric and non-symmetric solutions of the Heat Equation in the Hyperbolic Space

1. Il Problema

L'articolo studia il comportamento asintotico a lungo termine ($t \to \infty$) delle soluzioni positive dell'equazione del calore nell'iperspazio iperbolico $\mathbb{H}^d$ (con $d \ge 2$):
$$\begin{cases} u_t = \Delta_{\mathbb{H}^d} u, & x \in \mathbb{H}^d, \, t > 0 \\ u(0, \cdot) = u_0, & x \in \mathbb{H}^d \end{cases}$$
Il focus è su famiglie diverse di dati iniziali $u_0$:

1. Dati radiali: $u_0$ dipende solo dalla distanza geodetica da un polo fisso.
2. Dati orosferici: $u_0$ dipende solo dalla distanza (con segno) da una specifica orosfera.
3. Dati generali: $u_0 \in L^1(\mathbb{H}^d)$ senza particolari simmetrie.

Contesto e Sfide:
A differenza del caso euclideo, dove le soluzioni convergono alla soluzione fondamentale centrata nel baricentro della massa iniziale, nell'iperspazio iperbolico il comportamento è più complesso. La geometria non compatta e la curvatura negativa influenzano la diffusione. Risultati precedenti (es. Vázquez, Papageorgiou) avevano stabilito la convergenza verso profili asintotici, ma spesso senza fornire tassi di convergenza espliciti o, nel caso di dati non radiali, senza identificare completamente il profilo asintotico in $L^1$ per dati generali. Inoltre, la massa della soluzione tende a concentrarsi su sfere geodetiche a distanza $(d-1)t$ dall'origine, rendendo l'analisi classica inadeguata.

2. Metodologia

Il cuore metodologico del lavoro è l'adattamento delle stime di entropia (entropy estimates) a varietà Riemanniane non compatte con curvatura non banale.

• Cambiamento di Variabili e Coordinate:
  Gli autori introducono un cambiamento di variabili temporale e spaziale per "seguire" la massa della soluzione che si sposta verso l'infinito.

  • Tempo: $t(\tau) = e^\tau - 1$.
  • Spazio (radiale): $r(\tau, \rho) = (d-1)e^\tau + e^{\tau/2}\rho$.
    Questo sposta il sistema di riferimento verso la "sfera critica" dove risiede la maggior parte della massa.
    Vengono utilizzate coordinate polari (ellittiche) per soluzioni radiali e coordinate paraboliche (orosferiche) per soluzioni orosferiche.

• Funzione Ausiliaria e Entropia Relativa:
  Si definisce una funzione ausiliaria $v(\tau, \rho, \theta)$ che incorpora il fattore di volume della metrica iperbolica ($\sinh^{d-1}(r)$ o $e^{(d-1)r}$) per lavorare in uno spazio di misura standard.
  Viene introdotta una funzionale di entropia di Kullback-Leibler (entropia relativa):
  $$H(\tau) = \int v \ln\left(\frac{v}{F}\right) d\mu$$
  dove $F$ è un "profilo di equilibrio transitorio" (transient equilibrium) che dipende dal tempo e dalla geometria.

• Disuguaglianze Log-Sobolev:
  Per dimostrare il decadimento dell'entropia, gli autori utilizzano disuguaglianze Log-Sobolev unidimensionali adattate alle funzioni di equilibrio transitorio. Questo permette di ottenere un'equazione differenziale per $H(\tau)$ della forma $\partial_\tau H \le -\lambda H + \text{errori}$, risolvibile tramite il Lemma di Gronwall.

• Gestione dei Dati Non Simmetrici:
  Per il caso generale, viene introdotta una funzione di memoria $\Phi(\theta)$, che cattura la distribuzione angolare della massa iniziale. Il profilo asintotico non è più una semplice soluzione fondamentale, ma una combinazione ponderata di soluzioni fondamentali direzionali.

3. Risultati Chiave

A. Soluzioni Radialmente Simmetriche

Per dati iniziali $u_0(r)$, la soluzione converge alla soluzione fondamentale $P_d(t, x)$ centrata nell'origine.

• Risultato: Viene stabilito un tasso di convergenza preciso in norma $L^1$:
  $$\|u(t, \cdot) - M_R P_d(t, \cdot)\|_{L^1(\mathbb{H}^d)} \lesssim t^{-1/2}$$
  dove $M_R$ è la massa totale.
• Conseguenza: Grazie all'effetto di regolarizzazione $L^1-L^\infty$, si ottiene anche un tasso in $L^\infty$:
  $$\|u(t, \cdot) - M_R P_d(t, \cdot)\|_{L^\infty(\mathbb{H}^d)} \lesssim e^{-\lambda_1 t} t^{-2}$$
  dove $\lambda_1 = (d-1)^2/4$ è il primo autovalore dell'operatore di Laplace-Beltrami.

B. Soluzioni Orosferiche

Per dati iniziali dipendenti dalla distanza da un'orosfera, la soluzione converge a un profilo Gaussiano traslato.

• Risultato: Convergenza in $L^1$ con tasso $t^{-1/2}$ e in $L^\infty$ con tasso $t^{-1}$. Il profilo asintotico è una Gaussiana che si sposta linearmente con la velocità $(d-1)$.

C. Soluzioni Generali (Non Simmetriche)

Questo è il contributo principale del lavoro. Per $u_0 \in L^1(\mathbb{H}^d)$, il profilo asintotico non è universale ma dipende dalla distribuzione angolare della massa iniziale.

• Funzione di Memoria: Viene definita la funzione:
  $$\Phi(\theta_x) = \int_{\mathbb{H}^d} u_0(y) [\cosh(r_y) - \sinh(r_y) \langle \theta_x, \theta_y \rangle]^{-(d-1)} d\mu(y)$$
  Questa funzione rappresenta la "massa direzionale" pesata geometricamente.
• Convergenza Asintotica: La soluzione converge a:
  $$u(t, x) \sim \Phi(\theta_x) P_d(t, x)$$
  in norma $L^1$.
• Tasso di Convergenza: Sotto una condizione di crescita controllata sui momenti del dato iniziale (condizione (11) nel paper, che richiede un momento di ordine $\epsilon$), gli autori dimostrano il tasso di convergenza:
  $$\|u(t, x) - \Phi(\theta_x) P_d(t, x)\|_{L^1(\mathbb{H}^d)} \lesssim t^{-1/2}$$
  Questo migliora i risultati precedenti di Papageorgiou, che non fornivano tassi di convergenza espliciti per il caso $L^1$.

4. Contributi Principali

1. Tassi di Convergenza Espliciti: Per la prima volta, vengono forniti tassi di convergenza precisi ($t^{-1/2}$ in $L^1$) per soluzioni dell'equazione del calore in $\mathbb{H}^d$, sia nel caso simmetrico che in quello generale.
2. Identificazione del Profilo Asintotico in $L^1$: Si identifica chiaramente che per dati non radiali, il profilo non è una singola soluzione fondamentale, ma una sovrapposizione pesata da una funzione di memoria $\Phi(\theta)$ che conserva l'informazione sulla distribuzione iniziale della massa.
3. Adattamento delle Stime di Entropia: Il lavoro dimostra come adattare la teoria delle entropie (tipica dei processi di diffusione su spazi compatti o euclidei) a varietà Riemanniane non compatte con curvatura negativa, gestendo la dipendenza temporale dei profili di equilibrio e dei domini di integrazione.
4. Unificazione del Framework: Viene proposto un approccio unificato basato su strumenti di PDE (equazioni differenziali parziali) che evita l'uso complesso di analisi armonica geometrica, rendendo i risultati più accessibili alla comunità delle PDE.

5. Significato e Implicazioni

Questo studio colma un vuoto significativo nella teoria asintotica delle equazioni di diffusione su varietà non compatte.

• Teorico: Dimostra che la geometria iperbolica induce una "memoria" persistente nella soluzione: a differenza del caso euclideo dove la soluzione "dimentica" la forma iniziale e converge a una Gaussiana universale, in $\mathbb{H}^d$ la distribuzione angolare iniziale influenza permanentemente il profilo asintotico.
• Metodologico: Fornisce un modello per l'analisi di altri processi di diffusione su varietà non compatte, mostrando come le stime di entropia possano essere generalizzate oltre il contesto euclideo.
• Applicativo: I risultati sono rilevanti per la comprensione della diffusione in spazi con curvatura negativa, con potenziali applicazioni in fisica matematica, cosmologia e teoria dei grafi (dove i grafi iperbolici modellano reti complesse).

In sintesi, il paper offre una descrizione completa e quantitativa di come il calore si diffonde nello spazio iperbolico, fornendo per la prima volta una stima precisa della velocità con cui le soluzioni raggiungono il loro stato asintotico, che è intrinsecamente legato alla geometria dello spazio e alla distribuzione iniziale della massa.