사인-고든 방정식

편집 역사 토론
계산물리 물리 마지막 수정: 2026-08-10 04:33:07

1. 개요[편집]

사인-고든 방정식
Sine-Gordon Equation
표준형$u_{tt} - u_{xx} + \sin u = 0$
대칭로렌츠 불변 — 특수상대론의 부스트가 그대로 대칭
기본 해킹크 / 안티킹크 / 브리더
보존에너지·운동량 + 정수 토폴로지 전하 + 무한개 보존량
실물비틀림 진자 사슬, 조셉슨 접합의 자속양자, 결정 전위
이름클라인-고든의 말장난(Klein → sine)

사인-고든 방정식(sine-Gordon equation)은 선형 클라인-고든 방정식의 질량항 uu를 sin⁡u\sin u로 바꾼 비선형 쌍곡형 편미분방정식으로, 표준형은 다음과 같다.

utt−uxx+sin⁡u=0u_{tt} - u_{xx} + \sin u = 0

이름부터가 농담이다. 클라인-고든(Klein-Gordon)의 Klein을 sine으로 바꿔치기한 말장난이며, 크루스칼이 퍼뜨린 것으로 알려져 있다.1 하지만 방정식 자체는 19세기 미분기하에서 이미 등장했다 — 3차원 공간에 놓인 일정한 음의 가우스 곡률 곡면(의구면, pseudospherical surface)의 가우스-코다치 방정식이 정확히 이 식이고, 뷔클룬트 변환·해의 중첩 공식 같은 도구가 그 시절에 이미 준비돼 있었다.

솔리톤 일반론과 적분가능계 이야기는 각 문서에 있으니 여기서는 반복하지 않는다. 사인-고든만의 고유한 지분은 세 가지다. 로렌츠 불변성, 정수로 양자화된 토폴로지 전하, 그리고 조셉슨 접합의 자속양자라는 실물 대응. KdV 방정식이나 비선형 슈뢰딩거 방정식에는 이 셋 중 어느 것도 없다.

2. 로렌츠 불변성[편집]

utt−uxxu_{tt} - u_{xx}는 1+1차원 달랑베르시안이고, sin⁡u\sin u는 미분을 포함하지 않는다. 따라서 방정식 전체가 로렌츠 변환

x′=γ(x−vt),t′=γ(t−vx),γ=11−v2x' = \gamma(x - vt), \qquad t' = \gamma(t - vx), \qquad \gamma = \frac{1}{\sqrt{1-v^2}}

에 대해 형태가 그대로 유지된다(c=1c=1로 규격화된 단위계). 이 사실이 실용적으로 뜻하는 바는 명확하다. 정지해 하나만 구하면 움직이는 해는 부스트로 공짜로 얻는다. 그리고 그렇게 얻은 해는 로렌츠 수축을 겪는다 — 빠른 킹크일수록 폭이 1/γ1/\gamma로 좁아지고 에너지는 γ\gamma배로 커진다. 상대론적 입자 흉내를 내는 것이 아니라 진짜로 상대론적 운동학을 따르는 것이다.

수치해석 입장에서 이건 곧바로 격자 요구조건이 된다. v=0.99v = 0.99짜리 킹크를 던지면 폭이 1/γ≈0.141/\gamma \approx 0.14로 줄어들므로 Δx\Delta x를 그만큼 줄여야 하고, 안 그러면 뒤에서 이야기할 격자 고착이 즉시 나타난다.

3. 킹크와 토폴로지 전하[편집]

경계에서 sin⁡u=0\sin u = 0, 즉 u→2πnu \to 2\pi n이 되어야 에너지가 유한하다. 진공이 2π2\pi 간격으로 무한히 많다는 뜻이며, 서로 다른 진공을 잇는 해가 킹크(kink)다.

u(x,t)=4arctan⁡ ⁣(exp⁡(±γ(x−vt)))u(x,t) = 4\arctan\!\Big(\exp\big(\pm\gamma(x - vt)\big)\Big)

++가 킹크, −-가 안티킹크다. 정지 킹크의 기울기는 ux=2 sech xu_x = 2\,\mathrm{sech}\,x로 sech\mathrm{sech} 꼴이고, 에너지 밀도 12ut2+12ux2+(1−cos⁡u)\tfrac12u_t^2 + \tfrac12u_x^2 + (1-\cos u)를 적분하면 정지 에너지가 정확히 88, 속도 vv면 E=8γE = 8\gamma가 된다. 딱 상대론적 입자의 E=mγE = m\gamma 꼴이라 킹크의 정지질량은 8이라고 말해도 된다.

여기서 결정적인 양이 토폴로지 전하다.

Q=u(+∞)−u(−∞)2π∈ZQ = \frac{u(+\infty) - u(-\infty)}{2\pi} \in \mathbb{Z}

