粒子在 n 維空間中「正定勢場」的定態 Schrödinger 方程

1 引言

在《量子力學》¹這本書中，其中一節的內容就是關於「一維諧振子」的定態 Schrödinger 方程²。內容如下﹕

一維諧振子的勢場 $V(x)=\frac{1}{2}m{\omega }^2{x}^2$
其定態波函數 $\psi (x)$ 滿足 $-\frac{{\hslash }^2}{2m}\frac{d^2\psi }{d{x}^2}+\frac{1}{2}m{\omega }^2{x}^2\psi =E\psi$
解方程後得 ${\psi }_n(x)=N_n\exp (-\frac{{\alpha }^2{x}^2}{2})H_n(\alpha x)$

$N_n$ 是歸一化因子 $\sqrt{\frac{\alpha }{{\pi }^{\frac{1}{2}}\cdot 2^n\cdot n!}}$
$\alpha$ 是位移無量綱化因子 $\sqrt{\frac{m\omega }{\hslash }}$
$H_n(x)$ 是 $n$ 階 Hermite 多項式 $(-1{)}^n\exp (x^2)\frac{d^n}{d{x}^n}\exp (-x^2)$
$E_n$ 是粒子能級分佈 $(n+\frac{1}{2})\hslash \omega$。

我基於一維諧振子的理論，延伸到二維諧振子，逐步延伸 n 維諧振子，然後由二維諧振子延伸到二維「正定勢場」下的波函數，最後擴展到 n 維空間「正定勢場」下的波函數。

2 二維諧振子

二維諧振子下的勢場
$$V(x,y)=\frac{1}{2}m{\omega }^2(x^2+y^2)$$
其定態波函數 $\psi (x,y)$ 滿足
$$-\frac{{\hslash }^2}{2m}{\nabla }^2\psi (x,y)+\frac{1}{2}m{\omega }^2(x^2+y^2)\psi (x,y)=E\psi (x,y)$$
定義算子
$$\widehat{H_x}=-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial x^2}+\frac{1}{2}m{\omega }^2{x}^2$$
$$\widehat{H_y}=-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial y^2}+\frac{1}{2}m{\omega }^2{y}^2$$
利用分離變量法，設 $\psi (x,y)=f(x)g(y)$，所以有
$$\begin{gathered} \hat{H}\psi =g(y)\widehat{H_x}f(x)+f(x)\widehat{H_y}g(y)=Ef(x)g(y) \\ f(x)\widehat{H_y}g(y)=(Ef(x)-\widehat{H_x}f(x))g(y) \\ \widehat{H_y}g(y)=E_{y}g(y)\text{且}\widehat{H_x}f(x))=(E-E_y)f(x)=E_{x}f(x) \end{gathered}$$
由一維諧振子的結論可得，$\widehat{H_x}$ 和 $\widehat{H_y}$ 的本征值都是 { E 1 , n } n ≥ 0 \{E_{1,n}\}_{n\geq0} {E1,n​}n≥0​，對應的本征函數是 { ψ n ( x ) } n ≥ 0 \{\psi_n(x)\}_{n\geq0} {ψn​(x)}n≥0​ 和 { ψ n ( y ) } n ≥ 0 \{\psi_n(y)\}_{n\geq0} {ψn​(y)}n≥0​。
現在 $f$ 和 $g$ 分別是$\widehat{H_x}$ 和 $\widehat{H_y}$ 的本征函數，所以 f ∈ { ψ n ( x ) } n ≥ 0 f \in \{\psi_n(x)\}_{n\geq0} f∈{ψn​(x)}n≥0​ 和 g ∈ { ψ n ( y ) } n ≥ 0 g \in \{\psi_n(y)\}_{n\geq0} g∈{ψn​(y)}n≥0​，即${\psi }_{m+n}(x,y)={\psi }_m(x){\psi }_n(y)$。與此同時，對應的能級是 $E_{2,N}=E_{mn}=E_{1,m}+E_{1,n}=(m+n+1)\hslash \omega$，簡並度是 $\left(\begin{array}{c}N+1\\ 1\end{array}\right)=N+1$。

3 $n$ 維諧振子

以下用數學歸納法來證明如下命題﹕

