Representation-induced superposition breakdown in linear physics

Este artículo demuestra que el principio de superposición puede fallar en medios multicapa al expresar campos como series infinitas de ondas evanescentes, y propone el uso de modos de flujo de potencia en una base ortonormal para restaurar la convergencia y garantizar la conservación de la energía en sistemas de ondas escalares, electromagnéticas y elásticas.

Lo esencial

🌊 El Problema: Cuando las Ondas se "Ahogan" en el Cálculo

Imagina que estás intentando predecir cómo se comporta la luz (o el sonido) cuando atraviesa varias capas de vidrio o plástico, como si fuera un sándwich de materiales.

En la física clásica, tenemos una regla de oro llamada Principio de Superposición. Es como decir: "Si lanzas una pelota y rebota en una pared, y luego lanzas otra, el resultado final es simplemente la suma de ambas". En sistemas lineales (como la luz en el vidrio), esto suele funcionar perfecto. Los científicos usan esta regla para sumar infinitos "rebotes" internos de la luz para calcular cuánto pasa a través del sándwich.

Pero aquí está el truco:
En ciertos casos (cuando hay muchas capas y la luz intenta atravesarlas en ángulos muy raros), esta suma infinita se rompe. En lugar de dar un número lógico, la matemática explota y da resultados infinitos o absurdos.

El artículo dice que esto no es un error de la computadora ni un fallo de la física real. El problema es que los científicos estaban usando la herramienta equivocada para sumar.

🧱 La Analogía: Contar con Monedas Rotas

Imagina que quieres sumar el valor de un montón de monedas para saber cuánto dinero tienes.

• El método antiguo (Modos de onda estándar): Estás usando monedas que, en teoría, valen 1 euro, pero algunas de ellas son "fantasma". Son monedas que existen en el papel pero que no tienen peso real ni valor en el mundo físico (en física, estas son las ondas evanescentes: existen, pero no transportan energía de un lado a otro).
• Cuando intentas sumar un montón infinito de estas "monedas fantasma" que no se pueden medir correctamente, tu calculadora se vuelve loca. El total se dispara al infinito.

El problema no es que el dinero no exista, sino que estás intentando sumar cosas que no se pueden normalizar (no se pueden pesar) con la balanza que estás usando.

💡 La Solución: Cambiar la Balanza (Modos de Flujo de Potencia)

Los autores, Michael Mazilu y Andriejus Demčenka, proponen una solución brillante: Cambia la forma en que cuentas.

En lugar de sumar "ondas" (que pueden ser fantasmales), proponen sumar energía real (flujo de potencia).

• La nueva herramienta: Imagina que en lugar de contar monedas, usas una balanza que solo pesa lo que realmente se mueve y transporta energía.
• El resultado: Al definir sus "monedas" (o modos de onda) basándose en cuánta energía real transportan, todas las monedas tienen un peso definido y positivo. De repente, la suma infinita deja de explotar y converge a un número lógico y estable.

🚦 La Analogía del Tráfico

Piensa en el tráfico en una autopista con varios peajes (las capas de material).

• El problema: Si intentas calcular cuántos coches pasan sumando cada coche individualmente, pero algunos "coches" son ilusorios (ondas que se desvanecen y no llegan a ningún lado), tu contador se desborda.
• La solución de los autores: En lugar de contar coches, cuentan cuánta gasolina se consume en total. Como la gasolina es algo real y medible, no importa cuántos "coches fantasma" haya en el cálculo; la cantidad total de gasolina siempre será finita y correcta.

🌟 ¿Por qué es importante esto?

1. No es un error numérico: Antes, muchos pensaban que si los cálculos fallaban, era culpa de la computadora o de que había que "ajustar" los números (una técnica llamada regularización). Este paper dice: "No, la física es correcta, pero nuestra forma de escribirla matemáticamente estaba mal elegida".
2. Funciona para todo: Esta idea no solo sirve para la luz (óptica), sino también para el sonido (acústica), las ondas sísmicas (terremotos) e incluso para la mecánica cuántica (donde la "energía" se convierte en "probabilidad").
3. Estabilidad: Al usar estos nuevos "modos de flujo de potencia", los cálculos nunca se rompen, incluso en situaciones extremas donde la luz o el sonido casi no pueden atravesar el material.

