tecktom.com

時間発展:ユニタリと Schrödinger 方程式

2026年10月17日9分で読めますAI生成教科書『虎の巻』第1部 第2章 基準:状態ベクトルと演算子の公理
目次
  1. 共通例題②:2準位系の Rabi 振動を基準の言葉で解く
  2. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第1部 第2章 §1.2.4 ★★★ 確率を保つ連続な時間発展はユニタリ群に限られ、生成子がハミルトニアンになる

定義 1.2.8(公理 Q5(時間発展))

測定しない間、時刻 00 の状態 ∣ψ(0)⟩| \psi(0) \rangle は

∣ψ(t)⟩=U^(t)∣ψ(0)⟩,U^(t)=e−iH^t/ℏ(1.2.8)| \psi(t) \rangle =\hat U(t)| \psi(0) \rangle ,\qquad \hat U(t)=\mathrm{e}^{-\mathrm{i}\hat Ht/\hbar}\tag{1.2.8}

に発展する。H^\hat H は自己共役なハミルトニアンである。微分すると Schrödinger 方程式

iℏ∂∂t∣ψ(t)⟩=H^∣ψ(t)⟩(1.2.9)\mathrm{i}\hbar\frac{\partial }{\partial t}| \psi(t) \rangle =\hat H| \psi(t) \rangle\tag{1.2.9}

になる。

この公理の形は、次の3つの要請から定まる。第一に、時間発展は状態から状態への写像で、遷移確率 ∣⟨ϕ∣ψ⟩∣2\lvert \langle \phi | \psi \rangle \rvert ^2 を保つ(測定をしない間に、確実に区別できた2つの状態が区別できなくなったり、その逆が起きたりしない)。Wigner の定理によれば、射線の間の写像で遷移確率を保つものは、ユニタリ演算子かそれに複素共役を合成した反ユニタリ演算子のどちらかで実現される。第二に、時間発展は合成できる:U^(t+s)=U^(t)U^(s)\hat U(t+s)=\hat U(t)\hat U(s)、U^(0)=1^\hat U(0)=\hat1。U^(t)=U^(t/2)2\hat U(t)=\hat U(t/2)^2 と書けるので、反ユニタリでも2回合成するとユニタリになり、U^(t)\hat U(t) はユニタリに限られる。第三に、tt について連続である。この3つを満たすユニタリの1径数群は、Stone の定理によって自己共役な生成子を1つ持ち、U^(t)=e−iH^t/ℏ\hat U(t)=\mathrm{e}^{-\mathrm{i}\hat Ht/\hbar} と書ける1M. H. Stone, Ann. Math. 33, 643 (1932)。下の脚注へ。生成子を ℏ\hbar で割った形に取ったものがハミルトニアンで、それが古典のハミルトニアンに対応することは第2部第3章と 解析力学から様々な理論分野を組み立てる 第4部第2章で示す。

H^\hat H の固有状態 H^∣n⟩=En∣n⟩\hat H| n \rangle =E_n| n \rangle で展開すると ∣ψ(t)⟩=∑ncne−iEnt/ℏ∣n⟩| \psi(t) \rangle =\sum_nc_n\mathrm{e}^{-\mathrm{i} E_nt/\hbar}| n \rangle で、各成分は位相だけが回る。固有状態1つだけの状態は全体の位相しか変わらないので、どの観測量の分布も時間によらない(定常状態)。異なる EnE_n の成分どうしの相対位相 e−i(En−Em)t/ℏ\mathrm{e}^{-\mathrm{i}(E_n-E_m)t/\hbar} が回ると、分布は時間で変わる。

共通例題②:2準位系の Rabi 振動を基準の言葉で解く

2準位原子を、遷移周波数 ω0\omega_0 に近い周波数 ωL\omega_{\mathrm L} のレーザーで駆動する。基底 ∣e⟩=(1,0)T| e \rangle =(1,0)^{\mathsf T}、∣g⟩=(0,1)T| g \rangle =(0,1)^{\mathsf T}(σz∣e⟩=+∣e⟩\sigma_z| e \rangle =+| e \rangle)を取り、レーザーとともに回る座標に移って回転波近似をすると、ハミルトニアンは時間によらない形

