薛丁格方程、量子井與穿隧：能階量子化與穿牆機制

從波函數到能階：薛丁格方程的物理骨架

一個被困在盒子裡的電子，為什麼只能擁有特定的能量？一道高於粒子動能的位能牆，為什麼粒子卻能「穿牆而過」？這些違反古典直覺的現象，全都源自同一條主控方程——薛丁格方程。本篇假設讀者已熟悉古典力學與基礎微積分，我們直接從非相對論性的量子描述切入。

量子態由波函數 $\psi(\mathbf{r},t)$ 完整描述，其模平方 $|\psi|^2$ 為機率密度。含時薛丁格方程為

$$ i\hbar\,\frac{\partial}{\partial t}\psi(\mathbf{r},t) = \hat{H}\,\psi(\mathbf{r},t), \qquad \hat{H} = -\frac{\hbar^2}{2m}\nabla^2 + V(\mathbf{r}). $$

其中 $\hat{H}$ 為哈密頓算符，第一項對應動能算符 $\hat{p}^2/2m$（透過正則量子化 $\hat{p}\to -i\hbar\nabla$ 得到），第二項為位能。當 $V$ 不顯含時間時，可分離變數 $\psi=\phi(\mathbf{r})e^{-iEt/\hbar}$，得到定態（時間無關）薛丁格方程

$$ -\frac{\hbar^2}{2m}\nabla^2\phi + V(\mathbf{r})\,\phi = E\,\phi. $$

這是一個本徵值問題：能量 $E$ 是哈密頓算符的本徵值，而邊界條件（波函數須單值、連續、平方可積）會把連續的數學解篩成離散的物理譜——這正是「量子化」的數學根源。

無限深量子井：能階量子化的最小模型

考慮一維無限深位能井，$V(x)=0$ 於 $0<x<L$，井外為 $+\infty$。井外波函數必須為零，故邊界條件為 $\phi(0)=\phi(L)=0$。井內方程化為

$$ \frac{d^2\phi}{dx^2} + k^2\phi = 0, \qquad k^2 = \frac{2mE}{\hbar^2}. $$

通解為 $\phi(x)=A\sin(kx)+B\cos(kx)$。由 $\phi(0)=0$ 得 $B=0$；再由 $\phi(L)=0$ 要求 $\sin(kL)=0$，即 $kL=n\pi$（$n=1,2,3,\dots$）。代回 $k$ 與 $E$ 的關係，得能階

$$ E_n = \frac{\hbar^2 k_n^2}{2m} = \frac{n^2\pi^2\hbar^2}{2mL^2}, \qquad n=1,2,3,\dots $$

歸一化條件 $\int_0^L |\phi|^2\,dx = 1$ 給出 $A=\sqrt{2/L}$，故本徵函數為 $\phi_n(x)=\sqrt{2/L}\,\sin(n\pi x/L)$。

值得注意的特徵：能量隨 $n^2$ 增長、與 $L^2$ 成反比（井越窄能階間距越大），且最低能量 $E_1>0$——這個非零的零點能正是不確定性原理 $\Delta x\,\Delta p\gtrsim\hbar/2$ 的直接後果：把粒子限制在尺度 $L$ 內，動量必有不可消去的散布。

有限位能牆與量子穿隧

當位能牆高度有限時，古典上禁止的區域內波函數不再硬性歸零，而是指數衰減——這開啟了穿隧的可能。考慮粒子能量 $E$ 入射一道矩形位能牆：$V(x)=V_0$ 於 $0<x<a$，且 $E<V_0$。在牆內，定態方程變為

$$ \frac{d^2\phi}{dx^2} = \kappa^2\phi, \qquad \kappa = \frac{\sqrt{2m(V_0-E)}}{\hbar}. $$

解為 $\phi(x)=Ce^{-\kappa x}+De^{+\kappa x}$。波函數不會在牆內變為零，而是以特徵長度 $1/\kappa$ 指數衰減；只要牆有限厚，另一側仍有非零振幅，粒子便有機率穿透。

在「厚而高」的牆（$\kappa a \gg 1$）近似下，穿透係數

$$ T \approx \frac{16E(V_0-E)}{V_0^2}\,e^{-2\kappa a} \;\sim\; e^{-2\kappa a}. $$

關鍵物理在指數項 $e^{-2\kappa a}$：穿隧機率對牆厚 $a$ 與「能量虧損」$\sqrt{V_0-E}$ 都呈指數敏感。這解釋了 α 衰變的壽命跨越數十個數量級、掃描穿隧顯微鏡（STM）對針尖—樣品間距達原子級的靈敏度，以及半導體穿隧二極體的運作。

定量小範例：電子穿隧的數量估計

