Global bifurcations and basin geometry of the nonlinear non-Hermitian skin effect

1. Problema e Contexto

O artigo investiga a interação entre não-Hermiticidade (sistemas abertos com ganho e perda) e não linearidade, focando especificamente no Efeito de Pele Não-Hermitiano (NHSE).

• Contexto: No NHSE linear, um grande número de autoestados do volume acumula-se nas bordas do sistema sob condições de contorno abertas, quebrando a correspondência volume-borda convencional.
• Desafio: Estudos anteriores sobre NHSE não linear (como o modelo de Hatano-Nelson não linear) focaram em bifurcações de Hopf supercríticas locais, onde um ponto fixo perde estabilidade e dá origem a um ciclo limite estável, levando a uma transição de pele para estados estendidos.
• Questão Central: O que acontece se a não reciprocidade dependente da amplitude for não monótona e saturada (ou seja, o ganho efetivo é suprimido em amplitudes grandes)? Quais estruturas globais no espaço de fases emergem e como isso afeta a coexistência de estados localizados (pele) e estendidos?

2. Metodologia

Os autores utilizam uma abordagem baseada em sistemas dinâmicos para analisar um modelo contínuo de Hatano-Nelson não linear:

1. Modelo Hamiltoniano:

   • Consideram uma equação de Schrödinger estacionária não linear com um termo não-Hermitiano modulada pela amplitude da onda: $F(|\psi|^2) = \gamma + a|\psi|^2 - b|\psi|^4$.
   • O termo $+a|\psi|^2$ aumenta a não reciprocidade em amplitudes moderadas, enquanto $-b|\psi|^4$ satura e suprime o ganho em amplitudes grandes.
   • Focam em estados estacionários reais com energia positiva $E > 0$.

2. Formulação como Fluxo no Espaço de Fases:

   • A equação diferencial parcial é tratada como uma equação diferencial ordinária (EDO) onde a coordenada espacial $x$ atua como um parâmetro de evolução temporal.
   • O sistema é mapeado em um fluxo planar $(\psi, v)$, onde $v = \partial_x \psi$.
   • Estados estacionários correspondem a trajetórias neste espaço de fases. A condição de contorno semi-infinita é parametrizada pela inclinação inicial $s = \partial_x \psi(0)$.

3. Análise de Bifurcação e Atratores:

   • Identificam pontos fixos e ciclos limites no espaço de fases e mapeiam sua estabilidade em função do parâmetro de não reciprocidade $\gamma$.
   • Correspondência Atrator-Localização:
     • Atratores na origem $(0,0)$ $\rightarrow$ Decaimento exponencial $\rightarrow$ Estado de Pele (Skin).
     • Ciclos limites estáveis $\rightarrow$ Oscilações limitadas com amplitude não nula $\rightarrow$ Estado Estendido.

4. Equação de Amplitude Média:

   • Derivam uma equação de amplitude média (método de média) para obter previsões analíticas fechadas para as ramificações dos ciclos limites e o limiar de bifurcação, validando os resultados numéricos.

5. Parâmetro de Ordem de Fração da Bacia:

   • Introduzem um novo parâmetro de ordem baseado na geometria das bacias de atração para quantificar a probabilidade de um estado inicial evoluir para o modo de pele versus o estado estendido.

3. Contribuições e Resultados Principais

A. Cenário de Bifurcação Global

O estudo revela um cenário global controlado por dois mecanismos de bifurcação distintos, gerando uma janela de coexistência finita ($\gamma_c < \gamma < 0$):

1. Bifurcação de Hopf Subcrítica ($\gamma = 0$): A origem perde estabilidade, mas um ciclo limite instável de pequena amplitude já existia para $\gamma < 0$.
2. Bifurcação Sela-Nó de Ciclos Limites (SNLC) em $\gamma = \gamma_c < 0$: Um par de ciclos limites (um estável externo e um instável interno) é criado.

Regimes Dinâmicos:

• $\gamma < \gamma_c$: Apenas a origem é atratora (apenas estados de pele).
• $\gamma > 0$: Apenas o ciclo limite estável é atrator (apenas estados estendidos).
• $\gamma_c < \gamma < 0$ (Janela de Coexistência): Ocorre bistabilidade. A origem (pele) e o ciclo limite externo (estendido) são ambos atratores estáveis. Eles são separados por uma separatriz formada pelo ciclo limite interno instável.

B. Mecanismo de Coexistência Distinto

Diferente de sistemas lineares onde a coexistência de estados localizados e estendidos depende de estruturas espectrais (como bordas de mobilidade), neste sistema não linear:

• Os dois estados coexistem para uma energia fixa $E$.
• A separação é determinada pela geometria da bacia de atração no espaço de fases, não por propriedades espectrais lineares.

C. Previsões Analíticas

A equação de amplitude média derivada fornece:

• Expressões fechadas para as amplitudes dos ciclos limites interno e externo.
• Uma previsão precisa para o limiar da bifurcação SNLC: $\gamma_c = -a^2 / (8b)$.
• Os resultados analíticos mostram excelente concordância com a continuação numérica direta do sistema.

D. Parâmetro de Ordem e Transição de Fase

• Os autores definem $p_{skin}(\gamma)$, a fração de trajetórias que terminam na origem (pele) para uma distribuição de inclinações iniciais.
• Na bifurcação SNLC ($\gamma = \gamma_c$), observa-se um salto do tipo primeira ordem em $p_{skin}$. Isso ocorre porque, no momento da criação do ciclo estável externo, uma fração finita de condições iniciais é imediatamente capturada por ele, alterando abruptamente a probabilidade de encontrar o estado de pele.

E. Fenômenos Físicos na Janela de Coexistência

1. Transientes Espaciais de Longa Duração: Trajetórias iniciadas próximas à separatriz (ciclo limite instável) podem "sombrear" esse ciclo instável por longas distâncias espaciais antes de eventualmente colapsar para o atrator estável. Isso cria estados transitórios que parecem estendidos por muito tempo, mas que eventualmente decaem (ou vice-versa), um fenômeno ausente em sistemas lineares.
2. Histerese: Devido à bistabilidade, o estado final do sistema depende da história de varredura do parâmetro $\gamma$:
   • Ao diminuir $\gamma$ (de estendido para pele), a transição ocorre em $\gamma_c$.
   • Ao aumentar $\gamma$ (de pele para estendido), a transição ocorre apenas em $\gamma = 0$ (onde a origem perde estabilidade).
   • Isso cria um laço de histerese e um efeito de memória.

4. Significado e Impacto

Este trabalho avança significativamente a compreensão dos sistemas não-Hermitianos não lineares ao:

• Demonstrar que a geometria global dos atratores e bacias é uma ferramenta poderosa e complementar aos conceitos espectrais lineares para classificar estados estacionários.
• Revelar que a saturação não linear pode induzir regimes de bistabilidade complexa com coexistência de modos de pele e estendidos, algo não previsto em modelos puramente lineares ou com não linearidades simples.
• Fornecer uma estrutura teórica para entender fenômenos de memória (histerese) e transientes de longa duração em dispositivos ópticos, de circuitos ou de átomos frios que exibem efeitos de pele não lineares.
• Estabelecer uma ponte entre a teoria de bifurcações de sistemas dinâmicos (Hopf subcrítico, SNLC) e a física da matéria condensada não-Hermitiana.

Em resumo, o artigo propõe uma nova lente para entender a localização em sistemas não-Hermitianos não lineares, onde a história do sistema e a geometria do espaço de fases são tão importantes quanto os parâmetros do Hamiltoniano.