Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Кашурников Сверхтекучест и бозе-конденсация 2008

.pdf
Скачиваний:
117
Добавлен:
16.08.2013
Размер:
5.25 Mб
Скачать

известным в теории солитонов, естественно предположить, что солитоны могут быть и в конденсате. Здесь целесообразно рассмотреть основы теории бозе-конденсата, помещенного в ловушку.

4.1.1. Уравнение Гросса – Питаевского

Многочастичный гамильтониан, описывающий парно взаимодействующие бозоны, находящиеся во внешнем

потенциальном

 

поле

 

 

 

 

,

 

в

представлении

вторичного

 

 

 

записываетсяœ˜®как [96, 97]

 

 

 

квантования

 

 

 

 

 

 

 

¨

'RÆ)

 

Ý

 

b ¨ 'RÆ)W h'RÆ) b

 

 

 

= ‘ ’RÆ h

M

'RÆ

)

VW +

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

œ˜®

 

 

 

 

h'RÆ)

‘ ’RÆ’Rì hM'RÆ)hM'RÆì)¨'RÆ W RÆì)h'RÆì)h'RÆ)& 'F( )

где

и

+

h

M

'RÆ)

 

операторы

бозонного поля

(полевые

 

 

 

 

бозонов.

 

¨'RÆ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

операторы),

уничтожающие и

рождающие бозон

в

точке ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

взаимодействие

потенциал

 

 

 

 

 

описывает

двухчастичное

 

 

Приближение среднего поля (приближение БоголюбоваhX'RÆ& !) ) состоит в представлении полевого оператора , взятого в гейзенберговском представлении, как суммы оператора рождения частиц в основном состоянииhX'RÆ& !) =и1Xво'RÆ&всех!) bостальныхhX ì'RÆ& !)( состояниях: 'F(+)

Предполагая, что макроскопическое1X'RÆ&число!) частиц находится в

основном состоянии, 1оператор'RÆ& !) заменяется числовой

комплексной функциейhX'RÆ& !) :1'RÆ& !) b hX ì'RÆ& !)& 'F(*)

так что 1'RÆ& !) = •hX'RÆ& !)w. Эта функция называется волновой функцией конденсата и имеет смысл параметра порядка. Она характеризует поведение недиагональных элементов

131

при увеличении расстояния

 

 

 

 

_ì'RÆ& RÆ

A !)

= •h

M

'RÆ& !)h'RÆ

& !)w

одночастичной матрицы плотности

 

g

ì

 

 

 

 

ì

 

и уничтожения возбужденийR = gRÆ W RÆ

 

 

hX'RÆ& !)

 

 

 

 

>

 

!

 

1

 

 

 

 

между точками рождения

|

_'RÆ& RÆìA

) =

Á'RÆ&

! 1

'RÆì&

!

)( 'F(F)

 

 

 

 

)

 

 

 

 

квантового поля

 

 

 

 

:

 

 

Ùž“

Свойство матрицы плотности (4.4) является математическим выражением существования недиагонального дальнего порядка в

системе.

Уравнение

движения

для полевого

оператора

h'RÆ)

можно

 

 

 

Ž

 

 

 

 

 

записать, исходя из уравнения Гейзенберга

 

что дает

Ž

k Ž! h'RÆ) =

Ú & h'RÆ)Ü& 'F(f)

 

 

 

Ý

 

b ¨

'RÆ)W h'RÆ) b

 

 

k Ž! h'RÆ

)

= VW +

 

 

 

 

 

 

 

œ˜®

 

 

b ‘ ’RÆì hM'RÆì)¨'RÆ W RÆì)h'RÆì)h'RÆ)( 'F( )

1Если'RÆ& !теперь) = •hX'заменитьRÆ& !)w операторы h'RÆ) их средними значениями (с учетом малых флуктуаций плотности конденсата), пренебрегая всеми корреляторами, то получится

уравнение движенияŽ классического поля

