Dexter's story

This series is divergent; therefore, we may be able to do something with it. -- Oliver Heaviside

$\frac{1}{r}$꼴을 갖는 Coulomb potential은 IR 발산이 있는 것으로 유명하다. 좀 더 구체적으로 말하자면, 학부 역학 수준에서 계산할 수 있는 궤도방정식을 풀어 얻는 Rutherford scattering의 미분단면적(differential cross-section)을 계산할 경우 다음과 같은 $\sin^{-4} (\theta/2)$의 꼴을 갖는다는 것이 알려져 있다.

$$ \frac{d\sigma}{d\Omega} \propto \frac{1}{\sin^4 (\theta/2)} $$

이 식을 적분하여 얻는 총산란단면적(total cross-section)은 발산한다.

$$ \sigma_{\text{tot}} = \int \frac{d \sigma}{d \Omega} d \Omega \propto \int \frac{d(\cos \theta)}{\sin^4 (\theta/2)} \to \infty$$

양자역학에서 Coulomb potential이 주어졌을 때의 산란문제를 풀 때도 이 성질과 관련된 현상이 나타난다. Griffiths 양자역학에서는 Coulomb potential을 Yukawa potential의 질량이 없는 극한으로 생각하기 때문에 등장하지 않지만 Landau 3권이나 교수님 세대의 메인 레퍼런스(...)란 느낌이 있는 Shiff책을 뒤적이다 보면 asymptotic region에서 파동함수가 평면파인 $e^{ikz}$로 수렴하는 것이 아니라 로그가 붙은 추가적인 위상항(phase factor)이 등장하는 것을 볼 수 있다.

$$ \psi \sim e^{ikz + (i/k) \log [k(r-z)]} $$

교재에서는 이런 Coulomb potential의 IR 발산에 대해 'Coulomb potential이 장거리 상호작용(long-range interaction)이기 때문에 발생한다'는 설명을 써놓지만, 구체적으로 무한원점에서 0으로 수렴하는 다른 potential들과 어떻게 다른지에 대해 설명하는 경우는 드물다^[각주:1]. 왜 이런 현상이 일어나는지 고전역학적으로 이해하는 것이 이 포스트의 목표.

---

Coulomb potential이 주어졌을 때 그 potential을 따라 움직이는 시험 입자(test particle)의 궤도방정식을 푸는 문제는 몇 안 되는 정확하게 풀 수 있는 고전역학 문제이다. 심지어 궤도방정식 위키백과 페이지가 있을 정도. 시간에 대한 거리의 미분방정식을 각도에 대한 거리의 미분방정식으로 바꾼 뒤 $u = 1/r$이란 변수변환으로 조화진동자 방정식으로 바꾸는 과정이나 이렇게 얻은 궤도방정식으로부터 충돌 파라메터(impact parameter)에 대한 산란각(scattering angle)의 방정식을 얻는 과정은 많은 교재에서 충분히 다루고 있으니 여기서는 생략하기로 하자^[각주:2].

여기서는 eikonal 근사의 변종으로 Coulomb potential에서의 산란을 풀어보자. Eikonal은 기하광학에서 빛의 경로를 계산하기 위해 쓰는데, WKB 근사라고 생각해도 좋다. 여담으로 eikonal은 해밀턴이 기하광학을 풀기 위한 수학적 기법을 다듬으면서 같은 기법이 고전역학에도 적용될 수 있음을 알아차리면서 현재의 해밀턴역학과 심플렉틱기하를 만들어내는 계기가 되었고, 슈뢰딩거의 파동방정식은 기하광학의 eikonal 방정식에서 영감을 얻었다고 한다.

고전역학이든 양자역학이든 산란 문제에서 eikonal 근사란 '직선 근사'라고 생각하면 된다^[각주:3]. 구체적으로 이야기한다면, 입자의 경로를 1) 아무런 산란이 없는 직선 경로에 2) 산란을 일으키는 포텐셜의 효과를 집어넣어 얼마나 직선 경로에서 벗어나는지 섭동계산으로 구하는 방법이 되겠다.

이제 Coulomb potential에서의 고전적인 산란 문제에 eikonal 근사를 적용해보자. Landau 1권에서는 뉴턴역학을 기반으로 eikonal 근사를 사용하지만 여기서는 해밀턴역학을 기반으로 eikonal 근사를 써보기로 한다^[각주:4]. 먼저 해밀토니안을 다음과 같이 적는다.

$$ H = \frac{p^2}{2} - \frac{k}{r} $$

해밀턴 운동방정식은 금방 적을 수 있다.

