카피차 진자

1. 개요

카피차 진자 Kapitza pendulum
구성	지지점을 연직으로 고속 진동시킨 강체 진자
현상	거꾸로 선 평형점이 안정해짐
이론	A. Stephenson (1908) · P. Kapitza (1951)
안정 조건	(Aω)² > 2gl
상한	A/l ≲ 0.45 (주기배가로 재불안정)
도구	빠른/느린 분리 · 유효 퍼텐셜 · 플로케 이론
같은 원리	파울 트랩 · 폰데로모티브 힘

카피차 진자(Kapitza pendulum)는 지지점을 연직 방향으로 충분히 빠르고 작게 흔들면 원래 불안정하던 거꾸로 선 평형점이 안정해지는 강체 진자로, 고주파 진동이 느린 운동에 남기는 평균 힘 — 즉 동역학적 안정화(dynamic stabilization)의 교과서 표본이다. 표트르 카피차가 1951년에 유효 퍼텐셜로 설명하고 실물 장치를 만들어 보이면서 이름이 붙었지만, 안정화 자체는 1908년 스티븐슨이 이미 수학적으로 지적했다.¹

직관을 정면으로 거스른다는 것이 이 계의 매력이다. 흔들기의 시간 평균 힘은 정확히 0이다. 위로 밀린 만큼 아래로 밀린다. 그런데도 진자는 선다. 답은 평균이 0인 것과 효과가 0인 것은 다르다는 데 있다 — 흔들리는 강제력이 자기가 만든 진동 응답과 상관되면서, 평균이 0이 아닌 잔여력을 남긴다. 이 잔여력이 만드는 것이 유효 퍼텐셜(effective potential)이고, 거기서 위쪽이 골짜기가 된다.

이 문서는 현상 자체와 그 정량적 조건, 그리고 근사가 어디서 깨지는지를 다룬다. 평균화 정리의 일반론과 2차 평균의 형식적 기계는 이중 평균법이, 선형화 후의 안정선도·플로케 판정은 마티외 방정식이 소스이므로 그쪽으로 넘긴다.

2. 운동방정식

길이 $l$, 질량 $m$의 진자를 생각한다. 지지점의 연직 위치를 $y_p(t) = A\cos\omega t$로 흔든다. 지지점에 붙은 (비관성) 좌표계에서 보면 관성력이 중력에 그대로 더해지므로, 유효 중력가속도가

$$g_{\text{eff}}(t) = g + \ddot y_p = g - A\omega^2\cos\omega t$$

로 시간에 따라 요동친다. 아래쪽 연직에서 잰 각을 $\theta$라 하면

$$\ddot\theta = -\frac{1}{l}\bigl(g - A\omega^2\cos\omega t\bigr)\sin\theta$$

가 전부다. 여기서 눈여겨볼 것이 강제 진폭에 $\omega^2$이 붙는다는 점이다. $A$가 아무리 작아도 $\omega$를 올리면 순간 가속도는 얼마든지 커진다. “작게, 그러나 빠르게”라는 처방이 모순이 아닌 이유다.

지배적인 무차원수는 두 개다.

$$\varepsilon = \frac{A}{l} \ll 1, \qquad \frac{\omega}{\omega_0} \gg 1, \quad \omega_0=\sqrt{g/l}$$

$\varepsilon$은 흔들기 진폭이 진자보다 훨씬 작다는 것, 두 번째는 흔들기가 진자의 고유진동보다 훨씬 빠르다는 것이다. 뒤에 보겠지만 두 조건이 다 필요하다 — 하나만 만족하면 안정화가 안 되거나, 됐다가 다시 깨진다.

3. 빠른 변수와 느린 변수

각을 두 조각으로 쪼갠다. 느린 평균 운동 $\Theta$와, 흔들기 주기에 맞춰 떠는 작은 진동 $\xi$.

$$\theta(t) = \Theta(t) + \xi(t), \qquad |\xi|\ll1, \quad \langle\xi\rangle = 0$$

$\langle\cdot\rangle$은 빠른 주기 $2\pi/\omega$에 대한 평균이다.

빠른 조각. $\xi$의 방정식에서 $\ddot\xi$와 균형을 이룰 만큼 큰 항은 $\omega^2$이 붙은 강제항뿐이고, 그 계수 안의 $\Theta$는 이 시간척도에서 상수로 굳어 있다.

