Критическое поведение o(n)-симметричной модели с антисимметричным тензорным параметром порядка: четырёхпетлевое приближение Текст научной статьи по специальности «Физика»

УДК 53

Вестник СПбГУ. Физика и химия. 2018. Т. 5 (63). Вып. 1

Н. М. Лебедев, М. В. Компаниец

КРИТИЧЕСКОЕ ПОВЕДЕНИЕ 0(п)-СИММЕТРИЧНОЙ МОДЕЛИ С АНТИСИММЕТРИЧНЫМ ТЕНЗОРНЫМ ПАРАМЕТРОМ ПОРЯДКА: ЧЕТЫРЁХПЕТЛЕВОЕ ПРИБЛИЖЕНИЕ*

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

При помощи ренормгруппового подхода рассматривается критическое поведение 0(га)-симметричной модели ф4 с антисимметричным тензорным параметром порядка. Ранее данная модель изучалась в рамках трёхпетлевого приближения посредством методов е-раз-ложения с пересуммированием конформ-борелевским методом и псевдо-е-разложения. Было показано, что результаты данных подходов качественно согласуются друг с другом, однако полученные значения критических индексов могут довольно существенно отличаться. Более того, оставался открытым вопрос о чувствительности даже качественных результатов по отношению к учёту следующих порядков теории возмущений. В настоящей работе ренормгрупповые функции и критические индексы модели вычислены с четырёхпетлевой точностью в рамках обоих подходов. Их анализ показывает, что результаты применения различных подходов качественно согласуются и между собой и с результатами трёхпетлево-го анализа. Библиогр. 29 назв. Табл. 3.

Ключевые слова: ренормализационная группа, критическое поведение.

Для цитирования: Лебедев Н. М., Компаниец М. В. Критическое поведение 0(га)-сим-метричной модели с антисимметричным тензорным параметром порядка: четырёхпетлевое приближение // Вестник СПбГУ. Физика и химия. 2018. Т. 5 (63). Вып. 1. С. 20-31. https://doi.org/10.21638/11701/spbu04.2018.103

N. M. Lebedev, M. V. Kompaniets

CRITICAL BEHAVIOUR OF A O(n)-SYMMETRIO MODEL WITH ANTISYMMETRIC TENSOR ORDER PARAMETER: FOUR-LOOP APPROXIMATION

St. Petersburg State University,

7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation

The critical behavior of the 0(ra)-symmetric model with an antisymmetric tensor order parameter is studied by means of the renormalization group approach. Previously, this model was studied in the framework of the three-loop approximation by means of e-expansion with con-formal Borel resummation technique and pseudo-e-expansion. It was shown that the results of these approaches are in qualitative agreement with each other, however, the numerical values of the critical indices can differ quite significantly. Moreover, the question of the sensitivity of even qualitative results to the accounting of the higher orders of perturbation theory remained open. In the present paper, the RG functions and critical indices of the model are calculated up to four-loop accuracy within both approaches. Their analysis shows that the results of applying different approaches are qualitatively consistent both with each other and with the results of the three-loop analysis. Refs 29. Tables 3.

Keywords: renormalization group, critical behavior.

For citation: Lebedev N. M., Kompaniets M. V. Critical behaviour of a 0(ra)-symmetric model with antisymmetric tensor order parameter: Four-loop approximation. Vestnik SPbSU.

Physics and Chemistry. 2018. Vol. 5 (63), iss. 1. P. 20-31. https://doi.org/10.21638/11701/spbu04.2018.103

* Н. М. Лебедев выражает благодарность РФФИ за финансовую поддержку в рамках гранта 16-32-00086 мол-а.

© Санкт-Петербургский государственный университет, 2018

Введение. В работах [1-3] метод теоретико-полевой ренормгруппы применялся к 0(п)-симметричной модели ф4 с вещественным антисимметричным тензорным параметром порядка. Такая модель может иметь отношение к переходам между немати-ческой, холестерической и голубой фазами жидких кристаллов [4-9], переходу в дисси-метричную фазу в ферроэластиках [10-14] и (в комплексной модификации) к переходу в сверхпроводящее состояние в системах фермионов с высшими спинами или дополнительными степенями свободы [15-17]. При п = 2 и 3 модель соответственно сводится к хорошо изученным однозарядным скалярной и 0(3)-векторной ф4-моделям; при п ^ 4 это истинно двухзарядная модель.

В работе [1] ренормгрупповой (РГ) анализ исследуемой модели был выполнен в рамках однопетлевого приближения посредством метода е-разложения. Было показано, что только при п = 4 в модели имеются фиксированные точки с вещественными координатами: одна инфракрасно-устойчивая (ИК-устойчивая) и три неустойчивых. При п > 4 нетривиальные фиксированные точки в модели отсутствуют, а РГ-потоки уходят либо на бесконечность, либо в нефизическую область.

