Open access peer-reviewed chapter - ONLINE FIRST

Magnetic Solitons in Optical Lattice

By Xing-Dong Zhao

Submitted: November 2nd 2018Reviewed: February 27th 2019Published: April 16th 2019

DOI: 10.5772/intechopen.85492

Downloaded: 82

Abstract

In this chapter, we discuss the magnetic solitons achieved in atomic spinor Bose-Einstein condensates (BECs) confined within optical lattice. Spinor BECs at each lattice site behave like spin magnets and can interact with each other through the static magnetic dipole-dipole interaction (MDDI), due to which the magnetic soliton may exist in blue-detuned optical lattice. By imposing an external laser field into the lattice or loading atoms in a red-detuned optical lattice, the light-induced dipole-dipole interaction (LDDI) can produce new magnetic solitons. The long-range couplings induced by the MDDI and ODDI play a dominant role in the spin dynamics in an optical lattice. Compared with spin chain in solid material, the nearest-neighbor approximation, next-nearest-neighbor approximation, and long-range case are discussed, respectively.

Keywords

  • spinor Bose-Einstein condensates
  • spin wave
  • optical lattice
  • magnetic soliton

1. Introduction

Soliton can be classified into different species, such as matter-wave soliton, magnetic soliton, optical soliton, and so on [1, 2]. The magnetic solitons, which describe the localized magnetization, are very important excitations in the Heisenberg spin chain in solid system in condensed matter [3]. Because of the defect and impurity in solid material, it is difficult to perform experimental observation and manipulation in these systems. In addition, magnetic solitons originate from the Heisenberg-like short-range exchange interaction between electrons in solid state systems, so the theoretical models and treatment are limited to only the approximation of nearest-neighbor interaction, it is disadvantageous to the study of long-range interaction-induced dynamics.

On the other hand, the spinor ultracold atoms in an optical lattice offer a pure and well-controlled platform to study spin dynamics. If the lattice trapping atoms is deep enough, the spinor atoms undergo a ferromagnetic-like phase transition that leads to a macroscopic magnetization of the condensates array; then, the individual sites become independent with each other. Spinor BECs at each lattice site behave like spin magnets and can interact with each other through the static MDDI, which can cause the ferromagnetic phase transition and the spin wave excitation. It is atomic spin chain in optical lattice [4, 5, 6, 7]. By tuning the system parameters, magnetic soliton can be produced; however a key difference of the atomic spin chain from solid-state one is that we are facing a completely new spin system in which the long-range nonlinear interactions play a dominant role. To demonstrate the high controllability, an external laser can be imposed into the lattice and induce light-induced dipole-dipole interaction (LDDI) to produce new magnetic solitons. In fact, the LDDI can be induced by loading the atoms in a red-detuned optical lattice as used in [5]. However, the method used here has good feasibility, and some novel physical phenomena can be observed.

In order to make comparison with the Heisenberg spin chain, we used the nearest-neighbor approximation and the next-nearest-neighbor approximation in our discussion. In addition, the two approximations are effective and reasonable with some appropriate parameters chosen in experiments.

In this chapter, we are going to be seeing how magnetic soliton in atomic spin chain is, which will be theoretically described. We first discuss the formation of atomic spin chain in optical lattice in Section 2. Then we shift our focus to the condition for magnetic soliton in Section 3. Section 4 presents the conclusion and outlook.

2. Atomic spin chain

2.1 Atomic spin chain with magnetic dipole: dipole interaction

The spinor BECs can be trapped within one-, two-, and three-dimensional optical lattices with different experimental techniques; however, we only take one-dimensional case as an example. As depicted in Figure 1, the atomic gases with F=1are trapped in a one-dimensional optical lattice, which is formed by two π-polarized laser beams counter-propagating along the y-axis. The light-induced lattice potential takes the form:

VLr=ULexpr2/WL2cos2kLyE1

Figure 1.

