Эффект Зеемана в модифицированной модели Гросса--Невё в (2+1)-мерном пространстве-времени с компактификацией Текст научной статьи по специальности «Математика»

Эффект Зеемана в модифицированной модели Гросса-Невё в (2+1)-мерном пространстве-времени с компактификацией

Р. Н. Жоховa, В.Ч. Жуковскийb, П. Б. Колмаковс

Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра теоретической физики. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.

E-mail: a rzn95@mail.ru, b zhukovsk@phys.msu.ru, cpavel.b.kolmakov@yandex.ru

Статья поступила 04.03.2015, подписана в печать 20.04.2015.

В рамках модифицированной модели Гросса-Невё рассмотрено влияние магнитного поля, температуры и химического потенциала на фазовую структуру системы в (2+1)-мерном пространстве-времени с одним компактифицированным пространственным измерением (цилиндр). Рассмотрен эффект Зеемана, индуцированный внешним магнитным полем, параллельным оси цилиндра. Построены графики намагниченности системы, создаваемой за счет эффекта Зеемана. Указана возможность обобщения полученных результатов на углеродные нанотрубки.

Ключевые слова: модель Гросса-Невё, химический потенциал, эффект Зеемана.

УДК: 530.145. PACS: 11.10.Kk, 04.60.Kz, 11.10.Wx.

Введение

Модель Гросса-Невё (ГН) [1, 2] — одна из простейших, но в то же время крайне эффективных моделей четырехфермионного взаимодействия, которая будет рассматриваться в настоящей работе. С одной стороны, модель ГН нашла свое применение в качестве простой модели, демонстрирующей эффекты стандартной модели, такие как асимптотическая свобода, спонтанное нарушение симметрии и т. д. С другой стороны, известно, что многие результаты, полученные в модели ГН, могут быть использованы в физике конденсированного состояния вещества. Так, лагранжиан, подобный (1 + 1)-мерной модели ГН, возникает при рассмотрении одномерных полимеров (таких как полиацетилен [3, 4]), а лагранжиан, подобный (2+1)-мерной модели ГН, возникает при рассмотрении двумерных полимеров (таких как графен — см., например, [5]). Таким образом, особый интерес представляет рассмотрение модели ГН, используемой для описания подобных систем, находящейся под влиянием таких параметров, как температура, химический потенциал и внешние поля [5-9].

Если число измерений В = d + 1 больше, чем В = 1 + 1 (d > 1), то модель ГН является непере-нормируемой в рамках стандартной теории возмущений. Однако было показано, что трехмерная модель Гросса-Невё может быть перенормирована с использованием непертурбативного метода 1/Ы -разложения [10].

В настоящей работе модель ГН исследуется в (2+1)-мерном пространстве-времени с одним компактифицированным пространственным измерением (цилиндр). Такая модель в физике конденсированного состояния вещества могла бы описывать полимерные трубки (такие, как углеродные нанотрубки, которые могут быть рассмотрены как трубки графе-на). Рассматриваемый цилиндр будет вложен в обычное плоское четырехмерное пространство-время, при этом будет учитываться воздействие внешнего магнитного поля, направленного вдоль оси цилиндра на спин электронов (эффект Зеемана), и обусловленная

этим явлением намагниченность. Следует отметить, что в теориях, связанных с графеном, несмотря на введение уравнения Дирака для четырехкомпонент-ных фермионов (см. дальнейшие определения), их компоненты не связаны с реальным спином [11]. По этой причине при введении эффекта Зеемана реальный спин будет рассматриваться как дополнительное квантовое число (см. раздел 4), так же как это сделано, например, в работе [5] для плоского пространства.

В настоящей работе не будет рассматриваться вызываемый магнитным полем эффект Ааро-нова-Бома, влияние которого без учета эффекта Зеемана рассматривалось, в частности, в работах [8, 9, 12]. Одновременный учет эффектов Зеемана и Ааронова-Бома является предметом последующих публикаций.

1. Эффективный потенциал плоской трехмерной модели Гросса-Невё

Действие модели Гросса-Невё в плоском пространстве записывается в виде [1]

]Т^д^г + §

s =

d x

a=1

(1)

