ТЕОРЕТИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ ФИЗИКА
УДК 530.145
ПЕРЕНОРМИРОВКА ЭНЕРГИИ КАЗИМИРА С ЛОГАРИФМИЧЕСКИМИ РАСХОДИМОСТЯМИ
И.Ю. Малахов, К. А. Свешников, П. К. Силаев
(кафедра квантовой теории и физики высоких энергий) E-mail: [email protected]
Предложена процедура перенормировки вакуумной энергии поля для того случая, когда энергии уровней известны только численно. При этом не требуется ни явного вида трансцендентного уравнения для определения уровней, ни аналитического выражения для их асимптотики. Процедура дает однозначный ответ даже при наличии логарифмической расходимости в выражении для вакуумной энергии.
Введение
Еще в 1948 г. Казимир [1] нашел поправки к энергии макроскопической системы, обусловленные вакуумными флуктуациями квантованного электромагнитного поля. Эффект Казимира является наблюдаемым квантовым эффектом, играющим существенную роль как на микро-, так и на макроуровнях в различных областях физики [2, 3]. Между тем вычисление вакуумной энергии квантованного поля при наличии нетривиальных граничных условий (далее — энергии Казимира), за исключением простейших случаев свободных полей, заключенных в полостях с плоской границей, представляет собой весьма серьезную проблему.
В частности, простейшей системе — свободному полю в шаре, т. е. системе, в которой можно выписать явные трансцендентные уравнения для спектра, было посвящено достаточно большое количество работ: электромагнитное поле внутри проводящей сферы было исследовано в [4-6]; фермионы, заключенные внутри сферы с помощью условия запирания, принятого в моделях кварковых мешков — в [7-11]; эффекты, возникающие из-за ненулевой массы фермионов, были исследованы в [12, 13]. Тем не менее первые аналитические результаты в достаточно замкнутой форме были получены лишь в работе [14] для массивного скалярного поля и в [15] для фермионов. Следует подчеркнуть, что эти результаты существенным образом опираются на возможность выписать явное трансцендентное уравнение для спектра: при наличии явных уравнений для уровней возможен переход от сумм, содержащих неявные выражения, к интегралам, которые содержат только явные выражения [16].
Одной из наиболее принципиальных трудностей при перенормировке энергии Казимира является неоднозначность, связанная с логарифмическими расходимостями. Действительно, в отсутствие логарифмически расходящихся членов (например, в случае безмассовых полей в областях с плоской
границей), энергия определяется единственным размерным параметром L — характерным линейным размером. Второй размерный параметр задачи — параметр регуляризации а, который мы выберем также с размерностью длины. Поэтому в отсутствие логарифмической расходимости «минимальное вычитание», состоящее в отбрасывании в выражении для энергии сингулярных по а членов, является не только естественной, но и хорошо обоснованной процедурой. Из очевидных размерных соображений любое слагаемое, пропорциональное (s > 0), будет зависеть от L как т.е. будет пропорционально
неотрицательной степени L. Тем самым мы можем нормировать окончательный ответ в точке L = оо, в которой энергия Казимира должна обращаться в нуль, и вычесть все сингулярные слагаемые. После вычитания в пределе а ^ 0 останется только слагаемое ca°/L (с правильной зависимостью от L), которое и составляет окончательный ответ.
При наличии логарифмической расходимости вычитание становится неоднозначным: чтобы перенормировать слагаемое са° 1п(а/L)/L, из него следует вычесть са° \n(da/L)/L, где d — произвольная константа. В результате окончательный ответ будет содержать неопределенность, связанную с d. Нормировка окончательного ответа при L ^ оо в данном случае невозможна, поскольку неопределенность имеет вид const /L, т. е. неопределенность возникает как раз в той константе, которую нам необходимо найти.
Для свободных полей в шаре логарифмическая расходимость возникает и в безмассовом случае как следствие кривизны границы. При этом неоднозначность в процедуре перенормировки обходится с помощью следующего приема [3, 9]: к «внутренней» задаче (рассматриваются поля внутри полости) добавляется «внешняя» задача (рассматриваются поля вне полости). В указанном случае логарифмические расходимости внешней и внутренней задач компенсируются, и процедура перенормировки становится
однозначной. Однако в таком подходе речь идет уже не о перенормировке вакуумной энергии полей внутри полости, а о перенормировке вакуумной энергии модифицированной задачи, когда поля существуют как внутри, так и вне полости.
Целью настоящей работы является разработка метода, который позволил бы находить энергию Казимира для задач, в которых из-за вычислительных проблем явное аналитическое выражение для энергетических уровней отсутствует, т. е. энергии уровней могут быть определены только численно. Предлагаемая процедура позволяет также обойти неоднозначность, связанную с наличием логарифмической расходимости, и выделить искомую конечную часть именно в исходном (а не в модифицированном) выражении для вакуумной энергии.
Одномерный случай
Простейшей системой, в которой энергия Казимира содержит логарифмическую расходимость, является одномерное массивное скалярное поле на отрезке длины Ь с нулевыми граничными условиями на концах интервала. Исходное выражение для энергии Казимира в этом случае имеет вид
| ОО | оо
= 2 = 2 ^(т/1)2 + т2. (1)
п=1 п=1
Прямолинейный подход [1-3] к перенормировке этой величины состоит в том, что из суммы (1) вычитается энергия, которая была бы у поля, занимающего отрезок Ь, но не ограниченного на этом отрезке. В силу тривиального вида контрчленов такая процедура для одномерного массивного поля оказывается вполне эффективна и приводит к правильному результату. Однако уже для поля в сферической полости эти рассуждения оказываются недостаточными для последовательной перенормировки [3, 14]. Поэтому мы не будем вычислять разность энергий двух различных систем (т.е. величину, физический смысл которой по меньшей мере неочевиден), а рассмотрим перенормировку непосредственно выражения (1).
Регуляризация суммы (1) должна проводиться с помощью параметра а с размерностью длины, входящего в обрезающую функцию Р(ашп):
, оо
= (2) п=1
Асимптотика шп при п » 1 содержит линейно расходящееся слагаемое жп/Ь и логарифмически расходящееся слагаемое т2Ь/(2жп), поэтому разложение
(г)
£саа в ряд по степеням а будет иметь следующую структуру:
£«2 + <'\—т + с\т2 Ып(Ь / а) + с0Ц- +
аг а1 Ь
+ исчезающие при а члены.
Любое вычитание в данном случае будет приводить к неопределенности, поэтому можно попытаться избавиться от процедуры вычитания. Будем вычислять „р rW я /)2cW.
не C-cas , л 0LoCaa ■
Это выражение регулярно при а ^ 0, однозначно
(г)
определяется из £ca,s и поэтому не зависит от способа разбиения последней на сингулярную и регулярную части. Располагая <9|£cas, можно восстановить и саму искомую конечную вакуумную энергию поля
(г)
Ecaa, поскольку из регулярности по а и однозначности следует, что
^ificas = cas> (3)
f г)
а начальные условия для вычисления Е¿as вытекают из очевидных физических соображений: при любой
(г) (г)
разумной регуляризации как Е^ш, так и оьЕсш должны исчезать при L ^ оо:
E&(L^oo) = dLE&(L^oo)=0.
Собственно энергия Казимира поля после этого на-
(г)
ходится через предельный переход: Eca& = lim Е¿а! •
а—г-0
Формализуем теперь изложенную процедуру таким образом, чтобы она не опиралась на явное аналитическое выражение для спектра. Итак, мы располагаем некоторым спектром шп, для которого энергия Казимира содержит логарифмическую расходимость. Прежде всего найдем такой параметр ¿¿, что
оо
£ (4)
п=1
не будет содержать логарифмическую расходимость. Для массивного одномерного поля ¡л будет равно массе поля т. В более нетривиальном случае, например для скалярного поля в сферической полости, это будет некоторый параметр с размерностью массы, характеризующий итоговый (суммарный по всем угловым моментам) коэффициент при логарифмической расходимости.
Далее «добавляем» к полю массу М и исследуем зависимость энергии Казимира от этой добавленной массы в диапазоне от М = 0 до М = оо. Удобно ввести параметр М. такой, что М2 = М2+/л2, и рассматривать зависимость энергии Казимира от М:
, оо п=1
в диапазоне от М = ц до М = оо. Для этого находим (численно, поскольку по условию аналитическое выражение для уровней отсутствует) при все-
возможных значениях М следующую величину:
к м)
д2
¿-cas (
М
М
= д2
м
, оо £
п=1
fi*
м2
1 х
х F
(5)
Здесь, в отличие от (2), в аргумент обрезающей функции .Р(ж) подставлена безразмерная величина,
так что а также безразмерно. Переход к £^(М)/М обусловлен тем, что эта величина становится регулярной уже после двукратного дифференцирования
по М. Для самой £сЙ(М) число дифференцирований пришлось бы увеличить до трех, что крайне нежелательно с вычислительной точки зрения. Для дальнейшего также будет существенно, что после такого перехода в преобразованном спектре шп = \/(ш^ — р?)/М2 + 1 масса М играет роль эффективной «длины», единица — роль эффективной массы, а при М = ¡л восстанавливается исходный спектр шп, деленный на ¡л.
Выражение (5) построено таким образом, что все расходящиеся члены исчезают при дифференцировании. Действительно, при п ^ оо
■Г
UJi,
•г
М
м2
м
■Г
поэтому в сумме, которую порождает первое слагаемое, логарифмическая расходимость отсутствует согласно определению ¡л через условие (4), а остальные расходимости пропорциональны (М/а)2/М и (М/а)1 /М и исчезают при двукратном дифференцировании. Второе слагаемое порождает логарифмическую расходимость, но коэффициент при ней пропорционален М1, так что она тоже исчезает при двукратном дифференцировании.
Таким образом, выражение (5) будет регулярным при а ^ 0. Тогда искомая конечная вакуумная
энергия Есш{М) находится через соотношение, аналогичное (3), именно
д2м (еШ.М)/м) = д2м (£«(М)/М) , (6) с принятой нами точкой нормировки при М = оо
Е<&(М оо) = дмЕ^{М оо) = 0. (7)
Действительно, поле с добавленной нами бесконечной массой М должно иметь нулевую вакуумную энергию. Следовательно, вычислив из (6) вторую
производную &мЕсы}(М) для всех М е [fi, оо], через начальные условия (7) можно затем восстановить
искомое Eca_s = lim ЕсХ{М = fi), соответствующее
а—s-0
исходному спектру.
Реализация этой процедуры для массивного одномерного скалярного поля при т = ц = 1 на отрезке
длиной Ь с нулевыми граничными условиями на концах отрезка дает
L 10-в 10"5
д%{2 Е) -0.5235988 • 1018 -0.5235988 • 1015
2 Е -0.2617989 • 10е -0.2617944 • 10®
L 10-4 М-3
д'1(2 Е) -0.5235988 • 1012 -0.5235986 • 109
2 Е -0.2617494 • 104 -0.2613008 • 103
L 0.01 0.1
д%{2 Е) -0.5235829 • 10е -0.5220101 • 103
2 Е -0.2569007 • 102 -2.182667
L 1.0 10.0
д'1(2 Е) -0.3897729 -7.893493 • Ю-10
2 Е -4.666378 • 10"2 0.0
Приведенные данные получены с помощью регу-ляризирующих функций F(x) = ехр(^ж), F(x) = = ехр(^ж2), F(x) = ехр(^ж3), ..., F(x) = = ехр(^ж6), F(x) = ехр(^2сЬ(ж) + 2). Для всех функций ответы совпадают с избыточной для наших целей точностью до восьми знаков.
Как и следовало ожидать, вторая производная энергии гораздо быстрее выходит на безмассовый предел дрЕ0 = —тг/(6L3) и ^0.52359878/L3 при стремлении L ^ 0, чем сама энергия выходит на свой безмассовый предел 2Eq = —7г/(12L) и «^0.26179939/L.
Трехмерный случай
Изложенный в предыдущем разделе метод безусловно работает для одномерного массивного поля. Чтобы этот способ оказался эффективным и в трехмерном случае, количество дифференцирований следует увеличить на два. Действительно, в трехмерном случае основной сингулярный член будет пропорционален объему и тем самым имеет вид c_4L3/a4, и чтобы от него избавиться, необходимо взять четвертую производную di.
Логарифмический по а член при таком количестве дифференцирований также заведомо исчезнет. Преобразование спектра построено как раз таким образом, чтобы коэффициент при логарифмической расходимости содержал только положительные степени М. Отрицательная степень М 1 умножается на сумму \/w2 — /л2, которая не содержит логарифмической расходимости согласно определению ¡л через условие (4).
Для безмассового скалярного поля в шаре единичного радиуса R = 1 получаем
F(x) ß2 ß 2EaiB/ß
ехр(—ж3) 0.01895799 0.1376880 0.0550299
ехр(—ж5) 0.01895777 0.1376872 0.0550527
ехр(—ж7) 0.01895777 0.1376872 0.0550531
e-2ch(i)+2 0.01895774 0.1376871 0.0550455
Найденные значения четвертой производной для разных .Р(ж) совпадают с точностью до четырех первых знаков, причем наибольшее отклонение дает .Р(ж) = ехр(^ж3). Соответственно точность в определении Есш = 3.790 • 10"3 также можно оценить как четыре верных знака. Точность не слишком высока просто потому, что диапазон использованных значений регуляризирующего параметра а соответствует приблизительно двумстам 5-уровням и соответственно максимальному значению орбитального момента I примерно 600.
Следует отметить, что в ходе вычислений мы нашли не только энергию Казимира в безмассовом случае М = 0, но и энергию Казимира для всех возможных масс М в диапазоне от нуля до «эффективной» бесконечности:
Из приведенных результатов следует, что при таком способе перенормировки вакуумная энергия не имеет инфракрасной особенности при М ^ 0, в отличие от ответа, который приведен в работах [3, 14]. С точки зрения общей теории для полей, локализованных в конечном объеме, появление инфракрасных особенностей при М ^ 0 в энергии вакуума вообще достаточно проблематично. На самом деле появление таких особенностей в [3, 14] может быть связано с тем обстоятельством, что даже при наличии «естественной» точки нормировки М ^ оо в вычитательной процедуре все же остается определенный произвол. Именно, к окончательному выражению для энергии всегда может добавиться функция, которая достаточно быстро убывает при М ^ оо, но при этом имеет логарифмическую особенность при М ^ 0. Использованная нами процедура перенормировки никаких вычитаний и, следовательно, никаких неоднозначностей не содержит: при дифференцировании исключаются только положительные степени М. которые заведомо должны отсутствовать в ответе для энергии Казимира. При этом во всем диапазоне изменения М (от ¡л до оо) ответ для производной конечен и при М ^ оо убывает быстрее, чем М 5 (хотя и медленнее, чем экспонента, на что было впервые указано в работе [14]), так что при М = ц, т.е. при М = 0, мы получаем конечный ответ.
Заключение
Таким образом, в работе предложен метод численного нахождения энергии Казимира, который не содержит процедуры вычитания и, следовательно,
не требует нормировки для устранения неопределенностей в вычитательной процедуре. Подчеркнем, однако, что это один из возможных способов перенормировки.
Недостатком метода является необходимость вычислений с повышенным числом верных знаков (в частности, спектр системы следует находить с высокой точностью). В то же время явным достоинством метода является его идеологическая и техническая тривиальность, позволяющая применять его в тех случаях, когда аналитическое выражение для спектра (и аналитическое выражение для уравнений на уровни спектра) отсутствует. Более того, нет необходимости вычислять приближенные асимптотические выражения для уровней.
Предложенный метод применим, когда отсутствует явное трансцендентное уравнение для уровней. Кроме того, в случае трехмерных задач его целесообразно применять и для относительно простых спектров. Это связано с тем, что в трехмерных задачах (даже в случае сферической симметрии) суммируется как минимум двойной ряд, и даже для тривиальных систем его перенормировка представляет собой довольно громоздкую задачу [14].
Литература
1. Casimir H.B.G. // Proc. Kon. Nederl. Akad. Wet. 1948. 51. P. 793; Physica. 1953. 19. P. 846.
2. Ициксон К., Зюбер Ж.-Б. Квантовая теория поля. Т. 1. М., 1978.
3. Bordag М., Mohideen U., Mostepanenko V.M. // Phys. Reports. 2001. 353. P. 1.
4. Boyer Т.Н. // Phys. Rev. 1968. 174. P. 1764.
5. Balian R., Duplantier R. // Ann. Phys. 1978. 112. P. 165.
6. Milton K.A., De Raad L.L., Schwinger J. // Ann. Phys. 1978. 115. P. 388.
7. Milton К.A. // Phys. Rev. 1980. 22D. P. 1441, 1444.
8. Brown O.E., Jackson A.D., Rho M., Vento V. // Phys. Lett. 1984. B140. P. 285.
9. Vepstas L., Jackson A.D. // Phys. Reports. 1990. 187. P. 109.
10. De Francia M., Falomir H., Santangelo E.M. // Phys. Rev. 1992. 45D. P. 2129.
11. Hosaka A., Toki H. // Phys. Rep. 1996. 277. P. 65.
12. Baacke J., Igarashi Y. // Phys. Rev. 1983. 27D. P. 460.
13. Blau S.K., Wisser M., Wipf А. // Nucí. Phys. 1988. B310. P. 163.
14. Bordag M., Elizalde E., Kirsten К., Leseduate S. // Phys. Rev. 1997. 56D. P. 4896.
15. Elizalde E., Bordag M., Kirsten К. // J. Phys. 1998. A31 P. 1743.
16. Bordag M., Geyer В., Kirsten К., Elizalde E. // Comm. Math. Phys. 1996. 179. P. 215.
Поступила в редакцию 06.10.04
м 0 0.1 0.2 0.5
2 Е(М) 7.580 • Ю-3 7.372 • Ю-3 6.761 -10"3 4.236 • Ю-3
М 1.0 2.0 5.0 10.0
2 Е(М) 1.649 • Ю-3 4.272 • Ю-4 9.843 • Ю-5 1.368 • 10"5