Localized spatially nonlinear matter waves in atomic-molecular Bose-Einstein condensates with space-modulated nonlinearity

Abstract

The intrinsic nonlinearity is the most remarkable characteristic of the Bose-Einstein condensates (BECs) systems. Many studies have been done on atomic BECs with time- and space- modulated nonlinearities, while there is few work considering the atomic-molecular BECs with space-modulated nonlinearities. Here, we obtain two kinds of Jacobi elliptic solutions and a family of rational solutions of the atomic-molecular BECs with trapping potential and space-modulated nonlinearity and consider the effect of three-body interaction on the localized matter wave solutions. The topological properties of the localized nonlinear matter wave for no coupling are analysed: the parity of nonlinear matter wave functions depends only on the principal quantum number n, and the numbers of the density packets for each quantum state depend on both the principal quantum number n and the secondary quantum number l. When the coupling is not zero, the localized nonlinear matter waves given by the rational function, their topological properties are independent of the principal quantum number n, only depend on the secondary quantum number l. The Raman detuning and the chemical potential can change the number and the shape of the density packets. The stability of the Jacobi elliptic solutions depends on the principal quantum number n, while the stability of the rational solutions depends on the chemical potential and Raman detuning.

Localized spatially nonlinear matter waves in atomic-molecular Bose-Einstein condensates with space-modulated nonlinearity

OPEN

A

Yu-QinYao, Ji Li, Wei Han, Deng-ShanWang & Wu-Ming Liu

The intrinsic nonlinearity is the most remarkable characteristic of the Bose-Einstein condensates (BECs) systems. Many studies have been done on atomic BECs with time- and space- modulated nonlinearities, while there is few work considering the atomic-molecular BECs with space-modulated nonlinearities. Here, we obtain two kinds of Jacobi elliptic solutions and a family of rational solutions of the atomic- of three-body interaction on the localized matter wave solutions. The topological properties of the localized nonlinear matter wave for no coupling are analysed: the parity of nonlinear matter wave functions depends only on the principal quantum number n, and the numbers of the density packetsfor each quantum state depend on both the principal quantum number n and the secondary quantum number l. When the coupling is not zero, the localized nonlinear matter waves given by the rational function, their topological properties are independent of the principal quantum number n, only depend on the secondary quantum number l. The Raman detuning and the chemical potential can change the number and the shape of the density packets. The stability of the Jacobi elliptic solutions depends on the principal quantum number n, while the stability of the rational solutions depends on the chemical potential and Raman detuning.

It is known that the precise control of untracold atomic systems have brought the realization of Bose-Einstein condensates (BECs) and Fermi gases. An important challenge is to produce and control more complicated molecular systems because of their potential applications for the tests of fundamental physics and for the dris of fundamental constants. To date, several atomic-molecular conversion schemes have been provided110. Among

them, Feshbach resonance5,6 and photoassociation7,8 are two main techniques to produce cole molecules from an atomic BECs. In real experiment, the cold molecules can be produced from a Fermi gas of atoms11,12 or an atomic BECs based on Feshbach resonance, Raman photoassociation or stimulated Raman adiabatic passage1315. For

example, a two-photon stimulated Raman transition in a 87Rb BECs has been used to produce 87Rb2 molecules in a single rotational-vibrational state16, where the input Raman laser pulse couples the molecular levels and reduces spontaneous emission. There is a nonlinear resonant transfer between atoms and molecules, as well as term proportional to the densities in the coupled atomic-molecular BECs. This type of soliton solutions have been studied in the nonlinear optics17,18 and in the problem of the self-localization of impurity atoms BECs19. The parametric solitons have been investigated in ref. 20. The coherent dynamics of this coupled atomic-molecular BECs have also been studied, which shows very rich behaviors, such as exact dark states solution21,22, crystallized and amorphous vortices23, Rabi oscillations24 and so on.

The intrinsic nonlinearity is the most remarkable characteristic of the BECs systems. In the past years, many interesting experiments, for example, the sonic-analogue of black holes, could be explored with spatial modulation of the interatomic interaction on short length scales. In refs 25 and 26, a promising technique (optical

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Feshbach resonance) is proposed to control the scattering length. With the development of this topic, a successful control of a magnetic Feshbach resonance of alkali-metal atoms was illustrated in ref. 27. In ref. 28, submicron control of the scattering length has been demonstrated by applying a pulsed optical standing wave to a BECs of ytterbium (174Yb) atoms. In recent research, the Nonlinear Schrdinger equation (NLSE) or the Gross-Pitaevskii equation (GPE) with spatially dependent cubic and quintic nonlinearities can be applied to the pulse propagation on optical ber29, photonic crystals30, and the study of BECs31,32. The wide localized soliton solutions, the wide vector solutions, the dark soliton solutions and so on have been worked out3336. The localized nonlinear waves in quasi-two-dimensional BECs with spatially modulated nonlinearity and in two-component BECs with time- and space- modulated nonlinearities are constructed3739. However, there is few work considering the two-dimensional atomic-molecular BECs with space-modulated nonlinearities.