Однако в работе [2] было продемонстрировано, что данные результаты неустойчивы по отношению к учёту следующих порядков теории возмущений. В частности, непосредственный учёт трёхпетлевых поправок привёл к смене ИК-устойчивой в главном приближении фиксированной точки на седловидную, в то время как одна из седловидных точек, наоборот, стала ИК-устойчивой после учёта асимптотического характера е-разложений. В итоге в работе [2] был сделан вывод, что для достоверного определения типа перехода необходимо исследовать не только устойчивость полученных результатов по отношению к учёту следующих порядков теории возмущений, но также и их зависимость от выбора конкретной схемы ренормировки.

С этой целью в работе [3] РГ-функции и критические индексы модели были вычислены в рамках подхода ренормгруппы в фиксированной размерности пространства (в реальном пространстве) и метода псевдо-е-разложения (т-разложения) с трёхпет-левой точностью. Анализ полученных выражений позволил подтвердить качественные выводы работы [2], однако оставил открытым вопрос об их устойчивости по отношению к учёту следующих порядков теории возмущений.

В настоящей работе РГ-функции и критические индексы модели вычислены с че-тырёхпетлевой точностью в рамках обоих ранее исследованных подходов.

Известно, что т-разложение является эффективным методом изучения характеристик критического поведения различных систем. Более того, данный подход позволяет получать надёжные численные значения критических индексов путём непосредственного суммирования т-разложений без применения сложных техник пересуммирования, чувствительных к особенностям асимптотического поведения рядов теории возмущений

В то же время ряды теории возмущений, полученные в рамках е-разложения, демонстрируют существенно более быстрый рост коэффициентов. Поэтому в рамках данного подхода численные значения критических индексов были получены путём процедуры пересуммирования конформ-борелевским методом, учитывающей параметры ассимп-тотики е-разложения, вычисленные в работе [2].

Модель. В данной работе рассматривается модель вещественного антисимметричного тензорного поля ф = ф^(х) второго ранга (ф^ = -фы, = 1,...,п) в ¿-мерном евклидовом векторном пространстве. Модель задаётся функционалом действия:

[18-24].

6'(Ф) = ^ МФ(-32 + ш2)ф) - ^ (МФ2))2 - ^ МФ4).

1

(1)

4!

4!

Здесь и далее интегрирование по ¿-мерному пространству, а также суммирование по повторяющимся индексам подразумевается; д2 обозначает оператор Лапласа; т2 — отклонение температуры (или её аналога) от критического значения; дю, д20 — константы взаимодействия.

Для обеспечения стабильности модели необходимо, чтобы неквадратичная по полям часть функционала (1) была отрицательно определённой на любой конфигурации поля ф. Данное требование приводит к ограничениям на значения, которые могут принимать константы взаимодействия в случае чётных и нечётных значений п соответственно [1]:

2д10 + 920 > 0, пдю + д20 > 0; (2)

2дю + 920 > 0, (п - 1)дю + 920 > 0. (3)

Стандартный анализ симметрий функционала (1) и канонических размерностей полей и параметров модели показывает, что поверхностные ультрафиолетовые (УФ) расходимости, требующие устранения путём перенормировки, содержатся только в двух-и четырёхточечных 1-неприводимых функциях Грина. Как следствие, модель является мультипликативно ренормируемой с ренормированным действием:

БпШ = - Мф(-^52 + г2т2)ф) - ^з(Мф2))2 - ^4Мф4). (4)

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

