Four-wave mixing in spin-orbit coupled Bose-Einstein condensates

TL;DR

This paper explores four-wave mixing in spin-orbit coupled Bose-Einstein condensates, identifying conditions for different configurations, demonstrating unique dynamics through simulations, and highlighting potential experimental applications in quantum physics.

Contribution

It derives phase matching conditions for FWM in spin-orbit coupled BECs and demonstrates possible configurations and dynamics through numerical simulations.

Findings

Identified four distinct FWM configurations based on spinor states.

Discovered a unique configuration with probe waves slower than the pump.

Showed simultaneous occurrence of two FWM processes.

Abstract

We describe possibilities of spontaneous, degenerate four-wave mixing (FWM) processes in spin-orbit coupled Bose-Einstein condensates. Phase matching conditions (i.e., energy and momentum conservation laws) in such systems allow one to identify four different configurations characterized by involvement of distinct spinor states in which such a process can take place. We derived these conditions from first principles and then illustrated dynamics with direct numerical simulations. We found, among others, the unique configuration, where both probe waves have smaller group velocity than pump wave and proved numerically that it can be observed experimentally under proper choice of the parameters. We also reported the case when two different FWM processes can occur simultaneously. The described resonant interactions of matter waves is expected to play important role in the experiments of BEC with artificial gauge fields. Beams created by FWM processes are important source of correlated particles and can be used the experiments testing quantum properties of atomic ensembles.

Tables

Table 1. Table 1: Possible configurations of degenerate FWM processes (positive and negative q 𝑞 q are included). The first column are numbers identifying configurations, which corresponds to the specific choice of the spectrum branches for the pump and probe waves, indicated in second, third, and fourth columns, respectively. In the last column we show the maximal number, N max subscript 𝑁 max N_{\rm max} , of q 𝑞 q values solving Eq. ( 12 ).

Configuration	$s_1$	$s_2$	$s_3$	$N_{\max }$
1	$1$	$-1$	$-1$	2
2	$1$	$1$	$-1$	2
3	$1$	$-1$	$1$	2
4	$-1$	$-1$	$-1$	4

Equations

$$i{\partial }_t\mathbf{\Psi }=H\mathbf{\Psi }+\frac{1}{2}G(\mathbf{\Psi })\mathbf{\Psi },$$

$$H=\frac{1}{2}\left(-{\partial }_x^2+\mathbf{\Omega }\cdot \sigma -i\alpha {\sigma }_x{\partial }_x\right)$$

$$G(\mathbf{\Psi})=\left(\begin{array}[]{c c}g_{1}|\Psi_{1}|^{2}+g|\Psi_{2}|^{2}&0\\ 0&g|\Psi_{1}|^{2}+g_{2}|\Psi_{2}|^{2}\\ \end{array}\right)$$

$${\mathbf{\Psi }}_{\pm }(x,t)=e^{ikx-i{\mu }_{\pm }(k)t}{\psi }_{\pm }(k),$$

$${\mu }_{\pm }(k)=\frac{k^2}{2}\pm \frac{\epsilon (k)}{2}.$$

$$\epsilon (k)=\sqrt{{\Omega }_z^2+{\tilde{\Omega }}^2(k)}$$

$$\mathbf{\psi}_{\pm}(k)=\frac{1}{\sqrt{\widetilde{\Omega}^{2}+(\varepsilon(k)\mp\Omega_{z})^{2}}}\left(\begin{array}[]{c}\widetilde{\Omega}\\ \pm\varepsilon(k)-\Omega_{z}\end{array}\right).$$

$$2k_1=k_2+k_3,$$

$$2{\mu }_{{\nu }_1}(k_1)={\mu }_{{\nu }_2}(k_2)+{\mu }_{{\nu }_3}(k_3).$$

$$\Pi ={\int }_{-\infty }^{\infty }{\mathbf{\Psi }}^{†}\left(-i{\partial }_x+\frac{\alpha }{2}{\sigma }_x\right)\mathbf{\Psi }dx$$

$$k_2=k_1+q,k_3=k_1-q.$$

$$2q^2=2s_1\epsilon (k_1)-s_2\epsilon \left(k_1+q\right)-s_3\epsilon \left(k_1-q\right),$$