Schematic diagram of spinor BECs confined in a one-dimensional optical lattice. A strong static magnetic field B is applied to polarize the ground-state spin orientations of the atomic chain along the quantized z-axis.

where r=x2+z2,kL=2π/λLis the wave number, WLthe width of the lattice beams, and the potential depth UL=Ω2/6Δwith Ωbeing the Rabi frequency and Δbeing laser detuning.

It has been theoretically and experimentally demonstrated that the condensate trapped in optical lattice undergoes a superfluid-Mott-insulator transition as the depth of the lattice wells is increased. If the lattice potential is deep enough, the condensate is divided into a set of separated small condensates equally located at each lattice site, which forms the atomic chain. The condensate localized in each lattice site is approximately in its electronic ground state, while the two lattice laser beams are detuned far from atomic resonance frequency, and the laser detuning is defined as Δ=ωLωa, ωLbeing laser frequency and ωaatomic resonant frequency. We classify the optical lattice into two categories: red-detuned optical lattice (Δ<0) and blue-detuned one Δ>0. For the former, the condensed atoms are trapped at the standing wave nodes where the laser intensity is approximately zero; the light-induced interactions are neglected in this case. In order to realize ferromagnetic phase transition, a strong static magnetic field Bis applied to polarize the ground-state spin orientations of the atomic chain along the quantized z-axis. The Hamiltonian including the long-range MDDI between different lattice sites takes the following form [7, 8, 9]:

HM=mdrψ̂m22M2+VLψ̂m+m,m,n,ndrdrψ̂mψ̂mψ̂nψ̂nVmmnncolrr+12drdrψ̂mψ̂mψ̂nψ̂nV1rrF̂mnF̂mn12drdrψ̂mψ̂mψ̂nψ̂nV2rr×F̂mnrrF̂mnrr+m,ndrψ̂mμBmnψ̂nE2

where ψ̂mrtm=±10are the hyperfine ground-state atomic field operators. The total angular momentum operator F̂for the hyperfine spin of an atom has components represented by 3×3matrices in the f=1mf=msubspace. The final term in Eq. (2) represents the effect of polarizing magnetic field, and μis magnetic dipole moment. The interaction potential includes the two-body ground-state collisions and the magnetic dipole-dipole interactions, which can be described by the potentials:

Vmmnncolrr=λsδmnδmn+λaF̂mnF̂mnδrr,E3
V1rr=μ0γB24π1rr3,E4
V2rr=3μ0γB24π1rr5.E5

In collision interaction potential in Eq. (3), λsand λaare related to the s-wave scattering length for the spin-dependent interatomic collisions. In the magnetic dipole-dipole interaction potential Eqs. (4) and (5), μ0is the vacuum permeability; the parameter γB=μBgFis the gyromagnetic ratio with μBbeing the Bohr magneton and gFthe Landé gfactor.

Under tight-binding condition and single-mode approximation, the spatial wave function ϕmrof the m-th condensate is then determined by the Gross-Pitaevskii (GP) equation [5]:

22M2+Vmr+λsNm1ϕmr2ϕmr=μmϕmr,E6

with

Vmr=ULexprL2/W2cos2kLy+,0<y<λL/2,E7

and Nmis the number of condensed atoms at the m-th site; μmis the chemical potential. In this case, the individual condensate confined in each lattice behaves as spin magnets in the presence of external magnetic fields; such spin magnets form a 1D coherent spin chain with an effective Hamiltonian [5, 6]:

Hm=mλaŜm2γBŜmBnmJmnmagŜmŜn,E8

where we have defined collective spin operators Ŝmwith components Ŝmxyzand

λa=12λadrϕmr4.E9

The parameter Jmnmagdescribes the strength of the site-to-site spin coupling induced by the static MDDI. It takes the form:

Jmnmag=μ0γB216π2drdrR1r5ϕmr2ϕnrr2,E10

