Two-dimensional paradigm for symmetry breaking: the nonlinear Schrödinger equation with a four-well potential

We study the existence and stability of localized modes in the two-dimensional (2D) nonlinear Schrödinger/Gross-Pitaevskii (NLS/GP) equation with a symmetric four-well potential. Using the corresponding four-mode approximation, we trace the parametric evolution of the trapped stationary modes, starting from the linear limit, and thus derive a complete bifurcation diagram for families of the stationary modes. This provides the picture of spontaneous symmetry breaking in the fundamental 2D setting. In a broad parameter region, the predictions based on the four-mode decomposition are found to be in good agreement with full numerical solutions of the NLS/GP equation. Stability properties of the stationary states coincide with those suggested by the corresponding discrete model in the large-amplitude limit. The dynamics of unstable modes is explored by means of direct simulations. Finally, in addition to the full analysis for the case of the self-attractive nonlinearity, the bifurcation diagram for the case of self-repulsion is briefly considered too.