🌀 nonlinear sciences

Asymptotic Analysis of discrete nonlinear localised modes in a Kagome lattice

This paper investigates nonlinear localized modes in a Kagome lattice with Klein-Gordon interactions by employing multiple scales asymptotic methods to derive various nonlinear Schrödinger reductions, including a novel coupled system at the flat band-upper band intersection, which is then analyzed via Lie symmetries and numerical simulations to reveal diverse solitary wave solutions.

Original authors: Jonathan AD Wattis, Pilar R Gordoa, Andrew Pickering

Published 2026-05-12

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Jonathan AD Wattis, Pilar R Gordoa, Andrew Pickering

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a giant, intricate dance floor made of triangles and hexagons, known as a Kagome lattice. On this floor, there are dancers (the nodes) connected by springs. Some springs connect dancers within the same small group, while others connect dancers to their neighbors in the next group over.

This paper is about what happens when these dancers start moving, but with a twist: the springs aren't perfectly stiff. They get a little "squishy" or "stiff" depending on how hard the dancers push (this is the nonlinear part). The researchers wanted to understand how waves of movement travel across this specific dance floor and if they could get stuck in one spot, forming a "breather" (a localized pulse of energy that stays put).

Here is the breakdown of their findings using simple analogies:

1. The Three "Floors" of Movement

When the dancers move, they can do so in three different "modes" or frequencies. The researchers found that the relationship between how fast they move and how they move (the dispersion relation) creates three distinct "surfaces" or layers of possibility:

The Flat Band: Imagine a perfectly flat, level floor. No matter how the dancers move, the energy doesn't spread out; it stays exactly where it is. This is a "flat" band because the frequency doesn't change with the direction of movement.
The Acoustic and Optical Bands: These are like two hilly landscapes. Usually, these hills are separate. However, under certain conditions, the peaks and valleys of these hills can touch.
- Dirac Points: Sometimes, the two hills touch at a single sharp point, forming a cone shape (like two ice cream cones touching tip-to-tip). This is called a Dirac point.
- Tangential Meeting: Sometimes, the flat floor touches the top of the highest hill smoothly, like a ball resting gently on a table.

2. The "Zoom Lens" Technique

The dancers are moving very slowly and gently (small amplitude). To understand the complex choreography, the researchers used a mathematical "zoom lens" called asymptotic analysis.

Instead of tracking every single tiny wiggle of every dancer, they looked at the "envelope" or the overall shape of the wave.
This allowed them to simplify the messy, complex rules of the dance floor into a much cleaner set of rules known as the Nonlinear Schrödinger (NLS) equations. Think of this as translating a complex, chaotic jazz improvisation into a simple, readable sheet of music.

3. The New Discovery: A Coupled Dance

Most of the time, the researchers found that the wave behavior could be described by a single, standard equation (like a solo dancer). However, they discovered something new and exciting at the point where the Flat Band meets the top of the Optical Hill (the tangential meeting).

Here, the wave behavior couldn't be described by just one dancer. It required a coupled system—two dancers who are constantly influencing each other.

They derived a new, novel system of equations that describes how these two "dancers" (mathematical variables) interact.
This system is unique and hadn't been seen before in this specific context.

4. Finding the "Vortex" Solitary Wave

Using advanced mathematical tools (Lie symmetries), the researchers looked for special solutions to this new coupled system. They found a solution that looks like a vortex.

Imagine a whirlpool in a pond, but instead of water, it's a pulse of energy spinning in a specific pattern across the lattice.
They created a computer simulation to watch this "whirlpool" in action. The simulation showed that this energy pulse is robust. Even after a long time (over 80 cycles of movement), it didn't fall apart or spread out; it stayed together, just like a stable whirlpool.

5. The Stability Question

While the wave looked stable in the simulation, the researchers are honest about a remaining mystery. In physics, there is a famous rule (the Vakhitov-Kolokolov criterion) that often predicts when such waves will collapse. Their math suggests this new wave might be unstable according to that rule, but because the simulation showed it lasting a long time, the final verdict on its stability is still "open." It's like seeing a spinning top that looks like it should fall over, but keeps spinning for a surprisingly long time.

Summary

In short, the paper takes a complex, triangular grid of connected springs, simplifies the math to find the "rules of the game," and discovers that at a specific intersection of energy levels, the rules change. Instead of a simple wave, you get a complex, interacting pair of waves that can form a stable, spinning "vortex" of energy. They proved this mathematically and showed it works in a computer simulation.

Technical Summary: Asymptotic Analysis of Discrete Nonlinear Localised Modes in a Kagome Lattice

Problem Statement
This paper investigates small-amplitude, weakly nonlinear oscillations within a discrete Kagome lattice governed by Klein-Gordon interactions. The model consists of a lattice with hexagonal symmetry composed of triangles and hexagons, where each unit cell contains three nodes (A, B, and C). Each node possesses a single scalar degree of freedom. The system features linear nearest-neighbor interactions with distinct coupling strengths within the unit cell ( $\gamma$ ) and between cells ( $\lambda$ ), alongside a nonlinear on-site potential characterized by quartic anharmonicity ( $\beta$ ). The primary objective is to derive the governing equations for envelope solutions (breathers) and to analyze the dispersion relation, specifically focusing on regimes where dispersion surfaces meet, such as Dirac points and tangential intersections involving a flat band.

