水素原子におけるシュレーディンガー方程式の解

本項、**水素原子におけるシュレーディンガー方程式の解**（すいそげんしにおけるシュレーディンガーほうていしきのかい）では、ハミルトニアンが
と書ける二粒子系の時間非依存なシュレーディンガー方程式の厳密解を解く（式中の記号の意味は後述）。

物理学的にはこれは、
* 質量の正の電荷をもつ粒子と質量が負の電荷を持つ粒子がクーロン力により結合している状況において
* 外力は働いておらず、
* 相対論的効果を考えない量子力学の範囲内で、
* 時間に依存しない定常状態の
粒子の波動関数を決定する事を意味する。正の電荷をもつ粒子と負の電荷がそれぞれ陽子と電子だとすればこの系は**[[水素原子]]**に相当するが、一般の価数の原子核を持つ1電子系多価イオン（**水素様原子**）の系も同一の方程式から解を導ける。この方程式は様々な教科書で取り上げられている。

なお、微細構造、超微細構造、ラムシフトなどの効果は、いずれも相対論的な量子力学を必要とする為、本項の対象外である。

## シュレーディンガー方程式
本項の目的は、時間非依存なシュレディンガー方程式
でハミルトニアンが
と書ける場合の厳密解を求める事である。ここで

ここで
:\boldsymbol{x}_0=(x_0,y_0,z_0)
:\boldsymbol{x}_1=(x_1,y_1,z_1)

はの元であり、、、は正の定数であり、
:\Delta_j =
{\partial^2\over \partial x_j{}^2}
+{\partial^2\over \partial y_j{}^2}
+{\partial^2\over \partial z_j{}^2}、
であり、は換算プランク定数である。

前述した物理的状況においては、２つの粒子の電荷をそれぞれとし、真空の誘電率をとすれば、
:Q={e_1e_2\over 4\pi \varepsilon_0}

であるが、本項では一般の正の定数に対して解を導くので、必ずしもが上述の形である事を仮定しない。

### 重心系への還元
()、()により定義される方程式は、重心系に書き直す事により、より簡単な式に還元できる。２つの粒子の重心
:\boldsymbol{c}={m_0\boldsymbol{x}_0+m_1\boldsymbol{x}_1\over m_0+m_1}

と２つの粒子の位置の差
:\boldsymbol{x}=\boldsymbol{x}_1-\boldsymbol{x}_0

と換算質量

:\mu=\frac{m_0m_1}{m_0+m_1}

を使うと、ハミルトニアン()は
:\hat{H}=
-{\hbar^2 \over 2 (m_0+m_1)}\Delta_{\boldsymbol{c}}
 -{\hbar^2 \over 2 \mu} \Delta_{\boldsymbol{x}}
 - {Q \over |\boldsymbol{x}|} 

と書ける<sup>H13</sup>。

このハミルトニアンは
　　…()|border colour=#50C878|background colour=#ECFCF4}}

=
-
{\hbar^2 \over 2 \mu}\Delta_{\boldsymbol{x}}
- {Q \over |\boldsymbol{x}|}
　　…()|border colour=#50C878|background colour=#ECFCF4}}
の和である。
  
()のハミルトニアンはよく知られた自由粒子のハミルトニアンであり、その連続スペクトルは
:\sigma_c(\hat{H}_c)=[0,\infty)

であり、点スペクトルは
:\sigma_c(\hat{H}_c)=\emptyset

である<sup>H13</sup>。したがって後は非自明な部分である()のスペクトルを求めれば良いことになる<sup>新井</sup>。そこで以下()のみ焦点を当てる。