킹크는 Q=+1Q=+1, 안티킹크는 Q=−1Q=-1, 진공은 00. 이 값은 연속 변형으로 바뀔 수 없다 — 바꾸려면 경계에서 무한한 에너지를 써야 하기 때문이다. 즉 QQ는 뇌터 정리가 주는 대칭 기반 보존량이 아니라 경계조건의 위상적 성질에서 오는 보존량이고, 그래서 방정식에 감쇠나 강제항을 붙여 적분가능성을 완전히 망가뜨려도 여전히 보존된다. 킹크가 그렇게 질긴 이유가 이것이다.2

4. 브리더[편집]

킹크와 안티킹크는 QQ가 +1+1과 −1-1이라 서로 소멸할 수 있지만, 그 대신 속박 상태를 이룰 수도 있다. 그것이 브리더(breather)다.

u(x,t)=4arctan⁡ ⁣[1−ω2ω⋅sin⁡(ωt)cosh⁡ ⁣(1−ω2 x)],0<ω<1u(x,t) = 4\arctan\!\left[\frac{\sqrt{1-\omega^2}}{\omega}\cdot \frac{\sin(\omega t)}{\cosh\!\big(\sqrt{1-\omega^2}\,x\big)}\right], \qquad 0 < \omega < 1

공간에 국소적이면서 내부 진동수 ω\omega로 진동한다. 에너지는

Ebreather=161−ω2E_{\text{breather}} = 16\sqrt{1-\omega^2}

이고, ω→0\omega \to 0이면 16=8+816 = 8 + 8, 즉 정지한 킹크와 안티킹크 한 쌍의 에너지로 수렴한다. ω\omega가 1에 가까울수록 에너지가 작아지므로 그 차이가 곧 속박 에너지다. 진동수 스펙트럼이 연속이라는 점도 중요하다 — 이 계에는 임의로 약하게 묶인 브리더가 존재하며, 그래서 초기 조건에 잡음을 조금만 섞어도 브리더가 잘 만들어진다.

브리더는 사인-고든의 적분가능성이 얼마나 미묘한 것인지 보여주는 지표이기도 하다. ϕ4\phi^4 모형처럼 비슷하게 생긴 비적분 계에서는 브리더가 엄밀한 의미로 존재하지 않고, 대신 아주 오래 사는 준안정 상태(oscillon)로만 나타나 서서히 방사파를 흘리며 붕괴한다.

5. 충돌 — 위상 이동만 남는다[편집]

킹크끼리, 킹크와 안티킹크, 킹크와 브리더 — 어떤 조합으로 충돌시켜도 사인-고든에서는 각자의 속도와 형태가 정확히 복원되고 위치의 어긋남(위상 이동)만 남는다. 킹크끼리는 서로 밀치고 킹크-안티킹크는 서로 당기지만, 그 상호작용이 산란 후 아무 흔적도 남기지 않는다. 방사파를 조금도 흘리지 않는다는 것이 핵심이다.

대조군을 보면 이게 얼마나 특별한지 알 수 있다. ϕ4\phi^4 모형(utt−uxx−u+u3=0u_{tt}-u_{xx} - u + u^3 = 0)에서 킹크와 안티킹크를 충돌시키면 임계 속도 아래에서 서로 붙잡혔다가 튕겨나가는데, 그 결과가 초기 속도에 따라 프랙탈 구조를 이루는 공명 창(resonance window) 으로 나온다 — 킹크 내부 진동 모드로 에너지를 주고받기 때문이다. 사인-고든에는 그런 내부 모드 자체가 없다. “충돌시켜 보면 안다”는 솔리톤의 판별법이 가장 극적으로 드러나는 한 쌍이다.

6. 물리적 실체 — 진자 사슬에서 자속양자까지[편집]

사인-고든이 다른 적분가능 방정식과 구별되는 점은 손으로 만질 수 있는 모형이 여럿 있다는 것이다.

  • 비틀림 진자 사슬: 팽팽한 강선에 진자들을 촘촘히 매달면, 각 진자의 각도 uju_j는 이웃과의 비틀림 결합(→ uxxu_{xx})과 중력 복원력(→ sin⁡u\sin u)을 동시에 받는다. 연속극한이 정확히 사인-고든이고, 킹크는 사슬을 따라 굴러가는 “한 바퀴 꼬임”이다. 대학 데모 장비로 실제로 만들어 굴린다.
  • 프렌켈-콘토로바 모형: 주기적인 기판 퍼텐셜 위에 스프링으로 연결된 원자 사슬. 여기서 킹크는 여분의 원자 하나가 만드는 전위(dislocation)이고, 킹크의 이동이 곧 소성 변형의 미시적 메커니즘이다. 사인-고든은 이 이산 모형의 연속극한이다.
  • 긴 조셉슨 접합: 두 초전도체 사이 위상차 φ(x,t)\varphi(x,t)가 감쇠·구동항이 붙은 사인-고든을 따른다. 여기서 킹크는 접합 안을 달리는 플럭손(fluxon)이고, 그것이 나르는 자속은 정확히 자속양자 Φ0=h/2e≈2.07×10−15 Wb\Phi_0 = h/2e \approx 2.07\times10^{-15}\ \mathrm{Wb} 하나다. 토폴로지 전하 QQ가 정수인 것이 자속이 양자화되어 있다는 실험 사실과 같은 이야기가 된다. 플럭손 왕복을 이용한 발진기는 실제 소자로 만들어졌다.3

