Кашурников Сверхтекучест и бозе-конденсация 2008
.pdfизвестным в теории солитонов, естественно предположить, что солитоны могут быть и в конденсате. Здесь целесообразно рассмотреть основы теории бозе-конденсата, помещенного в ловушку.
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Æ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
операторы), |
уничтожающие и |
рождающие бозон |
в |
точке , |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
взаимодействие |
||
потенциал |
|
|
|
|
|
описывает |
двухчастичное |
|
|
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Æ |
µ+ |
gÝ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" |
|
RÆ |
|
|
1 RÆ |
|
|
Удобно ввести нормированные переменные и переписать это уравнение в безразмерной форме. Выберем в качестве нормированных переменных\ = ð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 + |
|
|
gÝ1 |
|
W |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
"$ |
|
|
A |
|
|
|
g |
|
|
|
|
'V |
|
) |
k |
|
|
Á |
|
W¨®Ùè¦'RÆ)gÁ1g |
W + g1g |
|
|
|
|
U |
||||||||||||
c |
d |
|
|
Žö |
W ö |
Žö |
|
|
|
Žö |
e W |
|
|
|
|
||||||||||||
|
|
= + dö |
|
|
Ž\ |
Ž\ |
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