실니코프 분기

계산물리 수치해석 물리 마지막 수정: 2026-08-12 04:16:28

1. 개요[편집]

실니코프 분기 Shilnikov bifurcation
제시	L. P. 실니코프, 1965
무대	3차원 이상 흐름의 안장-초점
고유값	γ > 0 (실수) · −ρ ± iω (ρ > 0)
안장 지표	δ = ρ / γ
δ < 1	가산 무한 개의 말굽 → 카오스
δ > 1	주기궤도 하나만 갈라져 나옴
여차원	1 (매개변수 하나로 도달 가능)

실니코프 분기는 3차원 이상 흐름에서 안장-초점(saddle-focus) 평형점에 호모클리닉 궤도가 붙었을 때, 그 평형점의 고유값 비율이 특정 부등식을 만족하면 궤도의 임의의 근방에 무한히 많은 스메일 말굽이 나타나는 대역 분기다. 레오니트 파블로비치 실니코프가 1965년에 증명했고, 결론을 한 줄로 줄이면 이렇다 — 고유값만 재면 카오스의 존재가 결정된다.¹

이 정리가 충격적인 이유는 요구 조건이 터무니없이 싸기 때문이다. 카오스를 보이려면 보통 궤도를 오래 적분하고 랴푸노프 지수를 재고 푸앵카레 단면을 뜯어봐야 한다. 그런데 실니코프는 평형점 하나의 야코비안 고유값 세 개와, 그 평형점에 붙은 고리 하나만 있으면 된다고 말한다. 게다가 안장-초점에 붙은 호모클리닉은 여차원 1이라 매개변수 하나만 조율하면 만나는 흔한 사건이다. “호모클리닉 고리 하나가 곧 카오스의 씨앗”이라는 이 바닥의 관용구는 정확히 이 정리를 가리킨다.

호모클리닉 궤도 자체의 정의, 횡단 교차와 스메일-버코프 정리, 멜니코프 방법, 궤도의 수치 계산법은 호모클리닉 궤도 문서가 다룬다. 여기서는 안장-초점이라는 특수한 무대에서만 벌어지는 일에 집중한다.

2. 무대: 안장-초점[편집]

$\dot{\mathbf x} = \mathbf f(\mathbf x, \mu)$ 의 평형점 $\mathbf x^*$ 에서 야코비안 고유값이

\lambda_1 = \gamma > 0, \qquad \lambda_{2,3} = -\rho \pm i\omega, \quad \rho > 0,\ \omega \neq 0

이면 이 평형점을 안장-초점이라 한다. 불안정 다양체 $W^u$ 는 1차원 곡선이고, 안정 다양체 $W^s$ 는 2차원인데 그 위에서 궤도가 나선을 그리며 감겨 들어온다. 시간을 뒤집은 배치(실수 고유값이 음수, 복소 쌍이 양수)도 똑같이 취급된다.

여기서 안장 지표(saddle index)를 정의한다.

\delta \;=\; \frac{|\mathrm{Re}\,\lambda_{2,3}|}{\lambda_1} \;=\; \frac{\rho}{\gamma}

$\delta$ 는 “빨려 들어가는 속도 대 튕겨 나가는 속도”의 비다. $\delta > 1$ 이면 수축이 팽창을 이기고, $\delta < 1$ 이면 팽창이 이긴다. 이 부호 하나가 전부를 가른다.

3. 실니코프 정리[편집]

$\mathbf x^*$ 가 안장-초점이고 여기에 호모클리닉 궤도 $\Gamma$ 가 붙어 있으며 $\delta < 1$ 이면, $\Gamma$ 의 임의로 작은 근방 안에 가산 무한 개의 말굽이 존재한다. 따라서 그 근방에는 임의로 긴 주기의 주기궤도가 무한히 많고, 비가산 개의 비주기 궤도가 있으며, 초기조건 민감성이 성립한다.

$\delta > 1$ 이면 정반대다. 호모클리닉이 깨질 때 극한주기궤도가 딱 하나 갈라져 나오고, 근방에는 그것 말고 아무 일도 없다. 평면계의 얌전한 호모클리닉 분기와 질적으로 같아진다.

주의할 점 두 가지.

정리가 주는 것은 카오스 불변집합이지 카오스 끌개가 아니다. 말굽은 측도 0의 안장형 집합이라, 실니코프 조건을 만족한다고 해서 화면에 이상 끌개가 뜬다는 보장은 없다. 실제로 관측되는 것이 과도 카오스로 끝나는 경우도 흔하다. 관측 가능한 끌개가 되려면 대역적 재주입 구조가 추가로 필요하다.
호모클리닉이 정확히 존재하는 매개변수는 측도 0이다. 정리는 $\mu = \mu_c$ 한 점에 대한 진술이지만, 실용적으로 의미가 있는 이유는 무한히 많은 말굽 중 상당수가 $\mu$ 를 조금 흔들어도 살아남기 때문이다. 말굽은 균등 쌍곡이라 구조적 안정성을 갖는다. 그래서 “호모클리닉 근처의 매개변수 구간 전체가 카오스적”이라는 실험적 관측이 정당화된다.