7 - 7 - 7 - 71/2 (к\

^91 — 92 — ТРГ? — '

Уравнение ренормгруппы получается стандартным способом на основе произвольности ренормировочной массы и для к-точечной 1-неприводимой функции Грина имеет вид [25, 26]

(|д| + Р1д31 + Р2д32 - ут2Вт2 - куф)Г^ = 0. (6)

При этом РГ-функции задаются стандартными определениями [25, 26]:

Ь = В11п в = Вг = 1, 2, (7)

где -О^ = хдх при фиксированных затравочных параметрах.

Инфракрасное асимптотическое поведение функций Грина определяется набором ИК-притягивающих (ИК-устойчивых) неподвижных точек соответствующих уравнений РГ, координаты которых ищутся из условия одновременного зануления всех в-функций:

&(д*)=0. (8)

Тип неподвижной точки определяется при помощи матрицы

ЫгН = двг/дди\д=д, . (9)

Неподвижная точка является ИК-притягивающей, если матрица ю положительна в этой точке, т. е. положительна вещественная часть всех её собственных чисел. Существование ИК-притягивающей неподвижной точки предполагает скейлинговое поведение для всех функций Грина, описываемое критическими индексами, в частности индексом Фишера:

П = 2уф(д1*,д2*). (10)

В работе используются две различные схемы ренормировки: схема минимальных вычитаний MS и схема ренормировки в фиксированной размерности пространства, аналогичная используемой в работе [27].

В рамках обоих подходов константы ренормировки были явно вычислены в форме рядов по степеням gi, g2 с трёхпетлевой точностью. Каждую диаграмму исследуемой модели, дающую вклад в Zi, можно представить в виде произведения соответствующей диаграммы скалярной ф4 модели и дополнительного тензорного множителя, возникающего как результат свёртки тензорных частей пропагатора и вершинных множителей (подробное объяснение см. в работах [1, 2]). В случае подхода фиксированной размерности пространства значения диаграмм скалярной модели были взяты в работе [27] для d =2 ив работе [19] для d = 3. В случае схемы минимальных вычитаний использовались значения, взятые в работе [26]. Вычисление тензорных свёрток было произведено с помощью программы FORM [28, 29].

Разложения для РГ-функций (7) были получены по известным рядам для констант ренормировки (подробное объяснение см. в работах [2, 3]). Полученные выражения приводятся и обсуждаются в нижеследующих разделах.

Стоит также отметить, что модель с дополнительными соотношениями g2 =0 для любого n или соотношением 4gi + 3g2 = 0 для n = 4 является замкнутой по отношению к процедуре ренормировки. Следствием этого является тот факт, что для координат неподвижных точек, обозначаемых в дальнейшем как A и C, соответствующие соотношения g2* =0 и 4gi* + 3g2* = 0 должны выполняться точно, а по отношению к соответствующим однозарядным моделям эти точки должны быть ИК-притягивающими. Данный факт можно использовать для дополнительной проверки полученных в данной работе результатов.

е-разложение. В рамках подхода е-разложения для всех РГ-функций были получены полностью аналитические выражения, однако в силу их громоздкости и того факта, что после применения процедуры пересуммирования значения критических индексов могут быть найдены лишь численно, здесь мы тоже приводим только численные значения коэффициентов разложения. Для краткости мы приводим впервые вычисленные коэффициенты четырёхпетлевого приближения, в то время как ранее известные с трёхпетлевой точностью члены разложения обозначены как pi3)(gi, g2),

Р2; (gi

,g2), уф (gi,g2) и могут быть найдены в работе [2] по приведённым в ней ф°рму-лам (22)-(24) соответственно. Также для удобства в этом разделе всюду сделана замена gi,2 ^ gi,2/(8n2).

р1 =в13) (g1; g2) + (0,000010n6 - 0,000030n5 - 0,0643234n4 + 0,128698n3 - 2,144418n2 + + 2,080064n - 12,557898)g5 + (-0,000246n5 + 0,000614n4 - 0,623339n3 + + 0,934395n2 - 14,215095n + 6,951835)g4g2 + (-0,012994n4 + 0,025988n3 -

- 3,180664n2 + 3,167670n - 18,497684)g3g2 + (-0,000042n5 + 0,000106n4 -

- 0,125824n3 + 0,188630n2 - 6,142543n + 3,039837)g2g| + (-0,000969n4 + + 0,001937n3 - 0,339387n2 + 0,338281n - 2,275556)g1g4 + (-0,002666n3 +

+ 0,003963n2 - 0,151275n + 0,074705)g5 + O(g6), (11)

e2 =e23)(g1,g2) + (0,000173n6 - 0,000519n5 - 0,009335n4 + 0,019536n3 - 2,212644n2 + + 2,202790n - 34,982149)g4g2 + (0,001444n5 - 0,003609n4 - 0,239656n3 +

+ 0,363094«2 - 23,046496«. + 11,462612)^3^2 + (0,003395«4 - 0,006790«3 -

- 2,713640«2 + 2,717035« - 18,023949)^2^3 + (-0,058531«3 + 0,087865«2 -

- 3,043559«. + 1,507386)^1^4 + (-0,000168«4 + 0,000317«3 - 0,069960«2 +

+ 0,070323« - 0,454975)д5 + 0(д6), (12)

УФ =Уф3)(д1, д2) + (0,000007«6 + 0,000023«5 + 0,000219«4 - 0,000475«3 + 0,004340«2 -

- 0,004099« + 0,012056)д4 + (-0,000060«5 + 0,000151«4 + 0,002050«3 -

- 0,003225«2 + 0,025198« - 0,012056)д3д2 + (-0,000090«4 +

+ 0,000181«3 + 0,009087«2 - 0,009178« + 0,013744)д2д2 + (0,000422«3 -

- 0,000633«2 + 0,006420« - 0,003105)д1д| + (0,000002«4 - 0,000004«3 +

+ 0,000339«2 - 0,000337« + 0,000640)д4 + 0(д5). (13)

По известным разложениям РГ-функций координаты трёх нетривиальных критических точек, обнаруженных авторами работы [1], и соответствующие им критические индексы были вычислены в форме е-разложения. В приведённых ниже выражениях сделана подстановка « = 4, поскольку в работе [1] установлено, что нетривиальные неподвижные точки с вещественными координатами могут существовать только в этом случае.

Точка А:

д1 д2 ю1 Ю2 П

-0,178258е'

4

0,428571е + 0,209912е2 - 0,100507е3 0,

е - 0,489795е2 + 0,948832е3 - 2,521871е4, -0,142858е - 0,151603е2 + 0,117489е3 - 0,331731е4, 0,020408е2 + 0,014889е3 - 0,006576е4.

