Two-component gap solitons in two- and one-dimensional Bose-Einstein condensates

We introduce two- and one-dimensional (1D) models of a binary BEC (Bose-Einstein condensate) in a periodic potential, with repulsive interactions. We chiefly consider the most fundamental case of the interspecies repulsion with zero intraspecies interactions. The same system may also model a mixture of two mutually repulsive fermionic species. Existence and stability regions for gap solitons (GSs) supported by the interplay of the interspecies repulsion and periodic potential are identified. Two-component GSs are constructed by means of the variational approximation (VA) and in a numerical form. The VA provides an accurate description for the GS which is a compound state of two tightly bound components, each essentially trapped in one cell of the periodic potential. GSs of this type dominate in the case of intragap solitons, with both components belonging to the first finite bandgap of the linear spectrum (only this type of solitons is possible in a weak lattice). Intergap solitons, with one component residing in the second bandgap, and intragap solitons which have both components in the second gap, are possible in a deeper periodic potential, with the strength essentially exceeding the recoil energy of the atoms. Intergap solitons are, typically, bound states of one tightly and one loosely bound component. In this case, results are obtained in a numerical form. The number of atoms in experimentally relevant situations is estimated to be similar to 5000 in the 2D intragap soliton, and similar to 25 000 in its intergap counterpart; in 1D solitons, it may be up to 10(5). For 2D solitons, the stability is identified in direct simulations, while in the 1D case it is done via eigenfrequencies of small perturbations, and then verified by simulations. In the latter case, if the intragap soliton in the first bandgap is weakly unstable, it evolves into a stable breather, while unstable solitons of other types (in particular, intergap solitons) get completely destroyed. The intragap 2D solitons in the first bandgap are less robust, and in some cases they are completely destroyed by the instability. Addition of intraspecies repulsion to the repulsion between the components leads to further stabilization of the GSs.