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

Enunciado do Problema
O método das múltiplas escalas (MMS) é uma técnica de perturbação padrão utilizada para resolver equações diferenciais autônomas que apresentam escalas de tempo amplamente separadas, como osciladores fracamente não lineares. Uma instrução fundamental do MMS é tratar a escala de tempo rápida ($\theta$) e a escala de tempo lenta ($\tau$) como variáveis independentes durante a derivação, apesar de ambas originarem de uma única variável física $t$. Embora essa independência formal resulte em excelentes aproximações numéricas, ela carece de uma justificativa rigorosa: a solução verdadeira $Y(t; \epsilon)$ possui invariância de translação temporal, o que implica que um deslocamento na variável temporal deve induzir um deslocamento acoplado em ambas as escalas. O artigo investiga a validade da suposição de que $\theta$ e $\tau$ são estritamente independentes, examinando especificamente se um deslocamento na fase rápida $\theta$ pode induzir um deslocamento não trivial e exponencialmente pequeno na escala lenta ou uma correção de fase dependente da escala lenta. Os autores focam na equação de Van der Pol como um modelo prototípico para bifurcações de Hopf e oscilações fracamente não lineares.

Metodologia
Os autores empregam uma análise assintótica rigorosa que estende o procedimento padrão do MMS para ordens arbitrariamente altas para expor as limitações da suposição de independência. A metodologia procede da seguinte forma:

1. Expansão de Alta Ordem: A solução é expandida em potências do parâmetro pequeno $\epsilon$. A análise rastreia os coeficientes da expansão até ordens elevadas $j \gg 1$.
2. Ansatz de Divergência e Fatorial-sobre-Potência: Os autores demonstram que a série assintótica gerada pelo MMS é divergente. Eles utilizam um ansatz de "fatorial-sobre-potência" (inspirado na teoria de Dingle) para caracterizar o comportamento dos coeficientes $A_{j,k}$ conforme $j \to \infty$. Isso revela que a série é dominada por singularidades no plano complexo $\tau$, especificamente em $\tau_n = (2n+1)i\pi$.
3. Fenômeno de Stokes e Truncamento Ótimo: A análise identifica linhas de Stokes emanando dessas singularidades. Ao truncar a série divergente de forma ótima próximo ao seu termo de menor valor, os autores isolam o resto exponencialmente pequeno $R$. Este resto representa a contribuição "além de todas as ordens" (beyond-all-orders) que o MMS padrão negligencia.
4. Correspondência e Alternância de Stokes (Stokes Switching): O resto é analisado através do cruzamento da linha de Stokes com o eixo real. Os autores fazem a correspondência entre o comportamento de termos tardios da série e o comportamento local próximo à singularidade para determinar o acionamento de termos exponencialmente pequenos. Isso envolve resolver equações linearizadas para o resto para encontrar a forma específica do deslocamento de fase e da translação da escala lenta.
5. Verificação Numérica: As previsões teóricas são validadas através de simulações numéricas da equação de Van der Pol. Os autores computam o deslocamento de fase do ciclo limite em relação à fase inicial para vários valores de $\epsilon$ e comparam a amplitude desse deslocamento contra as leis de escala analíticas.

Principais Contribuições e Resultados
O principal resultado do artigo é a derivação de um deslocamento de fase $\Delta\theta$ não trivial e exponencialmente pequeno que quebra a suposta independência das escalas lenta e rápida.

• Correção à Simetria de Transformação: O artigo invalida a transformação padrão do MMS $y(\theta, \tau) \to y(\theta + \phi, \tau)$ como uma descrição completa da simetria do sistema. Em vez disso, a transformação correta inclui uma correção exponencialmente pequena:
  $$y(\theta, \tau) \to y(\theta + \phi + \Delta\theta(\tau, \phi, \epsilon), \tau + \Delta\tau(\tau, \phi, \epsilon))$$
  onde $\Delta\theta$ e $\Delta\tau$ são funções exponencialmente pequenas de $\epsilon$.
• Forma Explícita do Deslocamento de Fase: Para a equação de Van der Pol, o deslocamento de fase $\Delta\theta$ no limite de longo tempo ($\tau \to \infty$) é encontrado como:
  $$\Delta\theta \sim 4\pi\Lambda e^{-2\pi/\epsilon} \cos(2\phi)$$
  Aqui, $\Lambda \approx 0.023080$ é uma constante numérica derivada das relações de recorrência, e $\phi$ é a fase inicial. O fator $e^{-2\pi/\epsilon}$ indica que o efeito é exponencialmente pequeno, desaparecendo mais rápido do que qualquer potência de $\epsilon$.
• Papel da Não Linearidade e Harmônicos: A magnitude do efeito é determinada pelo espaçamento dos harmônicos gerados pela não linearidade. Como a equação de Van der Pol contém uma não linearidade cúbica, apenas harmônicos ímpares são gerados, levando a um espaçamento harmônico de $\Delta k = 2$ no espaço de Fourier. Isso resulta no fator exponencial $e^{-2\pi/\epsilon}$. Os autores observam que, se não linearidades quadráticas estivessem presentes (gerando harmônicos pares com $\Delta k = 1$), o fator seria $e^{-\pi/\epsilon}$, tornando o efeito significativamente maior.
• Confirmação Numérica: Simulações numéricas confirmam a existência do deslocamento de fase e sua dependência da fase inicial $\phi$ como $\cos(2\phi)$. A amplitude da oscilação no deslocamento de fase condiz com a escala prevista $\Delta(\epsilon) \approx 0.58 e^{-2\pi/\epsilon}$ para valores de $\epsilon$ até 0.35.

Significância e Alegações
O artigo afirma elucidar o mecanismo preciso pelo qual as escalas lenta e rápida interagem em problemas diferenciais autônomos, um fenômeno frequentemente negligenciado na teoria de perturbação padrão.

• Validade do MMS: Os autores argumentam que a suposição de independência de escala no MMS não é fundamentalmente "errada", mas sim "benigna" em muitos contextos. A interação entre as escalas é exponencialmente fraca e ocorre apenas nas proximidades das linhas de Stokes (onde o crescimento da amplitude é mais rápido). Longe dessas linhas, as escalas efetivamente se desacoplam, justificando o sucesso prático do MMS.
• Origem Não Linear: O estudo estabelece que o acoplamento entre as escalas é um efeito puramente não linear. Em sistemas lineares, a solução é uma combinação simples de exponenciais onde as escalas permanecem estritamente fatorizadas; a quebra da independência surge especificamente da geração não linear de harmônicos superiores.
• Generalidade: Embora o cálculo seja realizado para a equação de Van der Pol, os autores sugerem que o resultado é essencialmente geral. A forma da correção depende dos harmônicos específicos gerados pela não linearidade, implicando que o método pode ser aplicado a outras bifurcações de Hopf e pacotes de ondas fracamente não lineares.
• Assintótica Exponencial: O artigo serve como uma demonstração detalhada de como aplicar a assintótica exponencial (especificamente o truncamento ótimo e a alternância de Stokes) a problemas de múltiplas escalas, fornecendo um modelo para analisar fenômenos semelhantes de "além de todas as ordens" em outros sistemas dinâmicos.