где N — число ароматов фермионов, а 7-матрицы определяются приводимым представлением:

7° =

0 ^ У = 0 -Л т2 = 0 -°2) {¡2 0) Ц 0 \а2 0 )' . (2) где а' — матрицы Паули, а матрица ¡2 — единичная матрица размера 2 х 2. В теориях, описывающих явления в графене, такое или схожие представления 7-матриц возникают естественным образом [13] при записи уравнения Гамильтона для квазичастиц, чем и обусловлен выбор этого представления в настоящей работе. При этом в теориях, связанных с гра-феном, N приобретает смысл числа спиновых состояний квазичастиц, а четыре компоненты ферми-

она ф не связаны с реальным спином [11], по этой причине в дальнейшем при расчете эффектов взаимодействия спина фермионов с магнитным полем различные спиновые состояния будут введены путем добавления нового квантового числа [5]. Киральная матрица 75 определяется стандартным образом:

75 =

Л2 (Л \0 -2/

(3)

Вводя вспомогательное поле о(х) = Nф(х)ф(х) и используя преобразование Стратоновича-Хаббар-да [14], можно привести действие к виду

5 =

¿3х

- - N ,

ф^Буф + ффа + — а

уг ' тт" ' 20' При этом производящий функционал

(4)

2 =

йфйфехр (¿Б [ф,Щ) = exp(-NVVeff),

где V — объем пространства-времени, с использованием приближения среднего поля а(х) = сопзЦх) = = а, являющегося ведущим порядком 1/N -разложения [15], позволяет определить эффективный потенциал Кей:

Veff = ^ - 2

20

с13р 1п [р2 + О2] .

(5)

В частности, а = 0 соответствует периодическим граничным условиям, а а = 1/2 — антипериодическим. Фактически, если разложить ф(х0, х1, х2) в ряд Фурье по х2, условие (6) означает, что ф(хо, х1, х2) = Е„ е2п1х2(п+а)/ь фп ^ ф(хо, х1, х2 + V) = = Е„ е2ж1х2(п+а)/Ь е'2п("+а) фп(х0, х1) = ф(х0, х1, х2)в2ж1а. В приложении к физике полимеров фаза а характеризует различные способы сворачивания листа полимера в трубку, приводящие к различным свойствам получающихся трубок; так, в частности, «металлическим» углеродным нанотрубкам соответствует а = 0, а «полупроводниковому» случаю соответствует а = ±1/3 [18]).

Рассмотрение модели при конечной температуре предполагает также компактификацию временного измерения; таким образом, вводя разложение Фурье также и по компактифицированному пространственному измерению, мы можем записать эффективный потенциал в виде

а

^ = ^ - 2

20

йр1 1 1 ^ V

1= — со

Ь

1п

2П) (1+1/2)2 + (2П) (п+а)2+р2+о5

(7)

2. Трехмерная модель Гросса-Невё на цилиндре

Рассмотрим цилиндр как (2+1)-мерную брану (с одним компактифицированным пространственным измерением), погруженную в плоское (3+1)-мер-ное пространство-время («балк»). Фермионы на бране подвергаются воздействию магнитного поля Н = И0гг постоянного и однородного во всем (3+1)-мерном балке. Здесь и далее в балке используются декартовы (х, у, г) и цилиндрические (р, ф, г) координаты, а на цилиндре — декартовы координаты (х1, х2), при этом компактифицированное измерение соответствует координате х2. Пусть радиус цилиндра равняется Я, а координатная ось г плоского (3+1)-мерного пространства совпадает с осью цилиндра: е1 = ег, е2 = вф. Векторный потенциал электромагнитного поля в балке записывается в виде А = 2вф, что приводит к потенциалу А = §е2 на цилиндре. Такой потенциал должен быть включен в ковариантную производную и приводит к возникновению эффекта Ааронова-Бома, однако в настоящей работе влияние данного эффекта рассматриваться не будет (в зависимости от геометрических параметров цилиндра и эффективных характеристик фермионов такое допущение может быть оправданным; в частности, магнитный поток через сечение цилиндра пропорционален квадрату радиуса цилиндра, и при достаточно малом радиусе даже значительные магнитные поля не создают существенного магнитного потока).