$$\ddot\xi \simeq \frac{A\omega^2}{l}\cos(\omega t)\sin\Theta \qquad\Longrightarrow\qquad \xi(t) = -\frac{A}{l}\cos(\omega t)\,\sin\Theta$$

두 번 적분하면서 $\omega^2$이 정확히 상쇄된다. 즉 떨림의 크기는 진동수와 무관하게 $\varepsilon\sin\Theta$ 정도다. $\Theta=0$이나 $\pi$, 즉 연직에 정확히 서 있으면 떨림이 아예 없다는 것도 읽힌다.

느린 조각. 원래 식의 $\sin\theta$를 $\sin\Theta + \xi\cos\Theta$로 전개하고 한 주기 평균한다. $\langle\cos\omega t\rangle = 0$이라 강제항의 1차 성분은 통째로 죽고, 강제와 떨림의 상관항만 살아남는다.

$$\bigl\langle \cos(\omega t)\,\xi(t)\bigr\rangle = -\frac{A}{l}\sin\Theta\,\bigl\langle\cos^2\omega t\bigr\rangle = -\frac{A}{2l}\sin\Theta$$

$\langle\cos^2\rangle = 1/2$라는 초등적인 사실 하나가 이 현상 전체를 떠받친다. 결과는

$$\ddot\Theta = -\frac{g}{l}\sin\Theta - \frac{A^2\omega^2}{2l^2}\sin\Theta\cos\Theta$$

새로 생긴 둘째 항은 $A^2\omega^2$에 비례한다. 강제력의 평균이 아니라 제곱의 평균이 들어왔다는 것이 요점이며, 평균이 0인데 효과가 0이 아닌 이유가 정확히 이것이다. 형식적으로 이것은 2차 평균화가 남기는 항이고, 왜 1차에서는 아무것도 안 나오는지에 대한 정직한 정리가 이중 평균법에 있다.

4. 유효 퍼텐셜과 안정 조건

위 방정식은 다음 퍼텐셜 안에서의 보통 운동이다.

$$U_{\text{eff}}(\Theta) = -mgl\cos\Theta + \frac{mA^2\omega^2}{4}\sin^2\Theta$$

$\ddot\Theta = -\,U_{\text{eff}}'(\Theta)/(ml^2)$을 확인하면 바로 맞는다.² 이 꼴은 하전입자가 비균질 고주파 전자기장에서 받는 폰데로모티브 퍼텐셜 $\langle f^2\rangle/(2m\omega^2)$와 정확히 같은 구조다 — 진동하는 힘이 센 쪽이 언덕이 된다.

첫째 항은 익숙한 중력 퍼텐셜로 $\Theta=0$이 골, $\Theta=\pi$가 마루다. 둘째 항은 $\sin^2$이라 $\Theta = \pi/2$에서 최대이고 연직 방향 양쪽($0,\pi$) 모두에서 0이다. 즉 흔들기는 수평 방향에 언덕을 쌓아 올린다. 그 언덕이 충분히 높으면 원래 마루였던 $\Theta=\pi$가 양옆이 더 높아지면서 상대적으로 골짜기가 된다.

정량적으로는 $\Theta=\pi$에서의 2계 도함수 부호다.

$$U_{\text{eff}}''(\pi) = mgl\bigl(-1\bigr) + \frac{mA^2\omega^2}{2} > 0 \qquad\Longleftrightarrow\qquad \boxed{\;(A\omega)^2 > 2gl\;}$$

$A\omega$는 지지점의 속도 진폭이다. 조건이 진폭이나 진동수 하나가 아니라 그 곱으로 나온다는 것이 실험적으로 중요하다 — 진폭을 반으로 줄이면 진동수를 두 배로 올려야 본전이다. 흔히 인용되는 무차원 형태는 $\varepsilon^2 (\omega/\omega_0)^2 > 2$, 즉

$$\frac{\omega}{\omega_0} > \frac{\sqrt2}{\varepsilon}$$

$\varepsilon = 0.1$이면 고유진동수의 14배 이상으로 흔들어야 한다는 뜻이다. 길이 20 cm 진자($\omega_0\approx7\,$rad/s)에 진폭 2 cm면 대략 100 rad/s, 16 Hz 근방. 탁상에서 재현 가능한 숫자라는 점이 이 데모가 100년 넘게 살아남은 이유다.

