본문으로 건너뛰기
426 / POISSON BRACKET LAB
먼저 알아두면 좋아요 · 정준 bracket·외적·직교 뉴턴 운동
단원 목록 (600) · 각운동량의 Poisson 대수와 보존
이해하는 속도로 읽으세요.개념도이론실험확인 문제
BUILD YOUR INTUITION

먼저 그림으로 이해하기

한 번에 한 단계씩 살펴보세요.

각운동량의 대수는 함수의 관계입니다. 실제 보존은 선택한 H와 토크가 결정하며 전체 공간과 한 궤적의 판정도 다릅니다.

Full 3D Cartesian canonical angular algebra {Lx,Ly}=Lz differs from time evolution {Li,H}=torque. Global component conservation depends on stiffness symmetry;selected zero-coordinate orbits can have extra constant components. Restricting z=pz=0 before differentiation changes the functions and phase space. Rotation variations are per angle,not physical velocities/forces. Three independent analytic flows;no new 3D integrator.

실제 공간의 xy 투영 -1.31101.311-1.2191.404x · my · m파랑 실제 궤적 · 주황 조회 위치t=0 sq0=(1, 1, 1) mv0=(1, 0, 1) m/s
세 축의 해석 운동을 해당 두 공간 축으로 투영합니다. 투영의 교차는 전체 위상 상태의 교차를 뜻하지 않습니다.
실제 공간의 yz 투영 -1.2190.092721.404-1.3111.311y · mz · m파랑 실제 궤적 · 주황 조회 위치t=0 sq0=(1, 1, 1) mv0=(1, 0, 1) m/s
세 축의 해석 운동을 해당 두 공간 축으로 투영합니다. 투영의 교차는 전체 위상 상태의 교차를 뜻하지 않습니다.
실제 공간의 zx 투영 -1.2190.092721.404-1.2191.404z · mx · m파랑 실제 궤적 · 주황 조회 위치t=0 sq0=(1, 1, 1) mv0=(1, 0, 1) m/s
세 축의 해석 운동을 해당 두 공간 축으로 투영합니다. 투영의 교차는 전체 위상 상태의 교차를 뜻하지 않습니다.
세 각운동량 성분 -1.201.202t · sLx, Ly, Lz · Js파랑 x · 주황 y · 초록 zL=(1, 0, -1) Jsglobal=1, 1, 1selected=1, 1, 1
전역 표시는 강성 대칭에서, 선택 궤적 표시는 영좌표 조건도 포함하여 판정합니다. 한 시점의 토크0과 다릅니다.
실제 토크와 각운동량 변화 -0.0600.0602t · staux, tauy, tauz · N m파랑 x · 주황 y · 초록 ztaux=0 N mtauy=0 N mtauz=0 N m
{Li,H}=τi이며 직접 미분한 v×p+x×F와 같습니다. 에너지 보존만으로 각운동량 보존을 판단하지 않습니다.
선택 성분의 서로의 bracket -1.1-0.50.102t · s{Lz, Lx/Ly/Lz} · Js파랑 Lx · 주황 Ly · 초록 Lzselected generator=Lzrow=(0, -1, 0) JsH와의 bracket이 아닙니다.
같은 대수 관계를 다른 H에서도 사용합니다. 서로 비가환이면서 H와는 가환하여 보존될 수 있습니다.
조회 상태의 각운동량 bracket 행렬LxLxLyLyLzLz0-101010-10행 Li · 열 Lj · 단위 Js{Lx,Ly}=Lz; {Ly,Lz}=Lx; {Lz,Lx}=LyL=(1, 0, -1) Js한 상태의 영행렬은 전역 영대수가 아닙니다.
전체 여섯 변수에서 미분한 뒤 상태를 대입합니다. 평면 제한을 먼저 대입한 함수는 다른 대상을 뜻합니다.
비등방 운동의 기계 에너지 -0.47012.3515.17102t · sH · JH=4.701 Jkx=2.467 N/mky=2.467 N/mkz=2.467 N/m
ki=m(πfi/4)²입니다. 진동수 코드를 고정하고 질량을 바꾸면 강성도 함께 바뀝니다.
Hamilton 변화율과 토크의 잔차 -5.33e-1605.33e-1602t · srate minus torque · N m파랑 x · 주황 y · 초록 zsample max=4.44e-16 N m해석식의 수치 평가 차이입니다.
토크와 각운동량의 실제 시간 미분을 별도로 대조합니다. 영잔차는 한 궤적의 보존을 뜻하지 않습니다.
각운동량 대수의 항등식 잔차 -6.00e-1906.00e-1902t · smutual bracket discrepancy · Jsselected generator=Lzsample max=0 Js동일한 단위의 성분만 비교합니다.
선택한 bracket 행을 순환 각운동량 성분과 대조합니다. 표시 최대값은 표본 진단입니다.
회전 생성자의 위치 변화 -1.201.2xyz항 또는 성분{x,Li} · mgenerator=Lz각도 매개변수 단위당 변화실제 속도가 아닙니다.
{x,Li}=ei×x입니다. 회전축 선택은 생성자 그림만 바꾸며 물리적 궤적은 바꾸지 않습니다.
회전 생성자의 운동량 변화 -0.10.51.1xyz항 또는 성분{p,Li} · kg m/sgenerator=Lz각도 매개변수 단위당 변화실제 힘이 아닙니다.
{p,Li}=ei×p입니다. 위치 변화와 운동량 변화를 단위가 다른 두 그림으로 나눕니다.
전체 3차원 각운동량 대수·토크·회전 생성자의 수치 평가
quantitynumerical valueunit
query time0s
mass1kg
initial x1m
initial y1m
initial z1m
initial vx1m/s
initial vy0m/s
initial vz1m/s
stiffness x2.467N/m
stiffness y2.467N/m
stiffness z2.467N/m
query x1m
query y1m
query z1m
query px1kg m/s
query py0kg m/s
query pz1kg m/s
Lx1Js
Ly0Js
Lz-1Js
bracket Lx,Ly-1Js
bracket Ly,Lz1Js
bracket Lz,Lx0Js
torque x0N m
torque y0N m
torque z0N m
Hamilton rate Lx0N m
Hamilton rate Ly0N m
Hamilton rate Lz0N m
energy4.701J
rotation position x-1m per unit angle
rotation position y1m per unit angle
rotation position z0m per unit angle
rotation momentum x0kg m/s per unit angle
rotation momentum y1kg m/s per unit angle
rotation momentum z0kg m/s per unit angle
rotation Jacobi example residual0m

