Two-dimensional intra-band solitons in lattice potentials with local defects and self-focusing nonlinearity Jianhua Zeng1, ∗ and Boris A. Malomed2, † 1

arXiv:1305.1834v1 [physics.optics] 8 May 2013

State Key Laboratory of Low Dimensional Quantum Physics, Department of Physics, Tsinghua University, Beijing 100084, China 2 Department of Physical Electronics, School of Electrical Engineering, Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel It is commonly known that stable bright solitons in periodic potentials, which represent gratings in photonics/plasmonics, or optical lattices in quantum gases, exist either in the spectral semi-infinite gap (SIG) or in finite bandgaps. Using numerical methods, we demonstrate that, under the action of the cubic self-focusing nonlinearity, defects in the form of “holes” in two-dimensional (2D) lattices support continuous families of 2D solitons embedded into the first two Bloch bands of the respective linear spectrum, where solitons normally do not exist. The two families of the embedded defect solitons (EDSs) are found to be continuously linked by the branch of gap defect solitons (GDSs) populating the first finite bandgap. Further, the EDS branch traversing the first band links the GDS family with the branch of regular defect-supported solitons populating the SIG. Thus, we construct a continuous chain of regular, embedded, and gap-mode solitons (“superfamily”) threading the entire bandgap structure considered here. The EDSs are stable in the first Bloch band, and partly stable in the second. They exist with the norm exceeding a minimum value, hence they do not originate from linear defect modes. Further, we demonstrate that double, triple and quadruple lattice defects support stable dipole-mode solitons and vortices. OCIS codes: 190.5530, 020.1475, 230.5298, 350.7420.

I.

INTRODUCTION

A powerful toolbox for controlling the dynamics of photonic, plasmonic, and atomic waves is based on the use of periodic potentials, which can be induced by photonic lattices and crystals in optics [1], by gratings built into plasmonic waveguides (see original works [2–13] and review [14]), and by optical lattices (OLs) in Bose-Einstein condensates (BECs) [15–17] or degenerate Fermi gases [18]. Acting in the combination with material nonlinearity, periodic potentials give rise to various species of solitons. These include ordinary ones, existing in the semi-infinite gap (SIG) of the corresponding linear spectrum, and gap solitons populating finite bandgaps, see reviews [16, 17, 19, 20] and book [21]. Extended localized states, in the form of truncated nonlinear Bloch waves, which appear in finite bandgaps, have been observed and modeled too [22–24]. On the other hand, a specific species of solitary waves is known in the form of embedded solitons, which exist, counter-intuitively, inside spectral Bloch bands populated by radiation modes [25]-[39]. Although embedded solitons cannot exist, generically, in continuous families, due to the resonance with radiation waves, isolated embedded modes are possible, for which the rate of the decay into radiation vanishes. Under more specific conditions, continuous families of embedded solitons were constructed too [34–37]. Very recently, the existence of solitons of this type was demonstrated in a model combining the quadratic (χ(2) ) nonlinearity and a complex

∗ †

Email: [email protected] Email: [email protected]

lattice potential, whose imaginary part is subject to the condition of the PT symmetry, representing symmetrically placed and mutually balanced local gain and loss [38].While embedded solitons have been studied in some detail in basic one-dimensional (1D) models, it is still an open question to find them in Bloch bands of 2D periodic potentials (stable intra-band solitons were reported in a 2D potential which combines a periodic lattice in one direction and a harmonic-oscillator trap acting perpendicular to the lattice [39]). In addition to the studies of ideal lattices, it has been demonstrated that lattice defects may also be used in optical devices, such as integrated circuits [40–42], microcavities [43], and defect-mode lasers [44]. Point and line defects help trapping and guiding light flows in structured photonic media [45]-[53]. In particular, the bandgap guidance of light by defects in 1D and 2D photonic lattices has been realized in Refs. [54–56]. Defectaffected transmission of light was also considered in 2D photonic quasicrystals [57, 58]. The dynamics of solitons under the action of defects in discrete nonlinearSchr¨odinger lattices has been studied too [59]-[64]. The dynamics of BEC trapped in OL potentials may also be strongly affected by defects [65]-[72]. The above-mentioned settings are based on the cubic nonlinearity. Solitons supported by local defects in the form of the localized χ(2) nonlinearity have been studied in Refs. [73–75]. Very recently, solitons pinned to local PT -symmetric defects, inserted into the 1D medium with the uniform χ(2) nonlinearity, were considered too [76]. On the other hand, a specific manifestation of the transmission of light [77]-[84] and matter waves [85]-[88] controlled by random sets of defects is the Anderson localization in the transverse plane, which does not require

