파울 트랩

물리 전자기학 계산물리 마지막 수정: 2026-09-02 04:39:08

1. 개요[편집]

파울 트랩 Paul Trap
고안	볼프강 파울 (1950년대, 본 대학)
노벨 물리학상	1989 (파울 · 데멜트 · 램지)
가두는 원리	RF 사중극장의 시간 평균 유효 퍼텐셜
지배 방정식	마티외 방정식
제어 변수	안정 파라미터 $(a,q)$
운동 구조	느린 영년 운동 + 구동 주파수 마이크로모션
대표 응용	사중극 질량분석기 · 이온 트랩 양자컴퓨터 · 광시계

가둘 수 없다면 흔들어라. 안장 위의 구슬을 떨어뜨리지 않는 방법은 안장을 돌리는 것뿐이다.

파울 트랩(Paul trap)은 시간에 따라 진동하는 무선주파수(RF) 사중극 전기장으로 하전입자를 3차원 공간에 가두는 장치다. 정전기장만으로는 원리적으로 불가능한 일을 시간 의존성 하나로 뚫었고, 그 대가로 갇힌 이온은 느린 영년 운동과 구동 주파수로 떨리는 마이크로모션이라는 두 겹의 운동을 하게 된다.

볼프강 파울이 1950년대에 사중극 질량분석기를 개발하며 세운 원리이고, 한스 데멜트의 페닝 트랩과 함께 1989년 노벨 물리학상을 받았다. 오늘날 이 장치는 질량분석 도구를 넘어 양자컴퓨터의 대표 플랫폼이자 가장 정확한 시계의 심장이 되었다.

동역학의 수학은 이미 심위키 안에 있다. 운동방정식이 마티외 방정식으로 정확히 환원되고 안정 영역이 어떻게 생기는지는 그쪽이, 빠른 진동이 남기는 유효 퍼텐셜이라는 개념 자체는 카피차 진자가 담당한다. 이 문서는 왜 이 장치가 필요했고, 그 수학이 실물 트랩에서 어떤 공학 문제로 나타나며, 이온 사슬과 시계와 양자 게이트가 어디서 갈라지는지를 다룬다.

2. 출발점 — 언쇼 정리[편집]

전하가 정전기장 안에서 안정하게 갇히려면 퍼텐셜 에너지 $U = e\Phi$ 에 국소 극소점이 있어야 한다. 그런데 전하가 없는 영역에서 퍼텐셜은 라플라스 방정식을 만족한다.

\nabla^2\Phi = \frac{\partial^2\Phi}{\partial x^2}+\frac{\partial^2\Phi}{\partial y^2}+\frac{\partial^2\Phi}{\partial z^2} = 0

세 축 방향 곡률의 합이 0이라는 뜻이다. 셋이 전부 양수일 수는 없으므로 극소점이 존재할 수 없고, 임계점은 언제나 안장점이다. 이것이 언쇼 정리(1842)이며, 정전기장만으로 하전입자를 가두는 장치는 존재하지 않는다는 결론이다. 같은 이유로 영구자석만으로 정지한 자석을 공중에 띄울 수도 없다.

빠져나가는 길은 세 갈래뿐이다.

시간에 따라 장을 흔든다 — 파울 트랩. 안장의 방향을 주기적으로 뒤집으면, 각 순간에는 한 방향으로 밀리지만 시간 평균이 가두는 쪽으로 남는다.
정자기장을 보탠다 — 페닝 트랩. 축 방향은 정전 사중극으로 가두고, 반경 방향의 발산은 강한 축 방향 자기장이 로런츠 힘으로 감아 돌린다.
유도된 쌍극자나 반자성을 쓴다 — 언쇼 정리의 전제(고정된 점전하) 자체를 벗어나는 광집게·반자성 부상 등.

3. RF 사중극과 유효 퍼텐셜[편집]

2차원 사중극 전극에 직류 $U$ 와 진폭 $V$ ·각진동수 $\Omega$ 의 RF 를 함께 걸면 전극 사이 퍼텐셜은

\Phi(x,y,t) = \big(U - V\cos\Omega t\big)\,\frac{x^2-y^2}{r_0^2}

꼴이다. $x$ 축이 가두는 순간에 $y$ 축은 밀어내고, 반주기 뒤에는 역할이 뒤바뀐다. 어느 순간에도 두 방향을 동시에 가두지 못한다는 것이 언쇼 정리의 현장 확인이다.

