Snakes and ghosts in a parity-time-symmetric chain of dimers

We consider linearly coupled discrete nonlinear Schrödinger equations with gain and loss terms and with a cubic-quintic nonlinearity. The system models a parity-time (PT)-symmetric coupler composed by a chain of dimers. We study uniform states and site-centered and bond-centered spatially localized solutions and present that each solution has a symmetric and antisymmetric configuration between the arms. The symmetric solutions can become unstable due to bifurcations of asymmetric ones, that are called ghost states, because they exist only when an otherwise real propagation constant is taken to be complex valued. When a parameter is varied, the resulting bifurcation diagrams for the existence of standing localized solutions have a snaking behavior. The critical gain and loss coefficient above which the PT symmetry is broken corresponds to the condition when bifurcation diagrams of symmetric and antisymmetric states merge. Past the symmetry breaking, the system no longer has time-independent states. Nevertheless, equilibrium solutions can be analytically continued by defining a dual equation that leads to ghost states associated with growth or decay, that are also identified and examined here. We show that ghost localized states also exhibit snaking bifurcation diagrams. We analyze the width of the snaking region and provide asymptotic approximations in the limit of strong and weak coupling where good agreement is obtained.