Eigenvalue stability of Hermitian and normal matrices

1. Problema e Contesto

Il lavoro si inserisce nella teoria delle perturbazioni degli operatori hermitiani e normali, un campo con una lunga storia che risale a Rellich e Kato. Il problema centrale affrontato dagli autori riguarda la regolarità e la stabilità delle mappe che associano a una famiglia di matrici i suoi autovalori.

In particolare, il documento indaga le proprietà di continuità della mappa caratteristica (o "characteristic map"), definita come la composizione di una famiglia di matrici con la mappa degli autovalori ordinati (nel caso hermitiano) o non ordinati (nel caso normale).
Mentre è noto che la mappa degli autovalori è Lipschitziana rispetto alla norma degli operatori (Teorema di perturbazione di Weyl), la questione aperta riguardava il comportamento di questa mappa quando agisce per composizione (superposizione) su spazi di funzioni con regolarità Sobolev ($W^{1,q}$) o Lipschitz ($C^{0,1}$).

Le domande chiave sono:

• La mappa caratteristica è continua quando si considera la topologia di $W^{1,q}$?
• La continuità vale anche nel caso Lipschitz ($q=\infty$)?
• Come si comportano gli autovalori di matrici normali (che sono complessi e non ordinabili globalmente) in spazi Sobolev?

2. Metodologia

Gli autori adottano un approccio analitico e geometrico sofisticato, combinando strumenti di analisi funzionale, teoria delle matrici e geometria algebrica reale.

• Decomposizione Blocco-Diagonale Uniforme: Un ingrediente fondamentale è un "lemma di splitting" per matrici triangolari complesse. Questo permette di dimostrare che, localmente, le matrici hermitiane e normali possono essere diagonalizzate unitariamente in modo uniforme, purché non tutti gli autovalori coincidano. Questo riduce il problema alla stabilità di blocchi di dimensioni inferiori, permettendo un'argomentazione per induzione sulla dimensione $d$ della matrice.
• Funzioni Sobolev a Valore Multi-valore (Almgren): Per il caso delle matrici normali, gli autovalori non possono essere ordinati globalmente in modo continuo. Gli autori utilizzano la teoria delle funzioni Sobolev a valori in spazi metrici di tuple non ordinate (spazi di Almgren, $A_d(\mathbb{C})$), definendo la regolarità Sobolev tramite un embedding bi-Lipschitz in uno spazio euclideo.
• Analisi della Convergenza: La dimostrazione della continuità si basa su:
  • Stime puntuali delle derivate delle curve di matrici.
  • Teoremi di convergenza (Vitali, Egorov) per gestire la convergenza quasi ovunque e l'integrabilità uniforme delle derivate.
  • Argomenti di "sezionamento" per estendere i risultati dal caso monodimensionale ($m=1$) a quello multidimensionale ($m > 1$).

3. Risultati Principali

A. Matrici Hermitiane

Per le matrici hermitiane $d \times d$, la mappa che ordina gli autovalori $\lambda^\uparrow$ induce una mappa caratteristica $E$.

• Continuità in $W^{1,q}$ ($1 \le q < \infty$): La mappa $E: W^{1,q}(U, \text{Herm}(d)) \to W^{1,q}(U, \mathbb{R}^d)$ è continua. Questo è un risultato forte, poiché la composizione con funzioni Lipschitziane non è generalmente continua su spazi Sobolev multidimensionali.
• Caso Lipschitz ($q = \infty$): La mappa $E: C^{0,1}(U, \text{Herm}(d)) \to C^{0,1}(U, \mathbb{R}^d)$ non è continua rispetto alla topologia naturale di $C^{0,1}$. Tuttavia, è continua se il codominio è equipaggiato con la topologia di traccia di $W^{1,q}$ (per ogni $q < \infty$) o con la topologia di Hölder $C^{0,\alpha}$ per $\alpha < 1$.
• Convergenza delle derivate: Se una successione di matrici $A_n$ converge a $A$ in $C^{0,1}$, allora le derivate degli autovalori $\partial_j E(A_n)$ convergono a $\partial_j E(A)$ quasi ovunque e in $L^q$ per ogni $q < \infty$.

B. Matrici Normali

Per le matrici normali, gli autovalori sono trattati come funzioni a valori in $A_d(\mathbb{C})$ (tuple complesse non ordinate).

• Continuità in $W^{1,q}$: La mappa caratteristica $E_u: W^{1,q}(U, \text{Norm}(d)) \to W^{1,q}(U, A_d(\mathbb{C}))$ è continua.
• Caso Lipschitz: Analogamente al caso hermitiano, la continuità non vale nella topologia $C^{0,1}$, ma vale in $W^{1,q}$ e in $C^{0,\alpha}$ per $\alpha < 1$.
• Parametrizzazioni: Il lavoro dimostra che, anche se non esiste una parametrizzazione continua globale degli autovalori per famiglie multidimensionali, esistono parametrizzazioni Sobolev che ereditano la regolarità della famiglia di matrici.

C. Ottimalità e Controesempi

Gli autori dimostrano che i risultati sono ottimali fornendo controesempi espliciti:

• La mappa non è uniformemente continua su $C^{0,1}$.
• Non esiste continuità uniforme nemmeno su sottospazi compatti per certi intervalli di esponenti di Hölder.
• Viene mostrato che la convergenza delle derivate non è uniforme in generale (anche se vale quasi ovunque).

D. Applicazioni

I risultati teorici hanno diverse applicazioni concrete:

1. Valori Singolari: La stabilità degli autovalori si estende immediatamente ai valori singolari delle matrici rettangolari.
2. Numero di Condizionamento: Viene dimostrata la stabilità del numero di condizionamento $\kappa(A)$ in spazi Sobolev, sotto l'ipotesi che la matrice sia non singolare.
3. Area Superficiale: Viene analizzata la stabilità dell'area superficiale del grafico degli autovalori.
4. Operatori Compatti: I risultati sono estesi agli operatori compatti autoaggiunti su spazi di Hilbert, permettendo di trattare spettri infiniti discreti in famiglie Lipschitziane.

4. Significato e Contributi

Questo lavoro rappresenta un avanzamento significativo rispetto alla letteratura precedente (inclusi i lavori recenti di Parusiński e Rainer su polinomi iperbolici) per diversi motivi:

1. Indipendenza dalla Dimensione: A differenza dei risultati sui polinomi, dove la regolarità Sobolev degli autovalori dipende dalla dimensione $d$ (es. $W^{1,q}$ con $q < d/(d-1)$), per le matrici hermitiane e normali la stabilità vale per tutti $1 \le q < \infty$, indipendentemente dalla dimensione $d$.
2. Globalità: Mentre per i polinomi iperbolici la regolarità è spesso locale, per le matrici gli spazi Sobolev sono definiti globalmente.
3. Unificazione: Il trattamento unificato dei casi hermitiano e normale, utilizzando la teoria di Almgren per il caso complesso, fornisce un quadro coerente per la stabilità spettrale.
4. Impatto sulla Fisica e Analisi Numerica: La comprensione della regolarità degli autovalori è cruciale per l'analisi di sistemi dinamici, teoremi adiabatici in fisica quantistica (dove la convergenza delle derivate è essenziale) e per la stabilità numerica di algoritmi che dipendono da matrici dipendenti da parametri.

In sintesi, il paper stabilisce che la mappa degli autovalori è "più regolare" di quanto ci si aspetterebbe dalla teoria generale della composizione di funzioni Lipschitziane, fornendo un quadro robusto per l'analisi di perturbazioni di operatori in spazi di funzioni.