🔬 optics

Tutorial of Fourier and Hankel transforms for ultrafast optics

Esta tutorial aclara conceptos erróneos en óptica ultrarrápida sobre señales analíticas y convenciones de transformadas de Fourier, ofrece derivaciones completas y discute múltiples aspectos de la transformada, además de presentar una nueva y superior implementación numérica de alta precisión para la transformada de Hankel.

Autores originales: Yi-Hao Chen

Publicado 2026-02-18

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Yi-Hao Chen

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

¡Hola! Imagina que este artículo es como un manual de instrucciones para cocineros de luz ultrarrápida.

El autor, Yi-Hao Chen, está un poco frustrado porque muchos científicos que estudian pulsos de luz que duran una fracción de un segundo (ultrafastos) están usando las herramientas matemáticas incorrectas, como si intentaran medir la temperatura de un horno con una regla. Esto lleva a errores que a veces ni siquiera se notan porque los resultados "parecen" suaves y bonitos.

Aquí te explico los 5 puntos clave del artículo usando analogías cotidianas:

1. La "Señal Analítica": Separar lo real de lo imaginario

El problema: En la vida real, la luz es una onda que sube y baja (como las olas del mar). Pero en matemáticas, es más fácil trabajar con números complejos (que tienen una parte "real" y una parte "imaginaria").
La analogía: Imagina que tienes una canción de rock. La parte "real" es la melodía que escuchas. Pero para analizarla mejor, necesitas separar las frecuencias graves de las agudas.
La lección: El autor dice que muchos olvidan que la luz real tiene dos "caras" (frecuencias positivas y negativas). Si solo miras una cara, pierdes información. Él nos enseña a usar la "Señal Analítica", que es como poner unos lentes especiales que solo te dejan ver la mitad útil de la onda, descartando el "ruido" redundante. Si no haces esto, tus cálculos de cómo la luz cambia de color o intensidad pueden salir mal (¡y a veces te falta un factor de 2 en la respuesta!).

2. El Transformado de Fourier: La traductora de idiomas

El problema: El Transformado de Fourier es la herramienta que convierte una señal de "tiempo" (cómo cambia la luz segundo a segundo) a "frecuencia" (qué colores tiene). Pero hay muchas formas de escribir esta fórmula (diferentes constantes, signos, etc.).
La analogía: Es como si en un país todos hablaran inglés, pero algunos escribieran las cartas de izquierda a derecha y otros de derecha a izquierda. Si le das una carta escrita de un modo a alguien que lee al revés, entenderá lo contrario de lo que dijiste.
La lección: El autor advierte que en óptica usamos una convención específica (como escribir de izquierda a derecha). Si mezclas las reglas de los matemáticos con las de los físicos, puedes terminar pensando que un pulso de luz se está acelerando cuando en realidad se está frenando. Es vital usar siempre el mismo "idioma" (convención) en toda la simulación.

3. El "DFT" y el efecto de los bordes (Aliasing)

El problema: Las computadoras no pueden ver el infinito; tienen una ventana de tiempo limitada. Si la señal es más larga que la ventana, la computadora la "corta" y crea ilusiones ópticas.
La analogía: Imagina que estás filmando un coche que pasa muy rápido, pero tu cámara solo tiene un marco de 1 metro. Si el coche es más largo, verás una parte al principio y otra al final, y pensarás que hay dos coches separados. O peor, si el coche da vueltas, parecerá que aparece y desaparece mágicamente.
La lección: Para evitar estas "fantasmas" (llamados aliasing), el autor explica cómo usar el "relleno con ceros" (zero-padding). Es como alargar la película para que el coche quepa entero sin cortarse, evitando que la computadora invente datos que no existen.

4. Análisis Tiempo-Frecuencia: El mapa de la luz

El problema: A veces la luz no es un color fijo; cambia de color mientras avanza (como un arcoíris que se desplaza).
La analogía: Si escuchas una canción, puedes saber qué notas hay (frecuencia) o cuándo suenan (tiempo), pero ¿cómo sabes qué nota suena exactamente en el segundo 3.5? Un análisis normal te da una de las dos cosas, no las dos a la vez.
La lección: El autor presenta una técnica llamada "Espectrograma de Alta Resolución". Imagina que tienes muchas cámaras con diferentes lentes: unas enfocan el tiempo (rápidas pero borrosas en color) y otras enfocan el color (precisas pero lentas en tiempo). Esta técnica combina todas las fotos para crear un mapa perfecto que te dice exactamente qué color tiene la luz en cada milisegundo, sin perder detalle.

5. La Transformada de Hankel: El giro de la luz

