Borel-Padé exponential asymptotics for the discrete nonlinear Schrödinger model with next-to-nearest neighbour interactions

In the present work we study discrete nonlinear Schr{ö}dinger models combining nearest (NN) and next-nearest (NNN) neighbor interactions, motivated by experiments in waveguide arrays. While we consider the more experimentally accessible case of positive ratio $μ$ of NNN to NN interactions, we focus on the intriguing case of competing such interactions $(μ<0)$, where stationary states can exist only for $-1/4 < μ< 0$. We analyze the key eigenvalues for the stability of the pulse-like stationary (ground) states, and find that such modes depend exponentially on the coupling parameter , with suitable polynomial prefactors and corrections that we analyze in detail. Very good agreement of the resulting predictions is found with systematic numerical computations of the associated eigenvalues. This analysis uses Borel-Padé exponential asymptotics to determine Stokes multipliers in the solution; these multipliers cannot be obtained using standard matched asymptotic expansion approaches as they are hidden beyond all asymptotic orders, even near singular points. By using Borel-Padé methods near the singularity, we construct a general asymptotic template for studying parametric problems which require the calculation of subdominant Stokes multipliers.