Квантова теорія поля на гратці

Теорія поля на ґратці це - галузь теоретичної фізики, зокрема квантової теорії поля, що вивчає теорії поля визначені на просторових або часо-просторових ґратках.

Теоретичний опис основних складових матерії та взаємодія між ними базується на квантовій теорії поля. Основними елементами теорії поля є поля. Вони є функціями $\phi$ кожної точки $x$ часу-простору $\phi (x)\ .$

У випадку класичних теорій поля, $\phi (x)$ зазвичай є елементом дійсного чи комплексного многовиду скінченної розмірності, котрий в більшості випадків є лінійним простором. Наприклад:

• дійсне скалярне поле $\phi (x)\in \mathbf {R} \ ,$
• комплексне скалярне поле $\phi (x)\in \mathbf {C} \ ,$
• $n$-векторне поле $\phi (x)\in \mathbf {R} ^{n}\ ,$
• безмасове (фотонне) поле $A_{\mu }(x)\in \mathbf {R} \ ,$ де $\mu =0,1,2,3$ лоренцівські індекси,
• поля Янга-Мілса $A_{\mu }(x)=\sum _{b}A_{\mu }^{b}(x)T_{b}\ ,$ компонентами якого є елементи алгебри Лі компактної групи Лі з генераторами $T_{b}\ .$

З іншого боку, в операторному формалізмі квантової теорії поля є операторами, що діють в гільбертовому просторі. (Якщо точніше, квантові поля $\phi (x)$ є операторозначними розподілами, що означає, що їхні інтеграли $\int f(x)\phi (x)dx$ з відповідними тестовими функціями $f(x)$ є операторами.)

Фізичний зміст теорії поля суттєво залежить від її Лагранжіану ${\mathcal {L}}(\phi (x),\partial ^{n}\phi (x))\ ,$ котрий є функцією $\phi (x)$ і її похідних. Лагранжіан визначає польові рівняння, які складають взаємодії. Якщо сила взаємодії визначається малим параметром $g\ ,$ отже можливо приблизно обчислити фізичні величини до задовільної точності в рамках теорії збурень, завдяки розкладу в степеневий ряд по $g\ .$ Цей метод також застосовується, наприклад, в Квантовій електродинаміці (КЕД), де величина взаємодії пропорційна сталій тонкої структури $\alpha \approx 1/137\ ,$ і багато спостережуваних величин можуть бути отримані у вигляді степеневого ряду по $\alpha \ .$ Але існують випадки, де теорія збурень виявилася незастосовною для обчислень фізичних величин. Найбільш вагомим прикладом є низькоенергетичний режим Квантової хромодинаміки (КХД) - теорія сильної взаємодії елементарних частинок.

Не лише Квантова хромодинаміка, а й інші складові Стандартної Моделі елементарних частинок та її розширення складають низку проблем, що не розв'язуються теорією збурень. Один з важливих кроків до вирішення цієї проблеми був зроблений Вільсоном у 1974 (Wilson, 1974). Він запропонував формулювання Квантової хромодинаміки на просторово-часовій ґратці, що дозволило обійтися без теорії збурень. Ця дискретизація буде описана детальніше далі. Такий підхід дозволив сформулювати ряд чітко визначених математичних проблем які, в принципі, мають розв'язок. Також треба звернути увагу на те, що впровадження просторово-часової ґратки може розглядатися як початкова точка математично чіткого формулювання Квантової теорії поля, відомої як Конструктивна квантова теорія поля.

В сучасній квантовій теорії поля впровадження просторово-часової ґратки є іншим підходом на противагу операторному формалізму. Це теорія поля на ґратці. Її основними складовими є:

• функціональні інтеграли,
• Евклідова теорія поля і
• просторово-часова дискретизація полів.

Теорія поля на ґратці виявилася успішною для обчислень фізичних величин у випадку, де теорія збурень була незастосовною. В цій статті наведено загальний огляд формулювання та методів Теорії поля на ґратці. Основні моменти проілюстровані на прикладі скалярної теорії поля.

Квантова теорія поля у формалізмі функціонального інтегралу

Формулювання Квантової теорії поля через функціональний інтеграл є узагальненням квантово-механічного інтегралу по шляхам. В квантовій механіці амплітуда переходу точкової частинки в одновимірному просторі визначена, як:

$\langle x'|\mathrm {e} ^{-\mathrm {i} HT}|x\rangle ,$

де $|x\rangle$ власний (невідповідний) стан оператора координати і $H$ гамільтонів оператор.

Path of a particle

Амплітуда переходу може бути записана через інтеграл по шляхам, як

$\langle x'|\mathrm {e} ^{-\mathrm {i} HT}|x\rangle =\int \!{\mathcal {D}}x\ \mathrm {e} ^{\mathrm {i} S},$

де інтегрування відбувається по всім можливим траєкторіям $x(t)$ від $x$ до $x'$ протягом часового інтервалу $T\ ,$ див. <figref>Particle-path.png</figref>, і

$S=\int _{0}^{T}\!dt\,L(x,{\dot {x}})$

є класичною дією для такої траєкторії.

Формально, міра інтегралу по шляхам визначається, як

${\mathcal {D}}x\equiv \prod _{t}dx(t)$

з точністю до нормалізаційного фактору. Для частинки в 3-х вимірному просторі цей вираз узагальнюється до траєкторій $x_{i}(t)\ ,$ де $i$=1,2,3, і

${\mathcal {D}}x=\prod _{t}\prod _{i}dx_{i}(t).$

