Eigenvalue stability of Hermitian and normal matrices

Resumen Técnico: Estabilidad de Valores Propios en Matrices Hermíticas y Normales

1. Planteamiento del Problema

El problema central de este trabajo es analizar la continuidad de la aplicación que asigna a una familia de matrices sus valores propios. Específicamente, se estudia la aplicación característica (o mapa de características), definida por la superposición con el mapa de valores propios ordenados (en el caso hermítico) o no ordenados (en el caso normal).

• Contexto: La teoría de perturbaciones de operadores hermíticos y normales tiene una larga historia (Rellich, Kato), donde se sabe que los valores propios ordenados forman un mapa Lipschitz en el espacio de matrices.
• La Pregunta: ¿Es esta aplicación característica continua cuando se considera en espacios de funciones con regularidad Sobolev ($W^{1,q}$) o Lipschitz ($C^{0,1}$)?
• Desafío: Aunque la aplicación es Lipschitz en el nivel de matrices individuales, la superposición con funciones de Sobolev no siempre preserva la continuidad en la topología de Sobolev, especialmente para $q=\infty$ (caso Lipschitz) o en dimensiones superiores ($m \ge 2$). Además, para matrices normales, los valores propios son números complejos que no admiten un orden global coherente, requiriendo un tratamiento como funciones multivaluadas.

2. Metodología y Herramientas Principales

Los autores emplean una combinación de análisis funcional, teoría de matrices y geometría algebraica real:

1. Teoremas de Perturbación Clásicos: Se basan fundamentalmente en el teorema de perturbación de Weyl y el teorema de Löwner para matrices hermíticas, y la extensión de Hoffman-Wielandt para matrices normales. Estos establecen que los valores propios son Lipschitz con constante 1 respecto a la norma de Frobenius.
2. Descomposición Unitaria Uniforme (Lema de Desglose): Una pieza clave es un lema de desglose (splitting lemma) para matrices triangulares complejas. Esto permite una diagonalización unitaria de bloques localmente uniforme para matrices hermíticas y normales siempre que no todos los valores propios coincidan.
   • Esto reduce el problema global a problemas locales en subespacios de dimensión menor, facilitando un argumento inductivo sobre el tamaño de la matriz $d$.
3. Funciones Sobolev Multivaluadas: Para el caso de matrices normales, los autores utilizan la teoría de funciones Sobolev con valores en el espacio métrico de $d$-tuplas no ordenadas de números complejos ($A^d(\mathbb{C})$), desarrollada por Almgren y posteriormente por De Lellis y Spadaro. Esto permite tratar los valores propios como funciones bien definidas en espacios de Sobolev sin necesidad de elegir una rama continua global (que a menudo no existe).
4. Análisis de Convergencia: Uso intensivo del teorema de convergencia de Vitali, desigualdades de Morrey y técnicas de "seccionamiento" (reducir el caso multivariable al caso univariable a lo largo de líneas paralelas a los ejes).

3. Contribuciones y Resultados Principales

El artículo establece resultados precisos sobre la continuidad de la aplicación característica $E$ (hermítica) y $E_u$ (normal):

A. Caso Hermítico (Valores Propios Reales Ordenados)

• Continuidad en $W^{1,q}$ ($1 \le q < \infty$): Se demuestra que la aplicación $E: W^{1,q}(U, \text{Herm}(d)) \to W^{1,q}(U, \mathbb{R}^d)$ es continua. Esto significa que si una sucesión de familias de matrices $A_n$ converge a $A$ en la norma $W^{1,q}$, entonces sus valores propios ordenados convergen en la misma norma.
• No Uniformidad: La aplicación no es uniformemente continua en $W^{1,q}$.
• Caso Lipschitz ($q = \infty$):
  • La aplicación $E: C^{0,1}(U, \text{Herm}(d)) \to C^{0,1}(U, \mathbb{R}^d)$ no es continua en la topología Lipschitz natural.
  • Sin embargo, es continua si el espacio objetivo se equipa con la topología de traza de la inmersión en $W^{1,q}$ (para cualquier $q < \infty$).
  • Es continua hacia espacios de Hölder $C^{0,\alpha}$ para $0 < \alpha < 1$, pero no para $\alpha = 1$.