設電子（$m=9.11\times10^{-31}\ \text{kg}$）能量 $E=1\ \text{eV}$，入射高 $V_0=2\ \text{eV}$、厚 $a=0.5\ \text{nm}$ 的位能牆，估算穿透係數。

步驟一：計算 $V_0-E = 1\ \text{eV} = 1.602\times10^{-19}\ \text{J}$。

步驟二：求衰減常數 $\kappa$。

$$ \kappa = \frac{\sqrt{2m(V_0-E)}}{\hbar} = \frac{\sqrt{2(9.11\times10^{-31})(1.602\times10^{-19})}}{1.055\times10^{-34}}. $$

分子 $=\sqrt{2.919\times10^{-49}}=5.40\times10^{-25}$，故 $\kappa = 5.40\times10^{-25}/1.055\times10^{-34}\approx 5.12\times10^{9}\ \text{m}^{-1}$。

步驟三：計算指數。$2\kappa a = 2\times(5.12\times10^{9})\times(0.5\times10^{-9}) \approx 5.12$，故 $e^{-2\kappa a}=e^{-5.12}\approx 5.97\times10^{-3}$。

步驟四：前因子 $16E(V_0-E)/V_0^2 = 16\times1\times1/4 = 4$。

合併得 $T \approx 4\times 5.97\times10^{-3} \approx 2.4\times10^{-2}$，約 2.4%。把牆加厚到 $a=1\ \text{nm}$，$2\kappa a\approx10.2$，$T$ 驟降到約 $4\times e^{-10.2}\approx1.5\times10^{-4}$——厚度加倍，穿透機率掉了兩個數量級，這正是指數律的威力。

深入探討（研究所視角）

把薛丁格方程放回更大的理論架構，能看出它在量子力學中的位置與限制。從哈密頓表述出發：古典哈密頓量 $H(q,p)=p^2/2m+V(q)$ 透過正則量子化 $\{q,p\}\to\frac{1}{i\hbar}[\hat q,\hat p]$、$[\hat x,\hat p]=i\hbar$，自然導出 $\hat H$；而薛丁格方程不過是態向量在希爾伯特空間中由么正演化算符 $U(t)=e^{-i\hat H t/\hbar}$ 驅動的微分形式。等價地，拉格朗日表述經費曼路徑積分 $\langle x_f|U|x_i\rangle = \int \mathcal{D}[x]\,e^{iS[x]/\hbar}$ 把振幅寫成對所有路徑的相位求和，作用量 $S=\int L\,dt$。在 $\hbar\to0$ 的古典極限下，穩相近似挑出使 $\delta S=0$ 的古典軌跡，呼應哈密頓—雅可比理論——這也給穿隧一個優雅的詮釋：穿隧對應「虛時間」（instanton）下的古典運動，衰減指數 $e^{-2\kappa a}$ 正是 WKB 近似中作用量積分 $\exp\!\big(-\frac{1}{\hbar}\int\sqrt{2m(V-E)}\,dx\big)$ 的展現。

WKB（Wentzel–Kramers–Brillouin）半古典方法是連結量子與古典的橋樑：對緩變位能，波函數寫成 $\phi\sim \exp(iS/\hbar)$ 並逐階展開 $\hbar$，在轉折點 $E=V(x)$ 處以 Airy 函數做連接公式，便能不解全方程而得到束縛態的 Bohr–Sommerfeld 量子化條件 $\oint p\,dx = (n+\tfrac12)2\pi\hbar$，以及一般位能牆的穿隧率。

往場論延伸，單粒子薛丁格方程的局限在於粒子數固定、且非洛倫茲協變。要描述粒子產生與湮滅、或高速情形，需把波函數提升為場算符 $\hat\psi(\mathbf{r})$ 並施加二次量子化的對易（玻色子）或反對易（費米子）關係——非相對論多體系統由此進入福克空間。相對論層面則由 Klein–Gordon 方程（自旋 0）與 Dirac 方程（自旋 1/2）接手；Dirac 方程的線性化要求自然導出自旋與反粒子，並在非相對論極限下還原薛丁格方程加上自旋—軌道與 Pauli 項。

這些表述在當代研究中交織成形。量子井模型是半導體異質結構、量子點與量子級聯雷射的理論基石；共振穿隧的相干輸運由 Landauer–Büttiker 形式描述，把穿透係數 $T$ 直接連到電導 $G=(2e^2/h)\sum_n T_n$。而路徑積分與 instanton 技術，則一路延伸到場論中的真空衰變、規範理論的非微擾效應與凝態系統的拓樸相變——一條看似樸素的二階微分方程，撐起了從原子到場的整片理論版圖。