Dinámica de la propagación de solitones: bifurcación, caos y perspectivas cuantitativas sobre la ecuación de Camassa–Holm modificada

Ondas que se niegan a romper

Imagine una ola de mar que viaja kilómetros sin perder su forma, que se desliza junto a otras olas como si no pasara nada. Estas olas tenaces, llamadas solitones, aparecen no solo en el agua, sino también en plasmas, fibras ópticas e incluso en sistemas mecánicos. Este artículo explora cómo se desplazan dichas ondas y cómo a veces se tornan caóticas en un modelo matemático ampliamente utilizado para ondas en el agua, revelando patrones que podrían ayudar a ingenieros a predecir y controlar mejor comportamientos complejos de ondas en la naturaleza y la tecnología.

Un plano moderno para ondas en aguas someras

El estudio se centra en la ecuación de Camassa–Holm modificada (MCH), un modelo potente para ondas en canales de aguas someras y entornos físicos relacionados. Versiones anteriores de esta familia de ecuaciones ayudaron a explicar las sorprendentes “peakons”: ondas solitarias con una cresta afilada y puntiaguda que imitan el rompimiento real de las olas más de cerca que los modelos clásicos de libro de texto. A lo largo de los años, los investigadores han ajustado estas ecuaciones para captar comportamientos más ricos, desde pulsos suaves en forma de campana hasta ondas que se vuelven más empinadas y rompen. Aun así, obtener muchas soluciones exactas y matemáticamente limpias ha seguido siendo difícil, limitando nuestra capacidad para entender todas las formas posibles de las ondas y su estabilidad.

Una nueva herramienta para construir formas de onda exactas

Para abordar este desafío, los autores usan un esquema analítico refinado llamado método modificado de (G′/G)-expansión (MG′/GE). En términos sencillos, convierten la ecuación original de las ondas en espacio y tiempo en una única “coordenada viajera” que se mueve con la onda. Esto transforma una ecuación en derivadas parciales complicada en una ecuación diferencial ordinaria más manejable. El método MG′/GE, a continuación, supone una forma en serie flexible para la onda y determina los coeficientes equilibrando términos y resolviendo un conjunto de ecuaciones algebraicas. Este marco es versátil: ajustando unos pocos parámetros puede generar muchos tipos diferentes de soluciones dentro de una misma receta unificada, en lugar de necesitar un nuevo truco para cada nueva forma de onda.

Un zoológico de solitones: desde pulsos suaves hasta picos singulares

Usando este método, el artículo descubre alrededor de treinta soluciones de onda viajera distintas de la ecuación MCH. Entre ellas hay solitones brillantes (picos aislados sobre un fondo plano), solitones oscuros (hondos localizados en un nivel por lo demás uniforme) y solitones “singulares” más exóticos en los que la altura de la onda se vuelve extremadamente empinada o efectivamente no acotada en un punto. Hay solitones singulares simples y dobles, así como configuraciones múltiples brillantes, oscuras y singulares. Algunas soluciones se expresan mediante funciones hiperbólicas (ondas que parecen jorobas aisladas), otras mediante funciones trigonométricas (ondas más oscilatorias) y otras mediante formas racionales (con transiciones más abruptas). Superficies 3D detalladas, mapas de contorno, gráficos de densidad y de evolución temporal muestran cómo se desplazan estas estructuras, interactúan y concentran energía en espacio y tiempo.

Cuando el orden se vuelve caos

Más allá de enumerar formas de onda, los autores investigan cuán estables son esos patrones y cómo se comporta el sistema cuando se perturba suavemente. Reescriben la ecuación de la onda viajera como un sistema dinámico de dos variables y analizan sus puntos fijos, o estados de equilibrio, usando herramientas como matrices Jacobianas y valores propios. A medida que cambia un parámetro clave de velocidad, el sistema sufre una bifurcación de horquilla (pitchfork): un único equilibrio se divide en tres, algunos estables y otros inestables. Retratos en el plano de fases trazan las posibles trayectorias que puede seguir el sistema, mientras que diagramas de bifurcación muestran cómo el comportamiento a largo plazo cambia con los parámetros. El equipo añade luego diferentes tipos de “forzamientos” dependientes del tiempo —como términos seno, coseno, gaussianos e hiperbólicos— y sigue el movimiento resultante con retratos de fase, secciones de Poincaré, series temporales e ideas al estilo de los exponentes de Lyapunov. Dependiendo del forzamiento, el sistema puede asentarse en ciclos regulares, derivar hacia un movimiento cuasi-periódico tipo toro o volverse inestable y no acotado, ofreciendo una guía visual clara de cómo trenes de ondas estructurados pueden inclinarse hacia comportamientos complejos o caóticos.

Por qué importan estos resultados

Para lectores no especialistas, la conclusión es que este estudio ofrece una especie de “mapa y caja de herramientas” para una ecuación de ondas ampliamente utilizada. Los autores muestran cómo un único método analítico puede producir un rico catálogo de formas solitónicas exactas, confirmar que muchas de ellas son estables frente a pequeñas perturbaciones y señalar cuándo la dinámica subyacente es propensa a volverse irregular o caótica. Dado que las mismas estructuras matemáticas aparecen en ingeniería costera, comunicaciones por fibra óptica, dispositivos de plasma y otras tecnologías, estos hallazgos pueden ayudar a los investigadores a diseñar sistemas que aprovechen solitones robustos para transportar energía e información o que eviten regímenes de onda destructivos. El trabajo también prepara el terreno para extensiones futuras a situaciones más realistas, como materiales con memoria, influencias aleatorias o dimensiones superiores.

Cita: Alam, M.N., Ahmed, S.F., Ismael, H.F. et al. Dynamics of soliton propagation: bifurcation, chaos, and quantitative insights into the modified Camassa–Holm equation. Sci Rep 16, 7588 (2026). https://doi.org/10.1038/s41598-026-37010-2

Palabras clave: solitones, ondas en aguas someras, dinámica no lineal, caos y bifurcación, ecuación de Camassa–Holm