where μ0is the vacuum permeability, R1=r23y2. It will be seen from this that the MDDI is determined by many parameters; however, we need to choose a parameter which is easy controllable in the experiment. The realistic interaction strengths of the MDDI and LDDI coupling coefficients have been calculated in ref. [4], in which we find the dependency of the transverse width.

Here, we only consider the ferromagnetic condensates where the spins of atoms at each lattice site align up along the direction of the applied magnetic field (the quantized z-axis) in the ground state of the ferromagnetic spin chain. Because of the Bose-enhancement effect, the spins are very large in these systems. As a result, at low temperature the spins can be treated approximately as classical vectors, and we can replace Ŝnwith its average value Sn[1]. From the Hamiltonian (8), we can derive the Heisenberg equation of motion for operator Sn±=Snx±iSnyi2=1(normalized by ):

idSn+dt=γBBzSn++jn2JnjmagSjzSn+SnzSj+.E11

Equation (11) determines the existence and the propagation of spin waves. From this equation, we can obtain the existence and dynamical evolving of the nonlinear magnetic soliton. In Sec.3, we attempt to show the long-range and controllable characteristics of MDDI and find how this long-range nonlinear dipole-dipole interaction affects the generation of the magnetic soliton. If the static magnetic field Bis strong enough, we have:

Snz=S2Sn+Sn=S2Sn+Sn+.E12

2.2 Atomic spin chain with light-induced dipole-dipole interaction

In fact, the site-to-site interaction can also be tuned by laser field, and the exchange interaction induced by the laser fields also plays an important role in atomic spin chain in optical lattice. We consider that an external modulational laser field is imposed into the system along the y-axis. This external laser field is properly chosen so that a new dipole-dipole interaction will be introduced into the system, i.e., the so-called LDDI. Its Hamiltonian takes the form [10]:

Hopt=mnmJmnoptŜmŜn++Sm+Sn.E13

Here, the coefficient Jmnoptdescribes the site-to-site spin coupling induced by the LDDI, which has the concrete form:

Jmnopt=Q144kL3drdrfcrexpR2WL2×coskLycoskLyye+1Wre1×ϕmr2ϕnrr2.E14

In above integral we define Q=γΩ12/Δ12to describe the magnitude of the LDDI; here γis the spontaneous emission rate of the atoms, and Ω1is the Rabi frequency; Δ1is the optical detuning for the induced atomic transition. A cutoff function fcr=expr/Lcis introduced to describe the effective interaction range of the LDDI, Lcbeing the coherence length. R2=r2+rr2with the transverse coordinate r=x2+z2. The unit vectors have the form:

e0=ez,E15
e±1=ex±iey,E16

and eαα=xyzwhich are the usual geocentric Cartesian frame. The tensor Wrdescribes the spatial profile of the LDDI and has the form:

Wr=3413cos2θsinξξ2+cosξξ3sin2θcosξξ,E17

where we define ξ=kLrr, and θis the angle between the dipole moment and the relative coordinate rr. Being different from MDDI in Eq. (8), it shows that the LDDI in this system leads to the spin coupling only in the transverse direction, i.e., the direction perpendicular to the external magnetic field, as shown in Eq. (13). The total Hamiltonian of the system can be written as [10].

H=Hm+Hopt.E18

Function (11) now has the form:

idϕndt=γBBzϕn+j=n±12SJnjmag1ϕj2ϕn+j=n±14SJnj1ϕn2ϕj.E19

The interaction coefficient Jijsatisfies Jij=Jijopt+12Jijmag.

Figure 2 shows the two spin coupling coefficients varying with the lattice sites and the transverse width of the condensate. It is clear that the spin coupling coefficients are sensitive to the variation of the transverse width Wof the condensate and exhibit the similar long-range behavior.

Figure 2.

The spin coupling coefficients Jmnz=Jmnmag and Jij as a function of the lattice site index for different transverse widths of the condensate in the MDDI and LDDI dominating optical lattice [4].

3. Magnetic soliton in atomic spin chain