조건을 넘으면 거꾸로 선 평형 주위의 미소 진동 각진동수도 바로 나온다.

$$\omega_{\text{eff}}^2 = \frac{A^2\omega^2}{2l^2} - \frac{g}{l}$$

영상에서 거꾸로 선 진자를 손가락으로 툭 밀면 느릿하게 되돌아오는데, 그 느릿한 진동수가 이것이다. 흔들기 진동수와는 완전히 다른 시간척도라 눈으로도 구분된다.

5. 평형점 지형이 통째로 바뀐다

퍼텐셜이 바뀌었으니 평형점 목록 자체가 달라진다. $U_{\text{eff}}'=0$은

$$\sin\Theta\Bigl(mgl + \frac{mA^2\omega^2}{2}\cos\Theta\Bigr) = 0$$

이고, $\Theta = 0,\pi$ 외에

$$\cos\Theta_* = -\frac{2gl}{A^2\omega^2}$$

라는 새 평형점이 등장한다. 우변의 절댓값이 1보다 작아야 하므로 존재 조건이 정확히 $(A\omega)^2>2gl$ — 안정화 조건과 같은 순간이다. 이 $\Theta_*$는 불안정(안장)이고, $\Theta=\pi$의 흡인 유역 경계 역할을 한다.

즉 문턱을 넘는 순간 벌어지는 일은 이렇다. 불안정하던 $\Theta=\pi$가 안정으로 바뀌면서, 그 대가로 불안정 평형 한 쌍이 $\pi$에서 갈라져 나온다. 대칭성을 가진 분기 이론의 갈래질(pitchfork) 그 자체다. 구동을 더 세게 하면 $\cos\Theta_*\to0$, 즉 $\Theta_*\to\pi/2$로 벌어지므로 거꾸로 선 상태의 유역은 문턱 바로 위에서 0이었다가 상반구 전체로 자란다. “안정해졌다”와 “웬만큼 밀어도 버틴다”는 별개 질문이고, 후자는 문턱보다 한참 위에서야 참이 된다.

6. 근사는 어디서 깨지는가

평균화는 공짜가 아니다. $(A\omega)^2>2gl$은 최저차 결과이고, 세 방향에서 교정을 받는다.

위쪽 진폭 한계. 거꾸로 선 위치 근방에서 선형화하면 ($y = \theta-\pi$)

$$\ddot y - \frac{g}{l}y + \frac{A\omega^2}{l}\cos(\omega t)\,y = 0$$

인데, $\tau=\omega t/2$로 재척도하면 마티외 방정식 표준형의

$$a = -\frac{4g}{l\omega^2}, \qquad |q| = \frac{2A}{l} = 2\varepsilon$$

이 된다. 안정 영역은 인스-스트럿 선도의 첫 번째 흰 띠, 즉 $a_0(q)$와 $b_1(q)$ 사이다. 아래쪽 경계 $a > a_0(q) \approx -q^2/2$가 방금의 $(A\omega)^2>2gl$이고 — 위쪽 경계 $a < b_1(q)$가 평균화 이론이 아예 모르는 조건이다. 빠른 구동 극한($a\to0^-$)에서 그 경계는 $q\lesssim0.908$이므로

$$\varepsilon = \frac{A}{l} \lesssim 0.45$$

가 된다. 너무 세게 흔들면 다시 넘어진다. 이때의 불안정은 플로케 승수가 $-1$을 통과하는 주기배가라, 진자가 그냥 쓰러지는 것이 아니라 구동 진동수의 절반으로 좌우로 펄럭이기 시작한다(flutter mode). 실험에서 확인된 현상이며, 이 상한이 존재한다는 사실 자체가 유효 퍼텐셜 그림의 한계를 보여 준다.

아래쪽 경계의 고차 보정. 같은 선도에서 $a_0(q) = -q^2/2 + 7q^4/128 - \cdots$까지 쓰면 조건이

$$(A\omega)^2\Bigl(1 - \frac{7}{16}\varepsilon^2\Bigr) > 2gl$$

로 조금 빡빡해진다. $\varepsilon=0.2$면 보정이 1.8 %, $\varepsilon=0.4$면 7 %다. 정밀 실험에서는 이 차이가 보인다.