2 nonlinearity. More recently, a great deal of attention has been drawn by soliton dynamics in nonlinear lattices, i.e., structures periodically modulating the local nonlinearity [89–91]. In this context, defects can be naturally considered as well. In particular, experimental techniques have been developed for the local control of local nonlinearity in BEC (via the Feshbach resonance [92, 93]), and in photonic crystals [94, 95], which allows one to design desirable nonlinear defects. In this work we report numerical analysis of 2D solitons in lattice potentials with local defects, under the action of the self-focusing nonlinearity. The model is formulated in Section II. The search for such soliton modes is suggested by previously published results demonstrating the existence of linear defect modes in photonic crystals [40]-[53], including 2D modes [96], as well as in ultrasonic crystals [97] and solid-state lattices (see Refs. [98–101] and references therein), although the solitons that we present below are not actually related to linear defect modes. In Section III, we report numerical results which demonstrate the existence of a new type of embedded solitons, alias intra-band solitons, viz., continuous families of 2D localized modes (unlike isolated solutions found previously in typical embedded-soliton models), which are embedded into the first and second Bloch bands of the linear spectrum, in terms of the propagation constant, and pinned to the defect, as concerns the spatial location of the modes. We name them embedded defect solitons (EDSs). It is relevant to stress that, for modes pinned to local lattice defects, it makes sense to identify the location of their propagation constant with respect to the bandgap spectrum of the infinite uniform OL, as the latter determines the spectrum of radiation waves into which the localized mode may decay. As mentioned above, the EDSs do not emerge from a continuation of linear localized defect modes, as they exist with the norm exceeding a finite minimum value (and the propagation constant of linear defect modes cannot be located inside a Bloch band). By means of systematic direct simulations, the EDSs are found to be fully stable in the first band and in a part of the second. Together with the regular defect solitons populating the SIG, and gap defect solitons (GDSs) nested in the first finite bandgap (in the spectral slot between the first and second Bloch bands), which are pinned to the same defect, the intra-band EDSs form a continuous “superfamily” threading the entire band structure. The stability of the GDSs supported by the self-attractive nonlinearity is a remarkable finding, as the previous analysis of gap solitons in defect-free 2D lattices was usually restricted to the case of the self-repulsion, assuming that they would be unstable in the opposite case [19, 102, 103] (gap solitons bifurcating under the action of the self-attractive nonlinearity from edges of the adjacent Bloch bands into the first finite bandgap were constructed in Ref. [104], but their stability was not studied there; in fact, it was recently found that, while 2D fundamental gap solitons

FIG. 1. (Color online) Contour plots of the lattice potential with the single (a), double (b), triple (c), and quadruple (d) defects, at ε = 2.

are completely unstable in the first finite bandgap under the self-attractive nonlinearity, a family of dipole-mode gap solitons are stable in a part of the first bandgap in the same case, provided that the depth of the lattice potential exceeds a certain threshold value [105]). Further, in Section IV we demonstrate the existence of stable EDS and GDS modes with the dipole structure, in the lattice with double and triple local defects, and of vortices supported by quadruple defects. The paper is concluded by Section V. II.

THE MODEL