En Resumen

El artículo nos enseña que la forma en que representamos la realidad importa. A veces, la física funciona perfectamente, pero si intentamos describirla con las "reglas de suma" incorrectas (usando componentes que no se pueden medir), todo se desmorona.

Al cambiar a una representación basada en energía real y conservada (como el flujo de potencia), recuperamos la estabilidad y la lógica, asegurando que la suma de infinitos rebotes siempre tenga sentido. Es como pasar de intentar contar sombras a contar la luz real que las produce.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "Representation-induced superposition breakdown in linear physics" (Ruptura de la superposición inducida por la representación en física lineal), escrito por Michael Mazilu y Andriejus Demčenka.

1. El Problema: Ruptura de la Superposición en Medios Multicapa

El principio de superposición es fundamental en la teoría de ondas lineales, permitiendo que las soluciones se expresen como combinaciones lineales de funciones base. Sin embargo, los autores identifican una falla crítica en sistemas multicapa (con tres o más interfaces) cuando los campos se expanden en series infinitas de ondas evanescentes e inhomogéneas utilizando bases convencionales (como amplitudes de ondas planas estándar).

• La Paradoja: Aunque las ecuaciones gobernantes son lineales, las series de dispersión múltiple (como las fórmulas de Airy o expansiones de matriz de transferencia) pueden divergir para ciertos rangos de parámetros (específicamente cerca de ángulos críticos o en configuraciones con modos evanescentes).
• Causa Raíz: La divergencia no es un artefacto numérico ni una violación de la linealidad, sino una consecuencia de la representación matemática. Las ondas evanescentes no transportan flujo de energía neta a través de las capas y, por lo tanto, no pueden normalizarse correctamente con respecto a la energía conservada en la base convencional. Esto rompe la unitariedad de la matriz de dispersión y hace que los autovalores del operador de dispersión excedan la unidad, provocando la divergencia de la serie.
• Contexto: Este problema afecta a la óptica (interferencia multicapa), ondas elásticas y sistemas cuánticos, donde la descomposición modal estándar falla en preservar la conservación de energía y la interpretabilidad física.

2. Metodología: Modos de Flujo de Potencia (Power Flux Modes)

Para resolver la divergencia, los autores proponen un cambio fundamental en la base de representación del campo, pasando de las amplitudes de campo estándar a modos propios del flujo de potencia.

• Definición de la Nueva Base: Se construyen modos que son ortonormales con respecto a la componente normal del vector de Poynting (para electromagnetismo) o el flujo de potencia mecánica (para elasticidad).
  • La base se define mediante la normalización de los modos de campo eléctrico/magnético (o desplazamiento/tensión) de tal manera que cada modo contribuya con una cantidad finita y físicamente significativa de energía.
  • Para ondas evanescentes e inhomogéneas, esto implica una combinación específica de modos conjugados que elimina la interferencia destructiva en la medida de energía.
• Reformulación de Operadores:
  • Matriz de Dispersión (Interfaz): En esta nueva base, la matriz de dispersión en las interfaces es unitaria independientemente del valor del vector de onda ($k_z$), garantizando la conservación de la energía.
  • Matriz de Propagación: Los autovalores de la matriz de propagación dentro de las capas están acotados ($\le 1$). Para ondas propagantes y evanescentes son iguales a 1 (en medios sin pérdidas), y menores que 1 para ondas inhomogéneas.
• Análisis Espectral: Se demuestra analíticamente que, en la base de flujo de potencia, el radio espectral del operador de dispersión global (producto de matrices de propagación y dispersión) nunca excede la unidad, lo que garantiza la convergencia de la serie de Neumann (expansión de Airy).

