院試hub

量子調和振動子生成消滅演算子・エネルギー準位・エルミート多項式

量子調和振動子 H^=p^2/2m+mω2x^2/2\hat H = \hat p^2/2m + m\omega^2\hat x^2/2 のエネルギー固有値は En=ℏω(n+12)E_n = \hbar\omega(n + \frac{1}{2})(n=0,1,2,…n = 0, 1, 2, \dots)で、準位は等間隔に並び、最も低い状態にも零点エネルギー ℏω/2\hbar\omega/2 が残る。固有値は、消滅演算子 a^\hat a と生成演算子 a^†\hat a^\dagger の交換関係 [a^,a^†]=1[\hat a, \hat a^\dagger] = 1 だけから求まる。固有関数はガウス関数にエルミート多項式を掛けた形で、基底状態は ψ0∝e−mωx2/2ℏ\psi_0 \propto e^{-m\omega x^2/2\hbar} である。

量子力学 の第2単元(全4単元)。先に読む単元:シュレディンガー方程式と井戸型ポテンシャル

量子力学を出題する院試の解答・解説PDF

大学・研究科ごとに、過去問の解答と解説を年度別にPDFにしている。収録年度・見本・価格を各商品ページで確認できる。

科目から解答PDFを探す

ハミルトニアンと無次元化

定義

量子調和振動子

質量 mm、角振動数 ω\omega の1次元調和振動子のハミルトニアンは

H^=p^22m+12mω2x^2,[x^,p^]=iℏ\hat H = \frac{\hat p^2}{2m} + \frac{1}{2}m\omega^2\hat x^2, \qquad [\hat x, \hat p] = i\hbar

である。長さの単位 x0=ℏ/(mω)x_0 = \sqrt{\hbar/(m\omega)} を使って ξ=x/x0\xi = x/x_0 とおくと、時間に依存しないシュレディンガー方程式は

d2ψdξ2+(ε−ξ2)ψ=0,ε=2Eℏω\frac{d^2\psi}{d\xi^2} + (\varepsilon - \xi^2)\psi = 0, \qquad \varepsilon = \frac{2E}{\hbar\omega}

になる。古典的な振動は振動の単元で扱う。

生成消滅演算子による固有値の求め方

定義

生成演算子と消滅演算子

a^=mω2ℏ(x^+imωp^),a^†=mω2ℏ(x^−imωp^)\hat a = \sqrt{\frac{m\omega}{2\hbar}}\left(\hat x + \frac{i}{m\omega}\hat p\right), \qquad \hat a^\dagger = \sqrt{\frac{m\omega}{2\hbar}}\left(\hat x - \frac{i}{m\omega}\hat p\right)

と定めると、[x^,p^]=iℏ[\hat x, \hat p] = i\hbar から [a^,a^†]=1[\hat a, \hat a^\dagger] = 1、H^=ℏω(a^†a^+12)\hat H = \hbar\omega\left(\hat a^\dagger\hat a + \frac{1}{2}\right) が出る。N^=a^†a^\hat N = \hat a^\dagger\hat a を数演算子という。

定理・法則

調和振動子のエネルギー固有値

条件

[a^,a^†]=1[\hat a, \hat a^\dagger] = 1 が成り立ち、状態のノルムが負にならない(⟨ψ∣a^†a^∣ψ⟩=∥a^ψ∥2≥0\langle\psi|\hat a^\dagger\hat a|\psi\rangle = \|\hat a\psi\|^2 \ge 0)。1次元の系を考える。

主張

N^\hat N の固有値は n=0,1,2,…n = 0, 1, 2, \dots で、

En=ℏω(n+12)E_n = \hbar\omega\left(n + \frac{1}{2}\right)

各準位は縮退せず、固有状態は ∣n⟩=(a^†)nn!∣0⟩|n\rangle = \frac{(\hat a^\dagger)^n}{\sqrt{n!}}|0\rangle と書ける。

要点

[N^,a^]=−a^[\hat N, \hat a] = -\hat a から、a^\hat a は N^\hat N の固有値を 11 下げる。固有値は負になれないので、下げる操作は a^∣0⟩=0\hat a|0\rangle = 0 で止まり、固有値は非負の整数に限られる。a^ψ0=0\hat a\psi_0 = 0 は1階の微分方程式で、解は定数倍を除いて1つなので縮退はない。

はしご演算子の公式

はしご演算子の作用
a^∣n⟩=n ∣n−1⟩,a^†∣n⟩=n+1 ∣n+1⟩\hat a|n\rangle = \sqrt{n}\,|n - 1\rangle, \qquad \hat a^\dagger|n\rangle = \sqrt{n + 1}\,|n + 1\rangle
位置と運動量
x^=ℏ2mω(a^+a^†),p^=imℏω2(a^†−a^)\hat x = \sqrt{\frac{\hbar}{2m\omega}}(\hat a + \hat a^\dagger), \qquad \hat p = i\sqrt{\frac{m\hbar\omega}{2}}(\hat a^\dagger - \hat a)
期待値
⟨n∣x^∣n⟩=0,⟨n∣x^2∣n⟩=ℏmω(n+12)\langle n|\hat x|n\rangle = 0, \qquad \langle n|\hat x^2|n\rangle = \frac{\hbar}{m\omega}\left(n + \frac{1}{2}\right)