$$0\le q^2\le {\epsilon }_{+}(k_1-q)\le {\epsilon }_{+}(k_1+q),$$

$$Q^3-\left(1+4s_1\sqrt{{\tilde{\omega }}^2+{\omega }_z^2}\right)Q^2$$

$${\Delta }_{s_1}={\tilde{\omega }}^2\left[15{\tilde{\omega }}^2+4s_1\sqrt{{\tilde{\omega }}^2+{\omega }_z^2}\left({\tilde{\omega }}^2+{\omega }_z^2+3\right)-12{\omega }_z^2-4\right].$$

$$v_{\pm }(k)\equiv \frac{\partial {\mu }_{\pm }(k)}{\partial k}=k\pm \frac{\tilde{\Omega }}{2\sqrt{{\tilde{\Omega }}^2+{\Omega }_z^2}},$$

$$\mathbf{\Psi }(x,t=0)=e^{-x^2/w^2}\sum _{j=1}^3{A}_j{\psi }_{s_j}(k_j)e^{i{k}_{j}x}.$$

$$k_2=k_1+q_1,k_3=k_1-q_1$$

$$k_4=k_1+q_2,k_5=k_1-q_2.$$

$$\eta (k_3)=\frac{\tilde{N}(k_3)}{N}\cdot 100\%,$$

Full text

Keywords: Spin-orbit coupling, four-wave mixing, Bose-Einstein condensate.

1 Introduction

Traditionally four-wave mixing (FWM) processes are associated with photon interactions via a non-linear polarization. It is a third-order parametric process in which two particles (from two so-called writing or pump beams - one particle from each beam) are annihilated when passing through a non-linear medium, and at the same time two new particles (constituting probe and signal beams) are generated. In optics, FWM is commonly associated with the third order Kerr nonlinearity. The phenomenon is ubiquitous (see e.g.,  [1, 2]) and its applications are very widespread, including: fiber optic communication (very often not welcome), wavelength conversion, parametric amplification, optical regeneration, optic phase conjugation and correction of the aberration of images.

FWM can be observed also for massive particles as it was predicted and observed in cold atomic gasses two decades ago [3]-[5] (see also [6]). In this case resonantly interacting particles are neutral atoms rather than photons. Flexibility of control of trapping potentials, as well as of nonlinear interactions, in atomic systems open interesting perspectives of managing both momentum conservation and energy conservation laws through the interplay of additional linear and nonlinear potentials. This issue was already explored in Refs. [7]-[11]. In this paper we explore a similar idea of controlling FWM processes through artificially created gauge potentials.

In order to consider FWM in a specific medium, one has to identify the characteristic eigenmodes: the elementary solutions to the linearized equations of motion in the form of plane waves, identified by their wavevectors and frequencies which satisfy the so-called phase matching conditions, that are equivalent to momentum and energy conservation laws. These are often quite demanding constraints depending on the particular form of dispersion relation characteristic for the system under investigation. For instance, in one dimension they cannot be satisfied for a system of cold atoms obeying parabolic dispersion relation and confined to (quasi-)one dimension. The situation can be improved by artificial change of the dispersion law using linear optical lattices [7, 8, 9] or by the manipulation of the wavenumbers of the matter waves involved in the process by means of nonlinear lattices [10, 11]. These modifications introduce the internal texture to the propagation medium, making it inherently inhomogeneous.

If a system has a spinor nature, i.e., consists of two subsystems, an alternative way to manipulate the linear properties of the medium, even preserving homogeneity, is to employ coherent coupling of the constituents. In optics, for example, one can satisfy the matching conditions for the FWM of light propagating in homogeneous parity-time ($\mathcal{PT}$) - symmetric coupled waveguides with gain and losses [12]. One can also consider the matching conditions, and thus observation of the FWM in $\mathcal{PT}$-symmetric optical lattices which are available experimentally [13]. And yet another interesting application was demonstrated in multi-component vector solitons consisting of two perpendicular FWM dipole components created by electromagnetically induced gratings [14].

Similar situation naturally occurs for spinor Bose-Einstein, where coupling between two atomic states by means of the spin-orbit coupling (SOC) allows one to manipulate the dispersion relation in the presence of external potential. This idea becomes attractive, since using various experimental techniques, spin-orbit-coupled Bose-Einstein condensates (SOC-BECs) of hyperfine states of 87Rb atoms has been created [15, 16, 17, 18, 19, 20]. Notice that atoms of s and p bands of the static lattice were considered as pseudospins [21].

