Аннотация:
Излагаются основные физические и математические идеи (П. Кюри, Дарбу, Пуанкаре, Дирак), приведшие к понятию магнитного заряда, общая конструкция магнитных лапласианов для магнитных монополей на римановых многообразиях и результаты Ю. А. Кордюкова и автора по квазиклассическому приближению для собственных сечений этих операторов.
Данная статья посвящена схеме квазиклассического квантования магнитных монополей, предложенной Кордюковым и автором в работе [1], примыкающей к статьям этих же авторов по формулам следа для магнитных лапласианов [2], [3]. Она является расширенным изложением доклада автора на международной конференции по математической физике, посвященной столетию со дня рождения В. С. Владимирова, в январе 2023 г. в Математическом институте им. В. А. Стеклова РАН.
Мы подробно останавливаемся на предыстории вопроса.
В разделе 2 обсуждаются работы П. Кюри, Дарбу и Пуанкаре, в которых с различных точек зрения рассматривались магнитные монополи в классической электродинамике. В разделе 3 излагается вывод условия квантования электрических зарядов, ради которого Дирак и ввел магнитные монополи в квантовой механике. Эта идея оказалась плодотворной, современное состояние ее развития изложено в монографии [4]. Общая конструкция магнитных лапласианов представлена в разделе 4.
В разделах 5, 6 мы излагаем на физическом уровне конструкцию квазиклассического квантования магнитного лапласиана. Она основана на расширении многомерного метода Вентцеля–Крамерса–Бриллюэна (ВКБ) (канонического оператора Маслова) на случай операторов, действующих на сечениях нетривиальных расслоений. Такое изложение делает конструкцию прозрачной, а за строгими математическими деталями мы отсылаем к работе [1].
2. Кюри, Пуанкаре и магнитный монополь
По-видимому, впервые с физической точки зрения магнитные монополи обсуждались Кюри. В краткой заметке [5] он привел некоторые физические следствия их существования. В 1894 г. Кюри пытался использовать эти выводы в своих экспериментах по обнаружению монополей, которые не привели к успеху. Тем не менее он указал, что1[x]1 “… au point de vue de l’énergétique, au point de vue de la symétrie, on peut concevoir sans absurdité les courants de magnétisme et les charges de magnétisme libre. Il serait certes téméraire d’induire de là que ces phénom‘enes existent réellement. Si cependant il en était ainsi, ils devraient satisfaire aux conditions que nous avons énoncées”.“с точки зрения энергетики, с точки зрения симметрии, можно представить себе без абсурда токи магнетизма и свободные магнитные заряды. Конечно, было бы опрометчиво делать из этого выводы, что эти явления действительно существуют. Но если бы это было так, они должны удовлетворять указанным нами условиям”.
В 1896 г. Пуанкаре предпринял попытку найти объяснение результатам экспериментов Биркеланда, норвежского физика и его бывшего ученика, с трубками Крукса (электрическими разрядными трубками). Обнаруженные в них катодные лучи были тогда предметом интенсивного исследования. Пуанкаре предложил математическую модель обнаруженного Биркеландом поведения катодных лучей [6]. Он начал со следующего предположения2[x]2 “… nous ecrirons les equations du rayon cathodique, en l’assimilant a une particule materielle en mouvement rapide, chargee d’electricite …”.: “напишем уравнения катодного луча, уподобляя его быстро движущейся материальной электрически заряженной частице”, и далее3[x]3 “Supposons un seul pole magnetique, que nous prendrons pour l’origine …”.: “Предположим, что есть одиночный магнитный полюс, который мы примем за начало координат”.
После этого Пуанкаре выписал уравнения движения r(t)=(x(t), y(t), z(t)) заряженной частицы в магнитном поле с особенностью в начале координат:
где r=|r|, а λ – постоянная, зависящая от интенсивности магнитного поля и “природы катодных лучей”4[x]4 “c’est-a-dire, dans l’hypothese de Crookes, de la masse de la particule materielle en mouvement et de sa charge electrique” – пер.: “то есть согласно гипотезе Крукса, от массы движущейся материальной частицы и ее электрического заряда”..
Обозначая через ˙r и ¨r скорость и ускорение частицы, перепишем эти уравнения в виде
¨r=λr3[r×˙r].
Из уравнений (1) Пуанкаре выводит закон сохранения энергии
|˙r|2=C=const
и важное соотношение
r2=Ct2+2Bt+A,
где постоянные A, B и C зависят от начальных условий. Действительно,
d2r2dt2=d2⟨r,r⟩dt2=2⟨˙r,˙r⟩+2⟨r,¨r⟩=2|˙r|2=2C,
где через ⟨⋅,⋅⟩ мы обозначаем скалярное (евклидово или эрмитово) произведение векторов. Прямым следствием основных уравнений также является то, что вектор [r×˙r]+(λ/r)r тоже постоянный, и мы получаем еще три закона сохранения:
[r×˙r]+λrr=v=(abc),
где вектор v является постоянным вдоль траекторий. Отсюда немедленно следует, что
⟨v,r⟩=ax+by+cz=λr,
т. е. частица движется по конусу вращения с вершиной в начале координат. Так как вектор ускорения ¨r в каждой точке конуса, кроме его вершины, перпендикулярен вектору скорости ˙r и образующей r этого конуса, то траектория является геодезической.
Таким образом, траектории движения частицы массы m и с зарядом q в магнитном поле (монополя)
H=H0r3r,H0=−mcλq,
где c – скорость света в вакууме, являются геодезическими на конусах (3), где векторы v определяются начальными данными согласно (2).
Заметим, что это объяснение Пуанкаре зависит от предположения Крукcа о том, что катодные линии являются потоками массивных заряженных частиц5[x]5 “si l’hypothese de Crookes n’est pas vraie, il semble bien que tout se passe comme se elle l’etait” – пер.: “если гипотеза Крукса неверна, то кажется, что все остается по-старому”.. Более того, он подчеркивает, что6[x]6 “Est-ce que la theorie est incomplete, parce que nous avons suppose un pole magnetique unique”. “теория неполна, поскольку мы предполагаем существование одиночного магнитного заряда”.
Уравнения (1), как оказалось, уже рассматривались в 1878 г. Дарбу в краткой заметке [7], не предполагавшей существование магнитного монополя, но тоже связанной с теорией магнетизма. Дарбу решил задачу о нахождении равновесия гибкого и нерастяжимого невесомого провода, проводящего ток и находящегося под действием полюса магнита7[x]7Аннотация к [7]: “Trouver la figure d’équilibre d’un fil flexible et inextensible non pesant, traversé par un courant et sounis à l’influence du póle d’un aimant”.. Дарбу свел ее к уравнениям (1) и описал возникающие кривые как геодезические на конусах вращения, совершенно аналогично тому, как это позднее сделал Пуанкаре (он видимо не был знаком с работой Дарбу и не ссылался на нее).
Через год после работы Пуанкаре, изучая катодные лучи, Томсон показал, что они являются потоками заряженных частиц, которые были названы электронами [8].
3. Монополь Дирака
В работе Дирака [9], занимающей чуть более 12 страниц, первая страница посвящена развитию математических основ современной физики и приведенные на ней рассуждения до сих пор не потеряли актуальность, вторая – обсуждению существования странных частиц, которые имеют ту же массу, что и электрон, но отрицательную энергию. Было отмечено, что их можно трактовать как дырки в ненаблюдаемом полностью заполненном распределении состояний с отрицательной энергией. Как отмечает Дирак8[x]8 “A hole, if there were one, would be a new kind of particle, unknown to experimental physics, having the same mass and opposite charge to an electron. We may call such a particle an anti-electron. We should not expect to find any of them in nature, on account of their rapid rate of recombination with electrons, but if they could be produced experimentally in high vacuum they would ne quite stable and amenable to observation”., “Дырка, если бы нашлась хоть одна, была бы частицей нового сорта, не известной экспериментальной физике, имеющей такую же массу, что и электрон, и противоположный заряд. Можно назвать такую частицу антиэлектроном. Не следует надеяться обнаружить какие-нибудь из них в природе из-за большой скорости их рекомбинации с электронами, и все же если бы они были получены экспериментально, то в высоком вакууме были бы вполне стабильными и поддающимися наблюдению ”.
Существование таких частиц вытекало из знаменитого уравнения Дирака, выведенного в 1928 г. Через год после выхода статьи [9], в 1932 г. эти частицы были открыты Андерсоном и известны сейчас как позитроны.
Далее Дирак переходит к основной цели статьи, состоящей в том, чтобы9[x]9 “to put forward a new idea which is in many respects comparable with this one about negative energies. It will be concerned essentially, not with electrons and protons, but with the reason for the existence of a smallest electric charge”.“выдвинуть новую идею, которая во многих отношениях сравнима с идеей об отрицательных энергиях. Она будет, по существу, относиться не к электронам и протонам, а к причине существования наименьшего электрического заряда”.
Пусть волновая функция ψ(x,y,z,t) задает движение частицы. Представим ее в виде
ψ=Aeiγ,
где A и γ – вещественнозначные функции. Если предположить, что функция ψ нормированная, она определяет состояние с точностью до произвольного постоянного множителя, равного по модулю нулю. Значит, можно предположить, что γ не имеет определенного значения в точке и для двух точек разность фаз определена по отношению к кривой, их соединяющей.
Для пары волновых функций φ и ψ модуль величины
⟨φ|ψ⟩=∫ˉφψdxdydz
имеет физический смысл и, как отмечает Дирак, мы должны предположить, что10[x]10 “The change in phase of a wave function round any closed curve must be the same for all the wave functions”.“изменение фазы волновой функции при обходе любой замкнутой кривой должно быть одинаковым для всех волновых функций”.
Достичь выполнения этих требований и принципа суперпозиции можно следующим образом. Пусть
ψ=ψ1eiβ,
где ψ1 – волновая функция с фазой, определенной в каждой точке, а неопределенность фазы ψ выражается в множителе eiβ. Предположим, что фаза не определена в каждой точке, но ее производные
ϰx=∂β∂x,ϰy=∂β∂y,ϰz=∂β∂z,ϰt=∂β∂t
определены, но при этом не обязаны удовлетворять условиям интегрируемости
∂ϰx∂y=∂ϰy∂x,….
Если ψ удовлетворяет какому-то уравнению, содержащему операторы импульса и энергии p и W, то ψ1 удовлетворяет тому же уравнению, в котором p и W заменены на p+hϰ и W−hϰ0. Таким образом, переход от ψ к ψ1 сводится к включению электромагнитного поля с потенциалом
A=ℏceϰ,
где h – постоянная Планка, c – скорость света в вакууме, e – элементарный заряд (положительный заряд, равный по модулю заряду электрона). При этом Дирак отмечает, что он рассматривает частицу с зарядом q=−e. Тем самым общая формула для A принимает вид
A=−ℏcqϰ,
где q – заряд частицы, для которой ψ является волновой функцией.
Дирак замечает, что связь между неинтегрируемостью фазы и электромагнитным полем является проявлением вейлевского принципа калибровочной инвариантности. Однако есть два дополнительных обстоятельства:
1) фаза всегда определена с точностью до числа, кратного 2π;
2) изменения фаз различных волновых функций вдоль замкнутых кривых могут быть различными и различаться на величины, кратные 2π.
Если первое утверждение очевидно, то остановимся на проявлениях второго утверждения. Изменение фазы вдоль замкнутой кривой τ в (x,y,z)-пространстве согласно теореме Стокса равно
∫τ(ϰxdx+ϰydy+ϰzdz)=∫Γ(curlϰ,dS),
где dS – элемент площади поверхности Γ, ограничивающей τ. Если в некоторой точке P мы имеем ψ(P)≠0, то при стягивании замкнутых кривых, обходящих эту точку, изменение фазы стремится к нулю. Если же ψ(P)=0, то мы можем только утверждать, что
eℏc∫Γ(H,dS)=2πnτ,
где H=(ℏc/e)curlϰ – магнитное поле, τ=∂Γ.
Если мы имеем замкнутую поверхность Γ такую, что она разбивается на малые области, содержащие по одному нулю функции ψ, мы, складывая предыдущее соотношение по всем таким областям, получим
eℏc∫Γ(H,dS)=2πn,n∈Z,
т. е. поток магнитного поля такой системы квантуется и кратен 2π. При N≠0 область, ограниченная поверхностью Γ, должна содержать магнитные заряды.
Рассмотрим простейший случай одноточечного магнитного заряда, находящегося в начале координат, и предположим, что электрическое поле равно нулю. Тогда магнитное поле имеет вид
H=gDrr3,
где r – радиус-вектор точки, r=|r|. Рассмотрим сферические координаты r, φ, θ. В качестве вектор-потенциала этого магнитного поля можно взять
ϰθ=ϰr=0,ϰφ=12rtgθ2.
В этом случае с помощью разделения переменных представим ψ1 в виде
и назовем оператор ΔL магнитным лапласианом на двумерной единичной сфере. Общее определение этого оператора и смысл символа L мы изложим в разделе 4.
В данном случае общее решение волнового уравнения было получено Таммом11[x]11Дирак ссылается на эту работу, вышедшую через месяц после [9]. [10]. Он показал, что собственные значения EN, N=0,1,…, имеют вид
EN=N2+2N+12
с кратностью 2N+2. Базис собственных “функций” при наименьшем собственном значении E=1/2 задан “функциями”
Sa=cosθ2,Sb=sinθ2e−iφ.
Причем Sa непрерывна всюду, а Sb имеет особенность при θ=π и ее фаза изменяется при обходе этой точки по малому контуру на 2π.
Базис (8) был получен Таммом [10]. В работе Дирака [9] он был воспроизведен с опечаткой: был неверный знак при φ в определении Sb. Приведенные ниже формулы (15) показывают, что Sa и Sb как сечения линейного расслоения с c1=−1 не имеют особенностей.
Однако основной целью Дирака был важный физический вывод. Из условия квантования (5) следует, что
egD=ℏc2n,n∈Z.
В § 1 в [9] e обозначает “наименьший” электрический заряд, но в § 3 в [9] в качестве e взят заряд произвольной частицы. Поэтому (9) есть условие квантования зарядов электрических частиц. Как отметил Дирак12[x]12 “… if there exists any monopole at all in the universe, all electric charges would have to be such that e times this monopole strength is equal to 12nℏc”. [11], “… если существует всего лишь один монополь во всей Вселенной, то все электрические заряды вынуждены быть такими, чтобы величина e, умноженная на заряд этого монополя, была равна 12nℏc ”.
Монополь Дирака до сих пор не обнаружен. Оценим его физические характеристики.
Поскольку в системе СГС постоянная тонкой структуры имеет вид
α=e2ℏc≈1137,
мы перепишем при n=1 условие квантования (5) в виде
gD≈1372e.
Мы не можем оценить массу монополя (как и электрона) из каких-либо теоретических предположений. Но если, например, предположить, что классические радиусы электрона и магнитного монополя совпадают:
re=rD,
где
re=e24πε0mec2,rD=g2D4πμ0mDc2
и ε0 и μ0 – электрическая и магнитная постоянные, то
mD=g2De2me=1(2α)2me≈4692me≈2.4 ГэВ.
4. Магнитный лапласиан
Мы привели выше рассуждения Дирака, чтобы показать, как он подошел с физической точки зрения к такому важному понятию, как U(1)-расслоение, и к условию квантования классов Черна таких расслоений. Объяснение работы Дирака с помощью этих топологических понятий было дано в работах [12], [13].
Мы рассмотрим более общую ситуацию, когда система определена на d-мерном дифференцируемом многообразии M. Все объекты на этом многообразии – формы, связности, метрики и т. д. – мы предполагаем дифференцируемыми необходимое число раз.
Под магнитным полем понимается замкнутая 2-форма F(0):
F(0)=∑F(0)jkdxj∧dxk,dF(0)=0.
Вектор-потенциал этой формы определяется как
A(0)=(A(0)k),∂A(0)k∂xj−∂A(0)j∂xk=F(0)jk.
Мы предполагаем, что c=ℏ=1, как принято в естественной системе единиц.
Если форма F(0) некогомологична нулю, то вектор-потенциал определяется только локально, в областях, на которые F(0) ограничивается до точной формы.
Пусть на M задана риманова метрика gjkdxjdxk (здесь и далее мы подразумеваем суммирование по верхним и нижним повторяющимся индексам), с помощью которой определяется гамильтониан (кинетическая энергия) движения частиц на M:
H=m2gjkpjpk,
где gjk – тензор, обратный к gjk:gjkgkl=δjl, p=(p1,…,pd)⊤ – импульс частицы. Далее для простоты положим m=1. Перепишем гамильтониан в виде
2H=1√gpj√ggjkpk,g=det(gjk),
и заменим импульс на отвечающий ему при квантовании оператор импульса,
pj→−i∂∂xj,
получив оператор Лапласа–Бельтрами как результат квантования гамильтониана свободной частицы:
Δ=−1√g∂∂xj√ggjk∂∂xk.
Этот оператор задается ковариантным выражением, т. е. инвариантен относительно замен координат.
Предположим теперь, что частица заряжена и имеет заряд q=Ze. В этом случае включение взаимодействующего с частицей магнитного поля F(0) состоит в следующей замене импульса:
pj→pj−qA(0)j.
Для краткости введем обозначения
F=qF(0),A=qA(0),dA=F.
При квантовании это приводит к замене оператора Лапласа–Бельтрами на магнитный лапласиан:
Δ=−1√g(∂∂xj−iAj)√ggjk(∂∂xk−iAk).
Рассуждения Дирака о неоднозначности функций ψ переформулируются следующим образом:
который является первым классом Черна расслоения L.
Заметим, что если условие квантования не выполняется, то A не задает связности в каком-либо линейном расслоении.
Магнитный лапласиан
Δ=−1√g∇j√ggjk∇k
инвариантен относительно калибровочных преобразований:
˜Δ(η⋅ψ)=η⋅Δψ,
но его определение существенно зависит от выбора связности A с заданной формой кривизны (“магнитным полем”) F.
Предположим, что связности A и A′ задают операторы Δ и Δ′, отвечающие одному и тому же магнитному полю F и действующие на сечениях одного и того же расслоения L. Тогда их разность Q=A′−A является замкнутой 1-формой: dQ=q(F−F)=0.
Если замкнутая форма Q=A′−A реализует нулевой класс одномерных когомологий:
[A′−A]=0∈H1(M;R),
то A′−A=df, где f:M→R – некоторая функция. Тогда согласно (13) из
Δψ=Eψ
следует, что
Δ′(e−ifψ)=E(e−ifψ),
и тем самым спектр магнитного лапласиана не зависит от выбора связности qA.
Если расслоение L тривиально, то, выбрав его тривиализацию L=M×C, можно определить связность на нем с помощью 1-формы A=Akdxk на многообразии M.
Пример “Монополь Дирака” на двумерной сфере. Рассмотрим единичную двумерную сферу M=S2 в трехмерном пространстве и ограничим на нее 2-форму, отвечающую магнитному полю (6). Мы получим
F(0)=gDsinθdθ∧dφ.
Условие квантования дает
2egD∈Z.
В области, полученной из сферы выкалыванием нижнего полюса (θ=π), вектор-потенциал магнитного поля F(0) можно задать в виде
A(0)θ=0,A(0)φ=gD(1−cosθ).
Умножим A(0) на заряд частицы
q=Ze,Z∈Z.
При
egD=12,q=Ze,Z=−1(электрон),
подставляя выражение для A=qA(0) в (12), получим в точности оператор Δ из уравнения (7).
Согласно монографии [10] наименьшее собственное значение для этого оператора равно λ=1/2. Базис собственных сечений можно выбрать в виде
Sa=cosθ2,Sb=sinθ2e−iφ.
Мы уже упоминали замечание Дирака о том, что Sb имеет особенность при θ=0. Как показали впервые Ву и Янг [13], эти “функции” являются сечениями линейного расслоения L и не имеют особенностей.
U(1)-расслоения над двумерной сферой классифицируются первым классом Черна. Рассмотрим стереографическую проекцию единичной сферы x2+y2+z2=1 на плоскость (x,y), которую мы отождествим с комплексной прямой с координатой z=x+iy. Эта проекция сохраняет ориентации и отождествляет двумерную сферу с комплексной проективной прямой CP1=C∪{∞}. Комплексное линейное расслоение L с c1(L)=n над CP1 обозначается через O(n).
Любое U(1)-расслоение, т. е. комплексное одномерное расслоение со структурной группой U(1), над CP1 задается функцией сцепления.
Пусть [u:v] – однородные координаты на CP1: точка на CP1 задается парой координат [u:v], которые не равны одновременно нулю и определены с точностью до умножения их на ненулевую постоянную, т. е. [u:v]∼[λu:λv], λ≠0. Сфера представляется как объединение перекрывающихся областей
Ua={u≠0},Ub={v≠0}
или как склейка двух дисков
Da={[1:z]∈Ua,|z|⩽
по общей границе. Так как диски D_a и D_b стягиваемы, то ограничения любого линейного расслоения на них тривиальны. Два тривиальных расслоения D_a\times\mathbb C и D_b\times\mathbb C склеиваются по границам в линейное расслоение над \mathbb CP^1. Склейка задается функцией сцепления \gamma(z), z=e^{i\varphi}, |\gamma(z)|=1, по формуле
и топологический тип расслоения определяется гомотопическим классом отображения \gamma\colon S^1\to S^1. Любое отображение из S^1 в S^1 гомотопно отображению вида \varphi\to e^{ik\varphi}, где k\in\mathbb Z.
Приведем в качестве примера тавтологическое расслоение \mathcal O(-1). Так как точке [u:v]\in\mathbb C P^1 отвечает прямая с направляющим вектором (u,v) в \mathbb C^2, то все такие прямые образуют линейное расслоение над \mathbb C P^1. На D_a расслоение тривиализуется D_a\times\mathbb C, и при этом паре координат (z,\lambda) отвечает точка с координатами (\lambda,\lambda z). Аналогично при тривиализации D_b\times\mathbb C имеет место соответствие (w,\mu)\to(\mu w,\mu). Эти расслоения склеиваются по границам D_a и D_b по очевидным правилам
Для двойственного расслоения \mathcal O(1)=\mathcal O(-1)^\ast функция сцепления суть e^{-i\varphi}. Для тензорных степеней этих расслоений \mathcal O(-1)^k и \mathcal O(1)^k, k>0, функции сцепления суть e^{-ik\varphi} и e^{ik\varphi}. Очевидно, что для тривиального расслоения \mathbb CP^1\times\mathbb C в качестве функции сцепления можно положить \gamma=1. Так как
и c_1(L\otimes L')=c_1(L)+c_1(L') для любой пары линейных расслоений L и L', то c_1(\mathcal O(k))=k, и для любого линейного расслоения L над \mathbb CP^1 функцию сцепления можно взять в виде
и не имеют особенностей. Они является сечениями тавтологического расслоения \mathcal O(-1) над \mathbb CP^1.
Явные формулы для монопольных гармоник – собственных сечений магнитных лапласианов, отвечающих постоянным магнитным полям на двумерной сфере, – выведены в работе [13].
Очевидно, что если \Delta^L\psi=E\psi, то, комплексно сопрягая аналитическое выражение для \psi, мы получим сечение \bar\psi сопряженного расслоения L^\ast=\overline L и \Delta^{L^\ast}\bar\psi=E\bar\psi. Физически это трактуется как обращение заряда частицы.
В квантовой механике метод ВКБ позволяет строить (квазиклассические) приближения собственных функций волновых операторов. Так как при этом “постоянная Планка” \hbar является малым параметром, по степеням которого строятся приближения, в дальнейшем мы не будем переходить в естественную систему измерений и полагать ее равной единице.
Напомним метод ВКБ для одномерного уравнения Шредингера
Если мы имеем движение точки в потенциальной яме \{U(x)\leqslant E\}, ограниченной парой точек, в которых U(x)=E, то вне этой ямы импульс p становится чисто мнимым. Хотелось бы с помощью выбора постоянных C_1 и C_2 достичь такой ситуации, чтобы вне ямы эта формула давала экспоненциально затухающие решения. Это можно сделать при выполнении условия Бора–Зоммерфельда
где интеграл взят по замкнутой траектории частицы в потенциальной яме. Естественно, это возможно не при всех значениях энергии E.
Метод канонического оператора Маслова, предложенный в [14] (см. также [15]), является многомерным вариантом метода ВКБ. В случае многомерного оператора Шредингера
он состоит в следующем. Пусть M – риманово многообразие размерности d с метрикой g_{jk}, а \Delta – оператор Лапласа–Бельтрами (11). Рассмотрим на кокасательном расслоении T^\ast M со стандартной симплектической структурой
Пусть \Lambda – лагранжево подмногообразие в T^\ast M, инвариантное относительно гамильтонова потока и с инвариантной мерой d\mu на \Lambda. Напомним, что подмногообразие T^\ast M лагранжево, если оно d-мерно и ограничение симплектической формы на него равно нулю.
Выберем покрытие \{V_\alpha\} подмногообразия \Lambda односвязными открытыми множествами (картами) с локальными координатами вида x^{j_1},\dots,x^{j_m}, p_{k_1},\dots,p_{k_{d-m}}. Выберем на \Lambda разбиение единицы \{\phi_\alpha\}: \phi_\alpha\colon\Lambda\to\mathbb R, 0\leqslant\phi_\alpha\leqslant 1, \sum\phi_\alpha\equiv 1, \phi_\alpha=0 вне V_\alpha.
Канонический оператор действует на дифференцируемых функциях u(s) на \Lambda и переводит их в функции на M. Он определяется для каждой карты V_\alpha с каноническими координатами (x^I,p_I): K^\hbar_\Lambda(V_\alpha)(\phi_\alpha u)(x), и общий оператор имеет вид
где \pi\colon\Lambda\to M – проекция \Lambda на M, \pi(y)=x, а \sqrt{g}\,dx – форма объема на M (g=\det(g_{jk}), dx=dx^1\wedge\dotsb\wedge dx^d). Здесь
– действие, которое получается интегрированием вдоль пути в V_\alpha\subset\Lambda из некоторой начальной точки y_0 в y. Так как область V_\alpha односвязна, а d(\sum p_k\,dx^k)=\Omega_0 и \Omega_0|_\Lambda=0, значение интеграла не зависит от выбора пути.
Сравним (16) и (19). Пусть \gamma – периодическая траектория частицы в одномерном потенциальном поле. Она является лагранжевым подмногообразием в T^\ast\mathbb R. Инвариантная мера на ней имеет вид
Замкнутая траектория \gamma с выкинутыми точками перегиба, в которых p=0, разбивается на два интервала, которые проецируются на интервал. Они задают две карты на \gamma, проходятся по отношению к x в разных направлениях, и отсюда возникают знаки \pm при \int p\,dx. Условие (17) заменяется в многомерной ситуации следующим условием квантования: для любой замкнутой кривой \gamma на \Lambda имеет место
где \mu(\gamma) – индекс Маслова кривой \gamma. В одномерной ситуации (21) сводится к (17).
Квазиклассические собственные функции оператора \widehat H строятся по лагранжевым подмногообразиям, лежащим на уровне энергии E и удовлетворяющим условию квантования (21). Для их построения канонический оператор применяется к функции u\equiv 1 на \Lambda:
Движение частицы в магнитном поле описывается гамильтоновой системой на T^\ast M с гамильтонианом H(x,p)=|p|^2/2 по отношению к скрученной симплектической структуре [16]:
Если форма F точна, то к системе можно применить канонический оператор, как это сделано в разделе 5. Только надо рассматривать лагранжевы подмногообразия по отношению к скрученной симплектической структуре и определять действие как
В разделе 4 при изложении общей схемы квантования мы предполагали \hbar=1. Здесь же у нас \hbar – малый параметр, и мы должны рассматривать выражения вида
так как при \hbar=1/N выражение (23) задает сечение расслоения L^N, т. е. N-й тензорной степени расслоения L.
Схема квазиклассического квантования магнитного монополя состоит в следующем [1]:
1) выбор лагранжева (по отношению к скрученной форме \Omega) подмногообразия \Lambda\subset T^\ast M, на котором гамильтониан равен постоянной: H\equiv E/2 (в отличие от рассмотренных ранее операторов Шредингера, здесь мы опускаем множитель 1/2 перед \Delta);
2) выбор покрытия \Lambda такими картами, что над их проекциями на M расслоение L тривиально;
3) для каждой карты V выбор ограничения на нее связности A_V=A_{V,k}\,dx^k;
4) далее построение канонического оператора по обычной схеме с заменой везде p_k на p_k-A_k, в частности для карты, проецирующейся на область из M, оператор K^\hbar_\Lambda(V) строится по формуле (19) с заменой действия S на (22), для других карт это потребует некоторых модификаций (см. [1]).
Для \hbar=1/N также должно выполняться дополнительное условие квантования на \Lambda (оно зависит от N=\hbar^{-1}): для любой замкнутой кривой \gamma на \Lambda
Заметим, что собственные сечения и их квазиклассические приближения являются сечениями одного и того же расслоения.
Поиск серий почти собственных значений сводится теперь к нахождению лагранжевых многообразий, удовлетворяющих условию квантования (25). В работе [1] это было сделано для монополя Дирака.
Гамильтонова система, отвечающая монополю Дирака на двумерной сфере (см. разделы 3, 4), описывает движение частицы с зарядом q=Ze на S^2 во внешнем магнитном поле F^{(0)}=\frac{1}{2}\sin\theta\,d\theta\wedge d\varphi. Эта система интегрируема, и два ее первых интеграла имеют вид
т. е. приближенный спектр сдвигается на постоянную 1/4 с сохранением кратности.
7. Заключительные замечания
Как мы уже отмечали, почти собственные числа находятся по инвариантным торам, удовлетворяющим условиям квантования (25). Выписать отвечающие им почти собственные сечения тоже несложно. Недавно явные формулы для случая монополя Дирака были получены Кордюковым и автором [17].
В [18] предложен иной подход к построению почти собственных значений – асимптотическое квантование. Оно сопоставляет символам на произвольном симплектическом многообразии операторы, которые действуют на так называемых пучках волновых пакетов на данном многообразии. Обсуждение асимптотического квантования магнитного лапласиана дано в работе [19], где подчеркнуто, что пучки волновых пакетов являются более сложными объектами, чем U(1)-расслоения, однако локально в координатах действия асимптотически квантованного лапласиана и самого лапласиана почти совпадают. В [19] найдены серии асимптотических собственных чисел для магнитных монополей на замкнутых гиперболических поверхностях в постоянных магнитных полях. Они строятся по лагранжевым торам, удовлетворяющим условиям квантования. Формулы для волновых пакетов (“асимптотических собственных функций”) пока не выписывались.
На отличие асимптотического квантования от квазиклассического указывает отсутствие в нем аналогов условия квантования постоянной Планка (24) и то, что настоящие собственные функции и их асимптотические приближения являются сечениями разных расслоений.
Благодарности
Автор благодарит Ю. А. Кордюкова за полезные обсуждения.
Конфликт интересов
Автор заявляет, что у него нет конфликта интересов.
Список литературы
1.
Ю. А. Кордюков, И. А. Тайманов, “Квазиклассическое приближение для магнитных монополей”, УМН, 75:6(456) (2020), 85–106
2.
Ю. А. Кордюков, И. А. Тайманов, “Формула следа для магнитного лапласиана”, УМН, 74:2(446) (2019), 149–186
3.
Y. A. Kordyukov, I. A. Taimanov, “Trace formula for the magnetic Laplacian on a compact hyperbolic surface”, Regul. Chaotic Dyn., 27:4 (2022), 460–476
4.
Ya. M. Shnir, Magnetic Monopoles, Theoretical and Mathematical Physics, Springer, Berlin, 2005
5.
P. Curie, “Sur la possibilité d'existence de la conductibilité magnétique et du magnétisme libre”, J. Phys. Theor. Appl., 3:1 (1894), 415–417
6.
H. Poincaré, “Remarques sur une expérience de M. Birkeland”, C. R. Acad. Sci., 123 (1896), 530–533
7.
G. Darboux, “Problème de mécanique”, Bulletin des Sciences Mathématiques et Astronomiques, Serie 2, 2:1 (1878), 433–436
8.
J. J. Thomson, “Cathode rays”, Philos. Mag., 44:269 (1897), 293–316
9.
П. А. М. Дирак, “Квантованные сингулярности в электромагнитном поле”, Собрание научных трудов, т. II, Квантовая теория (научные статьи 1924–1947 гг.), Физматлит, М., 2003, 388–398
10.
И. Е. Тамм, “Обобщенные шаровые функции и волновые функции электрона в поле магнитного полюса”, Собрание научных трудов, Наука, М., 1975, 186–195
11.
П. А. М. Дирак, “Концепция монополя”, Собрание научных трудов, т. III, Квантовая теория (научные статьи 1948–1984 гг.), Физматлит, М., 2004, 189–200
12.
T. T. Wu, C. N. Yang, “Concept of nonintegrable phase factors and global formulation of gauge fields”, Phys. Rev. D, 12:12 (1975), 3845–3857
13.
T. T. Wu, C. N. Yang, “Dirac monopole without strings: monopole harmonics”, Nucl. Phys. B, 107:3 (1976), 365–380
14.
В. П. Маслов, Теория возмущений и асимптотические методы, Изд-во Моск. ун-та, M., 1965
15.
В. П. Маслов, М. В. Федорюк, Квазиклассическое приближение для уравнений квантовой механики, Наука, М., 1976
16.
С. П. Новиков, “Гамильтонов формализм и многозначный аналог теории Морса”, УМН, 37:5(227) (1982), 3–49
17.
Ю. А. Кордюков, И. А. Тайманов, “Квазиклассическое приближение монопольных гармоник”, Матем. заметки, 114:6 (2023), 848–862
18.
М. В. Карасев, В. П. Маслов, Нелинейные скобки Пуассона. Геометрия и квантование, Наука, М., 1991
19.
Й. Брюнинг, Р. В. Некрасов, А. И. Шафаревич, “Квантование периодических движений на компактных поверхностях постоянной отрицательной кривизны в магнитном поле”, Матем. заметки, 81:1 (2007), 32–42
Образец цитирования:
И. А. Тайманов, “Геометрия и квазиклассическое квантование магнитных монополей”, ТМФ, 218:1 (2024), 149–167; Theoret. and Math. Phys., 218:1 (2024), 129–144
И. А. Тайманов, “Функции Флоке–Блоха на неодносвязных многообразиях, потоки Ааронова–Бома и конформные инварианты погруженных поверхностей”, Геометрия, топология, математическая физика, Сборник статей. К 85-летию академика Сергея Петровича Новикова, Труды МИАН, 325, МИАН, М., 2024, 297–308; I. A. Taimanov, “Floquet–Bloch Functions on Non-simply Connected Manifolds, the Aharonov–Bohm Fluxes, and Conformal Invariants of Immersed Surfaces”, Proc. Steklov Inst. Math., 325 (2024), 280–291