El problema: Cuando la luz viaja por un tubo o un láser, suele tener forma redonda (simetría radial). Usar herramientas cuadradas (como el Fourier normal) para cosas redondas es ineficiente, como intentar medir el perímetro de una pizza usando una regla cuadrada.
La analogía: Si quieres cortar una pizza, no usas un cuchillo que va en línea recta de lado a lado; usas un movimiento circular.
La lección: El autor introduce la Transformada de Hankel, que es la herramienta matemática diseñada específicamente para cosas redondas. Además, propone una versión mejorada y más rápida (llamada FHATHA) que corrige errores de energía.

El truco: A veces, al hacer los cálculos rápidos, la "pizza" pierde un poco de queso (energía). El autor sugiere un pequeño "ajuste final" para asegurar que la cantidad de queso sea la misma al principio y al final, o bien, ignorar ese ajuste si el error es tan pequeño que no importa, ahorrando tiempo de computadora.

En resumen

Este artículo es una guía para que los científicos de la luz no se pierdan en su propio laberinto matemático.

Usa los lentes correctos (Señal Analítica).
Habla el mismo idioma (Convención de Fourier).
No cortas la película (Evita el aliasing).
Mira el mapa completo (Espectrogramas de alta resolución).
Usa herramientas redondas para cosas redondas (Transformada de Hankel mejorada).

El autor incluso ha dejado el código de sus recetas (sus programas de computadora) en internet para que cualquiera pueda cocinar (simular) luz perfecta sin quemarse los dedos. ¡Es una guía práctica para que la física funcione de verdad!

A continuación presento un resumen técnico detallado del artículo "Tutorial of Fourier and Hankel transforms for ultrafast optics" de Yi-Hao Chen, estructurado según los puntos solicitados.

1. Problema Identificado

El artículo aborda una serie de malentendidos y deficiencias críticas en la simulación numérica de la óptica ultrarrápida, específicamente en el uso de las transformadas de Fourier y Hankel. Los problemas principales identificados son:

Ambigüedad en la Señal Analítica: Existe confusión frecuente sobre la descomposición de campos reales en señales analíticas (frecuencias positivas) y sus conjugados complejos. La omisión de la señal analítica correcta lleva a errores en la conservación de frecuencias (especialmente frecuencias negativas) y a factores de error de 2 en cálculos no lineales (ej. dispersión Raman, mezcla de cuatro ondas).
Convenciones de Transformada de Fourier: La falta de una convención unificada y la selección arbitraria de constantes en las transformadas de Fourier (FT) generan ecuaciones engañosas. El uso incorrecto de la convención (física vs. matemática) en simulaciones numéricas puede invertir el chirp de un pulso o alterar la evolución de la fase, llevando a resultados físicos incorrectos que no se detectan fácilmente con verificaciones superficiales (como la suavidad del perfil temporal).
Limitaciones de la Transformada de Hankel Numérica: Para sistemas con simetría radial (comunes en óptica de campo completo), la Transformada de Hankel es esencial. Sin embargo, los métodos existentes, como la Transformada de Hankel Rápida de Alta Precisión (FHATHA), sufren de problemas de conservación de energía y errores significativos al omitir el centro espacial ( $r=0$ ), lo cual es crítico para perfiles gaussianos.
Falta de Recursos Integrales: No existía una introducción completa que derivara las formulaciones generales de las relaciones de Fourier y Hankel, incluyendo sus constantes y su implementación discreta (DFT/DHT) con unidades físicas correctas.

2. Metodología

El autor emplea un enfoque pedagógico y derivativo riguroso, combinando teoría analítica con implementación numérica práctica:

Derivación Analítica Parametrizada: Se definen las transformadas de Fourier y Hankel con constantes genéricas ( $C_F, C_{IF}, C_H, C_{IH}$ ) y signos ( $c_s$ ) para demostrar cómo las convenciones afectan las ecuaciones físicas. Se derivan teoremas de convolución, Parseval y relaciones de autocorrelación bajo estas convenciones.
Análisis de Señal Analítica: Se establece la necesidad de partir siempre de una señal real y descomponerla en su señal analítica (parte de frecuencia positiva) antes de aplicar aproximaciones de envolvente lentamente variable (SVEA) o ecuaciones de propagación.
Desarrollo de Algoritmos Numéricos:
- Se analiza el DFT (Transformada Discreta de Fourier), abordando el uso de fftshift/ifftshift, el desenrollado de fase (phase unwrapping) y la convolución circular vs. lineal (incluyendo el relleno de ceros o zero-padding para evitar aliasing).
- Se propone una mejora al algoritmo FHATHA (Fast Hankel Transform of High Accuracy) mediante la inclusión de la contribución del centro espacial ( $r=0$ ) y una estrategia de restauración de energía.
- Se compara el FHATHA mejorado con la Transformada Discreta de Hankel (DHT) teórica.
Validación: Se utilizan ejemplos numéricos (pulsos chirpeados, generación de supercontinuo, propagación en fibras y telescopios) para verificar la precisión de las derivaciones y la estabilidad energética de los nuevos esquemas.

3. Contribuciones Clave

A. Clarificación de la Señal Analítica y SVEA

