Dexter's story

대학원 고전역학에서 다룰만한 내용으로 교수님과 이야기하다가 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 \]

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