query time s = 0

고정 기준 예제의 그림과 수치입니다. 각 단계에서 같은 그림의 의미를 차례로 살펴보세요. 고정된 기준 예제이며 아래 실험의 설정과 독립적입니다.
1 / 4

전체 여섯 변수에서 미분한다

세 정준 쌍을 유지하고 외적 함수의 기울기를 계산한 뒤 조회 상태를 넣습니다.

Li=(x×p)i

실험에서 확인하기 →
이 단원 목차

대수 관계와 보존 법칙을 나누기

두 관측량이 모두 보존되면서 서로의 Poisson bracket은 0이 아닐 수 있다. Bracket을 단순히 무엇이 변하는지 판정하는 연산으로만 읽으면 모순처럼 느껴진다. 그러나 물리적 시간 변화는 특정 해밀토니안이 생성한다. 다른 두 함수 사이의 bracket은 다른 질문에 답한다.

이를 3차원 각운동량으로 확인한다. 같은 각운동량 함수를 등방 진동자, 일부 축만 대칭인 진동자, 비등방 진동자에 적용한다. 각운동량 사이의 대수 관계는 유지되지만 해밀토니안에 따른 보존 성질은 달라진다. 실제 궤적은 세 축의 독립적인 직교 해석 흐름으로 계산하며 새로운 3차원 수치 적분기를 도입하지 않는다.