$$ \dot{\vec{r}} = \{ H , \vec{r} \} = \vec{p} \,,\, \dot{\vec{p}} = \{ H , \vec{p} \} = - \frac{k \vec{r}}{r^3} $$

이 역학계의 산란문제를 eikonal 근사로 푸는 것은 다음과 같은 ansatz를 이용해 섭동전개 파라메터 $k$에 대해 푸는 것으로 생각할 수 있다.

$$ \vec{p} = \vec{p}_0 + k \vec{p}_1 (t) + k^2 \vec{p}_2 (t) + \cdots \,,\, \vec{r} = \left( \vec{b} + \vec{p}_0 t \right) + k \vec{r}_1 (t) + k^2 \vec{r}_2 (t) + \cdots $$

여기서 $\vec{p}_0$는 asymptotic region에서의 운동량이고, $\vec{b}$는 충돌 파라메터의 역할을 한다. 이렇게 해석하려면 $\vec{b} \cdot \vec{p}_0 = 0$이란 조건을 추가로 얹어주는 것이 좋다. 섭동이 없는 원래 경로에서 시간 $t$의 원점을 재정의하는 것으로 이 조건을 맞출 수도 있고.

이제 위의 방정식을 풀어보자. 방정식을 풀려면 경계조건을 줘야 하는데, 가장 먼저 생각할 수 있는 경계조건은 다음 경계조건이다.

$$\vec{r}_{i>0} (-\infty) = \vec{p}_{i>0} (-\infty) = 0$$

언듯 보기에는 문제가 없는 경계조건으로 보인다. $t = -\infty$는 산란이 일어나기 한참 전의 과거이므로 섭동이 없는 원래 경로와 일치해야 한다는 직관과도 맞고. 하지만 이 경계조건은 절대로 맞춰줄 수 없다. Coulomb potential의 꼬리가 너무 길기 때문. 우선 이 문제를 무시하고 그냥 방정식을 풀어보자.

$\vec{p}_1$에 대한 운동방정식은 다음과 같이 주어진다.

$$ k \dot{\vec{p}}_1 (t) = - \frac{k (\vec{b} + \vec{p}_0 t)}{(b^2 + p_0^2 t^2)^{3/2}} $$

이 식에 처음 얹은 경계조건을 넣고 풀면 다음과 같은 답을 얻는다.

$$ \vec{p}_1 (t) = - \int_{-\infty}^t \frac{\vec{b} + \vec{p}_0 \tau}{(b^2 + p_0^2 \tau^2)^{3/2}} d\tau = -\frac{1}{ab^2} \left[ \left( 1 + \frac{at}{\sqrt{1 + a^2 t^2}} \right) \hat{b} - \frac{\hat{a}}{\sqrt{1 + a^2 t^2}} \right] $$

쌍곡함수로 변수변환을 하면 적분을 쉽게 할 수 있다. 문제를 풀 때 새로 정의한 변수들은 다음과 같다.

$$ \hat{b} := \frac{\vec{b}}{b} \,,\, \vec{a} := \frac{\vec{p}_0}{b} \,,\, \hat{a} := \frac{\vec{a}}{a} = \frac{\vec{p}_0}{p_0} $$

$k^1$ 차수에서 운동량 변화는 단순히 $\vec{p}_1 (+\infty)$를 읽어내면 된다.

$$\Delta \vec{p}_1 := \vec{p}_1 (+\infty) = - \frac{2 \hat{b}}{ab^2} = - \frac{2 \vec{b}}{p_0 b^2}$$

마찬가지로 $k^2$ 차수에서 운동량 변화는 $\vec{p}_2 (+\infty)$를 읽어내면 되는데, $\vec{p}_2$는 $\vec{r}_1$에 대한 해가 있어야 풀 수 있다^[각주:5].

$$k^2 \dot{\vec{p}}_2 = - k^2 \left[ \frac{\vec{r}_1}{r_0^3} - \frac{3 \vec{r}_0 (\vec{r}_0 \cdot \vec{r}_1)}{r_0^5} \right]$$

따라서 $\vec{r}_1(t)$를 풀어야 한다. 우선 식을 적어보자.

$$\vec{r}_1 (t) = \int_{-\infty}^{t} \vec{p}_1 (\tau) d\tau = - \frac{1}{ab^2} \int_{-\infty}^{t} \left[ \left( 1 + \frac{a\tau}{\sqrt{1 + a^2 \tau^2}} \right) \hat{b} - \frac{\hat{a}}{\sqrt{1 + a^2 \tau^2}} \right] d\tau$$