In this paper, we investigate the nonlinear matter waves in the two-dimensional atomic-molecular Bose-Einstein condensates with space-modulated nonlinearities, which can be described by the coupled GP equations with space-modulated nonlinearities. We work out three kinds of localized nonlinear wave solutions for both the attractive spatially inhomogeneous interactions and the repulsive ones by using the similarity transformation40. Our results show that the topological properties of the localized nonlinear matter waves given by the Jacobi elliptic function can be described by the principal quantum number n and the secondary quantum number l, while the topological properties of the localized nonlinear matter waves given by the rational function are independent of the principal quantum number n, only depend on the secondary quantum number l. The Jacobi elliptic solution is linearly stable only for the principal quantum number n = 1, while the stability of the rational form solutions depends on the chemical potential and Raman detuning.

Results

The coupled Gross-Pitaevskii equation with space-modulated nonlinearity. In real experiment, the coherent free-bound stimulated Raman transition can cause atomic BECs of 87Rb to produce a molecular BECs of 87Rb2. If the molecular spontaneous emission and the light shi eect can be ignored41,42, according to the mean eld theory, the coupled atomic-molecular BECs4244 with three-body interaction term can be written as

=

i t M g g V

i t M g g V

a

a

2 2 ,

2

2 , (1)

+ + + +

+

a a am m a a ext

a

+ + + +

=

2 2 2 2 4 ( )

2 2 2 2 4 ( )

2

a a m

m

m

am a m m m m ext

m

+ +

m a m

where = +

2 2

2

2 2 i(i=a, m) denotes the macroscopic wave function of atomic condensate and molecular condensate respectively,

= =

,

x y

g g

,

a

2 2

,

g = +

am

2

, a and m represent respectively the

cubic and quintic nonlinearity the strengths of interaction, = =

+

V i a m

( , )

ext

4 4

a a

m m

aM m

a M

a

M M M M

4

/ ( )

ama m a m

2 2 2

are the trapping potentials, Ma(Mm) is the mass of atomic (molecule), is the parametric coupling coefficient which describes the conversions of atoms into molecules due to stimulated Raman transitions. The parameter characterizes Raman detuning for a two photon resonance16,41,45. Integrating along the transverse coordinates, the above equations for the wave functions i(i=a, m) in dimensionless form can be written as the coupled GP equations

=

i M x y

( ) ( )2

i

i t g g

x y

+ + + + +

+

a

( )

2

m

2 2 ,

4 ( )

2 , (2)

a a am m a a a a m

=

i t g g x y

+ + + + +

2 2 2 4 2 2 2

2 2 2 4 2 2 2

2

am a m m m m

m a m

+ +

The unit of length, time and energy correspond to

. .

m um s

/( ) ( 1 07 ), ( 1 6 10 )

1 3

and , respectively. In this paper, we use the parameters of atomic-molecular BECs of 87Rb system with Mm=2Ma=2m (m=144.421027Kg), gm=2ga (aa=101.8aB), where aB is the Bohr radius.

Now we consider the spatially localized stationary solution

= =

a m of ((2)) with i(x, y) (i = a, m) being a real function for lim|x|,|y|i(x, y) = 0. This maps (2) into the following coupled equations

2

2

x y e x y e

( , ) , ( , )

i t

i t

a a

m m

2

2

+

a a

a a am m a a a a

a a a m

m m

m m am m a m m m

a m m a

3 2 5 2 2 2

+

+ + =

x y g g x y

x y g g x y

2 2 2 ( )

2 2 2 0,

4 4 4 ( )

8 4 2 2 0, (3)

2

2

+ + + + +

+ =

2

2

3 2 5 2 2 2

2

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

where a, m are chemical potentials. In order to solve the above equations, we take the similarity transformation

= =

x y U X x y V X x y

( , ) ( ( , )), ( ( , )), (4)

a a m m

to transform (3) to the ordinary dierential equations (ODEs)

+ + + + =

+ + + + =

3

U b U b UV b UV b U

V b U V b V b U b V

0,

5 0, (5)

XX

XX

11

12

2 13 14

5

2

3

2

21

22

23

24

where bij, i=1, 2, j= 1, 2, 3 are constants. Substituting (4) into (3) and letting U(X), V(X) to satisfy (5), we obtain a set of partial dierential equations (PDEs). Solving this set of PDEs, we have

g b

F e g

= = =

= =

=

2 , 4 , 2 ,

2 , 2 ,

2 ,

3

2 ,

11

2

+ + + + + +

+ + + +

( 4 ) 22

2

( 4 ) 12

2

x xy y

m

x xy y

am

( 4 )

x xy y

a

2 2 2 2 2 2

2 2 2 2

bF e g

bF e

bF e

14

4

( 6 ) 24

4

( 6 )

x xy y