При указанном способе компактификации мы можем рассмотреть нетривиальные граничные условия с произвольной фазой а, что может быть записано как

ф(х0, хь х2 + V) = в2таф(х0, хь х2). (6)

где мы использовали обозначения ру = (р0, р1, р2), Ь = 2пЯ (длина окружности компактифицированного измерения), р2,п = (2^) (п + а) и учли конечную температуру Т = 1/в путем использования техники разложения по мацубаровским частотам [16] р0Д = ^ = (2п/в)(/ + 1/2) (I = 0, ±1, ±2,...).

Известно, что углеродные нанотрубки, проявляющие «металлические» свойства, для которых выполняются периодические граничные условия а = 0, достаточно устойчивы к нарушениям кристаличе-ской решетки и сохраняют свои «металлические» свойства даже в экстремальных условиях [17]. По этой причине в дальнейшем мы будем рассматривать только случай а = 0, однако большая часть результатов может быть обобщена и на случай а = 0.

Следует также отметить некоторые формальные свойства эффективного потенциала (7): эффективный потенциал является четной функцией а, кроме того, он периодичен по а с периодом 1, что позволяет ограничить исследуемую область интервалом а е [0,1/2], и далее с учетом свойств четности и периодичности возможно получить значение эффективного потенциала для любого действительного значения а, используя формулы ^(а,...) = ^(1+а,...) и ^(а,...) = Veff(-а,...).

Далее после добавления химического потенциала у в эффективный потенциал (7) (способы введения химического потенциала подробно изложены, например, в [16]) и в предположении а = 0, которое обсуждалось выше, эффективный потенциал может быть записан в виде

О2 2 ^й(о, а, Т, у) = 20 - Ж £

/ЗЬ ^ £

2п

х

X

2

п= —оо

X

х 1п

(I ('+1)-"')

2

2

+ (2?) + *+а

(8)

3. Фазовая структура модели с учетом химического потенциала и температуры

После выделения в (8) слагаемого, не зависящего от а, эффективный потенциал может быть переписан в виде

^ то

Уей= 20- вь £ £

1= — оо п=-то '

dpl 2п

1п

1 +

((I +1) - '»2 + (Ь )2(п+ф)2+р?

(9)

С использованием формулы

то ( Ь2 \ ,1:0+«гЬ '

то

\2 , „2 | и2

V (, + а)2 + а2 у = (1п ((, + а)2 + а2 + Ь2) - 1п ((I + а)2 +

а)2 + а2)) =

dт 1п 1 +

(1 + т 2 + а2)

+

+ 1п

1 - 2 cos(2nа)e-2^Va?+Ь? + е^^^ 1 - 2 ^^па^-2^ + е-4^

, (10)

которая может быть найдена, например, в [19], после ряда упрощений эффективный потенциал переписывается как

^ = К(0) + %) + V), (11)

где

к(0) = 20 - 2

