푸앵카레 단면

1. 개요[편집]

푸앵카레 단면(Poincaré section)은 연속시간 흐름의 궤적이 위상공간에 놓인 어떤 초곡면을 정해진 방향으로 가로지르는 순간의 점만 골라 찍은 것이고, 그렇게 찍힌 점에서 다음 점으로 가는 대응 규칙이 되돌이 사상(first-return map, 푸앵카레 사상)이다. 앙리 푸앵카레가 3체 문제를 붙들고 있던 1890년대에 고안했으며, 미분방정식을 푸는 대신 그 해의 모양을 보겠다는 정성적 동역학의 출발점이다.

핵심은 두 가지다. 첫째, 차원이 하나 준다 — $n$ 차원 흐름이 $n-1$ 차원 사상이 된다. 둘째, 시간이 사라진다 — 미분방정식이 이산 사상으로 바뀌면서, 흐름을 따라가는 지루한 나선 감기가 점 몇 개로 압축된다. 로렌츠 끌개를 아무리 3D로 회전시켜 봐야 실뭉치처럼 보일 뿐이지만, 단면을 하나 끼워 넣는 순간 거의 1차원짜리 곡선 위의 점들이 튀어나오고 그제야 구조가 읽힌다.¹

2. 횡단 단면과 되돌이 사상[편집]

흐름 $\dot{\mathbf{x}} = \mathbf{F}(\mathbf{x})$, $\mathbf{x}\in\mathbb{R}^n$ 에 대해 초곡면을 $\Sigma = \{\mathbf{x} : g(\mathbf{x}) = 0\}$ 로 잡는다. 이때 반드시 요구되는 것이 횡단성 조건이다.

$$\nabla g(\mathbf{x}) \cdot \mathbf{F}(\mathbf{x}) \neq 0 \qquad \text{on } \Sigma$$

벡터장이 단면에 접하면 안 된다는 뜻이다. 접하는 순간 궤적이 단면을 스치고 되돌아가 버려서 “가로질렀다”는 판정 자체가 무의미해지고, 뒤에 나올 사상이 특이해진다. 실무에서 단면을 잘못 고르면 대개 이 조건이 국소적으로 깨진 것이다.

$\Sigma$ 위의 점 $\mathbf{p}$ 에서 출발한 궤적이 같은 부호 방향으로 처음 다시 $\Sigma$ 에 닿는 시각을 되돌이 시간 $\tau(\mathbf{p})$ 라 하면, 되돌이 사상은

$$P(\mathbf{p}) = \boldsymbol{\phi}_{\tau(\mathbf{p})}(\mathbf{p}), \qquad P: \Sigma \to \Sigma$$

이다. 되돌이 시간이 점마다 다르다는 것이 이 정의의 특징이자 골칫거리다. 방향을 하나로 못 박는 것도 필수인데, 양방향 교차를 다 찍으면 서로 다른 두 장의 사상이 겹쳐 그림이 뭉개진다.²

3. 주기궤도는 사상의 고정점[편집]

가장 유용한 성질은 이것이다. 흐름의 주기궤도가 $\Sigma$ 를 $m$ 번 가로지른다면, 그 궤도는 $P$ 의 주기-$m$ 점에 대응한다. 특히 한 번만 가로지르면 고정점 $P(\mathbf{p}^*) = \mathbf{p}^*$ 다.

덕분에 “주기해를 찾아라”라는 경계값 문제가 ”$P(\mathbf{p}) - \mathbf{p} = 0$ 을 풀어라”라는 비선형 대수방정식으로 바뀐다. 뉴턴-랩슨법을 얹으면 그것이 곧 슈팅법 기반 주기해 추적이고, 파라미터를 조금씩 밀면서 이 해를 따라가는 것이 수치 연속법이다. 분기 이론에서 주기궤도 가지를 그리는 도구가 정확히 이 조합이다.

4. 플로케 승수와 안정성[편집]

안정성은 사상의 야코비안 $DP(\mathbf{p}^*)$ 의 고윳값이 결정한다. 이 고윳값을 흐름 쪽 언어로 옮긴 것이 플로케 승수다. 궤도를 따라 변분방정식

$$\dot{\Phi} = D\mathbf{F}(\mathbf{x}(t))\,\Phi, \qquad \Phi(0) = I$$

를 한 주기 $T$ 만큼 적분해 얻은 $M = \Phi(T)$ 가 모노드로미 행렬이고, 그 고윳값 $\mu_i$ 가 플로케 승수다. 자율계에서는 흐름 방향에 대응하는 승수가 항상 정확히 1이며(궤도를 따라 밀어도 같은 궤도), 이 자명한 승수를 제거한 나머지 $n-1$ 개가 곧 $DP$ 의 고윳값이다. 전부 $|\mu|<1$ 이면 궤도가 점근 안정.

승수가 단위원을 뚫고 나가는 방식이 다음에 벌어질 일을 알려 준다.

