Asymptotic Analysis of discrete nonlinear localised modes in a Kagome lattice

Resumo Técnico: Análise Assintótica de Modos Localizados Não Lineares Discretos em uma Rede Kagome

Enunciação do Problema
Este artigo investiga oscilações de pequena amplitude e fracamente não lineares dentro de uma rede Kagome discreta governada por interações de Klein-Gordon. O modelo consiste em uma rede com simetria hexagonal composta por triângulos e hexágonos, onde cada célula unitária contém três nós (A, B e C). Cada nó possui um único grau de liberdade escalar. O sistema apresenta interações lineares entre primeiros vizinhos com constantes de acoplamento distintas dentro da célula unitária ($\gamma$) e entre células ($\lambda$), juntamente com um potencial não linear no sítio caracterizado por anarmonicidade quártica ($\beta$). O objetivo principal é derivar as equações governantes para soluções de envelope (breathers) e analisar a relação de dispersão, focando especificamente em regimes onde superfícies de dispersão se encontram, como pontos de Dirac e interseções tangenciais envolvendo uma banda plana.

Metodologia
Os autores empregam análise assintótica de múltiplas escalas para reduzir as equações de movimento discretas a equações diferenciais parciais (EDPs) contínuas.

1. Análise de Dispersão: O sistema linearizado é analisado para determinar a relação de dispersão, identificando três bandas: uma banda plana (frequência independente do número de onda) e duas bandas acústica e óptica dependentes do número de onda. A análise identifica pontos críticos onde essas bandas se encontram, incluindo pontos de Dirac (onde $\gamma = \lambda$) e interseções tangenciais (especificamente no máximo global onde a banda óptica encontra a banda plana em $k=l=0$).
2. Expansão Assintótica: Um ansatz de múltiplas escalas é introduzido, onde a solução é expandida em potências de um parâmetro pequeno $\epsilon$. As variáveis espaciais e temporais são reescalonadas em variáveis lentas ($x, y, \tau, T$).
3. Redução para Sistemas NLS: Ao aplicar a alternativa de Fredholm em ordens sucessivas de $\epsilon$, os autores derivam condições de solvabilidade.
   • Para pontos genéricos e extremos simples de banda, o sistema reduz-se a uma única equação de Schrödinger Não Linear (NLS) 2+1 dimensional.
   • Para o caso específico da interseção tangencial entre a banda plana e a banda óptica (Caso 7-8), o núcleo de ordem principal é bidimensional. Isso exige uma redução novel para um sistema acoplado de equações NLS 2+1 dimensões.
4. Análise de Simetria: O sistema NLS acoplado derivado é analisado usando simetrias de Lie de ponto para encontrar reduções de similaridade que transformam as EDPs em Equações Diferenciais Ordinárias (EDOs).
5. Simulação Numérica: As soluções de onda solitária resultantes são simuladas em uma rede discreta para testar sua robustez e compará-las com as previsões assintóticas.

Principais Contribuições e Resultados

• Caracterização da Relação de Dispersão: O estudo confirma a existência de três superfícies de dispersão. Detalha as condições sob as quais pontos de Dirac se formam (requerendo $\gamma = \lambda$) e identifica um único ponto de encontro tangencial entre a banda plana e a banda óptica na origem da zona de Brillouin.
• Derivação de um Novo Sistema NLS Acoplado: A contribuição teórica central do artigo é a derivação de um novo sistema de equações NLS acopladas 2+1 dimensões (Equação 2.47) surgindo especificamente da interseção tangencial das bandas plana e óptica. Este sistema envolve dois envelopes complexos, $F$ e $P$, e inclui termos de derivadas cruzadas ($F_{xy}, P_{xy}$) e acoplamentos não lineares cúbicos específicos.
• Quantidades Conservadas: Os autores identificam quantidades conservadas para este novo sistema acoplado, correspondentes à norma $L^2$ (carga/massa) e energia.
• Reduções de Similaridade e Ondas Solitárias: Usando análise de simetria de Lie, o sistema acoplado de EDPs é reduzido a EDOs. Os autores encontram uma redução de similaridade específica levando a uma solução de onda solitária tipo vórtice. Diferentemente dos solitons de Townes padrão encontrados em equações NLS únicas, essas soluções exibem uma estrutura de campo vetorial onde os três componentes da rede ($q_A, q_B, q_C$) são rotações uns dos outros por $120^\circ$.
• Soluções Aproximadas: No limite fracamente não linear, uma solução analítica aproximada é construída usando uma expansão em função de Bessel ($J_1$) com amplitude variando lentamente, derivada via uma condição de secularidade.
• Validação Numérica: Simulações em uma rede $160 \times 80$ demonstram que o modo derivado é robusto, persistindo por mais de 80 oscilações da onda portadora. A distribuição de energia permanece largamente estável, embora os autores observem que o decaimento algébrico da solução ($u \sim r^{-2}$) requer redes muito grandes para uma convergência numérica perfeita.

Significado e Afirmações
O artigo afirma fornecer uma derivação assintótica rigorosa de modos localizados não lineares em redes Kagome, estendendo trabalhos anteriores sobre redes quadradas, triangulares e de favo de mel.

• Novidade: Os autores enfatizam que o sistema NLS acoplado derivado para a interseção banda plana/óptica é um "sistema novel" não derivado ou estudado anteriormente nesta forma específica.
• Mecanismo: O trabalho destaca o comportamento distinto dos modos em interseções de bandas tangenciais em comparação com pontos de Dirac isolados. Enquanto pontos de Dirac levam a dinâmicas semelhantes a ondas planas ou reduções NLS padrão, a interseção tangencial produz uma estrutura acoplada mais complexa.
• Generalidade: Os autores sugerem que fenômenos semelhantes (encontros tangenciais de superfícies de dispersão) podem ocorrer em outras geometrias de rede, como redes "super-Kagome" ou redes hexagonais modificadas, implicando que os sistemas NLS acoplados derivados podem ser relevantes para uma classe mais ampla de sistemas físicos.
• Limitações: Os autores afirmam modestamente que a estabilidade das soluções tipo vórtice derivadas dentro da rede discreta original permanece uma questão em aberto. Embora o critério de Vakhitov-Kolokolov sugira instabilidade para a aproximação NLS contínua, a estabilidade no sistema discreto não é definitivamente resolvida. Além disso, a análise atual é restrita a incógnitas escalares; trabalhos futuros são propostos para abordar sistemas mecânicos com deslocamentos vetoriais e interações do tipo FPUT.

O artigo conclui que, embora os modos derivados sejam robustos em simulações de curto prazo, sua estabilidade de longo prazo e a natureza precisa de seu decaimento na rede discreta requerem investigação adicional.