그런데 이온이 충분히 무겁거나 $\Omega$ 가 충분히 빠르면 사정이 달라진다. 이온은 구동에 실시간으로 따라가지 못하고 작은 진폭으로 떨기만 하는데, 그 떨림이 자기가 놓인 장의 세기와 상관되면서 평균이 0이 아닌 잔여력을 남긴다. 이 잔여력이 만드는 것이 폰데로모티브 퍼텐셜, 또는 의사퍼텐셜(pseudopotential)이다. 진폭 $E_0(\mathbf r)$ 의 진동 전기장 안에서

\Psi(\mathbf r) = \frac{e^2 E_0(\mathbf r)^2}{4 m \Omega^2}

이 유효 퍼텐셜이 되고, 위 사중극에서 $E_0 = 2V\,r/r_0^2$ 이므로

\Psi(r) = \frac{e^2 V^2}{m\Omega^2 r_0^4}\,r^2 = \tfrac12 m\,\omega_r^2 r^2, \qquad \omega_r = \frac{\sqrt2\,eV}{m\,\Omega\,r_0^2}

라는 진짜 조화 우물이 된다. 시간 평균이 0인 힘이 0이 아닌 효과를 남긴다는 이 구조는 카피차 진자의 거꾸로 선 진자와 문자 그대로 같은 물리이며, 볼프강 파울이 1989년 노벨상 강연에 카피차식 데모 장치를 들고 나온 것도 그 때문이다.

유용한 어림 하나. 아래에서 정의할 $q$ 를 쓰면 우물의 깊이는 전극 반경에서

\Psi(r_0)/e = \frac{qV}{4}

볼트다. 전형적인 $q\approx0.3$ , $V\approx300\ \mathrm{V}$ 면 약 22 eV — 상온 열에너지( $\approx0.026$ eV)의 1000배 가까이 되므로, 한 번 잡힌 이온은 좀처럼 안 나간다. 실험실에서 이온 하나를 몇 달씩 가둬 두는 것이 가능한 이유다.

4. 마티외 방정식과 안정 영역[편집]

의사퍼텐셜은 근사다. 정확한 운동방정식은 $\xi = \Omega t/2$ 로 두면 정확히 마티외 꼴이 된다.

\frac{d^2 u}{d\xi^2} + \big(a_u - 2q_u\cos 2\xi\big)\,u = 0, \qquad a_x = \frac{8eU}{m r_0^2\Omega^2},\quad q_x = \frac{4eV}{m r_0^2\Omega^2}

$y$ 방향은 부호가 뒤집혀 $(a_y,q_y) = (-a_x,-q_x)$ 다. 두 방향이 동시에 안정해야 하므로 실제 작동 영역은 두 안정 영역의 교집합이고, 그래서 $(a,q)$ 평면의 첫 안정 영역이 삼각형 비슷한 모양이 된다. 특성곡선과 혀 구조, 사중극 질량분석기가 이 삼각형의 꼭짓점을 어떻게 써먹는지는 마티외 방정식이 다룬다.

이온 트랩 실무에서 기억할 숫자는 몇 개 없다.

RF 만 걸면( $a=0$ ) 조건은 $q < 0.908$ 하나로 줄어든다. 실제 트랩은 여유를 두고 $q \approx 0.1\sim0.5$ 에서 돌린다.
$|a|\ll1$ , $q\ll1$ 인 깊은 안정 영역에서 영년 진동수는

\omega_u \simeq \frac{\Omega}{2}\sqrt{a_u + \frac{q_u^2}{2}}

$a=0$ 이면 $\omega_r = q\Omega/(2\sqrt2)$ 로, 앞의 의사퍼텐셜 결과와 정확히 일치한다. 의사퍼텐셜 근사는 마티외 전개의 $q$ 2차 항인 셈이다.

$q$ 가 0.4를 넘어가면 이 근사의 오차가 눈에 띄게 커지고, 0.908 근처에서는 아예 무의미해진다.

5. 영년 운동과 마이크로모션[편집]

낮은 차수 해를 적으면 트랩 안 이온의 운동 구조가 한눈에 보인다.

u(t) \simeq u_0\cos(\omega_u t)\Big[1 - \frac{q_u}{2}\cos\Omega t\Big]