k Ž! 1'RÆ) = VW + Ý b ¨œ˜®'RÆ)W1'RÆ) b

b ‘ ’RÆì 1Á'RÆì)¨'RÆ W RÆì)1'RÆì)1'RÆ)( 'F(,)

В разреженном газе холодных атомов можно ограничиться только

взаимодействиями,

ì)

обусловленными столкновениями между

атомами, положив

$

ì)

F>

ì)

Параметр

 

¨'RÆназываетсяW RÆ = "длинойp'RÆ W RÆ

рассеяния=

p'RÆ-волны;W RÆ ( – масса'F(2)

 

Этот параметр зависит от деталей атомного потенциала и

атома.

 

 

 

 

 

 

 

ø

 

 

 

 

 

132

 

 

может быть положительным= +(,fилиcDотрицательным=.f(,,Напомним,cD что для найдено, что , для , для

= W (Ff cD

.

В таком приближении локального взаимодействия между бозонами уравнение движения для волновой функции конденсата

Ž

 

 

 

 

 

 

 

принимает вид

 

Ý b ¨

'RÆ)W1'RÆ) b " 1'RÆ)g 1'RÆ)( 'F(O)

k Ž! 1'RÆ

)

= VW +

 

 

 

 

œ˜®

$g

 

Это и есть уравнение Гросса – Питаевского [98, 99].

Уравнение Гросса – Питаевского может быть представлено в

гамильтоновой форме

k Ž1 = p!&Á & 'F( /)

 

!&

 

 

 

 

 

 

 

 

 

 

"$

 

 

 

 

 

 

 

 

 

 

 

Гросса – Питаевского,

 

 

 

 

где

 

 

 

функционал энергииŽ!

p1

b + g1'RÆ)g

¶( 'F( )

!& = ‘ ’RÆ

µ+

1'RÆ)g

b ¨œ˜®'RÆ)g1'RÆ)g

быть

 

 

 

¨ 'RÆ)

 

 

 

 

 

 

 

Потенциал

 

 

, удерживающий атомы конденсата, может

 

 

весьма

œ˜®произвольным. Часто в качестве функции

 

 

 

 

 

 

 

 

 

 

 

 

 

приближениеœ˜®

берется параболический потенциал, как хорошее

 

 

¨

'RÆ)

реального потенциала ловушки. Действительно, вблизи минимума потенциальной ямы ¨или'RÆ)около максимума потенциального барьера разложение œ˜® по степеням отклонения от точки экстремума во втором порядке дает параболическую зависимость

от координат. Частными случаями его являются: сферически симметричный потенциал

несимметричные

¨œ˜®'R) =

 

 

ð@R A 'F( +)

 

 

потенциалы, например

для

сигароподобной

ловушки он имеет форму

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

œ˜®

'Â& Y& ”) =

 

 

 

b Y

 

 

¨

+

»ðZ

 

 

b ð[

 

¼( 'F( *)

 

 

133

 

 

 

 

 

 

 

 

Непараболический удерживающий потенциал соответствует случаю, когда два облака конденсата разделены потенциальным барьером конечной высоты и слабо связаны из-за туннелирования атомов сквозь барьер. Такой потенциал можно аппроксимировать

следующим выражением:

 

b Y

 

b ð[

 

b ß'”

 

W é

 

¼( 'F( F)

¨

'Â& Y& ”) = + »ðZ

 

)

 

 

)

œ˜®

 

 

 

 

 

 

 

 

 

 

 

 

Другой пример непараболического потенциала ловушки – потенциал оптической решетки – будет рассмотрен позже.

4.1.2.Бозе-эйнштейновская конденсация в квазиодномерной параболической ловушке. Одномерное уравнение Гросса – Питаевского

Предположим, что потенциал, удерживающий конденсат, описывается формулой (4.13). В этом случае уравнение Гросса – Питаевского записывается как

Ž

 

Ž

 

 

 

 

 

 

 

 

 

 

 

 

 

k Ž! 1'RÆ) b +

Ž” 1'RÆ) W

+

»ð

Z

 

b Y

) b ð[

¼1'RÆ) W

 

 

$g

1

'

)g

 

'

)

= /( 'F( f)

 

 

W"

 

 

 

1

 

 

Удобно ввести нормированные переменные и переписать это уравнение в безразмерной форме. Выберем в качестве нормированных переменных\ = ðZ!A

] = Âá)A

 

 