3.1 Magnetic soliton under the MDDI

To confirm that the magnetic soliton can be formed in this system, we first consider the case that other interactions between spins are absent in this section. Namely if the solitons indeed exist, they can be generated by the MDDI alone Jij=Jijmag/2. In the usual way, we set:

Sn±=Sn±/S,E20
Snz=Snz/S,E21

here Sis the magnitude of spin for site nand only depends on the number of atoms trapped in each site. Under the external magnetic field Bconsidered here, we make an assumption that the spin deviation from the quantized z-axis is considerable, which is reasonable when the number of atoms trapped in each site is large; it is achievable according to present experiments. We also assume that the magnitude of the spin deviation snwill extend over a large number of lattice sites and vary slowly with site index n. Considering the spatial symmetry, we introduce a staggered variable:

ϕn=1nsn+,E22

where ϕnis complex; to find the existence conditions of the magnetic soliton modes, we need firstly to get the nonlinear Schrödinger equation (NLSE) for ϕnin continuum form, which is expected to be qualitatively correct for solitons in discrete limit.

3.1.1 Nearest-neighbor approximation with the MDDI

Here we firstly consider only the nearest-neighbor (NN) interactions. In other words, we assume that the on-site spin couples only to the spins at its two neighboring sites, and the longer range couplings are artificially cut off. In above equation, for small retortion from the ground state, we have:

1ϕj2112ϕj2.E23

Using the slow varying envelope approximation or long-wavelength limit, we have:

d2ϕytdy2y=ynϕn+1+ϕn12ϕna2.E24

where ais the lattice constant. Then we can obtain the continuum limit NLSE for ϕyt, namely

idt=2SJ01maga2d2ϕdy2+βϕ4SJ01magϕ2ϕ,E25

where a=λL/2is the lattice spacing between two adjacent sites, with λLbeing the wave length of the lattice laser beams, and the coefficient β=γBBz+8SJ01mag. This NLSE is integrable, and it has bright soliton solution if the effects of dispersion and nonlinearity can cancel each other. In other words, the condition for Eq. (25) to have one bright soliton solution is that the coefficient of the second-order spatial derivative term and the coefficient of the cubic nonlinear term are of the same sign [11]. The criterion for the existence of solitons requires:

fW=8S2a2J01mag2>0.E26

Since the spin coupling parameter J01magis positive, this condition cannot be satisfied, that is to say, there is no one soliton solution in atomic spin chain under the NN approximation. It drives us to in-depth consideration; we should consider more complex interactions, i.e., more site-to-site interactions would be considered.

3.1.2 Next-nearest-neighbor approximation with the MDDI

Even though the long-range effect of the dipole-dipole interaction is obvious, in this section, we consider a case slightly beyond the NN approximation, i.e., the next-nearest-neighbor (NNN) approximation, which is used in some solid-state systems. Under this approximation, each on-site spin is coupled to the four near-neighbor spins, two on its right side and the other two on its left side, respectively. In this case, we set:

αj=J0jmag/J01,E27
A2=J02/J01,E28
η=γBBz+4SJ01mag+2J01+J02mag2J02.E29

Additionally, the second-order spatial derivative is introduced:

d2ϕytdy2y=ynϕn+2+ϕn22ϕn2a2.E30

Then, we can easily get:

iytdt=4SJ0114A2a2d2ϕytdy2+ηϕyt2SJ01α1+A2α2+22A2ϕyt2ϕyt,E31

where

η=γBBz+4SJ01mag+2J01+J02mag2J02.E32

Because the LDDI is absent in this case, we have:

J01opt=0,E33
J01=12J01mag,E34
α1=α2=2.E35

Then, Eq. (31) can be rewritten as [10]:

iytdt=2SJ01mag14A2a2d2ϕytdy2+ηϕyt4SJ01magϕyt2ϕyt,E36