이탈 방식	분기	단면에서 보이는 변화
$\mu \to +1$	순환의 안장-마디	두 고정점이 충돌해 소멸
$\mu \to -1$	주기배가	점 하나가 번갈아 찍히는 두 점으로
복소켤레 쌍	나이마크-사커	고정점이 작은 닫힌 곡선으로

$\det M = \exp\!\int_0^T \nabla\cdot\mathbf{F}\,dt$ 라는 리우빌 공식도 유용한 검산이다. 소산계면 이 값이 1보다 작아 사상이 넓이를 줄이고, 보존계라면 정확히 1이라 $P$ 가 넓이 보존 사상이 된다. 계산한 승수 중에 1에 충분히 가까운 것이 없다면 적분 정확도부터 의심해야 한다.

5. 스트로보스코픽 단면[편집]

주기 $T$ 로 강제되는 비자율계 $\dot{\mathbf{x}} = \mathbf{F}(\mathbf{x}, t)$, $\mathbf{F}(\mathbf{x},t+T)=\mathbf{F}(\mathbf{x},t)$ 에서는 훨씬 편한 선택지가 있다. 시간을 상태변수로 승격시켜 확장 위상공간 $(\mathbf{x}, \theta)$, $\theta = t \bmod T$ 를 만들면 계가 자율계가 되고, $\theta = \text{const}$ 인 초평면이 자연스러운 단면이 된다. 결과적으로 매 $T$ 마다 상태를 한 장씩 찍는 것이 전부다 — 회전하는 물체에 스트로보스코프를 비추는 것과 정확히 같아서 이 이름이 붙었다.

되돌이 시간이 상수 $T$ 로 고정되므로 교차 검출도, 보간도 필요 없다. 그래서 강제 더핑 진동자, 강제 진자, 주기 강제를 받는 유동 혼합 문제에서는 단면 계산이 사실상 공짜다. 시간주기 유동장 속 입자 궤적을 시간 $T$ 마다 찍어 라그랑주 카오스와 혼합 영역을 읽는 것이 이 기법의 유체 쪽 대표 용례다.

6. 단면 위의 그림 읽기[편집]

찍힌 점들의 모양만 봐도 동역학의 정체가 거의 드러난다.

• 유한 개의 점 — 주기해. 점이 $m$ 개면 되돌이 사상 기준 주기 $m$, 즉 주기배가를 $\log_2 m$ 번 겪은 상태.³
• 닫힌 곡선 — 준주기해. 위상공간의 2-토러스가 단면에 남긴 자국이다. 곡선 위에서의 이동각을 한 바퀴 $2\pi$ 로 나눈 회전수가 무리수면 궤도가 곡선을 조밀하게 채우고, 유리수 $p/q$ 면 $q$ 개의 점으로 떨어진다(주파수 잠김).
• 프랙탈 점구름 — 카오스. 로렌츠 끌개의 단면은 얼핏 매끄러운 곡선처럼 보이지만, 확대하면 칸토어 집합 같은 층상 구조가 끝없이 나온다. 부피가 0이면서 차원이 정수가 아닌 프랙탈이라는 사실이 여기서 눈으로 확인된다.
• 면을 균일하게 뒤덮는 점 — 에르고딕한 카오스 바다. 보존계에서 흔하다.

정량 지표도 그대로 이어진다. 흐름의 랴푸노프 지수 스펙트럼에는 흐름 방향에 해당하는 $\lambda=0$ 이 반드시 하나 들어 있는데, 단면을 취하면 그 성분이 제거된다. 즉 단면 사상의 지수 집합 = 흐름의 지수 집합에서 0 하나를 뺀 것이며, 사상 쪽에서 $\lambda>0$ 이 남으면 카오스다.

7. 해밀토니안계 — KAM 토러스와 섬 사슬[편집]

해밀토니안 역학에서 이 도구는 특히 강력하다. 자유도 2인 계는 위상공간이 4차원이지만 에너지 보존으로 3차원 등에너지면에 갇히고, 단면을 취하면 2차원 넓이 보존 사상이 된다. 종이 한 장에 계의 모든 것이 담기는 셈이다.

적분가능계라면 궤도가 불변 토러스 위에 있고, 단면에는 동심원 비슷한 닫힌 곡선들만 나온다. 여기에 작은 섭동을 가하면 어떻게 되는가 — 이것이 KAM 정리의 무대다. 회전수가 유리수에서 충분히 멀다는 디오판토스 조건을 만족하는 토러스는 살짝 일그러진 채 살아남고, 단면에서 KAM 곡선으로 보인다. 반면 회전수가 유리수인 공명 토러스는 부서지는데, 푸앵카레-비르호프 정리에 따라 그 자리에 짝수 개의 주기점이 남는다. 절반은 타원형(주변에 작은 닫힌 곡선을 거느린 섬), 절반은 쌍곡형(안장)이다.

