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

{{Pathnav|[[物理学]]|[[量子力学]]|frame=1}}

本項、'''水素原子におけるシュレーディンガー方程式の解'''（すいそげんしにおけるシュレーディンガーほうていしきのかい）では、[[ハミルトニアン]]が{{Equation box 1|cellpadding
|border
|indent=:|equation=<math>\hat{H}=-{\hbar^2 \over 2 m_0}\Delta_0 -{\hbar^2 \over 2 m_1} \Delta_1 - {Q \over |\boldsymbol{x}_0-\boldsymbol{x}_1|}</math>|border colour=#50C878|background colour=#ECFCF4}}
と書ける二粒子系の時間非依存な[[シュレーディンガー方程式]]の厳密解を解く（式中の記号の意味は後述）。

物理学的にはこれは、
* [[質量]]{{Math|''m''{{sub|0}}}}の正の[[電荷]]をもつ粒子と質量が{{Math|''m''{{sub|1}}}}負の電荷を持つ粒子が[[クーロン力]]により結合している状況において
* 外力は働いておらず、
* 相対論的効果を考えない量子力学の範囲内で、
* 時間に依存しない定常状態の
粒子の[[波動関数]]を決定する事を意味する。正の電荷をもつ粒子と負の電荷がそれぞれ陽子と電子だとすればこの系は'''[[水素原子]]'''に相当するが、一般の価数の[[原子核]]を持つ1電子系[[多価イオン]]（'''水素様原子'''）の系も同一の方程式から解を導ける。この方程式は様々な教科書で取り上げられている<ref name="harashima">原島鮮「初等量子力学」裳華房</ref><ref name="shimizu">清水清孝「シュレーディンガー方程式の解き方教えます」共立出版</ref><ref name="kondoh">近藤保、真船文隆「量子化学」裳華房</ref>。

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

== シュレーディンガー方程式 ==
本項の目的は、時間非依存なシュレディンガー方程式{{Equation box 1|cellpadding
|border
|indent=:|equation=<math>\hat{H}\psi(\boldsymbol{x}_0,\boldsymbol{x}_1)= E \psi(\boldsymbol{x}_0,\boldsymbol{x}_1)</math>　　…({{EquationRef|X1}})|border colour=#50C878|background colour=#ECFCF4}}
でハミルトニアンが{{Equation box 1|cellpadding
|border
|indent=:|equation=<math>\hat{H}=-{\hbar^2 \over 2 m_0}\Delta_0 -{\hbar^2 \over 2 m_1} \Delta_1 - {Q \over |\boldsymbol{x}_0-\boldsymbol{x}_1|}</math>  …({{EquationRef|X2}})|border colour=#50C878|background colour=#ECFCF4}}
と書ける場合の厳密解を求める事である。ここで

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

は{{Math|'''R'''{{sup|3}}}}の元であり、{{Math|''m''{{sub|0}}}}、{{Math|''m''{{sub|1}}}}、{{Mvar|Q}}は正の定数であり、
:<math>\Delta_j =
{\partial^2\over \partial x_j{}^2}
+{\partial^2\over \partial y_j{}^2}
+{\partial^2\over \partial z_j{}^2}</math>、
であり、{{Mvar|&hbar;}}は[[ディラック定数|換算プランク定数]]である。

前述した物理的状況においては、２つの粒子の電荷をそれぞれ{{Math|''e''{{sub|1}}, &minus;''e''{{sub|2}}}}とし、[[真空の誘電率]]を{{math|&epsilon;{{sub|0}}}}とすれば、
:<math>Q={e_1e_2\over 4\pi \varepsilon_0}</math>

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

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

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

と[[換算質量]]

:<math>\mu=\frac{m_0m_1}{m_0+m_1}</math>

を使うと、ハミルトニアン({{EquationNote|X2}})は
:<math>\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}|}</math> 

と書ける[[#H13|<sup>H13</sup>]]{{Rp|page=416}}。

このハミルトニアンは
{{Equation box 1|cellpadding
|border
|indent=:|equation=<math>\hat{H}_c=-{\hbar^2 \over 2 (m_0+m_1)}\Delta_{\boldsymbol{c}}</math>　　…({{EquationRef|Y1}})|border colour=#50C878|background colour=#ECFCF4}}

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

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

である[[#H13|<sup>H13</sup>]]{{Rp|page=235}}。したがって後は非自明な部分である({{EquationNote|Y2}})のスペクトルを求めれば良いことになる[[#新井|<sup>新井</sup>]]{{Rp|page=263}}。そこで以下({{EquationNote|Y2}})のみ焦点を当てる。

=== 無次元化 ===
適切な値{{Math|''a''{{sub|0}}}}と定数{{Math|''E''{{sub|a}}}}を選び、長さとエネルギーをそれぞれ{{Math|''a''{{sub|0}}}}、{{Math|''E''{{sub|a}}}}が1となるように座標変換
:<math>(x',y',z')= (x/a_0, y/a_0, z/a_0)</math>
:<math>E'=E/E_a</math>
してやると、({{EquationNote|Y2}})のハミルトニアンに関する時間非依存なシュレディンガー方程式は{{Equation box 1|cellpadding
|border
|indent=:|equation=<math>
-{1 \over 2 }\Delta_{\boldsymbol{x}'}\psi
- {\psi \over |\boldsymbol{x}'|} = E' \psi
</math>　　…({{EquationRef|A1}})|border colour=#50C878|background colour=#ECFCF4}}
と[[無次元量|無次元化]]される[[#SO96|<sup>SO96</sup>]]{{Rp|2.1.1節}}

簡単な計算により{{Math|''a''{{sub|0}}}}、{{Math|''E''{{sub|a}}}}の具体的な値は
:<math>a_0={\hbar^2 \over \mu Q}</math>、<math>E_a={\mu Q^2\over \hbar^2}</math>　　…({{EquationRef|A2}})
である事が分かる。

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

== 求解 ==

本節では({{EquationNote|Y2}})のハミルトニアンを無次元した{{Equation box 1|cellpadding
|border
|indent=:|equation=<math>
\hat{H}_{\boldsymbol{x}'}=
-
{\hbar^2 \over 2 \mu}\Delta_{\boldsymbol{x}'}
- {Q \over |\boldsymbol{x}'|}
</math>  …（{{EquationNote|A1}}、再掲）|border colour=#50C878|background colour=#ECFCF4}}

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

=== 求解の方針 ===
({{EquationNote|A1}})を解く基本的アイデアは、無次元化した座標系{{math|1=(''x''&prime;,''y''&prime;,''z''&prime;) = (''x''/''a''{{sub|0}}, ''y''/''a''{{sub|0}}, ''z''/''a''{{sub|0}})}}を[[球面座標系|球面座標]]{{Math|(''r''&prime;,''&theta;'',''&phi;'')}}に変換するというものだが、直接球面座標を用いると、計算が複雑になる。そこで計算を楽にするため、以下の事実に着目する。

({{EquationNote|A1}})のハミルトニアンは球対称なポテンシャルを持っており、しかもラプラシアンは回転不変である事が知られているので、({{EquationNote|A1}})のハミルトニアンは回転不変である。よって({{EquationNote|A1}})のハミルトニアンは[[軌道角運動量]]演算子<math>(\hat{L}_x,\hat{L}_y,\hat{L}_z)</math>と可換である：
:<math>[\hat{H}_{\boldsymbol{x}'},\hat{L}_x]
=[\hat{H}_{\boldsymbol{x}'},\hat{L}_y]
=[\hat{H}_{\boldsymbol{x}'},\hat{L}_z]
=0</math>  
よって特に、軌道角運動量演算子の自乗{{Math|{{hat|'''''L'''''}}{{sup|2}}}}とも可換である：
:<math>[\hat{H}_{\boldsymbol{x}'},\hat{\boldsymbol{L}^2}]=0</math>  
よって{{Math|{{hat|''H''}}{{sub|'''''x'''''&prime;}}}}は{{Math|{{hat|'''''L'''''}}{{sup|2}}}}と同時対角化できるはずである<ref group="注" name=":0">厳密にいうと、量子力学で扱わねばならない無限次元の線形代数においては、２つの作用素が同時対角化可能であること('''強可換性''')は一般には交換子が0になる事（'''可換性'''）よりも強い条件である[[水素原子におけるシュレーディンガー方程式の解#.E6.96.B0.E4.BA.95|<sup>新井</sup>]]{{Rp|page=179}}。したがって可換性から同時対角化可能性を結論付けるのは本当は正しい推論ではない。したがってここはあくまで、交換子が0になってるため同時対角化可能で「あろう」という推測の元、発見的解法を試みたと解釈すべきである。</ref>、さらに
:<math>[\hat{\boldsymbol{L}^2},\hat{L}_z]=0</math>  
である事から、{{Math|{{hat|''H''}}{{sub|'''''x'''''&prime;}}, {{hat|'''''L'''''}}{{sup|2}}, {{hat|''L''}}{{sub|''z''}}}}の３つを同時対角化できるはずである<ref group="注" name=":0" />。

そこでまず、{{Math|{{hat|'''''L'''''}}{{sup|2}}, {{hat|''L''}}{{sub|''z''}}}}の同時固有関数を求め、これを利用して{{Math|{{hat|''H''}}{{sub|'''''x'''''&prime;}}}}の固有関数を求める。

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

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

</math>　　 …({{EquationRef|W1}})

である。ラプラシアンを[[球面座標系|球面座標]]{{Math|(''r''&prime;,''&theta;'',''&phi;'')}}で書き表し、動径方向と球面方向にわけると、
: <math>\Delta_{\boldsymbol{x}'}={1\over r'^2}(\Delta_{r'}+\Delta_S)</math>　   …({{EquationRef|W2}})
と書ける[[水素原子におけるシュレーディンガー方程式の解#.E6.AD.A6.E8.97.A411-15|<sup>武藤11-15</sup>]]{{Rp|page=13}}。ここで
: <math>\Delta_{r'}=\frac{\partial}{\partial r'}\left(r'^2\frac{\partial}{\partial r'}\right)</math>  、<math>\Delta_S=-{1\over \hbar^2}\hat{\boldsymbol{L}^2}</math>　　　…({{EquationRef|W3}})
であり[[水素原子におけるシュレーディンガー方程式の解#.E6.AD.A6.E8.97.A411-15|<sup>武藤11-15</sup>]]{{Rp|page=13}}、{{Math|{{hat|'''''L'''''}}{{sup|2}}}}は[[軌道角運動量]]演算子の自乗である。({{EquationNote|W1}})のラプラシアンを極座標表示した上で({{EquationNote|W1}})に({{EquationNote|B1}})の波動関数を代入すると、({{EquationNote|B1}})が{{Math|{{hat|'''''L'''''}}{{sup|2}}/''&hbar;''{{sup|2}}}}の固有値{{Math|''&ell;''(''&ell;''+1)}}に対応する固有関数であった事から、
: <math>{1\over 2r'^2}\left(\Delta_{r'}-\ell(\ell+1)
+2r+2r^2E'\right)
 R(r')
P_\ell^{|m|}(\cos \theta) \, e^{im\phi}=0
</math>　

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

束縛状態では{{Mvar|E}}は負の値しか取らないので、記号を簡単にするため
:<math>n:={1\over \sqrt{-2E'}}, \rho:={2r' \over n}</math>　　　…({{EquationRef|W4}})
と定義し[[#原94|<sup>原94</sup>]]{{Rp|page=103-105}}、{{Mvar|R}}を{{Mvar|ρ}}の関数とみなすと、
:<math>\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

</math>　…({{equationref|W5}})
が成立する[[#石川15|<sup>石川15</sup>]]{{Rp|page=154}}。

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

({{EquationNote|W5}})の方程式を解くことで各{{Math|1=''n'' = 0, 1, 2, &hellip;}}に対し、（({{EquationNote|A2}})の{{Math|''E''{{sub|a}}}}を単位とした）エネルギー
:<math>E'_n= -\frac{1}{2n^2}</math>                   ...({{EquationRef|B2}})
に対する解が見つかる[[#新井|<sup>新井</sup>]]{{Rp|page=285-286}}。{{Math|''E'''{{sub|''n''}}}}に対応する固有関数は、

<math>
\begin{cases}
0\le \ell \le n-1\\
|m|\le\ell
\end{cases}
</math>　　　　　　　　　　　　　　　　　　　　　…({{EquationRef|B3}})

に対してのみ存在し、そのときの{{Math|''R''(''r'')}}は[[ラゲールの陪多項式|ラゲールの陪関数]]
:<math>
R_{n,\ell}(\rho)=
 \exp \left( - \rho\right)
\rho^\ell L^{2\ell+1}_{n+\ell}\left(2\rho\right)
</math>　　　　　×規格化定数　　　　　…({{EquationRef|B4}})
に一致する。ここで
:<math>
\rho={r'\over n}
</math>、
:<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)

</math>
である。

=== 規格化定数 ===
3次元空間における体積要素{{Math|1=d''V'' = d''x''&prime;d''y''&prime;d''z''&prime;}}は動径方向の線素{{Math|d''r''&prime;}}と球面方向の面素{{Math|1=d''S'' = sin''&theta;''d''&theta;''d''&phi;''}}を用いて
: <math>\operatorname{d}V
=
r'^2
\operatorname{d}r'
\operatorname{d}S
</math>
と書けるので、({{EquationNote|B1}})における{{Mvar|ψ}}のノルム
: <math>\|\psi\|{}:=
\sqrt{\int_0^{\pi} |\psi(x',y',z')|^2\operatorname{d}V}
</math>　　　
も
: <math>\|\psi\|{}:=\|R\|_r\|Y\|_S

</math>　　                         　…({{EquationRef|M1}})
と「変数分離」する。ここで

: <math>Y(\theta,\varphi)=P_\ell^{|m|}(\cos \theta) \, e^{im\phi}

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

</math>　　　　…({{EquationRef|M3}})
({{EquationNote|M2}})のノルムを1にする規格化定数の値は「[[軌道角運動量]]」の項目に書いてあり、
:<math> (-1)^{(m+|m|)/2}
         \sqrt{ \frac{2\ell+1}{4\pi} \frac{(\ell-|m|)!}{(\ell+|m|)!} \,}
 
</math>
である[[#原94|<sup>原94</sup>]]{{Rp|page=97}}。

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

== 結論 ==
無次元化した({{EquationNote|A1}})をベースにしたこれまでの議論を通常の単位系に戻すことで以下の結論が得られる。
:<math>a_0={\hbar^2 \over \mu Q}</math>
とし、{{Math|''n''>0}}を自然数、{{Math|''&ell;'', ''m''}}を以下を満たす整数とする：

<math>
\begin{cases}
0\le \ell \le n-1\\
|m|\le\ell
\end{cases}
</math>　　　　　　　　　　…('''{{EquationNote|B3}}'''、再掲)　　　　　

このとき({{EquationNote|Y2}})のハミルトニアンはエネルギー
:<math>E_n= -\frac{\mu Q^2}{2\hbar^2n^2}</math>
に対し、
:<math>\hat{\boldsymbol{L}^2}\psi =\hbar^2 \ell(\ell+1)\psi</math>  
:<math>\hat{L}_z\psi =m\hbar \psi</math>  
を満たす固有関数
:<math>
\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^{|m|}(\cos \theta) \, e^{im\phi}
</math> 　×(規格化定数)　　…({{EquationRef|B5}})
を持つ。ここで
: <math>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}
</math>、 
: <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)

</math> 
: <math>
\rho_n={r \over na_0}
</math>
であり、規格化定数は
:<math> (-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)!} }
</math>
である。

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

</math>  s.t. <math>
\sum_{n,\ell,m} a_{n,\ell,m}{}^2 <\infty

</math> 
の形に書けるのは、{{Math|{{hat|''H''}}{{sub|'''''x'''''}}}}の負のスペクトルに対応するベクトルだけで、正のスペクトルに対応するベクトルはこの方法では表記できない[[#H13|<sup>H13{{Rp|page=399}}</sup>]]。

=== 量子数 ===
ハミルトニアン({{EquationNote|Y2}})の固有関数({{EquationNote|B5}})に登場する２つの変数は以下のように呼ばれる：
* {{Mvar|n}}は'''[[主量子数]]'''と呼ばれ、{{Math|{{hat|''H''}}{{sub|'''''x'''''}}}}のエネルギー固有値の大きさを司っている。
* {{Mvar|&ell;}}は[[方位量子数|'''軌道角運動量量子数''']]（方位量子数）と呼ばれ、{{Math|{{hat|'''''L'''''}}{{sup|2}}}}の固有値の大きさを司っている。
* {{Mvar|m}}は'''[[磁気量子数]]'''（軌道磁気量子数）と呼ばれ、{{Math|{{hat|''L''}}{{sub|''z''}}}}の固有値の大きさを司っている。

なお、{{Math|''n''&minus;''&ell;''&minus;1}}は、動径方向の波動関数の節の数を表している{{要出典|date=2017年8月}}。

==== 化学的意味 ====
{{See also|電子配置}}

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

=== リュードベリ定数 ===
{{See also|水素スペクトル系列}}
エネルギー準位が{{Math|''E''{{sub|n}}}}にある電子がエネルギー準位が{{Math|''E''{{sub|n&prime;}}}}に落ちると、
:<math>E_{n}-E_{n'}= \frac{\mu Q^2}{2\hbar^2}\left({1\over n'^2}-{1\over n^2}\right)</math>
のエネルギーが
:<math>E_{n}-E_{n'}= {\hbar c \over \lambda} </math>
を満たす波長{{Mvar|λ}}の光となって放出される。したがって
:<math>{1\over \lambda}= \frac{\mu Q^2}{2\hbar^3c}\left({1\over n'^2}-{1\over n^2}\right)</math>

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

== ({{EquationNote|W5}})の解 ==
本節の目的は、微分方程式({{EquationNote|W5}})を解き、({{EquationNote|B3}})、({{EquationNote|B4}})を導出することである。

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

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

==== {{Mvar|&rho;}}が十分小さい場合の({{EquationNote|W5}})の近似解 ====

({{EquationNote|W5}})における{{Mvar|R}}の係数は{{Mvar|&rho;}}が十分小さいところでは{{Math|''&ell;''(''&ell;''+1)}}と近似できるので、({{EquationNote|W5}})は
:<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>]]{{Rp|page=155}}。

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

{{Indent|<math>R(\rho)=\rho^\ell</math>   ・・・({{equationref|C1}})}}

の形で書ける。

==== {{Mvar|&rho;}}が十分大きい場合の({{EquationNote|W5}})の近似解 ====

({{EquationNote|W5}})式を{{Math|''ρ''{{sup|2}}}}で割った上で{{Math|''&rho;'' → &infin;}}の極限をとることで、{{Mvar|&rho;}}が十分大きいところでは({{EquationNote|W5}})は

{{Indent|<math>\frac{\operatorname{d}^2R}{\operatorname{d}\rho^2} - \frac{1}{4}R=0</math>}}

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

{{Indent|<math>R(\rho)=\mathrm{e}^{- \rho / 2}</math>  ・・・({{equationref|C2}})}}

である。

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

{{Indent|<math>R(\rho) = \rho^\ell u(\rho) \mathrm{e}^{- \rho / 2}</math>　　　…({{equationref|C3}})}}

の形に変数変換する。一般に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>

を満たすので、({{EquationNote|W5}})の第一項、および第二項は、

:<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>　　

である。上式を({{EquationNote|W5}})に代入すると、すべての項に{{Math|e{{sup|&minus;{{sfrac|''&rho;''|2}}}}}}が掛かっていることがわかる。よって各項を{{Math|e{{sup|&minus;{{sfrac|''&rho;''|2}}}}}}で割った上で式を整理して、

{{Indent|<math>
\rho^\ell\frac{\operatorname{d}^2u(\rho)}{\operatorname{d}\rho^2}+\left\{2\rho^{\ell-1}(\ell+1)-\rho^\ell\right\}\frac{\operatorname{d}u(\rho)}{\operatorname{d}\rho}+(n-\ell-1)\rho^{\ell-1}u(\rho) = 0
</math>}}

を得る。この式の両辺を{{Math|''&rho;''{{sup|''&ell;''&minus;1}}}}で割ると、

{{Indent|<math>
\rho\frac{\operatorname{d}^2u(\rho)}{\operatorname{d}\rho^2}+(2\ell+2-\rho)\frac{\operatorname{d}u(\rho)}{\operatorname{d}\rho}+(n-\ell-1)u(\rho) = 0
</math>}}

となる[[#石川15|<sup>石川15</sup>]]{{Rp|page=155}}。こうして得た式(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>     …({{Equationref|C4}})

ラゲールの陪方程式の解{{Math|''u''(''&rho;'')}}は'''[[ラゲールの陪多項式]]'''（ラゲール陪関数）と呼ばれる形の定数倍になることが知られている。ラゲールの陪多項式（ラゲール陪関数）{{Math|''L''{{su|p=''m''|b=''k''}}(''&rho;'')}}は下記のように定義される。
:<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}
</math>

ここで、{{Mvar|k}}は
:<math>
0\le k \le m
</math>   …({{Equationref|C5}})

を満たす整数である。

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

ラゲール陪多項式の係数の条件式({{EquationNote|C4}})から、<math>hoge\ell</math>は
:<math>
0\le \ell \le n-1
</math>   …({{Equationref|C7}})

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

==  規格化定数({{EquationNote|M5}})の導出 ==
規格化定数を{{Mvar|C&prime;}}とすると、規格化条件
:<math>
\int_{0}^{\infty} |R(r')|^2 r'^2 \operatorname{d}r' = 1
</math> 

は、({{EquationNote|W2}})、({{EquationNote|W4}})より、

:<math>
\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}
</math> …({{equationref|D1}})

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

:<math>\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}</math>

ので、後者の式を({{EquationNote|D1}})に対して用いる事で、
: <math>\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)!}

</math>

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

が得られる。

なお、無次元化する前のハミルトニアン({{EquationNote|Y2}})に対する規格化定数は、変数変換
:<math>
\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'
</math>   

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

となる[[#原94|<sup>原94</sup>]]{{Rp|page=106}}。

== 具体的な値 ==

水素原子の波動関数の{{Math|1=''&ell;''  = 0~3}}における角因子は以下のようになる。ここで{{Math|&Theta;(''&theta;'')}}、{{Math|&Phi;(''&phi;'')}}はそれぞれ、動径方向の関数
: <math>Y(\theta,\varphi)=P_\ell^{|m|}(\cos \theta) \, e^{im\varphi}

</math>　
の右辺の積の第一成分と第二成分を規格化したものである。なお、{{Math|&Phi;(''&phi;'')}}の指数関数の虚数部分はオイラーの公式により一対の{{Math|&Phi;(''&phi;'')}}関数の一次結合で書き換えられる。
{| class="wikitable" style="text-align:left; white-space:nowrap;"
|+ 
! {{mvar|&ell;}} || {{mvar|m}} || {{Math|&Phi;(''&phi;'')}} || {{Math|&Theta;(''&theta;'')}} ||  || {{Math|&Phi;(''&phi;'')&Theta;(''&theta;'')}}（極座標） || {{Math|&Phi;(''&phi;'')&Theta;(''&theta;'')}}（直交座標） || 記号
|-
| 0 || 0 || <math>\frac{1}{\sqrt{2\pi}}</math> || <math>\frac{1}{\sqrt{2}}</math> ||  || <math>\frac{1}{2\sqrt{\pi}}</math> || <math>\frac{1}{2\sqrt{\pi}}</math> || <math>\mbox{s} \,</math>
|-
| 1 || 0 || <math>\frac{1}{\sqrt{2\pi}}</math> || <math>\sqrt{\frac{3}{2}}\cos\theta</math> ||  || <math>\frac{1}{2}\sqrt{\frac{3}{\pi}}\cos\theta</math> || <math>\frac{1}{2}\sqrt{\frac{3}{\pi}}\frac{z}{r}</math> || <math>\mbox{p}_{z} \,</math>
|-
| 1 || +1 || <math>\frac{1}{\sqrt{2\pi}}\exp(i\phi)</math> || <math>\frac{\sqrt{3}}{2}\sin\theta</math> || rowspan="2" | <math>\Bigg\{</math> || <math>\frac{1}{2}\sqrt{\frac{3}{\pi}}\sin\theta\cos\phi</math> || <math>\frac{1}{2}\sqrt{\frac{3}{\pi}}\frac{x}{r}</math> || <math>\mbox{p}_{x} \,</math>
|-
| 1 || -1 || <math>\frac{1}{\sqrt{2\pi}}\exp(-i\phi)</math> || <math>\frac{\sqrt{3}}{2}\sin\theta</math> || <math>\frac{1}{2}\sqrt{\frac{3}{\pi}}\sin\theta\sin\phi</math> || <math>\frac{1}{2}\sqrt{\frac{3}{\pi}}\frac{y}{r}</math> || <math>\mbox{p}_{y} \,</math>
|-
| 2 || 0 || <math>\frac{1}{\sqrt{2\pi}}</math> || <math>\frac{1}{2}\sqrt{\frac{5}{2}}(3\cos^2\theta-1)</math> ||  || <math>\frac{1}{4}\sqrt{\frac{5}{\pi}}(3\cos^2\theta-1)</math> || <math>\frac{1}{4}\sqrt{\frac{5}{\pi}}\frac{2z^2-x^2-y^2}{r^2}</math> || <math>\mbox{d}_{3z^2-r^2}</math>
|-
| 2 || +1 || <math>\frac{1}{\sqrt{2\pi}}\exp(i\phi)</math> || <math>\frac{\sqrt{15}}{2}\sin\theta\cos\theta</math> || rowspan="2" | <math>\Bigg\{</math> || <math>\frac{1}{2}\sqrt{\frac{15}{\pi}}\sin\theta\cos\theta\cos\phi</math> || <math>\frac{1}{2}\sqrt{\frac{15}{\pi}}\frac{zx}{r^2}</math> || <math>\mbox{d}_{zx} \,</math>
|-
| 2 || -1 || <math>\frac{1}{\sqrt{2\pi}}\exp(-i\phi)</math> || <math>\frac{\sqrt{15}}{2}\sin\theta\cos\theta</math> || <math>\frac{1}{2}\sqrt{\frac{15}{\pi}}\sin\theta\cos\theta\sin\phi</math> || <math>\frac{1}{2}\sqrt{\frac{15}{\pi}}\frac{yz}{r^2}</math> || <math>\mbox{d}_{yz} \,</math>
|-
| 2 || +2 || <math>\frac{1}{\sqrt{2\pi}}\exp(2i\phi)</math> || <math>\frac{\sqrt{15}}{4}\sin^2\theta</math> || rowspan="2" | <math>\Bigg\{</math> || <math>\frac{1}{4}\sqrt{\frac{15}{\pi}}\sin^2\theta\cos2\phi</math> || <math>\frac{1}{4}\sqrt{\frac{15}{\pi}}\frac{x^2-y^2}{r^2}</math> || <math>\mbox{d}_{x^2-y^2}</math>
|-
| 2 || -2 || <math>\frac{1}{\sqrt{2\pi}}\exp(-2i\phi)</math> || <math>\frac{\sqrt{15}}{4}\sin^2\theta</math> || <math>\frac{1}{4}\sqrt{\frac{15}{\pi}}\sin^2\theta\sin2\phi</math> || <math>\frac{1}{2}\sqrt{\frac{15}{\pi}}\frac{xy}{r^2}</math> || <math>\mbox{d}_{xy} \,</math>
|-
| 3 || 0 || <math>\frac{1}{\sqrt{2\pi}}</math> || <math>\frac{1}{2}\sqrt{\frac{7}{2}}(5\cos^3\theta-3\cos\theta)</math> ||  || <math>\frac{1}{4}\sqrt{\frac{7}{\pi}}(5\cos^3\theta-3\cos\theta)</math> || <math>\frac{1}{4}\sqrt{\frac{7}{\pi}}\frac{z(2z^2-3x^2-3y^2)}{r^3}</math> || <math>\mbox{f}_{z(5z^2-3r^2)}</math>
|-
| 3 || +1 || <math>\frac{1}{\sqrt{2\pi}}\exp(i\phi)</math> || <math>\frac{1}{4}\sqrt{\frac{21}{2}}(5\cos^2\theta-1)\sin\theta</math> || rowspan="2" | <math>\Bigg\{</math> || <math>\frac{1}{4}\sqrt{\frac{21}{2\pi}}(5\cos^2\theta-1)\sin\theta\cos\phi</math> || <math>\frac{1}{4}\sqrt{\frac{21}{2\pi}}\frac{x(5z^2-r^2)}{r^3}</math> || <math>\mbox{f}_{x(5z^2-r^2)}</math>
|-
| 3 || -1 || <math>\frac{1}{\sqrt{2\pi}}\exp(-i\phi)</math> || <math>\frac{1}{4}\sqrt{\frac{21}{2}}(5\cos^2\theta-1)\sin\theta</math> || <math>\frac{1}{4}\sqrt{\frac{21}{2\pi}}(5\cos^2\theta-1)\sin\theta\sin\phi</math> || <math>\frac{1}{4}\sqrt{\frac{21}{2\pi}}\frac{y(5z^2-r^2)}{r^3}</math> || <math>\mbox{f}_{y(5z^2-r^2)}</math>
|-
| 3 || +2 || <math>\frac{1}{\sqrt{2\pi}}\exp(2i\phi)</math> || <math>\frac{\sqrt{105}}{4}\cos\theta\sin^2\theta</math> || rowspan="2" | <math>\Bigg\{</math> || <math>\frac{1}{4}\sqrt{\frac{105}{\pi}}\cos\theta\sin^2\theta\cos2\phi</math> || <math>\frac{1}{4}\sqrt{\frac{105}{\pi}}\frac{z(x^2-y^2)}{r^3}</math> || <math>\mbox{f}_{z(x^2-y^2)}</math>
|-
| 3 || -2 || <math>\frac{1}{\sqrt{2\pi}}\exp(-2i\phi)</math> || <math>\frac{\sqrt{105}}{4}\cos\theta\sin^2\theta</math> || <math>\frac{1}{4}\sqrt{\frac{105}{\pi}}\cos\theta\sin^2\theta\sin2\phi</math> || <math>\frac{1}{2}\sqrt{\frac{105}{\pi}}\frac{xyz}{r^3}</math> || <math>\mbox{f}_{xyz} \,</math>
|-
| 3 || +3 || <math>\frac{1}{\sqrt{2\pi}}\exp(3i\phi)</math> || <math>\frac{1}{4}\sqrt{\frac{35}{2}}\sin^3\theta</math> || rowspan="2" | <math>\Bigg\{</math> || <math>\frac{1}{4}\sqrt{\frac{35}{2\pi}}\sin^3\theta\cos3\phi</math> || <math>\frac{1}{4}\sqrt{\frac{35}{2\pi}}\frac{x(x^2-3y^2)}{r^3}</math> || <math>\mbox{f}_{x(x^2-3y^2)}</math>
|-
| 3 || -3 || <math>\frac{1}{\sqrt{2\pi}}\exp(-3i\phi)</math> || <math>\frac{1}{4}\sqrt{\frac{35}{2}}\sin^3\theta</math> || <math>\frac{1}{4}\sqrt{\frac{35}{2\pi}}\sin^3\theta\sin3\phi</math> || <math>\frac{1}{4}\sqrt{\frac{35}{2\pi}}\frac{y(3x^2-y^2)}{r^3}</math> || <math>\mbox{f}_{y(3x^2-y^2)}</math>
|}

[[原子番号]]{{Mvar|Z}}の水素様原子の動径関数は以下のようになる。
:<math>\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}</math>

{|
|-
| [[ファイル:R-1s-orbital.png|class=skin-invert-image|200px]] ||  ||  || 
|-
| 1s軌道の動径関数 ||  ||  || 
|-
| [[ファイル:R-2s-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R-2p-orbital.png|class=skin-invert-image|200px]] ||  || 
|-
| 2s軌道の動径関数 || 2p軌道の動径関数 ||  || 
|-
| [[ファイル:R-3s-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R-3p-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R-3d-orbital.png|class=skin-invert-image|200px]] || 
|-
| 3s軌道の動径関数 || 3p軌道の動径関数 || 3d軌道の動径関数 || 
|-
| [[ファイル:R-4s-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R-4p-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R-4d-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R-4f-orbital.png|class=skin-invert-image|200px]]
|-
| 4s軌道の動径関数 || 4p軌道の動径関数 || 4d軌道の動径関数 || 4f軌道の動径関数
|}

動径関数を2乗し{{Math|''r''{{sup|2}}}}を掛けた動径分布{{Math|''r''{{sup|2}}''R''(''r''){{sup|2}}}}は、核の中心からのある距離における電子の存在確率に相当する。
{|
|-
| [[ファイル:R2r2-1s-orbital.png|class=skin-invert-image|200px]] ||  ||  || 
|-
| 1s軌道の動径分布 ||  ||  || 
|-
| [[ファイル:R2r2-2s-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R2r2-2p-orbital.png|class=skin-invert-image|200px]] ||  || 
|-
| 2s軌道の動径分布 || 2p軌道の動径分布 ||  || 
|-
| [[ファイル:R2r2-3s-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R2r2-3p-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R2r2-3d-orbital.png|class=skin-invert-image|200px]] || 
|-
| 3s軌道の動径分布 || 3p軌道の動径分布 || 3d軌道の動径分布 || 
|-
| [[ファイル:R2r2-4s-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R2r2-4p-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R2r2-4d-orbital.png|class=skin-invert-image|200px]] || [[ファイル:R2r2-4f-orbital.png|class=skin-invert-image|200px]]
|-
| 4s軌道の動径分布 || 4p軌道の動径分布 || 4d軌道の動径分布 || 4f軌道の動径分布
|}

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

== 脚注 ==
=== 注釈 ===
<references group="注"/>
=== 出典 ===
{{reflist}}

== 参考文献 ==
* 書籍
** [新井97] {{Cite book|ref=A97|author=[[新井朝雄]]|title=ヒルベルト空間と量子力学|series=共立講座２１正規の数学１６|date=1997/1/25|year=|publisher=[[共立出版]]}}
** [原94] {{Cite book|和書|ref=原94|author=原康夫|authorlink=原康夫|title=５ 量子力学|series=岩波基礎物理シリーズ|date=1994/6/6|year=1994|publisher=[[岩波書店]]|ISBN=978-4000079259}}
** [H13] {{Cite book|ref=H13|author=Brian C.Hall|title=Quantum Theory for Mathematicians|series=Graduate Texts in Mathematics 267|date=2013/7/1|year=|publisher=[[Springer]]}}
** [SO96] {{Cite book|ref=SO96|author=Attila Szabo, Neil S. Ostlund|title=Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory|series=Dover Books on Chemistry|date=1996/7/2|publisher=Dover Publications|ISBN=978-0486691862}}
*** 邦訳：{{Cite book|author=A. ザボ, N.S. オストランド|title=新しい量子化学―電子構造の理論入門〈上〉、〈下〉|date=1996/7/2|year=1987/07|publisher=[[東京大学出版会]]|translator=大野公男, 望月祐志, 阪井健男}}
* レクチャーノート
** [武藤11-15]{{Cite web|和書|url=http://www.th.phys.titech.ac.jp/~muto/lectures/QMII11/QMII11_chap15.pdf|title=第１５章	中心力ポテンシャルでの束縛状態|accessdate=2017年8月13日|author=武藤一雄|format=pdf|work=[http://www.th.phys.titech.ac.jp/~muto/index.html 量子力学第二 平成２３年度　学部　５学期 ]|publisher=[[東京工業大学]]|ref=武藤11-15}}
** [石川15] {{Cite web|和書|url=https://www.sci.hokudai.ac.jp/~ishikawa/ryoushirikigaku4.pdf|title=量子力学|accessdate=2017年8月13日|author=石川健三|format=pdf|publisher=[[北海道大学]]理学部|ref=石川15|date=2015年1月21日}}

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

== 外部リンク ==
* [http://keisan.casio.jp/exec/system/1176445657 水素原子の電子分布の計算]

{{DEFAULTSORT:すいそけんしにおけるしゆれていんかほうていしきのかい}} 
[[Category:計算化学]]
[[Category:量子力学]]
[[Category:微分方程式]]
[[Category:水素|しゆれていんかほうていしきのかい]]
[[Category:電子軌道]]