Degenerate parabolic equations in divergence form: fundamental solution and Gaussian bounds

1. Il Problema e il Contesto

Il lavoro si concentra sullo studio di equazioni paraboliche degeneri del secondo ordine in forma di divergenza, definite sull'intero spazio-tempo $\mathbb{R} \times \mathbb{R}^n$. L'operatore considerato è della forma:
$$H u := \partial_t u - \omega^{-1} \text{div}_x (A(t, x) \nabla_x u)$$
dove:

• $A(t, x)$ è una matrice con coefficienti complessi, misurabili e dipendenti dal tempo.
• $\omega(x)$ è un peso spaziale, indipendente dal tempo, appartenente alla classe di Muckenhoupt $A_2(\mathbb{R}^n)$.
• La degenerazione dell'ellitticità è governata dal peso $\omega$: l'operatore soddisfa una condizione di ellitticità uniforme rispetto alla misura pesata $d\omega = \omega(x)dx$.

L'obiettivo principale è stabilire l'esistenza, l'unicità e le proprietà di regolarità della soluzione fondamentale (o nucleo di calore generalizzato) $\Gamma(t, x; s, y)$ associata a questo operatore, in particolare cercando di ottenere stime gaussiane superiori e inferiori.

2. Metodologia

L'autore adotta un approccio analitico moderno che combina la teoria degli operatori in spazi di Hilbert pesati con tecniche classiche di regolarità per equazioni paraboliche.

• Struttura Funzionale: Il lavoro si basa su spazi di Sobolev pesati $H^1_\omega(\mathbb{R}^n)$ e sulla teoria degli operatori accoppiati. Viene utilizzata una tecnica di embedding di tipo Lions per gestire la derivata temporale e le identità energetiche.
• Costruzione della Soluzione Fondamentale: L'esistenza e l'unicità della famiglia di operatori $\Gamma(t, s)$ sono ottenute estendendo i risultati di [AB25] (Auscher-Baadi), utilizzando forme sesquilinee coercitive e misurabili. Non si assume regolarità sui coefficienti $A$, che possono essere solo misurabili.
• Strumenti Chiave:
  • Disuguaglianza di Caccioppoli: Derivata dalla forma divergente, permette di controllare l'energia locale delle soluzioni deboli.
  • Stime di decadimento fuori diagonale: Dimostrazione che la soluzione fondamentale decade esponenzialmente lontano dalla diagonale spaziale.
  • Metodo di Moser: Utilizzato per collegare le stime $L^2$-$L^\infty$ delle soluzioni deboli alle proprietà del nucleo.
  • Disuguaglianza di Harnack: Applicata nel caso di coefficienti reali per ottenere stime inferiori.

3. Contributi Chiave e Risultati Principali

Il risultato centrale del paper è il Teorema 1.2, che stabilisce un'equivalenza fondamentale e fornisce stime gaussiane complete.

A. Equivalenza tra Stime Gaussiane e Stime di Moser

Il paper dimostra che l'esistenza di una soluzione fondamentale con stima gaussiana superiore è equivalente al soddisfacimento delle stime di Moser $L^2$-$L^\infty$ per le soluzioni deboli locali sia dell'operatore $H$ che del suo aggiunto $H^*$.
La stima gaussiana superiore assume la forma:
$$|\Gamma(t, x; s, y)| \leq \frac{K_0}{\sqrt{\omega(B(x, \sqrt{t-s}))} \sqrt{\omega(B(y, \sqrt{t-s}))}} e^{-k_0 \frac{|x-y|^2}{t-s}}$$
Questa equivalenza generalizza i risultati classici di Hofmann e Kim (per coefficienti non pesati) al caso degenerato pesato, senza richiedere regolarità sui coefficienti.

B. Il Caso di Coefficienti Reali

Nel caso specifico in cui la matrice $A$ abbia coefficienti reali:

1. Le stime di Moser $L^2$-$L^\infty$ sono automaticamente soddisfatte (grazie alla teoria di regolarità nota per soluzioni deboli di equazioni reali).
2. Di conseguenza, le stime gaussiane superiori sono sempre valide.
3. Utilizzando la disuguaglianza di Harnack (nota per coefficienti reali), l'autore deriva anche le stime gaussiane inferiori:
   $$\tilde{K}_0 \frac{e^{-\tilde{k}_0 \frac{|x-y|^2}{t-s}}}{\sqrt{\omega(B(x, \sqrt{t-s}))} \sqrt{\omega(B(y, \sqrt{t-s}))}} \leq \Gamma(t, x; s, y)$$
   Questo fornisce una stima completa (superiore e inferiore) del nucleo di calore generalizzato.

C. Proprietà Aggiuntive

Il paper stabilisce diverse proprietà fondamentali della soluzione fondamentale:

• Proprietà di Conservazione: La soluzione fondamentale preserva le costanti (in senso locale pesato), ovvero $\Gamma(t, s)1 = 1$.
• Identità di Chapman-Kolmogorov: Valide per la famiglia di operatori.
• Decadimento $L^2$: Stime di decadimento esponenziale per la soluzione fondamentale su insiemi disgiunti.

4. Significato e Impatto

Questo lavoro è significativo per diversi motivi:

1. Generalizzazione: Estende la teoria classica delle equazioni paraboliche (unweighted) e i risultati di Aronson/Nash al contesto degenerato pesato con coefficienti complessi e dipendenti dal tempo.
2. Rimozione di Ipotesi di Regolarità: A differenza di approcci precedenti che richiedevano coefficienti lisci o perturbazioni piccole di coefficienti reali, questo metodo funziona per coefficienti puramente misurabili, basandosi esclusivamente sulle proprietà strutturali dell'operatore e sulle stime di Moser.
3. Unificazione: Fornisce un quadro unificato che collega la teoria degli operatori, le stime di regolarità delle soluzioni deboli e le proprietà asintotiche del nucleo di calore.
4. Applicabilità: I risultati sono rilevanti per l'analisi di potenze frazionarie di equazioni paraboliche, diffusione anomala e nuclei di calore per equazioni di Schrödinger con potenziali singolari, come menzionato nell'introduzione.

In sintesi, Baadi dimostra che la struttura geometrica e analitica imposta dal peso $A_2$ e dalla forma divergente è sufficiente a garantire l'esistenza di un nucleo di calore ben comportato (con stime gaussiane), a patto che le soluzioni deboli godano di proprietà di regolarità locale (stime di Moser), che sono automatiche nel caso reale.