여섯 변수의 직교 위상 공간

위치를 x=(x,y,z), 운동량을 p=(px,py,pz), 양의 질량을 m으로 쓰자. 정준 쌍은 세 개다. 부호 약속은

{F,G}=∑j=13(∂F∂xj∂G∂pj−∂F∂pj∂G∂xj)\{F,G\}=\sum_{j=1}^3\left( \frac{\partial F}{\partial x_j}\frac{\partial G}{\partial p_j} -\frac{\partial F}{\partial p_j}\frac{\partial G}{\partial x_j}\right)

이다. 세 위치 성분은 모두 m, 세 운동량 성분은 모두 kg m/s다. 각운동량은 kg m²/s, 즉 Js이므로 각운동량 둘의 bracket도 Js다. 각운동량과 에너지의 bracket은 J이며, 여기서는 각도의 SI 무차원성에 따라 토크 N m로 해석한다.

단위가 같은 각운동량 세 성분을 하나의 축에서 비교할 수는 있다. 그러나 위치와 운동량을 그대로 합친 크기를 만들지는 않는다. (x,y), (y,z), (z,x) 그림은 명시한 두 공간 축으로의 투영이며, 여섯 변수의 위상 궤적과 다르다. 투영 그림의 교차만 보고 전체 상태 궤적이 자기 자신과 만난다고 결론 내리지 않는다.

각운동량을 먼저 함수로 미분하기

L=x×p를 성분으로 쓰면

Lx=ypz−zpy,Ly=zpx−xpz,Lz=xpy−ypxL_x=yp_z-zp_y,\qquad L_y=zp_x-xp_z,\qquad L_z=xp_y-yp_x

이다. Lx의 위치 기울기는 (0,pz,−py), 운동량 기울기는 (0,−z,y)다. 선택한 축의 단위 벡터를 e_i라 쓰면

∇xLi=p×ei,∇pLi=ei×x\nabla_xL_i=p\times e_i,\qquad \nabla_pL_i=e_i\times x

로 정리된다.

중요한 순서는 전체 여섯 변수에서 함수를 미분한 다음 조회 상태를 대입하는 것이다. 어떤 궤적에서 우연히 z와 pz가 0이라고 원래 함수의 그 변수에 대한 미분까지 사라지는 것은 아니다. 평면 운동의 직관과 전체 3차원 함수의 도함수를 구분한다.

한 bracket에서 순환 관계 얻기

Lx와 Ly의 기울기를 정의에 대입한다. 세 정준 쌍의 항을 더하면 xp_y−yp_x가 남는다. 즉

{Lx,Ly}=Lz,{Ly,Lz}=Lx,{Lz,Lx}=Ly\{L_x,L_y\}=L_z,\qquad \{L_y,L_z\}=L_x,\qquad \{L_z,L_x\}=L_y

이며 인자를 뒤집은 관계는 반대칭성으로 얻는다. 전체 행렬은

(0Lz−Ly−Lz0LxLy−Lx0)\begin{pmatrix} 0&L_z&-L_y\\ -L_z&0&L_x\\ L_y&-L_x&0 \end{pmatrix}

이다. 조회한 위상점의 L에 따라 행렬 값도 달라진다. L=0인 한 상태에서는 표시한 모든 성분이 0이지만 대수 자체가 전체 위상 공간에서 영인 것은 아니다.

이 유도에는 아직 H가 등장하지 않았다. 따라서 이 식만으로 어떤 각운동량 성분이 물리적 시간에 따라 변하는지는 알 수 없다. 서로의 대수와 H에 대한 변화율을 나눠 읽어야 한다.

Hamiltonian을 선택하고 토크 복원하기

세 방향 강성을 독립적인 비음수로 두자.

H=px2+py2+pz22m+kxx2+kyy2+kzz22.H=\frac{p_x^2+p_y^2+p_z^2}{2m} +\frac{k_xx^2+k_yy^2+k_zz^2}{2}.