m

x xy y

a

bF e

F c KummerM y x y x

c KummerU y x y x

X erf x y

2

a

( )

2 ( )

2 ,

1

+

a

2

2

( )

2 ( ),

3

2 ,

=

1

2 ( )

2 ,

+

=

m a

1 2 , (6)

where

( ) =

erf s e d

2

2 is called an error function,

KummerM ( ) y

s

y x 2

3 2

, , ( )

a 2 and

( ) 2

x

KummerU ( ) y

0

y x 2

3 2

a 2 are solutions of the the ordinary dierential equation

+ =

F Y YF Y F Y

( ) ( ) ( ) 0, (7)

YY Y a

where Y = y x. Specially, when =

, , ( )

( ) 2

x

3

a , the KummerM function can be simplied as exponential function

e

y x

( )2

2 . In this case, the interactions become as

g bc y x e g

b

c y x e

= =

11

x xy y

m

x xy y

2 ( ) , 4 ( ) ,

2 ( ) ,

+ + + +

+ +

a

2 2

( ) 22

1

2 2

( )

1

2 2 2 2

2 2

g bc y x e

=

12

x xy y

am

2 2

( )

1

which is experimentally feasible due to the exible and precise control of the scattering lengths achievable in BECs with magnetically tuning the Feshbach resonances5,6,28.

Rational solution of the atomic-molecular BECs with three-body interaction. When the coupling = 0, (4) and (5) gives the rational formal solution of (2)

=

=

e e

xy i t

a

a

,

m

2

g g g g g c erf x y g g g

e e

