고전 및 라그랑주 역학
Advanced Mechanics
2단계: 변분법과 라그랑주·해밀턴 역학
1부: 이론적 기초 — "왜 뉴턴을 버려야 하는가?"
1단계에서 우리는 F = ma, 각운동량, 에너지 보존, 강체 역학까지 뉴턴 역학의 전반을 다뤘다. 그런데 솔직하게 말해보자. 뉴턴의 방법은 어떤 경우에는 지독하게 불편하다. 이중 진자(double pendulum)를 생각해보자. 두 개의 막대가 끝에 연결된 단순해 보이는 장치인데, 뉴턴식으로 풀려면 두 질점에 작용하는 장력, 중력을 각각 벡터로 분해하고, 제약 조건(constraint)을 처리하기 위한 연립 미분방정식을 만들어야 한다. 미지수가 넘쳐나고 수식은 복잡해진다. 마찰이나 경사면이 추가되면 더 나빠진다.
여기서 근본적인 질문이 생긴다. 왜 힘(force)이라는 벡터 개념을 꼭 써야 하는가? 힘은 방향이 있어서 좌표계를 어떻게 잡느냐에 따라 표현이 달라진다. 그러나 에너지(energy)는 스칼라(scalar)다. 방향이 없다. 좌표계가 바뀌어도 운동 에너지 T나 퍼텐셜 에너지 V의 값은 변하지 않는다. 18세기의 수학자·물리학자들, 특히 조제프-루이 라그랑주(Joseph-Louis Lagrange)는 이 아이디어를 밀어붙였다. "힘을 쓰지 말고, 에너지만으로 운동을 기술할 수 있지 않을까?" 그 결과가 **라그랑주 역학(Lagrangian Mechanics)**이다.
[노트 기록] 뉴턴 역학 vs. 라그랑주 역학의 핵심 차이: 뉴턴은 벡터(힘) 중심, 라그랑주는 스칼라(에너지) 중심. 뉴턴은 각 질점의 힘을 분석하고, 라그랑주는 시스템 전체의 에너지를 분석한다.
그런데 에너지만으로 운동을 기술하려면, 먼저 좌표 체계를 새롭게 정의해야 한다. 이것이 **일반화 좌표(Generalized Coordinates)**의 개념이다. 단진자를 예로 들면, 줄에 매달린 추의 위치를 x, y 좌표 두 개로 표현할 수도 있다. 하지만 줄의 길이가 고정되어 있으므로 x² + y² = L²이라는 제약 조건이 생기고, 실제 자유도(degree of freedom)는 1개뿐이다. 이때 그냥 각도 θ 하나만 써도 추의 위치를 완벽하게 기술할 수 있다. 이 θ가 일반화 좌표다. 직교 좌표일 필요도 없고, 차원이 같을 필요도 없다. 시스템의 자유도를 가장 자연스럽게 포착하는 어떤 수(또는 수들의 집합)도 일반화 좌표가 될 수 있다. 이중 진자라면 두 각도 θ₁, θ₂가 일반화 좌표가 되어 복잡한 벡터 분석 없이 시스템을 기술한다.
[노트 기록] 일반화 좌표 qᵢ는 시스템의 자유도(DOF, Degree of Freedom) 수만큼 필요하다. DOF = (전체 좌표 수) - (제약 조건 수). 이중 진자: 4개 좌표(x₁,y₁,x₂,y₂) - 2개 제약(각 막대 길이 고정) = 2 DOF → q₁=θ₁, q₂=θ₂.
이제 더 근본적인 질문을 던져볼 시간이다. 공이 A 지점에서 B 지점으로 날아갈 때, 왜 공은 그 특정한 경로를 선택하는가? 다른 경로로 가면 안 되는가? 이 질문에 답하는 것이 **최소 작용 원리(Principle of Least Action)**이고, 이것을 수학적으로 다루는 도구가 **변분법(Calculus of Variations)**이다. 미적분학이 함수의 극값(최솟값, 최댓값)을 찾는 것이라면, 변분법은 범함수(Functional)의 극값을 찾는 것이다. 범함수란 "함수를 입력받아 숫자를 출력하는 것"이다. 예를 들어 "두 점 사이의 모든 가능한 경로들 중에서 어느 것의 길이가 가장 짧은가?"라는 질문에서, 길이(숫자)를 경로(함수)로부터 계산하는 것이 바로 범함수다.
변분법을 직관적으로 이해하기 위해 빛의 굴절을 생각해보자. 빛이 공기에서 물로 들어갈 때 왜 꺾이는가? 페르마(Fermat)는 "빛은 가장 빠른 경로를 선택한다"라고 했다(페르마의 원리). 공기보다 물에서 빛이 느리므로, 직선으로 가는 것보다 공기 중에서는 더 많이 이동하고 물속에서는 덜 이동하는 경로가 시간적으로 더 빠를 수 있다. 빛은 마치 "모든 가능한 경로를 계산해보고" 그 중 시간이 최소인 경로를 선택하는 것처럼 행동한다. 역학에서 최소 작용 원리도 이와 완벽하게 유사한 구조를 가진다.
2부: 본 내용 — 변분법에서 해밀턴 역학까지
2-1. 변분법과 오일러-라그랑주 방정식
**작용(Action)**이란 무엇인가? 라그랑지안(Lagrangian) L을 T - V, 즉 운동 에너지에서 퍼텐셜 에너지를 뺀 값으로 정의하자. (왜 T + V가 아니라 T - V인지는 나중에 명백해진다.) 시간 t₁에서 t₂까지의 **작용 S(Action)**는 다음과 같다.
S = ∫[t₁ to t₂] L(q, q̇, t) dt
이 S는 "시스템이 어떤 경로를 따라 움직였느냐"에 따라 다른 숫자를 내놓는다. 경로가 다르면 S도 다르다. 그래서 S는 경로라는 함수를 받아 숫자를 돌려주는 범함수다. 이제 해밀턴의 원리(Hamilton's Principle)가 등장한다.
[노트 기록] 해밀턴의 원리(최소 작용 원리): 실제 자연이 선택하는 경로는, 고정된 양 끝점(초기 상태와 최종 상태)을 연결하는 모든 가능한 경로들 중에서 작용 S를 극값(stationary value)으로 만드는 경로다. 수학적으로: δS = 0.
"δS = 0"에서 δ는 변분(variation)이라는 기호다. 미분이 "함수의 작은 변화에 대한 출력 변화"를 다룬다면, 변분은 "경로 자체의 작은 변화에 대한 S의 변화"를 다룬다. q(t)라는 실제 경로에서 아주 조금 벗어난 경로 q(t) + εη(t)를 생각하자. 여기서 η(t)는 임의의 변화 함수이고, 양 끝점에서는 η(t₁) = η(t₂) = 0이다(시작점과 끝점은 고정). S의 변분이 0이 되어야 한다는 조건을 수학적으로 처리하면, 모든 일반화 좌표 qᵢ에 대해 다음 방정식이 나온다.
오일러-라그랑주 방정식: d/dt (∂L/∂q̇ᵢ) - ∂L/∂qᵢ = 0
이것이 라그랑주 역학의 핵심이다. 이 방정식을 한번 음미해보자. ∂L/∂q̇ᵢ는 L을 일반화 속도로 편미분한 것인데, 이것을 일반화 운동량(generalized momentum) pᵢ라고 부른다. ∂L/∂qᵢ는 L을 일반화 좌표로 편미분한 것이고, 이것은 **일반화 힘(generalized force)**에 해당한다. 따라서 오일러-라그랑주 방정식은 사실 "일반화 운동량의 시간 변화율 = 일반화 힘"이라는 뜻이고, 이것은 뉴턴의 F = ma의 에너지 언어로의 번역이다.
단진자로 이것을 직접 확인해보자. 길이 L인 단진자에서 일반화 좌표 q = θ(각도)로 놓으면, 운동 에너지는 T = ½m(Lθ̇)²이고 퍼텐셜 에너지는 V = -mgL cosθ(최저점을 기준 0으로 잡을 수도 있지만, 지지점 기준으로 쓰면 -mgL cosθ)다. 따라서 L = T - V = ½mL²θ̇² + mgL cosθ. 이제 오일러-라그랑주 방정식을 적용하면: ∂L/∂θ̇ = mL²θ̇이므로 d/dt(mL²θ̇) = mL²θ̈. 그리고 ∂L/∂θ = -mgL sinθ. 따라서 mL²θ̈ - (-mgL sinθ) = 0, 즉 mL²θ̈ + mgL sinθ = 0, 곧 θ̈ + (g/L) sinθ = 0이 나온다. 이것은 우리가 뉴턴 방식으로 장력과 중력을 분해해서 얻는 바로 그 방정식이다. 그러나 라그랑주 방식에서는 장력이라는 제약력(constraint force)을 단 한 번도 고려하지 않았다! 이것이 라그랑주 역학의 강력함이다.
[노트 기록] 오일러-라그랑주 방정식 유도 핵심 단계: L 정의 → 각 일반화 좌표 qᵢ에 대해 ∂L/∂q̇ᵢ를 시간 미분, ∂L/∂qᵢ를 계산 → 두 값의 차가 0. 제약력은 신경 쓸 필요 없다.
2-2. 해밀턴 역학 체계와 위상 공간
라그랑주 역학은 일반화 좌표 qᵢ와 일반화 속도 q̇ᵢ로 시스템을 기술한다. 해밀턴(William Rowan Hamilton)은 이것을 더 우아하게 바꿨다. "속도 q̇ᵢ 대신 일반화 운동량 pᵢ = ∂L/∂q̇ᵢ를 독립 변수로 쓰면 어떨까?" 이 전환은 **르장드르 변환(Legendre Transform)**을 통해 이루어지고, 그 결과물이 해밀토니안(Hamiltonian) H다.
H(q, p, t) = Σ pᵢq̇ᵢ - L(q, q̇, t)
많은 경우(특히 퍼텐셜이 속도에 의존하지 않고 좌표 변환이 시간에 무관할 때), H는 정확히 전체 역학적 에너지 T + V와 같다. 라그랑지안 L = T - V와 비교해보라. L은 에너지의 "차이"를 다루고, H는 에너지의 "합"을 다룬다. 이 두 관점이 쌍을 이루는 것이 아름답다.
해밀토니안을 사용하면 운동 방정식이 **해밀턴 방정식(Hamilton's Equations)**으로 바뀐다.
q̇ᵢ = ∂H/∂pᵢ ṗᵢ = -∂H/∂qᵢ
[노트 기록] 라그랑주 방정식은 qᵢ에 대한 2차 미분방정식(가속도 θ̈가 나오므로) n개다. 해밀턴 방정식은 qᵢ, pᵢ에 대한 1차 미분방정식 2n개다. 수식이 두 배로 늘어났지만, 각각이 더 단순하고 대칭적이다. 이 대칭성이 나중에 엄청난 수학적 도구들을 불러온다.
이제 **위상 공간(Phase Space)**을 이해할 시간이다. 1단계에서 2D 강체 물리 엔진을 만들 때, 각 물체의 상태를 "위치 + 속도"로 표현했을 것이다. 위상 공간은 이 아이디어의 정확한 형식화다. n개의 자유도를 가진 시스템에서, (q₁, ..., qₙ, p₁, ..., pₙ)으로 이루어진 2n차원 공간을 위상 공간이라 한다. 시스템의 순간 상태는 이 공간의 한 점으로 완벽하게 기술되고, 시간이 지나면 이 점은 위상 공간에서 하나의 **궤적(trajectory)**을 그린다.
단진자를 다시 예로 들자. q = θ, p = mL²θ̇인 2차원 위상 공간을 생각해보자. 에너지가 낮으면(작은 진폭 진동) 궤적은 닫힌 타원처럼 그려진다. 에너지가 충분히 높으면(추가 계속 회전) 궤적은 끝없이 이어지는 파도 모양이 된다. 그리고 그 경계, 즉 정확히 추가 꼭대기(θ = π)에 극히 불안정하게 서 있는 경우에는 궤적이 **분리선(separatrix)**이라는 특별한 곡선을 그린다. 에너지가 같은 모든 상태는 위상 공간에서 같은 곡선 위에 있다. 이것은 해밀토니안 H = 상수가 위상 공간에서 등에너지 곡면(hypersurface)을 정의하기 때문이다.
[노트 기록] 류빌 정리(Liouville's Theorem): 위상 공간에서 초기 상태들의 집합이 이루는 부피는 시간이 지나도 변하지 않는다. (증명은 해밀턴 방정식의 발산(divergence)이 0임을 보이면 된다.) 이것은 해밀턴 역학 시스템이 "비압축성(incompressible)"임을 의미한다. 마치 물이 관 속을 흐를 때 부피가 보존되는 것처럼.
2-3. 정준 변환과 포아송 괄호
위상 공간에서 좌표 변환을 해도 해밀턴 방정식의 형태가 유지되는 변환을 **정준 변환(Canonical Transformation)**이라 한다. 왜 이것이 중요한가? 어떤 문제를 풀 때, 좌표계를 영리하게 바꾸면 해밀토니안이 엄청나게 단순해질 수 있기 때문이다. 마치 경사면 문제를 풀 때 경사를 따라 x축을 잡으면 계산이 쉬워지는 것처럼.
정준 변환의 핵심을 이해하는 도구가 **포아송 괄호(Poisson Bracket)**다. 두 관측량 f, g에 대해:
{f, g} = Σᵢ (∂f/∂qᵢ · ∂g/∂pᵢ - ∂f/∂pᵢ · ∂g/∂qᵢ)
포아송 괄호로 해밀턴 방정식을 다시 쓰면: q̇ᵢ = {qᵢ, H}, ṗᵢ = {pᵢ, H}. 그리고 임의의 관측량 f의 시간 변화율은 ḟ = {f, H} + ∂f/∂t가 된다. f가 명시적으로 시간에 의존하지 않는다면, {f, H} = 0일 때 f는 보존된다! 이것이 보존량의 기준이다. 정준 변환은 포아송 괄호를 보존하는 변환으로 정의된다. {qᵢ, pⱼ}의 기본 포아송 괄호 구조가 변환 후에도 유지되면, 그 변환은 정준 변환이다.
[노트 기록] 기본 포아송 괄호: {qᵢ, qⱼ} = 0, {pᵢ, pⱼ} = 0, {qᵢ, pⱼ} = δᵢⱼ (δᵢⱼ는 i=j일 때 1, 아닐 때 0인 크로네커 델타). 이 구조가 보존되어야 정준 변환이다.
**해밀턴-야코비 이론(Hamilton-Jacobi Theory)**은 정준 변환의 궁극적 응용이다. 이 이론에서는 해밀토니안 자체를 0으로 만드는 정준 변환을 찾는데, 그러면 운동 방정식이 trivial하게 풀린다. 이 이론은 양자역학의 슈뢰딩거 방정식과 깊이 연결되어 있어서, 고전 역학이 양자 역학으로 넘어가는 다리 역할을 한다.
2-4. 뇌터 정리 — 대칭성과 보존량의 심오한 연결
이제 이 단계의 가장 아름다운 결과로 온다. 에미 뇌터(Emmy Noether, 1918년)가 증명한 **뇌터 정리(Noether's Theorem)**는 현대 물리학의 가장 심오한 통찰 중 하나다.
뇌터 정리: 작용 S에 연속적인 대칭성이 있으면, 반드시 대응하는 보존량이 존재한다.
"대칭성"이란 무엇인가? 물리 법칙(라그랑지안)이 어떤 변환에 대해 변하지 않는 것이다. 구체적인 예를 보자. 시간 병진 대칭(time translation symmetry): 라그랑지안이 t를 명시적으로 포함하지 않는다면, 시스템은 시간을 어제 실험하든 오늘 실험하든 같은 방식으로 행동한다. 뇌터 정리에 의하면 이 대칭성에 대응하는 보존량은 **에너지(해밀토니안 H)**다. 공간 병진 대칭(spatial translation symmetry): 라그랑지안이 어떤 좌표 qᵢ를 명시적으로 포함하지 않으면(∂L/∂qᵢ = 0), 그 qᵢ에 대응하는 일반화 운동량 pᵢ = ∂L/∂q̇ᵢ가 보존된다. 이것은 오일러-라그랑주 방정식에서 직접 보인다. 회전 대칭(rotational symmetry): 공간이 회전에 대해 균일하다면, 이에 대응하는 보존량은 각운동량이다.
[노트 기록] 뇌터 정리 대응 관계 (반드시 외울 것):
- 시간 균일성 ↔ 에너지 보존
- 공간 균일성 ↔ 선운동량 보존
- 공간 등방성(회전 대칭) ↔ 각운동량 보존 이 세 가지는 단순한 법칙이 아니라, 우주가 가진 대칭성의 필연적 결과다.
뇌터 정리를 라그랑주 역학으로 증명해보자. 만약 좌표 qᵢ가 라그랑지안에 나타나지 않는다면(이런 좌표를 순환 좌표(cyclic coordinate) 혹은 **무시 좌표(ignorable coordinate)**라 한다), 오일러-라그랑주 방정식은 d/dt(∂L/∂q̇ᵢ) = ∂L/∂qᵢ = 0이 되어, pᵢ = ∂L/∂q̇ᵢ = 상수가 된다. pᵢ가 보존된다. 이것이 뉴턴 역학에서 "외력이 없으면 운동량 보존"이라는 사실의 더 심오한 버전이다. "외력이 없다"는 것은 결국 "공간이 그 방향으로 균일하다", 즉 대칭성이 있다는 뜻이기 때문이다.
이 사상은 입자물리학까지 이어진다. 전하 보존 법칙은 **게이지 대칭(gauge symmetry)**에 대응하고, 색전하(color charge) 보존은 SU(3) 대칭에 대응한다. 우주의 모든 보존 법칙은 그 배후에 어떤 대칭성을 가진다. 뇌터 정리는 이것을 수학적으로 보장한다.
3부: 프로젝트 — 이중 진자 동역학 분석
이제 배운 것을 직접 적용할 시간이다. 아래 문제들은 정답 없이 제시된다. 각 문제에서 어떤 개념이 필요한지 먼저 스스로 생각하고, 막히는 부분이 있으면 어느 개념이 필요한지 역방향으로 추적해보라. 총 풀이 시간은 40분 내외를 목표로 한다.
[프로젝트 A] 이중 진자의 라그랑지안 유도 (약 15분)
질량 m₁, m₂이고 길이 l₁, l₂인 두 개의 막대가 연결된 이중 진자를 생각하자. 첫 번째 막대는 고정점에 매달리고, 두 번째 막대는 첫 번째 막대의 끝에 매달린다. 중력 가속도는 g, 일반화 좌표는 θ₁(첫 번째 막대와 수직선의 각도), θ₂(두 번째 막대와 수직선의 각도)다.
문제 A-1. 각 질점의 위치를 직교 좌표 (x₁, y₁), (x₂, y₂)로 표현하되, θ₁, θ₂, l₁, l₂로만 나타내라. (힌트: 고정점을 원점으로 잡고, y축을 아래 방향으로 정의하면 편하다. 아니면 위로 정의해도 된다. 어떤 선택이든 일관성을 유지해라.)
문제 A-2. 각 질점의 속도의 제곱(v₁², v₂²)을 θ₁, θ₂, θ̇₁, θ̇₂, l₁, l₂로 표현하라. (힌트: vᵢ² = ẋᵢ² + ẏᵢ²이다. A-1에서 구한 위치를 시간 미분하면 속도가 나온다. 계산이 복잡하므로 차근차근 전개하라. cos(θ₁-θ₂) 항이 나타날 것이다.)
문제 A-3. 이중 진자의 라그랑지안 L = T - V를 완성하라. 운동 에너지 T와 퍼텐셜 에너지 V를 각각 써라.
문제 A-4. 오일러-라그랑주 방정식을 θ₁과 θ₂에 대해 각각 적용하여 이중 진자의 운동 방정식 두 개를 유도하라. (경고: 계산이 상당히 복잡하다. ∂L/∂θᵢ와 ∂L/∂θ̇ᵢ를 신중하게 계산하라. 특히 d/dt[∂L/∂θ̇ᵢ]에서 곱의 미분 규칙을 주의하라.)
[프로젝트 B] 해밀토니안과 에너지 보존 (약 10분)
프로젝트 A에서 유도한 라그랑지안을 사용한다.
문제 B-1. 이중 진자의 일반화 운동량 p₁ = ∂L/∂θ̇₁, p₂ = ∂L/∂θ̇₂를 계산하라.
문제 B-2. 이 시스템의 라그랑지안이 시간 t를 명시적으로 포함하는지 확인하라. 뇌터 정리에 의하면, 시간 병진 대칭이 있을 때 무엇이 보존되는가? 이중 진자의 해밀토니안 H = T + V가 운동 중에 보존될 것이라고 예상하는가? 그 이유는 무엇인가?
문제 B-3. (심화) 이중 진자에서 각운동량이 보존되지 않는 이유를 뇌터 정리의 관점에서 설명하라. 어떤 대칭성이 깨져 있는가?
[프로젝트 C] 위상 공간 분석 (약 15분)
이제 수치적(numerical) 관점으로 넘어간다. 실제로 이중 진자를 시뮬레이션하지 않아도, 다음 사고 실험을 통해 위상 공간의 본질을 이해할 수 있다.
문제 C-1. 단진자(single pendulum)를 생각하자. 위상 공간은 (θ, p) 2차원 평면이다. 다음 세 경우에 위상 공간 궤적의 형태와 특징을 (수식 없이) 서술하라: (i) 에너지가 작아서 작은 각도로만 진동하는 경우, (ii) 에너지가 매우 커서 추가 계속 회전하는 경우, (iii) 정확히 불안정 평형점(θ = π, 추가 꼭대기)에 대응하는 에너지를 가지는 경우. 각 경우의 궤적을 직접 스케치해보라.
문제 C-2. 단진자의 에너지는 E = p²/(2mL²) - mgL cosθ로 주어진다. 이 등에너지 곡선(E = constant)이 (θ, p) 평면에서 어떻게 생겼는지 분석하라. 특히, E = 2mgL(추가 꼭대기까지 정확히 닿을 수 있는 임계 에너지)일 때 궤적이 특별한 이유를 설명하라.
문제 C-3. 이중 진자는 4차원 위상 공간 (θ₁, θ₂, p₁, p₂)를 가진다. 이 4차원 공간을 직접 시각화할 수는 없지만, **푸앵카레 단면(Poincaré Section)**이라는 기법이 있다. θ₂ = 0이고 ṗ₂ > 0인 순간마다 (θ₁, p₁) 평면에 점을 찍는 방식이다. 에너지가 낮을 때(규칙적 운동)와 에너지가 높을 때(카오스적 운동)의 푸앵카레 단면은 각각 어떤 모습일지 예측하고, 그 이유를 위상 공간의 논리로 설명하라.
[프로젝트 D] 보너스 — 뇌터 정리 적용 (추가 도전)
문제 D-1. 균일한 자기장 B (z 방향) 속에서 움직이는 전하 q를 가진 입자를 생각하자. 이 경우 라그랑지안은 L = ½mv² + q(v · A)로 주어지며, A = ½(-By, Bx, 0)는 벡터 퍼텐셜이다. 이 시스템이 z축 회전에 대해 대칭인지 확인하고, 뇌터 정리에 의해 어떤 양이 보존되는지 일반화 좌표 (r, φ, z)를 사용하여 구하라.
마무리: 다음 단계를 향해
이번 2단계에서 우리는 "왜"라는 질문에서 출발했다. 왜 뉴턴 역학이 부족한가, 왜 자연은 특정 경로를 선택하는가, 왜 에너지가 보존되는가. 변분법은 이 "왜"에 수학적 언어를 부여했고, 라그랑지안은 에너지로 운동 방정식을 도출했으며, 해밀토니안과 위상 공간은 시스템의 상태와 진화를 기하학적으로 그려냈다. 그리고 뇌터 정리는 물리학의 가장 깊은 통찰—대칭성이 보존량을 만든다—을 명확하게 제시했다.
3단계에서는 이중 진자에서 예고된 **카오스(chaos)**의 세계로 들어간다. 뇌터 정리가 보존량의 존재를 보장하는데도, 왜 이중 진자는 예측 불가능해지는가? 그리고 N개의 물체가 서로 중력으로 당기는 계에서 수치 적분은 어떻게 작동하는가? 오늘 배운 위상 공간과 해밀토니안 역학은 그 카오스 이론의 직접적인 토대가 된다. 프로젝트 C-3에서 예측한 푸앵카레 단면의 모습이 3단계에서 검증될 것이다.