원 구속과 Dirac 위상 기하

위치 구속은 운동량도 제한한다

원 위에 구속된 입자는 임의의 법선 운동량을 가질 수 없습니다. 위치 방정식만 맞춘다고 일관된 위상 공간 상태가 되는 것은 아닙니다. 위치 관계를 시간 내내 유지하려면 운동량도 접선 방향이어야 합니다.

이 단원에서는 두 구속에 대한 Dirac bracket을 유도합니다. 원래의 정준 bracket과 계수2의 구속 구조를 비교하며 세 질문을 나눕니다. 여기서 bracket은 정규적인가요? 현재 상태는 원하는 구속면 위에 있나요? 선택한 Hamiltonian은 어떤 관측량을 보존하나요? 의도적으로 잘못된 입력에서는 답이 다르게 나옵니다.

모델은 정확한 순간 벡터장과 구조 항등식을 계산합니다. 수치 궤적을 적분하지 않으며, 기준 원을 그렸다고 운동을 시뮬레이션한 것도 아닙니다.

위상 공간·단위·구속 정하기

정준 좌표 순서는 $z=(q_x,q_y,p_x,p_y)$이고 q는 m, p는 kg m/s 단위입니다. 이 순서의 정준 행렬은

$$J=\begin{pmatrix}0&I\\-I&0\end{pmatrix},\qquad \{f,g\}=\nabla f^T J\nabla g$$

입니다. 양의 목표 반지름 R에 대해

$$\phi_1=\frac{q\cdot q-R^2}{2},\qquad \phi_2=q\cdot p$$

를 정의합니다. 첫째 단위는 m², 둘째는 Js입니다. 동일한 수치 허용오차로 합쳐 평가할 측정값 둘이 아니라 서로 다른 조건입니다. 원하는 물리적 위상 공간은 두 값이 동시에0인 집합입니다. 정준 자유 속도 $\dot q=p/m$로 첫 구속을 미분하면 $\dot\phi_1=\phi_2/m$이므로 운동량 조건이 필요합니다.

질량은 정수 kg 설정을 사용합니다. 반지름·직교 위치·운동량·외력 코드는 각각 m, kg m/s, N 단위의 입력/2입니다. 목표 반지름은 양수로 유지합니다. 위치0은 물리적 원 밖에 있는 특이 상태를 명시적으로 검사하기 위해 선택할 수 있습니다.

구속 행렬을 만들고 역행렬 구하기

구속을 대입하기 전에 bracket을 미분합니다. $s=q\cdot q$라 두면

$$C_{ab}=\{\phi_a,\phi_b\} =\begin{pmatrix}0&s\\-s&0\end{pmatrix},\qquad C^{-1}=\begin{pmatrix}0&-1/s\\1/s&0\end{pmatrix}$$

입니다. q가 비영이면 역행렬이 존재합니다. 행렬식은 m⁴ 단위의 $s^2$입니다. 목표 원에서는 $s=R^2$이므로 두 구속은 제2종 구속 쌍입니다. 역행렬이 존재하는 이 구속 bracket 행렬과, 뒤에서 구하는 계수2의 주변 Dirac 행렬은 서로 다른 대상입니다.

0인 칸에도 행과 열의 물리량에서 정해지는 단위가 있습니다. 예를 들어 C11은 m² s/kg, C22는 Js, 비영인 비대각 칸은 m²입니다. 역행렬의 비대각 칸은1/m²입니다. 모든 칸에 하나의 크기 색상 눈금을 적용하면 서로 다른 단위의 수치를 섞게 됩니다.

역행렬로 보정 bracket 정의하기

Dirac bracket은

$$\{f,g\}_D=\{f,g\}-\{f,\phi_a\}(C^{-1})^{ab}\{\phi_b,g\}$$

입니다. G를 두 구속의 기울기를 담은2×4 행렬이라 하면 코드는

$$D=J-JG^T C^{-1}GJ$$

로 전체 행렬을 직접 계산합니다. 그림을 보고 개별 성분의 부호를 추측하지 않습니다. 곱하면 $GD=0$이므로 정규 영역에서 두 구속은 모든 매끄러운 관측량과의 Dirac bracket이0입니다. 이 퇴화한 Poisson 구조의 Casimir라는 뜻입니다.

반대칭이라는 사실만으로 Poisson bracket이 되지는 않습니다. 독립 생성기는 수치를 대입하기 전에 일반 상태의 행렬을 미분하여 독립적인 좌표 Jacobi 삼중항 네 개를 기호적으로 검사합니다. 유한한 표본 몇 개가 맞는다는 결과를 이 기호 검사의 대체물로 제시하지 않습니다.