d3p ( 2 2\ 1 (а3 а2а0

(2п)3 1п (р2 + = Л У -

/а3

и

у = А е^ А

-Ьап

п= 1

п=1

= 2(^ (е-Ьа) + (е-)) ,

(12)

(13)

V) = - л Е

dpl 2п

х 1п( 1 + 2 сЬ(Ди) ехр(-в£п,р1) ехр(-2в£п,р1)) -

- - У

/ЗЬ ^

2п

1п (1 + е-в(£п+Р1 )) +

)+,п(

1 +е

-в(Е-р1)

(14)

в [21]) характеризует границы области применимости модели, поскольку в рамках теории возмущений трехмерная модель Гросса-Невё является непере-нормируемой. Величина а0 является равновесным значением конденсата а при отсутствии температуры, компактификации и внешних полей [2]. При построении графиков будет использоваться характерный масштаб температур в = 1/Тс = 2|§| 1п 2/п, соответствующий температуре, при которой восстанавливается симметрия в плоской модели без внешних полей, а также совпадающий с ним в системе единиц Н = с = 1 масштаб длины для окружности компактифицированного измерения Ьс = рс.

Дальнейшие расчеты и построение графиков производились численно. График на рис. 1 демонстрирует фазовую структуру модели в плоскости (Т, /) для Ь = 1.2Ьс. Фазовая структура схожа со структурой, полученной в [22] для двумерной модели Гросса-Не-вё, что может быть интерпретировано как проявление размерной редукции из асимметричной фазы трехмерной модели Гросса-Невё, существующей при достаточно низкой температуре и а = 0 [2].

где и — специальные функции — полилогарифмы [20], Бп,р1 = у/р2 + (2Ь)2 + а2 и £±,1 = ^ ± ± /, а также введены обозначения а0 = , g = 1/0 - п/Л <0, где параметр обрезания импульса в плоской модели Л (достаточно подробно процедура перенормировки модели изложена, например,

0.2 0.4 0.6 0.8 1.0

Рис. 1. Фазовая диаграмма модели с учетом температуры и химического потенциала при Ь = 1.2Ьс

Области I и III на рисунке соответствуют фазе нарушенной симметрии, для которой а = 0, различие этих областей заключается в том, что в области I существует только один минимум эффективного потенциала , являющийся нетривиальным а = 0, в области III, однако существует два минимума — глобальный минимум а = 0 и локальный минимум при а = 0. Области II и IV являются областями с ненарушенной симметрией, при этом в области II существует только тривиальный минимум а = 0, в то время как в области IV присутствует также второй, нетривиальный минимум а = 0, являющийся, однако, только локальным. Как следствие линия ВЕ является кривой фазового перехода первого рода, а линия АВ — кривой фазового перехода второго рода (это различие будет наглядно продемонстрировано в разделе 4). Линии ВВ и ВС не являются кривыми фазовых переходов, а соответствуют линиям, при пересечении которых появляется либо исчезает локальный минимум.

х

2

а

2

е

2

п

х

х

п= —то

X

4. Магнитные свойства модели

Как указывалось выше, при рассмотрении эффектов взаимодействия описываемой системы с магнитным полем мы введем дополнительное квантовое число, соответствующее проекции спина фермионов на направление магнитного поля. Причиной такого подхода является связь модели с теориями, описывающими явления в графене, где четыре компоненты спинора ф не описывают реального спина частиц, а конденсат а является следствием нарушения кристаллической решетки и не связан с четырехфер-мионным взаимодействием. При этом, однако, эффективный потенциал таких теорий совпадает или частично совпадает (при добавлении дополнительных конденсатов) с (11)—(14) (см., например, [5]). Таким образом, для рассмотрения эффекта Зеемана в эффективный потенциал следует добавить стандартный член, по форме совпадающий с химическим потенциалом, но зависящий от направления проекции реального спина частиц на направление магнитного поля:

а2 1 тото

Keff = 2G - ßiYl Е Е

х ln

/2п

ve

s=t,| 1=—то п= — то 2

+

dpi 2п

(/+2)

(2P) +p2+^2

(15)

1

к-п=-^EE -2ПК1+2Ch(^)x

dpi

= -- У

ßL ^

s=t,| п=—то '

x exp(-ßEn,p1) + exp(-2ßEnp)) =

dp1 2n

ln

+ ln 1 +e

-e(En,p1);

(1 + e—^ )+t) +

^ + ln + e

-ß(En,p1 ) +

4 +

+ ln ^ 1 +e—ß(En p1

, (16)

M = -

dve

eff

dH '

Для рассматриваемой модели это означает

где введены обозначения у^ = у + 5у = у + ^УвИ, у^ = у — 5у = у — &т уВИ, в которых ув = в/2т — магнетон Бора и gL — эффективный множитель Ланде.

По аналогии с формулой (14) для учета эффекта Зеемана требуется заменить УуТ) на

M=

dVeff dVeff(CT, ß, L, у, 5у) dVeff dSy

^eff

dH

dH

дбу dH '

(17)

(18)

Согласно данному выше определению,

дёу = gLУв дИ 2 .

Можно ожидать, что эффект Зеемана должен повлиять на восстановление киральной симметрии модели при характерном значении 5у ~ о0, где о0 — опреденный выше конденсат в плоской модели без внешних полей. Магнитное поле Иу, создающее соответствующее значение 5у, может быть оценено как