The main goal of this paper is to show that with properly adjusted SOC, one can satisfy the phase matching conditions for a homogeneous one-dimensional SOC-BECs. Importantly, the SOC properties in atomic systems are highly adjustable [22, 23, 24, 25] and the matching conditions reported below are experimentally feasible. Inter-atomic interactions are also tunable, most commonly by Feshbach resonance; see Ref. [26] for the observation of Fesbach resonances for SOC fermions, or Ref. [27] for observation of partial waves, with nonlinear interactions controlled by the dressing technique. Additional degrees of freedom in manipulating the effective dispersion relation, and thus, of the matching conditions [see (8) and (9) below], may be reached by using moving lattices [17, 21].

Understanding dynamics of FWM process is important in the context of creation od pairs of correlated particles, which are crucial in the ultra-precise metrology. Particles undergoing elastic collisions may create entanglement and it can be of practical use in the experiment. Here we mention the experiments [30, 31, 32], study of Cauchy Schwarz inequality [33, 34], and twin beam generation [35]. Also, recently there was a prediction [36] and its experimental confirmation [37], of the Bell inequality is violated in a many-body system of massive particles with spinor condensates.

In the context of present considerations we would like to mention two recent experimental achievements. In the first a stripe phase with supersolid properties in SOC-BECs has been observed [28]. In this case spin flip process with a momentum transfer has been realized and observed using Bragg scattering. In another experiment radio-frequency (rf) photons were dressed with tunable recoil momentum by combining rf pulses with an oscillating magnetic force. This leads to a new application of Floquet engineering: periodically driven systems can have time-averaged properties which cannot be achieved with constant fields and in our opinion hold a promise of mixing different waves [29].

Due to the spinor nature of the one-dimensional (1D) SOC-BEC, it is characterized by two branches of the dispersion relation. As a result, the matching conditions can be readily satisfied, as we shall see below. Moreover, unlike in the case BECs without SOC, where additional external potentials, like linear [7, 8, 9] or nonlinear [10, 11] lattices are employed for ensuring matching conditions, now one can find a diversity of distinct FWM processes, where the interacting waves, as well as waves generated may represent different spinor states at the same values of parameters of the system (similarly to the FWM with laser pulses reported in [12]). It is a goal of the present study to identify the FWM processes available in the FWM in SOC-BECs and show they can be efficient enough to be observed.

The organization of the paper is as follows: in Sec. 2 we discuss phase matching conditions in a quasi-1D SOC-BEC, and we focus on the degenerate case, where the central pump wave serves as a source for stimulated enhancement of two probe waves. Here we borrow the terminology from optics but we make no distinction between signal and probe beams. Next we identify four possible configurations of FWM and in Sec. 3 we perform feasibility study using real time simulations for all predicted configurations. The outcomes are summarized in Conclusion.

2 Phase matching conditions

2.1 General relations

Let us consider a quasi-1D SOC-BEC which is described by a two-component order parameter $\mathbf{\Psi}(x,t)=(\Psi_{1}(x,t),\Psi_{2}(x,t))^{\mathrm{T}}$ (hereafter $\mathrm{T}$ stands for transposition). The dynamics of the spinor $\mathbf{\Psi}(x,t)$ is governed by the coupled Gross-Pitaevskii equations (GPEs):

[TABLE]

where

[TABLE]

is the linear mean field Hamiltonian of the two component BEC, $\alpha$ is the SOC strength, $\mathbf{\Omega}$ is the vector of the Zeeman coupling (we admit external magnetic field), $\mathbf{\sigma}=(\sigma_{x},\sigma_{y},\sigma_{z})$ is the vector of Pauli matrices $\sigma_{x,y,z}$. Physical interpretation of these parameters depends on the particular realization. For instance in the experiment performed in an 87Rb Bose-Einstein condensate, a pair of Raman lasers created a momentum sensitive coupling between two internal atomic states [15]. This SOC was equivalent to that of an electronic system with equal contributions of Rashba and Dresselhaus couplings, and with a uniform magnetic field $B$ in the $(x,z)-$plane. In materials the SOC is due to intrinsic properties, which are largely determined by the specific material and the details of its growth. In these and other proposed schemes $\mathbf{\Omega}$ may have different physical interpretation, including Rabi frequencies of the dressing laser fields. It has to be taken into account when one consider tuning and range of vector $\mathbf{\Omega}$. Here we will call it Zeeman coupling, and $\Omega_{z}$ Zeeman splitting.