H^=−ℏΔ2σz−ℏΩ2σx(1.2.10)\hat H=-\frac{\hbar\Delta}{2}\sigma_z-\frac{\hbar\Omega}{2}\sigma_x\tag{1.2.10}

になる。Δ=ωL−ω0\Delta=\omega_{\mathrm L}-\omega_0 は離調、Ω\Omega は Rabi 周波数(実数に取る)である。導出は 光学と量子エレクトロニクス 第18章で行った。ここでは (1.2.10) をそのまま受け取り、公理 Q5 だけで ∣ψ(0)⟩=∣g⟩| \psi(0) \rangle =| g \rangle からの時間発展を解く。

一般化 Rabi 周波数 Ω′≔Ω2+Δ2\Omega'\coloneqq\sqrt{\Omega^2+\Delta^2} と単位ベクトル n=(Ω,0,Δ)/Ω′\bm{n}=(\Omega,0,\Delta)/\Omega' を使うと、(1.2.10) はH^=−ℏΩ′2n⋅σ\hat H=-\frac{\hbar\Omega'}{2}\bm{n}\cdot\bm{\sigma} と書ける。例題 1.2.2 の (n⋅σ)2=1^(\bm{n}\cdot\bm{\sigma})^2=\hat1 を使って指数関数を級数で展開すると、偶数次と奇数次がそれぞれcos⁡\cos と sin⁡\sin にまとまり

U^(t)=exp⁡(iΩ′t2n⋅σ)=cos⁡Ω′t2 1^+isin⁡Ω′t2 n⋅σ(1.2.11)\hat U(t)=\exp\Bigl(\mathrm{i}\frac{\Omega't}{2}\bm{n}\cdot\bm{\sigma}\Bigr)=\cos\frac{\Omega't}{2}\,\hat1+\mathrm{i}\sin\frac{\Omega't}{2}\,\bm{n}\cdot\bm{\sigma}\tag{1.2.11}

を得る。∣g⟩| g \rangle から ∣e⟩| e \rangle への振幅は、⟨e∣1^∣g⟩=0\langle e | \hat1 | g \rangle =0、⟨e∣σx∣g⟩=1\langle e | \sigma_x | g \rangle =1、⟨e∣σz∣g⟩=0\langle e | \sigma_z | g \rangle =0 より⟨e∣U^(t)∣g⟩=isin⁡(Ω′t/2) nx=i(Ω/Ω′)sin⁡(Ω′t/2)\langle e | \hat U(t) | g \rangle =\mathrm{i}\sin(\Omega't/2)\,n_x=\mathrm{i}(\Omega/\Omega')\sin(\Omega't/2)。Born の規則 (1.2.4) で2乗して

  Pe(t)=(ΩΩ′)2sin⁡2Ω′t2,Ω′=Ω2+Δ2  (1.2.12)\boxed{\;P_e(t)=\Bigl(\frac{\Omega}{\Omega'}\Bigr)^2\sin^2\frac{\Omega't}{2},\qquad \Omega'=\sqrt{\Omega^2+\Delta^2}\;}\tag{1.2.12}

となる。これが共通例題②の正答(第1部第1章の表)と同じ形の答えである。共鳴(Δ=0\Delta=0)では Pe=sin⁡2(Ωt/2)P_e=\sin^2(\Omega t/2) で、Ωt=π\Omega t=\pi(π\pi パルス)で確実に ∣e⟩| e \rangle へ移る。離調があると振動は速く(Ω′\Omega')、振幅は小さく((Ω/Ω′)2(\Omega/\Omega')^2)なる(図 1.2.2)。ほかの部では、同じ (1.2.12) を第2部第3章が Pauli 行列の Heisenberg 方程式で、第4部第5章が Bloch 方程式で、第6部第6章がスピンのコヒーレント状態の経路積分で出す。

図

図 1.2.2 式 (1.2.12) の Rabi 振動。離調 Δ\Delta を大きくすると周期 2π/Ω′2\pi/\Omega' が短くなり、最大値 Ω2/(Ω2+Δ2)\Omega^2/(\Omega^2+\Delta^2) が下がる。(1.2.10) の中では厳密解なので実線1種類で描いた。

例題 1.2.4 Rabi 振動の数値 ★★★

オリジナルRabi 周波数を Ω/2π=1.00 MHz\Omega/2\pi=1.00\,\mathrm{MHz} とする。

  1. (1) 共鳴(Δ=0\Delta=0)で ∣g⟩| g \rangle から ∣e⟩| e \rangle へ確実に移すのに要る時間(π\pi パルスの長さ)を求めよ。
  2. (2) 離調 Δ=3 Ω\Delta=\sqrt3\,\Omega のとき、PeP_e の最大値と、それに初めて達する時刻を求めよ。
  3. (3) (2) の時刻に ∣g⟩| g \rangle の確率を測る代わりに、U^(t)\hat U(t) がユニタリであることから Pe+Pg=1P_e+P_g=1 を確かめよ。
解答

(1) Ωt=π\Omega t=\pi なので t=π/Ω=π/(2π×1.00×106 s−1)=0.500 μst=\pi/\Omega=\pi/(2\pi\times1.00\times10^6\,\mathrm{s^{-1}})=0.500\,\mu\mathrm{s}。

(2) Ω′=Ω2+3Ω2=2Ω\Omega'=\sqrt{\Omega^2+3\Omega^2}=2\Omega なので最大値は (Ω/Ω′)2=1/4(\Omega/\Omega')^2=1/4。sin⁡2(Ω′t/2)=1\sin^2(\Omega't/2)=1 に初めてなるのは Ω′t=π\Omega't=\pi のときでt=π/(2Ω)=0.500 μs/2=0.250 μst=\pi/(2\Omega)=0.500\,\mu\mathrm{s}/2=0.250\,\mu\mathrm{s}。

(3) (1.2.11) から ⟨g∣U^∣g⟩=cos⁡(Ω′t/2)+isin⁡(Ω′t/2) nz⟨g∣σz∣g⟩=cos⁡(Ω′t/2)−i(Δ/Ω′)sin⁡(Ω′t/2)\langle g | \hat U | g \rangle =\cos(\Omega't/2)+\mathrm{i}\sin(\Omega't/2)\,n_z\langle g | \sigma_z | g \rangle =\cos(\Omega't/2)-\mathrm{i}(\Delta/\Omega')\sin(\Omega't/2)。よって Pg=cos⁡2(Ω′t/2)+(Δ/Ω′)2sin⁡2(Ω′t/2)P_g=\cos^2(\Omega't/2)+(\Delta/\Omega')^2\sin^2(\Omega't/2) で、Pe+Pg=cos⁡2+((Ω2+Δ2)/Ω′2)sin⁡2=1P_e+P_g=\cos^2+\bigl((\Omega^2+\Delta^2)/\Omega'^2\bigr)\sin^2=1。(2) の時刻では Pg=0+(Δ/Ω′)2=34P_g=0+(\Delta/\Omega')^2=\frac34 で、和は 11 である。

近似の裏付け

型 B(極限) | 小さいもの Ω/ω0→0\Omega/\omega_0\to0 かつ ∣Δ∣/ω0→0\lvert \Delta\rvert /\omega_0\to0(反回転の項を時間平均で捨てる) | 証明 条件付き(時間 1/Ω1/\Omega の桁の範囲で誤差 O(Ω/ω0)\mathcal{O}(\Omega/\omega_0)。周期的な係数の常微分方程式に対する平均化の定理による。判定は光学と量子エレクトロニクス 第18章に揃えた) | 行き先 近似の数理 第7章

1 M. H. Stone, Ann. Math. 33, 643 (1932)。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