4. 왜 안장 지표가 1보다 작으면 말굽이 생기는가[편집]

증명의 핵심은 되돌아옴 사상을 손으로 조립해 보는 것이다. 안장-초점 근방에 두 개의 작은 단면을 놓는다. $W^s$ 를 가로지르는 입구 단면 $\Sigma_{\mathrm{in}}$ 과, $W^u$ 를 가로지르는 출구 단면 $\Sigma_{\mathrm{out}}$ .

국소 통과. 선형화 좌표에서 나선 평면의 반지름은 $e^{-\rho t}$ 로 줄고, 불안정 방향 좌표 $z$ 는 $e^{\gamma t}$ 로 자란다. $z = z_0$ 에서 들어와 $z = h$ 에서 나간다면 체류 시간은

\tau = \frac{1}{\gamma}\ln\frac{h}{z_0}

이다. 그동안 반지름은 $e^{-\rho\tau} = (z_0/h)^{\delta}$ 배로 줄고, 위상각은 $\omega\tau = \frac{\omega}{\gamma}\ln(h/z_0)$ 만큼 돌아간다. 여기가 급소다. 안장에 가까이 들어올수록 체류 시간이 로그로 길어지고, 그만큼 위상각이 끝없이 감긴다.

대역 복귀. $\Sigma_{\mathrm{out}}$ 에서 호모클리닉 고리를 따라 한 바퀴 돌아 $\Sigma_{\mathrm{in}}$ 으로 돌아오는 사상은 매끄럽고 근사적으로 아핀이다. 이 둘을 합성하면 1차원으로 줄인 되돌아옴 사상이 다음 꼴이 된다.

z_{n+1} \;=\; \mu \;+\; A\, z_n^{\,\delta}\,\sin\!\left(\frac{\omega}{\gamma}\ln\frac{1}{z_n} + \phi\right)

$\mu$ 는 호모클리닉으로부터의 이탈량(분기 매개변수), $A, \phi$ 는 상수다. 이 사상의 그래프는 $z \to 0^+$ 에서 무한히 많은 골과 마루가 로그 주기로 쌓이는 모습을 한다. 연속한 극값들의 $z$ 값 비율은 $e^{-2\pi\gamma/\omega}$ 로 기하급수적으로 촘촘해진다.

기울기를 보면 결론이 나온다. $|dz_{n+1}/dz_n| \sim z_n^{\delta-1}$ 이므로

$\delta < 1$ — 기울기가 $z \to 0$ 에서 발산한다. 즉 0 근처의 무한히 많은 진동 구간마다 기울기 절댓값이 1을 넘는 단조 조각이 생기고, 각 조각이 자기 정의역을 통째로 덮으며 되돌아온다. 이것이 말굽의 조건이다. 조각이 무한히 많으니 말굽도 가산 무한 개.
$\delta > 1$ — 기울기가 0으로 수렴한다. 사상이 축소사상이 되어 고정점 하나로 얌전히 수렴한다. 고정점 = 극한주기궤도 하나.

“나선으로 감겨 들어오는데 감기는 것보다 튕겨 나가는 게 빠르다” 는 것이 실니코프 조건의 물리적 번역이다. 궤도가 되돌아올 때마다 조금씩 다른 위상으로 재주입되고, 그 위상 어긋남이 늘림-접힘을 만든다. 안장-초점이 아니라 실수 고유값만 있는 안장이면 나선이 없어 위상 어긋남 자체가 없고, 그래서 3차원이어도 카오스가 안 나온다.

5. 분기 다이어그램의 뱀[편집]

$\mu$ 를 $\mu_c$ 근처에서 쓸어 보면 주기궤도 가지가 로그 주기로 무한히 꺾이는 뱀 모양(Shilnikov wiggle)을 그린다. 주기 $T$ 는 안장 체류 시간이 지배하므로 $T \sim -\frac{1}{\gamma}\ln|\mu - \mu_c|$ 로 발산하고, 그 사이에 다음이 무한히 쌓인다.

안장-마디 분기와 주기배가 분기가 번갈아 나오는 접힘. 접힘의 진폭은 $\mu$ 축에서 $e^{-2\pi\rho/\omega}$ 비율로 줄어든다.
다중 펄스 호모클리닉 — 고리를 $N$ 번 감았다 돌아오는 $N$ -호모클리닉 궤도가 주 호모클리닉에 축적한다. 글렌디닝과 스패로우(1984)가 이 축적 구조를 체계적으로 정리했다.
각 접힘마다 주기배가 캐스케이드가 붙어 있어서, 분기 이론에서 말하는 파이겐바움 시나리오가 무한히 반복된다.

