\parindent=0pt
\pagestyle{empty}
\begin{document}
\section{Gamma equation}
\begin{equation}
\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt
\end{equation}
\begin{equation}
\int_0^\infty t^{z-1}e^{-\alpha t}\,dt=\frac{\Gamma(z)}{\alpha^z}
\end{equation}
\section{spherical harmonics}
\begin{eqnarray}
Y_{0,0}&=&\sqrt{\frac1{4\pi}}\\
Y_{1,0}&=&\sqrt{\frac3{4\pi}}\cos \theta\\
Y_{1,\pm1}&=&\mp\sqrt{\frac3{8\pi}}\sin \theta e^{\pm i\varphi}\\
Y_{2,0}&=&\sqrt{\frac5{16\pi}}(3\cos^2 \theta -1)\\
Y_{2,\pm1}&=&\mp\sqrt{\frac{15}{8\pi}}\sin \theta\cos \theta e^{\pm i\varphi}\\
Y_{2,\pm2}&=&\mp\sqrt{\frac{15}{32\pi}}\sin^2 \theta e^{\pm i2\varphi}
\end{eqnarray}
\section{H\_atom}
\subsection{$n=1, l=0$\hspace{1zw}(1s)}
\begin{equation}
\psi_\mathrm{1s}=2e^{-r}Y_{0,0}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{1s}|\psi_\mathrm{1s}\right>&=&4\int e^{-2r}|Y_{0,0}|^2\,dV\\
&=&4\int r^2e^{-2r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&4\frac{\Gamma(3)}{2^3}\\
&=&1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-2\int e^{-r}Y_{0,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)e^{-r}Y_{0,0}\,dV\\
&=&-2\int\left(r^2-2r\right)e^{-2r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-2\left[\frac{\Gamma(3)}{2^3}-2\frac{\Gamma(2)}{2^2}\right]\\
&=&-2\left(\frac14-\frac12\right)\\
&=&\frac12
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac1{\pi}\int e^{-r}\frac1{r}e^{-r}\,dV\\
&=&-4\int re^{-2r}\,dr\int\!\!\!\int |Y_{0,0}|^2\,d\theta d\varphi\\
&=&-4\frac{\Gamma(2)}{2^2}\\
&=&-1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&\frac12-1\\
&=&-\frac12
\end{eqnarray*}
\newpage
\subsection{$n=2,l=0$\hspace{1zw}(2s)}
\begin{equation}
\psi_\mathrm{2s}=\frac1{\sqrt 8}(2-r)e^{-\frac{r}2}Y_{0,0}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2s}|\psi_\mathrm{2s}\right>&=&\frac18\int (2-r)^2e^{-r}|Y_{0,0}|^2\,dV\\
&=&\frac18\int (r^4-4r^3+4r^2)e^{-r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&\frac18\left[\Gamma(5)-4\Gamma(4)+4\Gamma(3)\right]\\
&=&1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac1{16}\int (2-r)e^{-\frac{r}2}Y_{0,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)(2-r)e^{-\frac{r}2}Y_{0,0}\,dV\\
&=&-\frac1{16}\int (2-r)\left(-\frac14r^3+\frac52r^2-4r\right)e^{-r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin\theta\,d\theta d\varphi\\
&=&-\frac1{16}\int \left(\frac14 r^4-3 r^3+9r^2-8r\right)e^{-r}\,dr\\
&=&-\frac1{16}\left[\frac14\Gamma(5)-3\Gamma(4)+9\Gamma(3)-8\Gamma(2)\right]\\
&=&-\frac{6-18+18-8}{16}\\
&=&\frac18
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac18\int \frac1{r}(2-r)^2e^{-r}|Y_{0,0}|^2\,dV\\
&=&-\frac18\int (r^3-4r^2+4r)e^{-r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac18\left[\Gamma(4)-4\Gamma(3)+4\Gamma(2)\right]\\
&=&-\frac{6-8+4}{8}\\
&=&-\frac14
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}
\newpage
\subsection{$n=2,l=1$\hspace{1zw}(2p)}
\begin{equation}
\psi_\mathrm{2p}=\frac1{\sqrt{24}}re^{-\frac{r}2}Y_{1,m}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2p}|\psi_\mathrm{2p}\right>&=&\frac1{24}\int r^2e^{-r}|Y_{1,m}|^2,dV\\
&=&\frac1{24}\int r^4e^{-r}\,dr\int\!\!\!\int |Y_{1,m}|^2\,d\theta d\varphi\\
&=&\frac1{24}\Gamma(5)\\
&=&\frac1{24}\times 24\\
&=&1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac1{48}\int re^{-\frac{r}2}Y_{1,m}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)re^{-\frac{r}2}Y_{1,m}\,dV\\
&=&-\frac1{48}\int \left(\frac14r^4-2r^3\right)e^{-r}\,dr \int\!\!\!\int|Y_{1,m}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{48}\left[\frac14\Gamma(5)-2\Gamma(4)\right]\\
&=&-\frac{6-12}{48}\\
&=&\frac1{8}
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac1{24}\int re^{-r}|Y_{1,m}|^2\,dV\\
&=&-\frac1{24}\int r^3e^{-r}\,dr\int\!\!\!\int |Y_{1,m}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{24}\Gamma(4)\\
&=&-\frac14
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}
\newpage
\subsection{$n=3,l=0$\hspace{1zw}(3s)}
\begin{equation}
\psi_\mathrm{3s}=\frac2{\sqrt{6561\times 3}}\left(27-18r+2r^2\right)e^{-\frac{r}3}Y_{0,0}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3s}|\psi_\mathrm{3s}\right>&=&\frac4{6561\times 3}\int \left(27-18r+2r^2\right)^2e^{-\frac{2r}3}|Y_{0,0}|^2\,dV\\
&=&\frac4{6561\times 3}\int \left(729-972r+432r^2-72r^3+4r^4\right)r^2 e^{-\frac{2r}3}\,dr\int\!\!\!\int |Y_{0,0}|^2\,d\theta d\varphi\\
&=&\frac4{6561\times 3}\left[729\Gamma(3)\left(\frac32\right)^3-972\Gamma(4)\left(\frac32\right)^4+432\Gamma(5)\left(\frac32\right)^5-72\Gamma(6)\left(\frac32\right)^6+4\Gamma(7)\left(\frac32\right)^7\right]\\
&=&\frac4{3^9}\left(\frac{3^9}4-\frac{3^{10}}2+3^9\times 4-3^9\times 5+\frac52\times 3^9\right)\\
&=&1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac2{6561\times 3}\int \left(27-18r+2r^2\right)e^{-\frac{r}3}Y_{0,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)\left(27-18r+2r^2\right)e^{-\frac{r}3}Y_{0,0}\,dV\\
&=&-\frac2{6561\times 3}\int \left(27-18r+2r^2\right)\left(\frac29r^2-6r+39-\frac{54}{r}\right)e^{-\frac23r} r^2\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac2{6561\times 3}\int \left(\frac49r^6-16r^5+192r^4-4\times3^5r^3+25\times 3^4r^2-2\times 3^6r\right)e^{-\frac23r}\,dr\\
&=&-\frac2{6561\times 3}\left[\frac49\Gamma(7)\left(\frac32\right)^7-16\Gamma(6)\left(\frac32\right)^6+192\Gamma(5)\left(\frac32\right)^5\right.\\
&&{}-\left.4\times3^5\Gamma(4)\left(\frac32\right)^4+25\times 3^4\Gamma(3)\left(\frac32\right)^3-2\times 3^6\Gamma(2)\left(\frac32\right)^2\right]\\
&=&-\frac2{6561\times 3}\left(\frac52\times 3^7-10\times 3^7+16\times 3^7-\frac12\times 3^{10}+\frac{25}4\times 3^7-\frac12\times 3^8\right)\\
&=&\frac2{6561\times 3}\times\frac14\times 3^7\\
&=&\frac1{18}
\end{eqnarray*}
\newpage
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac4{6561\times 3}\int \frac1{r}\left(27-18r+2r^2\right)^2e^{-\frac23r}|Y_{0,0}|^2\,dV\\
&=&-\frac4{6561\times 3}\int \left(729r-972r^2+432r^3-72r^4+4r^5\right)e^{-\frac{2r}3}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&\frac4{6561\times 3}\left[729\Gamma(2)\left(\frac32\right)^2-972\Gamma(3)\left(\frac32\right)^3+432\Gamma(4)\left(\frac32\right)^4-72\Gamma(5)\left(\frac32\right)^5+4\Gamma(6)\left(\frac32\right)^6\right]\\
&=&\frac4{3^9}\left(\frac{3^8}4-3^8+2\times 3^8-2\times 3^8+\frac52\times 3^7\right)\\
&=&-\frac19
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}
\newpage
\subsection{$n=3,l=1$\hspace{1zw}{3p}}
\begin{equation}
\psi_\mathrm{3p}=\frac4{\sqrt{6561\times 6}}r\left(6-r\right)e^{-\frac{r}3}Y_{1,m}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3p}|\psi_\mathrm{3p}\right>&=&\frac{16}{6561\times 6}\int r^2\left(6-r\right)^2e^{-\frac23r}|Y_{1,m}|^2\,dV\\
&=&\frac{16}{6561\times 6}\int (r^6-12r^5+36r^4)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,m}|^2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 6}\left[\Gamma(7)\left(\frac32\right)^7-12\Gamma(6)\left(\frac32\right)^6+36\Gamma(5)\left(\frac32\right)^5\right]\\
&=&\frac{16}{6561\times 6}\left(\frac58\times 3^9-\frac52\times 3^8+3^8\right)\\
&=&1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac8{6561\times 6}\int r\left(6-r\right)e^{-\frac{r}3}Y_{1,m}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r\left(6-r\right)e^{-\frac{r}3}Y_{1,m}\,dV\\
&=&-\frac8{6561\times 6}\int (6r-r^2)\left(-\frac19r^4+\frac83r^3-12r^2\right)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,m}|^2\sin \theta \,d\theta\,d\varphi\\
&=&-\frac8{6561\times 6}\int\left(\frac19r^6-\frac{10}3r^5+28r^4-72r^3\right)e^{-\frac23r}\,dr\\
&=&-\frac8{6561\times 6}\left[\frac19\Gamma(7)\left(\frac32\right)^7-\frac{10}3\Gamma(6)\left(\frac32\right)^6+28\Gamma(5)\left(\frac32\right)^5-72\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac8{6561\times 6}\left(\frac58\times 3^7-\frac{25}4\times 3^6+7\times 3^6-3^7\right)\\
&=&-\frac1{18}
\end{eqnarray*}
\newpage
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac{16}{6561\times 6}\int r^3(6-r)^2e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,m}|^2\,d\theta\,d\varphi\\
&=&-\frac{16}{6561\times 6}\int (r^5-12r^4+36r^3)e^{-\frac23r}\,dr\\
&=&-\frac{16}{6561\times 6}\left[\Gamma(6)\left(\frac32\right)^6-12\Gamma(5)\left(\frac32\right)^5+36\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac{16}{6561\times 6}\left(\frac58\times 3^7-3^7+\frac12\times 3^7 \right)\\
&=&-\frac19
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}
\newpage
\subsection{$n=3,l=2$\hspace{1zw}(3d)}
\begin{equation}
\psi_\mathrm{3d}=\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}Y_{2,m}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d}|\psi_\mathrm{3d}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}|Y_{2,m}|^2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,m}|^2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac8{6561\times 30}\int r^2e^{-\frac{r}3}Y_{2,m}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}Y_{2,m}\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,m}|\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}
\newpage
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}|Y_{2,m}|^2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,m}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}
\subsubsection{$\left
\begin{eqnarray*}
\left
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}
\end{document}
