Discrete vector solitons in two-dimensional nonlinear waveguide arrays: solutions, stability, and dynamics.

We identify and investigate bimodal (vector) solitons in models of square-lattice arrays of nonlinear optical waveguides. These vector self-localized states are, in fact, self-induced channels in a nonlinear photonic-crystal matrix. Such two-dimensional discrete vector solitons are possible in waveguide arrays in which each element carries two light beams that are either orthogonally polarized or have different carrier wavelengths. Estimates of the physical parameters necessary to support such soliton solutions in waveguide arrays are given. Using Newton relaxation methods, we obtain stationary vector-soliton solutions, and examine their stability through the computation of linearized eigenvalues for small perturbations. Our results may also be applicable to other systems such as two-component Bose-Einstein condensates trapped in a two-dimensional optical lattice.