運動エネルギーと位置エネルギーの期待値はどちらも En/2E_n/2(ビリアル定理)。

不確定性
Δx Δp=ℏ(n+12)\Delta x\,\Delta p = \hbar\left(n + \frac{1}{2}\right)

基底状態 n=0n = 0 で最小値 ℏ/2\hbar/2 をとる。

固有関数とエルミート多項式

固有関数の公式

基底状態
ψ0(x)=(mωπℏ)1/4exp⁡(−mωx22ℏ)\psi_0(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left(-\frac{m\omega x^2}{2\hbar}\right)

a^ψ0=0\hat a\psi_0 = 0、つまり (x+ℏmωddx)ψ0=0\left(x + \frac{\hbar}{m\omega}\frac{d}{dx}\right)\psi_0 = 0 を解いて規格化したもの。

一般の固有関数
ψn(x)=12nn!(mωπℏ)1/4Hn(ξ) e−ξ2/2,ξ=mωℏ x\psi_n(x) = \frac{1}{\sqrt{2^nn!}}\left(\frac{m\omega}{\pi\hbar}\right)^{1/4}H_n(\xi)\,e^{-\xi^2/2}, \qquad \xi = \sqrt{\frac{m\omega}{\hbar}}\,x
エルミート多項式
Hn(ξ)=(−1)neξ2dndξne−ξ2,H0=1, H1=2ξ, H2=4ξ2−2, H3=8ξ3−12ξH_n(\xi) = (-1)^ne^{\xi^2}\frac{d^n}{d\xi^n}e^{-\xi^2}, \qquad H_0 = 1,\ H_1 = 2\xi,\ H_2 = 4\xi^2 - 2,\ H_3 = 8\xi^3 - 12\xi
偶奇性
ψn(−x)=(−1)nψn(x)\psi_n(-x) = (-1)^n\psi_n(x)

ポテンシャルが偶関数なので、固有関数は偶関数か奇関数になる。ψn\psi_n は nn 個の節をもつ。

級数解の方法では、ψ=h(ξ)e−ξ2/2\psi = h(\xi)e^{-\xi^2/2} とおいて h′′−2ξh′+(ε−1)h=0h'' - 2\xi h' + (\varepsilon - 1)h = 0 をべき級数で解く。級数が無限に続くと hh は eξ2e^{\xi^2} のように増えて規格化できないので、級数が有限で切れる条件 ε=2n+1\varepsilon = 2n + 1 が必要になる。これがエネルギーの量子化 En=ℏω(n+12)E_n = \hbar\omega(n + \frac{1}{2}) を与え、切れた多項式がエルミート多項式 HnH_n である。

はしご演算子と波動関数の計算例

計算例

基底状態の確かめ

ψ0(x)=Ce−αx2\psi_0(x) = Ce^{-\alpha x^2} がシュレディンガー方程式を満たす α\alpha と、そのときのエネルギー、規格化定数 CC を求めよ。

解答

ψ0′′=(4α2x2−2α)ψ0\psi_0'' = (4\alpha^2x^2 - 2\alpha)\psi_0 なので

−ℏ22mψ0′′+12mω2x2ψ0=[(12mω2−2ℏ2α2m)x2+ℏ2αm]ψ0-\frac{\hbar^2}{2m}\psi_0'' + \frac{1}{2}m\omega^2x^2\psi_0 = \left[\left(\frac{1}{2}m\omega^2 - \frac{2\hbar^2\alpha^2}{m}\right)x^2 + \frac{\hbar^2\alpha}{m}\right]\psi_0

右辺が Eψ0E\psi_0 になるには x2x^2 の係数が 00、つまり α=mω/2ℏ\alpha = m\omega/2\hbar であればよく、E0=ℏ2α/m=ℏω/2E_0 = \hbar^2\alpha/m = \hbar\omega/2。規格化は C2∫e−2αx2dx=C2π/(2α)=1C^2\int e^{-2\alpha x^2}dx = C^2\sqrt{\pi/(2\alpha)} = 1 から C=(mω/πℏ)1/4C = (m\omega/\pi\hbar)^{1/4}。

検算

E0=ℏω/2E_0 = \hbar\omega/2 は EnE_n の式の n=0n = 0 と一致する。⟨x2⟩=1/(4α)=ℏ/(2mω)\langle x^2\rangle = 1/(4\alpha) = \hbar/(2m\omega) で、はしご演算子で求めた期待値の n=0n = 0 とも一致する。

計算例

第1励起状態の波動関数

ψ1=a^†ψ0\psi_1 = \hat a^\dagger\psi_0 を計算し、ψ0\psi_0 と直交することを確かめよ。

解答

座標表示では a^†=12(ξ−ddξ)\hat a^\dagger = \frac{1}{\sqrt{2}}\left(\xi - \frac{d}{d\xi}\right)。e−ξ2/2e^{-\xi^2/2} に作用させると 12(ξ+ξ)e−ξ2/2=2 ξe−ξ2/2\frac{1}{\sqrt{2}}(\xi + \xi)e^{-\xi^2/2} = \sqrt{2}\,\xi e^{-\xi^2/2} なので

ψ1(x)=(mωπℏ)1/42 ξ e−ξ2/2\psi_1(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\sqrt{2}\,\xi\,e^{-\xi^2/2}

ψ0ψ1\psi_0\psi_1 は奇関数なので ∫ψ0ψ1 dx=0\int\psi_0\psi_1\,dx = 0。

検算

H1(ξ)=2ξH_1(\xi) = 2\xi を一般の式に入れると 12⋅2ξ=2 ξ\frac{1}{\sqrt{2}} \cdot 2\xi = \sqrt{2}\,\xi で一致する。a^†\hat a^\dagger の係数 0+1=1\sqrt{0 + 1} = 1 から、ψ0\psi_0 が規格化されていれば ψ1\psi_1 も規格化されている。

計算例

重ね合わせ状態の位置の期待値

時刻 00 に Ψ=(ψ0+ψ1)/2\Psi = (\psi_0 + \psi_1)/\sqrt{2} だった振動子について、⟨x⟩(t)\langle x\rangle(t) を求めよ。

解答

Ψ(t)=(ψ0e−iωt/2+ψ1e−3iωt/2)/2\Psi(t) = (\psi_0e^{-i\omega t/2} + \psi_1e^{-3i\omega t/2})/\sqrt{2}。x^\hat x の行列要素は ⟨0∣x^∣1⟩=⟨1∣x^∣0⟩=ℏ/2mω\langle 0|\hat x|1\rangle = \langle 1|\hat x|0\rangle = \sqrt{\hbar/2m\omega} で、対角要素は 00 なので

⟨x⟩(t)=12ℏ2mω(e−iωt+eiωt)=ℏ2mωcos⁡ωt\langle x\rangle(t) = \frac{1}{2}\sqrt{\frac{\hbar}{2m\omega}}\left(e^{-i\omega t} + e^{i\omega t}\right) = \sqrt{\frac{\hbar}{2m\omega}}\cos\omega t

検算

⟨x⟩\langle x\rangle は古典的な振動子と同じ角振動数 ω\omega で振動し、エーレンフェストの定理 d2dt2⟨x⟩=−ω2⟨x⟩\frac{d^2}{dt^2}\langle x\rangle = -\omega^2\langle x\rangle を満たす。t=0t = 0 の値 ℏ/2mω\sqrt{\hbar/2m\omega} は ∫x ψ0ψ1 dx\int x\,\psi_0\psi_1\,dx を直接計算しても得られる。

調和振動子でよくある誤り

  1. 1

    零点エネルギーを落として En=nℏωE_n = n\hbar\omega とする。基底状態でも E0=ℏω/2E_0 = \hbar\omega/2 が残り、不確定性関係 Δx Δp≥ℏ/2\Delta x\,\Delta p \ge \hbar/2 からエネルギーを 00 にはできない。

  2. 2

    a^†∣n⟩=n ∣n+1⟩\hat a^\dagger|n\rangle = \sqrt{n}\,|n + 1\rangle と係数を取り違える。上げるときは n+1\sqrt{n + 1}、下げるときは n\sqrt{n} で、a^†a^∣n⟩=n∣n⟩\hat a^\dagger\hat a|n\rangle = n|n\rangle になるかで確かめる。

  3. 3

    x^\hat x を a^\hat a、a^†\hat a^\dagger で書くときの係数 ℏ/2mω\sqrt{\hbar/2m\omega} の 22 を落とす。⟨0∣x^2∣0⟩=ℏ/(2mω)\langle 0|\hat x^2|0\rangle = \hbar/(2m\omega) になるかで確かめる。

  4. 4

    規格化定数の 1/2nn!1/\sqrt{2^nn!} を忘れる。∫Hn2e−ξ2dξ=2nn!π\int H_n^2e^{-\xi^2}d\xi = 2^nn!\sqrt{\pi} なので、この分を割る必要がある。

  5. 5

    ψn\psi_n の偶奇を逆に覚える。基底状態 n=0n = 0 は偶関数のガウス関数で、ψn(−x)=(−1)nψn(x)\psi_n(-x) = (-1)^n\psi_n(x) である。

量子調和振動子が出た院試の過去問

院試hubが年度別解説の大問を分類した集計から、この単元の内容を扱う出題を挙げる。問題文は各大学の公式過去問で確かめ、解き方は各年度の解説ページで読める。

量子力学の前後の単元

このページの本文を数式のLaTeXのまま読むには Markdown版。最終更新 2026-11-16。誤りに気づいたら 運営者 へ連絡してほしい。

量子力学を出題する院試の解答・解説PDF

大学・研究科ごとに、過去問の解答と解説を年度別にPDFにしている。収録年度・見本・価格を各商品ページで確認できる。

科目から解答PDFを探す