The $*$-algebra of unbounded GLT: construction and theoretical foundations

1. Il Problema e il Contesto

La teoria delle sequenze GLT (Generalized Locally Toeplitz) è uno strumento consolidato per l'analisi spettrale asintotica di matrici derivate dalla discretizzazione di Equazioni alle Derivate Parziali (PDE) su domini limitati (tipicamente ipercubi o sottoinsiemi regolari di $[0,1]^d$). La teoria associa a ogni sequenza di matrici un "simbolo" canonico (una funzione misurabile) che ne descrive la distribuzione dei valori singolari ed eigenvalori, preservando le operazioni algebriche.

Tuttavia, molte applicazioni fisiche e ingegneristiche (come le onde marine su domini aperti o problemi su domini non limitati) richiedono la discretizzazione su domini illimitati. Sebbene alcuni lavori recenti abbiano esteso la teoria GLT a domini limitati ma irregolari (GLT ridotti), mancava un quadro teorico rigoroso per domini illimitati di misura finita ($L^d(\Omega) < \infty$). La sfida principale risiede nel fatto che su domini illimitati, la struttura della griglia di discretizzazione e la dimensione delle matrici non seguono lo schema standard degli ipercubi, rendendo difficile definire un simbolo unico e un'algebra coerente.

2. Metodologia

Gli autori costruiscono la teoria delle GLT illimitate (unbounded GLT o uGLT) attraverso un approccio costruttivo basato su tre pilastri fondamentali:

1. Estensione agli ipercubi e operatori di restrizione/estensione:

   • Si definiscono algebre GLT su ipercubi generici $Q_{y,l}$.
   • Si introducono operatori lineari di restrizione ($R$) ed estensione ($E$) tra algebre su ipercubi diversi o tra ipercubi e domini limitati. Questi operatori preservano la struttura di *-algebra e la convergenza a.c.s. (approximating class of sequences).

2. GLT Ridotte su domini limitati:

   • Per un dominio limitato $\Omega$, si definisce l'algebra delle GLT ridotte come l'immagine dell'operatore di restrizione applicato a un GLT su un ipercubo contenente $\Omega$. Si dimostra che questa definizione è indipendente dalla scelta dell'ipercubo contenente.

3. Approccio tramite "Regular Exhaustion" e g.a.c.s.:

   • Per un dominio illimitato $\Omega$ con $L^d(\Omega) < \infty$ e $L^d(\partial\Omega)=0$, si utilizza una esaurizione regolare $\{\Omega_t\}_t$ (una famiglia crescente di domini limitati che coprono $\Omega$).
   • Si introduce il concetto di g.a.c.s. (generalized approximating class of sequences), una generalizzazione della classica a.c.s. che permette di approssimare una sequenza di matrici su $\Omega$ con sequenze di dimensioni diverse (definite sui $\Omega_t$) che convergono asintoticamente.
   • Una sequenza su $\Omega$ è definita unbounded GLT se può essere approssimata (in senso g.a.c.s.) da una sequenza di GLT ridotte definite sui domini $\Omega_t$ dell'esaurizione, i cui simboli convergono in misura a una funzione limite su $\Omega \times [-\pi, \pi]^d$.

3. Contributi Chiave e Risultati Principali

A. Costruzione dell'Algebra uGLT

Gli autori definiscono l'insieme $G_\Omega$ di tutte le sequenze uGLT su un dominio $\Omega$. La costruzione dimostra che $G_\Omega$ è una $\ast$-algebra su $\mathbb{C}$.

B. Proprietà Algebriche e Spettrali

Il Teorema 1.1 (e il Teorema 6.11) stabilisce che l'operatore di simbolo $\Phi_\Omega: G_\Omega \to L^0(\Omega \times [-\pi, \pi]^d)$ è un omomorfismo di $\ast$-algebre. Nello specifico, per sequenze $\{A_n\}, \{A'_n\} \in G_\Omega$ con simboli $f, g$:

• Coniugazione: $\{A_n^*\} \in G_\Omega$ con simbolo $\bar{f}$.
• Somma e Prodotto: $\{\alpha A_n + \beta A'_n\}$ e $\{A_n A'_n\}$ appartengono a $G_\Omega$ con simboli $\alpha f + \beta g$ e $fg$.
• Inverso di Moore-Penrose: Se $f \neq 0$ quasi ovunque, allora $\{A_n^\dagger\} \in G_\Omega$ con simbolo $f^{-1}$.

C. Isometria con le Funzioni Misurabili

Il risultato fondamentale è il Teorema 1.2 (e Teorema 6.16), che afferma che la mappa che associa a ogni sequenza uGLT il suo simbolo canonico è una isometria suriettiva tra lo spazio quoziente delle sequenze (rispetto alla distanza a.c.s.) e lo spazio delle funzioni misurabili $L^0(\Omega \times [-\pi, \pi]^d)$ (rispetto alla metrica della convergenza in misura).
Ciò significa che esiste una corrispondenza biunivoca perfetta tra le proprietà asintotiche delle matrici e le proprietà analitiche delle funzioni simbolo, anche su domini illimitati.

D. Unicità del Simbolo

Si dimostra che il simbolo canonico è unico a meno di uguaglianza quasi ovunque, indipendentemente dalla scelta dell'esaurizione regolare o della sequenza approssimante (Corollario 6.10).

E. Applicazioni Numeriche

La teoria è validata su un problema modello: la discretizzazione di un'equazione di Poisson con coefficiente variabile su un dominio illimitato bidimensionale ($\Omega = \{ (x,y) : x>0, y>0, y < g(x) \}$ con $g(x)$ decrescente).

• Vengono costruite matrici di Toeplitz illimitate e di campionamento diagonale illimitate.
• I test numerici confermano che la distribuzione spettrale degli autovalori delle matrici discretizzate converge alla distribuzione prevista dal simbolo teorico $f(x, \theta) = a(x)h(\theta)$, anche al variare della dimensione della griglia e del dominio di troncamento $t$.

4. Significato e Implicazioni

1. Generalizzazione della Teoria GLT: Questo lavoro estende la potenza della teoria GLT a una classe di problemi molto più ampia, permettendo l'analisi spettrale di PDE su domini non limitati ma di misura finita, un caso frequente in fluidodinamica e fisica matematica.
2. Strumento per Domini in Movimento: Gli autori sottolineano che questa struttura algebrica è il prerequisito fondamentale per studiare PDE su domini in movimento (es. onde marine, interfacce fluide). Poiché un dominio in movimento può essere visto come una successione di domini fissi, la teoria uGLT fornisce il linguaggio matematico per tracciare l'evoluzione spettrale della discretizzazione nel tempo.
3. Robustezza Strutturale: La dimostrazione che l'algebra è chiusa sotto operazioni fondamentali (incluso l'inverso generalizzato) garantisce che le tecniche di precondizionamento e analisi di stabilità sviluppate per i domini limitati possano essere estese a questi nuovi contesti.
4. Nuove Direzioni di Ricerca: Il lavoro apre la strada allo studio di classi specifiche di funzioni simbolo (es. $L^p$, Sobolev) all'interno del contesto uGLT e suggerisce possibili estensioni per gestire domini di misura infinita tramite griglie non uniformi.

In sintesi, il paper fornisce le fondamenta teoriche rigorose per trattare l'analisi spettrale asintotica di operatori discretizzati su domini illimitati, colmando un vuoto teorico significativo e offrendo nuovi strumenti per la simulazione numerica di problemi fisici complessi.