ハミルトニアン：$\hat{H}= -\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+ \frac{1}{2}m\omega^2 x^2$
エネルギー固有値：$E_n = \hbar\omega\!\left(n+\tfrac{1}{2}\right)$
波動関数：$\psi_n(x) = N_n\, H_n\!\left(\frac{x}{x_0}\right) e^{-x^2/2x_0^2}$
特性長：$x_0 = \sqrt{\dfrac{\hbar}{m\omega}}$
不確定性：$\Delta x \cdot \Delta p = \hbar\!\left(n+\tfrac{1}{2}\right)$
量子調和振動子シミュレーターとは
🙋
量子調和振動子って何ですか？普通のバネの振動と何が違うんですか？
🎓
ざっくり言うと、バネにくっついた粒子の「量子力学バージョン」だね。古典的なバネはどんなエネルギーでも取れるけど、量子の世界ではエネルギーが$E_n = \hbar\omega(n+\frac{1}{2})$って感じで、$n=0,1,2...$と離散的な値しか取れないんだ。このシミュレーターで上のスライダーを動かして量子数$n$を変えると、そのエネルギーに対応する「波動関数」の形がどう変わるか、すぐに確認できるよ。
🙋
え、$n=0$でもエネルギーが$\hbar\omega/2$ってあるんですか？ゼロじゃないんですね。それと、グラフで波動関数と一緒に確率密度も出てますけど、これは何を見てるんですか？
🎓
その通り！これが「ゼロ点エネルギー」で、量子力学の大きな特徴なんだ。粒子が完全に静止できないことを意味しているよ。確率密度$|\psi|^2$は、その位置$x$で粒子を見つける確率の高さを表している。例えば$n=0$の基底状態では、確率密度が釣鐘型で中心付近にいる確率が一番高い。でも$n$を大きくしていくと、確率密度の山が端の方に寄っていくのが分かるよね。これが古典的な振動（端で速度ゼロ、中心で速度最大）に近づく「対応原理」の一例だ。
🙋
なるほど！左下に表示されてる「不確定性」って、パラメータの質量$m$や角周波数$\omega$を変えるとどうなるんですか？
🎓
良いところに気が付いたね。不確定性原理$\Delta x \Delta p \geq \hbar/2$は、位置と運動量のばらつきの積には下限があるって原理だ。シミュレーターで質量$m$を大きく（重い粒子に）してみて。波動関数が狭くまとまって$\Delta x$が小さくなるけど、その分$\Delta p$が大きくなるのが分かるかな？逆に角周波数$\omega$を大きく（バネを硬く）すると、ポテンシャルの谷が狭くなるので、粒子の存在範囲$\Delta x$も自然と小さくなるんだ。実際に触って確かめてみよう！
よくある質問
波動関数のグラフで「実部」「虚部」「確率密度」を切り替えられますが、それぞれ何を表していますか？
「実部」と「虚部」は波動関数ψ(x)の複素数としての成分で、量子状態の位相情報を含みます。「確率密度」は|ψ(x)|²で、粒子が位置xに見つかる確率を表します。観測可能な物理量は確率密度です。
エネルギー準位のn（量子数）を大きくすると波動関数の形はどう変わりますか？
量子数nが大きくなるほど、波動関数の節（ゼロ交差点）の数がn個に増え、振幅が大きくなります。また、古典的な調和振動子の運動に近づき、端で確率密度が高くなる「古典的領域」が現れます。
質量mや角周波数ωを変えるとエネルギー準位はどう変化しますか？
エネルギー準位はE_n = ℏω(n+1/2)で与えられるため、ωを大きくすると準位間隔が広がり、mを変えても準位間隔は変わりません（ただし波動関数の空間的広がりはmに依存します）。
このシミュレーターで得られた結果は、現実の分子振動などに応用できますか？
はい。二原子分子の振動や結晶中の原子の格子振動（フォノン）は、低エネルギー領域では量子調和振動子で近似できます。ただし、実際の分子では非調和性が現れるため、高励起状態では注意が必要です。
実世界での応用
分子の振動解析： 二原子分子などで、原子間の結合をバネとみなす近似が広く使われます。赤外分光法で観測される分子の振動スペクトルは、このモデルに基づいて理解され、量子数$n$の遷移に対応する光子の吸収・放出として現れます。
固体物理学（格子振動）： 結晶中の原子は平衡位置の周りで振動しており、その量子化されたエネルギー量子は「フォノン」と呼ばれます。固体の比熱や熱伝導率などの物性は、この調和振動子の集合としてモデル化されて説明されます。
量子光学： 電磁場の量子論において、各周波数モードの電磁場は調和振動子と等価であることが知られています。レーザー光や量子もつれ状態の理解に、この基礎モデルが不可欠です。
CAE・分子動力学シミュレーション： 複雑な分子や材料のシミュレーションでは、原子間ポテンシャルの最も単純な近似として調和振動子ポテンシャルがよく用いられます。特に平衡位置付近の小さな振動を扱う際の出発点となります。
よくある誤解と注意点
まず、波動関数 $\psi(x)$ そのものは「確率の振幅」であり、直接観測できる物理量ではない点に注意だ。観測されるのは確率密度 $|\psi(x)|^2$ の方だよ。グラフ上で $\psi(x)$ が負の値を取っている部分を見て「確率が負？」と混乱する人が多いけど、それは振幅だから問題ない。確率密度は常に正またはゼロだ。
次に、パラメータ設定のコツ。質量 $m$ と角周波数 $\omega$ を同時に大きくすると、特性長 $x_0 = \sqrt{\hbar/(m\omega)}$ が急速に小さくなり、グラフ上で波動関数が非常に狭い範囲に縮んで見えなくなることがある。例えば、$m$ を陽子質量程度の $1.67 \times 10^{-27}$ kg、$\omega$ を $1 \times 10^{16}$ rad/s とすると、$x_0$ は約 $6 \times 10^{-13}$ m（フェムトメートルオーダー）となり、デフォルトの表示範囲では点にしか見えない。こういう時は、表示する $x$ の範囲を狭める（例えば -$2x_0$ から $2x_0$ まで）と形が確認しやすくなる。
もう一つ、実務的な落とし穴として「基底状態（n=0）は古典的には静止状態と誤解されがち」という点がある。確かに古典的な最低エネルギーは原点で静止している状態だが、量子の基底状態にはゼロ点エネルギーがあり、粒子は常に「ゆらぎ」を持っている。このゆらぎ（不確定性）を無視して古典的な設計をすると、極低温環境でのナノデバイスの挙動予測を大きく外す可能性があるんだ。