We provide a straightforward and numerically efficient procedure to perform local-density-$\text{approximation}+\text{Hubbard}$ I approximation $\left(\text{LDA}+\text{HIA}\right)$ calculations, including self-consistency over the charge density, within the full potential linearized augmented plane-wave (FP-LAPW) method. This implementation is all-electron, includes spin-orbit interaction, and makes no shape approximations for the charge density. The method is applied to calculate selected heavy actinides in the paramagnetic phase. The electronic structure and spectral properties of Am and Cm metals obtained are in agreement with previous dynamical mean-field theory $\left(\text{LDA}+\text{DMFT}\right)$ calculations and with available experimental data. We point out that the charge-density self-consistent $\text{LDA}+\text{HIA}$ calculations predict the $f$ charge on Bk to exceed the atomic integer $f^8$ value by 0.22.

• Received 12 March 2009

DOI:https://doi.org/10.1103/PhysRevB.80.085106