---
abstract: 'We have constructed a LDA+DMFT method in the framework of the iterative perturbation theory (IPT) with full LDA Hamiltonian without mapping onto the effective Wannier orbitals. We then apply this LDA+DMFT method to ferromagnetic bcc-Fe and fcc-Ni as a test of transition metal, and to antiferromagnetic NiO as an example of transition metal oxide. In Fe and Ni, the width of occupied 3d bands is narrower than those in LDA and Ni $6~{\mbox{eV}}$ satellite appears. In NiO, the resultant electronic structure is of charge-transfer insulator type and the band gap is $4.3~{\mbox{eV}}$. These results are in good agreement with the experimental XPS. The configuration mixing and dynamical correlation effects play a crucial role in these results.'
author:
- 'Oki Miura$^{1,2}$'
- 'Takeo Fujiwara$^{1,2,3}$'
title: ' Electronic structure and effects of dynamical electron correlation in ferromagnetic bcc-Fe, fcc-Ni and antiferromagnetic NiO '
---

Introduction
============

Much attention has been focused on strongly correlated electron systems, since these systems show anomalous physical properties such as various spin, charge and orbital order, metal-insulator transition and so on. The local density approximation (LDA) based on the density functional theory (DFT) is hardly applicable to these fruitful physical properties in strongly correlated electron systems. LDA+U method could not treat metallic phase near metal-insulator transition since it has been developed in order to discuss magnetic insulators. [@re:LDA+U2] Fluctuation of charge or/and spin densities plays an important role in low-energy excitations near the metal-insulator transition. To discuss realistic materials of strongly correlated electron systems, another method is needed for the “dynamical” electron correlation.

The dynamical mean field theory (DMFT) has been developed and applied to model systems of strongly correlated electrons, which led us to a unified picture of low- and high-energy excitations in anomalous metallic phase near metal-insulator transition. [@re:Rev-DMFT] DMFT is based on mapping of many electron systems in bulk onto single impurity atom embedded in effective medium, namely the single impurity problem. In this mapping procedure, the on-site dynamical correlation is included and the inter-atomic correlation is neglected. The combination of DMFT with LDA, called LDA+DMFT, has been developed in order to discuss realistic systems. [@re:Ans; @re:Lic] LDA or DFT meet their own new stage with the GW approximation, [@re:GW-Arya] which is based on the many body perturbation theory and the random phase approximation (RPA). The combination of DMFT with the GW approximation (GW+DMFT) was proposed to include both on-site and inter-atomic Coulomb interaction. [@re:GW]

To solve the mapped single impurity problem within DMFT, one can use several computational schemes such as the quantum Monte Carlo method (QMC), [@re:Rozen-QMC; @re:Geo-QMC; @re:QMC-difficult; @re:PQMC; @re:CTQMC] the iterative perturbation theory (IPT), [@re:Kaj-IPT; @re:Fuji; @re:miura; @re:MOIPT] the non-crossing approximation (NCA) [@re:NCA] and the exact diagonalization (ED). [@re:Geo-ED]

Hirsch-Fye QMC [@re:Hirsch-Fye-QMC] is widely used as a solver for the mapped single impurity problem within DMFT. However, it is not applicable in the low temperature limit. Moreover, Hirsch-Fye QMC has serious difficulty in application of multi-orbital systems with spin-flip and pair-hopping terms of the exchange interactions, since one cannot apply the Hubbard-Stratonovich transformation to these systems. [@re:QMC-difficult] To go beyond such problems of Hirsch-Fye QMC, several QMC methods have been proposed. Projective QMC [@re:PQMC] has been developed to calculate in the low temperature limit and Continuous Time QMC (CTQMC) [@re:CTQMC] to calculate multi-orbital systems with spin-flip and pair-hopping terms.

In spite of development of those QMC, the exact calculations such as QMC and ED can only be suitable for simple Hamiltonians with relatively small size because of their computational costs. Several LDA+DMFT methods with QMC adopt the projected effective Hamiltonian of Wannier-like functions with a few adopted bands for reducing computational costs. With the use of the effective Hamiltonian of Wannier-like functions in LDA+DMFT, only a few adopted bands are discussed with fixing the hybridization mixing and could not describe a possible change of hybridization due to local Coulomb interaction.

To carry out calculations in realistic materials with multiple-orbitals, approximate calculation schemes of IPT or NCA could be more suitable because of its efficient CPU-time, though IPT is not applicable to cases of large Coulomb interactions and NCA cannot yield the Fermi liquid behavior at low energies and in low temperature limit. [@re:NCA-error]

The IPT method was developed by Kajueter and Kotliar for non-degenerate orbital and the results show good agreements with the ED results. [@re:Kaj-IPT] However, it was reported that the quasiparticle peak in the IPT results tends to be sharper [@re:Rev-DMFT] and that the IPT results in La$_{1-x}$Sr$_x$TiO$_3$ does not reprduce quasiparticle peak at 1000K. [@re:Nekrasov-LaSrTiO3]

We generalized IPT for multi-orbital bands and applied this method to doubly degenerated ${\mbox{e}_{g}}$ bands and triply degenerated ${\mbox{t}_{2g}}$ bands on simple cubic and body-centered cubic lattices. [@re:Fuji; @re:miura] We also verified that this DMFT with IPT scheme is fairly applicable to arbitrary electron occupation cases with different Coulomb interaction $U$. The spectrum shows electron affinity and ionization levels with different electron configurations of different occupation numbers and the system becomes insulating state at an integer filling in sufficiently large $U$. [@re:Fuji; @re:miura] IPT for multi orbitals was also proposed by Laad, Craco and Müller-Hartmann and applied to several realistic materials. [@re:MOIPT]

In this paper, we apply this DMFT with IPT scheme to realistic materials, where we adopt full LDA Hamiltonian without reducing its size and IPT as a solver for the mapped single impurity problem. In the following, we notify this DMFT with IPT as “LDA+DMFT-IPT(1)”. The goal of the present paper is as follows:

\(i) To generalize LDA+DMFT-IPT(1) which is applicable to various realistic strongly correlated materials, both metallic and insulating, multi-atom (compound), spin-polarized and strongly hybridized cases between s, p and d-bands.

\(ii) To apply LDA+DMFT-IPT(1) to ferromagnetic bcc-Fe and fcc-Ni and to antiferromagnetic NiO. For ferromagnetic bcc-Fe and fcc-Ni, we will discuss whether LDA+DMFT-IPT(1) reproduces the accurate width of occupied 3d bands and the observed satellite of spectrum at $6~{\mbox{eV}}$ below the Fermi energy in Ni (“$6~{\mbox{eV}}$ satellite”). Application of LDA+DMFT-IPT(1) to ferromagnetic bcc-Fe and fcc-Ni is a test in transition metals since previous LDA+DMFT [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT] has been applied to those metals and reproduced those physical properties well. For antiferromagnetic NiO, we will discuss whether LDA+DMFT-IPT(1) reproduces the accurate band gap and the correct description of the charge-transfer insulator. These physical properties are not well reproduced in the previous LSDA calculation.

\(iii) To verify the applicability of LDA+DMFT-IPT(1) to various realistic strongly correlated materials. Actually, the physical properties mentioned in (ii) are reproduced fairly well by LDA+DMFT-IPT(1). We will discuss the origin of drastic changes of those physical properties and show the validity of LDA+DMFT-IPT(1) by comparing with other theoretical scheme including other LDA+DMFT method [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT; @re:Kunes-NiO-DMFT] and experiments. We will discuss the origin of nickel $6~{\mbox{eV}}$ satellite, which has been concluded as hole-hole scattering process, [@re:Ni-satellite] by using the spectrum of Ligand Field Theory (LFT) obtained in the procedure of LDA+DMFT-IPT(1). In bcc-Fe and fcc-Ni case, the narrowing of occupied 3d bands in LDA+DMFT-IPT(1) occurs due to the on-site dynamical electron correlation. In NiO case, the appearance of charge-transfer insulator and accurate band-gap in LDA+DMFT-IPT(1) is due to the change of hybridization between Ni-3d and O-2p bands caused by local Coulomb interaction of nickel 3d bands.

In Sec. \[sec:FORM\], the general formulation of LDA+DMFT-IPT(1) will be given. In Sec. \[sec:result-Fe-Ni\], we test LDA+DMFT-IPT(1) in systems of ferromagnetic bcc-Fe and fcc-Ni. In the energy spectra, one can see nickel $6~{\mbox{eV}}$ satellite caused by hole-hole scattering process and the narrowing of occupied 3d bands caused by on-site dynamical electron correlation. Section \[sec:result-NiO\] will be devoted to calculation of the electronic structures in antiferromagnetic NiO by LDA+DMFT-IPT(1). Local Coulomb interaction of nickel 3d bands enhances the hybridization between Ni-3d and O-2p bands and system becomes the charge-transfer insulator with accurate band gap. In both sections, we will present the energy spectra and ${\bm k}$-resolved spectrum and LFT spectra of single isolated ion. Section \[sec:CONC\] is the summary. We discuss the applicability of LDA+DMFT-IPT(1) to various realistic strongly correlated materials by comparing with other LDA+DMFT method and experiments.

