By means of the variational approximation (VA) and systematic simulations, we study dynamics and stability boundaries for solitons in a two-dimensional (2D) self-attracting Bose-Einstein condensate (BEC), trapped in an optical lattice (OL) whose amplitude is subjected to periodic time modulation (the modulation frequency $\ensuremath{\omega}$ may be in the range of several kHz). Regions of stability of the solitons against the collapse and decay are identified in the space of the model's parameters. A noteworthy result is that the stability limit may reach the largest (100%) modulation depth, and the collapse threshold may exceed its classical value in the static lattice (which corresponds to the norm of Townes soliton). The minimum norm
Nmin necessary for the stability of the solitons is identified too. It features a strong dependence on $\ensuremath{\omega}$ at low frequencies, due to a resonant decay of the soliton. Predictions of the VA are reasonably close to results of the simulations. In particular, the VA helps understand salient resonant features in the shape of the stability boundaries observed with the variation of $\ensuremath{\omega}$.
In a basic physical model where two-dimensional (2D) matter-wave solitons may be stable, namely, the Gross-Pitaevskii equation with the self-attractive nonlinearity and quasi-one-dimensional (1D) optical-lattice (OL) potential, we test robustness of the solitons against periodic time modulation of the OL strength. Stability diagrams for the 2D solitons are presented in the plane of the modulation depth and frequency. Basic features of the diagrams are explained with the help of the variational approximation for the stationary counterpart of the model. In the Bose-Einstein condensate of
7Li atoms, the stable 2D solitons may contain the number of atoms in the range of
104−−105, relevant values of the OL strength and modulation frequency being, respectively, $\ensuremath{\lesssim}5$ recoil energies and $\ensuremath{\lesssim}10\text{ }\text{kHZ}$. Head-on collisions between stable 2D solitons moving in the unconfined direction are studied in detail too, for velocities up to $\ensuremath{\sim}5\text{ }\text{cm}∕\text{s}$. A border between quasi-elastic collisions and merger of the solitons into a single localized state is identified. In some cases, the soliton produced by the merger is stable against collapse, which was not observed before in the static OL potential either.
We investigate stability of gap solitons (GSs) in the first two band gaps in the framework of the one-dimensional Gross-Pitaevskii equation, combining the repulsive nonlinearity and a moderately strong optical lattice (OL), which is subjected to ``management,'' in the form of time-periodic modulation of its depth. The analysis is performed for parameters relevant to the experiment, characteristic values of the modulation frequency being $\ensuremath{\omega}\ensuremath{\sim}2\ensuremath{\pi}\ifmmode\times\else\texttimes\fi{}20\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$. First, we present several GS species in the two band gaps in the absence of the management. These include fundamental solitons and their bound states, as well as a subfundamental soliton in the second gap, featuring two peaks of opposite signs in a single well of the periodic potential. This soliton is always unstable, and quickly transforms into a fundamental GS, losing a considerable part of its norm. In the first band gap, (stable) bound states of two fundamental GSs are possible solely with opposite signs, if they are separated by an empty site. Under the periodic modulation of the OL depth, we identify stability regions for various GS species, in terms of $\ensuremath{\omega}$ and modulation amplitude, at fixed values of the soliton's norm,
N. In either band gap, the GS species with smallest
N has a largest stability area; in the first and second gaps, they are, respectively, the fundamental GS proper, or the one spontaneously generated from the subfundamental soliton. However, with the increase of
N, the stability region of every species expands in the first gap, and shrinks in the second one. The outcome of the instability development is also different in the two band gaps: it is destruction of the GS in the first gap, and generation of extra side lobes by unstable GSs in the second one.