### 無次元化
適切な値と定数を選び、長さとエネルギーをそれぞれ、が1となるように座標変換
:(x',y',z')= (x/a_0, y/a_0, z/a_0)
:E'=E/E_a
してやると、()のハミルトニアンに関する時間非依存なシュレディンガー方程式は
と無次元化される<sup>SO96</sup>

簡単な計算により、の具体的な値は
:a_0={\hbar^2 \over \mu Q}、E_a={\mu Q^2\over \hbar^2}　　…()
である事が分かる。

#### ボーア半径・ハートリー
特に、陽子の質量が電子の質量より遥かに重いと仮定した場合の水素原子の系における、は
:Q={e^2\over 4\pi \varepsilon_0}
:\mu =m_1\left( 1+ {m_1\over m_0}\right)\approx m_1
より、
:a_0={4\pi \varepsilon_0\hbar^2 \over m_1e^2}、E_a={m_1e^4\over 16\pi^2 \varepsilon_0^2 \hbar^4}
である。ここでは電気素量である。この場合のを**[[ボーア半径]]**といい、を基準としたエネルギーの単位を**[[ハートリー]]**という<sup>SO96</sup>。

## 求解

本節では()のハミルトニアンを無次元した

のスペクトルを求める。なお、本節ではまず変数分離解を求めるが、後述するように実はこのハミルトニアンは変数分離解しか持たない。

### 求解の方針
()を解く基本的アイデアは、無次元化した座標系を球面座標に変換するというものだが、直接球面座標を用いると、計算が複雑になる。そこで計算を楽にするため、以下の事実に着目する。

()のハミルトニアンは球対称なポテンシャルを持っており、しかもラプラシアンは回転不変である事が知られているので、()のハミルトニアンは回転不変である。よって()のハミルトニアンは軌道角運動量演算子(\hat{L}_x,\hat{L}_y,\hat{L}_z)と可換である：
:[\hat{H}_{\boldsymbol{x}'},\hat{L}_x]
=[\hat{H}_{\boldsymbol{x}'},\hat{L}_y]
=[\hat{H}_{\boldsymbol{x}'},\hat{L}_z]
=0  
よって特に、軌道角運動量演算子の自乗とも可換である：
:[\hat{H}_{\boldsymbol{x}'},\hat{\boldsymbol{L}^2}]=0  
よってはと同時対角化できるはずである、さらに
:[\hat{\boldsymbol{L}^2},\hat{L}_z]=0  
である事から、の３つを同時対角化できるはずである。

そこでまず、の同時固有関数を求め、これを利用しての固有関数を求める。

### との同時固有関数
との同時固有関数の求め方は「軌道角運動量」の項目に書いてあるので、結論だけを言えば、に対し、
:\hat{\boldsymbol{L}^2}\psi =\hbar^2 \ell(\ell+1)\psi  
:\hat{L}_z\psi =m\hbar \psi  
を満たす固有関数が存在し、は極座標で
:
\psi(r',\theta,\varphi)= 
 R(r')
P_\ell^(\cos \theta) \, e^{im\phi}
 　　×(規格化定数)　…()
という形で書ける。ここでは任意の自乗可積分関数であり、はルジャンドルの陪多項式
: P_\ell{}^m (z) = {1\over 2^\ell \ell!}(1-z^2)^{m\over 2}{\operatorname{d}^{m+\ell}\over \operatorname{d}z^{m+\ell}}(z^2-1)^{\ell}
 
である<sup>新井</sup>。

### の決定
後はを決定するだけである。を決定するには()を()のハミルトニアンに入れてシュレディンガー方程式()を解けば良い。（）を式変形すると、
:
\left({1\over 2}\Delta_{\boldsymbol{x}'}+{1 \over |\boldsymbol{x}'|}+E'\right)\psi=0

　　 …()

である。ラプラシアンを球面座標で書き表し、動径方向と球面方向にわけると、
: \Delta_{\boldsymbol{x}'}={1\over r'^2}(\Delta_{r'}+\Delta_S)　   …()
と書ける<sup>武藤11-15</sup>。ここで
: \Delta_{r'}=\frac{\partial}{\partial r'}\left(r'^2\frac{\partial}{\partial r'}\right)  、\Delta_S=-{1\over \hbar^2}\hat{\boldsymbol{L}^2}　　　…()
であり<sup>武藤11-15</sup>、は軌道角運動量演算子の自乗である。()のラプラシアンを極座標表示した上で()に()の波動関数を代入すると、()がの固有値に対応する固有関数であった事から、
: {1\over 2r'^2}\left(\Delta_{r'}-\ell(\ell+1)
+2r+2r^2E'\right)
 R(r')
P_\ell^(\cos \theta) \, e^{im\phi}=0
　

すなわち
: \left(\Delta_{r'}-\ell(\ell+1)
+2r'+2r'^2E'\right)
 R(r)
=0
　

束縛状態ではは負の値しか取らないので、記号を簡単にするため
:n:={1\over \sqrt{-2E'}}, \rho:={2r' \over n}　　　…()
と定義し<sup>原94</sup>、をの関数とみなすと、
:\rho^2 \frac{\operatorname{d}^2R}{\operatorname{d}\rho^2} + 2\rho \frac{\operatorname{d}R}{\operatorname{d}\rho} 
- \left\{ \frac{\rho^2}{4} - n \rho + \ell(\ell+1) \right\} R = 0

　…()
が成立する<sup>石川15</sup>。

この方程式を解くのは複雑な計算を必要とするので後の章にまわし、ここでは結論のみを述べる。

()の方程式を解くことで各に対し、（()のを単位とした）エネルギー
:E'_n= -\frac{1}{2n^2}                   ...()
に対する解が見つかる<sup>新井</sup>。に対応する固有関数は、


\begin{cases}
0\le \ell \le n-1\\
|m|\le\ell
\end{cases}
　　　　　　　　　　　　　　　　　　　　　…()

に対してのみ存在し、そのときのはラゲールの陪関数
:
R_{n,\ell}(\rho)=
 \exp \left( - \rho\right)
\rho^\ell L^{2\ell+1}_{n+\ell}\left(2\rho\right)
　　　　　×規格化定数　　　　　…()
に一致する。ここで
:
\rho={r'\over n}
、
:
L^m_k(\rho) = \frac{\operatorname{d}^m}{\operatorname{d}\rho^m} e^\rho \frac{\operatorname{d}^k}{\operatorname{d}\rho^k} (e^{-\rho} \rho^k)


である。

### 規格化定数
3次元空間における体積要素は動径方向の線素と球面方向の面素を用いて
: \operatorname{d}V
=
r'^2
\operatorname{d}r'
\operatorname{d}S

と書けるので、()におけるのノルム
: \|\psi\|{}:=
\sqrt{\int_0^{\pi} |\psi(x',y',z')|^2\operatorname{d}V}
　　　
も
: \|\psi\|{}:=\|R\|_r\|Y\|_S

　　                         　…()
と「変数分離」する。ここで

: Y(\theta,\varphi)=P_\ell^(\cos \theta) \, e^{im\phi}

　　
であり、
: \|Y\|_S{}:=
\sqrt{\int_0^{\pi} |Y(\theta,\phi)|^2\operatorname{d}S}
　　　　…()
: \|R\|_r{}:=
\sqrt{\int_0^\infty |R(r)|^2r^2\operatorname{d}r}

　　　　…()
()のノルムを1にする規格化定数の値は「軌道角運動量」の項目に書いてあり、
: (-1)^{(m+|m|)/2}
         \sqrt{ \frac{2\ell+1}{4\pi} \frac{(\ell-|m|)!}{(\ell+|m|)!} \,}
 

である<sup>原94</sup>。

()のノルムを1にする規格化定数の値の計算は後述するが、結論から言えば規格化定数は
:\left({2\over n}\right)^{3/2}\sqrt{ \frac{(n-\ell-1)!}{2n(n+\ell)!} }　　　　　...()
である。

## 結論
無次元化した()をベースにしたこれまでの議論を通常の単位系に戻すことで以下の結論が得られる。
:a_0={\hbar^2 \over \mu Q}
とし、を自然数、を以下を満たす整数とする：


\begin{cases}
0\le \ell \le n-1\\
|m|\le\ell
\end{cases}
　　　　　　　　　　…(****、再掲)　　　　　

このとき()のハミルトニアンはエネルギー
:E_n= -\frac{\mu Q^2}{2\hbar^2n^2}
に対し、
:\hat{\boldsymbol{L}^2}\psi =\hbar^2 \ell(\ell+1)\psi  
:\hat{L}_z\psi =m\hbar \psi  
を満たす固有関数
:
\psi_{n,\ell,m}(r,\theta, \phi)      
=  R_{n,\ell}(r) Y_{lm}(\theta, \phi)=
 \exp ( - \rho_n )\rho_n{}^{\ell}
 L^{2\ell+1}_{n+\ell}(2\rho_n)
P_\ell^(\cos \theta) \, e^{im\phi}
 　×(規格化定数)　　…()
を持つ。ここで
: P_\ell{}^m (z) = {1\over 2^\ell \ell!}(1-z^2)^{m\over 2}{\operatorname{d}^{m+\ell}\over \operatorname{d}z^{m+\ell}}(z^2-1)^{\ell}
、 
: 
L^m_k(\rho) = \frac{\operatorname{d}^m}{\operatorname{d}\rho^m} e^\rho \frac{\operatorname{d}^k}{\operatorname{d}\rho^k} (e^{-\rho} \rho^k)

 
: 
\rho_n={r \over na_0}

であり、規格化定数は
: (-1)^{(m+|m|)/2}
 \sqrt{
 \left({2\over na_0}\right)^{3} 
 \frac{2\ell+1}{4\pi} \frac{(\ell-|m|)!}{(\ell+|m|)!} 
 \frac{(n-\ell-1)!}{2n(n+\ell)!} }

である。

以上では変数分離により発見的に解を求めたため、()、()に書いたものが解である事は間違いないものの、それ以外に解があるかどうかは不明である。しかし実は**これ以外に解がない事が知られている**<sup>H13</sup>。)=(0,\infty)
と点スペクトル
: \sigma_p(\hat{H}_{\boldsymbol{x}})=\left\{E_n \mid n=1,2,\ldots \right\}
を持ち、
E_nに対する固有関数は()、()で書かれた関数で貼られる次元空間である。|name=定理}}連続スペクトルに相当する部分は、物理的にいえば水素原子がイオン化している状態であり、したがって電子が陽子から逃れていってしまっている<sup>H13</sup>。なお、固有関数の和
:
\sum_{n,\ell,m} a_{n,\ell,m}\psi_{n,\ell,m}(r,\theta, \phi)

  s.t. 
\sum_{n,\ell,m} a_{n,\ell,m}{}^2  
の形に書けるのは、の負のスペクトルに対応するベクトルだけで、正のスペクトルに対応するベクトルはこの方法では表記できない<sup>H13</sup>。

### 量子数
ハミルトニアン()の固有関数()に登場する２つの変数は以下のように呼ばれる：
* は**[[主量子数]]**と呼ばれ、のエネルギー固有値の大きさを司っている。
* は'''軌道角運動量量子数'''（方位量子数）と呼ばれ、の固有値の大きさを司っている。
* は**[[磁気量子数]]**（軌道磁気量子数）と呼ばれ、の固有値の大きさを司っている。

なお、は、動径方向の波動関数の節の数を表している。

#### 化学的意味

3つの量子数のうち、には以下のような化学的意味がある：
* 主量子数  は電子殻の **K殻、L殻、M殻**、…に対応している。
* 方位量子数  は**[[s軌道]]、[[p軌道]]、[[d軌道]]、f軌道、g軌道**…に対応している。
水素原子において、s軌道, p軌道, d軌道, f軌道…のエネルギー準位は縮退している。これはエネルギー固有値が、となり、やに依存しないためである。なお、水素原子に磁場をかけると、これらのエネルギー準位は、スピン部分を無視して考えた場合、磁気量子数の違いにより分裂する(→ゼーマン効果)。電場をかけた場合も、シュタルク効果によって分裂する。このとき、異なるの軌道同士の線形結合をとった混成軌道がハミルトニアンの固有状態となる。

### リュードベリ定数

エネルギー準位がにある電子がエネルギー準位がに落ちると、
:E_{n}-E_{n'}= \frac{\mu Q^2}{2\hbar^2}\left({1\over n'^2}-{1\over n^2}\right)
のエネルギーが
:E_{n}-E_{n'}= {\hbar c \over \lambda} 
を満たす波長の光となって放出される。したがって
:{1\over \lambda}= \frac{\mu Q^2}{2\hbar^3c}\left({1\over n'^2}-{1\over n^2}\right)

水素原子の場合、すなわち
:Q={e^2\over 4\pi \varepsilon_0}
の場合の上式右辺の定数、もしくはその定数に対して近似
:\mu =m_1\left( 1+ {m_1\over m_0}\right)\approx m_1
を行ったときの値を**[[リュードベリ定数]]**という。

## ()の解
本節の目的は、微分方程式()を解き、()、()を導出することである。

### ラゲールの陪方程式にあてはめる

本節では式()をさらに式変形することで、()をラゲールの陪方程式（詳細後述）で書き表せる事を示す。ラゲールの陪方程式の解は特殊関数で書けることが知られているので、これにより式()が解けることになる。この目標に達するため、以下の3ステップを踏む。
* が十分小さいという条件下()の近似解を求める。
* **'が十分大きいという条件下()の近似解を求める。
* 上記2ステップの結論を参考にして、()の厳密解を変数変換し、()をラゲールの陪方程式に（近似なしで）変形する。

==== が十分小さい場合の()の近似解 ====

()におけるの係数はが十分小さいところではと近似できるので、()は
:<math>\rho^2 \frac{\operatorname{d}^2R}{\operatorname{d}\rho^2} + 2\rho \frac{\operatorname{d}R}{\operatorname{d}\rho} 
- \ell(\ell+1) R = 0

</math>
と近似できる[[#石川15|<sup>石川15</sup>]]。

この形の方程式は[[オイラーの微分方程式]]の解法に準ずる方法で解ける。その解は



の形で書ける。

==== が十分大きい場合の()の近似解 ====

()式をで割った上での極限をとることで、が十分大きいところでは()は



となる事がわかる。簡単な計算から上記の方程式の一般解は
:<math>R(\rho)=\mathrm{e}^{\rho\over 2}</math>、<math>\mathrm{e}^{-\rho\over 2}</math>
もしくはこれらの線形和である。は発散する不適切な解となるので、



である。

==== ()からのラゲールの陪多項式の導出 ====
()、()を参考に、()の厳密解R(&rho;)を



の形に変数変換する。一般に3つの関数の積の微分は公式
:<math>\begin{align}
(fgh)' &= f'gh + fg^'h + fgh' \\
(fgh)'' &= (f''gh + fg''h + fgh'')
+ 2 (f'g'h + fg'h' + f'gh')
\end{align}
</math>

を満たすので、()の第一項、および第二項は、

:<math>\begin{array}{lccl}
\displaystyle 
\rho^2 \frac{\operatorname{d}^2R}{\operatorname{d}\rho^2} &=& \frac{\operatorname{d}}{\operatorname{d}\rho}\left\{\rho^\ell u(\rho) \mathrm{e}^{-\frac{\rho}{2}}\right\}
& = & \ell\rho^{\ell-1}u(\rho)\mathrm{e}^{-\frac{\rho}{2}}
+ \rho^{\ell}\rho' \mathrm{e}^{-\frac{\rho}{2}}
- \frac{1}{2}\rho^\ell u(\rho) \mathrm{e}^{-\frac{\rho}{2}} \\

\displaystyle \rho \frac{\operatorname{d}R}{\operatorname{d}\rho} &=& \frac{\operatorname{d}^2}{\operatorname{d}\rho^2}\left\{\rho^\ell u(\rho) \mathrm{e}^{-\frac{\rho}{2}}\right\}
& = & \displaystyle \ell(\ell-1)\rho^{\ell-2} u(\rho) \mathrm{e}^{-\frac{\rho}{2}}
+ \rho^\ell \rho''\mathrm{e}^{-\frac{\rho}{2}}
+\frac{1}{4}\rho^\ell u(\rho)\mathrm{e}^{-\frac{\rho}{2}}\\

&&         & + & \displaystyle 2 \left\{ \ell\rho^{\ell-1}\rho'\mathrm{e}^{-\frac{\rho}{2}}
         - \frac{1}{2}\rho^{\ell}\rho'\mathrm{e}^{-\frac{\rho}{2}}
         - \frac{1}{2}\ell\rho^{\ell-1}u(\rho)\mathrm{e}^{-\frac{\rho}{2}}\right\}
\end{array}
</math>　　

である。上式を()に代入すると、すべての項にが掛かっていることがわかる。よって各項をで割った上で式を整理して、



を得る。この式の両辺をで割ると、



となる[[#石川15|<sup>石川15</sup>]]。こうして得た式(6.12)は下記の式(6.13)に示したラゲールの陪方程式（ラゲール陪関数）の形になっている。
:<math>
\rho\frac{\operatorname{d}^2u(\rho)}{\operatorname{d}\rho^2}+(m+1-\rho)\frac{\operatorname{d}u(\rho)}{\operatorname{d}\rho}+(k-m)u(\rho) = 0
</math>     …()

ラゲールの陪方程式の解は**ラゲールの陪多項式'''（ラゲール陪関数）と呼ばれる形の定数倍になることが知られている。ラゲールの陪多項式（ラゲール陪関数）は下記のように定義される。
:
L^m_k(\rho) = \frac{\operatorname{d}^m}{\operatorname{d}\rho^m} e^\rho \frac{\operatorname{d}^k}{\operatorname{d}\rho^k} (e^{-\rho} \rho^k) \quad \rho={2r\over n}


ここで、は
:
0\le k \le m
   …()

を満たす整数である。

よって()の解は
:
u(\rho) = L^{2\ell+1}_{n+\ell}(\rho)
   ×(規格化定数)
となる。これを変数変換の式()に代入して
: R(\rho) = \rho^\ell L^{2\ell+1}_{n+\ell}(\rho) \exp\left(-\frac{\rho}{2}\right)        ×(規格化定数)   …()
を得る。

ラゲール陪多項式の係数の条件式()から、hoge\ellは
:
0\le \ell \le n-1
   …()

を満たす整数でなければならない。

## 規格化定数()の導出
規格化定数をとすると、規格化条件
:
\int_{0}^{\infty} |R(r')|^2 r'^2 \operatorname{d}r' = 1
 

は、()、()より、

:
\begin{array}{lcl}
\displaystyle 1=\int_{0}^{\infty} |R(r')|^2 r'^2 \operatorname{d}r
 & = &\displaystyle  {n^3\over 8}\int_{0}^{\infty} R(\rho)^2 \rho^2 \operatorname{d}\rho\\ 
& = & \displaystyle  {n^3C'\over 8} \int_{0}^{\infty} \rho^{2\ell+2} \{L^{2\ell+1}_{n+\ell}(\rho)\}^2 \exp\left(-\rho\right) \operatorname{d}\rho\\
\\
\end{array}
 …()

ラゲールの陪多項式（ラゲール陪関数）は下記の直交性を満たすことが知られている

:\begin{align}
\int_{0}^{\infty} z^m \exp(-z) L_{k}^{m}(z) L_{\ell}^{m}(z) \operatorname{d}z &= \frac{(k!)^3}{(k-m)!} \delta_{k\ell} \\
\int_{0}^{\infty} z^{m+1} \exp(-z) \{ L_{k}^{m}(z) \}^2 \operatorname{d}z &= (2k+1-m) \frac{(k!)^3}{(k-m)!}
\end{align}

ので、後者の式を()に対して用いる事で、
: \displaystyle \int_{0}^{\infty} \rho^{2\ell+2} \{L^{2\ell+1}_{n+\ell}(\rho)\}^2 \exp\left(-\rho\right) \operatorname{d}\rho 
=\displaystyle (2n) \frac{[(n+\ell)!]^3}{(n-\ell-1)!}



これが()の左辺である1と等しいことから、規格化定数について解く事で
:C' = \left({2\over n}\right)^{3/2}\sqrt{ \frac{(n-\ell-1)!}{2n[(n+\ell)!]^3} }　　…()

が得られる。

なお、無次元化する前のハミルトニアン()に対する規格化定数は、変数変換
:
\int_{0}^{\infty} |R(r)|^2 r^2 \operatorname{d}r={1\over a_0}^3\int_{0}^{\infty} |R(r')|^2 r'^2 \operatorname{d}r'
   

の分だけ()のものとはずれるので、()に対する規格化定数は、
:C = \left({2\over na_0}\right)^{3/2}\sqrt{ \frac{(n-\ell-1)!}{2n[(n+\ell)!]^3} }　　…()

となる<sup>原94</sup>。

## 具体的な値

水素原子の波動関数のにおける角因子は以下のようになる。ここで、はそれぞれ、動径方向の関数
: Y(\theta,\varphi)=P_\ell^(\cos \theta) \, e^{im\varphi}

　
の右辺の積の第一成分と第二成分を規格化したものである。なお、の指数関数の虚数部分はオイラーの公式により一対の関数の一次結合で書き換えられる。

原子番号の水素様原子の動径関数は以下のようになる。
:\begin{align}
  R_{1 \mathrm s} &= 2 \left( \frac{Z}{a_0} \right)^{3/2} \exp \left( - \frac{Zr}{a_0} \right)
  \\
  R_{2 \mathrm s} &= \frac{1}{2 \sqrt{2}} \left( \frac{Z}{a_0} \right)^{3/2} \left( 2 - \frac{Zr}{a_0} \right)
  \exp \left( - \frac{Zr}{2a_0} \right) \\
  R_{2 \mathrm p} &= \frac{1}{2 \sqrt{6}} \left( \frac{Z}{a_0} \right)^{3/2} \frac{Zr}{a_0}
  \exp \left( - \frac{Zr}{2a_0} \right)
  \\
  R_{3 \mathrm s} &= \frac{2}{81 \sqrt{3}} \left( \frac{Z}{a_0} \right)^{3/2}
  \left( 27 - \frac{18 Zr}{a_0} + \frac{2 Z^2 r^2}{a_0^2} \right)
  \exp \left( - \frac{Zr}{3a_0} \right) \\
  R_{3 \mathrm p} &= \frac{4}{81 \sqrt{6}} \left( \frac{Z}{a_0} \right)^{3/2}
  \left( 6 - \frac{Zr}{a_0} \right) \frac{Zr}{a_0}
  \exp \left( - \frac{Zr}{3a_0} \right) \\
  R_{3 \mathrm d} &= \frac{4}{81 \sqrt{30}} \left( \frac{Z}{a_0} \right)^{3/2}
  \frac{Z^2 r^2}{a_0^2}
  \exp \left( - \frac{Zr}{3a_0} \right)
  \\
  R_{4 \mathrm s} &= \frac{1}{768} \left( \frac{Z}{a_0} \right)^{3/2}
  \left( 192 - \frac{144 Zr}{a_0} + \frac{24 Z^2 r^2}{a_0^2} - \frac{Z^3 r^3}{a_0^3} \right)
  \exp \left( - \frac{Zr}{4a_0} \right) \\
  R_{4 \mathrm p} &= \frac{1}{256 \sqrt{15}}
  \left(\frac{Z}{a_0} \right)^{3/2} \left(80 - \frac{20 Zr}{a_0} + \frac{Z^2 r^2}{a_0^2} \right)
  \frac{Zr}{a_0} \exp \left( - \frac{Zr}{4a_0} \right) \\
  R_{4 \mathrm d} &= \frac{1}{768 \sqrt{5}} \left(\frac{Z}{a_0} \right)^{3/2} \left(12 - \frac{Zr}{a_0} \right)
  \frac{Z^2 r^2}{a_0^2} \exp \left( - \frac{Zr}{4a_0} \right) \\
  R_{4 \mathrm f} &= \frac{1}{768 \sqrt{35}} \left(\frac{Z}{a_0} \right)^{3/2} \frac{Z^3 r^3}{a_0^3}
  \exp \left( - \frac{Zr}{4a_0} \right)
\end{align}

動径関数を2乗しを掛けた動径分布は、核の中心からのある距離における電子の存在確率に相当する。

詳しくは電子配置の項を参照のこと。

## 脚注
### 注釈

### 出典

## 参考文献
* 書籍
** [新井97] 
** [原94] 
** [H13] 
** [SO96] 
*** 邦訳：
* レクチャーノート
** [武藤11-15]
** [石川15] 

## 関連項目
* シュレーディンガー方程式
* 球面調和関数
* ラゲールの陪多項式
* 水素原子

## 外部リンク
* 水素原子の電子分布の計算