Analysis of multistability in discrete quantum droplets and bubbles

This study investigates the existence and stability of localized states in the discrete nonlinear Schrödinger equation with quadratic and cubic nonlinearities, describing the so-called quantum droplets and bubbles. As we vary a control parameter, those states exist within an interval known as the pinning region. 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 asymptotic and variational methods for the weak and strong coupling limits, respectively, to capture the pinning region's width. The width exhibits algebraic and exponentially small dependencies 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.