當 $n$ 維諧振子的勢場為 $V(\overrightarrow{r_n})=V(x_1,x_2,\cdots ,x_n)=\frac{1}{2}m{\omega }^2{\overrightarrow{r_n}}^2$，
本征波函數族為 { ∏ k = 1 n ψ N k ( x k ) ∣ ∑ k = 1 n N k = N } \{\prod\limits_{k=1}^{n}\psi_{N_k}(x_k)|\sum\limits_{k=1}^{n}N_k=N\} {k=1∏n​ψNk​​(xk​)∣k=1∑n​Nk​=N}，其中 ${\psi }_k$ 是一維諧振子各能級 $E_k$ 對應的定態波函數。與此同時，對應的本征能級是 $E_{n,N}=(\sum _{k=1}^n{N}_k+\frac{n}{2})\hslash \omega =(N+\frac{n}{2})\hslash \omega$，簡並度是$\left(\begin{array}{c}N+n-1\\ n-1\end{array}\right)$

(1.) 當 $n=1$ 時，命題等價於一維諧振子結論
(2.) 假設 $n=t-1$ 時命題成立，那麼當 $n=t$ 時，設 $\psi (\overrightarrow{r_t})=\psi (\overrightarrow{r_{t-1}},x_t)=f(\overrightarrow{r_{t-1}})g(x_t)$。
記
$$\begin{gathered} {\nabla }_{\overrightarrow{r_{t-1}}}^2=\sum _{i=1}^{t-1}\frac{{\partial }^2}{\partial x_i^2} \\ \hat{H_{\overrightarrow{r_{t-1}}}}=-\frac{{\hslash }^2}{2m}{\nabla }_{\overrightarrow{r_{t-1}}}^2+\frac{1}{2}m{\omega }^2{\overrightarrow{r_{t_1}}}^2 \\ \widehat{H_{x_t}}=-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial x_t^2}+\frac{1}{2}m{\omega }^2{x}_t^2 \end{gathered}$$
所以
$$\hat{H}=-\frac{{\hslash }^2}{2m}{\nabla }_{\overrightarrow{r_t}}^2+\frac{1}{2}m{\omega }^2{\overrightarrow{r_t}}^2=\left(-\frac{{\hslash }^2}{2m}{\nabla }_{\overrightarrow{r_{t-1}}}^2+\frac{1}{2}m{\omega }^2{\overrightarrow{r_{t_1}}}^2\right)+\left(-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial x_t^2}+\frac{1}{2}m{\omega }^2{x}_t^2\right)=\hat{H_{\overrightarrow{r_{t-1}}}}+\widehat{H_{x_t}}$$
利用分離變量法，設 $\psi (\overrightarrow{r_t})=f(\overrightarrow{r_{t-1}})g(r_t)$，類似二維諧振子般，可以得到如下關係式
$$\widehat{H_{x_t}}g(y)=E_{x_t}g(x_t)\text{且}\hat{H_{\overrightarrow{r_{t-1}}}}f(\overrightarrow{r_{t-1}}))=(E-E_{x_t})f(\overrightarrow{r_{t-1}})=E_{\overrightarrow{r_{t-1}}}f(\overrightarrow{r_{t-1}})$$
由 (1.) 知，$\widehat{H_{x_t}}$ 的本征值是 { E 1 , N } N ≥ 0 = { ( N + 1 2 ) ℏ ω } N ≥ 0 \{E_{1,N}\}_{N\geq0}=\{(N+\frac{1}{2})\hbar\omega\}_{N\geq0} {E1,N​}N≥0​={(N+21​)ℏω}N≥0​，對應的本征函數是 { ψ N ( x t ) } N ≥ 0 \{\psi_N(x_t)\}_{N\geq0} {ψN​(xt​)}N≥0​。
由歸納假設知，$\hat{H_{\overrightarrow{r_{t-1}}}}$ 的本征值是 { E t − 1 , N } N ≥ 0 = { ( N + t − 1 2 ) ℏ ω } n ≥ 0 \{E_{t-1,N}\}_{N\geq0}=\{(N+\frac{t-1}{2})\hbar\omega\}_{n\geq0} {Et−1,N​}N≥0​={(N+2t−1​)ℏω}n≥0​，每一個本征值 $E_{t-1,N}$ 對應的本征函數族是 { ∏ k = 1 t − 1 ψ N k ( x k ) ∣ ∑ k = 1 t − 1 N k = N } \{\prod\limits_{k=1}^{t-1}\psi_{N_k}(x_k)|\sum\limits_{k=1}^{t-1}N_k=N\} {k=1∏t−1​ψNk​​(xk​)∣k=1∑t−1​Nk​=N}
因此，$\hat{H}$ 的本征值便是 { E t , N } N ≥ 0 = { E t − 1 , N + E 1 , 0 } N ≥ 0 = { ( N + t − 1 2 + 0 + 1 2 ) ℏ ω } N ≥ 0 = { ( N + t 2 ) } N ≥ 0 \{E_{t,N}\}_{N\geq0}=\{E_{t-1,N}+E_{1,0}\}_{N\geq0}=\{(N+\frac{t-1}{2}+0+\frac{1}{2})\hbar\omega\}_{N\geq0}=\{(N+\frac{t}{2})\}_{N\geq0} {Et,N​}N≥0​={Et−1,N​+E1,0​}N≥0​={(N+2t−1​+0+21​)ℏω}N≥0​={(N+2t​)}N≥0​，對應的本征值函數族是 { ∏ k = 1 t ψ N k ( x k ) ∣ ∑ k = 1 t N k = N } \{\prod\limits_{k=1}^{t}\psi_{N_k}(x_k)|\sum\limits_{k=1}^{t}N_k=N\} {k=1∏t​ψNk​​(xk​)∣k=1∑t​Nk​=N}

