뇌터 정리

1. 개요[편집]

뇌터 정리(Noether’s theorem)는 물리계의 작용(action)이 연속 대칭을 가지면 그 대칭 하나마다 보존량 하나가 대응한다는 정리다. 1918년 에미 뇌터(Emmy Noether)가 논문 Invariante Variationsprobleme에서 증명했고, 그 뒤로 고전역학·장론·입자물리 전부를 관통하는 뼈대가 되었다.

이 정리가 대단한 이유는 보존법칙을 “관측된 사실”에서 “구조의 필연”으로 승격시켰다는 데 있다. 에너지가 보존되는 이유는 실험에서 그렇더라가 아니라, 물리 법칙이 어제와 오늘 같기 때문이다. 운동량이 보존되는 이유는 여기와 저기에서 법칙이 같기 때문이다.¹

수치해석·시뮬레이션 쪽 사람이 이 문서를 읽어야 할 이유는 따로 있다. 이산화는 대칭을 깨는 행위이기 때문이다. 격자를 놓는 순간 공간 병진 대칭이 깨지고, 시간 스텝 $h$ 를 고정하는 순간 시간 병진 대칭이 깨진다. 그러면 뇌터 정리에 따라 대응하는 보존량이 정확히 보존되지 않는다. 반대로 어떤 대칭을 이산 수준에서도 살려 두면 그 보존량은 반올림 오차 수준까지 정확히 보존된다. 이것이 베를레 적분 계열이 룽게-쿠타법보다 장시간 적분에서 강한 이유이자, 유한체적법이 산업 코드를 장악한 이유다.

대칭·정준구조 자체의 이론은 해밀토니안 역학, 라그랑주 역학, 포아송 괄호, 리우빌 정리 문서에 있으니 여기서는 반복하지 않는다. 이 문서는 대칭 ↔ 보존량 대응과 그것의 이산 버전에 집중한다.

2. 대칭과 보존량의 대응표[편집]

연속 대칭	보존량
시간 병진 (법칙이 시각에 무관)	에너지
공간 병진 (법칙이 위치에 무관)	선운동량
회전 (법칙이 방향에 무관)	각운동량
갈릴레이 부스트	질량중심의 등속 운동
전역 U(1) 위상 대칭	전하
스케일 대칭 (특수한 경우)	비리얼 유형 양

주의할 점 하나. 표의 마지막 줄 근처에 있는 게이지 이야기는 자주 뭉개져서 전달되는데, 전하 보존을 주는 것은 전역 U(1) 위상 대칭이고, 국소 게이지 대칭이 주는 것은 보존량이 아니라 운동방정식들 사이의 항등식(비앙키 항등식류)이다. 전자가 뇌터의 제1정리, 후자가 제2정리다.²

3. 유도 스케치[편집]

라그랑지안 $L(q, \dot q, t)$ 가 무한소 변환 $q \to q + \varepsilon K(q,t)$ 아래에서 불변이라고 하자(엄밀히는 전미분 $\mathrm{d}F/\mathrm{d}t$ 만큼의 차이는 허용된다). 변분은

$$\delta L = \frac{\partial L}{\partial q}\varepsilon K + \frac{\partial L}{\partial \dot q}\varepsilon \dot K$$

인데, 여기에 오일러-라그랑주 방정식 $\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial L}{\partial \dot q} = \frac{\partial L}{\partial q}$ 를 대입하면 두 항이 곱미분 형태로 합쳐진다.

$$\delta L = \varepsilon \frac{\mathrm{d}}{\mathrm{d}t}\!\left( \frac{\partial L}{\partial \dot q} K \right)$$

즉 $\delta L = \varepsilon \,\mathrm{d}F/\mathrm{d}t$ 라는 대칭 조건은 곧바로

$$\frac{\mathrm{d}}{\mathrm{d}t}\!\left( \frac{\partial L}{\partial \dot q} K - F \right) = 0$$