눈치가 빠른 분들은 알아차리셨겠지만, 이 정적분은 잘 정의되질 않는다. 두번째 항이 $\sim \tau^{-1}$의 꼴을 하고 있기 때문에 무한대에서 로그 발산이 있기 때문이다. 첫번째 항은 정적분으로 처리하고 두번째 항은 정적분을 포기하고 부정적분으로 처리할 경우 다음 식을 얻는다.

$$\vec{r}_1 (t) = - \frac{ e^{\sinh^{-1} (at)}}{a^2 b^2} \hat{b} + \left. \frac{\sinh^{-1}(at)}{a^2 b^2} \hat{a} \right|_{-\infty}^{t}$$

$x \in \mathbb{R}$일 때 $\sinh^{-1} x = \log (x + \sqrt{1+x^2})$이므로, 두번째 항의 발산은 예상대로 로그 발산임을 확인할 수 있다. 이 로그 발산은 다음과 같이 이해할 수 있다. Coulomb potential에서의 에너지 보존을 생각하면 무한대에서의 입자의 속력을 $v$라고 할 때 asymptotic region에서의 입자의 속력 $v$는 다음과 같다.

$$ \frac{v^2}{2} = \frac{v_0^2}{2} + \frac{k}{r} \Rightarrow v \sim v_0 + \frac{c}{r}$$

따라서 아무런 힘을 못 느끼고 $v_0$의 속력으로 이동하는 섭동이 없는 경로와 Coulomb potential의 영향을 받아 섭동이 있는 경로 사이의 변위(displacement)를 계산하면 다음과 같아진다.

$$ \Delta r \sim \int (v - v_0) dt \sim \int \frac{c}{r} dt \sim \int \frac{1}{dr/dt} \frac{c}{r} dr \sim \frac{c}{v_0} \log r $$

$r^{-1}$보다 빠르게 떨어지는 다른 potential의 경우 입자가 멀어져 가면서 potential로부터 받는 영향이 충분히 빠르게 줄어들어 섭동이 없는 경로와 potential의 영향을 받은 경로 사이의 변위가 일정하게 유지된다. 하지만 Coulomb potential의 경우 potential의 영향이 0으로 줄어드는 속도가 느려 아무리 멀어지더라도 변위의 차이가 계속 누적되는 것이다. 발산하는 총산란단면적이나 양자역학 산란 문제를 풀 때 평면파에 로그만큼의 위상항이 추가로 붙는 현상은 이 흔적이라고 이해할 수 있다.

---

여튼, $k^2$ 차수의 운동량 변화를 계산하는 문제로 돌아오자. 발산이 있으면 잡으면 되는 법이다.

가장 단순한 해법은 $t = - \infty$를 기준점으로 잡지 않고 $t = 0$를 기준점으로 잡는 것이다. 실제로 worldline quantum field theory(WQFT)를 도입해서 post-Minkowskian 계산을 하는 팀에서 이런 접근을 취하고 있는데, 이 접근법은 일관성이 있다는 장점이 있지만 asymptotic variable을 새로 계산해야 하는 번거로움이 있다.

다른 해법은 로그 발산을 미리 섭동계산의 경계조건에 반영하는 것이다. 구체적으로는 다음과 같이 로그 발산을 $\vec{r}_1^{(0)}$로 뽑아내고 $\vec{r}_1^{(1)}$에 대한 방정식을 푸는 것.

$$ \vec{r}_1 (t) = \vec{r}_1^{(0)} (t) + \vec{r}_1^{(1)} (t) \,,\, \vec{r}_1^{(0)} (t) = \frac{\sinh^{-1} (at)}{a^2 b^2} \hat{a} $$

로그 발산을 갖는 경계조건을 $\vec{r}_1^{(0)}$로 뽑아내었기 때문에 남는 경계조건은 $\vec{r}_1^{(1)} (-\infty) = 0$이 되며, $\vec{r}_1 (t)$는 다음과 같이 풀린다.

$$ \vec{r}_1 (t) = \vec{r}_1^{(0)} (t) + \vec{r}_1^{(1)} (t) = - \frac{ at + \sqrt{1 + a^2 t^2}}{a^2 b^2} \hat{b} + \frac{\log \left( at + \sqrt{1 + a^2 t^2} \right)}{a^2 b^2} \hat{a} $$

위 해를 $\vec{p}_2$에 대한 운동방정식에 집어넣으면 $k^2$ 차수의 운동량 변화를 구할 수 있다. 적분구간이 $(-\infty, +\infty)$로 대칭적이라는 것을 이용하면 식을 좀 다 단순화할 수 있다.