감쇠. 실물에는 마찰이 있다. 감쇠 $2\gamma\dot\theta$가 붙으면 $\theta = e^{-\gamma t}u$ 치환으로 무감쇠 마티외로 정확히 되돌아가고(마티외 방정식의 그 항등변환), 안정 영역이 $a$축으로 $\gamma^2$만큼 이동한다. 그런데 두 경계가 서로 다르게 반응한다는 것이 이 계의 재미있는 대목이다.

• 아래쪽 경계 $a_0$는 안 움직인다. $\gamma^2$만큼 밀려 들어간 만큼, $\theta$가 유계이려면 $u$가 $e^{\gamma t}$보다 느리게만 자라면 되므로 허용치도 딱 그만큼 늘어난다. 두 효과가 최저차에서 정확히 상쇄된다. 확인은 $q=0$에서 한 줄로 된다 — $\ddot\theta + 2\gamma\dot\theta + (g_{\!\ast}/l)\theta=0$의 특성근이 $-\gamma\pm\sqrt{\gamma^2 - g_{\!\ast}/l}$인데, 뒤집힌 평형($g_\ast<0$)이면 $\gamma$가 아무리 커도 한쪽 근이 항상 양수다. 마찰로는 넘어지는 것을 막지 못한다는, 생각해 보면 당연한 사실. 따라서 $(A\omega)^2>2gl$은 감쇠가 있어도 그대로 살아남는다.
• 위쪽 경계 $b_1$은 올라간다. 이쪽은 파라메트릭 공진형(주기배가) 불안정이라 성장률이 감쇠를 이겨야 켜진다. 즉 마찰이 있으면 $A/l$ 상한이 0.45보다 높아지고, 안정 띠가 위로 넓어진다.

덤으로 감쇠는 진짜 점근 안정을 만들어 준다. 무감쇠 마티외의 “안정”은 유계 진동일 뿐 수렴이 아니라서, 이상적인 무마찰 카피차 진자는 한 번 밀면 영원히 흔들린다. 실험 영상에서 진자가 결국 똑바로 멎는 것은 마찰 덕이다.

7. 왜 중요한가 — 동역학적 안정화

카피차 진자가 물리학 교과서에 남은 이유는 진자 자체보다 “고주파 흔들기로 불안정한 평형을 안정화한다”는 처방이 도처에서 재사용되기 때문이다.

• 파울 트랩과 이온 트랩. 라플라스 방정식 때문에 정전기장만으로는 하전입자를 3차원으로 가둘 수 없다 — 어느 방향으로든 안장점이 생긴다. 그 안장을 RF로 흔들어 가둔 것이 파울 트랩이고, 그 유효 퍼텐셜이 이 문서의 $U_{\text{eff}}$와 문자 그대로 같은 폰데로모티브 퍼텐셜이다. 볼프강 파울이 1989년 노벨상 강연에서 카피차 진자 데모 장치를 들고 나온 것도 우연이 아니다.
• 플라스마와 레이저 트랩. 비균질 고주파장 속 하전입자에 작용하는 폰데로모티브 힘, 광쌍극자 덫(optical dipole trap)의 갇힘, 자기 미러의 고주파 안정화가 전부 같은 계산이다.
• 유체 계면. 무거운 유체가 가벼운 유체 위에 놓인 레일리-테일러 불안정 배치도 용기를 연직으로 진동시키면 안정화된다. 뒤집힌 진자의 연속체 판이며, 조건도 유사한 꼴이다. 반대로 진폭을 더 올리면 패러데이 파가 터진다 — 안정선도의 위쪽 경계를 넘은 것이다.
• 진동 제어(vibrational control). 되먹임 없이 고주파 입력만으로 계를 안정화하는 제어 기법의 원형이다. 센서가 필요 없다는 것이 장점이고, 원하는 상태를 정확히 겨냥할 수 없다는 것이 단점이다. 되먹임으로 역진자를 세우는 것과는 철학이 정반대다.
• 다관절 확장. 애치슨과 멀린이 1993년에 보인 것이 압권이다. 관절이 여러 개인 사슬 진자도 지지점을 흔들면 통째로 곧게 선다. 인도의 로프 마술처럼 보인다고 해서 별명이 붙었고, 조건은 각 관절의 유효 퍼텐셜이 동시에 골이 되는 것이다.³

8. 수치적으로 다루기

