CC BY
6
УДК 532.517.45
Вестник СПбГУ. Физика и химия. 2018. Т. 5 (63). Вып. 1
С. Е. Воробьёва, Э. В. Иванова, В. Д. Серов
БОРЕЛЕВСКОЕ ПЕРЕСУММИРОВАНИЕ ДИНАМИЧЕСКОГО ИНДЕКСА г В МОДЕЛИ А КРИТИЧЕСКОЙ ДИНАМИКИ С УЧЁТОМ АСИМПТОТИКИ СИЛЬНОЙ СВЯЗИ
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
В данной работе рассмотрена модель А критической динамики, описывающая динамические характеристики критического поведения ферромагнетиков, а также некоторых других систем с несохраняющимся параметром порядка. Для теоретического описания таких моделей метод ренормализационной группы и е-разложения являются традиционными инструментами для построения ряда теории возмущений для критических показателей и других универсальных величин. Одной из проблем в данном формализме является асимптотическая природа рядов, получаемых в рамках е-разложения, поэтому для получения надёжных теоретических предсказаний в рамках данного формализма необходимо проведение борелевского пересуммирования. В настоящей работе такое пересуммирование выполнено для четырёхпетлевого е-разложения динамического критического индекса z, описывающего критическое замедление в системах типа ферромагнетиков. Решение этой задачи было стимулировано двумя факторами. Во-первых, недавно получены уточнённые четырёхпетлевые значения для коэффициентов е-разложения. Во-вторых, для статических моделей было показано, что включение асимптотики сильной связи в борелевское пересуммирование приводит к лучшей сходимости процедуры. Также было показано, что для задач критической статики учёт асимптотики сильной связи существенно улучшает предсказания старших членов ряда. Поскольку асимптотика сильной связи для данной модели неизвестна, данный параметр является свободным и фиксируется из условия наиболее точного предсказания последнего известного члена по предыдущим и скорости сходимости процедуры. Библиогр. 11 назв. Ил. 3.
Ключевые слова: ренормализационная группа, е-разложение, многопетлевые диаграммы, критические показатели, борелевское пересуммирование.
Для цитирования: Воробьёва, С. Е., Ива,ново, Э. В., Серов В. Д. Борелевское пересуммирование динамического индекса z в модели А критической динамики с учётом асимптотики сильной связи // Вестник СПбГУ. Физика и химия. 2018. Т. 5 (63). Вып. 1. C. 13—19. https://doi.org/10.21638/11701/spbu04.2018.102
S. Y. Vorobyeva, E. V. Ivanova, V. D. Serov
BOREL SUMMATION OF THE DINAMIC INDEX г
IN THE MODEL A OF CRITICAL DINAMICS WITH AN ACCOUNT
OF THE STRONG COUPLING ASYMPTOTICS
St. Petersburg State University,
7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
In this paper we consider the model A of critical dynamics which describes the dynamic characteristics of the critical behavior of ferromagnets as well as other systems with a nonconserved order parameter. For such a model the renormalization group method and the e-expansion are the common tool for constructing consistent perturbative expansion for critical exponents and other universal quantities. One of the problems in this formalism is the asymptotic nature of the series obtained within the e-expansion, so to get reliable theoretical predictions it is necessary to apply the Borel resummation technique. In this work we apply such a resummation technique to the four-loop e-expansion of the dynamical exponent z which describes the critical slowing down
© Санкт-Петербургский государственный университет, 2018
in the ferromagnet-type systems. The motivation of the paper is twofold: 1) improved values for the coefficients of the e-expansion have recently been obtained, 2) sit has been shown that incorporation of strong coupling asymptotics into the Borel resummation technique leads to better convergence of the resummation procedure for static models. In this investigation it was shown that taking into account the asymptotics of the strong coupling for critical static problems also improves the predictions of the higher terms of the series. Since the strong-coupling asymptotic for this model is unknown, this parameter is held free and adjusted in accordance with the condition of the most accurate prediction of the last known term from the previous ones and the rate of convergence of the procedure. Refs 11. Figs 3.
Keywords: renormalization group, e-expansion, multi-loop diagrams, critical exponents, Borel resummation.
For citation: Vorobyeva S. Y., Ivanova E. V., SerovV. D. Borel summation of the dynamic index г in the model A of critical dynamics with an account of the strong coupling asymptotics.
Vestnik SPbSU. Physics and Chemistry. 2018. Vol. 5 (63), iss. 1. P. 13-19. https://doi.org/10.21638/11r01/spbu04.2018.102
Одной из основных черт непрерывных фазовых переходов и критических явлений является их универсальность. Так, модель А критической динамики, до последнего времени использовавшаяся для описания критического замедления в ферромагнетиках [1], в настоящее время используется также для описания релаксационных процессов при структурных фазовых переходах [2]. В настоящей работе в рамках ренормгруппового подхода и е-разложения произведено пересуммирование четырёхпетлевых результатов для динамического критического индекса z, который определяется динамическим индексом у* и статическим индексом Фишера п:
z = 2 + у* - (1)
Здесь у* = у(и*) и п = 2уф(и*) — значения ренормгрупповых функций у(и) и Уф(и) в неподвижной точке и = и* — корне уравнения в(и*) = 0 , где в(и) — так называемая бета-функция. Величины у(и) и в(и) вычисляются в виде рядов по константе связи и (мы используем нормировку и, как в статье [4], эта нормировка отличается в 7/2 раза от принятой в работе [5]). Бета-функция динамической А модели совпадает со статической и известна в настоящее время с шестипетлевой точностью [5]. Для однокомпонентного поля она имеет вид
/ зи 17и2
|3(и) = и (-£+ — - + 4,06871ы3 - 16,97536ы4 + 89,01775ы5 -
- 543,37705и6 + 0(и7)^ , (2)
где е = 4 — d — отклонение размерности пространства d от критического значения dc =4.
Величина п(и) также известна с шестипетлевой точностью [7]:
п(и) = 0,02083333и2 — 0,007812500и3 + 0,02115885и4 —
— 0,06017453и5 + 0,2247499и6 + 0(иг). (3)
Рассматривая е как формально малый параметр, решая итерационно уравнение в(и*) = 0 и подставляя результат в (3), получаем е-разложение величины п(е):
п(е) = 0,0185185е2 + 0,01869е3 — 0,00832875е4 + 0,0256564е5 — 0,0812726е6 + O(er). (4)
Для функции y(g) максимальная достигнутая в настоящее время точность — четы-рёхпетлевое приближение — получена в работе [8]:
y(u) = C2U2 + С3М3 + с^и4. (5)
Погрешность расчёта коэффициента C4 в (5) в этой работе довольно велика — порядка 1%. Используя при вычислении фейнмановских диаграмм метод Sector Decomposition, нам удалось на два порядка увеличить точность расчётов. В результате были получены следующие значения коэффициентов в (5):
с2 = 0,07192051811,
с3 = -0,03552374(4), (6)
с4 = 0,0996502(26).
Соответствующий ряд для у*(е) имеет вид
у*(е) = 0,0319646747169е2 + 0,029726259(11)е3 - 0,0139079(5)е4 + 0(е5). (7)
Ряды (2)—(5), (7) являются асимптотическими и требуют пересуммирования для корректного определения значений критических показателей. На основе данных статьи [8] такое пересуммирование для показателя z было выполнено в работе [9], при этом сходимость итерационной процедуры оказалась плохой и значение индекса z было определено с большой погрешностью.
В настоящей работе мы используем для пересуммирования уточнённые значения (6) и модифицированный конформ-борелевский метод суммирования. Напомним его содержание. Пусть имеется функция A(g), для которой известны первые N коэффициентов An в разложении
^
A(g) = £ Angn, (8)
n=0
а также асимптотика высоких порядков (АВП) коэффициентов An:
An -► с (—a)nn! nb0. (9)
n—
Тогда конформ-борелевское представление функции A(g) имеет вид
А(д) = J™ dte-Hb (-^j B(w(gi)), (10)
где
цх)=(п)
у 1 + ах +1
В(х) = ^ Вп™п. (12)
п=0
Асимптотика высоких порядков (АВП) коэффициентов разложения функции (10) в ряд по д имеет вид (9) с заданным значением параметра а, заданное значение параметра 60 достигается путём выбора
3
Ъ = Ъ0 + -. (13)
Приближённое выражение )(д) для функции (10) получается путём учёта в (12) N первых членов суммы, коэффициенты Вп в которой выбираются из требования совпадения первых N членов разложения функции (10) с известными коэффициентами Ап в (8):
, XV N
„-иы
т(
,-Ж / . XV N
Вп
(14)
Параметр V в (10) контролирует асимптотику «сильной связи»:
А(д)
9У-
(15)
В обсуждаемой задаче параметр V неизвестен, и мы будем использовать его для улучшения сходимости процедуры суммирования при увеличении числа учитываемых слагаемых в (8).
Следуя работе [6], проиллюстрируем эффективность такого подхода на примере модельной функции
/>оо
I 2 4
А(д)= / йхв-х -дх . (16)
■ 'о
Для этой функции известны все коэффициенты разложения (8), а АВП имеет вид (9) с а = 4,Ъ0 = —1. Параметр V = —1/4. На рис. 1 показано сравнение точного значения А(1) = 0,6842134279 функции (16) при д =1 с результатом расчёта по формуле (10) с учётом различного числа слагаемых (горизонтальная ось) при различных значениях параметра V.
Рис. 1. Значение пересуммированного ряда А(м\д) функции (16) при д = 1 для различных N и V
На рис. 1 видно, что сходимость процедуры суммирования имеет место в некотором диапазоне значений параметра V, и скорость сходимости заметно уменьшается по мере удаления этого параметра от «истинного» значения V = —1/4.
Мы применили описанную процедуру для суммирования е-разложения динамического критического индекса г. Для этого из соотношений (1), (4) и (7) были найдены
Рис. 2. Пересуммированный критический индекс г при е =1 и различных N и V
2,025 2,020 ^2,015
II
С.0
и 2,010 2,005 2,000
V = 2,8
V = 2,7^^^^ *— ___к
X- ' •V = 2,6
<
V = 0
3
N
4
2,14 2,12 2,10
^ 2,08
II
2,06
N
2,04 2,02
Рис. 3. Пересуммированный крити- 2,00 ческий индекс г при е = 2 и различных N и V
четыре члена е-разложения величины г(е) = 2 + у(е) — п(е), и полученное разложение суммировалось согласно соотношениям (10)—(13). Параметры а и Ьо для п(е) выбирались, согласно [10], равными а = 1/3, Ьо = 3/2. В работе [11] показано, что такие же значения имеют эти параметры для величины -г(е). На рис. 2 показаны полученные в результате суммирования значения динамического критического индекса г для трёхмерного пространства с! = 3 (е =1) и двумерного пространства с! = 2 (е = 2) при различном числе учтённых членов разложения г(е) и значений параметра V.
Из сравнения рис. 1 и рис. 2 видно, что сходимость процедуры борелевского суммирования для г(е) вполне аналогична рассмотренной для модельной функции (16). Наибольшая скорость сходимости достигается при значении V = 2,6, что даёт следующие оценки для динамического индекса г для трёхмерного и двумерного случаев:
Число петель..... Две Три Четыре
е = 1 .............. 2,0210142 2,021218 2,0217974
е = 2.............. 2,1050275 2,10642878 2,1118359
Результаты работы показывают, что учёт параметра сильной связи V необходим для корректного определения динамического критического индекса г методом
V = 2,8 - - V =
X--- *
V = 2,6
< 1 1 Л 1
V =
V = 0:
конформ-борелевского суммирования, поскольку при значении v = 0, которое использовалось в работе [9], сходимость очень медленная и требуется учёт большого числа членов разложения z(e).
Исследования были проведены с использованием вычислительных ресурсов ресурсного центра «Вычислительный центр СПбГУ» (http://www.cc.spbu.ru/).
Литература
1. MarinelliM., MercuriF., BelangerD. P. Specific heat, thermal diffusivity, and thermal conductivity of FeF2 at the Neel temperature // Phys. Rev. (B). 1995. Vol. 51. P. 8897.
2. LivetF., FevreM., Beutier G., Sutton M. Ordering fluctuation dynamics in AuAgZn2 // Phys. Rev. (B). 2015. Vol. 92. 094102.
3. Антонов Н. В., Васильев А. Н. Критическая динамика как теория поля // Теор. мат. физика. 1984. Т. 60, № 1. C. 59-71.
4. АджемянЛ. Ц., Воробьёва, С. Е., КомпаниецМ. В. Представление несингулярными интегралами ß-функции и аномальных размерностей в моделях критической динамики // Теор. мат. физика. 2015. Т. 185, № 1. C. 3-11.
5. Kompaniets M. V., Panzer E. Minimally subtracted six loop renormalization of 0(ra)-symmetric ф4 theory and critical exponents // Phys. Rev. (D). 2017. Vol. 96. 036016.
6. Kompaniets M. V. Prediction of the higher-order terms based on Borel resummation with conformal mapping // J. Phys. Conf. Ser. 2016. Vol. 762. 012075.
7. Batkovich D. V., Chetyrkin K. G., Kompaniets M. V. Six loop analytical calculation of the field anomalous dimension and the critical exponent n in 0(ra)-symmetric ф4 model // Nucl. Phys. (B). 2016. Vol. 906. P. 147-167.
8. АджемянЛ. Ц., Новиков С. В., Сладкофф Л. Расчёт динамического индекса модели A критической динамики в порядке е4 // Вестник СПбГУ. Физика и химия. 2008. Вып. 4. C. 110-114.
9. Налимов М. Ю., Сергеев В. А., Сладкофф Л. Борелевское пересуммирование е-разложения динамического индекса г модели A ф4(0(га))-теории // Теор. мат. физика. 2009. Т. 159, № 1. C. 96-108.
10. Kleinert H., Schulte-Frohlinde V. Critical properties of ф4-theories. Freie Universität Berlin, 2001. 489 p.
11. Andreanov J., Honkonen J., Komarova M., Nalimov M. Large-order asymptotes for dynamic models // J. Phys. (A). 2006. Vol. 39. P. 1-10.
References
1. Marinelli M., Mercuri F., Belanger D. P. Specific heat, thermal diffusivity, and thermal conductivity of FeF2 at the Neel temperature. Phys. Rev. (B), 1995, vol. 51, pp. 8897.
2. Livet F., Fevre M., Beutier G., Sutton M. Ordering fluctuation dynamics in AuAgZn2. Phys. Rev. (B), 2015, vol. 92, 094102.
3. Antonov N. V., Vasil'ev A. N. Kriticheskaia dinamika kak teoriia polia [Critical dynamics as a field theory]. Teor. mat. fiziha. [Theor. Math. Phys.], 1984, vol. 60, no 1, pp. 59—71. (In Russian)
4. Adzhemyan L. Ts., Vorob'eva S. E., Kompaniets M.V. Predstavlenie nesinguliarnymi integralami ß-funktsii i anomal'nykh razmernostei v modeliakh kriticheskoi dinamiki [Representation of the ß-function and anomalous dimensions by nonsingular integrals in models of critical dynamics]. Teor. mat. fiziha. [Theor. Math. Phys.], 2015, vol. 185, no 1, pp. 3-11. (In Russian)
5. Kompaniets M. V., Panzer E. Minimally subtracted six loop renormalization of 0(ra)-symmetric ф4 theory and critical exponents. Phys. Rev. (D), 2017, vol. 96, 036016.
6. Kompaniets M.V. Prediction of the higher-order terms based on Borel resummation with conformal mapping. J. Phys. Conf. Ser., 2016, vol. 762, 012075.
7. Batkovich D. V., Chetyrkin K. G., Kompaniets M. V. Six loop analytical calculation of the field anomalous dimension and the critical exponent n in 0(ra)-symmetric ф4 model. Nucl. Phys. (B), 2016, vol. 906, pp. 147-167.
8. Adzhemyan L. Ts., Novikov S. V., Sladkoff L. Raschet dinamicheskogo indeksa modeli A kriticheskoi dinamiki v poriadke e4 [Calculation of the dynamical exponent in the model A of critical dynamics to order E4]. Vestnih SPbSU. Physics and Chemistry, 2008, iss. 4, pp. 110-114. (In Russian)
9. Nalimov M.Yu., Sergeev V. A., Sladkoff L. Borelevskoe peresummirovanie e-razlozheniia dinamich-eskogo indeksa z modeli A ф4(0(ra))-teorii [Borel resummation of the e-expansion of the dynamical exponent
г in model A of the ф4(0(га)) theory]. Teor. mat. fizika. [Theor. Math.. Phys.], 2009, vol. 159, no 1, pp. 96—108. (In Russian)
10. Kleinert H., Schulte-Frohlinde V. Critical properties of ф4-theories. Freie Universität Berlin, 2001. 489 p.
11. Andreanov J., Honkonen J., Komarova M., Nalimov M. Large-order asymptotes for dynamic models. J. Phys. (A), 2006, vol. 39, pp. 1-10.
Статья поступила в редакцию 20 декабря 2017 г.
Контактная информация
Воробьёва Светлана Евгеньевна — аспирантка; e-mail: [email protected] Иванова Элла Валерьевна — аспирантка; e-mail: [email protected] Серов Виталий Дмитриевич — студент; e-mail: [email protected]
Vorobyeva Svetlana Yevgenyevna — post-graduate student; e-mail: [email protected] Ivanova Ella Valeryevna — post-graduate student; e-mail: [email protected]
Serov Vitalii Dmitrievich, — student; e-mail: [email protected]