실제 힘은 (−kx x,−ky y,−kz z)다. 이제 각운동량과 H의 bracket을 계산하면

L˙x={Lx,H}=(ky−kz)yz,\dot L_x=\{L_x,H\}=(k_y-k_z)yz, L˙y=(kz−kx)zx,L˙z=(kx−ky)xy\dot L_y=(k_z-k_x)zx,\qquad \dot L_z=(k_x-k_y)xy

를 얻는다. 다른 경로로 L=x×p를 직접 시간 미분하면 v×p+x×F다. p=mv이므로 첫 항은 0이고 남은 것은 실제 토크다. 서로의 각운동량 bracket과 별개로 뉴턴 역학에서 같은 변화율을 검산할 수 있다.

강성이 모두 양수여도 각운동량이 보존되는 것은 아니다. 비등방 복원력은 일반적으로 위치 벡터와 나란하지 않다. 에너지를 보존하면서도 토크가 생겨 각운동량 벡터의 크기나 방향을 바꿀 수 있다. 에너지 오차만 작은 모델이 올바른 회전 성질까지 재현한다고 주장해서는 안 된다.

전체 공간에서 보존되는 성분

ky=kz이면 Lx의 변화율은 모든 위상점에서 0이다. 이 H에서는 Lx가 전역적으로 보존된다. 마찬가지로 Ly에는 kz=kx, Lz에는 kx=ky가 필요하다. 세 강성이 같으면 세 각운동량 모두 보존된다.

실습에는 이보다 약한 궤적별 판정도 있다. 토크 식의 두 위치 성분 중 하나가 궤적 전체에서 항상 0이면, 강성 차이가 있어도 선택한 토크는 계속 0이다. 이 분리된 선형 방정식에서 한 좌표가 항상 0인 조건은 그 좌표의 초기 위치와 초기 속도가 모두 0인 것이다.

현재 시각에서만 토크가 0인 것은 더 약한 사실이다. 지금 한 좌표가 0이어도 곧바로 비영이 될 수 있다. 따라서 함수 전체의 보존 항등식, 선택한 궤적 위의 일정함, 조회한 한 시점의 0을 각각 표시한다. 표본 곡선이 거의 평평해 보인다는 이유로 전역 보존 표시를 켜지는 않는다.

보존되지만 서로 가환하지 않는 예

m=1 kg, 진동수 코드 (2,2,2), 초기 위치 (1,1,1) m, 초기 속도 (1,0,1) m/s를 택하자. 초기 운동량도 (1,0,1) kg m/s이므로

L=(1,0,−1) JsL=(1,0,-1)\ \mathrm{Js}

이다. 세 축의 강성이 같아 각 성분은 모두 일정하다. 그런데 {Lx,Ly}=Lz=−1 Js로 0이 아니다. 보존을 판정하는 식은 {Lx,H}와 {Ly,H}이고, 둘 사이의 bracket은 다른 관계다.

이 코드에서 각 축의 Ω는 π/2 s⁻¹이다. 같은 진동수 덕분에 궤적은 초기 위치와 속도가 만드는 물리적 평면에 놓인다. 그럼에도 함수의 대수는 전체 3차원 위상 공간에서 계산한다. 편리한 운동 평면이 있다는 사실만으로 Lx와 Ly를 독립적인 평면 정준 운동량으로 바꾸어 해석하지 않는다.

각운동량은 0인데 토크는 비영인 시작

이번에는 진동수 코드 (0,1,2), m=1 kg, 초기 위치 (1,1,1) m, 초기 속도 0을 고르자. 초기 운동량이 0이어서 L과 서로의 bracket 행렬도 모두 0이다. 하지만 세 강성은 각각 0, π²/16, π²/4 N/m이고 초기 토크는

τ=(−3π216,π24,−π216) N m\tau=\left(-\frac{3\pi^2}{16},\frac{\pi^2}{4},-\frac{\pi^2}{16}\right) \ \mathrm{N\,m}