블록을 유도하고 운동량 부호 검사하기

위치·운동량 블록은 접선 투영자

$$P=I-\frac{qq^T}{s}$$

입니다. 기본 bracket들은

$$\{q_i,q_j\}_D=0,\qquad \{q_i,p_j\}_D=P_{ij},\qquad \{p_i,p_j\}_D=\frac{p_iq_j-q_ip_j}{s}$$

입니다. q-q 블록은 s/kg, q-p와 p-q 블록은 무차원, p-p 블록은 kg/s 단위입니다. 이는4차원 Euclidean 계량이 아닙니다. 행렬 계수는 의미가 있지만 단위가 다른 블록 크기를 더한 노름을 물리적 오차로 해석하면 안 됩니다.

$q=(R,0)$, $p=(0,P_t)$에서 부호를 검산하세요. $\{p_x,p_y\}_D=-P_t/R$입니다. 자유 운동에서는 $\dot p_x=-P_t^2/(mR)$로 안쪽 운동량 변화가 나옵니다. 반대 부호라면 바깥으로 밀어 구속의 미분 조건을 위반합니다. 정의·좌표 순서·이 검산이 사용할 부호 규약을 정합니다.

서로 다른 두 투영자 구분하기

2차원 P는 $P^2=P$이고 계수1이며 배치 벡터를 접선으로 투영합니다. 위상 공간 투영자는

$$\Pi=I-JG^TC^{-1}G$$

입니다. $\Pi^2=\Pi$, $G\Pi=0$, $D=\Pi J\Pi^T$를 만족하며 계수는2입니다. 위치와 운동량 단위가 다른 네 위상 성분에 Euclidean 내적을 임의로 넣었을 때의 직교 투영자라고 일반적으로 부를 수는 없습니다.

Dirac 행렬 D도 정규 주변 영역의4×4 행렬로는 계수2입니다. 정규2차원 구속 잎에 제한하면 그곳의 구조는 비퇴화합니다. 주변 공간에서 특이한 행렬이어도 적절한 축약 심플렉틱 위상 공간과 양립합니다. D에 무조건 계수4의 역행렬을 요구하지 말고 제한한 면의 차원을 확인하세요.

구속 Hamiltonian 벡터장 유도하기

일정한 직교 외력 F에 대해

$$H=\frac{p\cdot p}{2m}-F\cdot q$$

를 사용합니다. 보정된 벡터장은

$$\dot q=\frac{Pp}{m},\qquad \dot p=PF+\frac{p(q\cdot p)-q(p\cdot p)}{ms}$$

입니다. 두 구속 변화율은 정확히0입니다. 물리적인 영구속면에서는 $q\cdot p=0$이므로 속도가 p/m이고

$$\dot p=PF-\frac{p\cdot p}{mR^2}q$$

입니다. 접선 외력과 안쪽 구심 운동량 변화가 함께 들어갑니다. 이상적 반력은 $\dot p-F$이고 물리적 속도가 접선이므로 그 일률은0입니다. 코드는 두 구속이 모두 만족될 때만 물리적 반력을 제공합니다. 면 밖에서도 수학적 벡터장은 정의되지만 p가 질량×속도인 의도한 입자의 뉴턴 방정식으로 무작정 해석할 수 없습니다.

올바른 상태에서 정준 자유 운동과 비교하기

질량1 kg, 목표 반지름1 m, 위치 $q=(1,0)$ m, 운동량 $p=(0,1)$ kg m/s, 외력0을 사용합니다. 허용되는 위상 상태입니다.

정준 자유 벡터장은 $\dot q=(0,1)$ m/s, $\dot p=(0,0)$ N을 줍니다. 첫 구속의 순간 변화율은0이지만 둘째 변화율은1 J입니다. 이어지는 변화가 원과 양립하는 운동량 조건을 벗어나므로 초기의 접선 조건만으로는 충분하지 않습니다.

Dirac 벡터장은 같은 초기 속도와 $\dot p=(-1,0)$ N을 줍니다. 두 구속 변화율은0입니다. 행렬 계수는2이고 $\{q_x,p_x\}_D=0$, $\{q_y,p_y\}_D=1$, $\{p_x,p_y\}_D=-1$ kg/s입니다. 반력은 안쪽이며 일률은0입니다.

이 결과는 정확한 연속 벡터장에 관한 것입니다. 일반적인 유한 Euler 걸음은 입력 속도가 접선이어도 원을 벗어날 수 있습니다. 적분기는 별도의 구속·수렴 분석이 필요하며 연속 bracket 항등식이 수치 궤적을 인증하지는 않습니다.

보존과 초기 상태 수정 혼동하지 않기