3( ) ( 2 ) ( ( ) 4 ( )

, (8)

2 2

2 2

2 ( )

12 22 11 22 12

+ +

14 12 22

2

xy i t

a

where c is arbitrary constant and F is given in (6).

In order to investigate the topological properties of the exact spatially localized stationary solution (8), we plot their density distributions. In Fig.1, it can be observed that the energy packets are striped distribution, and the number of the energy stripes increases with the chemical potential a when is xed. It can also be seen that some zero points appear on the middle density stripe along line y=x when the number of the density stripes is odd.

Jacobi elliptic function solution. When the three-body eect is very weak and the coupling = 0, we have the following exact solutions of (5),

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 1. The density distributions |a|2 of the atomic-molecular BEC with three-body interaction term as the function of and a with =0.02, b11=3, b22=12, b12=1. The energy packets are striped distribution. (ac) show that the number of the energy stripes increases with chemical potential a when is xed. (df) illustrate that some zero points appear on the middle density stripe when the number of the density stripes is odd.

( ) , 2

2 ,

( ) , 2

2 (9)

=

U X c cn c X

V X c cn c X

=

1 0 0

2 0 0

or

( ) , 2

2 ,

( ) , 2

2 (10)

=

U X d sd d X

V X d sd d X

=

0 0

0 0

where c0, d0 are arbitrary constants,

= =

1 2

b bb b b b

b bb b b b

22 1211 22 12 21

11 2111 22 12 21 and cn, sd=sn/dn are Jacobi elliptic functions.

When imposing the bounded condition lim|x|,|y|i(x, y)= 0, we have =

,

2(2 1) 2

2 and =

( )

4 2

2 , where

( )

c

N K 0

d

NK 0

N is a natural number and

=

12 .

From Eqs(4), (6), (9) and (10), we obtain the Jacobi elliptic function solutions for the atomic-molecular BEC

(2)

K d

1 sin

( )

22 0

/2 1 2

2

( )

( )

n K Fe

cn nK erf x y e

2 2

xy

1

=

2 2

2

a

+

2 2 ( ) ,

22 ,

i t

a

m

n K Fe

cn nK erf x y e n N

xy

2 2

2

2

2 2

=

m

+

=

i t

2 2 ( ) ,

22 , ( 2 1), (11)

or

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

( )

( )

nK Fe

sd nK erf x y e

2 2

xy

=

1 2 2 2

a

+

2 2 ( ) ,

22 ,

i t

a

m

nK Fe

sd nK erf x y e n N

2 2

xy

=

2 2 2 2

2 2 ( ) ,

22 , ( 2 ) (12)

= .

i t

m

+

Here we discuss the existence regions of the spatially localized stationary solution (11) and (12) by assuming the two constraint conditions 1>0 and 2> 0. We have the eight cases of parameters b11, b12 and b22. According to the real experiment, we consider the following two cases:

(1) b22<b11<0 and >

b b b

2 12 11 22 .(2) b22>b11>0 and < <

b b b b

2

11 22 12 11.

These correspond to two cases of the intercomponent interaction parameters ga, gm and gam:

(a) gm>ga>0 and <

g g g

2 am a m .(b) gm<ga<0 and > >

g g g g

2

a m am a.

These are the regions that the exact spatially localized stationary solutions (11) and (12) exist. Now we only consider case (b), which denotes two self-attractive atom-atom interactions, two self-attractive molecular-molecular interactions, and attractive and repulsive atomic-molecular interactions. The other cases can be analysed in the same way.

In the following, we will see that the integer n and the number of the zero points of function F which equals to that of the KummerU and KummerM functions determine the topological properties of the atom and molecular packets, so we call n and l as the principal quantum number and the secondary quantum number, respectively. In order to investigate the topological properties of the exact spatially localized stationary solution (11) and (12), we plot their density distributions by manipulating the principal quantum number n when the secondary quantum number l is xed. In Fig.2, we analyse the atomic BEC when the secondary quantum number l is xed and the principal quantum number n is modulated. It is easy to see that the number of density packets for each quantum states is equal to 2n. And the number of density packets on each quantum states increases two by two when the principal quantum number n increases. The properties of the molecular BEC are similar to that of the atomic BEC. In Fig.3, we analyse the interactions of the atomic BEC and the molecular BEC when the secondary quantum number l is xed. It is shown that the interaction is stronger when N= 1 and becomes weaker with the increasing of N.

When the principal quantum number n is xed, we can adjust the secondary quantum number l to observe the properties of the atomic-molecular BEC. Figure4 demonstrates the density distributions of atomic-molecular BEC for dierent secondary quantum number l. It is easy to nd that the number of energy packets increases when l increases, and the number of the nodes for each quantum state equals to the secondary quantum number l. And some zero points appear on the middle density packets along the line y = x when the number of the secondary quantum number l is even. Figure5 demonstrates the interaction of the atomic BEC and molecular BEC when the principal quantum number n is xed. It is shown that the number of the atomic-molecular pair is the function of the secondary quantum number l and some zero points appear on the middle atomic-molecular pair along the line y=x when the number of the secondary quantum number l is even.

Now we analyse the eect of Raman detuning for the atomic-molecular BEC. From Fig.6, we can see that when < a and is xed, the number of the density packets increases one by one with the increasing of the chemical potential a. When a, there is only one density packets for each quantum states. The absolute of a aect the shape of the energy packets: when the absolute of a is small, the shape of the density packet is like circle, and when the absolute of a is larger, the shape of the density packet becomes narrow and long.

Rational formal solution. When the three-body eect is very weak and the coupling 0, (4) and (5) also gives the rational formal solution of (2)

=

13 + +

=

13 + +

= .

12 2

2 ( ) 27

,

b F

xy i t

a

a

( ) ( )

b c erf x y b

e e

2 2

2

2

11

b F

m

24

2 ( ) 27

,

13

xy i t

2

a

b c erf x y b

e e

2 2

2

2

11

28 (13)

bF e

13 ( 3 )

+ +

x xy y

13

2 2

where b13=b23, b11>0, c is arbitrary constant and F is given in (6).

In order to investigate the topological properties of the exact spatially localized stationary solution (13), we plot their density distributions by adjusting the secondary quantum number l. The secondary quantum number l

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

Figure 2. The density distributions |a|2 of the atomic BEC as the function of the principal quantum number n when the secondary quantum number is xed. The wave function a take the form in Eqs(11) and (12) with b11=3, b22=12, b12=1 and = 0.2. The number of density packets for each quantum states is equal to 2n. (ac) show the density distributions of the odd parity wave functions (11) for n= 1, 3, 5, respectively. (df) illustrate the density distributions of the even parity wave functions (12) for n= 2, 4, 6, respectively. The solution displayed in guer (a) is linear stable. The unit of the length is 1.07m.

is always zero for <a, and can be taken dierent values for a. In Fig.7, it can be observed that the energy packets are striped distribution, and the number of the energy stripes increases with the chemical potential a when the secondary quantum number l = 0 and is xed. When the secondary quantum number l 0, there is only one energy stripe and the energy stripe becomes more narrower with the increasing of the secondary quantum number l. It can also be seen that some zero points appear on the middle density stripe along line y=x when the number of the density stripes is odd. Figures1 and 7 show that the rational solution (8) and (13) have similar topological properties, which implies that three-body interaction doesnt hinder the formation of the localized nonlinear matter wave solutions.

Linear stability analysis. In the following, we analyse the linear stability of the solutions (11), (12) and (13) by using the linear stability analysis. A perturbed solution is constructed as46,47

i t i t i t

= + +

= + +

[ ( , ) ( , ) ( , ) ] ,

[ ( , ) ( , ) ( , ) ] (14)

a a

x y u x y e w x y e e

x y u x y e w x y e e

1 1

2 2

a

i t i t i t

m m

m

where |u1|1, |u2|1, |w1|1, |w2| 1 are small perturbation. Substituting this perturbed solution into (2) and neglecting the higher-order terms in u1, u2, w1 and w2, we obtain the eigenvalue problem

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

Figure 3. The density distributions |a|2+|m|2 of the atomic-molecular pair as the function of N when the secondary quantum number l=1. Here b11=3, b22=12, b12=1, and =0.2. (a) shows the interaction of the atomic-molecular pair for N=1. (b) shows the interaction of the middle atomic-molecular pair for N=2, andit also displays that the interactions are weaker than the interaction in (a). (af) illustrate that the interaction becomes weaker with the increasing of N. The unit of the length is 1.07m.

2 2

2 2

2

2

L g g g

g L g g

g g L g

g g g L

u w u w

u w u w

m a a a am a m am a m

a a m am a m am a m m

a am a m am a m m m

am a m am a m a m m

2

1

1

2

2

2

2

2

1

1

2

2

=

1

1

2

2

,

(15)

where

1

1 2 ( ) 2

2 2 2 2 2 2 2

L g g x y

L g g x y

= + + +

= + + + .

1 2 ( ) ,

14 ( ) 2

x y a a am m a

x y m m am a a

2 2 2 2 2 2 2

2

14 ( ) 2

Numerical experiments show that the eigenvalue of the eigenvalue problem (15) is real for n=1. This suggests that the localized nonlinear matter wave solution (11) is linearly stable for n= 1 and solution (12) is unstable. For the solution (13), it can be shown that the linear stability rests on the chemical potential a and the

Raman detuning (see Fig.8).

Discussion

In this paper, we focus on the analytic solutions of atomic-molecular BECs and the eects of the coupling and the Raman detuning on the atomic-molecular BECs. The system in this report is like the one in the ref. 43. Comparing to the atomic-molecular system given in the ref. 43, Gupta and Dastidar have proposed a more complicated model when they study the dynamics of atomic and molecular BECs of 87Rb in a spherically symmetric trap coupled by stimulated Raman photoassociation process in the ref. 42. In fact, the light shi eect in Gupta and Dastidars model almost has the same function as the Raman detuning term. So, it can be contributed to the Raman detuning term. Based on this reason, we dont consider the light shi eect and take the form of the atomic-molecular BECs system as the form in ref. 43.

In the ref. 43, they show that the coherent coupling between atoms and molecules changes the situation crucially and it is sensitive to the presence of vortices. For example, when the coupling is zero, each of the atoms and molecular BECs wave function forms an independent triangular vortex lattice, and a nonzero coupling proposes more dramatic changes. Our results show that the coupling can change the topological structure of

SCIENTIFIC REPORTS

7

www.nature.com/scientificreports/

Figure 4. The density distributions |a|2 of the atomic-molecular BEC as the function of the secondary quantum number l. Here =0.2, b11=3, b22=12, b12= 1. The number of the nodes for each quantum state equals to the secondary quantum number l and some zero points appear on the middle one when the number of the density packets is odd. (a1a4) show the density distributions of the odd parity wave functions (11) for n=1 and l= 1, 2, 3, 4, respectively. (b1b4) show the density distributions of the even parity wave functions (12) for n=2 and l= 1, 2, 3, 4, respectively. The solutions displayed in the rst and third gures on the upper row are linear stable. The unit of the length is 1.07m.

Figure 5. The density distributions |a|2+|m|2 of the atomic-molecular pair as the function of the secondary quantum number l with =0.2, b11=3, b22=12, b12=1. The number of the atomic-molecular pairs equals to l+1. (ac) show the density proles of the atomic-molecular pair for l is odd. (df) show the density proles of the atomic-molecular pair for l is even, and it also displays that some zero points appear on the middle one. The unit of the length is 1.07m.

the localized nonlinear wave of the atomic- molecular BECs. In the case of = 0, Figs14 illustrate that the topological structures depend on the principal quantum number n and the secondary quantum number l, and each density packet is like a circle and oval. When 0, Fig.6 display the density packets are striped distribution and their topological structures only reply on the secondary quantum number l and are independent on the principal quantum n.

In real experiment, spatial modulation of the interatomic interaction can be achieved. In the recent experiment28, the authors apply a pulsed optical standing wave to a BEC of ytterbium (174Yb) atoms and realize the submicron control of the scattering length. The experimental phenomena is well explained by the semi-classical

SCIENTIFIC REPORTS

8

www.nature.com/scientificreports/

Figure 6. The eect of Raman detuning for the atomic-molecular BEC. (ac) show that the number of the density packets increases with the chemical potential a when <a. (a,d) reveal that the number of the density packets dont depends on the chemical potential a and there is only one density packet for each quantum state when a, it also show the value of a eects the shape of the density packet. The solution displayed in gure (b) is linear stable. The unit of the length is 1.07m.

Figure 7. The density distributions |a|2 of the atomic-molecular BEC as the function of and a with =0.02, b11=3, b22=12, b12=1. The energy packets are striped distribution. (a1b3) show that the number of the energy stripes increases with chemical potential a when the secondary quantum number l=0. (c1c3) show that there is only one density stripe when the secondary quantum number l0. (b1,c3) illustrate that some zero points appear on the middle density stripe when the number of the density stripes is odd. The solutions displayed in gures (a1,b1) are linear stable.

theory of Bohn and Julienne47. In this paper, the interaction ga, gm, gam and the coherent coupling all depend on the spatial variables. Under that conditions, the stable exact solutions can be worked out for the rst time. The spatial modulation of the interaction can be realized by the above experiment, but there is no successful experiment for the spatial modulation of the coherent coupling. We hope that our research will stimulate the further research on the spatial modulation of the atomic- molecular BECs.

SCIENTIFIC REPORTS

9

www.nature.com/scientificreports/

Figure 8. Linear stability. Eigenvalue for dierent principal quantum numbers n with parameters b11=3,

b22=12, b12=1. (a1a3) show that the exact solution (11) is linearly stable only for n=1; (b1b3) show that the exact solution (12) is linearly unstable for all n; (c1c3) illustrate that the solution (13) are linearly stable in the two group parameters =0.01, m=0.1 and =0.1, m=0.2.

It is obvious that the Raman detuning term in the atomic- molecular BECs behaves just like the chemical potential to control the systems energy. In this paper, the results imply that a not only changes the altitude of the wave packets, but also changes the topological structures of the nonlinear waves. When a 0, the number of the energy packets changes with the chemical potential a. When a< 0, there is only one energy packet for each quantum state.

In summary, we have worked out three kinds of localized nonlinear matter wave solutions of the two-dimensional atomic-molecular BECs with space-modulated nonlinearity and considered the effect of three-body interaction on the localized nonlinear matter wave solutions. Our results show that the matter wave functions given by elliptic function have even parity for the even principal quantum number and odd parity for the odd one, the number of density packets for each quantum state is twice of the principal quantum number n, and the number of density packets increases two by two with the principal quantum number n. The number of the nodes equals to the secondary quantum number l. For the nonlinear matter wave given by rational function, the number of the energy stripes increases with the chemical potential a when the secondary quantum number l=0 and is xed. When the secondary quantum number l 0, there is only one energy stripe for each quantum state and the energy stripe becomes more narrower with the increasing of the secondary quantum number l. Odd (even) secondary quantum number l leads to even (odd) number of the energy packets (stripes). Some zero

SCIENTIFIC REPORTS

10

www.nature.com/scientificreports/

points appear on the middle energy packets (stripes) along line y=x for even secondary quantum number l. We also analyse the eect of Raman detuning for the atomic-molecular BECs. The value of a can change the number and shape of the energy packets (stripes). The stability of our solutions is analysed: the nonlinear matter wave solution (11) is linearly stable for the principal quantum number n= 1, the solution (12) is unstable, and the stability of the solution (13) rests on the chemical potential a and the Raman detuning . Our results are signi-cant to matter wave management in high-dimensional atomic-molecular BECs.

Methods

We use the coupled Gross-Pitaevskii equation to describe the atomic-molecular BECs. Taking into account the term responsible for the creation of molecules48, the Hamiltonian is taken as

g h m

a

2 2 4 2

2 [ ] (16)

[circumflexnosp3char][circumflexnosp3char]

[circumflexnosp4char][circumflexnosp4char]

[circumflexnosp3char][circumflexnosp3char] [circumflexnosp3char]

am a a m m m a a m a a

First, the coupled Gross-Pitaevskii equation is decomposed into two ODEs and a number of PDEs making use of the similarity transformation. Then we solve these ODEs and PDEs by using some solving techniques and some special functions, such as error function, KummerU function and Jacobi elliptic function. The nal interaction parameters are altered to = =

g g g g

2 ,

m a am a

13 and the chemical potential satises m=2a.

1. Regal, C. A., Ticknor, C., Bohn, J. L. & Jin, D. S. Tuning p-Wave Interactions in an Ultracold Fermi Gas of Atoms. Nature (London) 424, 47 (2003).

2. Strecher, K. E., Partridge, G. B. & Hulet, R. G. Conversion of an Atomic Fermi Gas to a Long-Lived Molecular Bose Gas. Phys. Rev. Lett. 91, 080406 (2003).

3. Xu, K. et al. Formation of Quantum-Degenerate Sodium Molecules. Phys. Rev. Lett. 91, 210402 (2003).4. Hodby, E. et al. Production Efficiency of Ultracold Feshbach Molecules in Bosonic and Fermionic Systems. Phys. Rev. Lett. 94, 120402 (2005).

5. Roberts, J. L. et al. Resonant magnetic el control of elastic scattering in cold 85Rb. Phys. Rev. Lett. 81, 5109 (1998).6. Donley, E. A., Claussen, N. R., Thompson, S. T. & Wieman, C. E. Atom-molecule coherence in a Bose-Einstein condensate. Nature (London) 417, 529 (2002).

7. Wynar, R. et al. Molecules in a Bose-Einstein Condensate. Science 287, 1016 (2000).8. Jones, K. M., Tiesinga, E., Lett, P. D. & Julienne, P. S. Ultracold photoassociation spectroscopy: Long-range molecules and atomic scattering. Rev. Mod. Phys. 78, 483 (2006).

9. Ling, H. Y., Pu, H. & Seaman, B. Creating a Stable Molecular Condensate Using a Generalized Raman Adiabatic Passage Scheme. Phys. Rev. Lett. 93, 250403 (2004).

10. Ling, H. Y., Maenner, P., Zhang, W. & Pu, H. Adiabatic theorem for a condensate system in an atom-molecule dark state. Phys. Rev. A 75, 033615 (2007).

11. Greiner, M., Regal, C. A. & Jin, D. S. Emergence of a molecular BoseCEinstein condensate from a Fermi gas. Nature (London) 426, 537 (2003).

12. Cubizolles, J. et al. Production of Long-Lived Ultracold Li2 Molecules from a Fermi Gas. Phys. Rev. Lett. 91, 2401401 (2003).13. Timmermans, E., Tommasini, P., Hussein, M. & Kerman, A. Feshbach resonances in atomic Bose-Einstein condensates. Phys. Rep. 315, 199 (1999).

14. Dunie, R. A. & Stoof, H. T. C. AtomCmolecule coherence in Bose gases. Phys. Rep. 396, 115 (2004).15. Volz, T. et al. Nat.Preparation of a quantum state with one molecule at each site of an optical lattice. Nat. Phys. 2, 692 (2006).16. Wynar, R. et al. Molecules in a Bose-Einstein Condensate. Science 287, 1016 (2000).17. Karpierz, M. A. Coupled solitons in waveguides with second- and third-order nonlinearities. Opt. Lett. 16, 1677 (1995).18. He, H., Werner, M. J. & Drummond, P. D. Simultaneous solitary-wave solutions in a nonlinear parametric waveguide. Phys. Rev. E 54, 896 (1996).

19. Sacha, K. & Tummermans, E. Self-localized impurities embedded in a one-dimensional Bose-Einstein condensate and their quantum uctuations. Phys. Rev. A 73, 063604 (2006).

20. Vaughan, T. G., Kheruntsyan, K. V. & Drummond, P. D. Three-dimensional solitons in coupled atomic-molecular Bose-Einstein condensates. Phys. Rev. A 70, 063611 (2004).

21. Winkler, K. et al. Atom-Molecule Dark States in a Bose-Einstein Condensate. Phys. Rev. Lett. 95, 063202 (2005).22. Dumlke, R. et al. Sub-natural-linewidth quantum interference features observed in photoassociation of a thermal gas. Phys. Rev. A 72, 041801(R) (2005).

23. Liu, C. F., Fan, H., Gou, S. C. & Liu, W. M. Crystallized and amorphous vortices in rotating atomic-molecular Bose-Einstein condensates. Scientic Report 4, 4224 (2014).

24. Ishkhanyan, A., Chemikov, G. P. & Nakamura, H. Rabi dynamics of coupled atomic and molecular Bose-Einstein condensates. Phys. Rev. A 70, 053611 (2004).

25. Fedichev, P. O., Kagan, Yu., Shlyapnikov, G. V. & Walraven, J. T. M. Inuence of Nearly Resonant Light on the Scattering Length in Low-Temperature Atomic Gases. Phys. Rev. Lett. 77, 2913 (1996).

26. Theis, M. et al. Tuning the Scattering Length with an Optically Induced Feshbach Resonance. Phys. Rev. Lett. 93, 123001 (2004).27. Dominik, M. B. et al. Control of a magnetic Feshbach resonance with laser light. Nat. Phys. 5, 339 (2009).28. Yamazaki, R., Taie, S., Sugawa, S. & Takahashi, Y. Submicron Spatial Modulation of an Interatomic Interaction in a Bose-Einstein Condensate. Phys. Rev. Lett. 105, 050405 (2010).

29. Mihalache, D., Mazilu, D., Bertolotti, M. & Sibilia, C. Exact solution for nonlinear thin-lm guided waves in higher-order nonlinear media. J. Opt. Soc. Am. B 5, 565 (1998).

30. Kivshar, Yu. S. & Agrawal, G. P. Optical solitons: From Fibers to Photonic Crystals (Academic Press, Boston, 2003).31. Pitaevskii, L. & Stringari, S. Bose-Einstein condensation (Clarendon Press, New York, 2003).32. Zhou, X. F., Zhou, J. & Wu, C. Vortex structures of rotating spin-orbit-coupled Bose-Einstein condensates. Phys. Rev. A. 84, 063624 (2011).

33. Arroyo Meza, L. E., de Souza Dutra, A. & Hott, M. B. Wide vector solitons in systems with time- and space-modulated nonlinearities. Phys. Rev. E. 88, 053202 (2013).

H d r hm

3 2 2 2 2

[circumflexnosp3char][circumflexnosp3char]

a a a m

[circumflexnosp4char][circumflexnosp4char]

g

a

=

+

[circumflexnosp3char][circumflexnosp3char][notdef] +

+ +

m

[circumflexnosp3char][circumflexnosp3char]

m m m

+ + +

.

g

References

SCIENTIFIC REPORTS

11

www.nature.com/scientificreports/

34. Belmonte-Beitia, J. & Cuevas, J. Symmetric and asymmetric bound states for the nonlinear Schrodinger equation with inhomogeneous nonlinearity. J. Phys. A: Math. Theor. 42, 165201 (2009).

35. Kengne, E., Vaillancour, R. & Malomed, B. A. Bose-Einstein condensates in optical lattices: the cubic-quintic nonlinear Schr? dinger equation with a periodic potential. J. Phys. B: At. Mol. Opt. Phys. 41, 205202 (2008).

36. Kengne, E. & Lakhssassi, A. Dynamics of solitary-wave structures in one-dimensional Gross-Pitaevskii equation with distributed coefficients. The European Physical Journal Plus 130, 197 (2015).

37. Wang, D. S, Hu, X. H. & Liu, W. M. Localized nonlinear matter waves in two-component Bose-Einstein condensates with time- and space-modulated nonlinearities. Phys. Rev. A. 82, 023612 (2010).

38. Wang, D. S., Hu, X. H., Hu, J. P. & Liu, W. M. Quantized quasi-two-dimensional Bose-Einstein condensates with spatially modulated nonlinearity. Phys. Rev. A. 81, 025604 (2010).

39. Sakaguchi, H. & Malomed, B. A. Two-dimensional solitons in the Gross-Pitaevskii equation with spatially modulated nonlinearity. Phys. Rev. E. 73, 026601 (2006).

40. Avelar, A. T., Bazeia, D. & Cardoso, W. B. Solitons with cubic and quintic nonlinearities modulated in space and time. Phys. Rev. E. 79, 025602(R) (2009); Juan Belmonte-Beitia, Vctor M. Prez-Garca, Vadym Vekslerchik, and Pedro J. Torres, Lie Symmetries and Solitons in Nonlinear Systems with Spatially Inhomogeneous Nonlinearities. Phys. Rev. Lett. 98, 064102 (2007).

41. McKenzie, C. et al. Photoassociation of Sodium in a Bose-Einstein Condensate. Phys. Rev. Lett. 88, 120403 (2002).42. Gupta, M. & Dastidar, K. R. Control of the dynamics of coupled atomic-molecular Bose-Einstein condensates: Modied Gross-Pitaevskii approach. Phys. Rev. A 80, 043618 (2009).

43. Woo, S. J., Park, Q. H. & Bigelow, N. P. Phases of Atom-Molecule Vortex Matter. Phys. Rev. Lett. 100, 120403 (2008).44. Tikhonenkov, I. & Vardi, A. Atom-molecule dephasing in an SU(1,1) interferometer based on the stimulated dissociation of a molecular Bose-Einstein condensate. Phys. Rev. A 80, 051604(R) (2009).

45. Heinzen, D. J., Wynar, R., Drummond, P. D. & Kheruntsyan, K. V. Superchemistry: Dynamics of Coupled Atomic and Molecular Bose-Einstein Condensates. Phys. Rev. Lett. 84, 5029 (2000).

46. Yang, J. K. Nonlinear Waves in Integrable and Nonintegrable Systems (Society for Industrial and Applied Mathematics, 2010).47. Afanasjev, V. V., Malomed, B. A. & Chu, P. L. Stability of bound states of pulses in the Ginzburg-Landau equations. Phys. Rev. E 56, 6020 (1997).

48. Bradley, A. S., Gardiner, C. W. & Davis, M. J. Bose-Einstein condensation from a rotating thermal cloud: Vortex nucleation and lattice formation. Phys. Rev. A 77, 033616 (2008).

Acknowledgements

Y.Q.Y. is supported by the NSFC under grant No. 11301179. W.M.L. is supported the NSFC under grants Nos 11434015, 61227902, 61378017, 11301179, 11271362 and 11375030, NKBRSFC under grants Nos 2012CB821305, SKLQOQOD under grants No. KF201403, SPRPCAS under grants Nos XDB01020300 and XDB21030300, Beijing Nova program No. Z131109000413029 and Beijing Finance Funds of Natural Science Program for Excellent Talents No. 2014000026833ZK19, W. Han is supported by the NSFC under grant No. 11547126.

Author Contributions

W.-M.L. conceived the idea and supervised the overall research. Y.-Q.Y., W.H. and. D.-S.W. performed the computations and writing the program for pictures. J.L. performed the analyse of linear stability. Y.-Q.Y. wrote the paper with helps from all other co-authors.

Additional Information

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Yao, Y.-Q. et al. Localized spatially nonlinear matter waves in atomic-molecular Bose-Einstein condensates with space-modulated nonlinearity. Sci. Rep. 6, 29566; doi: 10.1038/srep29566 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

12