3. Contribuciones Clave

1. Identificación de la Causa de la Divergencia: Se demuestra que la divergencia en series de dispersión múltiple es un problema de representación (falta de normalización adecuada de modos no transportadores de energía) y no un fallo de la física subyacente.
2. Construcción de una Base Unitaria: Desarrollo de una base de modos de flujo de potencia que restaura la unitariedad en las interfaces y acota los autovalores de propagación, eliminando la necesidad de regularizaciones ad hoc (como la suma de Cesàro o continuaciones analíticas) para obtener resultados físicos.
3. Generalización Multidisciplinaria: El marco teórico se aplica exitosamente a:
   • Óptica: Ondas electromagnéticas escalares y vectoriales (incluyendo materiales anisotrópicos).
   • Elastodinámica: Ondas acústicas y elásticas en medios sólidos estratificados.
   • Mecánica Cuántica: Se establece un paralelo con la conservación de la corriente de probabilidad y la unitariedad en matrices S cuánticas.
4. Resolución Numérica y Analítica: Se proporciona una formulación matricial recursiva que permite calcular campos reflejados y transmitidos en estructuras multicapa complejas sin inestabilidad numérica, incluso en regímenes de modos evanescentes fuertes.

4. Resultados

• Convergencia Garantizada: Las simulaciones numéricas muestran que, mientras que la expansión estándar de ondas planas diverge (los autovalores máximos $|\lambda| > 1$) en regiones de ángulos críticos y modos evanescentes, la expansión en la base de flujo de potencia mantiene $|\lambda| \le 1$, asegurando la convergencia de la serie.
• Conservación de Energía: En la base estándar, se observa un desequilibrio de potencia (no localidad) donde la potencia incidente parece depender de las propiedades de capas remotas debido a interferencias. En la base de flujo de potencia, la potencia incidente, reflejada y transmitida es ortogonal y se conserva estrictamente ($\Delta P = 0$).
• Validación con Pulsos Ultracortos: Se simuló la reflexión de pulsos ópticos ultracortos en una estructura de dos capas. La reconstrucción temporal usando modos estándar mostró divergencia y artefactos no físicos, mientras que el uso de modos de flujo de potencia reprodujo fielmente la dinámica del pulso con estabilidad.
• Aplicación a Materiales Anisotrópicos: En cristales birrefringentes, donde la mezcla de modos es compleja, la nueva base elimina la degeneración de mezcla que causaba divergencia en los modos estándar, manteniendo la causalidad y la dirección del flujo de energía.

5. Significado e Impacto

Este trabajo tiene implicaciones profundas para la física de ondas lineales:

• Fundamentos Teóricos: Corrige una limitación histórica en la teoría de dispersión múltiple y la óptica de multicapas, estableciendo que la elección de la base de representación es crítica para la validez física de la superposición en presencia de campos evanescentes.
• Herramienta de Diseño: Proporciona un marco robusto para el diseño de dispositivos fotónicos, sensores y estructuras elásticas donde los efectos de campo cercano y modos evanescentes son dominantes (ej. nanofotónica, acoplamiento cuántico, localización de Anderson).
• Eliminación de Regularizaciones: Ofrece una solución "intrínseca" al problema de divergencia, evitando técnicas de regularización matemática que carecen de interpretación física causal (como la adición de absorción ficticia y su extrapolación a cero).
• Universalidad: Al demostrar que el problema y la solución son análogos en electromagnetismo, elasticidad y mecánica cuántica, unifica el tratamiento de la dispersión en medios estratificados bajo el principio de conservación del flujo de energía.

En resumen, el artículo demuestra que la "ruptura" de la superposición es una ilusión matemática causada por una base de representación inadecuada, y que la restauración de la base basada en el flujo de potencia resuelve el problema, garantizando estabilidad, convergencia y conservación de energía en todos los sistemas de ondas lineales.