結合(1.) 和 (2.)，根據數學歸納原理可知，原命題對於任意正整數 $n$ 皆成立。

4 二維「正定勢場」下的定態 Schrödinger 方程

我先不引入「正定勢場」的概念，而是用另一個特殊勢場作為引入，求解其定態 Schrödinger 方程。
$$V(x,y)=\frac{1}{2}m{\omega }^2(x^2+xy+y^2)$$
考慮座標變換(讀者可以先自行思考)
$$\left(\begin{array}{c}u\\ v\end{array}\right)=\left(\begin{array}{cc}\frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}\\ \frac{1}{\sqrt{2}}& -\frac{1}{\sqrt{2}}\end{array}\right)\left(\begin{array}{c}x\\ y\end{array}\right)$$
用 $u$、$v$ 表示 ${\nabla }^2$﹕
$$\begin{gathered} \frac{\partial }{\partial x}=\frac{\partial }{\partial u}\frac{\partial u}{\partial x}+\frac{\partial }{\partial v}\frac{\partial v}{\partial x}=\frac{1}{\sqrt{2}}(\frac{\partial }{\partial u}+\frac{\partial }{\partial v}) \\ \frac{{\partial }^2}{\partial x^2}=\frac{1}{\sqrt{2}}(\frac{\partial }{\partial u}+\frac{\partial }{\partial v})\cdot \frac{1}{\sqrt{2}}(\frac{\partial }{\partial u}+\frac{\partial }{\partial v})=\frac{1}{2}(\frac{{\partial }^2}{\partial u^2}+\frac{{\partial }^2}{\partial v^2}+\frac{{\partial }^2}{\partial u\partial v}+\frac{{\partial }^2}{\partial v\partial u}) \end{gathered}$$
同理可得
$$\begin{gathered} \frac{{\partial }^2}{\partial y^2}=\frac{1}{2}(\frac{{\partial }^2}{\partial u^2}+\frac{{\partial }^2}{\partial v^2}-\frac{{\partial }^2}{\partial u\partial v}-\frac{{\partial }^2}{\partial v\partial u}) \\ {\nabla }^2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}=\frac{{\partial }^2}{\partial u^2}+\frac{{\partial }^2}{\partial v^2} \end{gathered}$$
簡記
$${\nabla }_{xy}^2={\nabla }_{uv}^2$$
用 $u$、$v$ 改寫勢場﹕
$$V=\frac{1}{2}m{\omega }^2(\frac{3}{2}{u}^2+\frac{1}{2}{v}^2)=\frac{1}{2}m(\sqrt{\frac{3}{2}}\omega {)}^2{u}^2+\frac{1}{2}m(\sqrt{\frac{1}{2}}\omega {)}^2{v}^2$$
記
$$\begin{gathered} \widehat{H_{x,\omega }}=-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial x^2}+\frac{1}{2}m{\omega }^2{x}^2 \\ {\psi }_{\omega ,n}(x)=\sqrt{\frac{\sqrt{\frac{m\omega }{\hslash }}}{{\pi }^{\frac{1}{2}}\cdot 2^n\cdot n!}}\exp (-\frac{\frac{m\omega }{\hslash }{x}^2}{2})H_n(\sqrt{\frac{m\omega }{\hslash }}x) \end{gathered}$$
所以在這個勢場下的哈密頓算符是
$$\hat{H}=\hat{H_{u,\sqrt{\frac{3}{2}}\omega }}+\hat{H_{v,\sqrt{\frac{1}{2}}\omega }}$$
如二維諧振子般利用分離變量法，可得本征波函數
$$\psi (x,y)={\psi }_{\sqrt{\frac{3}{2}}\omega ,m}(u){\psi }_{\sqrt{\frac{1}{2}}\omega ,n}(v)={\psi }_{\sqrt{\frac{3}{2}}\omega ,m}(\frac{x+y}{\sqrt{2}}){\psi }_{\sqrt{\frac{1}{2}}\omega ,n}(\frac{x-y}{\sqrt{2}})$$
對應的本征能級是
$$E_{m,n}=\left[(m+\frac{1}{2})\sqrt{\frac{3}{2}}+(n+\frac{1}{2})\sqrt{\frac{1}{2}}\right]\hslash \omega$$
顯然地，這個能級是非簡並的。

