author Frank Zirkelbach Mon, 18 Jun 2012 12:15:16 +0000 (14:15 +0200) committer Frank Zirkelbach Mon, 18 Jun 2012 12:15:16 +0000 (14:15 +0200)

index e4c2362..c4cf869 100644 (file)
@@ -169,12 +169,34 @@ HSC potential \ldots

\subsubsection{Fully separable form of the pseudopotential}

+KB transformation \ldots

+\subsection{Spin-orbit interaction}

-\subsection{Spin orbit interaction}
-
-
-\subsubsection{Perturbative treatment}
-
-\subsubsection{Non-perturbative method}
-
+Relativistic effects can be incorporated in the normconserving pseudopotential method up to but not including order $\alpha^2$ with $\alpha$ being the fine structure constant.
+With the solutions of the all-electron Dirac equations, the new pseudopotential reads
+\begin{equation}
+V(r)=\sum_{l,m}\left[
+\ket{l+\frac{1}{2},m+{\frac{1}{2}}}V_{l,l+\frac{1}{2}}(r)
+\bra{l+\frac{1}{2},m+{\frac{1}{2}}} +
+\ket{l-\frac{1}{2},m-{\frac{1}{2}}}V_{l,l-\frac{1}{2}}(r)
+\bra{l-\frac{1}{2},m-{\frac{1}{2}}}
+\right] \text{ .}
+\end{equation}
+By defining an averaged potential weighted by the different $j$ degeneracies of the $\ket{l\pm\frac{1}{2}}$ states
+\begin{equation}
+\bar{V}_l(r)=\frac{1}{2l+1}\left(
+l V_{l,l-\frac{1}{2}}(r)+(l+1)V_{l,l+\frac{1}{2}}(r)\right)
+\end{equation}
+and a potential describing the difference in the potential with respect to the spin
+\begin{equation}
+V^{\text{SO}}_l(r)=\frac{2}{2l+1}\left(
+V_{l,l+\frac{1}{2}}(r)-V_{l,l-\frac{1}{2}}(r)\right)
+\end{equation}
+the total potential can be expressed as
+\begin{equation}
+V(r)=\sum_l \ket{l}\left[\bar{V}_l(r)+V^{\text{SO}}_l(r)LS\right]\bra{l}
+\text{ ,}
+\end{equation}
+where the first term correpsonds to the mass velocity and Darwin relativistic corrections and the latter is associated with the spin-orbit (SO) coupling.