We consider possibilities to control dynamics of solitons of two types, maintained by the combination of cubic attraction and spin-orbit coupling (SOC) in a two-component system, namely, semi-dipoles (SDs) and mixed modes (MMs), by making the relative strength of the cross-attraction, γ, a function of time periodically oscillating around the critical value, γ = 1, which is an SD/MM stability boundary in the static system. The structure of SDs is represented by the combination of a fundamental soliton in one component and localized dipole mode in the other, while MMs combine fundamental and dipole terms in each component. Systematic numerical analysis reveals a finite bistability region for the SDs and MMs around γ = 1, which does not exist in the absence of the periodic temporal modulation ("management"), as well as emergence of specific instability troughs and stability tongues for the solitons of both types, which may be explained as manifestations of resonances between the time-periodic modulation and intrinsic modes of the solitons. The system can be implemented in Bose-Einstein condensates, and emulated in nonlinear optical waveguides.PACS numbers:
I. INTRODUCTIONAn active direction in the current work with Bose-Einstein condensates (BECs) in atomic gases is using them as a testbed for emulation of various effects originating in condensed-matter physics, which may be reproduced in a clean and easy-to-control form in ultracold bosonic gases [1][2][3]. In particular, a binary gas, with a pseudo-spinor twocomponent wave function, may emulate spin-orbit coupling (SOC) in semiconductors, i.e., the interaction between the electron's spin and its motion across the underlying ionic lattice [4,5], as first demonstrated in Ref. [6], see also reviews [7]- [11]. While most experimental works on the BEC simulation of SOC dealt with effectively one-dimensional (1D) settings, implementation of SOC in the quasi-2D geometry was reported too [12], making it relevant to consider 2D (and 3D) systems coupled by the spin-orbit interaction. In this way, SOC opens a straightforward way to the creation of topological modes characterized by vorticity, because linear operators accounting for the coupling of two components in the corresponding system of Gross-Pitaevskii equations (GPEs), see Eq. (1) below, generate vorticity in one component if the other one is taken in the zero-vorticity form.The SOC effect, being linear by itself, may be naturally combined with the intrinsic nonlinearity of bosonic gases, represented by cubic terms in the respective GPEs [13] (and/or by nonlocal cubic terms accounting for long-range interactions in BEC built of dipole atoms [14]). The interplay of SOC and nonlinearity makes it possible to predict a great variety of stable modes, including 1D and 2D solitons [15][16][17][18] and various nonlinear topological states in 2D, such as vortices and vortex lattices [19]-[28] and skyrmions [29]. In fact, the 2D and 3D SOC systems is one of the most prolific sources of nonlinear states with intrinsic topologic...