The nonlinearity $2\times 2$ matrix is given by

[TABLE]

with intra– and inter–component interactions $g_{1,2}$ and $g$, respectively, and we use the dimensionless units with $\hbar=m=1$.

To address the matching conditions for the FWM we start with the eigenmodes of the linear spectral problem, representing them in the form of the plane waves

[TABLE]

where $\mathbf{\psi}_{\pm}(k)=(\psi_{\pm}^{(1)}(k),\psi_{\pm}^{(2)}(k))^{\rm T}$ is a constant (i.e. $x$- and $t$-independent) spinor, $k$ is a mode wavenumber, $\mu_{\pm}(k)$ is its frequency, and $\pm$ indicate the upper (”$+$”) and lower (”$-$”) branches of the spectrum [i.e., $\mu_{-}(k)\leq\mu_{+}(k)$]. We will concentrate in the case of the magnetic field in the $(x,z)$ plane, i.e. $\mathbf{\Omega}=(\Omega_{x},0,\Omega_{z})$ [Without loss of generality we fix $\Omega_{x},\Omega_{z}\geq 0$], for which we compute the two branches of the dispersion relation

[TABLE]

Here

[TABLE]

with $\widetilde{\Omega}(k)=\alpha k+\Omega_{x}$, is the gap between the spectral branches at a given $k$. The lower (-) and the upper (+) branches of the respective eigenvectors are defined as

[TABLE]

It is to be emphasized that (4), as well as the matching conditions (8), (9), strictly speaking are valid only in linear regime, or approximately valid during the initial stages of the evolution. If however, the pulses with newly generated central frequencies separate in space from the initial pulse due to different group velocities, further generation of new harmonics is suppressed and one can observe propagation of quasi-monochromatic pulses. Confirmation of the described scenario requires numerical simulation of the full nonlinear problem. This approach is adopted below.

To reduce the number of parameters, here we investigate the degenerate FWM process, where two input spinor states are identical. We label their wavevectors by $k_{1}$ and using the optical terminology, we call them pump waves. Two spinors that are created in the FWM process with central wavevectors $k_{2}$ and $k_{3}$, will be referred to as probe waves. Respectively, the conservation of wavenumbers and frequencies of the pump and probe waves, are expressed in the from of phase-matching conditions

[TABLE]

The indexes $\nu_{j}$ ($j=1,2,3$) refer to either “$+$” or “$-$” branch of the spectrum.

Below we consider only the cases $k_{2,3}\neq k_{1}$, excluding the trivial case of self-phase modulation where all wavenumbers are equal. We note that the system (1) obeys gauge, rather than Gallilean, invariance, at which the generalized momentum

[TABLE]

is conserved: $d\Pi/dt=0$. For our consideration this means that the input wavenumber $k_{1}$ cannot be set arbitrarily to zero without changing the spinor eigenstates. It also means that at zero Zeeman field $\mathbf{\Omega}=0$, the linear Hamiltonian $H$ is gauge equivalent to the usual one-dimensional Schrödinger Hamiltonian $H_{0}=-\partial_{x}^{2}$ which does not support the matching conditions (8) and (9). In other words, while SOC controls the waves involved in resonant processes, the FWM itself requires nonzero Zeeman field.

For the following consideration it is convenient to rewrite (8) as

[TABLE]

Now matching condition for frequencies (9) can be rewritten in the form

[TABLE]