(14)

Точка В:

д1 д2 ю1 Ю2 П

0,705882е +0,818237е2 - 0,310641е3 + 0,626898е4, -0,705882е - 1,69021е2 + 0,205344е3 - 1,294462е4, е - 0,541522е2 + 1,019739е3 - 2,619199е4,

0,058823е - 0,106554е2 - 0,195408е3 0,020761е2 + 0,017295е3 — 0,003487е4.

0,040578е'

4

(15)

Точка С:

д1 д2 ю1 Ю2 П

0,818182е + 0,466567е2 - 0,258762е3 + 0,484615е4, — 1,09091е — 0,622089е2 + 0,345019е3 - 0,646153е4,

е - 0,570248е2 + 1,28289е3 - 3,78111е4, —0,090909е + 0,311796е2 - 0,230143е3 +

0,587198е'

4

0,020661е

2

0,018398е3 — 0,00744947е4.

(16)

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

(3)

Уф (д1,д2), могут быть найдены по приведённым в работе [3] формулам (15)-(17) соответственно. Кроме того, этом разделе так же, как и в предыдущем, всюду сделана замена д1;2 ^ д1г2/(8л2):

р1 =в3(д1, д2) + (0,000015п6 - 0,000045п5 - 0,006865п4 + 0,013805п3 - 0,239365п2 +

+ 0,232455п - 1,413744)д5 + (0,000261п5 - 0,000652п4 - 0,069723п3 + 0,105237п2 -

- 1,627854п + 0,796365)д4д2 + (-0,000349п4 + 0,000697п3 - 0,374186п2 + + 0,373837п - 2,099976)д3д| + (-0,000021п5 + 0,000052п4 - 0,013293п3 + + 0,019888п2 - 0,707143п + +0,350259)д2д| + (-0,000128п4 + 0,000257п3 -

- 0,036161п2 + 0,036070п - 0,264548)д 1 д4 + +(-0,000232п3 + 0,000351п2 -

- 0,016355п + 0,008115)д5 + 0(д6), (17)

р2 =в23)(д1, д2) + (-0,000055п6 + 0,000166п5 - 0,000115п4 - 0,000048п3 - 0,225369п2 + + 0,225421п - 3,883269)д4д2 + (-0,000299п5 + 0,000748п4 - 0,015290п3 + + 0,022187п2 - 2,531959п + 1,262307)д3д| + (-0,000411п4 + 0,000823п3 -

- 0,283863п2 + 0,283451п - 2,015581)д2д| + (-0,007615п3 + 0,011404п2 -

- 0,333321п + 0,164692)д1д4 + (-0,000024п4 + 0,000046п3 - 0,008252п2 +

+ 0,008230п - 0,052477)д5 + 0(д6), (18)

УФ =Уф3)(д1,д2) + (-0,00000021п6 + 0,00000062п5 + 0,00006381п4 - 0,00012861п3+

+ 0,00095355п2 - 0,00088912п + 0,00251239)д4 + (-0,00000166п5 + 0,00000415п4+ + 0,00051874п3 - 0,00078225п2 + 0,00528580п - 0,00251239)д3д2 + (0,00001333п4-

- 0,00002665п3 + 0,00199737п2 - 0,00198404п + 0,00279865)д2д| + (0,00012720п3-

- 0,00019080п2 + 0,00132477п - 0,00063058)д1д| + (0,00000121п4 - 0,00000242п3+ + 0,00007863п2 - 0,00007742п + 0,00012448)д4 + 0(д5). (19)

По известным разложениям РГ-функций координаты неподвижных точек и соответствующие им критические индексы были вычислены в форме т-разложения. Для этого в (17), (18) заменён линейный член д12 ^ тд^. В приведённых ниже выражениях сделана подстановка п = 4.

Точка А:

д1 = 0,428571т + 0,141237т2 +0,002840т3 + 0,002765т4,

д2 Ю1

Ю2 П

2

0,203958т3 — 0,147240т4,

т - 0,329554т -0,142857т - 0,101501т2 - 0,001258т3 - 0,023151т4, 0,012094т2 + 0,008979т3 + 0,003961т4.

Точка В:

д1 д2 ю1 Ю2 П

0,705882т + 0,547663т2 + 0,127564т3 - 0,045456т4, -0,705882т - 1,128978т2 - 0,512302т3 - 0,049290т4, т - 0,364091т2 + 0,183818т3 - 0,118377т4, 0,058823т — 0,070855т2 — 0,126034т3 - 0,064476т4,

0,012303т2 + 0,009984т3 + 0,005667т4.

3

4

(21)

Точка С:

д1 д2 ю1 Ю2 П

= 0,818182т + 0,313549т2 + 0,008818т3 + 0,007837т4,

= -1,09091т - 0,418065т2 - 0,011757т3 - 0,01045т4,

= т - 0,383226т2 + 0,272169т3 - 0,224950т4,