忽略能級簡並度的因素，我們去擴展這個問題﹕

如果勢函數 $V(x,y)$ 可以寫成
$$V(x,y)=\frac{1}{2}m{\omega }^2\left(\begin{array}{cc}x& y\end{array}\right)\left(\begin{array}{cc}a& b\\ b& c\end{array}\right)\left(\begin{array}{c}x\\ y\end{array}\right)$$
且當 $(x,y)\ne (0,0)$ 時 $V(x,y)>0$
那麼定態波函數 $\psi (x,y)$ 是否可以寫成 ${\psi }_{{\lambda }_1\omega ,m}(k_{1}x+b_{1}y)\cdot {\psi }_{{\lambda }_2\omega ,n}(k_{2}x+b_{2}y)$的形式呢？(其中，$k_1,b_1,k_2,b_2,{\lambda }_1,{\lambda }_2$ 都是常數)

注意到矩陣 $\left(\begin{array}{cc}a& b\\ b& c\end{array}\right)$ 是實對稱正定矩陣 (若讀者對正定矩陣不瞭解，可以參考 MIT 的 Introduction to Linear Algebra³ 書籍)，所以可以進行正交對角化。即存在一個正交矩陣 $Q=\left(\begin{array}{cc}{q}_{11}& q_{12}\\ q_{21}& g_{22}\end{array}\right)$ 和一個對角矩陣 $\Lambda =diag({\lambda }_1^2,{\lambda }_2^2)$，滿足
$$\begin{gathered} A=Q^T\Lambda Q \\ V(x,y)=\frac{1}{2}m{\omega }^2\left(\begin{array}{cc}x& y\end{array}\right)Q^T\Lambda Q\left(\begin{array}{c}x\\ y\end{array}\right)=\frac{1}{2}m{\omega }^2{\vec{q}}^T\Lambda \vec{q}=\frac{1}{2}m({\lambda }_1\omega {)}^2{u}^2+\frac{1}{2}m({\lambda }_2\omega {)}^2{v}^2 \end{gathered}$$
此時，$u=q_{11}x+q_{12}y$、$v=q_{21}x+q_{22}y$
下面證明﹕

$$\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}=\frac{{\partial }^2}{\partial u^2}+\frac{{\partial }^2}{\partial v^2}$$

證﹕

