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

本項、水素原子におけるシュレーディンガー方程式の解（すいそげんしにおけるシュレーディンガーほうていしきのかい）では、ハミルトニアンが
と書ける二粒子系の時間非依存なシュレーディンガー方程式の厳密解を解く（式中の記号の意味は後述）。
物理学的にはこれは、- 質量の正の電荷をもつ粒子と質量が負の電荷を持つ粒子がクーロン力により結合している状況において- 外力は働いておらず、- 相対論的効果を考えない量子力学の範囲内で、- 時間に依存しない定常状態の
粒子の波動関数を決定する事を意味する。正の電荷をもつ粒子と負の電荷がそれぞれ陽子と電子だとすればこの系は水素原子に相当するが、一般の価数の原子核を持つ1電子系多価イオン（水素様原子）の系も同一の方程式から解を導ける。この方程式は様々な教科書で取り上げられている。
なお、微細構造、超微細構造、ラムシフトなどの効果は、いずれも相対論的な量子力学を必要とする為、本項の対象外である。

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


はの元であり、、、Qは正の定数であり、
、
であり、&hbar;は換算プランク定数である。
前述した物理的状況においては、２つの粒子の電荷をそれぞれとし、真空の誘電率をとすれば、

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

**重心系への還元**
(X1)、(X2)により定義される方程式は、重心系に書き直す事により、より簡単な式に還元できる。２つの粒子の重心

と２つの粒子の位置の差

と換算質量

を使うと、ハミルトニアン(X2)は
 
と書けるH13page=416。
このハミルトニアンは
の和である。

(Y1)のハミルトニアンはよく知られた自由粒子のハミルトニアンであり、その連続スペクトルは

であり、点スペクトルは

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

**無次元化**
適切な値と定数を選び、長さとエネルギーをそれぞれ、が1となるように座標変換


してやると、(Y2)のハミルトニアンに関する時間非依存なシュレディンガー方程式は
と無次元化されるSO962.1.1節
簡単な計算により、の具体的な値は
、　　…(A2)
である事が分かる。

**ボーア半径・ハートリー**
特に、陽子の質量が電子の質量より遥かに重いと仮定した場合の水素原子の系における、は


より、
、
である。ここでeは電気素量である。この場合のをボーア半径といい、を基準としたエネルギーの単位をハートリーというSO962.1.1節。

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

**求解の方針**
(A1)を解く基本的アイデアは、無次元化した座標系を球面座標(r′,θ,φ)に変換するというものだが、直接球面座標を用いると、計算が複雑になる。そこで計算を楽にするため、以下の事実に着目する。
(A1)のハミルトニアンは球対称なポテンシャルを持っており、しかもラプラシアンは回転不変である事が知られているので、(A1)のハミルトニアンは回転不変である。よって(A1)のハミルトニアンは軌道角運動量演算子と可換である：
  
よって特に、軌道角運動量演算子の自乗とも可換である：
  
よってはと同時対角化できるはずである、さらに
  
である事から、の３つを同時対角化できるはずである。
そこでまず、の同時固有関数を求め、これを利用しての固有関数を求める。

**との同時固有関数**
との同時固有関数の求め方は「軌道角運動量」の項目に書いてあるので、結論だけを言えば、1=&ell; = 0, 1, 2, …, m = 0, &pm;1, &pm;2, … &pm;&ell;に対し、
  
  
を満たす固有関数ψが存在し、ψは極座標で
 　　×(規格化定数)　…(B1)
という形で書ける。ここでR(r)は任意の自乗可積分関数であり、はルジャンドルの陪多項式
  
である新井page=277。

**R(r′)の決定**
後はR(r′)を決定するだけである。R(r′)を決定するには(B1)を(Y2)のハミルトニアンに入れてシュレディンガー方程式(Y1)を解けば良い。（A1）を式変形すると、
　　 …(W1)
である。ラプラシアンを球面座標(r′,θ,φ)で書き表し、動径方向と球面方向にわけると、
 　   …(W2)
と書ける武藤11-15page=13。ここで
   、　　　…(W3)
であり武藤11-15page=13、は軌道角運動量演算子の自乗である。(W1)のラプラシアンを極座標表示した上で(W1)に(B1)の波動関数を代入すると、(B1)がの固有値&ell;(&ell;+1)に対応する固有関数であった事から、
 　
すなわち
 　
束縛状態ではEは負の値しか取らないので、記号を簡単にするため
　　　…(W4)
と定義し原94page=103-105、Rをρの関数とみなすと、
　…(W5)
が成立する石川15page=154。
この方程式を解くのは複雑な計算を必要とするので後の章にまわし、ここでは結論のみを述べる。
(W5)の方程式を解くことで各1=n = 0, 1, 2, …に対し、（(A2)のを単位とした）エネルギー
                   ...(B2)
に対する解が見つかる新井page=285-286。に対応する固有関数は、
　　　　　　　　　　　　　　　　　　　　　…(B3)
に対してのみ存在し、そのときのR(r)はラゲールの陪関数
　　　　　×規格化定数　　　　　…(B4)
に一致する。ここで
、

である。

**規格化定数**
3次元空間における体積要素1=dV = dx′dy′dz′は動径方向の線素dr′と球面方向の面素1=dS = sinθdθdφを用いて
 
と書けるので、(B1)におけるψのノルム
 　　　
も
 　　                         　…(M1)
と「変数分離」する。ここで
 　　
であり、
 　　　　…(M2)
 　　　　…(M3)
(M2)のノルムを1にする規格化定数の値は「軌道角運動量」の項目に書いてあり、

である原94page=97。
(M3)のノルムを1にする規格化定数の値の計算は後述するが、結論から言えば規格化定数は
　　　　　...(M5)
である。