= -0,090909т + 0,207585т2 + 0,009024т3 + 0,031832т4,

= 0,012244т2 + 0,010404т3 + 0,005027т4.

(22)

Случай й = 2. Для подхода фиксированной размерности пространства и случая с! = 2 были получены следующие выражения для РГ-функций. Как и ранее, мы приводим лишь впервые вычисленные коэффициенты четырёхпетлевого приближения, а трёхпетлевые разложения, обозначенные как в13)(д1,д2), в23)(д1, д2), уф3)(д1,д2), могут быть найдены по приведённым в работе [3] формулам (21)-(23) соответственно. Для удобства в этом разделе всюду сделана замена д1;2 ^ д1,2/(2п):

— =ЫЕ11£з1 (—0 000009??6 + 0.000026??5 - 0.034690??4 + 0,069336??3 - 1,042054??2 + 2 2

+ 1,007390« - 5,857337)д5 + (0,000073«5 - 0,000183«4 - 0,336953«3 + + 0,505613«2 - 6,816885« + 3,324168)д4д2 + (-0,007113«4 + 0,014225«3 -

- 1,655153«2 + 1,648041« - 8,572945)д3д| + (-0,000075«5 + 0,000188«4 -

- 0,070248«3 + 0,105184«2 - 2,971747« + 1,468349)д2д| + (-0,000667«4 + + 0,001333«3 - 0,167944«2 + 0,167304« - 1,053886)д1д4 + (-0,001393«3 +

+ 0,002092«2 - 0,069839« + 0,034566)д5 + 0(д6), (23)

^ =^2 + (—0,000080??6 + 0,000239??5 - 0,005155т?4 + 0,009912??3 - 1,043498т?2 +

+ 1,038582« - 15,990014)д4д2 + (-0,000313«5 + 0,000783«4 - 0,120510«3 + + 0,179981«2 - -10,6584217« + 5,299240)д3д| + (-0,000952«4 + 0,001904«3 -

- 1,314103«2 + 1,313151« - 8,165407)д2д3 + (-0,034892«3 + 0,052325«2 -

- 1,430708«. + 0,706584)^1^4 + (-0,000137«4 + 0,000272«3 - 0,036431«2 + + 0,036299«. - 0,211757)^5 + 0(д6), (24)

УФ =уф3)(д1 ,д2) + (—0,00000447«6 + 0,00001341«5 + 0,00037783«4 - 0,00077802«3 +

+ 0,00603383«2 - 0,00564259« + 0,01602439)д4 + (-0,00003576«5 + 0,00008939«4 + + 0,00320146«3 - 0,00489159«2 + 0,033685267« - 0,016024389)д3д2 + + (0,00001941«4 - -0,00003882«3 + 0,01267427«2 - 0,01265486« + + 0,01787080)д2д2 + (0,00074673«3 - 0,00112009)«2 + 0,00849283« -

- 0,00405973)д1д| + (0,00000706«4 - 0,00001412«3 + 0,00048436«2 -

- 0,00047730« + 0,00081960)д| + 0(д5). (25)

Как и в предыдущем разделе, по известным разложениям РГ-функций находим координаты трёх нетривиальных неподвижных точек и соответствующие им критические индексы в форме т-разложения для единственного нетривиального случая « = 4. Точка А:

д1

д2 ю1

Т

Ш2

Т

П

0,428571т + 0,239836т2 + 0,018647т3 + 0,008478т4,

0,

- 0,559618т2 + 0,539324т3 — 0,589932т4,

= -0,142857т - 0,171799т2 - 0,009016т3 - 0,085776т4, = 0,037432т2 + 0,039666т3 + 0,020168т4.

(26)

Точка В:

д1

д2 ю1

Т

Ш2

Т п

0,705882т + 0,926789т2 + 0,423377т3 - 0,222441т4, -0,705882т - 1,907940т2 - 1,580563т3 - 0,213826т4, т - 0,617968т2 + 0,448858т3 - 0,517774т4,

0,058823т - 0,119387т2 - 0,362933т3 - 0,311133т4, 0,038080т2 + 0,044797т3 + 0,034695т4.

(27)

Точка С:

д1

д2 ю1

Т

Ш2

Т п

0,818182т + 0,532022т2 + 0,066832т3 + 0,045014т4, -1,09091т - 0,709362т2 - 0,089109т3 - 0,060018т4, т - 0,650249т2 + 0,682279т3 - 0,893010т4,

-0,090909т + 0,350054т2 + 0,034867т3 + 0,117126т4, 0,037896т2 + 0,047027т3 + 0,028966т4.

(28)

Обсуждение результатов. По известным т- и е-разложениям координат неподвижных точек и критических индексов можно получить их численные значения. Для

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

Полученные таким путём численные значения индексов приведены в табл. 1-3.

Таблица 1

Сравнение критических индексов в размерности й = 2 для неподвижной точки Л при п = 4

Индекс й = = 3 й = = 2

е-разложение т-разложение е-разложение т-разложение

