自然単位系（atomic units）で計算しています。n を変えるとエネルギー準位と節の数が変わります。
—
古典反転点 x_c = √((2n+1)ℏ/mω)
上：放物線型ポテンシャル V(x)=½mω²x²（青）と離散エネルギー準位 E_0..E_n（黒水平線）／下：ψ_n(x)（緑）と |ψ_n(x)|²（赤）、縦破線＝古典反転点 ±x_c
1次元量子調和振動子のシュレディンガー方程式の解は、Hermite 多項式 $H_n(\xi)$ とガウス関数の積で書けます。
規格化された固有関数（$\xi = x\sqrt{m\omega/\hbar}$）：
$$\psi_n(x) = \frac{1}{\sqrt{2^n n!}}\left(\frac{m\omega}{\pi\hbar}\right)^{1/4} H_n(\xi)\,e^{-\xi^2/2}$$
Hermite 多項式の漸化式：
$$H_{n+1}(\xi) = 2\xi\,H_n(\xi) - 2n\,H_{n-1}(\xi),\quad H_0=1,\;H_1=2\xi$$
エネルギー準位と古典反転点：
$$E_n = \hbar\omega\left(n+\tfrac{1}{2}\right),\qquad x_c = \sqrt{\frac{(2n+1)\hbar}{m\omega}}$$
エネルギーは離散的に等間隔 ℏω で並び、基底状態 (n=0) でも零点エネルギー ℏω/2 が残ります。波動関数は古典反転点を越えて指数的にしみ出します。
量子調和振動子の波動関数シミュレーターとは
🙋
量子調和振動子って、要はバネにつながったボールの量子版ですよね？古典のバネと何が違うんですか？
🎓
ざっくり言うと、エネルギーが連続じゃなくて「飛び飛び」になるのが一番の違いだ。$E_n = \hbar\omega(n+1/2)$ で n=0,1,2,... と等間隔に並ぶ。シミュレーターで量子数 n を増やしてみて。エネルギー準位カードと上のグラフの黒水平線が ℏω ずつ上がっていくのが見える。デフォルト値（n=3, m=ω=ℏ=1）だと E_3 = 3.5 になる。
🙋
n=0 でもエネルギーがゼロにならないんですか？
🎓
そう、それが「零点エネルギー」だ。n=0 でも E_0 = ℏω/2 が残る。古典なら「バネをじっと静止させればエネルギーゼロ」だけど、量子では不確定性原理のせいで位置と運動量を同時にゼロにできない。だから絶対零度でも結晶格子は微振動を続けている。これが熱伝導や比熱の量子補正の出発点になる重要な性質だよ。
🙋
波動関数の緑のカーブを見ると、n=3 で 3 か所で符号が変わってますね。これって偶然ですか？
🎓
偶然じゃなくて、量子数 n の波動関数は必ず「節（ゼロ点）」を n 個持つ。これは Hermite 多項式 $H_n$ の零点数がちょうど n だからだ。n=0 は節ゼロ（純粋なガウス）、n=1 は中央に節 1 個、n=3 は節 3 個。シミュレーターの「節の数」カードがそれを示してる。一般的な「ノードが多いほど高エネルギー」というルールの教科書例だね。
🙋
古典反転点 ±x_c の縦破線の外側にも、波動関数が少しはみ出してますね。
🎓
いいところに気づいた。古典的には粒子は ±x_c の外には絶対行けない（運動エネルギーが負になるから）。でも量子では指数的に減衰しながらしみ出す。これが「トンネル効果」の原型だよ。半導体の量子井戸レーザーや走査トンネル顕微鏡（STM）、α崩壊なんかも、根本はこの「禁止領域へのしみ出し」が効いてる。シミュレーターで n を変えると、しみ出しの広さも変わるのが分かるよ。
よくある質問
なぜエネルギーが等間隔なのですか？
調和振動子のポテンシャル V(x) = ½mω²x² が放物線で、対称性が高いためです。代数的には生成・消滅演算子 a^†, a を使うと a^† が状態を1つ「上げる」働きをし、上げるたびに加わるエネルギーが常に ℏω と一定であることから等間隔が導かれます。水素原子のような 1/r ポテンシャルでは準位間隔が縮んでいくのと対照的です。
高い量子数 n では古典力学に近づきますか？
はい、これは「対応原理」と呼ばれます。n を大きくすると、確率密度 |ψ_n(x)|² の包絡線が古典振動子の存在確率（反転点付近で高く、中央で低い）に近づきます。シミュレーターで n=10 にして |ψ|² の赤いカーブを見ると、激しく振動するけれど包絡線として U 字型になっており、古典の確率分布に対応していることが分かります。
なぜ波動関数が古典反転点の外にも存在するのですか？
量子力学では波動関数は古典的に禁止された領域（E < V の領域）にも指数的に減衰しながら侵入できます。これは Schrödinger 方程式の解として、ポテンシャルが急に無限大にならない限り波動関数が連続でなければならないためです。この「しみ出し」がトンネル効果の起源で、α崩壊や走査トンネル顕微鏡 (STM) の原理になっています。
Hermite 多項式の代わりに数値的に解くことはできますか？
できます。Schrödinger 方程式を有限差分法や Numerov 法で離散化し、行列の固有値問題として解けば数値的にエネルギー準位と波動関数が得られます。調和振動子の場合は解析解（Hermite 多項式）が存在するので比較検証に便利ですが、より一般的なポテンシャル（モース型、非調和、井戸+障壁など）では数値解法が必須になります。
実世界での応用
分子の振動分光： 二原子分子（H₂、CO、HClなど）の伸縮振動は、平衡結合距離まわりのポテンシャルが近似的に放物線で表せるため、調和振動子モデルがそのまま使えます。赤外分光（IR）やラマン分光で観測される吸収・散乱線の位置が ℏω として直接読み取れ、結合の強さの指標になります。より正確には、解離を表すモース型ポテンシャルで補正します。
固体物理のフォノン： 結晶中の原子の格子振動は、隣接原子間のバネ的相互作用が連なった結合振動子系と見なせます。基準振動を取り出すと、それぞれが独立な調和振動子になり、量子化されたエネルギー量子が「フォノン」です。比熱の Einstein モデル・Debye モデル、熱伝導率、超伝導の電子-フォノン相互作用など、固体物理の根幹を支えます。
量子光学とレーザー： 電磁場の各モードも独立な調和振動子として量子化でき、各エネルギー量子が「光子」です。コヒーレント状態 |α⟩ はレーザー光の量子状態を記述し、スクィーズド状態は量子ノイズを低減する精密計測（重力波検出など）に使われます。本シミュレーターで n を変えると、レーザー光の Fock 状態 |n⟩ の波動関数を見ていることになります。
量子コンピュータと量子情報： イオントラップ型量子コンピュータでは、捕捉イオンの調和振動モードが情報伝達のバス（共通モード）として使われます。超伝導量子ビットでも LC 回路の量子化（電磁場の調和振動子化）が基礎で、本モデルの理解が量子ハードウェアの理解に直結します。
よくある誤解と注意点
最も多い誤解は、「波動関数 ψ そのものが物理的に観測できる量」だと考えてしまう ことです。実際に観測できるのは確率密度 |ψ|² であり、波動関数の位相や符号は単独では物理的意味を持ちません（干渉が起きたときに初めて相対位相が観測できる）。シミュレーターでも緑の ψ は符号を変えて振動しますが、赤の |ψ|² は常に非負で、これが「粒子がそこに存在する確率密度」を表します。波動関数の符号反転（ψ → −ψ）はまったく同じ物理状態を表します。
次に多いのが、古典反転点 x_c を「粒子が決して越えない壁」と思い込む 誤解です。古典力学ではそうですが、量子力学では指数的にしみ出す確率があります。シミュレーターで縦破線の外側を見ると、緑と赤のカーブがゼロにはなっておらず、急速に減衰しながら裾を引いていることが分かります。この「しみ出し」こそが量子トンネル効果の原型で、α崩壊や走査トンネル顕微鏡 (STM) や半導体デバイスの根本原理です。「禁止領域」は古典の概念で、量子では確率は小さくともゼロではありません。
最後に、調和振動子モデルが「常に成り立つ近似」だと過信する のは危険です。このモデルは平衡点まわりの小さな振動に対する2次のテイラー展開なので、振幅が大きくなると4次以上の補正（非調和性）が無視できなくなります。実際の分子振動では、量子数が高くなるにつれて準位間隔が縮み、最終的に解離します（モース型ポテンシャルが表現）。シミュレーターは純粋な調和振動子モデルで、現実の系を扱う際は非調和補正の必要性を意識してください。