^ = Yá)A

 

 

1'Â& Y& ”& !

= ”Þ)A

'

)

)

= •

V

 

 

_ ]& ^&

& \ ( 'F( )

В этих переменных уравнение Гросса – Питаевского выглядит так:

134

 

 

k

Ž_

b

 

âZ_ b

 

d

á)

e

Ž _

W

 

 

Ž\

+

+

Þ)

Ž

 

 

 

W

 

 

 

 

 

 

)_ W

"$

 

 

 

где

 

+ 'RZ b ;

 

ðZ g_g _ = /& 'F( ,)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

RZ

á) =

ðZ

A

 

 

 

 

 

 

 

 

 

= '] ðb[á)^ )VA

 

 

 

 

 

 

 

 

; = dðZÞ)e

 

A

 

 

 

 

 

 

 

 

 

âZ =

Ž

 

 

 

 

Ž

 

 

 

á)

 

 

 

 

Ž] b Ž^ ( 'F( 2)

Длина

характеризует поперечный размер облака конденсата,

 

который считается много меньшим продольного размера ловушки

Þ). Если ввести корреляционный радиус (healing length)

 

 

 

C = µ

+

 

 

V

= '2>•g g)‚V& 'F( O)

g

"

 

 

 

 

 

 

 

 

 

 

 

 

$g

в (4.17)$

 

á)

 

 

 

 

 

 

"

 

 

 

 

 

 

 

 

 

 

можно записать как

 

 

 

то константу нелинейности g

g

 

 

 

 

 

 

 

 

 

 

 

ðZ = `}' ) + d

 

e ( 'F(+/)

C

Если плотность атомов конденсата

 

определяется как интеграл

 

 

 

 

функцииè

конденсата, то параметр

 

от квадрата модуля волновой

 

 

 

 

 

 

 

 

 

 

связан с плотностью атомов соотношением

 

 

 

 

è = ‘g1'RÆ)g

’ R = • ‘g_'RÆ)g ’ R( 'F(+ )

 

 

 

 

 

_'RÆ)g

 

 

 

 

 

 

 

Выбрав нормировку функции

g

 

, можно

выразить

 

через

плотность атомов конденсата.

 

 

 

 

Если поперечный размер ловушки велик по сравнениюá ¬ с продольным и взаимодействие между атомами слабое, ) C, то структура облака конденсата в поперечном направлении

135

определяется самой ловушкой, и распределение в поперечном

Ž_

 

RZ

 

направлении задается линейным уравнением

k Ž\ b +

âZ_ W +

_ = /( 'F(++)

Эта задача напоминает поиск поперечного распределения электромагнитного поля в волоконном световоде. Таким образом, надо найти поперечные моды (решения линейного уравнения) и разложить искомое решение уравнения Гросса – Питаевского (4.17) по этим модам, полагая, что коэффициенты разложения являются функциями времени и продольной координаты.

 

 

 

_'RÆ& !

)

= ö'

)

 

' )

 

 

‚y€a

( 'F(+*)

Пусть

 

 

 

& !

S ]& ^

 

 

 

 

Žö á)

 

Ž ö ;

 

"$

 

 

 

 

 

 

Подстановка этого разложения в уравнение (4.17) дает

 

 

 

µk

Ž\

b + dÞ)e

Ž W

+

 

ö W ðZ gSög ö¶ S b

 

 

 

 

 

 

 

 

 

)

ö = /( 'F(+F)

 

 

 

b + ZS W RZ S b +ŠS

 

и используя краевые

âZS W RZ S b +ŠS = / 'F(+f)

Положив

 

 

 

 

 

 

 

 

 

 

= ]

 

b ^

 

 

 

 

 

 

 

 

 

S ž /

 

 

 

RZ

 

V

ž À

 

 

 

 

условия

 

 

 

при

 

 

 

ловушки.

,

можно найти поперечные моды сигарообразной '

 

 

)

 

Простое приближение соответствует выбору основной моды, для

которой

 

 

 

 

S$'RZ) = £

ÎÙ«

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b ( 'F(+ )

Умножив

(4.24)

на

Ž ö ;

и

 

интегрируя

по поперечным

 

Žö á)

 

 

"$•¸

 

 

 

 

 

 

 

уравнение$

 

 

 

 

 

координатам, получим

S 'RZ)

 

 

ö W ðZ gög ö = /&

 

 

k Ž\ b

+ d

 

e

Ž W +

 

 

'F(+,)

где

C

 

 

 

 

 

 

 

136

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S$

'RZ)RZ@RZ( 'F(+2)

 

 

 

 

 

 

 

 

 

 

¸ = +> ‘

 

 

 

 

 

 

 

 

 

 

 

 

 

$

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выбирая должным образом параметры

 

,

 

и

 

 

, можно привести

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

виду одномерного уравнения

уравнение (4.27) к каноническому

 

 

• ¸

 

 

Þ)

 

 

 

 

 

 

 

 

 

Žö

 

Ž ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гросса – Питаевского:

b

 

W

 

ö W Ugög ö = /( 'F(+O)

 

 

 

 

 

k

Ž\

 

 

Здесь

 

 

 

 

 

U = `}'"

