In the review by Oliviera et. al., I do not understand how the authors get equation 14: \begin{equation} E[\rho]=\sum_i^Mn_i\langle\Psi_i|-\frac{1}{2}\nabla^2+v_{ext}(\vec{r})+\frac{1}{2}\int\frac{\rho(\vec{r'})}{|\vec{r}-\vec{r'}|}|\Psi_i\rangle+E_{xc}+\frac{1}{2}\sum_\beta^N\sum_{\beta\neq\alpha}^N\frac{Z_\alpha Z_\beta}{|\vec{R_\alpha}-\vec{R_\beta}|}\tag{14.} \end{equation} To me the energy should be in this form:
$$E = \langle\psi|T + V_\mathrm{ext} + V_\mathrm{ee}|\psi\rangle$$
How did they manage to put some of term out of the bracket? I also have trouble retrieving (16) as they cite, are they using a Taylor expansion?
References:
- Oliveira, A. F.; Seifert, G.; Heine, T.; Duarte, H. A. Density-functional based tight-binding: an approximate DFT method. J. Braz. Chem. Soc. 2009, 20 (7), 1193–1205. DOI: 10.1590/S0103-50532009000700002.