Це, мабуть, найбільш інтуїтивна ілюстрація квантово-механічної амплітуди переходу. Вона визначається як інтеграл вкладів всіх можливих шляхів від початкової до кінцевої точок з ваговим множником визначеним класичною дією, котра відповідає кожному конкретному шляху.

Для більш детального і математично акуратного формулювання інтегралів по шляхам читача може зацікавити підручник (Glimm and Jaffe, 1987).

Формулювання квантової механіки через інтеграл по шляхам може бути використано для побудови теорії поля. Розглянемо скалярне поле $\varphi (x)\ ,$ де $x=({\vec {x}},t)$ позначає просторово-часові координати, а еволюція в часі $\varphi ({\vec {x}},t)$ визначається як

$\varphi ({\vec {x}},t)=\mathrm {e} ^{\mathrm {i} Ht}\,\varphi ({\vec {x}},t=0)\,\mathrm {e} ^{-\mathrm {i} Ht}.$

В теорії поля основними об'єктами досліджень є вакуумні значення дії (часовпорядкованих) польових операторів, функції Гріна:

$\langle 0|\varphi (x_{1})\varphi (x_{2})\dots \varphi (x_{n})|0\rangle ,\qquad t_{1}>t_{2}>\dots >t_{n}.$

Важливим прикладом є пропагатори $\langle 0|\varphi (x)\varphi (y)|0\rangle .$ Функції Гріна, по суті, містять всю фізичну інформацію. Як от елементи S-матриці пов'язані з функціями Гріна, наприклад, у випадку двочастинкового розсіяння можуть бути отримані як $\langle 0|\varphi (x_{1})\dots \varphi (x_{4})|0\rangle .$

Замість того, щоб розглядати квантову теорію поля у формулюванні функціонального інтегралу з самого початку, ми обмежимося переходом від квантово-механічного підходу до теорії поля шляхом розгляду аналогій. В цьому випадку основні змінні $x_{i}(t)$ переходять в поля $\phi ({\vec {x}},t)\ .$ Правила переходів:

${\begin{aligned}x_{i}(t)\quad &\longleftrightarrow \quad \phi ({\vec {x}},t)\\i\quad &\longleftrightarrow \quad {\vec {x}}\\\prod _{t,i}dx_{i}(t)\quad &\longleftrightarrow \quad \prod _{t,{\vec {x}}}d\phi ({\vec {x}},t)\equiv {\mathcal {D}}\phi \\S=\int \!dt\ L\quad &\longleftrightarrow \quad S=\int \!dt\,d^{3}x\ {\mathcal {L}},\end{aligned}}$

де $S$ класична дія.

Для скалярної дії можна розглянути наступний Лагранжіан з густиною ${\begin{aligned}{\mathcal {L}}&={\frac {1}{2}}\left(({\dot {\phi }}(x))^{2}-(\nabla \phi (x))^{2}\right)-{\frac {m_{0}^{2}}{2}}\phi (x)^{2}-{\frac {g_{0}}{4!}}\phi (x)^{4}\\&={\frac {1}{2}}(\partial _{\mu }\phi )(\partial ^{\mu }\phi )-{\frac {m_{0}^{2}}{2}}\phi (x)^{2}-{\frac {g_{0}}{4!}}\phi (x)^{4}.\end{aligned}}$

Маса $m_{0}$ і стала взаємодії $g_{0}$ мають індекс $0\ ,$ оскільки вони є ненормалізованими параметрами. Ця теорія лежить в основі моделі Хіггса-Юкави, де $\phi (x)$ поле Хіггса.

За аналогією до квантово-механічного інтегралу по шляхам, представлення функцій Гріна через функціональний інтеграл визначається як

$\langle 0|\varphi (x_{1})\varphi (x_{2})\dots \varphi (x_{n})|0\rangle ={\frac {1}{Z}}\int \!{\mathcal {D}}\phi \ \phi (x_{1})\phi (x_{2})\dots \phi (x_{n})\ \mathrm {e} ^{\mathrm {i} S}$

де $Z=\int \!{\mathcal {D}}\phi \ \mathrm {e} ^{\mathrm {i} S}.$ Цей інтеграл включає по всім конфігураціям класичних полів.

Як зазначалося вище, ми не вдавалися до виведення функціонального інтегралу, а лише сформулювали його за аналогією. Більше того, ми перейшли від квантово-механічної амплітуди переходу до функції Гріна, що, загалом, має дещо інший зміст.

Формулювання функціонального інтегралу викликає ряд запитань. По-перше, як утворюється проекція на основний стан $|0\rangle$? По-друге, ці інтеграли містять осцилюючі підінтегральні елементи за рахунок уявних експонент; як на рахунок збіжності? На останок, чи можливо їх обчислювати чисельно?

Далі буде викладено як впровадження уявного часу допомагає дати відповідь на ці запитання.

Евклідова теорія поля

В рамках квантової механіки теж можна сформулювати функції Гріна, Наприклад

$G(t_{1},t_{2})=\langle 0|X(t_{1})X(t_{2})|0\rangle ,\qquad t_{1}>t_{2},$

де $X(t)$ оператор координати в представленні Гейзенберга. Далі буде продемонстровано, що такі функції Гріна пов'язані з квантово-механічної амплітудою переходу в комплексному часі через аналітичне продовження. Розглянемо матричний елемент

$\langle x',t'|X(t_{1})X(t_{2})|x,t\rangle =\langle x'|\mathrm {e} ^{-\mathrm {i} H(t'-t_{1})}X\mathrm {e} ^{-\mathrm {i} H(t_{1}-t_{2})}X\mathrm {e} ^{-\mathrm {i} H(t_{2}-t)}|x\rangle$