위치는 $(1,0)$으로 두고 운동량을 $(1,1)$로 바꿉니다. 위치는 목표 원 위에 있지만 $\phi_2=1$ Js이므로 위상 상태는 잘못되었습니다. 그래도 bracket은 정규적이고 두 구속 변화율은0입니다.

Dirac 속도는 $(0,1)$ m/s이며 $p/m=(1,1)$ m/s와 다릅니다. 운동량 변화율은 $(-1,1)$ N입니다. 흐름은 현재의 비영인 phi 두 값을 유지할 뿐,0으로 초기화하지 않습니다. 코드가 입력 운동량을 접선 투영으로 몰래 바꾼 것도 아닙니다.

위치 오차도 같습니다. 목표 반지름1에서 $q=(1/2,0)$이면 phi 하나는−3/8 m²입니다. 정규 벡터장은 이 잘못된 수준을 보존합니다. Casimir 항등식은 현재 잎 안의 운동을 설명합니다. 원하는 영수준 잎에 도달하려면 올바른 초기 상태, 일관된 초기화 알고리즘 또는 별도로 명시한 투영 절차가 필요합니다.

에너지와 각운동량 검사 나누기

반대칭성 때문에 모든 정규 잎에서 $\dot H=\nabla H^TD\nabla H=0$입니다. 의도한 물리적 원이 아닌 잎에서도 성립합니다. 따라서 에너지 보존 하나로 초기 상태의 적합성을 인증할 수 없습니다.

각운동량 $L=q_xp_y-q_yp_x$에 대해서는

$$\dot L=q_xF_y-q_yF_x$$

입니다. 이 외력족에서는 외력0일 때 전역 회전 보존이 성립합니다. 일정한 외력이 비영이면 일반적으로 깨지지만 조회 상태의 토크는0일 수 있습니다. $q=(1,0)$, $p=(0,1)$, $F=(0,1)$ N에서는 각운동량 변화율이1 N m이면서 에너지는 보존됩니다. 구속 반력의 일률도 여전히0입니다.

올바른 구속, 정규 bracket, H의 보존과 L의 보존은 서로 다른 검사입니다. 외력 설정은 Hamiltonian을 바꾸며 bracket의 구성은 바꾸지 않습니다. 어떤 Hamiltonian과 영역에 대한 보존인지 명시하세요.

각도의 분기 절단과 bracket 특이성 구분하기

국소 각도 미분의 배치 기울기는 $(-q_y/s,q_x/s)$입니다. 각운동량과 결합하면 국소 정준 관계 $\{\theta,L\}_D=1$을 얻습니다.

하지만 원 전체에서 하나의 단일값 각도가 생기는 것은 아닙니다. 표시한 주값 atan2는 음의 x축에 절단이 있습니다. 그 축의 비영 점에서도 C는 가역이고 다른 국소 각도 좌표에서는 bracket1입니다. 선택한 주값 좌표의 미분만 사용할 수 없으므로 그 bracket 표시는 미정으로 남깁니다.

물리적 영구속 잎의 국소 정준 쌍에서는 $H_{\mathrm{red}}=L^2/(2mR^2)-R(F_x\cos\theta+F_y\sin\theta)$입니다. $\dot\theta=L/(mR^2)$와 주변 Dirac 벡터장과 같은 외부 토크를 줍니다. 코드는 허용되는 상태에서만 이 축약 에너지와 주변 H를 비교합니다. 국소 정준 쌍이 있다고 위상 공간의 원기둥 전체가 하나의 좌표 평면과 전역적으로 같아지지는 않습니다.

원점에서는 문제가 더 큽니다. C 자체가 특이하므로 각도 미분과 Dirac 역행렬이 모두 정의되지 않습니다. 정규 원 위에서 다른 좌표로 해결할 절단과 q0의 실제 역행렬 실패를 구분하세요. 작은 분모를 더하면 해결되는 것이 아니라 두 문제를 가립니다.

특이 상태와 기준 그림을 정직하게 해석하기

목표 반지름 설정은 항상 양수이므로 원점은 물리적 면 위에 올 수 없습니다. q0에서는 Dirac 기하·역행렬·구속 벡터장·축약 잎 차원을 제공하지 않습니다. 반면 주변 정준 J는 그곳에서도 정의되므로 정준 벡터장은 표시할 수 있습니다.

원점의 구속 기울기 계수는 운동량이 비영이면1, 운동량도0이면0일 수 있습니다. 구속 식이 두 개라고 적는 것만으로 정규성을 확인한 것은 아닙니다.