91 0,542 0,575 1,102 0,695

92 0 0 0 0

С01 0,781 0,727 1,347 0,780

С02 -0,226 -0,269 -0,487 -0,819

л 0,018 0,025 0,052 0,097

Таблица 2

Сравнение критических индексов в размерности й = 2 для неподвижной точки В при п = 4

Индекс й = = 3 й = = 2

е-разложение т-разложение е-разложение т-разложение

91 1,145 1,336 2,401 1,834

92 -1,747 -2,396 -3,954 -4,408

С01 0,755 0,701 1,277 0,626

С02 -0,055 -0,203 -0,198 -1,469

Л 0,017 0,028 0,045 0,118

Таблица 3

Сравнение критических индексов в размерности й = 2 для неподвижной точки С при п = 4

Индекс й = = 3 й = = 2

е-разложение т-разложение е-разложение т-разложение

91 1,027 1,148 2,019 1,462

92 -1,369 -1,531 -2,692 -1,949

С01 0,762 0,664 1,298 0,278

С02 0,077 0,158 0,245 0,822

Л 0,017 0,028 0,044 0,114

Из табл. 1-3 видно, что для всех трёх неподвижных точек результаты суммирования т-разложений качественно согласуются с результатами обработки е-разложений. Более того, они так же качественно согласуются с результатами, полученными ранее в рамках трёхпетлевого анализа [3].

Тем не менее значения критических индексов могут довольно существенно различаться в зависимости от выбора схемы ренормировки при переходе к следующему по-

рядку теории возмущений в рамках одной схемы. Особенно сильно различаются результаты в случае d =2. Более того, из табл. 2 видно, что результаты непосредственного суммирования т-разложений координат неподвижной точки B приводят к тому, что для неё перестают выполняться неравенства (2), (3), ив результате она покидает границу физической области параметров. Возможным объяснением данного факта может служить более сильная расходимость рядов теории возмущений в случае d =2, которая приводит к необходимости учитывать асимптотические свойства т-разложений уже на уровне трёх-четырёх петель. Этим же можно объяснить и более существенное расхождение значений индексов в сравнении со случаем d =3.

Таким образом, можно заключить, что предсказания относительно возможных типов критического поведения изучаемой модели, сделанные в работе [2], оказались устойчивыми не только по отношению к выбору схемы ренормировки, но и по отношению к учёту следующих порядков теории возмущений. Это, в свою очередь, позволяет надеяться, что учёт следующих порядков теории возмущений, вкупе с нахождением параметров асимптотики т-разложений и разработкой учитывающей их процедуры пересуммирования, позволит получить надёжные значения критических индексов модели.

Литература

1. AntonovN. V., Kompaniets M. V., LebedevN.M. Critical behaviour of the 0(га)-ф4 model with an antisymmetric tensor order parameter //J. Phys. (A). 2013 Vol. 46. 405002.

2. AntonovN. V., Kompaniets M. V., LebedevN.M. Critical behavior of the O(ra)^4-model with an antisymmetric tensor order parameter: Three-loop approximation // Theor. Math. Phys. 2017. Vol. 190, N 2. P. 204-216.

3. Лебедев Н.М., Компаниец М. В. Критическое поведение 0(га)-симметричной модели с антисимметричным тензорным параметром порядка: ренормгруппа в реальном пространстве // Вестник СПбГУ. Физика и химия. 2017. Т. 4 (62), вып. 4. C. 417-428.

4. BrazovskiS. A., DmitrievD. Phase transitions in cholesteric liquid crystals // Sov. Phys. JETP. 1975. Vol. 42, N 3. P. 497.

5. BrazovskiS. A., FilevV.M. Critical phenomena in cholesteric liquid crystals // Sov. Phys. JETP. 1978 Vol. 48, N 3. P. 497-573.

6. GrebelH., Hornreich R. M., Shtrikman S. Landau theory of cholesteric blue phases // Phys. Rev. (A).

1983. Vol. 28. P. 1114.

7. GrebelH., Hornreich R. M., Shtrikman S. Landau theory of cholesteric blue phases: The role of higher harmonics // Phys. Rev. (A). 1984. Vol. 30. P. 3264.

8. Hornreich R. M., Shtrikman S. Some open questions in cholesteric blue phases // Z. Phys. (B): Cond. Matt. 1987. Vol. 68. P. 369.

9. Belyakov V. A., Dmitrienko V. E. The blue phase of liquid crystals // Sov. Phys. Uspekhi. 1985. Vol. 28. P. 535.

10. Изюмов Ю. А., Сыромятников В. Н. Фазовые переходы и симметрия кристаллов. М.: Наука,

1984.

11. DoveM. T., Redfern S. A. T. Lattice simulation studies of the ferroelastic phase transitions in (Na,K)AlSi3Os and (Sr,Ca)Al2Si2Os feldspar solid solutions // Am. Mineral. 1997. Vol. 82. P. 8.