영년 운동(secular motion): 느린 진동 $\omega_u \ll \Omega$ . 우리가 “이온이 우물에 갇혀 있다”고 말할 때의 그 운동이고, 레이저 냉각으로 식히는 대상이며, 양자 게이트가 이용하는 자유도다.
마이크로모션(micromotion): 구동 주파수 $\Omega$ 로 떠는 작은 흔들림. 진폭이 영년 진폭에 비례하고 계수가 $q/2$ 다.

여기서 두 종류를 반드시 구분해야 한다.

고유 마이크로모션은 위 식 그대로 영년 진폭에 비례하므로, 이온을 RF 널(장이 0인 점 또는 선)에 정확히 놓고 잘 식히면 원리적으로 0으로 수렴한다.

초과 마이크로모션(excess micromotion)은 다르다. 표류 전하, 전극 가공 오차, 위상 불일치 같은 이유로 이온이 RF 널에서 밀려나 앉으면, 아무리 식혀도 그 자리의 RF 장 때문에 진폭 $\propto E_{\rm dc}/\Omega^2$ 짜리 떨림이 남는다. 이것이 정밀 측정의 주적이다.

초과 마이크로모션이 남기는 것	결과
2차 도플러(시간 팽창) 편이	광시계의 계통 오차 예산에서 가장 큰 항 중 하나
진동하는 AC 슈타르크 편이	전이 주파수가 조건에 따라 흔들림
반송파 세기 감소, 마이크로모션 측대역 출현	게이트 충실도 저하
도플러 냉각 효율 저하	이온이 원하는 온도까지 안 내려감

그래서 이온 트랩 실험의 일상 업무 중 하나가 마이크로모션 보정이다. DC 보정 전극으로 이온을 RF 널로 밀어 넣으면서, RF 위상과 상관된 형광 광자 계수·측대역 세기 비·분해된 측대역 스펙트럼 같은 진단으로 잔여량을 잰다. 광시계 논문의 불확도 표에서 “excess micromotion” 줄이 항상 눈에 띄는 위치에 있는 이유다.¹

6. 선형 파울 트랩과 이온 사슬[편집]

원조 3차원 파울 트랩은 쌍곡면 링 전극과 두 엔드캡으로 되어 있고, RF 널이 점 하나다. 이온을 하나 가두는 데는 완벽하지만 여러 개를 가두면 두 번째부터는 널 밖에 앉게 되어 초과 마이크로모션을 뒤집어쓴다. 그래서 오늘날 양자정보 실험의 표준은 선형 파울 트랩이다.

네 개의 막대 전극이 $x$ – $y$ 평면에서 RF 사중극을 만들고, 축( $z$ ) 방향은 정전 엔드캡으로 가둔다. 결정적인 차이는 RF 널이 점이 아니라 축을 따라 뻗은 선이라는 것 — 이온 여러 개를 그 선 위에 일렬로 세우면 전부 마이크로모션 없이 앉는다.

축 방향 구속이 조화 퍼텐셜 $\tfrac12 m\omega_z^2 z^2$ 이고 이온끼리는 쿨롱 반발을 하므로, 평형 배치는 두 힘의 균형으로 정해진다. 특성 길이

\ell = \left(\frac{e^2}{4\pi\varepsilon_0\, m\,\omega_z^2}\right)^{1/3}

가 자연 단위이고, 위치는 이 단위로 무차원화된 대수 방정식을 풀어 얻는다. 두 가지 성질이 실무를 지배한다.

간격이 균일하지 않다. 사슬 중앙이 촘촘하고 양 끝이 성기다. $N$ 이온 사슬의 최소 간격은 대략 $\ell N^{-0.56}$ 규모로 줄어들어, 사슬이 길어질수록 개별 이온에 레이저를 따로 쏘는 것(개별 주소 지정)이 어려워진다.
집단 진동 모드가 곧 데이터 버스다. 축 방향 정규 모드는 $N$ 개이며, 가장 낮은 것이 사슬 전체가 함께 움직이는 질량중심 모드( $\omega_z$ ), 다음이 좌우로 늘었다 줄었다 하는 호흡 모드( $\sqrt3\,\omega_z$ )다. 이 모드를 매개로 멀리 떨어진 두 이온의 내부 상태를 얽는 것이 이온 트랩 양자 게이트의 원리다.