이다. 각운동량은 즉시 변하기 시작한다. 초기 각운동량 값이나 초기 대수 행렬이 0인 것은 보존의 증거가 아니다. 이때 전체 에너지는 5π²/32 J로 유지된다.

일부 대칭을 비교하려면 코드 (4,1,1)을 사용한다. ky=kz여서 Lx는 전역 보존되지만 Ly와 Lz는 일반적으로 변한다. 단순히 “각운동량 보존”이라는 하나의 참·거짓 표시만 두면 어떤 성분이 어떤 조건에서 보존되는지 놓친다.

평면 제한을 먼저 대입하는 함정

초기 위치 (1,0,0) m와 속도 (0,1,0) m/s, 진동수 코드 (1,2,0)을 생각하자. 운동은 z=0 평면 안에 남는다. 궤적에서 Lx=Ly=0이고 초기 Lz=1 Js다. 처음에는 y=0이므로 z축 토크도 0이다. 그러나 곧 y가 변하고 kx와 ky가 다르므로 Lz는 보존되지 않는다.

전체 3차원 위상 공간에서는 {Lx,Ly}=Lz이므로 초기 상태에서 그 bracket은 1 Js다. 반면 함수에 z=pz=0을 먼저 넣어 버리면 Lx와 Ly는 작은 위상 공간의 영함수가 되어 그 평면 bracket은 0이다. 서로 다른 공간과 서로 다른 함수에 적용한 연산을 비교한 것이다.

기하적으로도 이유가 있다. Lx와 Ly가 생성하는 회전은 일반적으로 상태를 선택한 평면 밖으로 움직인다. z와 pz를 제거한 뒤에도 이 함수들이 원래의 공간 회전 생성자 역할을 그대로 유지할 수는 없다. 반면 Lz는 평면 안의 회전을 생성한다.

구속된 이론을 만들려면 올바른 축약 변수와 bracket을 함께 정해야 한다. 이번 예제는 두 연산의 불일치를 분명히 보여 주는 데 목적이 있다. 일반적인 축약 절차나 Dirac bracket 이론까지 완성했다고 주장하지 않는다. 그 조건은 이후 구속과 구조의 심화에서 다룬다.

각운동량이 회전을 생성하는 모습

축 i 주위의 작은 무차원 각도 ε에 대해

δx=ϵ{x,Li}=ϵ(ei×x),δp=ϵ{p,Li}=ϵ(ei×p)\delta x=\epsilon\{x,L_i\}=\epsilon(e_i\times x),\qquad \delta p=\epsilon\{p,L_i\}=\epsilon(e_i\times p)

다. 위치가 (1,0,0), 운동량이 (0,1,0)인 상태를 z축 주위로 돌리면 위치 변화 방향은 (0,1,0), 운동량 변화 방향은 (−1,0,0)이다. 이것은 회전 매개변수 단위당 변화이며 실제 속도나 실제 힘이 아니다. 생성자가 H인지 Li인지와 그 단위가 해석을 결정한다.

회전을 생성할 수 있는 것과 그 회전이 특정 H를 보존하는 것은 다르다. 비등방 H도 Li가 생성하는 회전을 적용할 수 있지만 일반적인 방향으로 돌린 뒤 H의 값은 바뀔 수 있다.

Jacobi 부호를 검사하는 예로 (x,Lx,Ly)를 사용한다. 순환 항은 {x,Lz}=−y, {Lx,−z}=+y, 0이며 길이 단위에서 상쇄된다. 이는 손으로 확인한 예제이지 모든 관측량에 대한 일반 증명을 대신하지 않는다.

조작 변수와 수치 결과 읽기

조작 변수는 질량, 세 축의 진동수 코드, 초기 위치 예제, 초기 속도 예제, 표시할 회전 축, 물리적 조회 시각이다. 초기 위치 선택지는 (1,0,0), (1,1,1), (1,−1,0.5), (0,0,0), (−1,0.5,1) m다. 초기 속도는 (0,1,0), (1,0,1), (−1,1,0.5), (0,0,0), (0.5,−1,1) m/s 중에서 고른다.

