On the independence of the slow and fast scales in multiple-scale expansions, with application to Van der Pol's equation

Planteamiento del Problema
El método de las múltiples escalas (MMS) es una técnica de perturbación estándar utilizada para resolver ecuaciones diferenciales autónomas que presentan escalas de tiempo ampliamente separadas, como los osciladores débilmente no lineales. Una instrucción fundamental del MMS es tratar la escala de tiempo rápida ($\theta$) y la escala de tiempo lenta ($\tau$) como variables independientes durante la derivación, a pesar de que ambas se originan de la única variable física $t$. Aunque esta independencia formal produce excelentes aproximaciones numéricas, carece de una justificación rigurosa: la solución verdadera $Y(t; \epsilon)$ posee invarianza de traslación temporal, lo que implica que un desplazamiento en la variable de tiempo induciría un desplazamiento acoplado en ambas escalas. El artículo investiga la validez de la suposición de que $\theta$ y $\tau$ son estrictamente independientes, examinando específicamente si un desplazamiento en la fase rápida $\theta$ puede inducir un desplazamiento no trivial y exponencialmente pequeño en la escala lenta o una corrección de fase dependiente de la escala lenta. Los autores se centran en la ecuación de Van der Pol como un modelo prototípico para las bifurcaciones de Hopf y las oscilaciones débilmente no lineales.

Metodología
Los autores emplean un análisis asintótico riguroso que extiende el procedimiento estándar de MMS a órdenes arbitrariamente altos para exponer las limitaciones de la suposición de independencia. La metodología procede de la siguiente manera:

1. Expansión de alto orden: La solución se expande en potencias del parámetro pequeño $\epsilon$. El análisis rastrea los coeficientes de la expansión hasta órdenes grandes $j \gg 1$.
2. Divergencia y el ansatz de "factorial sobre potencia": Los autores demuestran que la serie asintótica generada por MMS es divergente. Utilizan un ansatz de "factorial sobre potencia" (inspirado en la teoría de Dingle) para caracterizar el comportamiento de los coeficientes $A_{j,k}$ cuando $j \to \infty$. Esto revela que la serie está dominada por singularidades en el plano complejo $\tau$, específicamente en $\tau_n = (2n+1)i\pi$.
3. Fenómeno de Stokes y truncamiento óptimo: El análisis identifica líneas de Stokes que emanan de estas singularidades. Al truncar óptimamente la serie divergente cerca de su término más pequeño, los autores aíslan el resto exponencialmente pequeño $R$. Este resto representa la contribución "más allá de todos los órdenes" que el MMS estándar ignora.
4. Emparejamiento y cambio de Stokes (Stokes Switching): El resto se analiza a través del cruce de la línea de Stokes con el eje real. Los autores emparejan el comportamiento de términos tardíos de la serie con el comportamiento local cerca de la singularidad para determinar la activación de términos exponencialmente pequeños. Esto implica resolver ecuaciones linealizadas para el resto para encontrar la forma específica del desplazamiento de fase y la traslación de la escala lenta.
5. Verificación numérica: Las predicciones teóricas se validan mediante simulaciones numéricas de la ecuación de Van der Pol. Los autores computan el desplazamiento de fase del ciclo límite relativo a la fase inicial para varios valores de $\epsilon$ y comparan la amplitud de este desplazamiento contra las leyes de escala analíticas.

Contribuciones Clave y Resultados
El resultado principal del artículo es la derivación de un desplazamiento de fase $\Delta\theta$ no trivial y exponencialmente pequeño que rompe la supuesta independencia de las escalas lenta y rápida.

• Corrección a la simetría de transformación: El artículo invalida la transformación estándar de MMS $y(\theta, \tau) \to y(\theta + \phi, \tau)$ como una descripción completa de la simetría del sistema. En su lugar, la transformación correcta incluye una corrección exponencialmente pequeña:
  $$y(\theta, \tau) \to y(\theta + \phi + \Delta\theta(\tau, \phi, \epsilon), \tau + \Delta\tau(\tau, \phi, \epsilon))$$
  donde $\Delta\theta$ y $\Delta\tau$ son funciones exponencialmente pequeñas de $\epsilon$.
• Forma explícita del desplazamiento de fase: Para la ecuación de Van der Pol, el desplazamiento de fase $\Delta\theta$ en el límite de tiempo largo ($\tau \to \infty$) se encuentra como:
  $$\Delta\theta \sim 4\pi\Lambda e^{-2\pi/\epsilon} \cos(2\phi)$$
  Aquí, $\Lambda \approx 0.023080$ es una constante numérica derivada de las relaciones de recurrencia, y $\phi$ es la fase inicial. El factor $e^{-2\pi/\epsilon}$ indica que el efecto es exponencialmente pequeño, desvaneciéndose más rápido que cualquier potencia de $\epsilon$.
• Papel de la no linealidad y los armónicos: La magnitud del efecto está determinada por el espaciamiento de los armónicos generados por la no linealidad. Debido a que la ecuación de Van der Pol contiene una no linealidad cúbica, solo se generan armónicos impares, lo que resulta en un espaciamiento armónico de $\Delta k = 2$ en el espacio de Fourier. Esto da como resultado el factor exponencial $e^{-2\pi/\epsilon}$. Los autores señalan que si estuvieran presentes no linealidades cuadráticas (generando armónicos pares con $\Delta k = 1$), el factor sería $e^{-\pi/\epsilon}$, haciendo que el efecto fuera significativamente mayor.
• Confirmación numérica: Las simulaciones numéricas confirman la existencia del desplazamiento de fase y su dependencia de la fase inicial $\phi$ como $\cos(2\phi)$. La amplitud de la oscilación en el desplazamiento de fase coincide con la ley de escala predicha $\Delta(\epsilon) \approx 0.58 e^{-2\pi/\epsilon}$ para valores de $\epsilon$ de hasta 0.35.

Significado y Reivindicaciones
El artículo pretende dilucidar el mecanismo preciso mediante el cual las escalas lenta y rápida interactúan en problemas diferenciales autónomos, un fenómeno que a menudo se pasa por alto en la teoría de perturbaciones estándar.

• Validez de MMS: Los autores argumentan que la suposición de independencia de escalas en MMS no es fundamentalmente "errónea", sino más bien "benigna" en muchos contextos. La interacción entre escalas es exponencialmente débil y ocurre solo en la vecindad de las líneas de Stokes (donde el crecimiento de la amplitud es más rápido). Lejos de estas líneas, las escalas se desacoplan efectivamente, justificando el éxito práctico de MMS.
• Origen No Lineal: El estudio establece que el acoplamiento entre escalas es un efecto puramente no lineal. En sistemas lineales, la solución es una combinación simple de exponenciales donde las escalas permanecen estrictamente factorizadas; la ruptura de la independencia surge específicamente de la generación no lineal de armónicos superiores.
• Generalidad: Aunque el cálculo se realiza para la ecuación de Van der Pol, los autores sugieren que el resultado es esencialmente general. La forma de la corrección depende de los armónicos específicos generados por la no linealidad, lo que implica que el método puede aplicarse a otras bifurcaciones de Hopf y paquetes de ondas débilmente no lineales.
• Asintótica Exponencial: El artículo sirve como una demostración detallada de cómo aplicar la asintótica exponencial (específicamente el truncamiento óptimo y el cambio de Stokes) a problemas de múltiples escalas, proporcionando una plantilla para analizar otros fenómenos similares de "más allá de todos los órdenes".