Se demuestra que la señal analítica es la única forma correcta de representar el campo complejo en óptica no lineal para evitar la pérdida de componentes de frecuencia negativa.
Se distingue claramente entre la aproximación de envolvente lentamente variable (SVEA) tradicional (que asume un solo modo espacial y frecuencia) y la señal analítica con desplazamiento de frecuencia/constante de propagación. Se aclara que los parámetros de desplazamiento ( $\omega_0, \beta_0$ ) pueden elegirse arbitrariamente para optimizar la velocidad de simulación, no necesariamente deben coincidir con el centro del pulso.

B. Convenciones de Transformada de Fourier

Se establece que en óptica (convención física $e^{-i\omega t}$ ), la transformada hacia el dominio espectral debe corresponder a la DFT inversa matemática (usualmente ifft en MATLAB) y viceversa.
Se advierte que mezclar convenciones sin aplicar el conjugado complejo al señal resulta en la propagación de un pulso con el chirp invertido (ej. un pulso up-chirped se comporta como down-chirped bajo la convención errónea).
Se proporcionan fórmulas precisas para convertir unidades físicas (W, J, nJ/THz) entre el dominio continuo y discreto, resolviendo problemas comunes de escalado en simulaciones.

C. Análisis Tiempo-Frecuencia de Alta Resolución

Se introduce un método para generar espectrogramas de alta resolución combinando múltiples espectrogramas con diferentes tamaños de ventana deslizante. Esto supera la limitación de la resolución tiempo-frecuencia (principio de incertidumbre) al identificar los máximos y mínimos normalizados de la energía en cada punto, eliminando artefactos de suavizado.

D. Mejora de la Transformada de Hankel (FHATHA)

Nueva Esquema Numérico: Se propone por primera vez una modificación al FHATHA que incluye explícitamente el valor del campo en el centro espacial ( $r=0$ ) y en el centro del espacio de frecuencias ( $k=0$ ). Esto corrige errores significativos en perfiles gaussianos donde el campo es máximo en el centro.
Restauración de Energía: Se introduce un paso de escalado uniforme al final de la transformación para restaurar la energía total, mitigando la no conservación de energía inherente a la discretización exponencial del FHATHA original.
Estrategia Híbrida: Se sugiere un esquema de simulación donde la restauración de energía se aplica solo en las transformadas inicial y final (entrada/salida del dominio espacial), permitiendo el uso del FHATHA estándar (sin restauración) durante la propagación intermedia para mayor eficiencia computacional sin acumulación de error.

4. Resultados

Precisión Mejorada: Las simulaciones que utilizan la señal analítica correcta y las convenciones de Fourier unificadas eliminan discrepancias de un factor de 2 en efectos no lineales y previenen la evolución física incorrecta de pulsos chirpeados.
Conservación de Energía: La versión modificada del FHATHA (con inclusión del centro y restauración de energía) muestra una conservación de energía casi perfecta y una reducción drástica del error cuadrático medio (RMSE) en comparación con el FHATHA original, especialmente con un número bajo de puntos de muestreo ( $N_\perp$ ).
Eficiencia Computacional: El nuevo esquema permite simulaciones de propagación de pulsos ultrarrápidos en 3D con simetría radial que son computacionalmente eficientes (usando FFT) pero tan precisas como métodos más costosos, superando a la DHT pura en velocidad sin sacrificar la conservación de energía física.
Herramientas Disponibles: El autor proporciona código fuente en GitHub para la implementación de estas técnicas, incluyendo herramientas para propagación en guías de onda y no guiada, así como funciones para espectrogramas de alta resolución.

5. Significado

Este tutorial es fundamental para la comunidad de óptica ultrarrápida por varias razones:

Corrección de Errores Sistemáticos: Proporciona la base teórica necesaria para evitar errores sutiles pero catastróficos en la simulación de fenómenos no lineales, asegurando que los resultados numéricos reflejen la física real.
Unificación de Convenciones: Ofrece un marco estandarizado para el manejo de constantes y unidades en transformadas discretas, facilitando la reproducibilidad y la comparación entre diferentes grupos de investigación.
Avance en Métodos Numéricos: La mejora propuesta para la Transformada de Hankel llena un vacío crítico en la simulación de sistemas de campo completo con simetría radial, permitiendo simulaciones más rápidas y precisas de fenómenos como la autoenfoque, la generación de supercontinuo y la compresión de pulsos.
Recurso Educativo: Sirve como una referencia exhaustiva que conecta la teoría matemática abstracta con la implementación práctica en software, llenando un vacío en la literatura existente que a menudo asume conocimientos previos o ignora detalles críticos de la implementación numérica.

En resumen, el artículo no solo educa sobre los fundamentos matemáticos, sino que propone soluciones prácticas y verificadas para mejorar la fiabilidad y eficiencia de las simulaciones en óptica ultrarrápida moderna.