де $t'>t_{1}>t_{2}>t\ .$ Тепер нехай всі часи розглянемо суто комплексними

$t=-\mathrm {i} \tau ,$

впорядковані як $\tau '>\tau _{1}>\tau _{2}>\tau \ .$ Це приводить до виразу

$\langle x'|\mathrm {e} ^{-H(\tau '-\tau _{1})}X\mathrm {e} ^{-H(\tau _{1}-\tau _{2})}X\mathrm {e} ^{-H(\tau _{2}-\tau )}|x\rangle .$

Підставляючи набір власних станів енергій, отримуємо розклад оператора часової еволюції в комплексному часу

$\mathrm {e} ^{-H\tau }=\sum _{n=0}^{\infty }\mathrm {e} ^{-E_{n}\tau }|n\rangle \langle n|=|0\rangle \langle 0|+\mathrm {e} ^{-E_{1}\tau }|1\rangle \langle 1|+\dots ,$

де основний енергетичний стан віднормовано до $E_{0}=0\ .$ Для великих $\tau$ спрощується до проектора на основний стан. Отже, в границі $\tau '\rightarrow \infty$ і $\tau \rightarrow -\infty$ отримуємо матричний елемент

$\langle x'|0\rangle \langle 0|X\mathrm {e} ^{-H(\tau _{1}-\tau _{2})}X|0\rangle \langle 0|x\rangle ,$

і, аналогічно,

$\langle x'|\mathrm {e} ^{-H(\tau '-\tau )}|x\rangle \longrightarrow \langle x'|0\rangle \langle 0|x\rangle .$

Тому функції Гріна в комплексному часі,

$G_{E}(\tau _{1},\tau _{2})=\langle 0|X\mathrm {e} ^{-H(\tau _{1}-\tau _{2})}X|0\rangle ,$

можуть бути виражені як