이 대응이 왜 중요하냐면, 사인-고든 시뮬레이션의 검증이 “해석해와 맞나”에서 끝나지 않고 측정 가능한 물리량과 맞나로 이어지기 때문이다. 조셉슨 접합 쪽 응용에서는 감쇠 αut\alpha u_t와 바이어스 전류 −Γ-\Gamma가 붙어 적분가능성이 깨지고, 플럭손이 손실과 구동이 균형을 이루는 종단속도로 안착한다.

7. 수치해석[편집]

쌍곡형 2계 방정식이라 KdV 방정식이나 비선형 슈뢰딩거 방정식과 달리 강성이 심하지 않고, CFL 조건도 Δt≤Δx\Delta t \le \Delta x라는 평범한 형태다. 그 대신 이 방정식만의 함정이 있다.

  • 에너지 보존 도식: utt=uxx−sin⁡uu_{tt} = u_{xx} - \sin u에 표준 립프로그를 쓰고 sin⁡un\sin u^n을 명시적으로 넣으면 전체 에너지가 천천히 표류한다. 고전적인 처방은 비선형항을 이산 그래디언트로 쓰는 것이다. sin⁡u ⟶ cos⁡un−cos⁡un+1un+1−un\sin u \ \longrightarrow\ \frac{\cos u^{n} - \cos u^{n+1}}{u^{n+1} - u^{n}} 이렇게 하면 이산 에너지가 정확히 보존되는 암시적 도식이 되고, 대신 매 스텝 스칼라 비선형 방정식을 풀어야 한다. 장시간 브리더 적분처럼 위상이 중요한 계산에서는 이 비용을 낼 가치가 있다. 기하 수치적분 관점에서는 심플렉틱 적분기와 같은 가족이다.
  • 토폴로지 전하는 공짜로 보존된다: 이산 전하 Qh=12π∑j(uj+1−uj)Q_h = \frac{1}{2\pi}\sum_j (u_{j+1}-u_j)는 망원경 합이라 (uN−u0)/2π(u_N - u_0)/2\pi로 줄어든다. 경계값을 고정하는 한 어떤 도식을 쓰든 기계 정밀도로 보존된다는 뜻이다. 그러므로 QQ가 보존됐다는 것은 코드가 정확하다는 증거가 전혀 아니다 — 검증에 쓸 양은 에너지, 킹크 속도, 브리더 주기 쪽이다.
  • 페이얼스-나바로 장벽: 이산화는 연속 방정식의 병진 불변성을 깨뜨린다. 그 결과 킹크의 에너지가 격자에 대한 위치에 따라 미세하게 달라지는 주기적 퍼텐셜(PN 장벽)이 생기고, 폭에 비해 Δx\Delta x가 크면 느린 킹크가 격자에 고착되어 멈춰버린다. 물리적 현상이 아니라 순수한 이산화 인공물이며, 킹크 폭당 격자점 10개 이상을 국룰로 잡는 이유다. 앞서 말한 로렌츠 수축과 겹치면 빠른 킹크에서 특히 위험하다.4
  • 경계: 유한 구간에서 킹크를 오래 돌리려면 방사파가 벽에서 되튀지 않도록 흡수 경계나 충분히 긴 도메인이 필요하다. 되튄 잔물결이 브리더 진동수를 오염시켜 “브리더가 저절로 붕괴한다”는 가짜 결론을 만들기 쉽다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. “sine-Gordon”이 하도 자연스럽게 굳어져서, 정작 이 이름이 농담이라는 사실을 모르고 쓰는 사람이 더 많다. 물리학계 작명 센스의 대표작으로 quark, strangeness, penguin diagram과 함께 자주 소환된다. ↩

  2. 그래서 “토폴로지적으로 보호된다”는 표현을 쓴다. 다만 이 보호는 전하 총합에 대한 것이지 개별 킹크의 형태에 대한 것이 아니다. 감쇠를 걸면 킹크는 느려지고 뭉개지지만 QQ는 여전히 1이다. 보호받는 것과 멀쩡한 것은 다르다. ↩

  3. 킹크 하나 = 자속양자 하나라는 대응은 “토폴로지 전하가 정수”라는 수학적 진술에 실측 가능한 물리적 살을 붙여준다. 사인-고든이 다른 적분가능 방정식보다 응집물질 쪽에서 유난히 자주 인용되는 이유이기도 하다. ↩

  4. PN 장벽에 걸려 멈춘 킹크를 보고 “우리 계에 새로운 고착 전이가 있다”고 발표하는 사고가 실제로 있었다. Δx\Delta x를 절반으로 줄여봤을 때 사라지는 현상은 물리가 아니라 격자다. ↩