Integrable perturbation theory for dark solitons of the defocusing nonlinear Schrödinger equation

Título: Teoria de Perturbação Integrável para Solitons Escuros da Equação de Schrödinger Não Linear Defocante

Autores: Nicholas J. Ossi, Barbara Prinari e Jianke Yang.

1. Problema e Contexto

O artigo aborda a limitação fundamental na aplicação da teoria de perturbação baseada em expansões de autofunções (método de inversão espectral) a solitons escuros na equação de Schrödinger não linear (NLS) defocante com fundo não nulo.

• O Desafio: Diferente dos solitons brilhantes (que decaem no infinito), os solitons escuros existem sobre um fundo constante não nulo ($q_0 \neq 0$). Perturbações externas (como amortecimento, dissipação ou auto-empinamento) não apenas alteram lentamente os parâmetros do soliton (velocidade, amplitude, fase, posição), mas também geram uma "prateleira" de radiação (radiation shelf) que se expande ao redor do soliton.
• A Lacuna Teórica: Trabalhos anteriores (como os de Konotop, Vekslerchik, Chen, Huang e outros) tentaram estender a teoria de perturbação integrável para este cenário. No entanto, essas abordagens falharam em prever corretamente a evolução adiabática de todos os parâmetros do soliton e não conseguiam descrever a formação da prateleira de radiação.
• A Hipótese de Falha: A existência da prateleira levou alguns pesquisadores (como Ablowitz et al.) a conjecturarem que as autofunções quadradas associadas ao soliton formavam uma base insuficiente para descrever a solução perturbada, questionando a existência de uma relação de fechamento (closure relation) válida para este problema.

2. Metodologia

Os autores revisitam e reformulam a teoria de perturbação baseada em autofunções quadradas, introduzindo correções cruciais na análise matemática do espectro contínuo.

• Prova de Completude: O trabalho fornece uma prova rigorosa da completude das autofunções quadradas para potenciais gerais com condições de contorno não nulas.
• Tratamento das Singularidades: A inovação central reside no tratamento correto das singularidades dos dados de espalhamento nos pontos de ramificação do espectro contínuo, $z = \pm q_0$.
  • Em trabalhos anteriores, as contribuições desses pontos de ramificação foram negligenciadas ou tratadas incorretamente.
  • Os autores demonstram que é necessário isolar as contribuições dos polos em $\pm q_0$ da parte integral da relação de fechamento, tratando o restante como um valor principal de Cauchy.
• Relação de Fechamento (Closure Relation): Derivam uma relação de fechamento completa para o caso de um único soliton (1-soliton) que inclui:
  1. Um termo de integral de valor principal sobre o espectro contínuo.
  2. Termos de resíduo associados aos autovalores discretos (o modo de Goldstone translacional e o modo generalizado de escala).
  3. Crucialmente: Termos de resíduo adicionais provenientes das singularidades em $z = \pm q_0$, que modificam os modos discretos efetivos usados nas condições de ortogonalidade.
• Teoria de Perturbação de Múltiplas Escalas: Utilizam uma expansão perturbativa onde os parâmetros do soliton ($k_1, \Lambda_1, x_1, \sigma_1$) evoluem em uma escala de tempo lenta $T = \epsilon t$. A correção de primeira ordem é expandida na base completa de autofunções quadradas.
• Supressão de Crescimento Secular: Impõem condições de ortogonalidade para eliminar termos que crescem linearmente no tempo (secular growth) na correção de primeira ordem. Isso gera equações diferenciais para a evolução dos parâmetros do soliton.

3. Contribuições Principais

1. Correção da Relação de Fechamento: Demonstram que a relação de fechamento anterior estava incompleta devido à falta de tratamento das singularidades em $\pm q_0$. A inclusão correta desses termos leva a modos discretos generalizados ($\tilde{Z}_1, \tilde{\Upsilon}_1$) que são ortogonais à função de perturbação de forma correta.
2. Predição da Prateleira de Radiação (Shelf): Mostram que a correção de primeira ordem contém um polo devido às singularidades em $\pm q_0$. A análise assintótica dessa integral revela a formação de uma prateleira de radiação em ambos os lados do soliton.
   • Derivam fórmulas explícitas para a altura ($h_{\pm}$) e o gradiente de fase ($\phi_{\xi}$) da prateleira.
   • Demonstram que a existência da prateleira está intrinsecamente ligada à evolução da fase do soliton ($\sigma_1$). Se a perturbação não gera singularidades (ex: auto-empinamento), não há prateleira.