$G_{E}(\tau _{1},\tau _{2})=\lim _{\tau '\rightarrow \infty ,\ \tau \rightarrow -\infty }{\frac {\langle x'|\mathrm {e} ^{-H(\tau '-\tau _{1})}X\mathrm {e} ^{-H(\tau _{1}-\tau _{2})}X\mathrm {e} ^{-H(\tau _{2}-\tau )}|x\rangle }{\langle x'|\mathrm {e} ^{-H(\tau '-\tau )}|x\rangle }}.$

Тепер знаменник, як і чисельник, можуть бути сформульовані через інтеграли по шляхам. Відмінність випадку комплексного часу від дійсного полягає в тому, що тепер має використовуватися

$\langle x|\mathrm {e} ^{-H\Delta \tau }|y\rangle \approx {\sqrt {\frac {m}{2\pi \Delta \tau }}}\ \exp -\Delta \tau \left\{{\frac {m}{2}}\left({\frac {x-y}{\Delta \tau }}\right)^{2}+V(x)\right\}.$

Це приводить до формулювання через інтеграли по шляхам

$G_{E}(\tau _{1},\tau _{2})={\frac {1}{Z}}\int \!{\mathcal {D}}x\ x(\tau _{1})x(\tau _{2})\,\mathrm {e} ^{-S_{E}},$

де $\quad Z=\int \!{\mathcal {D}}x\ \mathrm {e} ^{-S_{E}}\ \,$ і

$S_{E}=\int \!d\tau \left\{{\frac {m}{2}}\left({\frac {dx}{d\tau }}\right)^{2}+V(x(\tau ))\right\}.$

Функції Гріна в дійсному часі, в яких було зацікавлення початково, може бути отримано з $G_{E}$ через аналітичне продовження, $G(t_{1},t_{2})=G_{E}(\mathrm {i} t_{1},\mathrm {i} t_{2})\ .$

Wick rotation from imaginary to real time coordinates

Аналітичне продовження має бути виконано таким чином, щоб всі часові аргументи чергувалися одночасно проти годинникової стрілки на комплексній $t$-площині. Це так зване обертання Віка, проілюстроване на <figref>Wickrotation.png</figref>.

Повертаючись до теорії поля, функції Гріна

$G(x_{1},\dots ,x_{n})=\langle 0|T\varphi (x_{1})\dots \varphi (x_{n})|0\rangle ,$

продовжені до комплексного часу, $t=-\mathrm {i} \tau \ ,$ є так званими функціями Швінгера.

$G_{E}(({\vec {x}}_{1},\tau _{1}),\dots ,({\vec {x}}_{n},\tau _{n}))=G(({\vec {x}}_{1},-i\tau _{1}),\dots ,({\vec {x}}_{n},-i\tau _{n})).$

По аналогії з випадком квантової механіки, їх формулювання через функціональні інтеграли

$G_{E}(x_{1},\dots ,x_{n})={\frac {1}{Z}}\int \!{\mathcal {D}}\phi \ \phi (x_{1})\dots \phi (x_{n})\,\mathrm {e} ^{-S_{E}}$

де $Z=\int \!{\mathcal {D}}\phi \ \mathrm {e} ^{-S_{E}}$ і ${\begin{aligned}S_{E}&=\int \!d^{3}xd\tau \left\{{\frac {1}{2}}\left({\frac {d\phi }{d\tau }}\right)^{2}+{\frac {1}{2}}(\nabla \phi )^{2}+{\frac {m_{0}^{2}}{2}}\phi ^{2}+{\frac {g_{0}}{4!}}\phi ^{4}\right\}\\&=\int \!d^{4}x\left\{{\frac {1}{2}}(\partial _{\mu }\phi )^{2}+{\frac {m_{0}^{2}}{2}}\phi ^{2}+{\frac {g_{0}}{4!}}\phi ^{4}\right\}.\end{aligned}}$

Як можна побачити з кінетичної частини $S_{E}\ ,$ метрика Міньковського простору

$-ds^{2}=-dt^{2}+dx_{1}^{2}+dx_{2}^{2}+dx_{3}^{2}$

переходить в

$d\tau ^{2}+dx_{1}^{2}+dx_{2}^{2}+dx_{3}^{2},$

що є Евклідовою метрикою. Саме тому такі об'єкти називають Евклідовими функціями Гріна $G_{E}$ і Евклідовими функціональними інтегралами. Вони є початковими елементами дослідження незбурюваних теорій поля.

Чи можливий подальший розгляд певної теорії поля аналітично при переході від дійсного часу до комплексного і навпаки залежить від ряду необхідних критеріїв. Для великого класу теорій поля ці умови були проаналізовані і сформульовані в роботах Osterwalder і Schrader, див. (Osterwalder and Schrader, 1973, 1975). Наприклад, Евклідова теорія поля має задовольняти так званий критерій reflection positivity щоб мати відповідну теорію поля у просторі Міньковського.

Оскільки $S_{E}$ дійсна, шукані інтеграли тепер дійсні і не мають неприємних осцилюючих складників. Більш того, оскільки $S_{E}$ обмежена знизу, фактор $\exp(-S_{E})$ в інтегралі теж обмежений. Сильно флуктуюючі поля дають велике значення Евклідової дії $S_{E}$, а отже подавляють фактор $\exp(-S_{E})\ .$ (Строго кажучи, це твердження не має сенсу якщо не брати до уваги ренормалізацію.) Це робить Евклідові інтеграли більш привабливими в порівнянні з Міньковськими аналогами.

Можна подумати, що в Евклідовому формалізмі всі величини нефізичні і немає можливості напряму отримати фізичні результати з Евклідових функцій Гріна. Але це не так. Наприклад, спектр в рамках теорії може бути отриманий наступним чином. Розглянемо вакуумне значення в наступному вигляді

$\langle 0|A_{1}\mathrm {e} ^{-H\tau }A_{2}|0\rangle ,$

де $A_{i}$ формуються з полів $\varphi \ ,$ наприклад як, $A=\varphi ({\vec {x}},0)$ або $A=\int \!d^{3}x\ \varphi ({\vec {x}},0)\ .$ Тепер підставимо набір власних станів енергій, аналогічно до раніше розглянутого випадку,

$\langle 0|A_{1}\mathrm {e} ^{-H\tau }A_{2}|0\rangle =\sum _{n}\langle 0|A_{1}|n\rangle \mathrm {e} ^{-E_{n}\tau }\langle n|A_{2}|0\rangle .$

У випадку неперервного спектру суму варто трактувати як інтеграл. З іншого боку формулювання вакуумного значення через функціональний інтеграл приводить до

${\frac {1}{Z}}\int \!{\mathcal {D}}\phi \ \mathrm {e} ^{-S_{E}}A_{1}(\tau )A_{2}(0)=\sum _{n}\langle 0|A_{1}|n\rangle \langle n|A_{2}|0\rangle \mathrm {e} ^{-E_{n}\tau }.$

Це схоже до проекції основного стану на початку цього розділу. Для великих $\tau$ найнижчий власний стан енергії домінуватиме в сумі. Отже можна отримати нижній регіон спектру з асимптотичної поведінки цього вакуумного значення. Обираючи $A_{1},A_{2}$ зручним чином, наприклад, для

$A\equiv A_{1}=A_{2}=\int \!d^{3}x\ \varphi ({\vec {x}},0),$

таким щоб $\langle 0|A|1\rangle \neq 0$ для одночастинкового стану $|1\rangle$ з нульовим моментом ${\vec {p}}=0$ і масою $m_{1}\ ,$ можна отримати

${\frac {1}{Z}}\int \!{\mathcal {D}}\phi \ \mathrm {e} ^{-S_{E}}A(\tau )A(0)=|\langle 0|A|1\rangle |^{2}\mathrm {e} ^{-m_{1}\tau }+\dots ,$

що робить можливим знаходження маси частинки.

Надалі ми продовжуватимемо працювати в Евклідовому просторі і знехтуємо індексом $E\ ,$ маючи на увазі $S\equiv S_{E}$ Евклідову дію.

Дискретизація на ґратці

Одне з основних питань все ще залишається відкритим: чи інтегрування нескінченної розмірності по всім конфігураціям класичних полів, наприклад

${\mathcal {D}}\phi =\prod _{x}d\phi (x),$

має сенс взагалі? Як воно визначене?

3-dimensional lattice

В квантовій механіці формулювання інтегралу по шляхам може бути виведено як границя дискретизації в часі. В теорії поля, поля залежать від чотирьох Евклідових координат замість однієї часової і трьох просторових координати. Ми можемо ввести дискретизований прострі-час у вигляді ґратки, наприклад гіперкубічної, визначеної як

$x_{\mu }=an_{\mu },\qquad n_{\mu }\in \mathbf {Z} ,$

див. <figref>3d-lattice.png</figref>.

Величина $a$ називається періодом ґратки. Варто звернути увагу, що, будучи розмірним параметром, період ґратки не є параметром дискретизованої теорії, котрий міг би бути використаний в комп'ютерній програмі для обчислення інтегралів по шляхам. Період ґратки у фізичних одиницях є обчислюваною величиною, що визначається динамікою моделі. Це буде детальніше описано в секції "Неперервна границя".

Скалярне поле

$\phi (x),\qquad x\in {\mbox{lattice}},$

тепер визначене лише на вузлах ґратки. Частинні похідні заміняються скінченними різницями

$\partial _{\mu }\phi \longrightarrow \Delta _{\mu }\phi (x)\equiv {\frac {1}{a}}(\phi (x+a{\hat {\mu }})-\phi (x)),$

а просторово-часові інтеграли заміняються сумуванням:

$\int \!d^{4}x\quad \longrightarrow \sum _{x}a^{4}\ .$

Дія дискретизованої $\phi ^{4}$-теорії може бути записана як

$S=\sum _{x}a^{4}\left\{{\frac {1}{2}}\sum _{\mu =1}^{4}(\Delta _{\mu }\phi (x))^{2}+{\frac {m_{0}^{2}}{2}}\phi (x)^{2}+{\frac {g_{0}}{4!}}\phi (x)^{4}\right\}.$

У функціональному інтегралі міра ${\mathcal {D}}\phi$ тепер включає лише вузли ґратки $x$. Отже набір дискретних змінних має інтегруватися. Якщо ґратка є скінченною тоді й інтеграл має скінченну розмірність.

Дискретизація простору-часу з використанням ґратки має один дуже важливий наслідок. Завдяки ненульовому періоду ґратки виникає відсічення в просторі імпульсів. Воно може спостерігатися при Фур'є перетворенні поля

${\tilde {\phi }}(p)=\sum _{x}a^{4}\ \mathrm {e} ^{-\mathrm {i} px}\ \phi (x).$

Фур'є образ поля періодичний в просторі імпульсів, щоб можна було визначити

$p_{\mu }\cong p_{\mu }+{\frac {2\pi }{a}}$

і обмежити імпульс до першої зони Брілюена

$-{\frac {\pi }{a}}\,<\,p_{\mu }\,\leq {\frac {\pi }{a}}.$

Обернене перетворення Фур'є, наприклад, визначається як

$\phi (x)=\int _{-\pi /a}^{\pi /a}{\frac {d^{4}p}{(2\pi )^{4}}}\ \mathrm {e} ^{\mathrm {i} px}\ {\tilde {\phi }}(p).$

Тут можна впізнати ультрафіолетове відсічення

$|p_{\mu }|\leq {\frac {\pi }{a}}.$

Отже теорії поля на ґратці регуляризується природнім чином.

Почнемо з розгляду скінченної ґратки для збереження чітких означень. Розглянемо гіперкубічну ґратку з періодом $L_{1}=L_{2}=L_{3}=L$ в кожній просторовій розмірності і з довжиною $L_{4}=T$ в Евклідовому часі,

$x_{\mu }=an_{\mu },\qquad n_{\mu }=0,1,2,\dots ,L_{\mu }-1,$

з скінченним об'ємом $V=L^{3}T\ .$ В скінченному об'ємі необхідно визначити крайові умови. Популярним є вибір періодичної крайової умови

$\phi (x)=\phi (x+aL_{\mu }\,{\hat {\mu }}),$

де ${\hat {\mu }}$ одиничний вектор в $\mu$-напрямку. Це означає, що й імпульс дискретний,

$p_{\mu }={\frac {2\pi }{a}}\,{\frac {l_{\mu }}{L_{\mu }}}\qquad {\mbox{with}}\ l_{\mu }=0,1,2,\dots ,L_{\mu }-1,$

а отже й інтегрування в імпульсному просторі також замінюється скінченним сумуванням

$\int \!{\frac {d^{4}p}{(2\pi )^{4}}}\ \longrightarrow \ {\frac {1}{a^{4}L^{3}T}}\sum _{l_{\mu }}.$

Таким чином всі функціональні інтеграли перетворилися в регуляризовані і скінченні вирази.

Звичайно, тепер необхідно відновити фізику в неперервному та нескінченному просторі-часі. Для цього необхідно перейти до границі нескінченного об'єму,

$L,T\longrightarrow \infty ,$

і до неперервної границі,

$a\longrightarrow 0.$

Побудова неперервної границі в теорії поля на ґратці зазвичай дуже нетривіальна задача.

Формулювання Евклідової квантової теорії поля на ґратці несе корисну аналогію зі статистичною механікою. Функціональний інтеграл має вигляд статистичної суми, а отже можна отримати наступні відповідності:

Евклідова теорія поля	Статистична механіка
Генеруючий функціонал $\int \!{\mathcal {D}}\phi \ \mathrm {e} ^{-S}$	Статистична сума $\sum \mathrm {e} ^{-\beta {\mathcal {H}}}$
Дія $S$	Гамільтоніан $\beta {\mathcal {H}}$
Маса $m$ $G\sim \mathrm {e} ^{-mt}$	Обернена кореляційна довжина $1/\xi$ $G\sim \mathrm {e} ^{-{\frac {x}{\xi }}}$

Ця формальна аналогія дозволяє використовувати добре відомі методі статистичної механіки в теорії поля і навпаки. Навіть термінологія обох дисциплін часто збігається. Наприклад, в теорії поля часто застосовується високотемпературний розклад або наближення середнього поля, а в статистичний механіці ренормалізаційна група.

Гамільтонова теорія поля на ґратці

Альтернативою Евклідовій теорії поля на ґратці є Гамільтонова теорія, запропонована Kogut і Susskind (Kogut and Susskind, 1975). В цьому формулюванні лише трьохвимірний простір координат дискретизується на ґратці, а інші координати залишаються неперервними. Більш того, час залишається дійсним і не продовжується до евклідового інтервалу. Гамільтонова теорія дозволяє застосовувати деякі аналітичні методи такі, як розклад сильної взаємодії або теорія збурень. На жаль, чисельний метод Монте-Карло незастосовний в рамках даної теорії через, що вона втратила свою популярність.

Калібрувальна теорія на ґратці

Теорії калібрувальних полів також можуть бути сформульовані на просторов-часовій ґратці. Більш детально це описано в окремій статті, тут лише буде наведено основні елементи для калібрувальної групи SU(N).

Link $b=\langle x,y\rangle$ between lattice points $x$ and $y$

Шляхи, що з'єднують найближчі точки на ґратці називаються лінками.

Кожному лінку

$b=\langle x+a{\hat {\mu }},x\rangle$

у напрямку ґратки ${\hat {\mu }}$ відповідає змінна лінку

$U(b)\equiv U(x+a{\hat {\mu }},x)\equiv U_{x\mu }\in {\rm {SU}}(N)$

Група цих змінних представляє калібрувальне поле. Калібрувальне поле дискретної алгебри Лі $A_{\mu }^{b}(x)$ може бути введено як

$U_{x\mu }\equiv \exp\{\mathrm {i} g_{0}aA_{\mu }^{b}(x)T_{b}\},$

де $T^{a}$ генератори калібрувальної групи, а $g_{0}$ стала зв'язку.

Plaquette $p$ in lattice directions $\mu$ and $\nu$

Найменші замкнуті шляхи га ґратці - плакетки, як показано на <figref>Lattice-plaquette.png</figref>.

Змінні плакеток

$U(p)=U_{x\mu \nu }\equiv U_{(x+a{\hat {\nu }})(-\nu )}U_{(x+a{\hat {\mu }}+a{\hat {\nu }})(-\mu )}U_{(x+a{\hat {\mu }})\nu }U_{x\mu }$

входять в дію Вільсона'

$S_{W}=-\sum _{p}{\frac {2}{g_{0}^{2}}}{\rm {Re}}({\rm {Tr}}(U(p))).$

В неперервній границі, де $a$ прямую до нуля, можна отримати

$S_{W}={\frac {g_{0}^{2}}{4}}\sum _{x}a^{4}F_{\mu \nu }^{b}F_{\mu \nu }^{b}+{\mathcal {O}}(a^{5}),$

що відтворює дію Янга-Мілса.

Інтеграл по всім конфігураціям калібрувальних полів на ґратці становить інтеграл по всім змінним лінків $U(b)\ .$ Отже, для вакуумного значення довільної спостережуваної величини $A$ можна записати

$\langle A\rangle ={\frac {1}{Z}}\int \!\prod _{b}dU(b)\ A\ \mathrm {e} ^{-S_{W}},$

де інтегрування $dU(b)$ для даного лінка $b$ має розумітися, як інваріантне інтегрування по групі многовиду нормованої як

$\int \!dU=1.$

Методи

В попередніх розділах були визначені функціональні інтеграли для теорій поля на ґратці. Іншим аспектом є методи обчислення цих інтегралів високої розмірності. Точні обчислення даних виразів в загальному випадку не є можливим. В цьому розділі будуть наведені деякі методи наближених обчислень функціональних інтегралів.

Теорія збурень

Хоча теорія поля на ґратці дає можливість вивчати елементи теорії без застосування теорії збурень, це все ще дуже цінний інструмент і для використання в теоріях на ґратках. Особливо для порівняння результатів збурених та незбурених теорій в межах де обидва методи можуть бути застосовні.

Застосування теорії збурень полягає в степеневому розкладі по сталій зв'язку. Ґратка забезпечує природне ультрафіолетове відсічення $\pi /a$ для всіх імпульсів. Також можна помітити, що пропагатори та вершини відрізняються від відповідних аналогів у неперервному випадку, що забезпечується особливою формою дії на ґратці. Наприклад, з'являється самодія глюона всіх порядків, не лише трьох- та чотирьох-глюонна вершини.

Розклад сильної взаємодії

Аналогія між Евклідовою теорією поля і статистичною механікою була встановлена. В статистичній механіці давно відомий метод високотемпературного розкладу. Для калібрувальних теорій на ґратці це розклад по степеням

$\beta \sim {\frac {1}{g_{0}^{2}}},$

що є малою величиною при великих значеннях сталої взаємодії $g_{0}\ .$ А отже це й справді розклад по сильній взаємодії. Больцманівський фактор розкладається як

$\exp {\left(\beta {\frac {1}{N}}{\rm {Re}}({\rm {Tr}}(U(p)))\right)}=1+\beta {\frac {1}{N}}{\rm {Re}}({\rm {Tr}}(U(p)))+\dots \ .$

Результуючий розклад можна зобразити діаграмою, аналогічно до фейнманівських діаграм теорії збурень. Але діаграмними елементами є плакетки $p$ на ґратці. Кожна додаткова степінь $\beta$ додає ще одну плакетку.

У випадку скалярного поля, відповідним методом є стрибкоподібне розширення параметрів, що відповідає розкладу за параметром $\kappa \ ,$ котрий є малим при великих значеннях маси $m_{0}\ .$

Розклад сильної взаємодії і метод стрибкоподібного розширення параметрів мають скінченний радіус збіжності на відміну від теорії збурень, яка, зазвичай, розбіжна і в кращому випадку асимптотична.

Інші аналітичні методи

Інші аналітичні методи доступні для приблизного обчислення функціональних інтегралів калібрувальних теорій на ґратці. Деякі з них:

• наближення середнього поля
• ренормалізаційна група
• ${\frac {1}{N}}$-розклад

Методи Монте-Карло

На скінченній ґратці обчислення вакуумних середніх значень потребує обчислення інтегралів скінченної розмірност. Це відразу підштовхує до застосування чисельних методів. Наївно буде запропонувати інтегрування квадратурами, адже типова ґратка з 40 вузлами матиме $4\cdot 40^{4}$ змінних лінків. Для калібрувальної групи SU(3) це становитиме 81,920,000 справжніх змінних. Це непідйомна задача для методу інтегрування квадратурами навіть в майбутньому. Отже необхідні статистичні методи. Утворення випадкової конфігурації калібрувальної ґратки виявилося дуже неефективним. Важливою є концепція вибірки важливого: для заданої дії ґратки $S$ квадратурні точки $x_{i}$ генеруються з імовірністю

$p(x_{i})\sim \exp\{-S(x_{i})\}.$

Це забезпечує велику кількість точок у необхідних межах інтегрування, значно підвищуючи точність.

У випадку калібрувальної теорії на ґратці квадратурнимим точками є конфігурації $U^{(i)}=\left\{U_{x\mu }^{(i)}\right\}\ .$ Вакуумні середні

$\langle 0|A|0\rangle ={\frac {1}{Z}}\int \!{\mathcal {D}}U\ A(U)\ \mathrm {e} ^{-S(U)}$

отримуються чисельно наближено як середнє

${\bar {A}}\equiv {\frac {1}{n}}\sum _{i=1}^{n}A(U^{(i)}).$

Метод Монте-Карло полягає в створенні послідовності конфігурацій $U^{(1)}\rightarrow U^{(2)}\rightarrow U^{(3)}\rightarrow \dots$ з відповідними ймовірностями статистичним шляхом. Це звісно є відбувається на комп'ютері. Оновлення це - крок , при якому змінюється лише єдина змінна лінку $U_{x\mu }$, тоді як свіп означає прохід через всю ґратку оновлюючи кожну змінну лінку. Поширена техніка для оновлення лінку є алгоритм метрополісу.

Важливою особливістю статистичного підходу до обчислень є існування статистичних похибок. Результат таких обрахунків зазвичай подається у вигляді

$\langle A\rangle ={\bar {A}}\pm \sigma _{\bar {A}},$

де розкид ${\bar {A}}$ зменшується зі зростанням числа $n$ конфігурацій

$\sigma _{\bar {A}}\sim {\frac {1}{n^{1/2}}}.$

Джерела похибок

Результати, отримані методом Монте-Карло, відрізняються від бажаних фізичних результатів наявність похибок різного характеру. Основними причинами похибок є

• наявність статистичних похибок: за рахунок скінченної кількості конфігурацій в симуляціях, $\sim 1/n^{1/2}\ ,$
• ефекти ґратки: за рахунок скінченного періоду ґратки $a\ ,$ часто $\sim a$ or $a^{2}\ ,$
• ефекти об'єму: за рахунок скінченного об'єму ґратки, часто $\sim 1/L\ ,$ $1/L^{2}\ ,$ or $\mathrm {e} ^{-mL}\ ,$
• великі кваркові маси: $m_{q}$ завеликі для Монте-Карло симуляцій,
• применшене наближення: $\det Q=1\ ,$ нехтування динамікою ферміонів.

Неперервна границя

Оскільки обчислення можна проводити лише на ґратці скінченного періоду, важливо мати можливість перейти до неперервної границі. Оскільки період ґратки є регулятором теорії, доцільно використати ренормалізаційну групу. Знаючи функціональну залежність сталої взаємодії $g_{0}$ від регулятора, іншими словами ренормалізаційне рівняння, можна варіювати сталу взаємодії теорії таким чином, щоб досягти неперервної границі.

2-dimensional lattices with increasing correlation lengths $\xi$

В неперервній границі період ґратки $a$ має прямувати до нуля, при цьому фізичні маси $m$ мають прямувати до скінченного значення. Період ґратки не є безрозмірною величиною, а отже необхідно зафіксувати масштаби мас $m\ ,$ наприклад маси деяких частинок і розглянути границю $am\rightarrow 0\ .$ Обернена величина

${\frac {1}{am}}\equiv \xi ,$

може вважатися кореляційною довжиною. В неперервній границі $\xi$ має прямувати до нескінченності, що називається критичною точкою теорії. В <figref>Lattice-spacings.png</figref> це проілюстровано на двовимірній ґратці з різними кореляційними довжинами.

В чистій калібрувальні теорії є єдина безрозмірна стала взаємодії $g_{0}$ і $am$ i справді функція $g_{0}\ .$ Для того щоб підійти до неперервної границі, необхідно варіювати $g_{0}$ таким чином, щоб $am\rightarrow 0\ .$ Як це відбувається визначається ренормалізаційною групою:

$-a{\frac {\partial g_{0}}{\partial a}}=\beta _{LAT}(g_{0})=-\beta _{0}g_{0}^{3}-\beta _{1}g_{0}^{5}+\dots ,$

де перший член розкладу

$\beta _{0}={\frac {11}{3}}\,N\,{\frac {1}{16\pi ^{2}}}.$

В збуреному режимі $g_{0}$ вказує на те, що при спаданні $am$ стала взаємодії $g_{0}$ також спадає, прямуючи до нуля. Тому неперервна границя пов'язана з границею

$g_{0}\rightarrow 0\qquad {\mbox{(continuum limit).}}$

Розв'язок ренормгрупового рівняння до другого порядку $g_{0}$ є

$a=\Lambda _{LAT}^{-1}\ \exp \left(-{\frac {1}{2\beta _{0}g_{0}^{2}}}\right)\ (\beta _{0}g_{0}^{2})^{-{\frac {\beta _{1}}{2\beta _{0}^{2}}}}\ \{1+{\mathcal {O}}(g_{0}^{2})\},$

де з'являється $\Lambda$-параметр ґратки $\Lambda _{LAT}$. Розв'язок для $g_{0}$ дає

$g_{0}^{2}={\frac {-1}{\beta _{0}\log {a^{2}\Lambda _{LAT}^{2}}}}+\dots ,$

що також демонструє стрімке спадання $g_{0}$ в неперервній границі:

$g_{0}^{2}\rightarrow 0\quad {\mbox{for}}\ a\rightarrow 0.$

Також можна помітити, що

$am=C\,\exp \left(-{\frac {1}{2\beta _{0}g_{0}^{2}}}\right)\cdot (\dots ),$

що свідчить про походження маси $m$ з незбуреної теорії.

Ці міркування, засновані на збуреній $\beta$-функції, підштовхують до гіпотези: неперервна границя калібрувальної теорії на ґратці має розглядатися при $g_{0}\rightarrow 0\ .$

Отже сценарій переходу до неперервної границі наступний. Обчислити маси у величинах ґратки, наприклад величинах $am\ ,$ і зменшуючи $g_{0}\ ,$ досягнути області, де безрозмірні величини $am$ поводяться як попередньому рівнянні, що називається асимптотичним масштабуванням.

Для відношення мас можна показати, що експоненційна залежність від $1/g_{0}^{2}$ зникає і, можливо, поблизу неперервної границі

${\frac {m_{1}}{m_{2}}}={\mbox{const.}}\times (1+{\mathcal {O}}(a^{p}))$

для деяких цілих $p\ .$ Така поведінка, $m_{1}/m_{2}\approx$ const., називається масштабуванням. В чисельних симуляціях масштабування різних фізичних величин було визначено для калібрувальних теорій поля, квантової хромодинаміки на ґратці та інших моделей, а от підтвердження асимптотичного масштабування набагато більш вимогливе.

Посилання

• http://www.scholarpedia.org/article/Lattice_quantum_field_theory
• Wilson, K (1974). "Confinement of Quarks". Phys. Rev. (D 10): 2445.
• Kogut, J; Susskind, L (1975). "Hamiltonian formulation of Wilson's lattice gauge theories". Phys. Rev. (D 11): 395.
• Osterwalder, K; Schrader, R (1975) [1973]. "Axioms for Euclidean Green's functions". Comm. Math. Phys. (31 (1973)). 83; 42 (1975): 281.

Рекомендовано

• Creutz, Michael (1985). Quarks, Gluons and Lattices. Cambridge: Cambridge University Press. Search this book on
• Montvay, Istvan; Münster, Gernot (1994). Quantum Fields on a Lattice. Cambridge: Cambridge University Press. Search this book on
• Rothe, H J (2005). Lattice Gauge Theories: An Introduction. Singapore: World Scientific. Search this book on
• Smit, J (2002). Introduction to Quantum Fields on a Lattice. Cambridge: Cambridge University Press. Search this book on
• Creutz, Michael, ed. (1992). Quantum Fields on the Computer. Singapore: World Scientific. Search this book on
• DeGrand, T; DeTar, C (2006). Lattice Methods for Quantum Chromodynamics. Singapore: World Scientific. Search this book on
• Gattringer, C; Lang, C B (2010). Quantum Chromodynamics on the Lattice: An Introductory Presentation. Berlin: Springer. Search this book on
• Glimm, James; Jaffe, Arthur (1987). Quantum Physics: A Functional Integral Point of View. Berlin: Springer. Search this book on

Зовнішні посилання

G. Münster, M. Walzl: Lattice gauge theory - a short primer, Proceedings of the Summer School on Phenomenology of Gauge Interactions, August 13-19, 2000, Zuoz, Switzerland, ed. D. Graudenz, V. Markushin. arXiv:hep-lat/0012005

R. Gupta: Introduction to Lattice QCD, Lectures given at the LXVIII Les Houches Summer School "Probing the Standard Model of Particle Interactions", July 28-Sept 5, 1997. arXiv:hep-lat/9807028

Kenway, R. D. (19 Oct 2000). "Lattice Field Theory". arXiv:hep-ph/0010219.

Morningstar, Colin (20 Feb 2007). "The Monte Carlo method in quantum field theory". arXiv:hep-lat/0702020.

U.-J. Wiese: An Introduction to Lattice Field Theory.

Також

• Path integral
• Gauge theories
• Lattice gauge theories
• Quantum electrodynamics
• Quantum chromodynamics
• Asymptotic freedom
• Lattice chiral fermions
• Renormalization
• Renormalization group

References

This article "Квантова теорія поля на гратці" is from Wikipedia. The list of its authors can be seen in its historical and/or the page Edithistory:Квантова теорія поля на гратці. Articles copied from Draft Namespace on Wikipedia could be seen on the Draft Namespace of Wikipedia and not main one.

📰 Article(s) of the same category(ies)

• AD Aerodrom Nikola Tesla
• Studypdf
• 高等华人
• Marpu Padmanabham
• Eurázsiai Gazdasági Közösség
• Samantha Gilabert
• محمد کاشانکی

---