$1/2 < \delta < 1$ 과 $\delta < 1/2$ 는 결이 또 다르다. 후자에서는 다양체 사이에 접점(tangency)이 지속적으로 생겨 뉴하우스 현상까지 겹친다는 것이 알려져 있고(오브샨니코프-실니코프), 이 영역에서는 “무한히 많은 안정 주기궤도가 공존하는” 병리적 상황도 배제할 수 없다.²

6. 시계열에서 알아보기[편집]

실니코프 카오스는 눈으로 알아볼 수 있는 시계열 지문을 가진다. 궤도가 안장-초점 근처로 들어가면 나선을 그리며 작은 진폭으로 여러 번 진동하다가, 불안정 방향으로 충분히 자라면 대역 고리를 한 바퀴 도는 큰 스파이크를 한 번 치고 다시 돌아온다. 즉 신호가 “작은 진동 $n$ 번 + 큰 스파이크 한 번”의 반복인데, $n$ 이 매번 다르다. 이 감김 횟수의 불규칙성이 무작위성의 원천이며, 신경계 버스팅에서 버스트당 스파이크 개수가 요동치는 현상과 같은 골격이다.

정량적 지문 셋을 순서대로 확인하면 된다.

체류 시간 히스토그램의 로그 주기 구조. 스파이크 간 간격 분포에 봉우리가 여러 개 나타나고, 그 간격이 나선 주기 $2\pi/\omega$ 다. 한 바퀴 더 감기면 그만큼 더 머문다는 사실이 히스토그램에 그대로 찍힌다.
되돌아옴 사상의 재구성. 푸앵카레 단면에서 연속한 교점을 $(z_n, z_{n+1})$ 로 찍는다. 앞 절의 $z^{\delta}\sin(\cdot\ln z)$ 그래프가 봉우리 여러 개가 원점 쪽으로 촘촘해지는 형태로 실제로 눈에 보이면 강한 증거다. 반대로 봉우리가 하나뿐인 단봉 사상이 나오면 그건 파이겐바움 계열이지 실니코프가 아니다.
주기의 로그 발산. 매개변수를 호모클리닉 쪽으로 밀 때 주기가 $-\ln|\mu-\mu_c|$ 로 자라면 호모클리닉, $|\mu-\mu_c|^{-1/2}$ 로 자라면 안장-마디 계열이다. 호모클리닉 궤도 문서의 스케일링 판별이 여기서도 1차 관문 역할을 한다.

실험 데이터만 있고 방정식이 없을 때도 이 셋은 쓸 수 있다. 지연 좌표로 상공간을 재구성한 뒤 단면을 잘라 되돌아옴 사상을 그리는 것이 화학·레이저·전기생리 실험에서 실니코프를 주장하는 표준 절차다.

7. 사촌들[편집]

안장-초점 호모클리닉 하나만 놓고 봐도 변주가 여럿이라, 이름을 섞어 쓰면 곧바로 오해가 생긴다.

구조	평형점의 고유값	결과
실니코프 (안장-초점 호모클리닉)	실수 하나 + 복소 쌍, $\delta < 1$	가산 무한 개 말굽
같은 구조, $\delta > 1$	실수 하나 + 복소 쌍, $\delta > 1$	극한주기궤도 하나
실수 안장 호모클리닉	실수 고유값만	평면형 대역 분기, 카오스 없음
SNIC (안장-마디 위 불변원)	실수 고유값 하나가 0	주기가 제곱근으로 발산, 카오스 없음

여차원 2로 올라가면 조합이 늘어난다. 안장-초점의 호모클리닉과 호프 분기가 같은 지점에서 만나는 실니코프-호프, 두 개의 안장-초점을 잇는 헤테로클리닉 사이클인 비코프 사이클이 대표적이다. 후자는 두 안장의 지표를 곱한 조건이 판정 기준이 되고, 매개변수 평면에서 실니코프 곡선들이 무한히 축적하는 복잡한 지도를 만든다. 셋 다 두 매개변수 지도를 그려야만 보이는 대상이라, 매개변수 하나만 쓸어보는 실험에서는 정체를 알 수 없다.

8. 실제로 어디서 나오나[편집]

추아 회로. 다이오드 하나 저항 둘 인덕터 하나 커패시터 둘로 만드는 추아 회로는 조각별 선형 3차원 계이고, 그 이중 스크롤 끌개는 원점 부근 안장-초점에 붙은 호모클리닉으로 설명된다. 추아·고무로·마쓰모토(1986)가 실니코프 정리를 적용해 이 회로의 카오스를 논증했다. 조각별 선형이라 각 영역의 고유값을 해석적으로 손으로 쓸 수 있어서 $\delta$ 를 종이 위에서 계산할 수 있다 — 실니코프 조건이 검증 가능한 형태로 확인되는 드문 사례다.

