15.04.2013 Views

Bose-Einstein Condensates in Rotating Traps and Optical ... - BEC

Bose-Einstein Condensates in Rotating Traps and Optical ... - BEC

Bose-Einstein Condensates in Rotating Traps and Optical ... - BEC

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

48 Groundstate of a <strong>BEC</strong> <strong>in</strong> an optical lattice<br />

the effects of harmonic trapp<strong>in</strong>g have been presented <strong>in</strong> [102]. In this thesis, also results for the<br />

density profile are <strong>in</strong>cluded <strong>and</strong> the dependence on lattice depth of the density at the center<br />

of the harmonic trap averaged over one lattice period is described. We briefly comment on<br />

the generalization to 2D lattices. Based on an ansatz for the order parameter valid <strong>in</strong> deep<br />

lattices, the effective coupl<strong>in</strong>g constant ˜g <strong>and</strong> its effect on the harmonically trapped system<br />

have been previously discussed <strong>in</strong> [70].<br />

5.1 Density profile, energy <strong>and</strong> chemical potential<br />

When look<strong>in</strong>g for stationary solutions, the GP-equation (3.14) takes the form<br />

<br />

− ¯h2 ∂<br />

2m<br />

2<br />

∂z2 + sER s<strong>in</strong> 2<br />

<br />

πz<br />

+ g |Ψ(z)|<br />

d<br />

2<br />

<br />

Ψ(z) =µΨ(z) . (5.1)<br />

The groundstate is given by the solution of this equation with the lowest energy.<br />

Let us rewrite Eq.(5.1) <strong>in</strong> a more convenient form. First, we <strong>in</strong>troduce the rescaled order<br />

parameter<br />

<br />

ϕ(z) =<br />

L2 Ψ(z) ,<br />

N<br />

(5.2)<br />

with L the transverse size of the system <strong>and</strong> N the number of particles per lattice site such<br />

that d/2<br />

dx |ϕ(z)|<br />

−d/2<br />

2 The GPE (5.1) then reads<br />

=1. (5.3)<br />

<br />

− ¯h2 ∂<br />

2m<br />

2<br />

∂z2 + sER s<strong>in</strong> 2<br />

<br />

πz<br />

+ dgn|ϕ(z)|<br />

d<br />

2<br />

<br />

ϕ(z) =µϕ(z) , (5.4)<br />

where we have <strong>in</strong>troduced the average density<br />

n = N<br />

. (5.5)<br />

dL2 Cast<strong>in</strong>g Eq. (5.4) <strong>in</strong> dimensionless form, <strong>in</strong> dimensionless form it is possible to identify the<br />

govern<strong>in</strong>g parameters of the problem. We shall measure length <strong>in</strong> units of d/π, momentum<br />

<strong>in</strong> units of the Bragg-momentum qB =¯hπ/d <strong>and</strong> energy <strong>in</strong> units of the recoil energy ER =<br />

¯h 2 π 2 /2md 2 . In this way, we obta<strong>in</strong> the dimensionless GPE<br />

<br />

− ∂2<br />

∂z 2 + s s<strong>in</strong>2 (z)+π gn<br />

|ϕ(z)|<br />

ER<br />

2<br />

<br />

ϕ(z) = µ<br />

ER<br />

ϕ(z) . (5.6)<br />

In contrast to the case of a s<strong>in</strong>gle particle discussed <strong>in</strong> chapter 4, there are now two govern<strong>in</strong>g<br />

parameters: the lattice depth s <strong>and</strong> the <strong>in</strong>teraction parameter gn/ER. The latter quantity can<br />

be changed by vary<strong>in</strong>g the average density or the scatter<strong>in</strong>g length by means of a Feshbach<br />

resonance.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!