$$ \Delta \vec{p}_2 = \int_{-\infty}^{+\infty} \left[ \frac{1}{a^2 b^5 (1 + a^2 \tau^2)} - 3 \frac{\sqrt{1 + a^2 \tau^2} - a\tau \log (a\tau + \sqrt{1 + a^2\tau^2})}{a^2 b^5 (1 + a^2 \tau^2)^{5/2}} \right] \hat{b} d\tau \\ - \int_{-\infty}^{+\infty} \left[ \frac{3a^2\tau^2}{a^2 b^5 (1 + a^2 \tau^2)^{5/2}} \right] \hat{a} d\tau $$

얼핏 봐서는 적분이 꽤 복잡하게 보이는데, 의외로 적분하고 나면 값 자체는 단순하다.

$$ \Delta \vec{p}_2 = - \frac{2 \vec{a}}{a^4 b^5} = - \frac{2 \vec{p}_0}{p_0^4 b^2}$$

$k$를 전부 살린 산란 후 운동량은 다음과 같은데

$$\vec{p} (+ \infty) = \left( 1 - \frac{2 k^2}{p_0^4 b^2} \right) \vec{p}_0 - \frac{2k}{p_0 b^2} \vec{b} + \mathcal{O}(k^3)$$

제곱해보면 $k^2$ 차수에서 에너지 보존이 성립한다는 것도 확인할 수 있다.

$$ \left| {\vec{p} (+ \infty)} \right|^2 = p_0^2 + \mathcal{O}(k^3) $$

대학원 고전역학에서 다룰만한 내용으로 교수님과 이야기하다가 Dirac bracket 이야기가 나와서 간단(?)하게 트위터에서 주절거렸던 내용을 정리. 해당 타래는 이것.

이왕 Dirac bracket 이야기가 나왔으니 간단하게 아이디어 설명.

뉴턴역학/라그랑주역학은 2계미분방정식을 이용하기 때문에 초기조건으로 함수의 값과 도함수값을 요구한다. 해밀턴역학은 2계미분방정식을 1계미분방정식으로 바꿔서 초기조건에 함수의 값만 필요하도록 바꾼 것으로 이해할 수 있다.

— Dexter Kim (@AstralDexter) January 29, 2019

모든 미분방정식은 충분한 숫자의 변수를 도입하는 것으로 1계미분방정식으로 만들 수 있다. 예컨대 $y''+y=0$이란 미분방정식이 있다면 $x=y'$이란 독립변수 $x$를 도입하여 $x'+y=0$으로 만들 수 있다. 해밀턴역학도 어떤 의미에서는 그런 접근의 연장선상에 놓여있다. 르장드르 변환과도 엮여있기 때문에 좀 복잡한 방식으로 이 과정을 이용하기는 하지만.

트윗 타래에서 설명했듯, 해밀턴역학에서 해밀토니안 함수는 위상공간 위에서의 흐름(flow)을 만들어내는 물체로 생각할 수 있다. 해밀토니안 함수와 그에 대응되는 흐름 혹은 벡터장을 연결해주는 역할을 하는 것이 포아송 괄호(Poisson bracket)이다. 연결 방법은 $H \to \{H,\bullet \}$. 물론 위상공간 위에서의 흐름을 만들어내는 해밀토니안이 실제 계의 동역학과 관계가 있어야 할 이유는 없다. 보다 추상적인 임의의 함수도 포아송 괄호를 통해 위상공간 위에서 흐름을 만들어낼 수 있으며, 일반적으로는 계의 보존량 $Q$를 이용해 이런 흐름을 만들어낼 때 $Q$를 대칭 생성자(symmetry generator)라고 부른다. 이쪽은 운동량 사상(moment map)과 연결되는 방향이지만 이 글의 주제에서는 벗어나니 다음 기회에^[각주:1].

임의의 함수는 포아송 괄호를 통해 위상공간 위에서의 벡터장과 대응될 수 있다.

위의 관점은 계의 모든 변수가 독립변수인 경우에는 문제 없이 적용이 가능하지만 계의 모든 변수가 독립변수가 아닌 경우, 즉 제약조건(constraint)이 존재하는 계의 경우에는 위의 관점을 적용하는데 무리가 있다. 이 경우 좌표를 새로 잘 정의해서 새 좌표에서는 모든 변수가 독립변수가 되도록 하는 것으로 위의 관점을 살려내는 방법이 있다. 물론 새 좌표를 찾는다는 것은 원칙상 가능하다는 뜻이고, 이 좌표를 찾는 일이 항상 쉬우리란 보장은 없다. 다른 방법은 디락의 디락 괄호(Dirac bracket)를 도입하는 것.