where $s_{j}=\pm 1$. Since each wave belongs to either upper or lower branch, these are eight different equations for given $k_{1}$ and $q$. However only four of them have nontrivial solutions. To justify this we first notice that if $s_{2}=s_{3}$, the condition (12) is symmetric under $q\leftrightarrow-q$ exchange. If however $s_{2}\neq s_{3}$, then (12) is symmetric with respect to simultaneous change $(s_{2},q)\leftrightarrow(s_{3},-q)$. This allows one to restrict the analysis to the case $q>0$. Next we use $\varepsilon(k)>0$ property and conclude that the case $(s_{1},s_{2},s_{3})=(-1,1,1)$ does not have solutions since the right hand side of Eq. (12) becomes negative. Let us now consider $k_{1}\geq 0$ (the case $k_{1}<0$ is fully analogous). For non-negative values of $k_{1}$ we find the following inequalities

[TABLE]

excluding the cases $(-1,1,-1)$ and $(-1,-1,1)$. Hence, the initial pulse from the lower brunch of the spectrum may originate degenerate FWM in processes involving modes from the lowest branch only (this is the configuration 4 in the Table 1 below). Finally, using inequalities (13) one can exclude also the case $(1,1,1)$, leaving only four possible configurations summarized in the Table. 1.

In the last column of the Table. 1 we list the maximal number of solutions for particular configuration. In what follows we present analytical and graphic considerations that led us to these counts.

2.2 Analysis of possible configurations

By straightforward algebraic manipulations we can eliminate square-root terms in the equation (12) (simple sequence of transfers and squaring). As a result all four phase matched processes listed in the Table 1 are determined by the following cubic equations

[TABLE]

where $\widetilde{\omega}=\widetilde{\Omega}/\alpha^{2}$, $\omega_{z}=\Omega_{z}/\alpha^{2}$ and $Q=q^{2}/\alpha^{2}$. For obvious reasons we are interested only in positive roots of (14) and exclude the root $Q=0$ (i.e. $q=0$) which does not correspond to FWM but to the self-phase modulation.

A number of real roots of Eq. (14) is determined by the sign of the discriminant

[TABLE]

If $\Delta_{s_{1}}>0$, there exists one real root. Three distinct real solutions exist if $\Delta_{s_{1}}<0$. At $\Delta_{s_{1}}=0$ all roots are real and at least one is multiple [38].

Now we inspect systematically configurations listed in Tab. 1 which manifest qualitatively different types of dynamics. Starting with the last one, we set $s_{1}=-1$. Now the discriminant (15) depends on the two parameters $\{\widetilde{\omega},\omega_{z}\}$ for different values of which $\Delta_{s_{1}}$ can acquire any sign or be zero. Analyzing Vieta formulae one can exclude the possibility of all three real roots being positive [38]. It means that in the configuration 4 (see the Table 1) there may exist either one or two real positive roots of Eq. (14). Taking into account that roots appear in pairs, $\pm q$, this corresponds to at most four possible arrangements allowed by the phase matching condition (12).

Next we turn to the configurations 1, 2, and 3 in the Table 1, by setting $s_{1}=+1$ in (14) and in (15). In this case we find that $\Delta_{s_{1}}<0$ for all values of the parameters $\{\widetilde{\omega},\omega_{z}\}$, i.e. all roots are real. Moreover, they are all positive as follows from the Vieta formulae. The three real positive roots of $Q$, i.e. six roots $q$, correspond to the three different configurations [notice that the difference among these configurations was removed upon squaring in obtaining Eq. (14)].

It is easy to establish one-to-one correspondence of the roots and the configurations. For if one root is common for two configurations we use equation (12) to show that one of the terms $\varepsilon_{s}(k_{1}\pm q)$ is equal to zero. Obviously, for each of the configurations with $s_{1}=+1$ always there is at least one positive root (all roots cannot be negative, because their product is positive). Let $0<q_{1},q_{2},q_{3}$ denote the positive roots, so that the configuration 1 has two symmetric roots, which we denote by $\pm q_{1}$ while the roots of the configurations 2 and 3 are given by ($q_{2}$, $-q_{3}$) and ($-q_{2}$, $q_{3}$), respectively.

These arguments are exemplified in Fig. 1 where we show graphical solutions of phase matching equation. In the panels (a), obtained for $\alpha=2$, and (b), obtained for $\alpha=10$, the dashed blue curve is the plot of the left hand side (LHS) of Eq. (12). Right hand side (RHS) of this equation is represented by red, pink and black lines correspond to configurations 1, 2 and 3 in Tab. 1. Crossings of solid and dashed lines yield the roots of the phase matching equations. Comparing the values of the LHS and RHS of Eq. (12) at $q=0$ and at $q\to\infty$, we conclude that each of the configurations 1, 2 and 3 must have at least one root for $q>0$.