**結論**
無次元化した(A1)をベースにしたこれまでの議論を通常の単位系に戻すことで以下の結論が得られる。

とし、n>0を自然数、&ell;, mを以下を満たす整数とする：
　　　　　　　　　　…(B3、再掲)　　　　　
このとき(Y2)のハミルトニアンはエネルギー

に対し、
  
  
を満たす固有関数
 　×(規格化定数)　　…(B5)
を持つ。ここで
 、 
  
 
であり、規格化定数は

である。
以上では変数分離により発見的に解を求めたため、(B3)、(B5)に書いたものが解である事は間違いないものの、それ以外に解があるかどうかは不明である。しかし実はこれ以外に解がない事が知られているH13page=399-400。連続スペクトルに相当する部分は、物理的にいえば水素原子がイオン化している状態であり、したがって電子が陽子から逃れていってしまっている#H13|H13page=399。なお、固有関数の和
  s.t.  
の形に書けるのは、の負のスペクトルに対応するベクトルだけで、正のスペクトルに対応するベクトルはこの方法では表記できない#H13|H13page=399。

**量子数**
ハミルトニアン(Y2)の固有関数(B5)に登場する２つの変数は以下のように呼ばれる：- nは主量子数と呼ばれ、のエネルギー固有値の大きさを司っている。- &ell;は軌道角運動量量子数（方位量子数）と呼ばれ、の固有値の大きさを司っている。- mは磁気量子数（軌道磁気量子数）と呼ばれ、の固有値の大きさを司っている。
なお、n−&ell;−1は、動径方向の波動関数の節の数を表しているdate=2017年8月。

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

**リュードベリ定数**
エネルギー準位がにある電子がエネルギー準位がに落ちると、

のエネルギーが

を満たす波長λの光となって放出される。したがって

水素原子の場合、すなわち

の場合の上式右辺の定数、もしくはその定数に対して近似

を行ったときの値をリュードベリ定数という。

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

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

**ρが十分小さい場合の(W5)の近似解**
(W5)におけるRの係数はρが十分小さいところでは&ell;(&ell;+1)と近似できるので、(W5)は

と近似できる石川15page=155。
この形の方程式はオイラーの微分方程式の解法に準ずる方法で解ける。その解は
の形で書ける。

**ρが十分大きい場合の(W5)の近似解**
(W5)式をで割った上でρ → ∞の極限をとることで、ρが十分大きいところでは(W5)は
となる事がわかる。簡単な計算から上記の方程式の一般解は
、
もしくはこれらの線形和である。は発散する不適切な解となるので、
である。

**(W5)からのラゲールの陪多項式の導出**
(C1)、(C2)を参考に、(W5)の厳密解R(ρ)を
の形に変数変換する。一般に3つの関数の積の微分は公式

を満たすので、(W5)の第一項、および第二項は、
　　
である。上式を(W5)に代入すると、すべての項にが掛かっていることがわかる。よって各項をで割った上で式を整理して、
を得る。この式の両辺をで割ると、
となる石川15page=155。こうして得た式(6.12)は下記の式(6.13)に示したラゲールの陪方程式（ラゲール陪関数）の形になっている。
     …(C4)
ラゲールの陪方程式の解u(ρ)はラゲールの陪多項式（ラゲール陪関数）と呼ばれる形の定数倍になることが知られている。ラゲールの陪多項式（ラゲール陪関数）は下記のように定義される。

ここで、kは
   …(C5)
を満たす整数である。
よって(C4)の解は
   ×(規格化定数)
となる。これを変数変換の式(C3)に代入して
         ×(規格化定数)   …(C6)
を得る。
ラゲール陪多項式の係数の条件式(C4)から、は
   …(C7)
を満たす整数でなければならない。

** 規格化定数(M5)の導出**
規格化定数をC′とすると、規格化条件
 
は、(W2)、(W4)より、
 …(D1)
ラゲールの陪多項式（ラゲール陪関数）は下記の直交性を満たすことが知られている

ので、後者の式を(D1)に対して用いる事で、
 
これが(D1)の左辺である1と等しいことから、規格化定数C′について解く事で
　　…(D2)
が得られる。
なお、無次元化する前のハミルトニアン(Y2)に対する規格化定数は、変数変換
   
の分だけ(D2)のものとはずれるので、(Y2)に対する規格化定数は、
　　…(D3)
となる原94page=106。

**具体的な値**
水素原子の波動関数の1=&ell;  = 0~3における角因子は以下のようになる。ここでΘ(θ)、Φ(φ)はそれぞれ、動径方向の関数
 　
の右辺の積の第一成分と第二成分を規格化したものである。なお、Φ(φ)の指数関数の虚数部分はオイラーの公式により一対のΦ(φ)関数の一次結合で書き換えられる。
原子番号Zの水素様原子の動径関数は以下のようになる。

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

**脚注**

**注釈**

**出典**

**参考文献**
 書籍
** 新井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 ISBN=978-0486691862
*** 邦訳：Cite book|author=A. ザボ, N.S. オストランド|title=新しい量子化学―電子構造の理論入門〈上〉、〈下〉|date=1996/7/2|year=1987/07|publisher=東京大学出版会|translator=大野公男, 望月祐志, 阪井健男- レクチャーノート
** 武藤11-15Cite web|和書|url=http://www.th.phys.titech.ac.jp/~muto/lectures/QMII11/QMII11_chap15.pdf|title=第１５章	中心力ポテンシャルでの束縛状態|accessdate=2017年8月13日|author=武藤一雄|format=pdf|work= 量子力学第二 平成２３年度　学部　５学期 |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日

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

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