잠시 원래 이야기에서 벗어나 역사적인 맥락을 살펴보면, 디락이 디락 괄호의 도입을 생각하게 된 이유는 양자전기역학이었다고 한다. 디락은 포아송 괄호를 교환자(commutator)로 교체하는 것으로 고전계를 양자화할 수 있다는 것을 발견했는데, 같은 방법을 전자기학에 적용하려니 뭔가 잘 안 맞는다는 것을 알게 된 것이다. 디락은 가우스 법칙에 의해 전자기장이 가질 수 있는 값에 제약이 생기는 것이 원인이라는 것을 알게 되었고, 제약조건이 있는 계의 포아송 괄호에 해당하는 물체를 어떻게 찾아낼 것인가를 고민한 결과 디락 괄호를 찾아내게 된다.

다시 원래 이야기로 돌아와서, 제약조건이 있다는 뜻은 전체 위상공간 중 그 부분집합에 해당하는 $f_i(\vec{p},\vec{q})=0$을 만족하는 $(\vec{p},\vec{q})$만 실제 계의 상태를 나타낸다는 관점으로도 이해할 수 있다. 일반적으로 해밀토니안에 의해 만들어지는 흐름은 이 제약조건을 만족하는 위상공간 속 부분다양체(submanifold) 위에서 출발하더라도 그 밖을 벗어나게 되리라고 예상할 수 있다.

해밀토니안에 의해 만들어지는 흐름(연두)은 제약조건을 만족하는 부분다양체(연파랑) 위에서 출발하더라도 그 부분다양체 위에서 움직이는 방향(녹색)과 그 부분다양체에서 벗어나는 방향(적색)을 모두 포함한다.

이제 문제는 포아송 괄호를 통해 얻은 해밀토니안 함수에 대응되는 흐름에서 제약조건을 만족하지 못하게 하는 방향의 흐름을 제거하는 것이다. 위의 그림에서 적색 화살표에 해당하는 성분을 제거하는 것이 목표인 셈. 이 목표는 제약조건을 만족하는 경우 0이란 값을 갖는 제약조건에 해당하는 함수 $f_i$들을 적당히 더하는 것으로 이루어진다. $f_i$에 의해 만들어지는 흐름 $\{f_i,\bullet\}$은 일반적으로 0이 아니기 때문. 수식으로 나타내면 다음과 같다.

\[ H \to \{ H, \bullet \}_{\text{Dirac}} = \{ H + c_i f_i , \bullet \} \]

이제 문제는 1. 충분한 숫자의 $f_i$를 찾아서 어떤 방향으로 벗어나더라도 벗어나는 방향을 제거할 수 있을 것 2. 계수들 $c_i$를 결정할 것 두가지로 나뉘게 된다. 첫번째 문제에 대한 답은 제약조건을 primary/secondary constraint와 1st class/2nd class constraint로 분류하는 과정과 관련이 있는데^[각주:2] 여기서는 일단 충분한 숫자의 $f_i$들을 구했다고 가정하기로 하자.

디락 괄호는 포아송 괄호에 보정을 가해서 제약조건을 만족시키도록 한 것으로 볼 수 있다.

계수들 $c_i$는 어떤 해밀토니안 함수를 통해 생성된 흐름이더라도 제약조건 $f_i$의 값을 0으로 유지시켜야 한다는 것으로부터 구할 수 있다. 따라서 다음 방정식의 해를 구해야 한다는 뜻이다.

\[ \forall i \,, \{ H, f_i \}_{\text{Dirac}} = 0 \]

이 문제는 다음 가설풀이(ansatz)를 적용해서 풀 수 있다. 이런 가설풀이를 도입하는 이유는 포아송 괄호의 성질들 중 필요한 성질들을 보존하기 위함인데, 그 이야기까지 하기에는 글이 너무 길어지므로 대충 넘어가기로 하자.

\[ c_i(H) = - \{ H, f_j \}M^{ji} \]

위의 가설풀이를 적용하면 이제 풀어야 할 방정식은 아래와 같이 바뀐다.

\[ \{ H, f_i \}_{\text{Dirac}} = \{ H, f_i \} - \{ H, f_k \} M^{kj} \{ f_j, f_i \} = 0\]

고맙게도 위 방정식은 단순한 역행렬 계산으로 풀 수 있다.

\[ M^{ij} \text{ is the solution to } M^{ij} \{ f_j, f_k \} = \delta^i_k \]

이 정도가 디락 괄호의 핵심적인 아이디어에 속한다.