12. Watson G. W., Parker S. C. Dynamical instabilities in a-quartz and a-berlinite: A mechanism for amorphization // Phys. Rev. (B). 1995. Vol. 52. P. 13306.

13. Goryainov S. V., Ovsyuk N. N. Twisting of a-quartz tetrahedra at pressures near the transition to the amorphous state //J. Exp. Theor. Phys. Lett. 1999. Vol. 69. P. 467.

14. Goryainov S. V., Ovsyuk N.N. Mechanism of the formation of a soft mode in ferroelastic phase transition // J. Exp. Theor. Phys. Lett. 2001. Vol. 73, iss. 8. P. 408-410.

15. NalimovM. Yu., Komarova M. V., Honkonen J. Temperature Green's functions in Fermi systems: The superconducting phase transition // Theor. Math. Phys. 2013. Vol. 176, N 1. P. 89-97.

16. Kalagov G. A., Kompaniets M. V., NalimovM. Yu. Renormalization-group investigation of a superconducting U(r)-phase transition using five loops calculations // Nucl. Phys. (B). 2016. Vol. 905. P. 16-44.

17. Kalagov G. A., Kompaniets M. V., NalimovM. Yu. Renormalization-group study of a superconducting phase transition: Asymptotic behavior of higher expansion orders and results of three-loop calculations // Theor. Math. Phys. 2014. Vol. 181, N 2. P. 1448-1458.

18. Orlov E. V., Sokolov A. I. Critical thermodynamics of two-dimensional systems in the vi-loop renormalization-group approximation // Phys. Solid State. 2000. Vol. 42, N 11. P. 2151-2158.

19. Nickel B. G., Meiron D. I., Baker G. A. Compilation of 2-pt and 4-pt graphs for continuous spin model // University of Guelph rep. 1977.

20. Baker G. A., Nickel B. G., Meiron D. I. Critical indices from perturbation analysis of the Callan — Sy-manzik equation // Phys. Rev. (B). 1978. Vol. 17. P. 1365-1374.

21. Le GuillouJ. C., Zinn-Justin J. Critical exponents from field theory // Phys. Rev. (B). 1980. Vol. 21, N 9. P. 3976-3998.

22. Guida R., Zinn-Justin J. Critical exponents of the N-vector model //J. Phys. (A). 1998. Vol. 31. P. 8103.

23. FolkR., Holovatch Yu., Yavorskii T. Pseudo-e expansion of six-loop renormalization-group functions of an anisotropic cubic model // Phys. Rev. (B). 2000. Vol. 62. P. 12195.

24. Holovatch Yu., Ivaneiko D., Delamotte B. On the criticality of frustrated spin systems with non-collinear order // J. Phys. (A). 2004. Vol. 37. P. 3569.

25. Zinn-Justin J. Quantum field theory and critical phenomena. Oxford: Clarendon, 1989.

26. VasilievA.N. The field theoretic renormalization group in critical behavior theory and stochastic dynamics. Boca Raton: Chapman & Hall/CRC, 2004.

27. AdzhemyanL. Ts., Kirienko Yu. V., Kompaniets M. V. Critical exponent n in 2D O(N)-symmetric ^4-model up to 6 loops // arXiv: 1602.02324.

28. Vermaseren J. A. M. New features of FORM // arXiv: math-ph/0010025.

29. KuipersJ., UedaT., Vermaseren J. A. M., Vollinga J. FORM version 4.0 // Comp. Phys. Comm. 2013. Vol. 184, N 5. P. 1453-1467.

References

1. Antonov N. V., Kompaniets M.V., Lebedev N. M. Critical behaviour of the O(ra)-^4 model with an antisymmetric tensor order parameter. J. Phys. (A), 2013, vol. 46, 405002.

2. Antonov N. V., Kompaniets M.V., Lebedev N.M. Critical behavior of the O(ra)-^4-model with an antisymmetric tensor order parameter: Three-loop approximation. Theor. Math. Phys., 2017, vol. 190, no 2, pp. 204-216.

3. Lebedev N.M., Kompaniets M. V. Kriticheskoe povedenie O(n)-simmetrichnoi modeli s antisimmet-richnym tenzornym parametrom poriadka: renormgruppa v real'nom prostranstve [Critical behaviour of a O(n)-symmetric model with an antisymmetric tensor order parameter: the real-space renormalization group]. Vestnik SPbSU. Physics and Chemistry, 2017, vol. 4 (62), iss. 4, pp. 417-428. (In Russian)

4. Brazovski S. A., Dmitriev D. Phase transitions in cholesteric liquid crystals. Sov. Phys. JETP.,1975, vol. 42, no 3, pp. 497.

5. Brazovski S.A., Filev V. M. Critical phenomena in cholesteric liquid crystals. Sov. Phys. JETP., 1978 vol. 48, no 3, pp. 497-573.

6. Grebel H., Hornreich R. M., Shtrikman S. Landau theory of cholesteric blue phases. Phys. Rev. (A), 1983, vol. 28, pp. 1114.