$$\begin{gathered} \frac{\partial }{\partial x}=\frac{\partial }{\partial u}\frac{\partial u}{\partial x}+\frac{\partial }{\partial v}\frac{\partial v}{\partial x}=q_{11}\frac{\partial }{\partial u}+q_{21}\frac{\partial }{\partial v} \\ \frac{{\partial }^2}{\partial x^2}=\frac{\partial }{\partial x}\cdot \frac{\partial }{\partial x}=(q_{11}\frac{\partial }{\partial u}+q_{21}\frac{\partial }{\partial v})(q_{11}\frac{\partial }{\partial u}+q_{21}\frac{\partial }{\partial v})=q_{11}^2\frac{{\partial }^2}{\partial u^2}+q_{21}^2\frac{{\partial }^2}{\partial u^2}+q_{11}{q}_{21}(\frac{{\partial }^2}{\partial u\partial v}+\frac{{\partial }^2}{\partial v\partial u}) \end{gathered}$$
同理可得
$$\frac{{\partial }^2}{\partial y^2}=q_{12}^2\frac{{\partial }^2}{\partial u^2}+q_{22}^2\frac{{\partial }^2}{\partial u^2}+q_{12}{q}_{22}(\frac{{\partial }^2}{\partial u\partial v}+\frac{{\partial }^2}{\partial v\partial u})$$
利用正交矩陣的定義
$$\begin{gathered} q_{11}^2+q_{12}^2=q_{21}^2+q_{22}^2=1 \\ {q}_{11}{q}_{21}+q_{12}{q}_{22}=0 \end{gathered}$$
可證得
$$\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}=\frac{{\partial }^2}{\partial u^2}+\frac{{\partial }^2}{\partial v^2}$$

利用上方結論，可以得出算符量
$$\begin{gathered} \hat{T}=\widehat{T_u}+\widehat{T_v} \\ \hat{V}=\widehat{V_{{\lambda }_1\omega }}+\widehat{V_{{\lambda }_2\omega }} \\ \hat{H}=\widehat{H_{{\lambda }_1\omega ,u}}+\widehat{H_{{\lambda }_2\omega ,v}} \end{gathered}$$
因此，本征波函數為
$$\psi (x,y)={\psi }_{{\lambda }_1\omega ,m}(u){\psi }_{{\lambda }_2\omega ,n}(v)={\psi }_{{\lambda }_1\omega ,m}(q_{11}x+q_{12}y){\psi }_{{\lambda }_2\omega ,n}(q_{21}x+q_{22}y)$$
對應的本征能級為
$$E_{m,n}=\left[(m+\frac{1}{2}){\lambda }_1+(n+\frac{1}{2}){\lambda }_2\right]\hslash \omega =\left(\begin{array}{cc}m+\frac{1}{2}& n+\frac{1}{2}\end{array}\right)\sqrt{\Lambda }\left(\begin{array}{c}1\\ 1\end{array}\right)\hslash \omega$$

5 正定勢場的定義

定義﹕ 對於任意向量 $\vec{r}=(x_1,x_2,\cdots ,x_n{)}^T\ne {\vec{0}}^T$，勢場
$$V(\vec{r})=\frac{1}{2}m{\omega }^2{\vec{r}}^{T}A\vec{r}>0$$
其中 $A$ 是一個實數矩陣，那麼勢場 $V$ 就是正定勢場。

一般而言，$A$ 會被寫成實對稱矩陣 (若非，則 $A$ 可用 $\frac{A+A^T}{2}$ 代替)，下一節的矩陣 $A$ 都會默認是實對稱矩陣。根據線性代數知識可知，$A$ 就是正定矩陣。

6 $n$ 維「正定勢場」下的定態 Schrödinger 方程

試着把二維情況推廣到 $n$ 維空間中，結論又是否會成立呢？
考慮
$$\begin{gathered} \vec{r}=(x_1,x_2,\cdots ,x_n{)}^T\ne {\vec{0}}^T \\ V(\vec{r})=\frac{1}{2}m{\omega }^2{\vec{r}}^T{A}_n\vec{r} \end{gathered}$$
$V$ 是一個正定勢場，尋找其本征波函數 $\psi (\vec{r})$。
首先，對 $A_n$ 進行正交對角化，同時改寫勢能表達式﹕
$$\begin{gathered} A_n=Q^T\Lambda Q \\ \Lambda =diag({\lambda }_1^2,{\lambda }_2^2,\cdots ,{\lambda }_n^2) \\ Q=(q_{ij}{)}_{n\times n}=\left(\begin{array}{c}{\overrightarrow{q_1}}^T\\ {\overrightarrow{q_2}}^T\\ ⋮\\ {\overrightarrow{q_n}}^T\end{array}\right) \\ V=\frac{1}{2}m{\omega }^2{\vec{r}}^T{Q}^T\Lambda Q\vec{r}=\frac{1}{2}m{\omega }^2{\vec{u}}^T\Lambda \vec{u}=\sum _{i=1}^n\frac{1}{2}m({\lambda }_i\omega {)}^2{u}_i^2 \\ (u_1,u_2,\cdots ,u_n{)}^T=\vec{u}=Q\vec{r} \end{gathered}$$
然後，我們證明一個結論﹕