을 뜻한다. 괄호 안이 보존량, 즉 뇌터 전하다. 병진 대칭이면 $K = 1$ 이라 보존량이 정준운동량 $\partial L/\partial \dot q$ 자체가 되고, 시간 병진이면 보존량이 해밀토니안이 된다. 장론으로 넘어가면 같은 논리가 국소 보존형 $\partial_\mu j^\mu = 0$ 을 낳는데, 이 발산 형태가 바로 다음 절에서 수치해석과 연결되는 고리다.

4. 이산 뇌터 정리와 변분 적분기[편집]

여기서 계산 쪽의 결정적인 아이디어가 나온다. 운동방정식을 이산화하지 말고, 작용을 이산화하자.

시간 격자 $t_k = kh$ 위에서 이산 라그랑지안 $L_d(q_k, q_{k+1}) \approx \int_{t_k}^{t_{k+1}} L \,\mathrm{d}t$ 를 정의하고, 이산 작용합 $S_d = \sum_k L_d$ 에 변분원리를 적용하면 이산 오일러-라그랑주 방정식이 나온다. 이렇게 만든 적분기를 변분 적분기(variational integrator)라 부르며, 여기에는 뇌터 정리의 이산 버전이 그대로 따라온다. $L_d$ 가 어떤 리 군 작용에 대해 불변이면 대응하는 이산 운동량 사상이 정확히 보존된다. 근사적으로가 아니라, 부동소수점 반올림을 빼면 정확히.

• $L_d$ 가 공간 병진 불변 → 총 선운동량이 매 스텝 정확히 보존
• $L_d$ 가 회전 불변 → 총 각운동량이 매 스텝 정확히 보존

실제로 베를레 적분(velocity-Verlet, 립프로그)은 사다리꼴 이산 라그랑지안에서 유도되는 가장 단순한 변분 적분기이며, 그래서 중심력장에서 각운동량을 정확히 지킨다. 4차 룽게-쿠타법은 국소 오차 차수는 훨씬 높지만 이런 구조적 보존은 하나도 하지 못한다. 정확도와 보존은 다른 축이라는 게 기하 수치적분(geometric numerical integration)의 핵심 교훈이다.³

5. 에너지만 예외인 이유 — 그림자 해밀토니안[편집]

위 목록에서 눈치챘겠지만 에너지가 빠져 있다. 이유는 정확히 뇌터적이다. 스텝 크기 $h$ 를 고정한 이산 시간축은 시간 병진 대칭을 깨기 때문에 에너지에 대응하는 정확한 이산 보존량이 없다.

그런데 심플렉틱 적분기로 궤도를 오래 돌려 보면 에너지가 표류하지 않고 작은 폭으로 진동만 한다. 이걸 설명하는 것이 후방 오차 해석(backward error analysis)이다. 심플렉틱 적분기가 만드는 이산 사상은, 원래 해밀토니안 $H$ 가 아니라

$$\tilde H = H + h^2 H_2 + h^4 H_4 + \cdots$$

라는 그림자 해밀토니안(shadow Hamiltonian)을 갖는 어떤 연속 해밀토니안 계의 정확한 흐름에 지수적으로 가깝다. 즉 이 적분기는 “내가 풀려던 문제를 대충 푸는” 것이 아니라 “이웃한 다른 해밀토니안 문제를 거의 정확히 푸는” 중이다. $\tilde H$ 는 그 계에서 진짜 보존되므로, 실제 $H$ 는 $\tilde H$ 주변에서 $O(h^2)$ 크기로 진동만 하고 장기 표류가 없다. 반면 비심플렉틱 적분기는 이런 그림자가 없어서 에너지가 시간에 비례해 줄줄 새거나 부푼다.

여기에 붙는 실무 함정도 두 개다. (1) 이 성질은 고정 스텝에서만 성립한다. 스텝 크기를 오차 추정으로 적응 조절하면 그 순간 심플렉틱 구조가 깨져 에너지가 표류하기 시작한다. (2) 서모스탯처럼 일부러 계에 에너지를 넣고 빼는 분자동역학 설정에서는 애초에 보존이 목표가 아니다. 이때 지켜야 할 것은 에너지가 아니라 표본이 따르는 분포다(리우빌 정리).