The fourth configuration ($s_{1}=s_{2}=s_{3}=-1$) is illustrated in panel (c). Here we observe three different possibilities: no positive solutions of the phase matching condition (12) (the light green curve does not cross the dash blue curve at $\alpha=5$); one positive root (the deep green and dash blue curves tangent to each other when $\alpha\approx 7.67$), and two positive solutions (the brown curve crosses the dash blue curve in two points, when $\alpha=9$).

Finally, in the panel (d) we varied $k_{1}$ (while holding $\alpha=9$) and shown the regions where the phase matching equation of the fourth configuration supports one positive (the black curves) and two positive solutions (the pink curves). In all panels (a)-(d) we fixed $\Omega_{x}=2.5$ and $\Omega_{z}=8$.

2.3 Matching of group velocities

While no matching conditions on the group velocities is imposed, for practical observation of different scenarios of FWM in numerical simulations, the issue of the group velocities (GVs)

[TABLE]

becomes relevant. On the one hand where all wavepackets move with respect to each other, it is important that the spinor involved in the process have similar values of GVs: otherwise fast separation of wavepackets in space may drastically reduce the conversion efficiency. On the other hand, GVs should have sufficient difference in order to observe spatial separation of the probe wavepackets. Thus in addition to solving the matching conditions we set a task of finding optimal conditions in the context of FWM numerical simulations (they are presented in the next section).

For each configuration listed in Table 1, at a given $k_{1}$ one can determine $q$, i.e., the wavenumbers $k_{2}$ and $k_{3}$, and consequently their GVs.

In Fig. 2 (a) we observe that the GV of the pump wavepacket, $v_{+}(k_{1})$, is close to either $v_{-}(k_{2})$ or $v_{-}(k_{3})$, almost for all $k_{1}$ except the vicinity of $k_{1}=0$. This means that to obtain clear separation of the wavepackets, generated in the FWM at relatively short time intervals, $k_{1}$ should be chosen close to zero. Then, the separation between velocities grows rapidly enough allowing direct observation of separated pulses. On the other hand, the time that pulses overlap is still long enough to generate a substantial four wave mixing signal.

Fig. 2 (b) shows the dependence of GVs of phase matched wavepackets versus pump momentum for configurations 2 and 3 from the Table 1. We again observe that in some regions GVs are close to each other or even coincide what does not allow observation of separation of the generated wavepackets from the initial one. However an interesting situation occurs in the vicinity of $k_{1}=0$. Here GVs of the second and third waves have bigger absolute values than $v_{+}(k_{1})$ and the same sign. In this case both created waves move faster that the initial wavepacket. We should emphasize that this is not common in the usual realizations of FWM and this is solely due to the SOC coupling. Note that in these configurations it is also possible to initiate FWM process where one of the velocities of the probe waves is smaller and one bigger than that of the pump.

The situation is more complicated for the case of the fourth configuration in Table 1, as shown in Fig. 2 (c). The (thick) red curve represents GV of the pump wave of the negative branch $v_{-}(k_{1})$ [see Eq. (16)]. The other (black) curves, that have forms of three loops, represent GVs of generated waves [$v_{-}(k_{2})$ and $v_{-}(k_{3})$] that correspond to other (non-trivial) solutions. Like in the previous cases, to reach significant separation of the pulses in the real space we choose $k_{1}$ in a region far from the crossing of the curves.

3 Numerical results

Equipped with the solutions of matching conditions and with the ideas of optimization the conversion efficiency in terms of the GVs we now turn to direct numerical simulations of the configurations of the FWM processes summarized in Table 1. In order to find favorable conditions to observe clear evidence of specific FWM process, first one has to select proper momentum $k_{1}$. Note that phase matching will automatically determine all participating wavepackets GVs as explained in the previous section. Then, appropriate initial widths and amplitudes of the pump and probe waves need to be adjusted to ensure long enough and strong enough nonlinear interaction.

In all simulations we use the wavepackets having equal widths and completely overlapping at $t=0$, i.e.,

[TABLE]