5у — сто - n/g, 2ctq

Hу —

gLyB

(19)

(20)

Формулы (12), (13), (16) позволяют рассчитать намагниченность, вызываемую эффектом Зеемана. Легко видеть, что из всех слагаемых, входящих в ^, только У(уТ) содержит зависимость от 5у. Для еще большего упрощения возможно отдельно провести вычисления для слагаемых У(уТ), соответствующих 5 и 5 =1, таким образом предварительно вычислив концентрацию частиц с различными направлениями спина:

дУ(

ns =

(ут)

ду8

= 1 ТО \dp1 (

= L ^ 2п I

п= — пп ^

sh(вУs)

2п Vch(вys) + ch(вEn,pl)

(21)

где (Еп,р1 = Еп,р1 ± уц. Из этой формулы можно сделать вывод, что значение эффективного потенциала остается неизменным при замене у ^ 5у, а следовательно, фазовая структура, полученная в разделе 3 на плоскости (Т, у), совпадает с фазовой структурой в плоскости (Т, 5у) при у = 0.

Рассмотрим намагниченность, обусловленную взаимодействием фермионов с внешним магнитным полем за счет эффект Зеемана. В общем случае намагниченность определяется как частная производная эффективного потенциала:

где производные берутся в точке стабильного минимума ст eq, для которой dVeff/дст|ст=ст = 0

и д2VVff/дст21ст >0.

1 ст —CTeq

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

дМ

М = ув(щ - пф), х = дн (22)

где формула для намагниченности фактически совпадает с (17), но использует промежуточный результат (21).

График на рис. 2 демонстрирует зависимость намагниченности, вызванной эффектом Зеемана, от величины магнитного поля при L = 0.7Lc < Lc, у = 0. На данном графике хорошо видно, что фазовый переход между кирально симметричной и несимметричной фазами является переходом первого рода, если в симметричной фазе сохраняется локальный нетривиальный минимум (переход из области III в область IV и обратно на рис. 1), и, напротив, переходом второго рода, если нетривиальный минимум исчезает при фазовом переходе (переход из области I в область II и обратно на рис. 1). При T = 1 Tc с ростом магнитного поля график пересекает кривую фазового перехода первого рода (BE на рис. 1), при T = 4 Tc — кривую фазового перехода второго рода

х

ст=ст

eq

)

ст=ст

eq

п=

М

Рис. 2. Намагниченность системы М как функция магнитного поля Н при различных температурах: Т = 3Тс (штриховая кривая), Т = |Тс (сплошная кривая), Т = 2Тс (точечная кривая)

(АВ на рис. 1), а при Т = 3Тс кривые фазовых переходов не пересекаются, поскольку симметрия оказывается восстановлена даже при Н = 0.

Заключение

В настоящей работе исследована роль эффекта Зеемана, вызванного внешним магнитным полем, и химического потенциала в модели Гросса-Невё в (2+1)-мерном пространстве-времени с одним компактифицированным пространственным измерением (цилиндр) с периодическими граничными условиями при конечной температуре. Получены явные выражения для эффективного потенциала (11)—(13), (16).

Построена фазовая диаграмма модели с учетом конечной температуры и химического потенциала (рис. 1), совпадающая также с фазовой диаграммой модели с учетом эффекта Зеемана при нулевом химическом потенциале. Построенная диаграмма принципиально совпадает с фазовой диаграммой двумерной модели Гросса-Невё [22], что является проявлением размерной редукции в трехмерной модели Гросса-Невё с компактификацией при периодических граничных условиях [2].

Получены формулы (21, 22) для расчета намагниченности, индуцированной эффектом Зеемана, вызванным магнитным полем, параллельным оси цилиндра. Построены графики зависимости намагниченности от приложенного магнитного поля при

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

Авторы искренне благодарят Д. Эберта за поддержку работы и многочисленные советы, а также участников научного семинара кафедры теоретической физики физического факультета МГУ за обсуждение работы и полезные замечания.

Список литературы

1. Gross D.J., Neveu A. // Phys. Rev. D. 1974. 10. P. 3235.

2. Bietenholz W., Gfeller A., Wiese U.-J. // JHEP. 2003. 10. P. 018.

3. Caldas H. // Nucl. Phys. B 2009. 807 [FS]. P. 651.

4. Roscher D., Braun J., Drut J.E. // Phys. Rev. A. 2014. 89. P. 063609.

5. Caldas H, Ramos R.O. // Phys. Rev. B. 2009. 80. P. 115428.

6. Novoselov K.S., Geim A.K., Morozov S.V. et al. // Nature. 2005. 438. P. 197.

7. Vshivtsev A.S., Zhukovsky V.Ch., Klimenko K.G., Mag-nitsky B.V. // Phys. Part. Nucl. 1998. 29. P. 523.

8. Ebert D., Klimenko K.G., Tyukov A.V., Zhukovsky V.Ch. // Phys. Rev. D. 2008. 78. P. 045008.

9. Ebert D., Khunjua T.G., Klimenko K.G., Zhukovsky V.Ch. // Int. J. Mod. Phys. A. 2012. 27. P. 1250162.

10. Rosenstein D, WarrB.J., Park S.H. // Phys. Rep. 1991. 205. P. 59.

11. Castro Neto A.H., Guinea F., Peres N.M.R. et al. // Rev. Mod. Phys. 2009. 81 . P. 109.

12. Жуковский В.Ч., Колмаков П.Б. // Вестн. Моск. ун-та. Физ. Астрон. 2013. № 4. С. 8 (Zhukovsky V.Ch., Kolmakov P.B. // Moscow University Phys. Bull. 2013. 68. N 4. P. 272).

13. Gusynin V.P., Sharapov S.G., Carbotte J.P. // Int. J. Mod. Phys. B. 2007. 21. P. 4611.

14. Стратонович Р.Л. // Докл. АН СССР. 1957. 157. С. 1097; Hubbard J. // Phys. Rev. Lett. 1959. 3. P. 77.

15. Rosenstein B, WarrB.J., Park S.H. // Phys. Rep. 1991. 205. P. 59.

16. Kapusta J.I, Gale C. Finite-temperature Field Theory: Principles and Applications. Cambridge, 2006.

17. Ando T. // J. Phys. Soc. Japan. 2005. 74. P. 777.

18. Elizalde E., Odintsov S.D., Saharian A.A. // Phys. Rev.D. 2011. 83. P. 105023.

19. SongD.Y. // Phys. Rev. D. 1993. 48. P. 3925.

20. Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. Т. 1. Элементарные функции. М., 1981.

21. Ramos R.O., Manso P.H.A. // Phys. Rev. D. 2013. 87. P. 125014.

22. Schoen V., Thies M. // At the Frontier of Particle Physics: Handbook of QCD: Boris Ioffe Festschrift. Vol. 3. World Scientific, 2001. P. 1945.

The Zeeman effect in a modified Gross-Neveu model in (2+1)-dimensional space-time with compactification

R.N. Zhokhova, V.Ch. Zhukovskyb, P. B. Kolmakovc

Department of Theoretical Physics, Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia.

E-mail: a rzn95@mail.ru, b zhukovsk@phys.msu.ru, cpavel.b.kolmakov@yandex.ru.

The influence of the magnetic field, temperature, and chemical potential on the phase structure of a system in (2+1)-dimensional space-time with one compactified spatial dimension (cylinder) was studied in the framework of the modified Gross-Neveu model. The Zeeman effect that is induced by an external magnetic field parallel to the cylinder axis is considered. The system magnetization that is created due to the Zeeman effect was plotted. The possible generalization of these results to carbon nanotubes is noted.

Keywords: Gross-Neveu model, chemical potential, Zeeman effect. PACS: 11.10.Kk, 04.60.Kz, 11.10.Wx.

Received 4 March 2015.

English version: Moscow University Physics Bulletin 4(2015).

Сведения об авторах

1. Жохов Роман Николаевич — аспирант; e-mail: rzn95@mail.ru.

2. Жуковский Владимир Чеславович — докт. физ.-мат. наук, профессор; тел.: (495) 939-31-77, e-mail: zhukovsk@phys.msu.ru.

3. Колмаков Павел Борисович — аспирант; e-mail: pavel.b.kolmakov@yandex.ru.