where η=γBBz+8SJ01mag. For two adjacent sites, A2=J02mag/J01magmeasures the relative strength of the NNN spin coupling to the NN spin coupling. The condition for Eq. (36) to have one soliton solution is

fW=8S2a2J01mag214A2>0;E37

we have A2>0.25. From the definition of A2, we can rewrite the condition as 0.25<A2<1, which agrees with the instability condition of the extended nonlinear band edge spin wave modes in the linear stability analysis in Ref. [4]. It means that one soliton solution can exist in atomic spin chain after the NNN interactions between spins are considered. To a certain extent, this result is very useful because the NNN approximation is good enough to describe this system under appropriate parameters and it provides us an opportunity to make comparison between the atomic spin chain and Heisenberg spin chain.

3.1.3 Long-range case with the MDDI

To study the effects of the long-range spin coupling on the nonlinear spin wave dynamics, we need consider the spin coupling at all sites through the whole optical lattice, though the MDDI decreases rapidly in some cases. We firstly use Eq. (22) to rewrite the nonlinear motion equation of ϕnas

idϕndt=γBBz+jn2SJnjmagϕnjnSJnjzϕj2ϕnjn4SJnjϕj1jn+jn2SJnjϕn2ϕj1jn.E38

Then we set:

ϕnϕyt,E39
JnjmagJmagyy,E40
JnjJyy,E41
γBBz+jn2SJnjmagωy.E42

Considering the discreteness of the optical lattice, we can treat the symbolic terms as below:

1jn=cosjnπ=cos2πλLyy.E43

After changing the sum to integration, Eq. (38) becomes:

iytdt=ωyϕyt2SλLϕytJmagyyϕ2ytdy8SλLJyyϕytcos2πλLyydy+4SλLϕ2ytJyyϕytcos2πλLyydy.E44

Denoting yy=ξ, the ϕytcan be expanded as:

ϕyt=ϕyt+ϕytyξ+12!2ϕyty2ξ2+.E45

Taking above series into Eq. (44), the integration variable is changed, i.e., dy, and using the assumption ϕyt/y1, finally we get [10]:

iytdt=2β12ϕyty2+ωy4β0ϕyt+2β0γϕyt2ϕyt,E46

where we have used Jij=12Jijmag, and the coefficients are defined as below:

βn=SλLJmagξξ2ncos2πλLξ,n=0,1,E47
γ=2SλLJmagξ.E48

The condition for Eq. (46) to have envelope soliton solution is

fW=2β12β0γ>0.E49

From the definitions of βnand γ, we can easily get 2β0<γ. Then Eq. (49) reduces to β1>0, i.e.,

Jmagξξ2cos2πλLξ>0.E50

This result is consistent with the instability condition of modulational instability of nonlinear coherent spin wave modes near the Brillouin zone boundary under the long-wavelength modulation [4]. There, the modulational instability criteria take the corresponding discrete form:

14A2+9A316A4+<0,E51

where Aj=J0jmag/J01mag( j=2,3) is the relative strength of the longitudinal spin coupling of site jto that of NN.

The function fWversus the transverse width Wof the condensate is plotted in Figure 3 for three different approximations. fW=0marked in dashed line in Figure 3 characterize critical value, above which the magnetic solitons can occur. It is clear that the NN approximation and the NNN approximation take effect for a small transverse width.

Figure 3.

The existence conditions for magnetic solitons in atomic spin chain dominated by the MDDI. The three lines correspond to the NN approximation (blue), the NNN approximation (green), and the long-range case, respectively. The critical value of the transverse width is W1 and W2 [10].

3.2 Magnetic soliton under the LDDI

When we choose red-detuned optical lattice or an external modulational laser field is imposed into the system, both the static MDDI and the laser-induced LDDI begin to take effect on the dynamics of the soliton excitation, and the criterion for the existence of solitons could be changed. Here we also consider three interaction distances as doing in Section 3.1.

3.2.1 Nearest-neighbor approximation with the LDDI

In this case, we can obtain the continuum limit NLSE for ϕytwith the form:

