We study nonlinear surface modes in two-dimensional anisotropic periodic photonic lattices and demonstrate that, in a sharp contrast to one-dimensional discrete surface solitons, the mode threshold power is lower at the surface, and two-dimensional discrete solitons can be generated easier near the lattice corners and edges. We analyze the crossover between effectively one-and two-dimensional regimes of the surface-mediated beam localization in the lattice.
PACS numbers:Surface modes have been studied in different branches of physics; in guided wave optics surface states were predicted to exist at interfaces separating periodic and homogeneous dielectric media [1]. The interest in studying surface waves has been renewed recently because the interplay of discreteness and nonlinearity can facilitate the formation of discrete surface solitons [2,3] at the edge of the waveguide array. That can be understood as the localization of a discrete optical soliton near the surface [4] for powers exceeding a certain threshold value, for which the repulsive effect of the surface is balanced. A similar effect of light localization near the edge of the waveguide array and the formation of surface gap solitons have been predicted and observed for defocusing nonlinear media [5,6].It is important to analyze how the properties of nonlinear surface waves are modified by the lattice dimensionality, and the first studies of different types of discrete surface solitons in two-dimensional lattices [7,8,9,10] revealed, in particular, that the presence of a surface increases the stability region for two-dimensional (2D) discrete solitons [10] and the threshold power for the edge surface state is slightly higher than that for the corner soliton [9].In this Letter we consider anisotropic semi-infinite twodimensional photonic lattices and study the crossover between one-and two-dimensional surface solitons emphasizing the crucial effect of the lattice dimensionality on the formation of surface solitons.We consider a semi-infinite 2D lattice [shown schematically in Fig.2(a) below], described by the system of coupled-mode equations for the normalized amplitudes u n,m [11,12],where ξ is the normalized propagation distance. We de- fine the lattice coupling as follows:where α characterizes the lattice anisotropy.Linear lattice waves of the form u n,m (ξ) = u 0 sin(kn) sin(qm) exp(iβξ) satisfy the dispersion relation β kq = 2(cos k + α cos q). In the nonlinear case, we look for localized stationary solutions of the form u n,m (ξ) = u n,m exp(iλξ), where the amplitudes u n,m are real, and λ is the nonlinear propagation constant. For a given λ, localized solutions are found in a 15 × 15 lattice by using the Newton-Raphson method.We calculate the power threshold P th that characterizes the discrete solitons in 2D lattices [13]. We study three different modes: corner [ Fig. 1(a)], edge surface [ Fig. 1(b,c)] and central [ Fig. 1(d)] localized modes. The corner and edge modes represent 2D surface localized