원식은 강성(stiff)이 아니지만 시간척도가 두 개라는 것이 함정을 만든다.

• 스텝 크기. 빠른 주기 $2\pi/\omega$를 최소 수십 스텝으로 나눠야 한다. 이걸 어기면 재앙이 아주 그럴듯한 얼굴로 온다 — 스텝이 빠른 진동을 언더샘플링하면 적분기가 강제항을 저주파 별칭으로 읽고, 문턱 아래에서도 진자가 얌전히 서 있는 가짜 안정화를 보여 준다. “돌려 보니 되던데요”가 가장 위험한 계 중 하나다. 구동 주기당 스텝 수를 두 배로 늘려도 답이 같은지 반드시 확인할 것.
• 적분기. 감쇠 없는 경우는 심플렉틱 적분기(예: 속도 벨레) 쪽이 장시간 에너지 표류가 없어 유리하고, 감쇠가 있으면 룽게-쿠타법 4차로 충분하다.
• 안정성 판정은 적분이 아니라 플로케로. “궤적이 안 넘어지네”는 유한 시간 관찰일 뿐이다. 선형화 방정식을 한 주기 $[0,2\pi/\omega]$ 동안 두 초기조건 $(1,0),(0,1)$으로 적분해 모노드로미 행렬 $M$을 만들고 $|\operatorname{tr}M|$을 재는 것이 정확한 판정이다. 무감쇠면 $\det M=1$이라 $|\operatorname{tr}M|<2$가 곧 안정이다(플로케 이론).
• 검증 루틴. 평균화가 맞는지 보려면 완전 비선형 방정식의 궤적에서 $2\pi/\omega$ 이동 평균을 뽑아 $\Theta(t)$를 얻고, 유효 퍼텐셜 방정식의 해와 겹쳐 그린다. 문턱 근처에서 두 곡선이 갈라지기 시작하는 $\varepsilon$이 곧 근사의 유효 한계다. 미소 진동 진동수 $\omega_{\text{eff}}$를 FFT로 재서 예측식과 맞춰 보는 것도 좋은 정량 검사다.

9. 관련 문서

• 이중 평균법 · 마티외 방정식 · 플로케 이론 · 힐 방정식
• 분기 이론 · 안장-마디 분기 · 흡인 유역 · 정규형 이론
• 파울 트랩 · 패러데이 파 · 레일리-테일러 불안정 · 역진자
• 심플렉틱 적분기 · 룽게-쿠타법 · 고속 푸리에 변환
• 더핑 방정식 · 극한 순환 · 동역학계 · 감쇠

10. Footnotes

1. A. Stephenson, On induced stability, Phil. Mag. 15 (1908) 233. 카피차의 1951년 논문(Dynamic stability of a pendulum when its point of suspension vibrates)이 이름을 가져간 것은 그가 유효 퍼텐셜이라는 설명 도구를 세웠기 때문이다. 스티븐슨은 “이 해가 유계다”를 보였고 카피차는 “왜 그런지”를 그림 한 장으로 말했다. 물리학사에서 이름은 대체로 후자에게 간다. ↩

2. 부호가 헷갈리기 쉬운 자리다. $U_{\text{eff}}$를 $ml^2$로 나눈 것이 각가속도의 퍼텐셜이므로 $\ddot\Theta = -U_{\text{eff}}'/(ml^2)$이고, $\tfrac{d}{d\Theta}\bigl(\tfrac{A^2\omega^2}{4}\sin^2\Theta\bigr) = \tfrac{A^2\omega^2}{4}\cdot2\sin\Theta\cos\Theta$를 $l^2$로 나누면 본문의 $\tfrac{A^2\omega^2}{2l^2}\sin\Theta\cos\Theta$가 나온다. 계수 2와 4가 어디서 왔는지만 놓치지 않으면 된다. ↩

3. D. J. Acheson & T. Mullin, Upside-down pendulums, Nature 366 (1993) 215. 애치슨은 임의의 $N$-링크 진자가 뒤집힌 상태로 안정화될 수 있음을 보였고, 멀린이 실험으로 확인했다. 영상으로 보면 스파게티가 접시 위에 서 있는 것처럼 보이는데, 물리학과 학부 데모 중 “이게 왜 되는지 설명해 보라”는 질문으로 대학원생을 가장 오래 붙잡아 두는 부류에 속한다. ↩