Here $A_{j}$ are initial amplitudes of the wave-packets with the central wavevectors $k_{j}$ of the spinors defined in accordance with Eq. (7).

For the FWM process corresponding to configuration 1, the initial state is formed with $A_{1}=1$, $A_{2}=0.2$, $A_{3}=0$ and $s_{j}$ are chosen according to Table 1: $s_{1}=1$, $s_{2}=-1$ and $s_{3}=-1$. In Fig. 3 we show an example of the FWM for this configuration. Due to the FWM process, by the expense of the highly populated initial state $A_{1}$ we observe strong amplification of the seed state and growth of the third matter wave with phase matched momentum $k_{3}$.

This process is depicted with snapshots at the beginning and end (i.e., at $t=0$ and $t=300$) of the simulations in Fig. 3, where main panels (a), (c) [(b), (d)] refer to the first [second] spinor component. In particular blue contours in panels (a), (b) represent initially overlapping pump ($k_{1}$) and probe ($k_{2}$) waves in the configuration space. They are fully separated after evolution time ($t=300$) due to the difference in GVs and new, clearly visible, wave of central momentum $k_{3}$ is generated. Panels (c), (d) show the corresponding features in the Fourier space. Here we distinct two waves as narrow blue peaks, at initial time and again three waves at the end of evolution. The most explicit feature is substantial broadening of all participating matter waves during the evolution. The inset in each panel shows full time evolution of modulus of the spinor components - (a) and (b) in the real space and (c) and (d) in the momentum space.

In the next two figures we illustrate the FWM process corresponding the second and third configuration from Table 1, with GV configuration shown in panel (b) of the Fig. 2. As mentioned above in these two configurations the exist two different roots ($q_{2}\neq q_{3}$) of phase matching condition (12). As one can see directly in panels (a) and (b) of Fig. 1 these configurations are related by the transformation $q_{1}\rightarrow-q_{2}$ and $q_{2}\rightarrow-q_{1}$, i.e. the analysis of second and third configuration are analogous.

Interestingly, in Fig. 4 both probe and created waves are generated in the same side of the pump wave. The evolution of the probe wavepackets in Fig. 5 looks qualitatively similar to that one shown in Fig. 3. However, since each newborn wavepacket bears a quasi-spin, the emergent spinors (more precisely the left propagating waves) are different in these cases.

Turning to the fourth in the Table (1) , we recall that in this case one can obtain up to five solutions from the phase matching condition (including the trivial case of $q=0$). We start with Fig. 6, where initial group velocities can be identified in the panel (c) of the Fig. 2 and are marked as dots on red and black curves. In principle, the dynamics presented in this case is very similar to that shown for the configurations 1 and 3, except that now different spinor states are involved (respectively the wavepackets bear different quasi-spins). Also closer look at Fig. 6 (c) and Fig. 6 (d), showing the spectra of the components reveals an interesting feature. Namely one can spot oscillations of the amplitude of pump wavepacket and we attribute them to the self phase modulation which was mentioned above as the trivial solution $Q=0$ (or $q=0$) of the Eq. (14).

For completeness we present the fourth configuration in the case when the phase matching allows simultaneously for two FWM processes. In this last simulation we used parameters corresponding to the region, where the phase matching equation supports four different root: $q_{1}$, $q_{2}$, $-q_{1}$ and $-q_{2}$. Having these roots in hand, we define momenta of the two sets of probe waves in the following way:

[TABLE]

In present case dynamics of the FWM process involves waves with the amplitudes: $A_{1}$ as a pump, and two pairs of probs $A_{2}$ and $A_{3}$, as well as $A_{4}$ and $A_{5}$ (by convention a mode with $Aj$ has momentum $k_{j}$, where $j=1,...,5$). In numerical simulations we initiate the dynamics by putting $A_{1}=1$, $A_{3}=A_{5}=0.2$, and $A_{2}=A_{4}=0$ [see Eq.(17); the values of the rest of parameters are listed in the Fig. 7 caption]. In the nonlinear evolution the pump is interacting with both sets of the probe wavepackets creating two new waves $A_{2}$ and $A_{4}$ with momenta $k_{2}$ and $k_{4}$ respectively, in two simultaneous FWM processes. Figure 7 shows two components of the spinor wavefunctions in the configuration [(a) and (b)] and momentum [(c) and (d)] spaces, where we can clearly see two sets of probe waves. This time we can not refer to Fig. 2, since the values of parameters were slightly different from those used for panel(c), but for the sake of clarifications we propose to look at panel (d) of the Fig. 1. Due to the strong spreading there is substantial overlap of wavepackets, but when we look at them in the momentum space all peaks can be clearly identified. Close inspection reveals also oscillations on the pump wave due to the self-phase modulation mentioned above.