각 강성은 ki=m(πfi/4)²다. 진동수 코드를 유지한 채 질량을 바꾸면 강성도 바뀌어 공간 궤적은 같지만 운동량, 각운동량, 에너지, 토크는 질량에 비례한다. 질량만 바꾸고 강성을 고정한 실험과 혼동하지 않는다.

모델은 두 초의 201시점을 수치 평가한다. Bracket·토크·에너지 잔차는 작은 값을 임의로 0으로 덮지 않는다. 최대 잔차는 표본 관측값이지 연속 시간의 보증 상한이 아니다. 독립적으로 생성한 기호 기울기와 60자리 궤적 미분을 대조하며, 위치·운동량 단위가 섞인 잔차 크기는 만들지 않는다. 공간 투영은 표시한 두 축의 단위를 유지하고 물리적 궤적과 위상 공간을 구분한다.

해설이 있는 연습 문제

  1. 좌표 기울기에서 {Lx,Ly}를 유도하자. 해설: 상태를 대입하기 전에 세 정준 쌍을 전개하면 xp_y−yp_x=Lz가 남는다. 인자를 바꾸면 −Lz다.
  2. 서로의 bracket이 비영인데 Lx와 Ly가 모두 보존될 수 있을까? 해설: 가능하다. 등방 H에서 {Lx,H}={Ly,H}=0이고 서로의 bracket은 비영일 수 있는 Lz다. 각각 다른 질문에 답한다.
  3. kx가 다르고 ky=kz이면 무엇이 전역 보존될까? 해설: 토크 계수가 ky−kz인 Lx다. 다른 성분도 특별한 선택 궤적에서는 일정할 수 있지만 이는 전역 판정과 구분한다.
  4. 초기 L=0인데 초기 변화율은 왜 0이 아닐 수 있을까? 해설: 코드 (0,1,2)의 정지 시작에서는 운동량이 0이지만 위치와 비등방 힘의 외적이 비영이다. 초기 토크가 즉각적인 각운동량 변화를 예측한다.
  5. 평면에서 두 bracket 계산이 달라지는 이유를 설명하자. 해설: 미분 전에 z=pz=0으로 제한하면 작은 위상 공간의 다른 함수가 된다. 이를 나중에 상태만 대입한 전체 여섯 변수의 대수와 동일시할 수 없다.
  6. {x,Lz}는 실제 속도인가? 해설: 회전 매개변수에 따른 위치 변화다. 실제 속도는 {x,H}다. 생성자의 의미와 단위를 먼저 확인해야 한다.

AI 연구와 읽을거리

학습한 동역학에서는 표현의 대수, 해밀토니안의 대칭, 선택한 데이터 궤적의 보존을 분리한다. 등방 궤적만 보고 세운 보존 가설이 비등방 계에서도 성립한다고 할 수 없다. 강성과 초기 조건을 개입하여 바꾸고 각운동량의 변화를 실제 토크와 대조한다. 에너지 보존 하나만 맞는지 보는 검증보다 구체적인 실패 원인을 구분할 수 있다.

각운동량 bracket과 생성자 관점의 읽을거리는 David Tong의 Hamiltonian Formalism 4.3.1·4.4.1절이다. 이 단원의 비등방 실험, 전체·제한 위상 공간 비교, 수치 예제는 직접 유도하고 독립 기준과 검산했다.

CHECK YOUR UNDERSTANDING

설명할 수 있나요?

실험 설정과 무관한 기준 예제를 묻습니다. 답을 고르고 설명을 확인해보세요.

1. Lx와 Ly가 모두 보존되면 서로의 bracket도 0인가요?
2. 평면 궤적에서 z=pz=0을 먼저 넣으면 3차원 대수가 유지되나요?

선택한 답과 학습 표시는 이 브라우저에 저장됩니다. 다시 풀거나 표시를 해제할 수 있습니다.