7. Grebel H., Hornreich R. M., Shtrikman S. Landau theory of cholesteric blue phases: The role of higher harmonics. Phys. Rev. (A), 1984, vol. 30, pp. 3264.

8. Hornreich R. M., Shtrikman S. Some open questions in cholesteric blue phases. Z. Phys. (B): Cond. Matt., 1987, vol. 68, pp. 369.

9. Belyakov V. A., Dmitrienko V. E. The blue phase of liquid crystals. Sov. Phys. Uspekhi, 1985, vol. 28, pp. 535.

10. Iziumov Iu.A., Syromiatnikov V.N. Fazovye perekhody i simmetriia kristallov [Phase transitions and symmetry of crystals]. Moscow, Nauka Publ., 1984. (In Russian)

11. Dove M. T., Redfern S. A. T. Lattice simulation studies of the ferroelastic phase transitions in (Na,K)AlSi3Os and (Sr,Ca)Al2Si2Os feldspar solid solutions. Am. Mineral, 1997, vol. 82, pp. 8.

12. Watson G.W., Parker S. C. Dynamical instabilities in a-quartz and a-berlinite: A mechanism for amorphization. Phys. Rev. (B), 1995, vol. 52, pp. 13306.

13. Goryainov S. V., Ovsyuk N. N. Twisting of a-quartz tetrahedra at pressures near the transition to the amorphous state. J. Exp. Theor. Phys. Lett., 1999, vol. 69, pp. 467.

14. Goryainov S. V., Ovsyuk N. N. Mechanism of the formation of a soft mode in ferroelastic phase transition. J. Exp. Theor. Phys. Lett., 2001, vol. 73, iss. 8, pp. 408-410.

15. Nalimov M. Yu., Komarova M. V., Honkonen J. Temperature Green's functions in Fermi systems: The superconducting phase transition. Theor. Math. Phys., 2013, vol. 176, no 1, pp. 89-97.

16. Kalagov G. A., Kompaniets M. V., Nalimov M. Yu. Renormalization-group investigation of a superconducting U(r)-phase transition using five loops calculations. Nucl. Phys. (B), 2016, vol. 905, pp. 16-44.

17. Kalagov G. A., Kompaniets M. V., Nalimov M. Yu. Renormalization-group study of a superconducting phase transition: Asymptotic behavior of higher expansion orders and results of three-loop calculations. Theor. Math. Phys., 2014, vol. 181, no 2, pp. 1448-1458.

18. Orlov E. V., Sokolov A. I. Critical thermodynamics of two-dimensional systems in the vi-loop renormalization-group approximation. Phys. Solid State, 2000, vol. 42, no 11, pp. 2151-2158.

19. Nickel B. G., Meiron D. I., Baker G. A. Compilation of 2-pt and 4-pt graphs for continuous spin model. University of Guelph rep., 1977.

20. Baker G. A., Nickel B. G., Meiron D. I. Critical indices from perturbation analysis of the Callan — Sy-manzik equation. Phys. Rev. (B), 1978, vol. 17, pp. 1365-1374.

21. Le Guillou J.C., Zinn-Justin J. Critical exponents from field theory. Phys. Rev. (B), 1980, vol. 21, no 9, pp. 3976-3998.

22. Guida R., Zinn-Justin J. Critical exponents of the N-vector model. J. Phys. (A), 1998, vol. 31, pp. 8103.

23. Folk R., Holovatch Yu., Yavorskii T. Pseudo-e expansion of six-loop renormalization-group functions of an anisotropic cubic model. Phys. Rev. (B), 2000, vol. 62, pp. 12195.

24. Holovatch Yu., Ivaneiko D., Delamotte B. On the criticality of frustrated spin systems with non-collinear order. J. Phys. (A), 2004, vol. 37, pp. 3569.

25. Zinn-Justin J. Quantum field theory and critical phenomena. Oxford, Clarendon, 1989.

26. Vasiliev A. N. The field theoretic renormalization group in critical behavior theory and stochastic dynamics. Boca Raton, Chapman & Hall/CRC, 2004.

27. Adzhemyan L. Ts., Kirienko Yu. V., Kompaniets M. V. Critical exponent n in 2D O(N)-symmetric 94-model up to 6 loops. arXiv: 1602.02324.

28. Vermaseren J. A.M. New features of FORM. arXiv : math-ph/0010025.

29. Kuipers J., Ueda T., Vermaseren J. A.M., Vollinga J. FORM version 4.0. Comp. Phys. Comm.., 2013, vol. 184, no 5, pp. 1453-1467.

Статья поступила в редакцию 4 декабря 2017 г.

Контактная информация

Лебедев Никита Михайлович — аспирант; e-mail: [email protected]

Компаниец Михаил Владимирович — доктор физико-математических наук, доцент; e-mail: [email protected]

Lebedev Nikita Mihailovich — post-graduate student; e-mail: [email protected]

Kompaniets Mikhail Vladimirivich — Doctor of Physics and Mathematics, Associate Professor; e-mail: [email protected]