Methodology
The authors employ multiple scales asymptotic analysis to reduce the discrete equations of motion to continuous partial differential equations (PDEs).

Dispersion Analysis: The linearized system is analyzed to determine the dispersion relation, identifying three bands: a flat band (frequency independent of wavenumber) and two wavenumber-dependent acoustic and optical bands. The analysis identifies critical points where these bands meet, including Dirac points (where $\gamma = \lambda$ ) and tangential intersections (specifically at the global maximum where the optical band meets the flat band at $k=l=0$ ).
Asymptotic Expansion: A multiscale ansatz is introduced where the solution is expanded in powers of a small parameter $\epsilon$ . The spatial and temporal variables are rescaled into slow variables ( $x, y, \tau, T$ ).
Reduction to NLS Systems: By applying the Fredholm alternative at successive orders of $\epsilon$ $ϵ$ , the authors derive solvability conditions.
- For generic points and simple band extrema, the system reduces to a single 2+1 dimensional Nonlinear Schrödinger (NLS) equation.
- For the specific case of the tangential intersection between the flat band and the optical band (Case 7-8), the leading-order kernel is two-dimensional. This necessitates a novel reduction to a coupled system of 2+1 dimensional NLS equations.
Symmetry Analysis: The derived coupled NLS system is analyzed using Lie point symmetries to find similarity reductions that transform the PDEs into Ordinary Differential Equations (ODEs).
Numerical Simulation: The resulting solitary wave solutions are simulated on a discrete lattice to test their robustness and compare them with the asymptotic predictions.

Key Contributions and Results

Dispersion Relation Characterization: The study confirms the existence of three dispersion surfaces. It details the conditions under which Dirac points form (requiring $\gamma = \lambda$ ) and identifies a unique tangential meeting point between the flat band and the optical band at the origin of the Brillouin zone.
Derivation of a Novel Coupled NLS System: The paper's central theoretical contribution is the derivation of a new system of coupled 2+1 dimensional NLS equations (Equation 2.47) arising specifically from the tangential intersection of the flat and optical bands. This system involves two complex envelopes, $F$ and $P$ , and includes cross-derivative terms ( $F_{xy}, P_{xy}$ ) and specific cubic nonlinear couplings.
Conserved Quantities: The authors identify conserved quantities for this new coupled system, corresponding to the $L^2$ norm (charge/mass) and energy.
Similarity Reductions and Solitary Waves: Using Lie symmetry analysis, the coupled PDE system is reduced to ODEs. The authors find a specific similarity reduction leading to a vortex-like solitary wave solution. Unlike standard Townes solitons found in single NLS equations, these solutions exhibit a vector field structure where the three lattice components ( $q_A, q_B, q_C$ ) are rotations of one another by $120^\circ$ .
Approximate Solutions: In the weakly nonlinear limit, an approximate analytical solution is constructed using a Bessel function expansion ( $J_1$ ) with a slowly varying amplitude, derived via a secularity condition.
Numerical Validation: Simulations on a $160 \times 80$ lattice demonstrate that the derived mode is robust, persisting for over 80 carrier wave oscillations. The energy distribution remains largely stable, though the authors note that the algebraic decay of the solution ( $u \sim r^{-2}$ ) requires very large lattices for perfect numerical convergence.

Significance and Claims
The paper claims to provide a rigorous asymptotic derivation of nonlinear localized modes in Kagome lattices, extending previous work on square, triangular, and honeycomb lattices.

Novelty: The authors emphasize that the coupled NLS system derived for the flat-band/optical-band intersection is a "novel system" not previously derived or studied in this specific form.
Mechanism: The work highlights the distinct behavior of modes at tangential band intersections compared to isolated Dirac points. While Dirac points lead to plane-wave-like dynamics or standard NLS reductions, the tangential intersection yields a more complex, coupled structure.
Generality: The authors suggest that similar phenomena (tangential meetings of dispersion surfaces) may occur in other lattice geometries, such as "super-Kagome" lattices or modified hexagonal lattices, implying that the derived coupled NLS systems could be relevant to a broader class of physical systems.
Limitations: The authors modestly state that the stability of the derived vortex-like solutions within the original discrete lattice remains an open question. While the Vakhitov-Kolokolov criterion suggests instability for the continuous NLS approximation, the stability in the discrete system is not definitively resolved. Furthermore, the current analysis is restricted to scalar unknowns; future work is proposed to address mechanical systems with vector displacements and FPUT-type interactions.

The paper concludes that while the derived modes are robust in short-term simulations, their long-term stability and the precise nature of their decay in the discrete lattice require further investigation.