若存在一個正交矩陣 Q = { q i j } n × n Q=\{q_{ij}\}_{n\times n} Q={qij​}n×n​，滿足 $\vec{u}=Q\vec{r}=\left(\begin{array}{c}{\overrightarrow{q_1}}^T\\ {\overrightarrow{q_2}}^T\\ ⋮\\ {\overrightarrow{q_n}}^T\end{array}\right)\vec{r}$，那麼
$${\nabla }_{\vec{r}}^2={\nabla }_{\vec{u}}^2(\sum _{i=1}^n\frac{{\partial }^2}{\partial x_i^2}=\sum _{i=1}^n\frac{{\partial }^2}{\partial u_i^2})$$
證﹕
$$\begin{gathered} \frac{\partial }{\partial x_i}=\sum _{k=1}^n\frac{\partial }{\partial u_k}\frac{\partial u_k}{\partial x_i}=\sum _{k=1}^n{q}_{ki}\frac{\partial }{\partial u_k} \\ \frac{{\partial }^2}{\partial x_i^2}=(\sum _{k=1}^n{q}_{ki}\frac{\partial }{\partial u_k})(\sum _{m=1}^n{q}_{mi}\frac{\partial }{\partial u_m})=\sum _{k=1}^n\sum _{m=1}^n{q}_{ki}{q}_{mi}\frac{{\partial }^2}{\partial u_k\partial u_m} \\ {\nabla }_{\vec{r}}^2=\sum _{i=1}^n\frac{{\partial }^2}{\partial x_i^2}=\sum _{i=1}^n\sum _{k=1}^n\sum _{m=1}^n{q}_{ki}{q}_{mi}\frac{{\partial }^2}{\partial u_k\partial u_m}=\sum _{k=1}^n\sum _{m=1}^n\left[\frac{{\partial }^2}{\partial u_k\partial u_m}\left(\sum _{i=1}^n{q}_{ki}{q}_{mi}\right)\right]=\sum _{k=1}^n\sum _{m=1}^n\left[\frac{{\partial }^2}{\partial u_k\partial u_m}\left({\overrightarrow{q_k}}^T\overrightarrow{q_m}\right)\right]=\sum _{k=1}^n\sum _{m=1}^n\left[\frac{{\partial }^2}{\partial u_k\partial u_m}{\delta }_{km}\right]=\sum _{k=1}^n\frac{{\partial }^2}{\partial u_k^2}={\nabla }_{\vec{u}}^2 \end{gathered}$$

因此，對應的算符量是
$$\begin{gathered} \widehat{T_i}=-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial u_i^2} \\ \widehat{V_{{\lambda }_i\omega ,i}}=\frac{1}{2}m({\lambda }_i\omega {)}^2{u}_i^2 \\ \widehat{H_{{\lambda }_i\omega ,i}}=\widehat{T_i}+\widehat{V_{\lambda \omega ,i}} \\ \hat{H}=\sum _{i=1}^n\widehat{H_{{\lambda }_i\omega ,i}} \end{gathered}$$
如第 3 節的「$n$ 維諧振子」般利用數學歸納法可得，本征波函數
$$\psi (\vec{r})=\prod _{i=1}^n{\psi }_{{\lambda }_i\omega ,m_i}({\overrightarrow{q_i}}^T\vec{r})$$
對應的本征能級 $\vec{m}=(m_i{)}_{n\times 1}^T$ 是
$$E_{\vec{m}}=\hslash \omega \sum _{i=1}^n(m_i+\frac{1}{2}){\lambda }_i=({\vec{m}}^T+(\frac{1}{2}{)}_{1\times n})\sqrt{\Lambda }(1{)}_{n\times 1}\hslash \omega$$

7 資料來源

---

1. 《量子力學》，錢伯初，高等教育出版社，2017年。 ↩︎

2. 一維諧振子https://en.wikipedia.org/wiki/Schr%C3%B6dinger_equation#Harmonic_oscillator ↩︎

3. Gilbert Strang, Introduction to Linear Algebra, fifth edition ↩︎