Analysis of multistability in discrete quantum droplets and bubbles

This study investigates the existence and stability of localized states in the discrete nonlinear Schr\"odinger (DNLS) equation with quadratic and cubic nonlinearities, describing the so-called quantum droplets and bubbles. Those states exist within an interval known as the pinning region, as we vary a control parameter. Within the interval, multistable states are connected through multiple hysteresis, called homoclinic snaking. In particular, we explore its mechanism and consider two limiting cases of coupling strength: weak (anti-continuum) and strong (continuum) limits. We employ an asymptotic and a variational method for the weak and strong coupling limits, respectively, to capture the pinning region's width. The width exhibits an algebraic and an exponentially small dependence on the coupling constant for the weak and strong coupling, respectively. This finding is supported by both analytical and numerical results, which show excellent agreement. We also consider the modulational instability of spatially uniform solutions. Our work sheds light on the intricate interplay between multistability and homoclinic snaking in discrete quantum systems, paving the way for further exploration of complex nonlinear phenomena in this context.
Comment: Chaos, Solitons and Fractals