プリセット
安定捕捉：(a, q) は第1安定領域内 — 粒子は中心付近に閉じ込められます
リアルタイム捕捉アニメーション
永年軌道
マイクロモーション
速度 v
力 F
等電位（瞬時）
理論・主要公式
ポールトラップの電位：
$$\phi(x,y,t)=\frac{[U-V_0\cos\Omega t]\,(x^2-y^2)}{2r_0^2}$$
各軸の運動はマシュー方程式に帰着する：
$$\frac{d^2u}{d\tau^2}+(a-2q\cos 2\tau)\,u=0,\quad \tau=\tfrac{\Omega t}{2}$$
無次元パラメータ：$a=\dfrac{4(q/m)U}{\Omega^2 r_0^2}$、$q=\dfrac{2(q/m)V_0}{\Omega^2 r_0^2}$。
$(a,q)$ が第1安定領域（$a=0$ で $q\lesssim 0.908$）の内側にあるとき軌道は有界＝捕捉。永年振動の角周波数は $\omega_{sec}\approx\tfrac{\beta}{2}\Omega$、$\beta\approx\sqrt{a+q^2/2}$（Dehmelt近似）。
検証：数値積分で $a=0$・$q=0.40$ は60周期以上有界、$q=0.95$ で発散（理論境界 $q\approx0.908$ と一致）。
荷電粒子トラップ（ポールトラップ）とは
🙋
ポールトラップって、どうやってイオンを宙に浮かせるんですか？普通の電場だと、プラス電極にマイナスイオンが吸い寄せられてしまう気がするんですが。
🎓
いいところに気が付いたね。大まかに言うと、電極のプラスとマイナスを非常に速く切り替えるんだ。そうすると、イオンは一箇所に落ち着く暇がなくて、中心付近で振動しながら捕らえられるんだよ。このシミュレーターの上の「AC電圧」と「周波数」のスライダーが、その切り替えの強さと速さを決めるパラメータだ。動かしてみると、粒子の動きが大きく変わるのがわかるよ。
🙋
え、そうなんですか！でも、高速で切り替わる電場の中、イオンがどう動くか想像つきません…。何か決まった動き方があるんですか？
🎓
その「決まった動き方」を表すのが「マシュー方程式」というんだ。これは、パラメータの組み合わせで、粒子の運動が安定（トラップされる）か不安定（吹っ飛ぶ）かが決まる、というルールを教えてくれる。シミュレーターの画面に描かれている縞模様の「安定領域」がその答えだ。パラメータを変えて、点（a, q）がその縞模様の中に入っているか外れているか、確かめてみて。
🙋
「DC成分」のスライダーは何をするんですか？ACだけで十分な気がするのですが…。
🎓
実はACだけだと、ある方向にだけ粒子を押さえつけることになるんだ。そこで、静的なDC電圧を重ねることで、別の方向の運動も安定化させている。例えば、自動車のサスペンションみたいに、バネ（DC成分）とダンパー（AC成分）を組み合わせて車体を安定させるイメージだね。DC成分をゼロにしたり大きくしたりして、粒子の軌道がどう変わるか観察してみよう。
よくある質問
まずAC電圧（q値）を下げてみてください。高すぎると不安定になります。次にDC電圧（a値）を0付近に設定し、周波数を調整してマシュー図の安定領域（a-q平面の三角形領域）に入るよう試行してください。
粒子の軌跡が発散せず、トラップ中心付近に留まり続ける状態が安定領域です。パラメータを変えながら、粒子が長時間（例：100周期以上）捕捉されるかどうかを観察してください。マシュー図の理論曲線を参考にすると効率的です。
実験前にAC電圧・周波数・DC成分の最適な組み合わせを探索するのに役立ちます。特に、マシュー方程式の安定領域を視覚的に確認できるため、実際の装置で粒子を捕捉するためのパラメータ調整の効率が向上します。
aはDC電圧による静的な閉じ込めの強さ、qはAC電圧による動的な閉じ込めの強さを表します。両者の組み合わせで粒子の軌道安定性が決まり、マシュー方程式の解が有界になる領域（安定領域）でのみ捕捉が可能です。
実世界での応用
質量分析計：ポールトラップは「イオントラップ質量分析計」の心臓部です。トラップに捕捉したイオンに特定の周波数の電場を加えると、質量に応じて共鳴して動きが大きくなり、検出器にぶつかります。この共鳴周波数から質量を特定するため、化学分析や薬物検査で広く使われています。
量子コンピュータ：捕捉した単一イオン（例えばYb+）を量子ビットとして利用します。レーザーでイオンの内部状態（スピン）を精密に制御し、量子情報処理を行います。極めて安定なトラップ環境が量子状態の長寿命化に不可欠です。
超高精度時計（イオン時計）：トラップされたイオンの超微細遷移周波数を「時計の振り子」として用います。宇宙航行やGPS衛星に搭載される原子時計よりもさらに高精度が求められる基礎物理実験の分野で活躍しています。
プラズマ物理・核融合研究：高温プラズマ中の荷電粒子を模擬し、その閉じ込めや加熱の基礎現象を研究するツールとしても用いられます。大型装置を建設する前に、小型トラップで原理検証を行うことがあります。
よくある誤解と注意点
まず、「AC電圧を上げれば上げるほど、より強く捕捉できる」と思いがちですが、これは誤解です。確かにAC電圧（qパラメータ）を上げると、粒子を中心に引き戻す力は強まりますが、同時に振動も激しくなります。ある閾値を超えると、粒子の振動が発散してあっという間に脱出してしまいます。例えば、DC成分(a)が0の状態で、qを0.9以上にすると、多くの場合で粒子は不安定になります。パラメータ調整は「力加減」が全てで、安定領域の「中」を狙うことが鉄則です。
次に、シミュレーションは理想環境であることを忘れないでください。実際の実験では、真空度が悪いと残留ガスとの衝突でイオンがエネルギーを失い、捕捉できなくなります。また、電極の形状のわずかな歪みや電源のノイズも無視できません。このツールで「理論上は捕捉できるはずなのに、実機ではうまくいかない」というギャップを想像しながら使うと、実務感覚が養えます。
最後に、「DC成分はオマケ」と考えないでください。DC電圧(a)は、x方向とy方向の捕捉の「バランス」を取る重要な役割があります。例えば、aを0から少し正の値（例えば0.1）にすると、安定領域の形が変わり、捕捉可能なqの範囲が狭まることがわかります。これは、特定の質量のイオンだけを選り分ける「質量選択性」の原理そのものです。DCとACは車の両輪だと心得ましょう。
具体的な計算例
本シミュレーターはスケール化した教育用モデルで、実電圧[V]や周波数[Hz]ではなく無次元パラメータを扱います。例として「安定捕捉」プリセット(V₀=45、Ω=4.0 rad/s、U=0、q/m=1.0)を選ぶと、粒子はトラップ中心付近で有界に振動し続け、ステータスカードの「捕捉中」が維持されます。ここから V₀ を 90 まで上げ Ω を 1.5 rad/s に下げると(「不安定」プリセット相当)、マシュー安定図上の動作点が安定領域から外れ、粒子が次々と「脱出済み」に変わります。DC電圧 U を加えると x・y 方向の閉じ込めバランスが変わり安定領域の形が変化します。表示される平均KEや捕捉数は相対的な傾向を見るためのもので、絶対値は実機の eV やイオン種を意味しません。