In general finding optimal conditions to observe effective FWM process is not an easy task especially taking into account that the problem is multi-parametric and the values of parameters have to be carefully selected. In order to support this statement we have chosen the Zeeman splitting $\Omega_{z}$ as a control parameter and in Fig. 8 we present the efficiency (unit of per cent) of FWM defined as

[TABLE]

where $N$ is total number of atoms in the system and the number of atoms corresponding to the generated wave $\widetilde{N}(k_{3})$ is evaluated after long enough evolution time, when the wave packets are well separated. The efficiency is calculated for the FWM conversion in the process identified above as configuration 1 for different for different types of the nonlinearity parameterized by the quantities $\overline{g}=(g_{1}+g_{2})/2$ and $\Delta g=(g_{1}-g_{2})/2$. In all of the calculations presented in the Fig. 8 we took $g=\overline{g}$.

In both panels of the Fig. 8 we vary $\Omega_{z}$ while keeping the rest of the parameters fixed (see Fig. 3 for the exact values of all parameters). In the left panel different curves correspond to different values of the ratio $\Delta g/g$ with $g=0.8$, while in the right panel $\Delta g$ is equal to zero and different curves were obtained choosing different values of the interaction strength $g$. In both panels we observe non-monotonic behavior of the efficiency on $\Omega_{z}$. At $g_{1}=g_{2}=g$ (the right panel and the red curve in the left panel) we find that the efficiency goes to zero in the limits $\Omega_{z}\to 0$ and $\Omega_{z}\to\infty$ having maximum at some finite value of the Zeeman splitting (notice that both these limits correspond to the integrable cases [39]). Such behavior may be explained by the fact that the case of all equal nonlinearities corresponds to the SU(2) symmetric (also known as Manakov [40]) nonlinearity. Thus at $\Omega_{z}=0$ the linear Hamiltonian (2) can be diagonalized by the global rotation. An approximate diagonalization can be also performed in the formal limit $\Omega_{z}\to\infty$ at $\Omega_{x}$ fixed, which after rescaling can be viewed as the limit of $\Omega_{x}\to 0$ at $\Omega_{z}$ fixed. Next we recall, that for a 1D NLS equation the matching conditions cannot be satisfied, this meaning that no FWM processes involving eigenstates from the same branch can be observed. On the other hand, the SU(2) nonlinearity does not support FWM processes where states from the different branches are involved. This can be easily seen if we rewrite the Manakov nonlinearity in the form $G(\bm{\Psi})\equiv\bm{\Psi}^{\dagger}\bm{\Psi}$. Indeed, since the states belonging to different branches of the spectrum are mutually orthogonal, this nonlinearity does not support mixing of mentioned states, i.e. the respective frequency conversion processes is not phase matched. On the other hand, when the nonlinear coefficients are not exactly equal (blue and green lines in the left panel of Fig. 8) we observe that even at $\Omega_{z}\to 0$ the efficiency does not vanish, although becomes very small, as shown in the inset in the left panel of Fig. 8.

4 Conclusions

In our study we analyzed the four-wave mixing process in Bose-Einstein condensates with spin-orbit coupling. We found all phase matched configurations for degenerate case where two identical initial states interact with two probe ones. We performed numerical simulations to illustrate the dynamics in which we seeded one of the probe and observed stimulated growth of the latter combined with resonant generation of extra waves. We found unique conditions where both probe waves have smaller group velocity than pump wave, and also reported the case when two FWM process can occur simultaneously.

5 Acknowledgments

The work was supported by the Polish National Science Centre 2016/22/M/ST2/00261 (M.T.), Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.01-2017.55 (N.V.H.), and Portuguese Foundation for Science and Technology (FCT) under Contract no. UIDB/00618/2020 (V.V.K.).

6 References