Formulation {#sec:FORM}
===========

Hamiltonian and Coulomb interactions
------------------------------------

We proceed with the Hubbard-type Hamiltonian as follows : $$\begin{aligned}
\! \! \! H &=& H_{{\mbox{\scriptsize{LDA}}}}^{dc} + H_{int}~,
\label{eqn:ham}\\
\! \! \!  H_{{\mbox{\scriptsize{LDA}}}}^{dc} \! \! 
&=& \! \! \! \! \! \! \! \! \! \! \! \sum_{j(m{\sigma})j'(m')} \! \! \! \! \! \! \! 
(h^{jj'({\mbox{\scriptsize{LDA}}})}_{mm'{\sigma}}
  -\Delta h^{{\sigma}}_{dc}\delta_{jj'}\delta_{mm'}) c_{jm{\sigma}}^{\dag}c_{j'm'{\sigma}}~,
\label{eqn:HLDAdc}\\
\! \! \!  H_{int} &=& \! \! {\frac{1}{2}} \! \! \sum_{j(m{\sigma})}\! \! 
U_{m_1m_2:m_3m_4}c_{jm_1{\sigma}_1}^{\dag}c_{jm_2{\sigma}_2}^{\dag}c_{jm_4{\sigma}_4}c_{jm_3{\sigma}_3}~. 
\nonumber \\
\label{eqn:int}\end{aligned}$$ The Hamiltonian $H$ is divided into two parts: the unperturbed part ($H^{dc}_{{\mbox{\scriptsize{LDA}}}}$) and the interaction part ($H_{int}$). $h^{jj'({\mbox{\scriptsize{LDA}}})}_{mm'{\sigma}}$ in Eq. (\[eqn:HLDAdc\]) is the full LDA Hamiltonian constructed by the Tight-Binding LMTO (TB-LMTO) method, [@re:LMTO] without projecting onto any kind of effective local orbitals. The use of original full LMTO Hamiltonian in LDA+DMFT-IPT(1) enables us to describe the change of hybridization caused by local Coulomb interaction. $\Delta h^{{\sigma}}_{dc}$ in Eq. (\[eqn:HLDAdc\]) is the double-counting term included in LDA Hamiltonian as averaged Coulomb and exchange interactions, because we need to subtract $\Delta h^{{\sigma}}_{dc}$ from $h^{jj'({\mbox{\scriptsize{LDA}}})}_{mm'{\sigma}}$ in order to construct the unperturbed part of the Hamiltonian $H^{dc}_{{\mbox{\scriptsize{LDA}}}}$. Note that an index $j$ in the summation of Eqs. (\[eqn:HLDAdc\]) and (\[eqn:int\]) runs over atomic sites, $\{m\}$ orbital indices and $\{{\sigma}\}$ spins.

$U_{m_1m_2:m_3m_4}$ in Eq. (\[eqn:int\]) is the on-site electron-electron interaction matrix and presented by using the Slater integrals $F^k$ as $$\begin{aligned}
U_{m_1m_2:m_3m_4} &=& \sum_{0\leq k\leq 2l}a_k(m_1m_2:m_3m_4)F^k ,
 \label{eqn:U}\\
a_k(m_1m_2:m_3m_4) &=& {\frac{4\pi}{2k+1}} \sum_{q=-k}^{k} \langle lm_1|Y_{kq}|lm_2\rangle \nonumber\\
&\times & \langle lm_3|Y_{kq}^{*}|lm_4\rangle,
 \label{eqn:ak}\end{aligned}$$ where $\{|lm\rangle\}$ are the basis set of complex spherical harmonics. Note that Eqs. (\[eqn:U\]) and (\[eqn:ak\]) are the general expression and invariant under any rotational operation. [@re:LDA+U2] For d-orbitals we only need $F^0, F^2$ and $F^4$. The averaged Coulomb and exchange parameters $U$ and $J$ for 3d-orbitals are defined in terms of the Slater integrals $F^0, F^2$ and $F^4$. Here we regard these equations as a definition of the Slater integrals $F^0, F^2$ and $F^4$ in terms of $U$ and $J$: $$\begin{aligned}
      U &=& F^0,  \label{eqn:F0}\\
      J &=& \frac{1}{14}(F^2+F^4),  \label{eqn:F2pF4} \end{aligned}$$ with a ratio of $F^2$ and $F^4$ for 3d orbitals as [@re:F2F4-ratio1] $$\begin{aligned}
F^4/F^2 &\sim& 0.625.  \label{eqn:F2F4}\end{aligned}$$ Therefore, once we obtain experimentally or assume the values $U$ and $J$ of 3d transition metal elements, the values of the Slater integrals $F^0$, $F^2$ and $F^4$ can be evaluated by using Eqs. (\[eqn:F0\]), (\[eqn:F2pF4\]) and (\[eqn:F2F4\]) and then each matrix element of the on-site electron-electron interaction $U_{m_1m_2:m_3m_4}$ can be evaluated by using Eq. (\[eqn:U\]) and (\[eqn:ak\]).

As discussed above, we treat the interaction Hamiltonian Eq. (\[eqn:int\]) for full 3d-orbitals on transition metals within the framework of DMFT. Thus, the present Hamiltonian in Eq. (\[eqn:int\]) does work for the cases with partially filled $\mbox{e}_{g}$ and $\mbox{t}_{2g}$-orbitals of ferromagnetic bcc-Fe, fcc-Ni and antiferromagnetic NiO. On the other hand, the case of triply degenerated ${\mbox{t}_{2g}}$-orbitals and doubly degenerated Hubbard model were treated within the framework of DMFT in Refs.  and , respectively.

Dynamical mean field theory
---------------------------

The matrices of the lattice Green’s function and the local Green’s function are presented on the basis of the non-orthogonal local base set as $$\begin{aligned}
&& [G({\bm k}, {i{\omega}_n})]^{-1}
\nonumber\\
&&\ \ \ \ =[({i{\omega}_n}+\mu)O({\bm k})-\{H^{dc}_{{\mbox{\scriptsize{LDA}}}}({\bm k})+{\Sigma}({i{\omega}_n})\}],
\label{eqn:Gkw}\\
&& G({i{\omega}_n})= {\frac{1}{V}}\int d{\bm k}~G({\bm k},{i{\omega}_n}) ,
\label{eqn:Gw}\end{aligned}$$ where $H^{dc}_{{\mbox{\scriptsize{LDA}}}}({\bm k})$, $O({\bm k})$ and $\mu$ are the Hamiltonian matrix of Eq. (\[eqn:HLDAdc\]), the overlap matrix, both in the ${\bm k}$-space, and the chemical potential, respectively. Here we neglect the ${\bm k}$-dependence of the self-energy ${\Sigma}$ within the framework of DMFT. The chemical potential $\mu$ is determined to satisfy Luttinger’s theorem. [@re:Mul] All the matrices in Eqs. (\[eqn:Gkw\]) and (\[eqn:Gw\]) have suffices $\{jm\sigma, j^\prime m^\prime \sigma^\prime\}$, where $j$, $m$ and $\sigma$ stand for atomic sites in a unit cell, orbitals and spins, respectively. 

We assume that the Green’s function in the effective medium $G^0({i{\omega}_n})$ may be expressed as the self consistent equation of DMFT: $$\begin{aligned}
G^0({i{\omega}_n})^{-1}=G({i{\omega}_n})^{-1}+{\Sigma}({i{\omega}_n})+\tilde{\mu}-\mu,~
\label{eqn:Gw-local}\end{aligned}$$ where $\tilde{\mu}$ is the chemical potential in the effective medium. 

Iterative Perturbation Theory {#sec:IPT}
-----------------------------

In this paper, IPT method is adopted as a solver for the mapped single impurity problem, generalized by Fujiwara [*et al.*]{} for multi-orbital bands on arbitrary electron occupation with different Coulomb $U$. [@re:Fuji; @re:miura] In the IPT scheme, the self-energy is determined with the interpolation scheme between the high frequency limit and the strong interaction limit (the atomic limit) by using the second order self-energy.

The second order self-energy is calculated as $$\begin{aligned}
\!\!\!\!\!\!\!\!\Sigma^{(2)}_{m{\sigma}}(\tau) \!\!\!\! &=& \!\!\!\! - \!\!\!\!\!\!\!\!\!\!\sum_{m_1m_2m_3{\sigma}'}\!\!\!\!\!\!\!(U_{m_1m_2:mm_3}-U_{mm_2:m_3m_1}\delta_{{\sigma}{\sigma}'})\nonumber\\
&\times& \!\!\!\! U_{m_1m_3:mm_2}G^0_{m_1{\sigma}}(\tau)G^0_{m_3{\sigma}'}(\tau)G^0_{m_2{\sigma}'}(-\tau)
\label{eqn:2nd-pertu-Sig}\\
\!\!\!\!\!\!\!\!\Sigma^{(2)}(i\omega_{n}) \!\!\!\! &=& \!\!\!\!\! \int_{0}^{\beta} \!\!\!\! d\tau 
e^{i\omega_{n}\tau}\Sigma^{(2)}(\tau) \ . \label{eqn:2nd-pertu-Sig-fourier}  \end{aligned}$$ By using the second order self-energy $\Sigma^{(2)}({i{\omega}_n})$,  the self-energy is then expressed in a matrix form as $$\begin{aligned}
\! \! \Sigma(i\omega_n) \! = \! \Sigma^{{\mbox{\scriptsize{HF}}}} \! \! \! 
+\! A\Sigma^{(2)}(i\omega_n)[1\! -\! B(i\omega_n)\Sigma^{(2)}(i\omega_n)]^{-1},
\label{eqn:sig-IPT}   \end{aligned}$$ where $\Sigma^{{\mbox{\scriptsize{HF}}}}$ is the Hartree-Fock energy. The matrices $A$ and $B(i\omega_n)$ are determined by requiring the self-energy ${\Sigma}({i{\omega}_n})$ to be exact in the high-frequency limit $({i{\omega}_n}\to \infty)$,  and in the atomic limit $(U\to \infty)$ [@re:Fuji] as $$\begin{aligned}
\!\!\!\!\!\! A \!\! &=& \!\! \lim_{i{\omega}_n\to\infty}\{\Sigma(i\omega_n)-\Sigma^{HF}\}\Sigma^{(2)}(i\omega_n)^{-1},
\label{eqn:Sig-IPT-A}\\ 
\!\!\!\!\!\! B(i{\omega}_n) \!\! &=& \!\! \Sigma^{at(2)}(i\omega_n)^{-1}-\{\Sigma^{at}(i\omega_n)-\Sigma^{HF}\}^{-1}A ,
\label{eqn:Sig-IPT-B} \end{aligned}$$ where $\Sigma^{at}(i\omega_n)$ and $\Sigma^{at(2)}(i\omega_n)$ in Eq. (\[eqn:Sig-IPT-B\]) are the self-energy and the second order self-energy in the atomic limit.  These values are calculated exactly within the framework of the ligand field theory (LFT), which is discussed in Section \[sec:LF\]. 