나선형 카오스. 뢰슬러 방정식처럼 “평면에서 나선을 그리며 커지다가 한 번 튕겨 되돌아오는” 구조의 계는 전형적인 실니코프 후보다. 나선 부분이 안장-초점의 $W^s$ , 튕김이 대역 복귀에 해당한다. 화학 진동 반응(벨루소프-자보틴스키 계열 모형), CO 산화 반응, 글로 방전, 광펌핑 레이저, 신경 모형의 버스팅에서도 같은 골격이 보고된다.

로렌츠 방정식은 실니코프가 아니다. 자주 헷갈리는 부분인데, 고전적 로렌츠 끌개를 만드는 $\rho \approx 13.926$ 의 호모클리닉은 실수 고유값만 가진 안장(원점)에 붙은 것이라 실니코프 시나리오가 아니라 다른 메커니즘(대역 폭발)이다. 다만 다른 매개변수 영역에서는 로렌츠 계에도 안장-초점 호모클리닉이 존재한다.

9. 수치적으로 판정하기[편집]

절차는 의외로 짧다. 궤도 계산 자체(절단 + 사영 경계조건 BVP, AUTO의 HomCont, 린의 방법)는 호모클리닉 궤도 문서에 있으니 여기서는 판정에 필요한 것만 적는다.

평형점을 다 찾고 고유값을 뽑는다. 실수 하나 + 복소 켤레 쌍의 조합이 있는지 본다. 없으면 실니코프는 애초에 후보가 아니다.
$\delta = \rho/\gamma$ 를 계산한다. 여기서 흔한 실수 — 복소 쌍의 절댓값이 아니라 실부를 써야 한다. $|\lambda_{2,3}| = \sqrt{\rho^2+\omega^2}$ 를 넣으면 $\omega$ 가 큰 계에서 $\delta$ 가 과대평가되어 멀쩡한 카오스를 “실니코프 아님”으로 오판한다.
호모클리닉이 실제로 있는지 확인한다. $\delta$ 만으로는 아무것도 증명되지 않는다. 고유값은 국소 정보이고, 호모클리닉의 존재는 대역 정보다. 이 단계를 건너뛰고 ” $\delta < 1$ 이니 실니코프 카오스”라고 쓴 논문이 실재하며, 그건 문 열쇠만 보여주고 문이 있다고 주장하는 격이다.³
수치 연속법으로 매개변수를 따라간다. 두 매개변수 평면에서 호모클리닉 곡선을 추적하면, $\delta = 1$ 인 선(벨류코프 경계)이 그 곡선을 가로지르는 지점이 나온다. 그 지점을 기준으로 한쪽은 얌전한 극한주기궤도, 다른 쪽은 카오스다. 이 교차점이 실니코프 해석의 실질적 결과물이다.

실무적 함정도 있다. 안장 근처 체류 시간이 로그로 발산하므로 적분 구간이 길어지고, 나선 방향의 강한 수축 때문에 슈팅 방식은 조건수가 지수적으로 무너진다. 그리고 $\delta$ 가 1에 가까우면( $0.9 \sim 1.1$ ) 판정이 사실상 불가능하다 — 그 정도 차이는 평형점 위치의 수렴 오차와 야코비안의 유한차분 오차에 묻힌다.

10. 관련 문서[편집]

러시아어 Шильников 를 옮긴 표기가 Shilnikov · Shil’nikov · Šilnikov 로 제각각이라, 검색할 때 셋 다 넣어야 논문이 다 나온다. 원 논문은 1965년 Soviet Math. Dokl. 에 실린 짧은 노트이며 제목이 “가산 개의 주기 운동이 존재하는 한 경우”다. 제목만 봐도 결론이 다 적혀 있는 보기 드문 논문. ↩
실니코프 조건을 “카오스가 보인다”로 번역하면 절반만 맞다. 정리가 보증하는 것은 말굽의 존재이고, 말굽은 대체로 끌개가 아니다. 화면에 안 나온다고 정리가 틀린 게 아니라, 당신이 보고 있는 것이 유역의 다른 부분일 뿐이다. 물론 이 변명은 리뷰어를 설득하지 못한다. ↩
반대 방향 오용도 있다. “우리 계에 안장-초점이 있고 카오스가 관측되므로 실니코프 메커니즘이다”는 삼단논법도 아니다. 3차원 카오스로 가는 길은 주기배가, 준주기 붕괴, 간헐성, 안장-마디 등 여럿이고 실니코프는 그중 하나다. 메커니즘을 주장하려면 호모클리닉 곡선을 실제로 연속해서 보여야 한다. ↩