3. Equações de Evolução Completa: Determinam as equações de evolução adiabática para os quatro parâmetros do soliton:
   • Velocidade ($k_1$)
   • Amplitude/Profundidade ($\Lambda_1$)
   • Posição do centro ($x_1$)
   • Fase global ($\sigma_1$)
   • Nota: A evolução da posição do centro e da fase requer uma condição de ortogonalidade específica que relaciona $\sigma_1$ e $x_1$, a qual estava ausente em trabalhos anteriores.
4. Validação Numérica: Todas as previsões teóricas são validadas através de simulações numéricas diretas da NLS perturbada, utilizando o método ETDRK4.

4. Resultados Chave

• Acordo com a Teoria Direta: Os resultados para amplitude, velocidade e profundidade do soliton concordam com a teoria de perturbação direta (baseada em leis de conservação perturbadas) de Ablowitz et al. [2] até a ordem $\mathcal{O}(\epsilon)$.
• Discrepância no Centro e Fase:
  • A teoria baseada em autofunções quadradas (neste trabalho) e a teoria direta de Ablowitz concordam em $\mathcal{O}(\epsilon)$, mas divergem na evolução de longo prazo do centro ($x_1$) e da fase ($\sigma_1$).
  • A divergência ocorre porque a evolução do centro depende de termos de ordem $\mathcal{O}(\epsilon^2)$ (derivados do Hamiltoniano em trabalhos anteriores) ou de uma dependência em uma escala de tempo ainda mais lenta ($\epsilon^2 t$).
  • O trabalho explica que a falha em capturar termos quadráticos na evolução do centro é uma limitação intrínseca da expansão de primeira ordem quando a velocidade é integrada no tempo para definir o referencial móvel.
• Casos de Aplicação:
  • Amortecimento Linear e Não Linear: Preveem corretamente a redução da amplitude e da velocidade, e a formação de prateleiras elevadas.
  • Dissipação: Preveem que a amplitude e velocidade permanecem constantes, mas o centro e a fase evoluem. A prateleira é deprimida.
  • Auto-empinamento (Self-steepening): Mostram que, para esta perturbação específica, as condições de ortogonalidade anulam a geração de prateleira (pois o produto interno com os modos em $\pm q_0$ é zero), e apenas a posição do centro se desloca linearmente. Isso contradiz trabalhos anteriores que previam evolução de amplitude.

5. Significado e Conclusão

Este trabalho resolve um problema aberto na teoria de perturbação de sistemas integráveis com fundo não nulo.

• Validação da Base de Autofunções: Refuta a conjectura de que as autofunções quadradas são uma base insuficiente. O problema não era a falta de completude, mas sim a aplicação incorreta da relação de fechamento (ignorando as singularidades de borda).
• Unificação de Fenômenos: Estabelece uma conexão direta entre a evolução dos parâmetros do soliton e a formação de estruturas de radiação (prateleiras), algo que teorias puramente baseadas em leis de conservação não descrevem tão naturalmente.
• Limitações e Futuro: O trabalho destaca que, para previsões de alta precisão em tempos longos ($t \sim 1/\epsilon$), é necessário incluir correções de segunda ordem ou escalas de tempo adicionais ($\epsilon^2 t$), uma limitação que afeta tanto a teoria de perturbação integrável quanto a direta.
• Impacto: A metodologia desenvolvida pode ser generalizada para sistemas acoplados (como o sistema Manakov) e sistemas discretos, abrindo novas vias de pesquisa para o estudo de ondas solitárias em meios não integráveis.

Em resumo, o artigo restaura a validade e o poder da teoria de perturbação baseada em autofunções quadradas para solitons escuros, corrigindo erros históricos na manipulação das singularidades espectrais e fornecendo uma descrição completa e verificável numericamente da dinâmica perturbada, incluindo a formação de prateleiras de radiação.