“…Our numerical computations show that, in absence of coupling between the modes, the bright soliton can be identified in the first gap for µ 1d + ν < µ b < µ 1d (i.e., µ b ∈ [0.6192,0.6755], and ν here as well as below denotes an appropriate shift), whereas bubble-type solutions also exist for µ 2u + ν < µ d < µ 2u (i.e., µ d ∈ [0.4618,0.5181]). Due to the nonlinear shift of the excited second band toward lower values of µ, the propagation constant of the bubble falls into the range of first gap of the linear band structure [43]. In turn, the dark soliton can be identified for lower values of the propagation constant, namely, for µ 3u + ν < µ d < µ 3u (i.e., µ d ∈ [0.3004,0.3567]).…”