= `}'

 

 

( 'F(*/)

 

 

 

 

 

 

Ž

 

 

 

$)

 

 

 

 

 

$)

 

 

 

 

Как исходное

уравнение

Гросса – Питаевского,

 

так и его

одномерная версия могут быть получены

в

форме уравнений

 

 

 

 

' d)

 

k

 

Á Ž

 

 

 

Ž

 

Á

 

 

 

 

 

 

 

Эйлера – Лагранжа вариационной задачи с лагранжианом

 

 

c

 

=

+ d1

Ž!

1 W

1 Ž! 1 e b +

 

 

1

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"$

 

 

A

 

 

 

g

 

 

 

'V

 

)

k

 

 

Á

 

®Ùè¦'RÆ)gÁ1g

W + g1g

 

 

 

 

U

c

d

 

 

Žö

W ö

Žö

 

 

 

Žö

e W

 

 

 

 

 

 

= +

 

 

Ž\

Ž\

e b eŽ

 

 

gög W

+ gög ( 'F(* )

Функционалы действия ¤ = ‘ c' d)Ú1Ü ’RÆ’! 'F(*+)

или

¤ = ‘ c'Vd)ÚöÜ ’\’ 'F(**)

инвариантны относительно сдвигов по оси времени и умножения полей на фазовый множитель. Но функционалы действия не инвариантны по отношению к трансляции вдоль осей координат из-за потенциала ловушки.

137

4.1.3. Одномерное нелинейное уравнение Шредингера

Предельный случай сильноанизотропной длинной сигарообразной ловушки может быть описан по аналогии с электромагнитными волнами в волоконных световодах, где профиль (необязательно параболический, но обычно цилиндрически симметричный) показателя преломления волокна играет роль удерживающего потенциала, тогда как в продольном направлении потенциальные стенки разнесены на огромное расстояние, и их влиянием можно пренебречь. Конденсат, сосредоточенный внутри такой вытянутой ловушки, может быть описан одномерным уравнением Гросса – Питаевского, где

линейное по полю слагаемоеŽö Ž отöсутствует:

k Ž\ b Ž W Ugög ö = /( 'F(*F)

Это известное одномерное нелинейное уравнение Шредингера, широко используемое в теории нелинейных волн в самых разнообразных областях физики. Опираясь на эту аналогию, можно ожидать существования квазиодномерных «волн вещества» в форме ярких солитонов, подобных временным солитонам в нелинейных оптических волокнах, если

взаимодействие атомов притягивающее (

 

 

). И действительно,

 

 

 

 

 

 

 

:,

где взаимодействие

яркие солитоны в конденсате атомов

/

 

независимо

 

 

 

 

 

притяжения,

были

атомов

имеет

характер

 

 

 

 

 

 

 

 

 

 

 

обнаружены

в

двух

экспериментах

[100, 101].

До

 

этого

в

конденсатах

с

отталкивающим

взаимодействием

(

,

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

которые

удалось

сформировать

темные

солитоны

[102, 103],

 

 

 

выглядят как провал в однородном фоне концентрации атомов.

Солитонные решения нелинейного уравнения Шредингера образуют множество,]& Q характеризуемое набором непрерывных параметров . Наиболее известным является фундаментальный солитон, имеющий следующий вид:

138

 

]

]

ö'\& ) = WU

æü V + ' W Q\)W å

å ij ´k

µQ

b µ] W Q \¶ ( 'F(*f)

+

 

+

Это решение нелинейного уравнения Шредингера иногда называют ярким (или светлым) солитоном (рис. 4.1), подчеркивая тем \самым& ž mÀто, что он располагается на нулевом фоне, т.е. исчезает при .

Рис. 4.1. Яркий солитон – решение нелинейного уравнения Шредингера

 

 

 

 

 

 

(4.35) с параметрами

] = + +, Q = f

 

 

 

 

 

 

 

Яркий солитон существует при

 

 

 

(параметр

 

определен в

 

 

 

 

 

 

 

 

 

 

 

 

между$

атомами конденсата$

 

 

, и

(4.8)). В случае отталкивания

 

"

: /

 

 

 

 

 

"

 

солитонных$

нелинейное

 

 

 

 

 

 

 

Шредингера

 

 

 

 

 

 

 

 

 

 

 

 

уравнение

 

 

 

 

 

 

не

 

имеет

 

 

 

"

\ /

 

 

 

 

 

gö'

 

 

,

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

!& ” ž mÀ

 

 

 

 

 

& \)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решений

на

нулевом фоне.

Если при

 

 

 

 

 

 

плотность

g

\)

g

ž gö

g

 

 

 

конечная и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö' &

 

 

 

 

 

постоянная, т.е. при

 

 

 

 

конденсата

 

 

 

 

 

 

 

 

 

 

 

\&

ž mÀ

 

 

 

 

 

 

 

 

 

 

 

 

 

существует

решение

нелинейного

 

уравнения

Шредингера, называемое$

]

 

 

 

]

 

 

 

 

 

 

 

 

 

 

 

 

темным солитоном (рис. 4.2):

 

 

 

 

 

 

 

 

ö'\&

 

) = m WU Ö }ü V + '

W Q\)W å

 

 

 

 

 

 

 

 

 

 

 

 

 

å ij ´W

k

 

W µ+] U b

Q

\¶ ( 'F(* )

 

 

 

 

 

 

 

+ µQ

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

139

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 4.2. Темный солитон – решение нелинейного] = + + Qуравнения= f Шредингера (4.36) с параметрами ,

Темные солитоны, как и светлые, упруго взаимодействуют между собой и устойчивы (хотя и неасимптотически) по отношению к слабым возмущениям. Столкновение двух фундаментальных солитонов описывается двухсолитонным решением нелинейного уравнения Шредингера [103].

Следующим шагом в сторону реальной ситуации является учет удерживающего потенциала в продольном направлении. Вместо нелинейного уравнения Шредингера рассматривается одномерное уравнение Гросса – Питаевского для конденсата в параболической ловушке (4.29).

Нелинейного уравнения Шредингера недостаточно для описания неоднородных состояний бозе-эйнштейновского конденсата, где ловушка может играть важную роль, но как модель, позволяющая предсказать ожидаемые явления, данное уравнение может оказаться полезным.

140