기준 원은 목표 위치 집합을 보여주기 위한 삼각함수 표본 그림입니다. 정확 유리수 궤적이나 구속 보존의 증명이 아닙니다. 구조의 계수·영점·힘·bracket은 정확 유리수로 계산하고 주값 각도와 원 그림은 수치 기하 보조로 구분합니다.

학습한 구속 모델에 구조 검사 적용하기

학습한 bracket에는 보기 좋은 궤적 하나보다 강한 검사가 필요합니다. 반대칭·구속 소거·정규 영역의 계수·Jacobi 항등식을 검사하세요. 그 뒤 선택한 Hamiltonian 벡터장, 초기 상태의 적합성, 국소 각도와 주변 좌표의 차이를 확인해야 합니다.

반대로 에너지만 일정하게 만드는 신경망은 잘못된 구속 잎을 유지할 수 있습니다. 한 표본의 접선 속도만 보면 둘째 운동량 구속을 놓칠 수 있습니다. 단위가 다른 행렬 블록을 하나의 노름으로 합치면 어느 항등식이 실패했는지 숨길 수 있습니다.

이 예제는 제2종 구속 쌍이며 국소 게이지 중복성이 아닙니다. 제1종 구속, 게이지 고정, 양자화와 연산자 순서 문제에는 추가 구성이 필요합니다. 원의 bracket은 그 주제로 들어가는 출발점이지 전체 분석을 구현했다는 주장이 아닙니다.

해설이 있는 연습문제

1. 빠진 조건. 위 물리적 자유 상태에서 정준 벡터장이 주는 두 구속 변화율을 비교하세요.

   해설: $\dot\phi_1=q\cdot p/m=0$이지만 $\dot\phi_2=p^2/m=1$ J입니다. 순간 위치 접선 조건이 운동량 조건까지 유지하지는 않습니다. Dirac 벡터장은 두 변화율을 모두0으로 만듭니다.

2. 운동량 bracket 부호. $q=(R,0)$, $p=(0,P_t)$에서 p-p bracket으로 법선 운동량 변화율을 구하세요.

   해설: $\{p_x,p_y\}_D=-P_t/R$입니다. $\partial H/\partial p_y=P_t/m$을 곱하면 안쪽인 $-P_t^2/(mR)$ N입니다. 반대 부호는 원의 일관성 조건을 위반합니다.

3. 잘못된 운동량 잎. 목표 반지름1, $q=(1,0)$, $p=(1,1)$에서 Dirac 구속 변화율0은 초기 상태를 수정했다는 증거인가요?

   해설: 아닙니다. phi 두는1 Js이며 비영인 그 값에 남습니다. 속도는 $(0,1)$ m/s로 p/m과 다릅니다. 정규적이지만 잘못된 잎에 있고 의도한 영구속면의 물리적 반력은 제공하지 않습니다.

4. 외부 토크와 에너지. 올바른 단위 상태와 $F=(0,1)$ N에서 각운동량 변화율과 반력 일률을 구하세요.

   해설: 토크는 $q_xF_y-q_yF_x=1$ N m입니다. 반력은 $(-1,0)$ N, 속도는 $(0,1)$ m/s이므로 일률은0입니다. L은 변해도 전체 H는 보존됩니다.

5. 좌표 경계. $q=(-1,0)$에서 주값 각도가 뛰면 Dirac bracket도 특이한가요?

   해설: 아닙니다. $q^2=1$이므로 C는 가역이고 축약 잎은2차원입니다. 다른 국소 좌표로 $\{\theta,L\}_D=1$을 얻습니다. 선택한 주값 분기의 미분만 절단에서 사용할 수 없습니다.

6. 실제 특이성. q0이고 p가 비영일 때 분모에 목표 반지름을 대신 넣어 bracket을 정의해도 되나요?

   해설: 아닙니다. C는 실제 q로 정해지며 행렬식이0입니다. s를 목표 $R^2$로 바꾼 행렬은 C의 역행렬이 아니며 일반적으로 Casimir 항등식을 잃습니다. 원점은 양의 반지름 구속면 밖에 있으므로 초기화를 수정해야 합니다.

읽을거리와 부호 규약 검산

Brian P. Dolan의 Constrained Dynamics in the Hamiltonian formalism 정의(12)와3·5절은 구속 역행렬 구성을 설명합니다. Hagen Kleinert와 Sergei V. Shabanov의 Proper Dirac Quantization of Free Particle on D-Dimensional Sphere 식(5)–(9)는 bracket 부호를 교차 확인하는 자료이며 여기서는 고전 구조만 사용합니다. 성분 부호는 이 단원에 명시한 정의에서 유도하고 Casimir·Jacobi·안쪽 힘으로 검산했습니다. 그 논문의 양자적 결론을 고전 예제로 구현하거나 추론하지 않습니다.