축 구속을 너무 세게 하거나 반경 구속을 풀면, 어느 임계 비율에서 일자 사슬이 지그재그로 접힌다(선형–지그재그 구조 전이). 이 전이가 한 트랩에 세울 수 있는 이온 수의 실질적 상한을 만들고, 그래서 규모를 키우려면 사슬을 길게 늘이는 대신 트랩을 여러 구역으로 나누고 이온을 구역 사이로 셔틀링하는 아키텍처로 간다.

7. 응용[편집]

질량분석. 원조이자 가장 널리 쓰이는 용도다. $(a,q)$ 가 둘 다 질량에 반비례하므로 $U/V$ 를 고정한 채 전압을 키우면 서로 다른 $m/z$ 가 $(a,q)$ 평면의 한 직선 위에 늘어서고, 그 직선을 안정 영역 꼭짓점에 아슬아슬하게 걸치면 좁은 질량 창만 살아남는다. 원리의 상세는 마티외 방정식에 있다. 3차원 파울 트랩을 이온 저장통으로 쓰고 불안정 경계를 스캔하며 이온을 순서대로 뱉어내는 이온 트랩 질량분석기도 널리 쓰인다.

양자 컴퓨팅. 1995년 치락과 촐러가 선형 트랩의 공통 진동 모드를 매개로 2큐비트 게이트를 만드는 방법을 제안했고, 같은 해 NIST의 와인랜드 그룹이 단일 이온에서 CNOT 을 시연했다. 이후 뫼르머–쇠렌센 게이트처럼 온도에 덜 민감한 방식이 나오면서 2큐비트 게이트 충실도가 99.9% 대에 이르렀다. 이온 큐비트의 장점은 명확하다 — 모든 이온이 자연이 만든 동일 제품이라 소자 편차가 없고, 결맞음 시간이 초 단위이며, 모든 이온쌍이 공통 모드를 통해 연결된다(전연결성). 단점도 명확하다 — 게이트 속도가 μs~수십 μs 규모로 초전도 큐비트보다 느리고, 사슬을 길게 늘일수록 모드가 촘촘해져 개별 모드를 골라 쓰기가 어려워진다. 양자 회로와 양자 오류 정정의 논리 큐비트 실험 상당수가 이 플랫폼에서 나왔다.

광시계. 갇힌 이온 하나를 레이저 냉각으로 운동 바닥상태 근처까지 식히면, 도플러 편이와 충돌 편이가 사라진 이상적인 분광 표본이 된다. Al⁺·Yb⁺·Sr⁺·Ca⁺ 기반 광시계가 $10^{-18}$ 대의 상대 불확도에 도달해 있고, 초의 재정의 후보로 논의된다. 앞서 말한 초과 마이크로모션 보정이 이 불확도 예산의 핵심 항목이다.

기타. 나노입자·먼지 입자를 가두는 거대 파울 트랩은 아예 눈에 보이는 크기의 입자로 마티외 안정도를 강의실에서 시연하는 데 쓰이고, 저온 화학·이온 이동도 측정·분자 이온 분광에도 널리 쓰인다.

8. 페닝 트랩과의 대비[편집]

같은 언쇼 문제를 정반대 방향으로 푼 형제 장치다.

항목	파울 트랩	페닝 트랩
구속 방식	RF 사중극 (시간 의존 전기장)	정전 사중극 + 강한 축 방향 자기장
운동 성분	영년 운동 + 마이크로모션	축 방향 · 수정 사이클로트론 · 자전(magnetron)
약점	RF 가열, 마이크로모션, 2차 도플러	자전 운동이 준안정(에너지 극대라 감쇠하면 밖으로 나감)
장비	자석 불필요, 광학 접근 자유로움	초전도 자석 필요, 접근 제한
대표 용도	이온 QC, 광시계, 사중극 질량분석	전자 $g$ 인자 정밀 측정, FT-ICR 질량분석, 반물질 연구

페닝 트랩에서 세 고유진동수는 $\omega_+\omega_- = \omega_z^2/2$ , $\omega_+ + \omega_- = \omega_c$ 로 얽혀 있고, 실무에서는 자기장의 비균일성에 둔감한 불변량 정리 $\omega_c^2 = \omega_+^2 + \omega_-^2 + \omega_z^2$ 를 써서 사이클로트론 진동수를 추출한다. 자전 운동이 준안정이라는 것이 이 장치의 고약한 점이다 — 냉각을 잘못 걸면 이온이 식으면서 밖으로 나간다.

