これらの導出はウェルトンの量子光学[1]による。
真空の電磁場のゆらぎは原子核の電位ポテンシャルに揺らぎを与え、電子の位置に揺らぎを与える。この揺らぎが準位のずれを引き起こす。電子の位置エネルギーの差は以下の式で表される。
![{\displaystyle \Delta V=V({\vec {r}}+\delta {\vec {r}})-V({\vec {r}})=\delta {\vec {r}}\cdot \nabla V({\vec {r}})+{\frac {1}{2}}(\delta {\vec {r}}\cdot \nabla )^{2}V({\vec {r}})+\cdots }](//wikimedia.org/api/rest_v1/media/math/render/svg/ff0033ab8c73508c8a171583e129342eea75d0ba)
揺らぎは等方的なので以下が成り立つ。
![{\displaystyle {\begin{aligned}\langle \delta {\vec {r}}\rangle _{\mathrm {vac} }&=0,\\\langle (\delta {\vec {r}}\cdot \nabla )^{2}\rangle _{\mathrm {vac} }&={\frac {1}{3}}\langle (\delta {\vec {r}})^{2}\rangle _{\mathrm {vac} }\nabla ^{2}.\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/116ff53470d36604fce99fde296729e13b6bcbf0)
よって
![{\displaystyle \langle \Delta V\rangle ={\frac {1}{6}}\langle (\delta {\vec {r}})^{2}\rangle _{\mathrm {vac} }\left\langle \nabla ^{2}\left({\frac {-e^{2}}{4\pi \epsilon _{0}r}}\right)\right\rangle _{\mathrm {at} }.}](//wikimedia.org/api/rest_v1/media/math/render/svg/c65f916fbcce34fd75cd0be5aac4e57834d8f1da)
いま、波数 k→ 周波数νの電場によりずれ(δr)k→ が生じるとする。このとき電子の運動方程式は
![{\displaystyle m{\frac {d^{2}}{dt^{2}}}(\delta r)_{\vec {k}}=-eE_{\vec {k}},}](//wikimedia.org/api/rest_v1/media/math/render/svg/5d00770fe2a1a1233f4a14ad34292258536979bc)
これは周波数ν がボーア軌道の周波数ν0 よりも大きいときのみ成り立つ。このため
。電場の揺らぎの周波数が軌道周波数よりも小さい場合、電子は電場に対して反応することができない。
ν で振動する電場に対しては
![{\displaystyle \delta r(t)\cong \delta r(0)e^{-i\nu t}+{\text{c.c.}},}](//wikimedia.org/api/rest_v1/media/math/render/svg/c9d8fc9986206bd996a764c84384e59da5d9e0c3)
であるため、
![{\displaystyle (\delta r)_{\vec {k}}\cong {\frac {e}{mc^{2}k^{2}}}E_{\vec {k}}={\frac {e}{mc^{2}k^{2}}}{\mathcal {E}}_{\vec {k}}\left(a_{\vec {k}}e^{-i\nu t+i{\vec {k}}\cdot {\vec {r}}}+h.c.\right)\qquad {\text{with}}\qquad {\mathcal {E}}_{\vec {k}}=\left({\frac {\hbar ck/2}{\epsilon _{0}\Omega }}\right)^{1/2},}](//wikimedia.org/api/rest_v1/media/math/render/svg/83a2a909ca99fcb1d0a700bad3c49aa7d149ff71)
ここで
は繰り込みに用いる体積 (水素原子を包む仮想的な箱の体積)である。すべての
について和をとると
![{\displaystyle {\begin{aligned}\langle (\delta {\vec {r}})^{2}\rangle _{\mathrm {vac} }&=\sum _{\vec {k}}\left({\frac {e}{mc^{2}k^{2}}}\right)^{2}\left\langle 0\left|(E_{\vec {k}})^{2}\right|0\right\rangle \\&=\sum _{\vec {k}}\left({\frac {e}{mc^{2}k^{2}}}\right)^{2}\left({\frac {\hbar ck}{2\epsilon _{0}\Omega }}\right)\\&=2{\frac {\Omega }{(2\pi )^{3}}}4\pi \int dkk^{2}\left({\frac {e}{mc^{2}k^{2}}}\right)^{2}\left({\frac {\hbar ck}{2\epsilon _{0}\Omega }}\right)&&{\text{since continuity of }}{\vec {k}}{\text{ implies }}\sum _{\vec {k}}\to 2{\frac {\Omega }{(2\pi )^{3}}}\int d^{3}k\\&={\frac {1}{2\epsilon _{0}\pi ^{2}}}\left({\frac {e^{2}}{\hbar c}}\right)\left({\frac {\hbar }{mc}}\right)^{2}\int {\frac {dk}{k}}\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/6fe2366ace730f6e9f49952d43226bc0c4d46ba4)
この積分は周波数の上限と下限を定めない限り発散する。下限は上に述べた
から
と求まる。また上限はコンプトン波長とし、
と求まる。これらの制限から積分の収束値が求まる。
.
原子軌道とクーロン場から
![{\displaystyle \left\langle \nabla ^{2}\left({\frac {-e^{2}}{4\pi \epsilon _{0}r}}\right)\right\rangle _{at}={\frac {-e^{2}}{4\pi \epsilon _{0}}}\int d{\vec {r}}\psi ^{*}({\vec {r}})\nabla ^{2}\left({\frac {1}{r}}\right)\psi ({\vec {r}})={\frac {e^{2}}{\epsilon _{0}}}|\psi (0)|^{2},}](//wikimedia.org/api/rest_v1/media/math/render/svg/337968c8dcdad059c778a900b544fa0b14cbe6a1)
また以下の等式により
![{\displaystyle \nabla ^{2}\left({\frac {1}{r}}\right)=-4\pi \delta ({\vec {r}}).}](//wikimedia.org/api/rest_v1/media/math/render/svg/9d1b0799bfb63998d2347db6c26bceab8583394b)
p 軌道では、波動関数は原点でゼロになる。このためエネルギーのずれは起こらない。しかし、s軌道は原点において
![{\displaystyle \psi _{2S}(0)={\frac {1}{(8\pi a_{0}^{3})^{1/2}}},}](//wikimedia.org/api/rest_v1/media/math/render/svg/0473920f6b152dcc1fb13246c911ef45b65f1111)
という値をもつ。ここで、ボーア半径
![{\displaystyle a_{0}={\frac {4\pi \epsilon _{0}\hbar ^{2}}{me^{2}}}.}](//wikimedia.org/api/rest_v1/media/math/render/svg/7e2b9af7128430d1ea9de5de73478b7a57c9986d)
を用いた。このため
.
したがって、位置エネルギーの差は
![{\displaystyle \langle \Delta V\rangle ={\frac {4}{3}}{\frac {e^{2}}{4\pi \epsilon _{0}}}{\frac {e^{2}}{4\pi \epsilon _{0}\hbar c}}\left({\frac {\hbar }{mc}}\right)^{2}{\frac {1}{8\pi a_{0}^{3}}}\ln {\frac {4\epsilon _{0}\hbar c}{e^{2}}}=\alpha ^{5}mc^{2}{\frac {1}{6\pi }}\ln {\frac {1}{\pi \alpha }},}](//wikimedia.org/api/rest_v1/media/math/render/svg/16ae6264925a0323ae2ecf120a9f2f8b8112334e)
ここで
は微細構造定数である。これによる周波数のシフトは1 GHz となる、実験で観測されたシフトと一致する。
ウェルトンによるラムシフトの導出はツィッターベヴェーグンクを用いたダーウィン項の導出との類似点がある[2]:80–81。