6. 왜 보존형 이산화를 고집하는가[편집]

같은 사고방식이 편미분방정식 쪽에도 그대로 이식된다. 유체·전자기·탄성의 지배 방정식은 대부분 발산(divergence) 형태

$$\frac{\partial u}{\partial t} + \nabla \cdot \mathbf{F}(u) = 0$$

으로 쓸 수 있는데, 이 형태 자체가 연속 대칭에서 나온 국소 보존법칙이다. 유한체적법은 이 식을 검사체적마다 적분한 뒤 면 플럭스를 이웃 셀과 하나의 값으로 공유하도록 이산화한다. 그러면 내부 면의 기여가 텔레스코핑으로 정확히 상쇄되어 전체 영역에서 질량·운동량·에너지가 반올림 수준까지 보존된다. 격자가 아무리 거칠어도 이 성질은 안 깨진다.

이게 취향 문제가 아니라는 것을 못 박은 것이 락스-벤드로프 정리다. 보존형이면서 일관된(consistent) 도식이 수렴한다면 그 극한은 반드시 약해(weak solution)이며, 따라서 충격파의 전파 속도가 랭킨-위고니오 관계를 만족한다. 뒤집으면, 비보존형으로 이산화한 도식은 그럴싸하게 수렴해도 충격파를 틀린 속도로 옮긴다. 압축성 유동과 리만 솔버 문서가 집요하게 플럭스 형태를 따지는 이유가 여기 있다.

비압축성 유동에서도 비슷한 이야기가 있다. 대류항을 왜대칭(skew-symmetric) 형태로 이산화하면 대류가 운동에너지를 생성하지도 소멸시키지도 않게 되어, 수치 점성이 물리 점성을 오염시키지 않는다. 직접수치모사나 대와류 모사처럼 에너지 수지 자체가 결과인 계산에서는 사실상 필수 조건이다.

정리하면 이렇다. 시뮬레이션 코드를 짤 때 물어야 할 질문은 “이 도식의 정확도 차수가 몇 차인가”만이 아니다. “내가 어떤 대칭을 깨고 있고, 그래서 어떤 보존량을 잃는가” 가 장시간·다물리 계산에서는 훨씬 중요한 질문이다. 뇌터 정리는 그 질문을 던지는 법을 알려준다.

7. 관련 문서[편집]

• 해밀토니안 역학 · 라그랑주 역학 · 포아송 괄호 · 리우빌 정리
• 베를레 적분 · 룽게-쿠타법 · 변분법
• 유한체적법 · 리만 솔버 · 충격파 · 압축성 유동
• 분자동역학 · 서모스탯 · 직접수치모사
• 적분가능계 · 게이지 이론 · 회전행렬

8. Footnotes[편집]

1. 뒤집으면 무서운 이야기도 된다. 팽창하는 우주는 시간 병진 대칭을 갖지 않으므로, 전체 에너지 보존이라는 개념이 일반적으로 성립하지 않는다. “에너지는 어디 갔나요”라는 질문에 우주론자가 “그 질문이 잘 정의되지 않습니다”라고 답하는 근거가 이것이다. ↩

2. 뇌터의 원 논문에 정리가 두 개 실려 있는데, 물리 교과서 대부분이 제1정리만 다루고 지나간다. 정작 제2정리는 게이지 이론의 구속 구조를 이해하는 열쇠라, 나중에 정준 양자화를 하려는 사람이 뒤늦게 되돌아와 읽는 부분이다. ↩

3. 그래서 태양계 장기 적분이나 분자동역학에서 “RK4가 2차 립프로그보다 좋겠지”라고 생각하고 갈아탔다가 며칠 뒤 에너지가 우상향하는 그래프를 보게 되는 것은 이 바닥의 통과의례다. 정확도는 높은데 궤도가 나선으로 풀려 나간다. ↩