The analysis is based on the 2D nonlinear Schr¨odinger/Gross-Pitaevskii equation for the mean-field BEC wave function, or the amplitude of the electromagnetic wave propagating in a photonic lattice, ψ(x, y, z), ˜ y)]VOL ψ − |ψ|2 ψ iψz = −(1/2)∇2 ψ − [1 − δ(x,

(1)

(in BEC, propagation distance z is replaced by time t), with the self-attractive nonlinearity and ∇2 = ∂x2 +∂y2 . In the absence of the defect, the lattice potential with depth 2ε and period normalized to be π is taken in the usual form, VOL = ε[cos(2x) + cos(2y)]. The defect is a “hole” of radius r0 = 1.2 in the potential, which is accounted ˜ y) = 1 at x2 + y 2 < r2 , and δ(x, ˜ y) ≡ 0 at for by δ(x, 0 2 2 2 x +y ≥ r0 in Eq. (1). Generic results can be adequately demonstrated for this shape of the defect with ε = 2, which is fixed below (in fact, fully stable soliton families could be found in the interval of 1.1 ≤ r0 ≤ 1.3, as well as for different shapes of the “hole”). The corresponding contour plots of the OL potential with the solitary or compound defects are displayed in Fig. 1. Stationary solutions to Eq. (1) with propagation constant −µ (or chemical potential µ, in terms of BEC)

3

(a)

2

0.5

(b)

2nd band

−2 0

200 400 Mode number (c)

−2 0 2

C

D

−0.5

1st band gap

−1

200 400 Mode number

A

−1.5

(d)

1

2

1st band 3

4

5

N

µ

µ

2

0

µ

0

B

0

µ

µ

2

0 −2 0

0 200 Mode number

−2 400 0

200 400 Mode number

FIG. 2. Linear spectra of the perfect OL (a), and OLs with the single (b), double (c), and triple (d) defects displayed in Fig. 1. Isolated points emerging in the second finite bandgap represent the corresponding linear defect modes.

are sought for as ψ = φ(x, y) exp(−iµz), with function φ(x, y) obeying the stationary equation, 2

2

µφ = −(1/2)∇ φ − [1 − δ(x, y)]VOL (x, y) φ − |φ| φ. (2)

While fundamental and dipole-mode solitons are described by real solutions for φ (x, y), vortices (that can be supported by quadruple defects, see below) naturally correspond to complex ones. Figure 2 displays the linear spectrum of Eq. (1), which has been produced by a numerical solution of the linearized version of Eq. (2). It is seen that, for given values of the parameters, a growing number of defectinduced isolated eigenvalues appear in the second finite bandgap with the increase in the size of the defect (single → double → triple), while the first bandgap remains unaffected. Below, we focus not on quite obvious quasilinear defect modes corresponding to these isolated eigenvalues, but rather on localized states which have no linear limit, but may be created by the defects in the first finite bandgap, as well as in the two Bloch bands adjacent to it. These states were constructed as numerical solutions of Eq. (2) by means of the Newton’s method in the domain of size 30 × 30, covered by a grid of 192 × 192 points. The initial Gaussian,  guess was
the isotropic
 φ (x, y) = A exp − x2 + y 2 / 2W 2 , with amplitude A, width W , and power (or number of atoms, in R Rtotal 2 the BEC), N ≡ φ2 (x, y) dxdy = π (AW ) . The stability of the solutions was subsequently tested in direct simulations of perturbed evolution in the framework of Eq. (1), using the same numerical domain with absorbing boundary conditions.

FIG. 3. The propagation constant (−µ) versus the total power norm, for the families of solitons pinned to the single defect [it is shown in Fig. 1(a)], which were found inside the first and second Bloch bands, as well as in the first finite bandgap between them. The solid and dashed curves indicate stable and unstable solutions, respectively. Typical examples of the pinned solitons, corresponding to the marked points, are shown below in Figs. 4 and 5. Some irregularity of the curves is explained by the fact that they were plotted with interval ∆µ = 0.1.

III.

RESULTS OF THE ANALYSIS FOR THE SINGLE DEFECT

As shown in Fig. 3, the numerical solution reveals the existence of continuous families of localized modes (fundamental EDSs), pinned to the single defect [the one shown in Fig. 1(a)], inside the first and second Bloch bands. The EDS families are continuously linked to the family of GDSs, attached to the same defect, whose propagation constant belongs to the first finite bandgap. The EDS family in the first Bloch band starts exactly at its boundary with the SIG, where it is linked to the family of stable regular solitons pinned to the defect. The solitons of the latter type are not presented here in detail, as their existence and properties are quite obvious. It is observed that the combined family of the pinned solitons in Fig. 3 exists above a finite minimum value of the total power, Nmin = 1.56, which exactly corresponds to the boundary between the first Bloch band and the first finite bandgap. This means, as said above, that the family of the localized modes supported by the defect does not emerge as a nonlinear continuation of any linear defect mode, as the latter would correspond to N → 0. Direct simulations demonstrates that the EDS family is completely stable in the first Bloch band, while the EDS becomes unstable in the second band beyond a certain critical point. It is worthy to note that the stability of the intra-band solitons supported by the defect in the first Bloch band agrees with the Vakhitov-Kolokolov (VK) criterion, which is relevant to localized modes in the case of the self-focusing nonlinearity [106], ∂µ/∂N < 0, while the EDS branch features ∂µ/∂N > 0 in the second band, where the modes are stable in a limited region (strictly speaking, this may be a region of a very weak instability, which is, nevertheless, tantamount to stabil-

4

FIG. 5. (Color online) An example of a stable gap soliton pinned to the single defect, found in the first bandgap (at point D marked in Fig. 3) under the self-focusing nonlinearity.

solitons in the first finite bandgap (in the present case, they all are pinned to the single defect).

FIG. 4. (Color online) Examples of the intra-band (embedded ) defect solitons, EDSs, found in the first (a) and second (b,c) Bloch bands, which correspond to points A and B,C, respectively, in Fig. 3. The shapes of the modes are shown by means of 3D images (left) and contour plots (right). The mode in panels (a) and (b) are stable, while the one in (c) is unstable.

ity in terms of a possible experiment). Of course, there is no direct proof of the applicability of the VK criterion to the present model. Typical examples of stable EDSs in the first and second Bloch bands, found at points marked (A) and (B) in Fig. 3, along with an example of an unstable EDS corresponding to point C, are displayed in Fig. 4. Throughout the first band, the shape of the solitons is very similar to that displayed in Fig. 4(a). On the other hand, the localized modes feature a more complex structure, with pronounced side peaks, in the second band [see Figs. 4(b,c)], in accordance with the fact that the shape of Bloch modes is more complex too in the same band. In fact, the extended tail of the soliton shown in Fig. 4(c) initiates the onset of the instability of this soliton. In direct simulations, unstable modes decay into radiation (not shown here in detail). To the best of our knowledge, the current model presents the first example of nonlinear localized modes found inside Bloch bands of the spectrum induced by the 2D lattice (in Ref. [39], which was dealing with 2D solitons of the embedded type, the lattice was onedimensional). In fact, the EDS family threading the first Bloch band plays the role of a missing link between the two well-known types of solitons existing in the lattice potential, viz., regular ones populating the SIG, and gap

As concerns the family of GDSs found inside the first finite bandgap, which links the intra-band EDS families populating the two adjacent Bloch bands (see Fig. 3), our numerical tests have demonstrated that these solitons, being pinned to the defect, are stable only in the case of the self-focusing nonlinearity [as set in Eq. (1)], on the contrary to gap solitons in perfect lattices, which are usually assumed to be stable solely in the case of the self-defocusing [19, 102]. The case of the self-defocusing nonlinearity was explored too, and it was concluded that the pinned GDSs are unstable in that case (not displayed here in detail). The latter finding may be explained by the fact that the effective mass of the gap soliton supported by the self-defocusing nonlinearity is negative, thus reversing the character of the soliton-defect interaction, and making unstable the bound state of the gap soliton pinned to the attractive defect. In line with this argument, additional analysis demonstrates that the local defect of the opposite sign, in comparison with the one considered here, gives rise to stable pinned states of gap solitons under the self-defocusing nonlinearity, and unstable states in the case of the self-focusing. The negative mass of gap solitons explains a number of other counterintuitive dynamical effects featured by these modes [107][110]. Numerical results also demonstrate that the selfdefocusing nonlinearity does not give rise to solitons that would be embedded into Bloch bands and pinned to the defect. A typical example of the stable pinned GDS, found in the first finite bandgap in the case of the self-focusing nonlinearity, which corresponds to point D in Fig. 3, is presented in Fig. 5. Finally, our numerical simulations suggest that stable GDSs pinned to the defect do not exist in the relatively narrow second finite bandgap (see Fig. 2 ), in agreement with the general trend of gap solitons to be unstable in the second bandgap [16, 17, 21].

5

FIG. 6. (Color online) Examples of stable dipole-mode solitons supported by the double defect, shown in Fig. 1(b), in the first Bloch band (left, with µ = −1.4, N = 3.9), and in the first finite bandgap (right, with µ = −0.3, N = 5.4). FIG. 8. (Color online) Examples of stable vortices with topological charge 1, supported by the quadruple defect inside the first Bloch band (top, µ = −1.5 and N = 8.8), and in the first bandgap (bottom, µ = −0.8 and N = 7.2 ). The contour plots show the phase distribution in the vortices. 0.5 2nd band

µ

0 −0.5

1st band gap

−1

FIG. 7. (Color online) Examples of stable dipole-mode solitons supported by the triple triangular defect, shown in Fig. 1(c), in the first Bloch band (left, µ = −1.45, N = 6.3), and in the first finite bandgap (right, µ = −0.5, N = 6.9).

IV.

RESULTS FOR MULTIPLE DEFECTS: STABLE DIPOLES AND VORTICES

We have extended the above analysis for the double, triple, and quadruple lattice defects, which are shown in Fig. 1(b,c,d). While the solitary defect may only support fundamental solitons of both the EDS and GDS types, as shown above, the double defects readily give rise to stable dipole-mode bound states of two solitons with opposite signs. Examples of such modes, found in the first Bloch band and in the first finite bandgap (i.e., of the EDS and GDS types, respectively), are shown in Fig. 6. The triangular defects [see Fig. 1(c)] also support stable triangularly shaped dipole modes, see examples displayed in Fig. 7. Further, the quadruple defects [see Fig. 1(d)] create stable solitary vortices, characteristic examples of which are shown in Fig. 8. As is typical for vortices supported by lattice potentials [19, 91, 111, 112], they are built as rhombic sets of four intensity peaks, with a nearly empty site in the center (therefore this type of vortices is often called onsite-centered), and phase shifts π/2 between the peaks, which corresponds to the global phase circulation of 2π (i.e., topological charge 1). As indicated in Fig. 9, direct simulations demonstrate that the

1st band

−1.5 6

8

10

12

14

N FIG. 9. Curve µ(N ) for the vortices with topological charge 1, pinned to the quadrupole lattice defect. The solid and dashed curves show stable and unstable solutions, respectively.

pinned vortices are unstable in the second Bloch band, being stable elsewhere.

V.

CONCLUSION

In this work, we have studied localized states which can be supported by single, double, triple, and quadruple defects of 2D lattice potentials in the self-focusing medium. The model describes lattices in nonlinear photonic media and BEC trapped in the OL. The first result is that the single defect supports families of intra-band (embedded ) solitons, in the first and second Bloch bands. These modes exist above the minimum value of the norm, i.e., they do not emerge from any linear defect mode. Direct simulations demonstrate stability of the family in the first Bloch band, and in a part of the second band. The branch of the stable pinned solitons embedded into the first band links families of regular and gap solitons

6 hosted by the semi-infinite and first finite bandgaps, respectively, and pinned to the same defect. The family of the gap solitons, pinned to the defect, is stable, despite the fact that the nonlinearity is self-attractive, while it is usually assumed that 2D fundamental gap solitons may be stable only under the self-repulsion (in uniform lattices). Finally, it has been found that the double and triple defects support stable dipole solitons, in the first Bloch band and in the first finite bandgap alike, and, similarly, the quadruple defect creates stable vortex solitons. This work can be naturally extended in other directions—in particular, to the model with a rectilinear defect, as an elongated version of the double one. Search for stable vortex solitons with higher topological charges, and, possibly, “supervortices” (vortex rings built of compact localized vortices) [113, 114] may be relevant too. On the other hand, it may be interesting to consider defects in complex lattices, whose imaginary part is subject

to the condition of the PT symmetry. 2D gap solitons in uniform PT -symmetric lattices were found recently, but they are unstable [115]. Stable 1D solitons pinned to defects in PT -symmetric lattices, embedded into the medium with the self-defocusing cubic nonlinearity, were reported very recently in Ref. [116]. In this connection, it is also relevant to mention that 1D and 2D gap solitons can be readily made stable in a more general GinzburgLandau system, which does not impose the balance condition on the gain and dissipation [117, 118].

[1] J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the Flow of Light (Princeton University Press: Princeton, 2008). [2] J. A. Porto, F. J. Garcia-Vidal, and J. B. Pendry, “Transmission Resonances on Metallic Gratings with Very Narrow Slits,” Phys. Rev. Lett. 83, 2845-2848 (1998). [3] F. J. Garcia-Vidal and L. Martin-Moreno, “Transmission and focusing of light in one-dimensional periodically nanostructured metals,” Phys. Rev. B 66, 155412 (2002). [4] S. A. Darmanyan and A. V. Zayats, “Light tunneling via resonant surface plasmon polariton states and the enhanced transmission of periodically nanostructured metal films: An analytical study,” Phys. Rev. B 67, 035424 (2003). [5] K. J. K. Koerkamp, S. Enoch, F. B. Segerink, N. F. van Hulst, and L. Kuipers, “Strong Influence of Hole Shape on Extraordinary Transmission through Periodic Arrays of Subwavelength Holes,” Phys. Rev. Lett. 92, 183901 (2004). [6] J. G. Rivas, M. Kuttge, P. H. Bolivar, H. Kurz, and J. A. Sanchez-Gil, “Propagation of Surface Plasmon Polaritons on Semiconductor Gratings,” Phys. Rev. Lett. 93, 256804 (2004). [7] A. Boltasseva, S. I. Bozhevolnyi, T. Nikolajsen, and K. Leosson, “Compact Bragg Gratings for Long-Range Surface Plasmon Polaritons,” J. Lightwave Tech. 24, 912-918 (2006). [8] Z. Han, E. Forsberg, and S. He, “Surface Plasmon Bragg Gratings Formed in Metal-Insulator-Metal Waveguides,” IEEE Phot. Tech. Lett. 19, 91-93 (2007). [9] J. Park, H. Kim, and B. Lee, “High order plasmonic Bragg reflection in the metal-insulator-metal waveguide Bragg grating,” Opt. Exp. 16, 413-425 (2008). [10] P. Genevet, J. P. Tetienne, E. Gatzogiannis, R. Blanchard, M. A. Kats, M. O. Scully, F. Capasso, “Large Enhancement of Nonlinear Optical Phenomena by Plasmonic Nanocavity Gratings,” Nano Lett. 10, 4880-4883

(2010). [11] G. D’Aguanno, N. Mattiucci, M. J. Bloemer, D. de Ceglia, M. A. Vincenti, and A. Al´ u, “Transmission resonances in plasmonic metallic gratings,” J. Opt. Soc. Am. B 28, 253-264 (2011). [12] D. de Ceglia, M. A. Vincenti, M. Scalora, N. Akozbek, and M. J. Bloemer, “Plasmonic band edge effects on the transmission properties of metal gratings,” AIP Advances 1, 032151 (2011). [13] I. Dolev, I. Epstein, and A. Arie, “Surface-Plasmon Holographic Beam Shaping ,” Phys. Rev. Lett. 109, 203903 (2012). [14] L. Dal Negro and S. V. Boriskina, “Deterministic aperiodic nanostructures for photonics and plasmonics applications,” Laser Phot. Rev. 6, 178-218 (2012). [15] C. J. Pethick and H. Smith, Bose-Einstein condensate in dilute gas (Cambridge University Press: Cambridge, 2008). [16] V. A. Brazhnyi and V. V. Konotop, “Theory of nonlinear matter waves in optical lattices,”Mod. Phys. Lett. B 18, 627-651 (2004). [17] O. Morsch and M. Oberthaler, “Dynamics of BoseEinstein condensates in optical lattices,” Rev. Mod. Phys. 78, 179-215 (2006). [18] H. T. C. Stoof, K. B. Gubbels, and D. B. M. Dickerscheid, Ultracold Quantum Fields (Springer: Dordrecht, 2009). [19] B. A. Malomed, D. Mihalache, F. Wise, and L. Torner, “Spatiotemporal optical solitons,” J. Optics B: Quant. Semicl. Opt. 7, R53-R72 (2005). [20] Y. V. Kartashov, V. A. Vysloukh, and L. Torner, “Soliton Shape and Mobility Control in Optical Lattices,” Progress in Optics 52, 63-148 (ed. by E. Wolf: North Holland,Amsterdam, 2009). [21] D. E. Pelinovsky, Localization in Periodic Potential: From Schr¨ odinger Operators to the Gross-Pitaevskii Equation (Cambridge University Press: Cambridge, 2011).

VI.

ACKNOWLEDGMENTS

J.Z. acknowledges support from the Natural Science Foundation of China (Project No. 11204151) and China Postdoctoral Science Foundation (2012M520238).

7 [22] T. Anker, M. Albiez, R. Gati, S. Hunsmann, B. Eiermann, A. Trombettoni, and M. K. Oberthaler, “Nonlinear Self-Trapping of Matter Waves in Periodic Potentials,” Phys. Rev. Lett. 94, 020403 (2005). [23] T. J. Alexander, E. A. Ostrovskaya, and Y. S. Kivshar, “Self-Trapped Nonlinear Matter Waves in Periodic Potentials,” Phys. Rev. Lett. 96, 040401 (2006). [24] F. H. Bennet, T. J. Alexander, F. Haslinger, A. Mitchell, D. N. Neshev, and Y. S. Kivshar, “Observation of Nonlinear Self-Trapping of Broad Beams in Defocusing Waveguide Arrays ,” Phys. Rev. Lett. 106, 093901 (2011). [25] J. Fujioka and A. Espinosa, “Soliton-Like Solution of an Extended NLS Equation Existing in Resonance with Linear Dispersive Waves,”J. Phys. Soc. Jpn. 66, 26012607 (1997). [26] J. Yang, B. A. Malomed, D. J. Kaup, “Embedded Solitons in Second-Harmonic-Generating Systems,”Phys. Rev. Lett. 83, 1958-1961 (1999). [27] A. R. Champneys and B. A. Malomed, “ Moving embedded solitons,”J. Phys. A: Math. Gen. 32, L547-L553 (1999). [28] A. R. Champneys, B. A. Malomed, J. Yang, and D. J. Kaup, “ Embedded solitons: solitary waves in resonance with the linear spectrum,” Physica D 152, 340354 (2001). [29] J. Atai and B. A. Malomed, “Solitary waves in systems with separated Bragg grating and nonlinearity,”Phys. Rev. E 64, 066617 (2001). [30] D. E. Pelinovsky and J. Yang, “ A normal form for nonlinear resonance of embedded solitons,”Proc. Roy. Soc. L. A 458, 1469-1497 (2002). [31] A. Espinosa-Cer´ on, J. Fujioka, and A. G´ omezRodr´ıguez, “ Embedded Solitons: Four-Frequency Radiation, Front Propagation and Radiation Inhibition,” Phys. Scripta 67, 314-324 (2003). [32] R. F. Rodr´ıguez, J. A. Reyes, A. Espinosa-Cer´ on, J. Fujioka, and B. A. Malomed, “ Standard and embedded solitons in nematic optical fibers,”Phys. Rev. E 68, 036606 (2003). [33] A. Gubeskys and B. A. Malomed, “Spontaneous soliton symmetry breaking in two-dimensional coupled Bose-Einstein condensates supported by optical lattices,”Phys. Rev. A 76, 043623 (2007). [34] J. Yang, “Dynamics of Embedded Solitons in the Extended Korteweg - de Vries Equations,” Stud. Appl. Math. 106, 337-365 (2001). [35] J. Yang and T. R. Akylas, “Continuous Families of Embedded Solitons in the Third-Order Nonlinear Schr¨ odinger Equation,” Stud. Appl. Math. 111, 359-375 (2003). [36] J. Yang, “Stable Embedded Solitons,” Phys. Rev. Lett. 91, 143903 (2003). [37] W. C. K. Mak, B. A. Malomed, and P. L. Chu, “Symmetric and asymmetric solitons in linearly coupled Bragg gratings,” Phys. Rev. E 69, 066610 (2004). [38] F. C. Moreira, V. V. Konotop, and B. A. Malomed, “Solitons in PT -symmetric periodic systems with the quadratic nonlinearity,” Phys. Rev. A 87, 013832 (2013). [39] J. Yang, “Fully localized two-dimensional embedded solitons,” Phys. Rev. A 82, 053828 (2010). [40] S. Noda, A. Chutinan, and M. Imada,“Trapping and emission of photons by a single defect in a photonic

[41]

[42]

[43]

[44]

[45]

[46]

[47]

[48]

[49]

[50]

[51]

[52]

[53] [54] [55] [56]

bandgap structure,” Nature (London) 407, 608-610 (2000). M. Qi, E. Lidorikis, P. T. Rakich, S. G. Johnson, J. D. Joannopoulos, E. P. Ippen, and H. I. Smith, “A threedimensional optical photonic crystal with designed point defects,” Nature (London) 429,538-542 (2004). S. A. Rinne, F. Garc´ıa-Santamar´ıa, and P. V. Braun, “Embedded cavities and waveguides in threedimensional silicon photonic crystals,” Nature Photon. 2, 52-56 (2008). P. R. Villeneuve, S. Fan, and J. D. Joannopoulos, “Microcavities in photonic crystals: Mode symmetry, tunability, and coupling efficiency,” Phys. Rev. B 54, 78377842 (1996). O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-Dimensional Photonic Band-Gap Defect Mode Laser,” Science 284, 1819-1821 (1999). E. Yablonovitch, T. J. Gmitter, R. D. Meade, A. M. Rappe, K. D. Brommer, and J. D. Joannopoulos, “Donor and acceptor modes in photonic band structure,” Phys. Rev. Lett. 67, 3380-3383 (1991). M. Bayindir, B. Temelkuran, and E. Ozbay, “TightBinding Description of the Coupled Defect Modes in Three-Dimensional Photonic Crystals,” Phys. Rev. Lett. 84, 2140-2143 (2000). V. Yannopapas, N. Stefanou, and A. Modinos, “Effect of Stacking Faults on the Optical Properties of Inverted Opals,” Phys. Rev. Lett. 86, 4811-4814 (2001). A. A. Sukhorukov and Y. S. Kivshar, “Nonlinear Localized Waves in a Periodic Medium ,” Phys. Rev. Lett. 87, 083901 (2001). M. Notomi, K. Yamada, A. Shinya, J. Takahashi, C. Takahashi, and I. Yokohama, “Extremely Large GroupVelocity Dispersion of Line-Defect Waveguides in Photonic Crystal Slabs,” Phys. Rev. Lett. 87, 253902 (2001). D. Englund, D. Fattal, E. Waks, G. Solomon, B. Zhang, T. Nakaoka, Y. Arakawa, Y. Yamamoto, and J. Vuˇckovi´c, “Controlling the Spontaneous Emission Rate of Single Quantum Dots in a Two-Dimensional Photonic Crystal,” Phys. Rev. Lett. 95, 013904 (2005). E. Gavartin, R. Braive, I. Sagnes, O. Arcizet, A. Beveratos, T. J. Kippenberg, and I. Robert-Philip, “Optomechanical Coupling in a Two-Dimensional Photonic Crystal Defect Cavity ,” Phys. Rev. Lett. 106, 203902 (2011). S. Lin, E. Chow, V. Hietala, P. R. Villeneuve, and J. D. Joannopoulos, “Experimental Demonstration of Guiding and Bending of Electromagnetic Waves in a Photonic Crystal,” Science, 282, 274-276 (1998). S. Ogawa, M. Imada, S. Yoshimoto, M. Okano, S. Noda, “Control of Light Emission by 3D Photonic Crystals,” Science 305, 227-229 (2004). J. Yang and Z. Chen, “Defect solitons in photonic lattices,” Phys. Rev. E 73, 026609 (2006). I. Makasyuk, Z. Chen, and J. Yang, “Band-Gap Guidance in Optically Induced Photonic Lattices with a Negative Defect,” Phys. Rev. Lett. 96, 223903 (2006). M. Rechtsman, A. Szameit, F. Dreisow, M. Heinrich, R. Keil, S. Nolte, and M. Segev, “Amorphous Photonic Lattices: Band Gaps, Effective Mass, and Suppressed Transport,” Phys. Rev. Lett. 106, 193904 (2011).

8 [57] Y. S. Chan, C. T. Chan, and Z. Y. Liu, “Photonic Band Gaps in Two Dimensional Photonic Quasicrystals,” Phys. Rev. Lett. 80, 956-959 (1998). [58] B. Freedman, G. Bartal, M. Segev, R. Lifshitz, D. N. Christodoulides, and J. W. Fleischer, “Wave and defect dynamics in nonlinear photonic quasicrystals,” Nature (London) 440, 1166-1169 (2006). [59] M. I. Molina and G. Tsironis, “Nonlinear impurities in a linear chain”, Phys. Rev. B 47, 15330-15333 (1993). [60] B. C. Gupta and K. Kundu, “Formation of stationary localized states due to nonlinear impurities using the discrete nonlinear Schr¨ odinger equation Phys. Rev. B 55, 894-905 (1997). [61] B. C. Gupta and K. Kundu, “Stationary localized states due to a nonlinear dimeric impurity embedded in a perfect one-dimensional chain”, Phys. Rev. B 55, 1103311036 (1997). [62] E. Bulgakov, K. Pichugin, and A. Sadreev, “Symmetry breaking for transmission in a photonic waveguide coupled with two off-channel nonlinear defects”, Phys. Rev. B 83, 045109 (2011). [63] V. A. Brazhnyi and B. A. Malomed, “ Localization and delocalization of two-dimensional discrete solitons pinned to linear and nonlinear defects”, Phys. Rev. E 83, 016604 (2011). [64] Phys. Rev. E 84, “Dragging two-dimensional discrete solitons by moving linear defects”,016608 (2011). [65] T. Mayteevarunyoo, B. A. Malomed, and G. Dong, “Spontaneous symmetry breaking in a nonlinear doublewell structure,” Phys. Rev. A 78, 053601 (2008). [66] N. Dror and B. A. Malomed, “Solitons supported by localized nonlinearities in periodic media,” Phys. Rev. A 83, 033828 (2011). [67] T. Mayteevarunyoo, B. A. Malomed, and A. Reoksabutr, “Spontaneous symmetry breaking of photonic and matter waves in two-dimensional pseudopotentials,” J. Mod. Opt. 58, 1977-1989 (2011) [68] V. A. Brazhnyi, V. V. Konotop, and V. M. P´erezGarc´ıa, “Defect modes of a Bose-Einstein condensate in an optical lattice with a localized impurity,” Phys. Rev. A 74, 023614 (2006). [69] V. A. Brazhnyi, V. V. Konotop, and V. M. P´erezGarc´ıa, “Driving Defect Modes of Bose-Einstein Condensates in Optical Lattices,” Phys. Rev. Lett. 96, 060403 (2006). [70] Y. V. Kartashov, V. V. Konotop, V. A. Vysloukh, and L. Torner, “Dissipative defect modes in periodic structures,” Opt. Lett. 35, 1638-1640 (2010). [71] S. Longhi, “Bloch oscillations in tight-binding lattices with defects,” Phys. Rev. B 81, 195118 (2010). [72] S. Longhi, “Dynamic reflectionless defects in tightbinding lattices,” Phys. Rev. B 84. 193105 (2011). [73] A. A. Sukhorukov and Y. S. Kivshar, “Spatial optical solitons in nonlinear photonic crystals,” Phys. Rev. E 65, 036609 (2002). [74] A. Shapira, N. Voloch-Bloch, B. A. Malomed, and A. Arie, “Spatial quadratic solitons guided by narrow layers of a nonlinear material,” J. Opt. Soc. Am. B 28, 1481-1489 (2011). [75] V. A. Brazhnyi and B. A. Malomed, “Symmetric and asymmetric localized modes in linear lattices with an embedded pair of χ(2) -nonlinear sites,” Phys. Rev. A, 86,013829 (2012).

[76] F. C. Moreira, F. Kh. Abdullaev, V. V. Konotop, and A. V. Yulin, “Localized modes in χ(2) media with PT -symmetric localized potential”, Phys. Rev. A 86, 053815 (2012). [77] D. S. Wiersma, P. Bartolini, A. Lagendijk, and R. Righini, “Localization of light in a disordered medium,” Nature 390, 671-673 (1997)¡£ [78] M. C. W. van Rossum and T. M. Nieuwenhuizen, “Multiple scattering of classical waves: microscopy, mesoscopy, and diffusion,” Rev. Mod. Phys. 71, 313371 (1999)¡£ [79] Y. A. Vlasov, M. A. Kaliteevski, and V. V. Nikolaev, “Different regimes of light localization in a disordered photonic crystal,” Phys. Rev. B 60, 1555-1562 (1999)¡£ [80] H. Cao, J. Y. Xu, E. W. Seelig, R. P. H. Chang, “Microlaser made of disordered media,” Appl. Phys. Lett. 76, 2997-2999 (2000)¡£ [81] C. Vanneste and P. Sebbah, “Selective Excitation of Localized Modes in Active Random Media,” Phys. Rev. Lett. 87, 183903 (2001). [82] M. Storzer, P. Gross, C. M. Aegerter, and G. Maret, “Observation of the Critical Regime Near Anderson Localization of Light,” Phys. Rev. Lett. 96, 063904 (2006). [83] T. Schwartz, G. Bartal, S. Fishman, and M. Segev, “Transport and Anderson localization in disordered two-dimensional photonic lattices,” Nature 446, 52-55 (2007). [84] Y. Lahini, A. Avidan, F. Pozzi, M. Sorel, R. Morandotti, D. N. Christodoulides, Y. Silberberg, “Anderson Localization and Nonlinearity in One-Dimensional Disordered Photonic Lattices,” Phys. Rev. Lett. 100, 013906 (2008). [85] U. Gavish and Y. Castin, “Matter-Wave Localization in Disordered Cold Atom Lattices,” Phys. Rev. Lett. 95, 020401 (2005). [86] J. Billy, V. Josse, Z. C. Zuo, A. Bernard, B. Hambrecht, P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer, and A. Aspect, “Direct observation of Anderson localization of matter waves in a controlled disorder,” Nature 453, 891-894 (2008). [87] G. Roati, C. D’Errico, L. Fallani, M. Fattori, C. Fort, M. Zaccanti, G. Modugno, M. Modugno, M. Inguscio, “Anderson localization of a non-interacting Bose-Einstein condensate,” Nature 453, 895-898 (2008). [88] L. Sanchez-Palencia and M. Lewenstein, “Disordered quantum gases under control,” Nature Physics 6, 8795 (2010). [89] Y. V. Kartashov, B. A. Malomed, and L. Torner, “Solitons in nonlinear lattices,” Rev. Mod. Phys. 83, 247-305 (2011). [90] J. Zeng and B. A. Malomed, “Stabilization of onedimensional solitons against the critical collapse by quintic nonlinear lattices,” Phys. Rev. A 85,023824 (2012). [91] J. Zeng and B. A. Malomed, “Two-dimensional solitons and vortices in media with incommensurate linear and nonlinear lattice potentials,” Phys. Scripta T149, 014035 (2012). [92] S. Ghanbari, T. D. Kieu, A. Sidorov, and P. Hannaford, “Permanent magnetic lattices for ultracold atoms and quantum degenerate gases,” J. Phys. B: At. Moil. Opt. Phys. 39, 847-860 (2006).

9 [93] A. Abdelrahman, P. Hannaford, and K. Alameh, “Adiabatically induced coherent Josephson oscillations of ultracold atoms in an asymmetric two-dimensional magnetic lattice,” Opt. Express 17, 24358-24370 (2009). [94] V. Berger, “Nonlinear Photonic Crystals,” Phys. Rev. Lett. 81, 4136-4139 (1998). [95] J. Hukriede, D. Runde, and D. Kip, “Fabrication and application of holographic Bragg gratings in lithium niobate channel waveguides,” J. Phys. D 36, R1-R16 (2003). [96] O. Painter, J. Vuˇckovi´c, and A. Scherer, “Defect modes of a two-dimensional photonic crystal in an optically thin dielectric slab,” J. Opt. Soc. Am. B 16, 275-285 (1999). [97] A. Khelif, A. Choujaa, B. Djafari-Rouhani, M. Wilm, S. Ballandras, and V. Laude, “Trapping and guiding of acoustic waves by defect modes in a full-band-gap ultrasonic crystal,” Phys. Rev. B 68, 214301 (2003). [98] S. R. Bickham and A. J. Sievers, “Intrinsic localized modes in a monatomic lattice with weakly anharmonic nearest-neighbor force constants,” Phys. Rev. B 43, 2339-2346 (1991). [99] H. R. Schober and B. B. Laird, “Localized low-frequency vibrational modes in glasses,” Phys. Rev. B 44, 67466754 (1991). [100] G. Kaczmarczyk, A. Kaschner, A. Hoffmann, and C. Thomsen, “Impurity-induced modes of Mg, As, Si, and C in hexagonal and cubic GaN,” Phys. Rev. B 61, 53535357 (2000). [101] G. Van de Walle and J. Neugebauer,“First-principles calculations for defects and impurities: Applications to III-nitrides,” J. Appl. Phys. 95, 3851-3879 (2004). [102] E. A. Ostrovskaya and Y. S. Kivshar, “Matter-Wave Gap Solitons in Atomic Band-Gap Structures,” Phys. Rev. Lett. 90, 160407 (2003). [103] E. A. Ostrovskaya and Y. S. Kivshar, “Photonic crystals for matter waves: Bose-Einstein condensates in optical lattices,” Opt. Express 12, 19-29 (2004). [104] Z. Shi and J. Yang, “Solitary waves bifurcated from Bloch-band edges in two-dimensional periodic media”, Phys. Rev. E 75, 056602 (2007). [105] N. Dror and B. A. Malomed, “Stability of dipole gap solitons in two-dimensional lattice potentials”, in: Spectral Analysis, Stability and Bifurcations in Nonlinear Physical Systems, edited by O. Kirillov and D. E. Peli-

novsky (Springer, in press). [106] M. Vakhitov and A. Kolokolov, “Stationary solutions of the wave equation in a medium with nonlinearity saturation,” Radiophys. Quantum. Electron. 16, 783-789 (1973). [107] H. Sakaguchi and B. A. Malomed, “Dynamics of positive- and negative-mass solitons in optical lattices and inverted traps,” J. Phys. B 37,1443-1459 (2004). [108] H. Sakaguchi and B. A. Malomed, “Two-dimensional loosely and tightly bound solitons in optical lattices and inverted traps,” J. Phys. B 37, 2225-2239 (2004). [109] M. Matuszewski, B. A. Malomed, and M. Trippenbach, “Competition between attractive and repulsive interactions in two-component Bose-Einstein condensates trapped in an optical lattice,” Phys. Rev. A 76, 043826 (2007). [110] S. K. Adhikari and B. A. Malomed, “Symbiotic gap and semigap solitons in Bose-Einstein condensates,” Phys. Rev. A 77, 023607 (2008). [111] B. B. Baizakov, B. A. Malomed, and M. Salerno, “Multidimensional solitons in periodic potentials,” Europhys. Lett. 63, 642-648 (2003). [112] J. Yang and Z. H. Musslimani, “Fundamental and vortex solitons in a two-dimensional optical lattice,” Opt. Lett. 28, 2094-2096 (2003). [113] H. Sakaguchi and B. A. Malomed, “Higher-order vortex solitons, multipoles, and supervortices on a square optical lattice,” Europhys. Lett. 72, 698-704 (2005). [114] R. Driben and B. A. Malomed, “Stabilization of twodimensional solitons and vortices against supercritical collapse by lattice potentials,” Eur. Phys. J. D 50, 317323 (2008). [115] J. Zeng and Y. Lan, “Two-dimensional solitons in PT linear lattice potentials,” Phys. Rev. E 85, 047601 (2012). [116] S. Hu and W. Hu, “Defect solitons in parity-time symmetric optical lattices with self-defocusing nonlinearity”, Opt. Commun. 294, 311-324 (2013). [117] H. Sakaguchi and B. A. Malomed, “Gap solitons in Ginzburg-Landau media,” Phys. Rev. E 77, 056606 (2008). [118] H. Sakaguchi and B. A. Malomed, “Two-dimensional dissipative gap solitons,” Phys. Rev. E 80, 026606 (2009).