쌍곡점의 안정 다양체와 불안정 다양체가 한 번 횡단 교차하면 무한히 많이 교차할 수밖에 없고, 그 결과 생기는 것이 푸앵카레가 “너무 복잡해서 그리려는 시도조차 하지 않겠다”고 적은 호모클리닉 얽힘이다. 단면에서는 섬 사슬 사이를 채우는 카오스 층으로 나타나며, 섬 하나를 확대하면 그 안에 또 섬 사슬이 나오는 자기유사 구조가 계속된다. 섭동을 키우면 KAM 곡선이 하나씩 무너지고 카오스 바다가 연결되어 전역 수송이 열린다. 에농-하일레스 계와 표준 사상이 이 시나리오의 교과서 예제이며, 가속기의 빔 동역학에서 동적 구경(dynamic aperture)을 정하는 계산도 결국 1턴 사상의 단면을 그리는 일이다.

8. 수치 구현에서 조심할 것[편집]

교차 검출 자체는 쉽다. 매 스텝 $g(\mathbf{x}_n)$ 의 부호 변화를 감시하고 $\nabla g\cdot\mathbf{F}$ 로 방향을 걸러내면 된다. 문제는 교차점을 얼마나 정확히 잡느냐다.

• 선형 보간은 적분기 차수를 깎는다. 4차 룽게-쿠타법으로 열심히 적분해 놓고 마지막 한 스텝을 두 점 사이 직선으로 이으면, 단면 위 좌표의 정확도가 1차로 떨어진다. 이 오차는 사상을 반복할수록 누적되어 끌개 모양을 흐리고, 회전수와 플로케 승수를 오염시킨다.

• 밀집 출력(dense output) 을 제공하는 적분기라면 그 연속 확장으로 $g=0$ 인 시각을 뉴턴-랩슨법으로 풀면 된다. 적분기 차수가 그대로 유지된다.

• 에농의 요령(1982)이 가장 우아하다. 부호가 바뀐 것을 확인한 뒤 마지막 한 스텝에서 독립변수를 시간 $t$ 에서 $g$ 자체로 바꿔치기한다.

  $$\frac{d\mathbf{x}}{dg} = \frac{\mathbf{F}(\mathbf{x})}{\nabla g\cdot\mathbf{F}(\mathbf{x})}$$

  이 방정식을 현재 $g$ 값에서 0까지 정확히 한 스텝 적분하면 교차점에 정확히 도달한다. 보간도, 반복도 필요 없고 적분기 차수도 보존된다. 분모가 횡단성 조건 그 자체이므로, 이 요령이 깨지는 자리는 애초에 단면 선택이 나쁜 자리다.

• 표본 수가 통계를 정한다. 단면 위의 점은 교차할 때만 하나씩 생긴다. 되돌이 시간이 긴 계에서 점 1만 개를 모으려면 궤적을 그만큼 길게 적분해야 하고, 그동안 카오스계에서는 궤적이 이미 참 해와 무관해져 있다. 그래도 통계량은 살아남는다는 것이 그림자 정리의 위안이다.

• 실험 데이터에서는 단면을 직접 정의할 수 없으므로, 관측 시계열의 국소 극대점을 되돌이 사상으로 삼거나(로렌츠계의 $z$ 최댓값 사상이 텐트 사상 모양이 되는 것이 유명하다) 힐베르트 변환으로 위상을 뽑아 $\theta = 0$ 을 단면으로 쓴다.

9. 관련 문서[편집]

• 카오스 이론 · 랴푸노프 지수
• 분기 이론 · 극한 순환
• 해밀토니안 역학 · 리우빌 정리 · 포아송 괄호
• 프랙탈 · 다중프랙탈
• 룽게-쿠타법 · 변분 적분기 · 뉴턴-랩슨법
• 자코비안 행렬 · 고유값 문제

10. Footnotes[편집]

1. 푸앵카레가 이 아이디어를 얻은 계기가 오스카르 2세 상 논문의 오류를 수습하는 과정이었다는 건 유명한 일화다. 이미 인쇄된 논문을 회수하고 자비로 다시 찍었는데, 그 수정판에서 호모클리닉 얽힘을 발견했다. 오타 수정하다 카오스를 발견한 셈. ↩

2. 단면 위치를 옮기면 사상 자체는 바뀌지만 고정점의 플로케 승수, 회전수, 랴푸노프 지수 같은 불변량은 그대로다. 그림이 예뻐지는 위치를 마음껏 찾아도 된다는 뜻이고, 실제로 논문 그림 하나 뽑으려고 단면을 열 번 옮기는 것이 이 바닥의 흔한 풍경이다. ↩

3. “단면에 점이 딱 세 개 찍혔으니 주기 3이다”라고 결론 내리기 전에, 그게 정말 수렴한 상태인지 확인해야 한다. 과도 구간을 충분히 버리지 않으면 아직 감겨 들어가는 중인 궤적이 유령 점을 뿌린다. 국룰은 앞쪽 절반을 통째로 버리는 것. ↩