iytdt=4SJ01a2d2ϕytdy2+γBBz+4SJ01mag+8SJ01ϕyt2SJ01mag+2J01ϕyt2ϕyt.E52

The existence condition of magnetic soliton now has the form

fW=8S2a2J01J01mag+2J01>0.E53

Obviously, magnetic soliton can be observed in this case, which is different from the case where only the MDDI plays role.

3.2.2 Next-nearest-neighbor approximation with the LDDI

In this case, the NLSE for ϕythas the form:

iytdt=4SJ0114A2d2ϕytdy2a2+ηϕyt4SJ0114A2ϕyt2ϕyt.E54

We assume that the intensity of the external laser is strong enough which corresponds to strong LDDI; then the terms like J0jmag/J0jhave been ignored. So for two adjacent sites, A2=J02/J01measures the relative strength of the NNN spin coupling to the NN spin coupling. The condition for Eq. (54) to have one soliton solution is

fW=16S2J01214A21A2>0.E55

3.2.3 Long-range case with the LDDI

In fact, once the external laser is imposed into the lattice, the LDDI dominates the system rapidly; the NLSE now takes the form [10]:

iytdt=2β12ϕyty2+ωy4β0ϕyt+2β0γϕyt2ϕyt,E56

where the coefficients are defined as below:

βn=2SλLJξξ2ncos2πλLξ,n=0,1,E57
γ=2SλLJmagξ.E58

The value of βndepends on both the MDDI and the LDDI. When a strong enough external modulated laser is imposed into the system, we have JmagξJξ, which corresponds to γ2β0. Thus the existence criteria of one soliton solution of Eq. (56) becomes β1β0<0, namely

fW=Jξξ2cos2πλLξJξcos2πλLξ<0E59

Similarly, Eq. (59) is consistent with the corresponding discrete form of MI condition of NCSW modes near the Brillouin zone boundary under the long-wavelength modulation [4]:

14A2+9A316A4+1A2+A3A4+<0,E60

where Aj=J0j/J01( j=2,3) is the relative strength of the transverse spin coupling of site jto that of NN.

We plot the existence regions of magnetic solitons in atomic spin chain dominated by the LDDI in Figure 4. The longitudinal coordinates stand for the intensity of the modulated laser. It is clear that the magnetic soliton can be achieved by tuning the transverse width and the modulated laser, which demonstrates the high controllability of optical lattice.

Figure 4.

The existence regions of magnetic solitons under the NNN approximation (a) and the long-range case (b) in atomic spin chain dominated by the LDDI, respectively. The critical value of the transverse width is Wa and Wb. The blank regions correspond to the existence of solitons for both cases [10].

4. Conclusion

In this chapter, we have shown the existence conditions of magnetic solitons that can occur in optical lattice. Compared to the more conventional solid-state magnetic materials, we discuss how the NN, the NNN, and the long-range interactions (the MDDI and LDDI) dominate the system. We also show that they can be tuned by the laser parameters and the shape of the condensate including the laser detuning, the laser intensity, and the transverse width of the condensate.

Besides studying the magnetic solitons, the atomic spin chain in optical lattice provides us with a useful tool to study the fundamental static and dynamic aspects of magnetism and lattice dynamics [12, 13, 14]. In experimental applications, the atomic spin chain has become potential candidates for multi-bit quantum computation due to their long coherence and controllability [13, 15, 16, 17]. The theoretical study of magnetic solitons in optical lattices will play a guiding role when the optical lattice is used in cold atomic physics, condensed matter physics, and quantum information.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant No. 11604086. Thanks to Ying-Ying Zhang for her contributions to physics discussion and English writing of the chapter.

Download

chapter PDF

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Xing-Dong Zhao (April 16th 2019). Magnetic Solitons in Optical Lattice [Online First], IntechOpen, DOI: 10.5772/intechopen.85492. Available from:

chapter statistics

82total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us