水素原子におけるシュレーディンガー方程式の解（すいそげんしにおけるシュレーディンガーほうていしきのかい）は、陽子1個と電子1個がクーロン力により結合している状態に対して、シュレーディンガー方程式を解くことによって求めた電子の波動関数である。陽子の代わりに、一般的な価数の原子核に置き換えた1電子系多価イオンは水素様原子と呼ばれており、本項の手法がそのまま適用できる。水素原子は最も単純な原子モデルであるため、水素原子を用いると波動関数や固有エネルギー、電子殻や化学結合など化学の基礎となる重要な要素の理論的背景の説明が比較的容易になる。このため、水素原子におけるシュレディンガー方程式は様々な教科書で取り上げられている。なお、微細構造、超微細構造、ラムシフトなどの効果は、電子-陽子間のクーロン相互作用に比べると極めて小さく、本項では取り扱わない。時間に依存しないシュレーディンガー方程式は以下の式で与えられる。水素原子のハミルトニアンは、∇を極座標におけるラプラシアンとすると、次のようになる。ここで、ポテンシャルV(r)は点電荷のつくるスカラーポテンシャルを用いることにより以下の式となる。また、μは換算質量であり、水素原子核の質量をM、電子の質量をmとすると、であるが、M»mであるからμ=mと近似でき、結局(1.3)式はとなる。よって、水素原子におけるシュレーディンガー方程式は、(1.1)式および及び(1.6)式より、すなわちとなる。(1.8)式を解くために、ψ(r,θ,φ)=R(r)Y(θ,φ)と書けると仮定して(1.8)式に代入すると次式を得る。Y(θ,φ)は球面調和関数という。両辺をR(r)Y(θ,φ)で割ると、となり、さらに両辺にformula_1を掛けて整理すると、となる。左辺はrのみの関数、右辺はθ、φのみの関数であるので、両辺が恒等的に等しくなるためには両辺が定数である必要がある。よって、とおける。(2.4)式は動径方程式、(2.5)式は角度方程式と呼ばれる。また、(2.5)式の両辺に-Y(θ,φ)sinθを掛けて整理すると、となり、さらにY(θ,φ)=Θ(θ)Φ(φ)と書けると仮定して(2.6)式に代入し、両辺をΘ(θ)Φ(φ)で割って整理すると、となる。左辺はθのみの関数、右辺はφのみの関数であるので、とおける。整理すると、水素原子のシュレディンガー方程式は変数分離の結果、下記3つの微分方程式(2.4, 2.9, 2.10)に帰着された。ただし、ψ(r,θ,φ)=R(r)Θ(θ)Φ(φ)である。(2.10)式を変形すると、であるので、その解はすぐに求まり、である（Aは定数）。また、波動関数はいたるところで連続であるので、Φ(φ+2π)=Φ(φ)が境界条件となり、すなわちが必要である。(3.4)式は、オイラーの公式より、と書けるので、m=0,±1,±2,… を得る。したがって、Φ(φ)は、となる。このmを磁気量子数（軌道磁気量子数）といい、量子化軸方向の軌道の角運動量の成分を指定する。より、であるから、規格化因子Aはとなる。ただし、ここで A は正の実数にとった。まずは下記の(2.9)を変形するところから始める。両辺にformula_2を掛けてこれを変形すると、と書ける。次に(4.1)をz=cosθ、Θ(θ)=P(z)として変数変換する。であるから、式変形するとを得る。(4.2)式において β=l(l+1) とおく。lを軌道角運動量量子数（方位量子数）と呼び、全軌道角運動量を指定する。この(4.3)式の解P(z)はルジャンドル陪関数と呼ばれ、極座標で微分方程式を解く際によく現れる形の関数である。次に、磁気量子数 m=0 とおく。この(4.4)式の形はルジャンドルの微分方程式として知られており、l=0,1,2...のときその解はルジャンドルの多項式(下記の式4.5)となる。また、ルジャンドル陪関数はとなることが知られている。よって、Θ(θ)=P(z)より、である（Aは規格化因子）。ルジャンドル陪関数は下記の関係を満たすことが知られている。よって、とするためには、sqrt{frac{(l-|m|)!}{(l+|m|)!}}P^{m}_{l}(z)\となることが解の条件である。規格化因子を求めると、解は下記のようになる。sqrt{l+frac{1}{2}}sqrt{frac{(l-|m|)!}{(l+|m|)!}}P^{m}_{l}(cos heta)ただし、formula_3(ルジャンドル多項式)である。最後に、動径方向の関数R(r)を解く。式(2.4)を下記のような流れで変形する。βをβ=l(l+1)で置き換えて移項すると、両辺にformula_4を掛けて、ここで、formula_5（ここでnは主量子数と呼ばれる）として変数変換する。ボーア半径aを用いれば、formula_6と書ける。式(5.2)の微分の項はであるから、式(5.3)を(5.2)に代入してここで、formula_7として変数変換する。ハートリーEを用いれば、formula_8と書ける。(5.5)を整理して、式(5.6)の両辺をformula_9で割って式(5.7)の解R(ρ)がラゲールの陪多項式（ラゲール陪関数）を用いて書けることを示す。両辺にformula_9をかけて、ρが十分小さいところではこの形の方程式はオイラーの微分方程式の解法に準ずる方法で解ける。その解はの形で書ける。また、ρが十分大きいところにおいて(5.7)はとなるので、formula_11とおいて(6.4)に代入するととなる。formula_12では発散する不適切な解となるので、である。こうして、ρが小さいところではformula_13、大きいところではformula_14となることがわかった。formula_15が十分に小さい場合・大きい場合のどちらにもあてはまらない場合、と書き、解R(ρ)をこれら3つの場合(式6.3, 6.5, 6.6)の積として下記のように表せると仮定する。この(6.7)を(5.7)に代入すると、を得る。3つの関数の積の微分はformula_16であるから、formula_17　　(6.9)(6.9)を(6.8)に代入すると、すべての項にformula_18が掛かっていることがわかる。よって各項をformula_18で割るとこれを整理して、(6.11)の両辺をformula_20で割ると、こうして得た式(6.12)は下記の式(6.13)に示したラゲールの陪方程式（ラゲール陪関数）の形になっている。ラゲールの陪方程式の解formula_21はラゲールの陪多項式（ラゲール陪関数）と呼ばれる形の定数倍になることが知られている。ラゲールの陪多項式（ラゲール陪関数）formula_22は下記のように定義される。よって式(6.13)の解はとなる。これを式(6.7)に代入してを得る。ラゲールの陪多項式（ラゲール陪関数）は下記の直交性を満たすことが知られている。よって、およびより(7.2)および(7.4)より、式(7.5)を規格化定数ｃについて解くと、解は下記のようになる。ただし、:formula_23はラゲールの陪多項式である。ψ(r,θ,φ)=R(r)Θ(θ)Φ(φ) として変数分離した関数のそれぞれの解はであるから、水素原子の波動関数ψ(r,θ,φ)は下記のようになる。整理して、を得る。(7.7)式から、formula_24の場合は、formula_25しか許されないことがわかるが、これはformula_26とならなければならないためである。ちなみに、formula_27は、動径方向の波動関数の節の数を表している。また、(4.8)式から、formula_25ならば、formula_29とならなければならない。この電子状態をs軌道（1s軌道）と呼ぶ。また、formula_30の場合には、同様に(7.7)式から、formula_31かformula_25でなければならないことがわかる。formula_31のとき、(4.8)式から、formula_34 , formula_29 , formula_36のいずれかとなる。この電子状態をp軌道（2p軌道）という。なお、formula_25ならばformula_38となり、この電子状態を2s軌道という。さらに、formula_39のとき、formula_40 , formula_31 , formula_25のいずれかとなるが、特に、formula_40のとき、formula_44 , formula_34 , formula_29 , formula_36 , formula_48のいずれかとなり、この電子状態をd軌道（3d軌道）という。一般に、formula_49のとき、1個のs軌道、3個のp軌道、5個のd軌道、7個のf軌道…のformula_50種類の電子軌道があり、全部でformula_51個の電子状態がある。水素原子において、s軌道,p軌道,d軌道,f軌道…のエネルギー準位は縮退している。これはエネルギー固有値が、formula_52となり、lやmに依存しないためである。なお、水素原子に磁場をかけると、これらのエネルギー準位は、スピン部分を無視して考えた場合、磁気量子数mの違いにより分裂する(→ゼーマン効果)。電場をかけた場合も、シュタルク効果によって分裂する。このとき、異なるlの軌道同士の線形結合をとった混成軌道がハミルトニアンの固有状態となる。水素原子の波動関数のl=0～3における角因子は以下のようになる。ここでformula_53の指数関数の虚数部分はオイラーの公式により一対のformula_53関数の一次結合で書き換えられる。原子番号Zの水素様原子の動径関数は以下のようになる。left(frac{Z}{a_0} ight)^{3/2}left(2-frac{Zr}{a_0} ight)exp(-Zr/2a_0)formula_55formula_56formula_57formula_58formula_59formula_60formula_61formula_62動径関数を2乗しrの2乗を掛けた動径分布は、核の中心からのある距離における電子の存在確率に相当する。詳しくは電子配置の項を参照のこと。

出典:wikipedia