The self-energy in Eq. (\[eqn:sig-IPT\]) also includes spin-flip and pair-hopping terms of the exchange interactions since the second order self-energy in Eq. (\[eqn:2nd-pertu-Sig\]) includes those terms. In addition, the IPT can be achieved within efficient CPU time since the self-energy in Eq. (\[eqn:sig-IPT\]) is obtained by only matrix calculation and hence the IPT is applicable to realistic materials with larger size by using full-LDA Hamiltonian. Form the discussion in this subsection, present IPT is much appropriate impurity solver and present IPT appropriately treats on-site dynamical correlation effects including both hole-hole and electron-electron scattering as well as electron-hole scattering.

In the following Sections (Secs. \[sec:result-Fe-Ni\] and \[sec:result-NiO\]), we assume the cubic symmetry around transition metal ions and, therefore, the self-energy can be assumed as $$\begin{aligned}
   \Sigma_{jm\sigma j'm'\sigma '}(i\omega_n)=
   \Sigma_{jm\sigma}(i\omega_n)\delta_{mm'}\delta_{jj'}\delta_{\sigma\sigma '}.
\label{eqn:sig-scalar}   \end{aligned}$$ For a cubic symmetry, the local Green’s function becomes diagonal after the ${\bm k}$-integration.

Calculation in the atomic limit based on the ligand field theory {#sec:LF}
----------------------------------------------------------------

The Hamiltonian mapped onto an isolated atom is defined as $$\begin{aligned}
 H_{atom}      &=& H^{dc}_{atom}+H_{int} \ ,
\label{eqn:HAM-atom}\\
 H^{dc}_{atom} &=& \sum_{(m{\sigma})}{\epsilon}_{mm'{\sigma}}c_{m{\sigma}}^{\dag}c_{m'{\sigma}} \ ,
\label{eqn:HAM-atom-dc}\\
 {\epsilon}_{mm'{\sigma}}&=&h^{({\mbox{\scriptsize{LDA}}})}_{mm'{\sigma}}-\Delta h^{{\sigma}}_{dc}\delta_{mm'} \ .
\label{eqn:ep-atom}\end{aligned}$$ The first term $H^{dc}_{atom}$ in Eq. (\[eqn:HAM-atom\]) is an unperturbed one-electron part. The Coulomb interaction $H_{int}$ is the same as Eq. (\[eqn:int\]). $h^{({\mbox{\scriptsize{LDA}}})}_{mm'{\sigma}}$ is the on-site element of LDA Hamiltonian and therefore $h^{({\mbox{\scriptsize{LDA}}})}_{mm'{\sigma}}$ includes the information of bulk matrix with electron transfer, crystal field and orbital hybridization. $\Delta h^{{\sigma}}_{dc}$ is the double counting term same as in Eq. (\[eqn:HLDAdc\]). Then one electron energy ${\epsilon}_{mm'{\sigma}}$ in Eq. (\[eqn:ep-atom\]) includes the crystal field splitting. Once we obtain ${\epsilon}_{mm'{\sigma}}$ and the Slater integrals, many-electron eigenstates of Eq. (\[eqn:HAM-atom\]) should be obtained within the framework of the ligand field theory (LFT). [@re:LFT] Full configuration interaction (CI) calculation is carried out by using the basis of all the multiplet of d electrons, the number of which is $2^{10}=1024$. Using Full-CI calculation yields an exact solution for local self-energy of Eq. (\[eqn:sig-IPT\]) in the atomic limit. The Green’s function, the self-energy and the second-order self-energy in the atomic limit are then defined as $$\begin{aligned}
G^{at}_{mm'{\sigma}}(i\omega_{n})
&=& {\frac{1}{Z}}\sum_{\eta \nu } 
{\frac{\langle \eta|c_{m{\sigma}}|\nu\rangle \langle \nu|c^{\dag}_{m'{\sigma}}|\eta\rangle}{i{\omega}_n +(E_{\eta}-\mu N_{\eta})-(E_{\nu}-\mu N_{\nu})}}
\nonumber\\
&\:\:& \:\: \:\: \:\: \:\:   \times \{e^{-\beta( E_{\eta}-\mu N_{\eta})}+e^{-\beta(E_{\nu}-\mu N_{\nu})}\}, 
\label{eqn:G-atom}\\
{\Sigma}^{at}_{mm'{\sigma}}(i{\omega})
&=& i{\omega}-({\epsilon}_{mm'{\sigma}}-\mu)-[G^{at}(i{\omega})]^{-1}_{mm'{\sigma}}, 
\label{eqn:Sig-atom}\\
\Sigma^{at(2)}_{m{\sigma}}(i\omega) 
&=& \sum_{m_1m_2m_3{\sigma}'}(U_{mm_2:m_1m_3}-U_{mm_2:m_3m_1}\delta_{{\sigma}{\sigma}'})\nonumber\\
\times && \! \! \! \! \! \! U_{m_1m_3:mm_2}\nonumber\\
\times && \! \! \! \! \! \! {\frac{ f_{m_2{\sigma}'}(1-f_{m_1{\sigma}}-f_{m_3{\sigma}'})+f_{m_1{\sigma}}f_{m_3{\sigma}'}}{i{\omega}- {\epsilon}_{m_1{\sigma}} - {\epsilon}_{m_3{\sigma}'} + {\epsilon}_{m_2{\sigma}'} +{\tilde{\mu}}}} \ ,
\label{eqn:sig2-atom}\end{aligned}$$ where $|\eta\rangle$, $E_{\eta}$, $N_{\eta}$, $f_{m{\sigma}}$ are the many-electron eigenstates in the Hilbert space of $2^{10}$ multiplets, the energy eigenvalue, the total electron number of the eigenstate $|\eta\rangle$ and the Fermi distribution function obtained from full-CI calculation within the framework of LFT.

One can assign, with a help of the atomic spectrum in LFT, the origin of peaks in LDA+DMFT spectra, e.g. the initial and final multiplets of the corresponding transitions. Generally, peak structures in LFT spectra would suggest the existence of strong scattering channels in LDA+DMFT spectra at around these particular energies.

From discussion in Sections \[sec:IPT\] and \[sec:LF\], the IPT is a appropriate approximation method and is worth being extended to the application of realistic materials in following four reasons: (1) Efficient computational cost (2) Applicability of realistic materials with larger size and more complicated hybridization by using full-LDA Hamiltonian \[This is followed by (1)\] (3) Care of spin-flip and pair-hopping terms of the exchange interactions in multi-orbital systems (4) Direct comparison of LFT spectra within the IPT method with the LDA+DMFT spectra and assignment of the origin of peaks in LDA+DMFT spectra

Energy spectrum and ${\bm k}$-resolved spectrum {#sec:Eq-Gw-Gkw}
-----------------------------------------------

The energy spectrum and the ${\bm k}$-resolved spectrum may be more physically understandable if we present the Green’s function in the orthogonal base, subject to the sum rule of the energy integration of the local Green’s function.  The lattice Green’s function and the local Green’s function in the orthogonal base are defined as $$\begin{aligned}
\!\!\!\!\!\!\!\!\![{\mathscr{G}}({\bm k}, {\omega})]^{-1}\!\!\!\! &=& \!\!\! O^{-\frac{1}{2}}({\bm k})G^{-1}({\bm k}, {\omega})O^{-\frac{1}{2}}({\bm k})
\label{eqn:Gkw-ortho-1}\\
&=&\!\!\! [({\omega}+\mu)
\nonumber\\
&-&\!\!\! O^{-\frac{1}{2}}({\bm k})
\{H^{dc}_{{\mbox{\scriptsize{LDA}}}}({\bm k})+{\Sigma}({\omega})\}O^{-\frac{1}{2}}({\bm k})],
\label{eqn:Gkw-ortho-2}\\
{\mathscr{G}}({\omega})&=& {\frac{1}{V}}\int d{\bm k}~{\mathscr{G}}({\bm k},{\omega}).~
\label{eqn:Gw-ortho-1}\end{aligned}$$ The energy spectrum will be shown with the imaginary part of ${\mathscr{G}}({\omega})$. ${\bm k}$-resolved spectrum will be also presented with an imaginary part of ${\mathscr{G}}({\bm k},{\omega})$, which may be compared with the angle-resolved photoemission spectra (ARPES).

Computational Details {#sec:Comp-Detals}
---------------------

The $8 \times 8 \times 8$ ${\bm k}$-points are used in the ${\bm k}$-integration within the whole Brillouin zone. The ${\bm k}$-integration is carried out by using a generalized tetrahedron method. [@re:Fuji]  The Padé approximation is adopted for analytic continuation of the Green’s function from the Matsubara frequencies ${i{\omega}_n}$ to the real ${\omega}$-axis.  We adopt $2^{11}=2048$ Matsubara frequencies.

Here, all the calculation were carried out at $T=1000~{\mbox{K}}$ to reduce the number of adopted Matsubara frequency. To consider temperature dependence of magnetic order transition, one should include entropy calculation as well as total energy calculation. [@re:Biermann-Ce] Here, we do not include entropy calculation within LDA+DMFT scheme. Thus temperature $T=1000{\mbox{K}}\sim 0.1{\mbox{eV}}$ causes the spectrum broadening by at most $0.1~{\mbox{eV}}$ and that do not change the magnetic order.

We first carry out LDA calculation to obtain input LDA Hamiltonian and then carry out DMFT calculation. We adopt the converged DMFT results as output. We do not feed back the density by DMFT to the LDA calculation.

In NiO case, on-site self-energy of nickel 3d bands on each sublattice is only included and inter-atom self-energy between nickel atoms on different sublattices are neglected. Thus, our calculation is “single-site” DMFT, but not cluster-DMFT. [@re:rev-cluster-DMFT]

Ferromagnetic bcc-Fe and fcc-Ni {#sec:result-Fe-Ni}
===============================

LSDA calculation has achieved a great success on overall understanding of the ground state electronic structures of crystalline transition metals, including bcc-Fe and fcc-Ni. However, LDA overestimates the width of occupied 3d bands for both ferromagnetic bcc-Fe and fcc-Ni and could not produce the observed satellite of spectrum at $6~{\mbox{eV}}$ below the Fermi energy in Ni (“$6~{\mbox{eV}}$ satellite”). On the other hand, LDA+DMFT calculation [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT] gives the width of occupied 3d bands for ferromagnetic bcc-Fe and fcc-Ni and the Ni $6~{\mbox{eV}}$ satellite, in good agreement with XPS experimental results. [@re:Fe-expt-XPS1; @re:Ni-expt-XPS1]

The aim of the calculation for ferromagnetic bcc-Fe and fcc-Ni in this section is to see whether LDA+DMFT-IPT(1) scheme can reproduce a reasonable width of occupied 3d bands and Ni $6~{\mbox{eV}}$ satellite. In fact, we will show that LDA+DMFT-IPT(1) presents the width of occupied 3d bands for bcc-Fe and fcc-Ni and the $6~{\mbox{eV}}$ satellite in fcc-Ni to be in good agreement with both experimental ARPES [@re:Fe-expt-XPS1; @re:Ni-expt-XPS1] and also with previous LDA+DMFT results. [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT]

Hamiltonian and $U$, $J$ values
-------------------------------

Lattice constants, atomic sphere radii in the TB-LMTO method and the averaged values of Coulomb and exchange interactions, $U$ and $J$, for ferromagnetic bcc-Fe and fcc-Ni to construct the Hamiltonian are summarized in Table \[tab:parameters\_and\_results\]. The values of $U$ and $J$ are obtained with the constraint-LDA calculation. [@re:Fe-Ni-UJ-paper] The basis set of the Hamiltonian consists of full valence bands of 4s, 4p and 3d orbitals, totally nine orbitals.

            lattice constant $a$ (Å)        atomic sphere radius $s_0$ (Å)             $U$ (eV)                     $J$ (eV)
  -------- ---------------------------- -------------------------------------- ------------------------- ------------------------------
   fcc-Ni             3.5233                            1.3768                  3.0[@re:Fe-Ni-UJ-paper]     0.9[@re:Fe-Ni-UJ-paper]
   bcc-Fe             2.8708                            1.4128                  2.0[@re:Fe-Ni-UJ-paper]     0.9[@re:Fe-Ni-UJ-paper]
    NiO               4.1948             1.2318 (Ni), 1.0699 (O), 0.8882 (ES)   7.0[@re:NiO-fujimori2]    0.9[@re:Anisimov-NiO-LSDA+U]

                                    LSDA    GW[@re:yamasaki-GW-TM]   LDA+DMFT       LDA+DMFT-IPT(1)        Expt.  
  ------------------------------- -------- ------------------------ ---------- -------------------------- ------- ------------------------------
      $\mu_{\rm spin}(\mu_B)$      bcc-Fe            2.17              2.31     2.28[@re:Fe-Ni-KKR+DMFT]   2.16    2.13[@re:Fe-Ni-expt-mag-mom]
                                   fcc-Ni            0.47              0.55     0.57[@re:Fe-Ni-KKR+DMFT]   0.47    0.57[@re:Fe-Ni-expt-mag-mom]
   $W_{\rm d, occ}({\mbox{eV}})$   bcc-Fe            3.7               3.4                                  3.4       3.3[@re:Fe-expt-XPS1]
                                   fcc-Ni            4.5               3.3                                  3.5       3.2[@re:Ni-expt-XPS1]

                                 LSDA   LSDA+U[@re:Anisimov-NiO-LSDA+U]   GW[@re:NiO-GW-Ferdi; @re:NiO-GW-Nohara]   QPscGW[@re:NiO-GW-Kotani]   LDA+DMFT[@re:Kunes-NiO-DMFT]   LDA+DMFT-IPT(1)             Expt.
  ----------------------------- ------ --------------------------------- ----------------------------------------- --------------------------- ------------------------------ ----------------- ----------------------------
   $E_{\rm gap}({\mbox{eV}}) $   0.2                  3.7                                   0.2                                4.8                          4.3                      4.3             4.3[@re:NiO-expt]
     $\mu_{\rm spin}(\mu_B)$     1.00                1.59                                  1.00                               1.72                          1.85                    1.00         1.64[@re:NiO-expt-mag-mom]

Energy Spectrum
---------------

Figure \[fig:dos-Fe-Ni\] shows the energy spectrum $-\mbox{Im} {\mathscr{G}}({\omega})$ obtained by LDA and LDA+DMFT-IPT(1) with the atomic spectrum by LFT. The occupied 3d-band becomes narrower in LDA+DMFT-IPT(1) than in LDA, narrower by $0.3~{\mbox{eV}}$ and $1.0~{\mbox{eV}}$ for bcc-Fe and fcc-Ni, respectively. Particularly, the spectra near the Fermi energy, in the region $-1.5~{\mbox{eV}}<{\omega}<0.0~{\mbox{eV}}$ for both bcc-Fe and fcc-Ni in LDA+DMFT-IPT(1), becomes narrower than that in the low energy region. These results are attributed to strong renormalization of the quasiparticle caused by the on-site dynamical electron correlation within the framework of DMFT. The precise width of occupied 3d bands is discussed in Sec. \[sec:Gkw-Fe-Ni-result\].

In addition, the satellite appears in LDA+DMFT-IPT(1) spectrum of fcc-Ni at $6~{\mbox{eV}}$ below the Fermi energy, in good agreement with experimental XPS results. [@re:Ni-expt-XPS1; @re:Ni-expt-spinXPS] Moreover, this satellite structure is strongly spin-dependent, much enhanced in the spectrum of the majority spin. These effects mainly come from the multiplet scattering of ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$, which is assigned by spin-dependent peak of LFT spectrum of ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$ at around $-5~{\mbox{eV}}$ in Fig. \[fig:dos-Fe-Ni\]-(2c). The satellite is due to hole-hole scattering process, [@re:Ni-satellite] and can appear only when one treats whole on-site dynamical correlation effects. It should be noted that $6~{\mbox{eV}}$ satellite is not observed in GW approximation (GWA) [@re:yamasaki-GW-TM] since GWA includes the dynamical correlation only within RPA (electron-hole excitations) and thus not the electron-electron and hole-hole scattering processes.

Let us focus on the atomic spectrum obtained by LFT in Fig. \[fig:dos-Fe-Ni\]. The atomic spectra are separately shown in Figs. \[fig:dos-Fe-Ni\]-(1c) (2c) for ionization spectra and Figs. \[fig:dos-Fe-Ni\]-(1d) (2d) for affinity spectra. The initial state is a mixture of ${\mbox{d}}^6$, ${\mbox{d}}^7$ in Fe and ${\mbox{d}}^8$, ${\mbox{d}}^9$ in Ni, since the total occupation numbers are non-integers, $n^{{\scriptsize{\mbox{Fe}}}}_{\rm tot}=6.829$ and $n^{{\scriptsize{\mbox{Ni}}}}_{\rm tot}=8.799$ from the results of bulk systems. The excited state is then a mixture of ${\mbox{d}}^5$, ${\mbox{d}}^6$, ${\mbox{d}}^7$, ${\mbox{d}}^8$ in Fe and ${\mbox{d}}^7$, ${\mbox{d}}^8$, ${\mbox{d}}^9$, ${\mbox{d}}^{10}$ in Ni. Small components appear in ionization spectra above $E_F$ and in affinity spectra below $E_F$. Moreover, characteristic atomic spectra at around $-0.5~{\mbox{eV}}$ in Fe and Ni originate from both ionization and affinity process; ${\mbox{d}}^6 {\rightarrow}{\mbox{d}}^7$ and ${\mbox{d}}^7 {\rightarrow}{\mbox{d}}^6$ in Fe in Figs. \[fig:dos-Fe-Ni\]-(1c)(1d) and ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^9$ and ${\mbox{d}}^9 {\rightarrow}{\mbox{d}}^8$ in Ni in Figs. \[fig:dos-Fe-Ni\]-(2c)(2d). These effects comes from the fact that higher energy occupied states are not fully occupied.

Characteristic structure in both atomic ionization and affinity spectra originate from various multiplet scattering; ${\mbox{d}}^7 {\rightarrow}{\mbox{d}}^6$ and ${\mbox{d}}^6 {\rightarrow}{\mbox{d}}^5$ at around $-4.0~{\mbox{eV}}$ and $-6.0~{\mbox{eV}}$ and ${\mbox{d}}^7 {\rightarrow}{\mbox{d}}^8$ and ${\mbox{d}}^6 {\rightarrow}{\mbox{d}}^7$ at around $1.5~{\mbox{eV}}$ in Fe of Figs. \[fig:dos-Fe-Ni\]-(1c)(1d) and ${\mbox{d}}^9 {\rightarrow}{\mbox{d}}^8$ and ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$ at around $-5.5~{\mbox{eV}}$ in Ni of Figs. \[fig:dos-Fe-Ni\]-(2c)(2d). This is due to the multiplet scattering caused by the exchange interaction $J$. The primary splitting of multiplets comes from the Coulomb interaction $U$, for example the transition spectrum of ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$ should locate in the lower energy region than that of ${\mbox{d}}^9 {\rightarrow}{\mbox{d}}^8$ of Ni spectra. Then additional multiplet splitting is caused by $J$. In LFT spectra in $-6.0~{\mbox{eV}}< {\omega}< 1.5~{\mbox{eV}}$, the energetic order of the spectral positions does not follow the above first simple rule and this fact implies that the scattering process by $J$ changes the multiplet spectra drastically. We conclude that the structures in LDA+DMFT-IPT(1) spectra is broadened and smoothed due to the multiplet scattering of ${\mbox{d}}^7 {\rightarrow}{\mbox{d}}^6$ and ${\mbox{d}}^6 {\rightarrow}{\mbox{d}}^5$ at $-6.0~{\mbox{eV}}< {\omega}< -2.0~{\mbox{eV}}$ and that of ${\mbox{d}}^7 {\rightarrow}{\mbox{d}}^8$ and ${\mbox{d}}^6 {\rightarrow}{\mbox{d}}^7$ at around $1.5~{\mbox{eV}}$ in bcc-Fe and that of ${\mbox{d}}^9 {\rightarrow}{\mbox{d}}^8$ and ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$ at around $-6.0~{\mbox{eV}}< {\omega}<-2.0~{\mbox{eV}}$ in fcc-Ni.

${\bm k}$-resolved spectrum and magnetic moment {#sec:Gkw-Fe-Ni-result}
-----------------------------------------------

Figure \[fig:Gkw-Fe-Ni\] shows the ${\bm k}$-resolved spectrum $- \mbox{Im}\,{\mathscr{G}}({\bm k},{\omega})$ by LDA+DMFT-IPT(1) with the LDA energy bands. LDA results overestimate the width of occupied 3d valence bands, which is defined to be the energy difference between the Fermi energy and the energy eigenvalue at P-point (Fe) or L-point (Ni), in comparison with experimental results. The width of occupied 3d valence bands and the magnetic moment in LDA+DMFT-IPT(1) are shown in Table \[tab:parameters\_and\_results\] in comparison with those by LDA and the experiments. [@re:Fe-Ni-expt-mag-mom; @re:Fe-expt-XPS1; @re:Ni-expt-XPS1] The valence band width of LDA+DMFT-IPT(1) is in reasonable agreement with experiments.

One can observe flat branches in bcc-Fe at around $-1.5~{\mbox{eV}}$ of majority spin and that at around $-2~{\mbox{eV}}$ of minority spin. These flat bands correspond to the local multiplet excitations of ${\mbox{d}}^7 {\rightarrow}{\mbox{d}}^6$. One can also observe flat branches in fcc-Ni at $-6~{\mbox{eV}}$ of both majority and minority spins. These flat bands at $-6~{\mbox{eV}}$ are due to the local multiplet excitations of ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$, where the intensity of the ${\bm k}$-resolved spectrum of majority spin is larger than that of minority spin. This intensity difference causes the strong spin dependence of $-6~{\mbox{eV}}$ satellite in the energy spectrum in Fig. \[fig:dos-Fe-Ni\]-(2b).

The ${\bm k}$-resolved spectrum in bcc-Fe is more diffusive than that in fcc-Ni at around $-3.0~{\mbox{eV}}\sim 0.0~{\mbox{eV}}$ in spite of smaller $U$ for bcc-Fe than fcc-Ni. This comes from that the scattering process by $J$ compared with $U$ is stronger in bcc-Fe than in fcc-Ni since $J/U$ in bcc-Fe is larger than that in fcc-Ni.

The magnetic moment $\mu_{\rm spin}$ of both bcc-Fe and fcc-Ni is almost the same as LDA result. $\mu_{\rm spin}$ of bcc-Fe in LDA+DMFT-IPT(1) result is also in good agreement with experiment, while that of fcc-Ni is slightly smaller than experimental result. This is an artifact due to insufficient number of ${\bm k}$-points to integrate the lattice Green’s function $G({\bm k}, {\omega})$ to obtain the local Green’s function $G({\omega})$ by using generalized tetrahedron method. The adopted value of the total number of the ${\bm k}$-points in LDA+DMFT-IPT(1) is much smaller than that in TB-LMTO method and the discrepancy of $\mu_{\rm spin}$ for fcc-Ni may be improved by increasing the total number of the ${\bm k}$-points in LDA+DMFT-IPT(1).

Comparison with previous LDA+DMFT results {#sec:Compare_LDA+DMFT_for_FeNi}
-----------------------------------------

Here, we compare the results for ferromagnetic bcc-Fe and fcc-Ni obtained by LDA+DMFT-IPT(1) with other previous LDA+DMFT [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT]. Lichtenstein [*et al.*]{} [@re:Fe-Ni-UJ-paper] has used the LDA+DMFT with QMC as an impurity solver. Minár [*et al.*]{} [@re:Fe-Ni-KKR+DMFT] has used the KKR+DMFT with perturbative SPTF (spin-polarized T-matrix+FLEX) as an impurity solver. The energy spectra in LDA+DMFT-IPT(1) show a good agreement with those LDA+DMFT calculations. The existence of Ni 6 eV satellite is also very similar to those LDA+DMFT calculations, presumably much better coincident with the position of the observed spectra. The magnetic moment of bcc-Fe is in good agreement with those LDA+DMFT results. Slightly smaller value of the magnetic moment of fcc-Ni in LDA+DMFT-IPT(1) than previous LDA+DMFT is not due to the use of perturbative IPT approach as an impurity solver but due to insufficient number of ${\bm k}$-points mentioned above. Thus, we can conclude that LDA+DMFT-IPT(1) reproduces reasonable results for ferromagnetic bcc-Fe and fcc-Ni and that LDA+DMFT-IPT(1) is applicable to realistic metallic materials in strongly correlated electron systems as well as other previous LDA+DMFT methods. [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT]

Antiferromagnetic NiO {#sec:result-NiO}
=====================

NiO is a type-II antiferromagnetic insulator with Néel temperature $523~{\mbox{K}}$. The experimentally observed band gap is $4.3~{\mbox{eV}}$.[@re:NiO-expt] Resonance photoemission experiments [@re:NiO-kyoumei] show that electronic structure of NiO should be of the charge-transfer type and the low energy satellite mainly consists of nickel 3d bands.

Various theoretical methods have been applied to NiO. The band gap $E_{\rm gap}$ and the magnetic moment $\mu_{\rm spin}$ are summarized in Table \[tab:NiO-past-result\]. LSDA calculation shows that NiO is Mott-Hubbard type insulator with a small band gap of $0.2~{\mbox{eV}}$ for antiferromagnetic phase, [@re:Terakura-NiO-LSDA] in which oxygen 2p bands are located at lower energy region than the occupied nickel 3d bands.

The LSDA+U method was applied and the resultant band gap is almost good agreement with experiment and a system becomes charge-transfer insulator. [@re:Anisimov-NiO-LSDA+U] However, bonding states of nickel ${\mbox{e}_{g}}$ bands are observed at around $8~{\mbox{eV}}$ below the Fermi energy and this causes less components of the occupied main peak and too much components of the occupied satellite peak, compared with experimental XPS spectrum. These problems mainly come from the static potential correction with orbital dependence in LSDA+U and this implies that one should include dynamical correlation effects.

GW approximation was applied and the resultant band gap still remains very small (0.2). [@re:NiO-GW-Ferdi; @re:NiO-GW-Nohara] Quasiparticle self-consistent GW (QPscGW) approximation [@re:NiO-GW-Kotani] was also applied and a system becomes charge-transfer insulator. However, the resultant band gap is overestimated (4.8) and the bonding states of nickel ${\mbox{e}_{g}}$ bands are observed at around $6~{\mbox{eV}}$ below the Fermi energy. These fact implies that the dynamical electron correlation should play a more crucial role in NiO than treated by RPA in GW approximation.

LDA+DMFT was applied to paramagnetic NiO [@re:Kunes-NiO-DMFT; @re:Ren-NiO-DMFT] and the resultant band gap is in good agreement with experimentally observed XPS result. [@re:NiO-expt] and the system becomes charge-transfer insulator. However, the top of valence bands mainly consists of nickel 3d bands, while the cluster-model CI calculation [@re:NiO-fujimori2; @re:NiO-vanElp] and oxygen x-ray absorption of Li$_x$Ni$_{1-x}$O [@re:NiO-Kuiper-expt] show that the top of valence bands is mainly based on oxygen 2p bands. This difference comes from that these LDA+DMFT calculations use projected effective Hamiltonian constructed by Wannier-like functions with fixing the hybridization mixing. Moreover, no calculation has been carried out for antiferromagnetic NiO by LDA+DMFT.

The target of LDA+DMFT calculation is to get electronic structures of antiferromagnetic NiO with (i) the accurate band gap, (ii) the correct description of the charge-transfer type insulator and (iii) the occupied main peak corresponding to oxygen 2p states and the occupied satellite peak corresponding to nickel 3d states.

Hamiltonian, the Coulomb and exchange interactions $U$ and $J$
--------------------------------------------------------------

The structure of NiO is not of the dense packing and we put empty atom spheres in vacant region of the lattice in LMTO formalism. The lattice constants, atomic sphere radii of each atom and the averaged values of Coulomb and exchange interactions, $U$ and $J$, are summarized in Table \[tab:parameters\_and\_results\]. The total number of atom spheres in an antiferromagnetic unit cell is eight; two Ni, two O and four empty atoms (ES). We use the full LDA Hamiltonian of the LMTO formalism and adopted muffin-tin orbitals are 4s, 4p and 3d in Ni and 2s, 2p in O and 1s, 1p in ES and, thus, the total number of basis is 42. $U$ refers to the experimental value of $U=7.0~{\mbox{eV}}$ [@re:NiO-fujimori2] and $J$ to the constraint-LDA calculation of $J=0.9.~{\mbox{eV}}$ [@re:Anisimov-NiO-LSDA+U] The difference of $U$-values in fcc-Ni and NiO is due to the difference of screening mechanism in metals and insulators.

The value of $U$ obtained by constraint LDA is $U=8.0~\mbox{eV}$, [@re:Anisimov-NiO-LSDA+U] which is slightly larger than experimental value. Within the framework of constraint LDA, all the screening channels are switch off. When we evaluate suitable value of $U$ for 3d-orbitals required in the framework of LDA+DMFT, only the screening channels associated with 3d electrons should be switch off and those associated with 4s and 4p electrons should remain. In this sense, a suitable value of $U$ is slightly smaller than that for constraint LDA. Thus, the adopted value of $U=7.0~\mbox{eV}$ in the present paper is reasonable within the framework of LDA+DMFT.

Energy Spectrum
---------------

Figure \[fig:dos-NiO\] shows the energy spectrum $-\mbox{Im} {\mathscr{G}}({\omega})$ of antiferromagnetic NiO by LDA+DMFT-IPT(1) with experimental XPS spectrum. [@re:NiO-expt] Experimental XPS spectrum mainly consists of three parts: a main peak at $4~{\mbox{eV}}$, a main peak at $-3.5~{\mbox{eV}}< {\omega}< -0.5~{\mbox{eV}}$ and a satellite peak at $-10.0~{\mbox{eV}}< {\omega}< -6.0~{\mbox{eV}}$ as in Fig. \[fig:dos-NiO\]-(4). With the cluster model CI calculation, [@re:NiO-fujimori2] these three structures are assigned to ${\mbox{d}}^9$, ${\mbox{d}}^8{\mbox{\underline{L}}}$ and ${\mbox{d}}^7$ final state, respectively, where ${\mbox{\underline{L}}}$ is a ligand hole created in oxygen 2p orbitals.

The energy spectrum by LDA+DMFT-IPT(1) reproduces these three structures fairly well, the positions of two main peaks are in good agreement with experimentally observed XPS result, but the position of satellite peak at $-10.0~{\mbox{eV}}< {\omega}< -6.0~{\mbox{eV}}$ shifts slightly upward in comparison with experimental results. [@re:NiO-expt] Appreciable component of oxygen 2p bands appears just at the top of the valence bands. The main peak of the occupied states at $-1.5$ eV and the satellite at $-6.0~{\rm eV}\sim -7.5~{\rm eV}$ originate from t$_{2g}$ orbitals with minority spin of Ni and e$_g$+t$_{2g}$ orbitals with majority spin of Ni, respectively. The main peak of conduction bands comes from the Ni-e$_g$ orbitals with minority spin. The spectrum of LDA+DMFT-IPT(1) in Fig. \[fig:dos-NiO\] shows that the electronic structure is of the charge-transfer insulator type, while the electronic structure in LDA is of Mott-Hubbard type. In fact, the hybridization mixing between oxygen 2p bands and nickel 3d bands is much enhanced in LDA+DMFT-IPT(1) in comparison with that in LDA, though the Coulomb matrix elements between Ni-3d and O-2p bands and among O-2p bands are not included.

On-site Coulomb interaction between nickel ${\mbox{e}_{g}}$ bands makes the occupied nickel ${\mbox{e}_{g}}$ bands shift to the lower energy side and the unoccupied nickel ${\mbox{e}_{g}}$ bands to the higher energy side. Due to the shift of the occupied nickel ${\mbox{e}_{g}}$ bands to lower energy side, the hybridization between Ni-${\mbox{e}_{g}}$ and Ni-${\mbox{t}_{2g}}$ increases and the occupied Ni-${\mbox{t}_{2g}}$ bands shift to the lower energy side. Since the occupied Ni-${\mbox{e}_{g}}$ and ${\mbox{t}_{2g}}$ bands shift to the lower energy side, the hybridization between Ni-3d and O-2p bands increases and whole bands are broadened in the region $-8.5~{\mbox{eV}}<{\omega}<-0.5~{\mbox{eV}}$. These effects enhance the hybridization mixing and change the character of NiO to be of the charge-transfer insulator type. The energy spectrum and the electronic structure by LDA+DMFT-IPT(1) is in good agreement with XPS experiments.

Atomic spectrum for single nickel ion obtained by LFT is shown in Fig. \[fig:dos-NiO\]-(3). Initial state is ${^3\mbox{A}_{2g}}({\mbox{d}}^8)$ of the electron configuration ${{\mbox{t}^{3}_{2g}}}_{{\uparrow}}{{\mbox{t}^{3}_{2g}}}_{{\downarrow}}{{\mbox{e}^{2}_{g}}}_{{\uparrow}}$. The excited states of atomic ionization spectra a1 and a2 are ${^4\mbox{T}_{1g}}({\mbox{d}}^7)$, b1 and b2 are ${^2\mbox{E}_{g}}({\mbox{d}}^7)$ and c1 and c2 are ${^2\mbox{T}_{1g}}({\mbox{d}}^7)$, respectively. This assignment is perfectly consistent with previous LFT calculation. [@re:NiO-fujimori2]

In the atomic calculation of LFT, a single nickel ion is considered and inter-atomic electron transfer is not allowed. Therefore, present atomic calculation does not include explicitly the ${\mbox{d}}^8{\mbox{\underline{L}}}$ final state. The initial state in the LFT spectra is fixed in ${\mbox{d}}^8$ since the total electron number of Ni ion is eight in NiO. Therefore, the unoccupied and occupied spectra in Fig. \[fig:dos-NiO\]-(3) originate from the transition of ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^9$ and ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$, respectively. Though the ${\mbox{d}}^8{\mbox{\underline{L}}}$ final state is not included in LFT, of Fig. \[fig:dos-NiO\]-(3), the occupied main peak in the spectrum of LDA+DMFT-IPT(1) is in good agreement with experimental XPS result. In fact, the initial states in LDA+DMFT-IPT(1) is a mixture of ${\mbox{d}}^8$ and ${\mbox{d}}^9{\mbox{\underline{L}}}$, since the electron occupation numbers of nickel 3d bands and oxygen 2p bands in LDA+DMFT-IPT(1) are 8.17 and 5.14, respectively. Therefore, we can assign, in the spectrum of LDA+DMFT-IPT(1), the unoccupied main peak and occupied satellite peak to be the ${\mbox{d}}^9,~{\mbox{d}}^7$ final state, respectively. The occupied main peak is also assigned to be the ${\mbox{d}}^8{\mbox{\underline{L}}}$ final state.

The splitting between ${\mbox{d}}^8{\mbox{\underline{L}}}$ and ${\mbox{d}}^7$ final state configurations is small and, in Fig. \[fig:dos-NiO\]-(4), the occupied satellite peak in LDA+DMFT-IPT(1) appear at higher energy region ($ -6 \sim -7$ eV) than that in experimental XPS result ($-8 \sim -9$ eV). This may be due to the fact that we do not include the ${\mbox{d}}^8{\mbox{\underline{L}}}$ final state in LFT calculation. The cluster model CI-calculation instead of LFT calculation would give more precise position of occupied satellite peak. 

Unoccupied ${\mbox{e}_{g}}$ bands of nickel majority spin locates at slightly higher energy region than that of minority spin and oxygen 2p bands in Fig. \[fig:dos-NiO\]-(1). It would not be the case, if we include the inter-atomic component of the self-energy or cluster CI-calculation.

The band gap $E_{\rm gap}$ and the magnetic moment $\mu_{\rm spin}$ are summarized in Table \[tab:NiO-past-result\], The band gap in LDA+DMFT-IPT(1) is in fairly good agreement with experimental result. On the other hand, the calculated magnetic moment is almost the same as the calculated result by LDA and much smaller than the observed one. The magnetic moment of antiferromagnetic NiO comes from the electron occupation of nickel ${\mbox{e}_{g}}$ bands since nickel ${\mbox{t}_{2g}}$ bands are almost fully occupied. The present result of smaller value of the magnetic moment in LDA+DMFT-IPT(1) than the experimental result comes from that unoccupied ${\mbox{e}_{g}}$ bands of nickel majority spin has some intensity. In the present calculation, we neglect the inter-atomic and inter-spin components of the self-energy and the discrepancy may be attributed to this approximation, which should be left as a future study.

${\bm k}$-resolved Spectrum
---------------------------

The ${\bm k}$-resolved spectrum $- \mbox{Im}\,{\mathscr{G}}({\bm k},{\omega})$ by LDA+DMFT-IPT(1) with both the energy bands by LDA (dashed lines) and the angle-resolved photoemission spectrum (dots) is shown in Fig. \[fig:Gkw-NiO\]. LDA band at $0.0~{\mbox{eV}}< {\omega}< 1.5~{\mbox{eV}}$ is shifted to the LDA+DMFT-IPT(1) band at $3.0~{\mbox{eV}}< {\omega}< 4.0~{\mbox{eV}}$. These spectra correspond to unoccupied nickel ${\mbox{e}_{g}}$ bands in Fig. \[fig:dos-NiO\]-(1). These unoccupied nickel 3d bands in LDA+DMFT-IPT(1) may be more localized than that in LDA.

In LDA band structure, the occupied bands mainly consists of two parts: (i) $-2.5~{\mbox{eV}}< {\omega}< 0.0~{\mbox{eV}}$ bands structure mainly based on Ni-3d bands and (ii) $-8.5~{\mbox{eV}}< {\omega}< -3.5~{\mbox{eV}}$ bands structure mainly based on O-2p bands.  In LDA+DMFT-IPT(1) results, (i) and (ii) have been shifted to lower and higher energy side, respectively. This shift of these results causes the increase of hybridization between Ni-3d and O-2p bands. Especially, strong hybridization occurs at L and $\Gamma$ points. 

We observe a broad and flat diffusive structure at $-5~{\mbox{eV}}$ in LDA+DMFT-IPT(1) spectrum, which corresponds to occupied satellite peak of Ni-${\mbox{t}_{2g}}$ bands and also Ni-${\mbox{t}_{2g}}$ state in LFT in Figs. \[fig:dos-NiO\]-(2)(3). Therefore, we conclude that Ni-${\mbox{t}_{2g}}$ bands may be more localized due to the strong scattering by $U$.

One can see flat branches at $-5.5eV< \omega <-5.0~{\mbox{eV}}$and at $ -6.5~{\mbox{eV}}< \omega < -6.0~{\mbox{eV}}$ in ${\bm k}$-resolved spectrum of LDA+DMFT-IPT(1). This corresponds to occupied satellite peak of Ni-${\mbox{e}_{g}}$ bands in Fig. \[fig:dos-NiO\]-(1).  This flat band implies that the occupied satellite peak of Ni-${\mbox{e}_{g}}$ bands is not due to the band structure but to the multiplet scattering, particularly ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$ assigned by the atomic calculation of LFT.

The ${\bm k}$-resolved spectrum at $-5.5~{\mbox{eV}}< {\omega}< -2.5~{\mbox{eV}}$ of LDA+DMFT-IPT(1) is in good agreement with experimental ARPES result [@re:NiO-ARPES] along $\Gamma$ to X point, which is due to the inclusion of dynamical correlation within DMFT, which causes strong hybridization between Ni-3d and O-2p bands.

Finally, we should mention that much enhanced intensity profile can be seen in $- \mbox{Im}\,G({\bm k},{\omega})$ but $- \mbox{Im}\,{\mathscr{G}}({\bm k},{\omega})$ becomes more diffusive due to overlap integrals.

Comparison with other calculations and experiments {#sec:Compare_with_others_NiO}
--------------------------------------------------

Here, we compare the results for antiferromagnetic NiO obtained by LDA+DMFT-IPT(1) with other previous calculations. LDA+DMFT-IPT(1) reproduces $4.3~{\mbox{eV}}$ band gap and the characteristics of the charge-transfer type insulator. Moreover, the results obtained by LDA+DMFT-IPT(1) shows no Ni-${\mbox{e}_{g}}$ bonding states below the Fermi energy, which is observed in LSDA+U [@re:Anisimov-NiO-LSDA+U] and QPscGW results. [@re:NiO-GW-Kotani] LDA+DMFT-IPT(1) also reproduces the top of valence bands with mainly based on oxygen 2p bands, while application of the LDA+DMFT method using Hamiltonian constructed by Wannier-like function to paramagnetic NiO [@re:Kunes-NiO-DMFT; @re:Ren-NiO-DMFT] shows that the top of valence bands is mainly based on nickel ${\mbox{e}_{g}}$ bands. Those effects obtained by LDA+DMFT-IPT(1) are in good agreement with experimental XPS results. Those drastic change in LDA+DMFT-IPT(1) is due to enhancement of hybridization between Ni-3d and O-2p bands caused by on-site Coulomb interaction of nickel 3d bands within DMFT scheme.

The assignment of final states of the main and satellite occupied peaks and the unoccupied main obtained in LDA+DMFT spectrum is consistent with that in the cluster-model CI calculation. [@re:NiO-fujimori2] This comes from the use of LFT spectrum within the IPT method and hence the assignment of the origin of peaks in LDA+DMFT spectra by using LFT spectra is a great advantage of the IPT method, compared with other impurity solvers.

However we still have several problems, slightly higher energy positions of the occupied satellite peak and the unoccupied main peak and smaller value of the magnetic moment. From discussion in the previous two subsections, these are mainly due to neglect of the inter-atomic and inter-spin components of the self-energy within DMFT scheme. Particularly, dynamical correlation of nickel 3d bands between two different sublattices and that of nickel 3d and oxygen 2p bands are not included since DMFT neglects all the inter-atomic dynamical correlation.

To improve those problems in antiferromagnetic NiO, one may adopt cluster DMFT [@re:rev-cluster-DMFT], which is based on mapping of many electron systems in bulk onto single cluster impurity problem with including the intra-cluster dynamical Coulomb interaction and neglecting the inter-cluster Coulomb interaction. In LDA+DMFT-IPT(1), the extension of LDA+DMFT to LDA+cluster DMFT is consistent with the extension of LFT of single nickel isolated ion to the cluster model CI-calculation of NiO$_6$ cluster similar to Ref.  and . The use of the cluster model CI-calculation will reproduce better positioning of occupied satellite peak and direct assignment of occupied main peak with ${\mbox{d}}^8{\mbox{\underline{L}}}$ final state.

From above discussion, we can conclude that LDA+DMFT-IPT(1) reproduces reasonable results for antiferromagnetic NiO within the DMFT scheme. and that LDA+DMFT-IPT(1) is applicable to realistic compound cases and insulating cases as well as LDA+DMFT with nonperturbative impurity solver.

Conclusion {#sec:CONC}
==========

In summary, we have proposed LDA+DMFT-IPT(1) where we use the full LDA Hamiltonian and IPT as an impurity solver. We then applied LDA+DMFT-IPT(1) to ferromagnetic bcc-Fe, fcc-Ni and antiferromagnetic NiO.

For bcc-Fe and fcc-Ni case, on-site dynamical correlation effect causes the narrowing of the width of occupied 3d bands. The multiplet scattering effect of ${\mbox{d}}^8 {\rightarrow}{\mbox{d}}^7$ assigned by the LFT spectrum causes the spin dependent $6~{\mbox{eV}}$ satellite for fcc-Ni. The energy spectra and magnetic moment in LDA+DMFT-IPT(1) show a good agreement with previous LDA+DMFT calculations [@re:Fe-Ni-UJ-paper; @re:Fe-Ni-KKR+DMFT].

For antiferromagnetic NiO case, on-site Coulomb interaction of nickel 3d bands enhances the hybridization between Ni-3d and O-2p bands and this causes the band gap of $4.3~{\mbox{eV}}$, the charge-transfer type insulator, and the top of valence bands with mainly oxygen 2p bands. Those effects are all in good agreement with experiments. The drastic change of hybridization between nickel 3d and oxygen 2p bands is due to the use of full-LDA Hamiltonian and on-site dynamical electron correlation of nickel site within DMFT. The successful application of LDA+DMFT-IPT(1) to antiferromagnetic NiO implies that LDA+DMFT-IPT(1) is fairly applicable to larger compound cases with more complicated hybridization between more than two atoms.

In addition, the assignment of characteristic peaks of NiO in LDA+DMFT spectrum is consistent with that in the cluster-model CI calculation. [@re:NiO-fujimori2] Thus, the assignment of the origin of peaks by using LFT spectra is a great advantage of the IPT method in understanding the existence of strong scattering channels in LDA+DMFT spectra.

Several remained problems in antiferromagnetic NiO, positions of the occupied satellite and the unoccupied main peak, would be solved by including the inter-atomic dynamical correlation. Thus, we can conclude that LDA+DMFT-IPT(1) is applicable to various realistic materials, such as both metallic and insulating cases, multi-atom (compound) cases, spin-polarized cases and strongly hybridized cases between s, p and d-bands.

[99]{}

A.I. Lichtenstein, V.I. Anisimov and J. Zaanen, Phys. Rev. B [**52**]{} R5467 (1995).

A. Georges, G. Kotliar, W. Krauth and M.J. Rozenberg, Rev. Mod. Phys. [**68**]{} 13 (1996).

V.I. Anisimov, A.I. Poteryaev, M.A. Korotin, A.O. Anokhin and G. Kotliar, J. Phys: Condens. Matter [**9**]{} 7539 (1997).

A.I. Lichtenstein and M.I. Katsnelson, Phys. Rev. B [**57**]{} 6884 (1998).

L. Hedin, Phys. Rev. [**139**]{} A796 (1965).

S. Biermann, F. Aryasetiawan and A. Georges, Phys. Rev. Lett. [**90**]{} 086402 (2003).

M.J. Rozenberg, X. Y. Zhang and G. Kotliar, Phys. Rev. Lett. [**69**]{} 1236 (1992).

A. Georges and G. Krauth, [Phys. Rev. Lett.]{} [**69**]{} 1240 (1992).

K. Held and D. Vollhardt, [Eur. Phys. J. B]{} [**5**]{} 473 (1998).

M. Feldbacher, K. Held, and F. F. Assaad, [Phys. Rev. Lett.]{} [**93**]{} 136405 (2003).

A. N. Rubtsov, V. V. Savkin and A. I. Lichtenstein, [Phys. Rev. B]{} [**72**]{} 035122 (2005).\
P. Werner, A. Comanac, L. de’ Medici, M. Troyer and A. J. Millis, [Phys. Rev. Lett.]{} [**97**]{} 076405 (2006).\
S. Sakai, R. Arita, K. Held and H. Aoki, [Phys. Rev. B]{} [**74**]{}, 155102 (2006).

H. Kajueter and G. Kotliar, [Phys. Rev. Lett.]{} [**77**]{} 131 (1996).

T. Fujiwara, S. Yamamoto and Y. Ishii, [J. Phys. Soc. Jpn.]{} [**72**]{} 777 (2003).

O. Miura and T. Fujiwara, J. Phys. Soc. Jpn. [**75**]{} 014703 (2006).

M. S. Laad, L. Craco, and E. Müller-Hartmann, [Phys. Rev. Lett.]{} [**91**]{} 156402 (2003).\
M. S. Laad, L. Craco, and E. Müller-Hartmann, [Phys. Rev. B]{} [**73**]{}, 045109 (2006).

Th. Pruschke and N. Grewe, [Z. Phys. B]{} [**74**]{} 439 (1989).

A. Georges and G. Krauth, [Phys. Rev. B]{} [**48**]{} 7167 (1993).

 J.E. Hirsch, Phys. Rev. B [**31**]{}, 4403 (1985);\
J. E. Hirsch and R.M. Fye, Phys. Rev. Lett. [**56**]{}, 2521 (1986).

Y. Kuramoto and H. Kojima, [Z. Phys. B]{} [**57**]{} 95 (1984).

I.A. Nekrasov, K. Held, N. Blümer, A.I. Poteryaev, V.I. Anisimov, and D. Vollhardt, [Eur. Phys. J. B]{} [**18**]{}, 55 (2000).

A.I. Lichtenstein, M.I. Katsnelson and G. Kotliar, Phys. Rev. Lett. [**87**]{} 067205 (2001).

J. Minár, L. Chioncel, A. Perlov, H. Ebert, M.I. Katsnelson, and A.I. Lichtenstein, Phys. Rev. B [**72**]{}, 045125 (2005).

J. Kuneš, V.I. Anisimov, A.V. Lukoyanov, and D. Vollhardt, Phys. Rev. B [**75**]{}, 165115 (2007).

A. Liebsch, Phys. Rev. B [**23**]{} 5203 (1981);\
D.R. Penn, Phys. Rev. Lett. [**42**]{}, 921 (1979).

O.K. Andersen and O. Jepsen, Phys Rev. Lett. [**53**]{} 2571 (1984);\
O.K. Andersen, O. Jepsen, and D. Glötzel, in [*Canonical Description of the Band Structures of Metals*]{}, Proceedings of the International School of Physics, “Enrico Fermi”, Course LXXXIX, Varenna, 1985, edited by F. Bassani, F. Fumi, and M.P. Tosi (North-Holland, Amsterdam, 1985), p. 59.

F.M.F. de Groot, J.C. Fuggle, B.T. Thole, and G.A. Sawatzky, [Phys. Rev. B]{} [**42**]{} 5459 (1990).

T. Pruschke and R. Bulla, [Eur. Phys. J. B]{} [**44**]{}, 217 (2005).

E. Müller-Hartmann, [Z. Phys. B-Condens. Matter]{} [**76**]{} 211 (1989).

Y. Tanabe and S. Sugano, J. Phys. Soc. Jpn. [**9**]{} 753,766 (1954). B. Amadon, S. Biermann, A. Georges, and F. Aryasetiawan, Phys. Rev. Lett. [**96**]{}, 066402 (2006).

T. Maier, M. Jarrell, T. Pruschke and M. H. Hettler, Rev. Mod. Phys. [**77**]{}, 1027 (2005).

H. Höchst, A. Goldmann and S. Hüfner, Z. Phys. B 24, 245 (1976) ;\
R.E. Kirby, E. Kisker, F.K. King and E.L. Garwin, Solid State Commun. [**56**]{}, 425 (1985) ;\
H. Höchst, S. Hüfner and A. Goldmann , Phys. Lett. 57A, 265 (1976).

H. Höchst, S. Hüfner and A. Goldmann, Z. Phys. B 26, 133 (1977) ;\
M. Lädeniemi, E. Ojala and M. Okoochi, Phys. Status Solidi B [**108**]{}, K61 (1981).

A. Fujimori and F. Minami, Phys. Rev. B [**30**]{}, 957 (1984).

V.I. Anisimov, J. Zaanen and O.K. Andersen, Phys. Rev. B [**44**]{}, 943 (1991).

T. Kinoshita, T. Ikoma, A. Kakizaki, T. Ishii J. Fujii, H. Fukutani K. Shimada, A. Fujimori, T. Okane and S. Sato, Phys. Rev. B [**47**]{}, 47, 6787(1993).

A. Yamasaki and T. Fujiwara, J. Phys. Soc. Jpn. [**72**]{} 607 (2003).

M. B. Stearns, [*in Magnetic Properties in Metals*]{}, edited by H.P.J. Wijn, Landolt-Börnstein, New Series, Vol. III, Pt. 19a (Springer, Berlin, 1984) ;\
D. Bonnenberg, K. A. Hempel and H. P. J. Wijn, [*ibid*]{}.

G.A. Sawatzky and J.W. Allen, Phys. Rev. B [**53**]{} 2339 (1984).

M.R. Thuler, R.L. Benbow, and Z. Hurych, Phys. Rev. B [**27**]{}, 2082 (1983).

F. Aryasetiawan and O. Gunnarsson, Phys. Rev. Lett. [**74**]{}, 3221 (1995).

Y. Nohara and T. Fujiwara, private communication.

S. V. Faleev, M. van Schilfgaarde and T. Kotani, Phys. Rev. Lett. [**93**]{}, 126406 (2004).

A.K. Cheetham and D. A. O. Hope, Phys. Rev. B [**27**]{}, 6964 (1983).

K. Terakura, T. Oguchi, A.R. Williams and J. Kuebler, Phys. Rev. B [**30**]{}, 4734 (1984).

X. Ren, I. Leonov, G. Keller, M. Kollar, I. Nekrasov, and D. Vollhardt, Phys. Rev. B [**74**]{}, 195114 (2006).

J. van Elp, H. Eskes, P. Kuiper, and G. A. Sawatzky, Phys. Rev. B [**45**]{}, 1612 (1992).

P. Kuiper, G. Kruizinga, J. Ghijsen, G. A. Sawatzky, and H. Verweij, Phys. Rev. Lett. [**62**]{}, 221 (1989).

Z.-X. Shen, R.S. List, D.S. Dessau, B.O. Wells, O. Jepsen, A.J. Arko, R. Barttlet, C.K. Shih, F. Parmigiani, J.C. Huang and P.A.P. Lindberg, Phys. Rev. B [**44**]{}, 3604 (1991).
