Snakes on Lieb lattice

We consider the discrete Allen–Cahn equation with cubic and quintic nonlinearity on the Lieb lattice. We study localized nonlinear solutions of the system that have linear multistability and hysteresis in their bifurcation diagram. In this work, we investigate the system’s homoclinic snaking, i.e. snaking-like structure of the bifurcation diagram, particularly the effect of the lattice type. Numerical continuation using a pseudo-arclength method is used to obtain localized solutions along the bifurcation diagram. We then develop an active-cell approximation to classify the type of solution at the turning points, which gives good agreement with the numerical results when the sites are weakly coupled. Time dynamics of localized solutions inside and outside the pinning region is also discussed.