• Convergencia Puntual: Bajo convergencia Lipschitz, las derivadas de los valores propios convergen casi por todas partes (c.p.t.) y en $L^q$ para cualquier $q < \infty$.

B. Caso Normal (Valores Propios Complejos No Ordenados)

• Generalización: Los resultados se extienden a matrices normales, donde los valores propios se interpretan como funciones Sobolev multivaluadas en el sentido de Almgren.
• Continuidad: La aplicación característica $E_u: W^{1,q}(U, \text{Norm}(d)) \to W^{1,q}(U, A^d(\mathbb{C}))$ es continua para todo $1 \le q < \infty$.
• Independencia de la Dimensión: A diferencia de los polinomios hiperbólicos (donde la regularidad Sobolev óptima depende del grado $d$), en el caso de matrices, la estabilidad en $W^{1,q}$ se mantiene para todo $q < \infty$ independientemente del tamaño $d$ de la matriz.

C. Aplicaciones Derivadas
Los autores demuestran que estos resultados tienen implicaciones directas en:

1. Valores Singulares: La estabilidad de los valores singulares de matrices rectangulares hereda las propiedades de continuidad de los valores propios de matrices hermíticas asociadas.
2. Números de Condición: Se prueba la continuidad del número de condición $\kappa(A) = \sigma_1(A)/\sigma_d(A)$ en familias de matrices en $W^{1,q}$, siempre que la matriz sea no singular.
3. Área de Superficie: Se establece la convergencia del área de la superficie de los gráficos de los valores propios bajo perturbaciones Lipschitz.
4. Operadores Compactos Autoadjuntos: Se extienden los resultados a familias de operadores compactos autoadjuntos en espacios de Hilbert, demostrando la estabilidad de sus valores propios ordenados.

4. Optimalidad y Contraejemplos

El artículo es riguroso al demostrar que sus resultados son óptimos mediante ejemplos explícitos (Sección 7):

• Se muestra que la continuidad falla en la topología Lipschitz ($C^{0,1}$) debido a la formación de "picos" en las derivadas de los valores propios cerca de puntos donde los valores propios coinciden (intersección de niveles).
• Se demuestra que la aplicación no es uniformemente continua ni en topologías Lipschitz ni Hölder, incluso en subconjuntos compactos, para ciertos rangos de parámetros.

5. Significado e Impacto

Este trabajo es fundamental por varias razones:

• Unificación: Proporciona un marco unificado para la estabilidad de valores propios en espacios de Sobolev, abarcando tanto el caso hermítico como el normal.
• Superioridad sobre Polinomios: Destaca que el problema de matrices es "más estable" que el problema de raíces de polinomios generales. Mientras que para polinomios la regularidad Sobolev de las raíces depende del grado ($W^{1,q}$ con $q < d/(d-1)$), para matrices la regularidad es $W^{1,q}$ para todo $q < \infty$.
• Herramientas para Análisis No Lineal: La introducción y uso de funciones Sobolev multivaluadas para valores propios de matrices normales abre nuevas vías para el análisis de sistemas dinámicos y ecuaciones diferenciales donde los espectros cruzan o se degeneran.
• Aplicaciones Físicas y Numéricas: Los resultados sobre números de condición y estabilidad de operadores compactos son relevantes para el análisis numérico de problemas mal condicionados y para la mecánica cuántica (teoremas adiabáticos).

En conclusión, Parusiński y Rainer establecen que, aunque la elección de valores propios individuales puede ser discontinua, la estructura global de los valores propios (ordenados o como conjunto) es robusta y continua en espacios de Sobolev, proporcionando una base teórica sólida para el análisis de perturbaciones en sistemas matriciales y operadores.