9. 수치적으로 다룰 때[편집]

파울 트랩은 겉보기에 단순한 상미분방정식이지만 적분할 때 함정이 있다.

시간 척도가 두 개다. 마이크로모션의 $\Omega$ 와 영년 운동의 $\omega \sim q\Omega/(2\sqrt2)$ 가 한 자릿수 이상 벌어진다. 정확한 궤적을 원하면 구동 주기당 수십~수백 스텝을 밟아야 하고, 그러면서 영년 주기의 수천 배를 적분해야 하므로 스텝 수가 순식간에 커진다. 룽게-쿠타법으로 오래 돌리면 에너지가 조용히 드리프트하므로, 장시간 궤적에는 심플렉틱 적분기나 시간 가역 도식을 쓰는 편이 안전하다.
의사퍼텐셜로 갈아타면 싸지지만 잃는 것이 있다. 구동을 평균해 버린 모형에서는 마이크로모션이 아예 존재하지 않고, 따라서 RF 가열도 안 보인다. 이온이 여럿인 구름에서 쿨롱 충돌이 마이크로모션 에너지를 영년 운동으로 옮기는 현상이 RF 가열인데, 이건 구동을 명시적으로 푸는 계산에서만 나온다. 파라미터 탐색은 의사퍼텐셜로, 최종 검증은 완전 적분으로 — 이게 국룰이다.
이온 사이는 쿨롱이라 $O(N^2)$ 다. 사슬 규모( $N\lesssim100$ )에서는 직접 합이 낫고, 수천 개짜리 구름에서는 N체 문제의 트리·다중극 기법이나 PIC 계열로 넘어간다.
검증 루틴. $(a,q)$ 를 첫 안정 영역 경계 근처로 밀어 가며 궤적의 유계성을 확인하면 안정 선도를 코드가 스스로 재현하는지 볼 수 있고, 이것이 트랩 적분기의 표준 회귀 시험이다. 경계 바로 바깥에서 지수적으로 발산하지 않으면 대개 스텝이 너무 크다.²

10. 여담[편집]

이름은 사람 이름이다. Wolfgang Paul, 파울. “폴 트랩”이라고 읽는 문헌도 흔하지만 독일식으로 파울이 맞고, 본인이 “나는 파울리(Pauli)가 아니라 파울”이라며 자기를 파울리의 배타 원리에 대비되는 존재로 농담 삼았다는 일화가 전해진다.³
갇힌 이온 하나에 레이저를 쏘면 산란광이 맨눈에 보인다. 1980년 데멜트 그룹이 찍은 단일 바륨 이온 사진은 “원자 하나를 눈으로 본” 최초 사례로 인용되며, 지금도 이온 트랩 실험실 문화에서 처음 이온을 잡은 날은 기념일 취급을 받는다.
마이크로모션이 골칫거리이기만 한 것은 아니다. 마이크로모션 측대역의 세기 비가 이온의 위치를 나노미터 정밀도로 알려 주는 진단 신호라, 보정 절차 자체가 이 “잡음”을 자로 쓰는 일이다.

11. 관련 문서[편집]

이온 트랩 대학원생의 통과의례가 “이온이 안 식는다”로 며칠 헤매다 결국 보정 전압이 어긋나 있었음을 발견하는 것이다. 표류 전하는 레이저가 전극에 산란되며 만든 광전자에서도 나오기 때문에, 어제 맞춰 놓은 보정값이 오늘 틀리는 일이 예사다. 자동 보정 루틴을 짜는 것이 사실상 필수 과제. ↩
안정 경계 근처는 마티외 특성곡선이 서로 접근하는 곳이라 해가 준주기적으로 아주 느리게 자란다. 짧게 적분하고 “안정하다”고 판정하면 틀린다. 실물 트랩에서도 같은 함정이 있어서, 경계 근처 파라미터로 잡은 이온은 몇 분 뒤에 조용히 사라지곤 한다. ↩
파울리 배타 원리가 “같은 자리에 둘이 못 있는다”면 파울 트랩은 “같은 자리에 붙잡아 둔다”는 대비인데, 정작 파울 트랩 안의 이온들은 쿨롱 반발로 서로를 밀어내며 결정 격자처럼 늘어선다. 농담의 완성도가 높다. ↩