New logarithmic power nonlinear Schrödinger equations with super-Gaussons

1. Problema e Contexto

A equação de Schrödinger não linear (NLS) é fundamental para modelar a propagação de ondas em meios dispersivos, com aplicações em óptica, condensados de Bose-Einstein (BECs) e física de plasmas. A versão clássica com não linearidade cúbica é amplamente estudada, mas a equação NLS logarítmica (log-NLS) ocupa um lugar especial devido às suas propriedades físicas e matemáticas únicas, como a separabilidade de subsistemas não interagentes e a existência de soluções localizadas exatas chamadas "Gaussons" (perfis gaussianos).

No entanto, o modelo log-NLS possui uma limitação significativa: suas soluções localizadas são estritamente gaussianas. Em muitos contextos físicos reais, os estados auto-aprisionados exibem perfis de "topo plano" (flat-top) ou super-gaussianos, caracterizados por uma amplitude quase constante no núcleo e decaimento abrupto nas bordas. O modelo log-NLS padrão não consegue capturar essa estrutura. O objetivo deste trabalho é superar essa limitação introduzindo uma nova classe de equações que admitam soluções exatas com perfis super-gaussianos.

2. Metodologia

O autor propõe e analisa uma nova equação denominada NLS logarítmico-potência (logp-NLS), definida pela equação:

$$i\partial_t \psi = -\Delta \psi - 4p^2 \left(\ln^2 |\psi|\right)^{\frac{p-1}{2p}} \left[\ln |\psi| + \left(1 + \frac{d-2}{2p}\right)\right] \psi$$

Onde $p \ge 1$ é um expoente que controla a "rigidez" da lei de compressibilidade não linear.

A metodologia empregada inclui:

• Derivação Analítica: A equação foi construída para admitir soluções estacionárias exatas na forma de super-Gaussons: $\psi(x, t) = e^{-r^{2p}} e^{i\theta}$. A forma da não linearidade foi determinada exigindo que o Laplaciano radial equilibre exatamente este ansatz.
• Estrutura Variacional e Hidrodinâmica: O modelo foi colocado no quadro Hamiltoniano padrão, derivando uma lei de pressão generalizada e uma lei de entalpia que dependem de uma potência logarítmica da densidade. Isso permite uma interpretação hidrodinâmica onde o meio possui uma compressibilidade efetiva que se torna infinitamente rígida no núcleo da soliton.
• Análise de Estabilidade Linear: Foi realizada uma análise de estabilidade linear perturbando a solução estacionária. O problema foi reduzido a um sistema de equações acopladas do tipo Schrödinger com potenciais anarmônicos e singulares.
• Simulações Numéricas: Simulações diretas em uma dimensão espacial foram realizadas para investigar a dinâmica não linear e colisões. Devido à singularidade da derivada da não linearidade em $|\psi|=1$, foi adotada uma estratégia de regularização (substituindo $\ln|\psi|$ por $\ln_\delta|\psi|$) e inicialização com amplitudes ligeiramente inferiores a 1 para garantir estabilidade numérica.

3. Principais Contribuições

• Introdução do Modelo logp-NLS: Proposição de uma generalização do log-NLS que introduz o parâmetro $p$, permitindo a transição contínua entre o Gausson gaussiano ($p=1$) e perfis de topo plano cada vez mais acentuados ($p > 1$).
• Soluções Exatas "Super-Gaussons": Demonstração da existência de soluções estacionárias exatas na forma de super-gaussianos, preenchendo a lacuna entre os modelos gaussianos e os perfis de topo plano observados experimentalmente.
• Interpretação Física da Rigidez: Estabelecimento de que o parâmetro $p$ governa a rigidez de uma lei de compressibilidade logarítmico-potência. A divergência da derivada da pressão próxima à densidade unitária explica a formação de um núcleo de densidade constante e bordas abruptas.
• Análise de Colisões Inelásticas: Investigação detalhada da dinâmica de colisões, revelando como o expoente $p$ e a velocidade de colisão influenciam a troca de energia entre modos translacionais e internos.

4. Resultados

• Estrutura das Soluções: À medida que $p$ aumenta, o núcleo da soliton torna-se mais plano e as bordas mais afiadas, aproximando-se de um perfil cilíndrico com bordas abruptas no limite $p \to \infty$.
• Estabilidade Linear: A análise espectral do operador linearizado mostra que todos os autovalores estão no eixo imaginário, indicando estabilidade linear para os Super-Gaussons considerados. O espectro é discreto, sem espectro contínuo.
• Dinâmica de Colisões:
  • As colisões são inelásticas, confirmando que o logp-NLS não é integrável (diferente da NLS cúbica).
  • Para baixas velocidades, a inelasticidade aumenta com $p$. A interação prolongada permite que os núcleos de topo plano se sobreponham, excitando modos internos de "respiração" (oscilação de largura) e radiação.
  • Para altas velocidades, a tendência se inverte: colisões com $p$ maior tornam-se relativamente mais elásticas. A rigidez extrema do núcleo (devido à alta $p$) impede a deformação durante o curto tempo de interação, fazendo com que as solitons se comportem como gotas quase incompressíveis que atravessam-se com menos perda de energia.
• Desafios Numéricos: Identificou-se que a singularidade da derivada da não linearidade em $|\psi|=1$ impõe restrições severas aos esquemas de integração temporal, exigindo regularização para simulações de longo prazo.

5. Significância e Implicações

Este trabalho expande significativamente o entendimento das equações de Schrödinger não lineares ao fornecer um modelo matemático rigoroso para solitons de topo plano.

• Teórico: O modelo conecta a física de solitons a leis de compressibilidade generalizadas, oferecendo uma estrutura Hamiltoniana para estados que antes eram apenas fenomenológicos.
• Aplicações Práticas: A capacidade de modelar perfis super-gaussianos é crucial para áreas onde tais formas são observadas, como:
  • Óptica Não Linear: Para pulsos de laser com perfis de intensidade uniformes.
  • Condensados de Bose-Einstein: Para descrever estados com densidade quase constante no núcleo.
  • Meios Granulares e Fluidos: Onde interações efetivas podem gerar forças logarítmicas ou de potência.

Em suma, o logp-NLS oferece uma ponte entre a teoria matemática de solitons e a realidade física de sistemas com respostas não lineares complexas, sugerindo novas direções para a investigação de estados localizados em meios não lineares.