Comparison Between the f-Electron and Conduction-Electron Density of States in the Falicov-Kimball Model at Low Temperature

TL;DR

This paper investigates whether the spectral gap at low temperature is the same for f-electrons and conduction electrons in the Falicov-Kimball model, finding strong evidence that they are indeed identical.

Contribution

The study provides new insights into the low-temperature spectral properties of the Falicov-Kimball model, specifically comparing the electron densities of states.

Findings

Spectral gaps for f-electrons and conduction electrons are the same as temperature approaches zero.

The f-electron density of states exhibits an orthogonality catastrophe singularity in the metallic phase.

The f-electron density of states is gapped in the insulating phase.

Abstract

The spinless Falicov-Kimball model is one of the simplest models of many-body physics. While the conduction-electron density of states is temperature independent in the normal state, the f-electron density of states is strongly temperature dependent---it has an orthogonality catastrophe singularity in the metallic phase and is gapped in the insulating phase. The question we address here is whether the spectral gap is the same for both electron species as T->0. We find strong evidence to indicate that the answer is affirmative.

Equations

$${\mathcal{H}}_{FK}=\sum _{ij}\left(-t_{ij}-\mu {\delta }_{ij}\right)c_i^{†}{c}_j^{}+E_f\sum _i{f}_i^{†}{f}_i^{}+U\sum _i{f}_i^{†}{f}_i^{}{c}_i^{†}{c}_i^{}$$

$${\rho }_{Z\to \infty }(ϵ)=\frac{\sqrt{4t^{\ast 2}-{ϵ}^2}}{2\pi t^{\ast 2}}$$

$$G_0(t)=-\frac{i}{\pi }{\int }_{-\infty }^{\infty }d\omega \text{Im}\left[{\mathcal{G}}_0(\omega )\right]e^{-i\omega t}[f(\omega )-\theta (t)]$$

$${\mathcal{G}}_0(\omega )=\frac{1}{\omega +i\delta +\mu -\lambda (\omega +i\delta )}$$

$$G_f^{>}(t,t^{\prime })$$

$$G_f^{<}(t,t^{\prime })$$

$$G_{f}^{c}(t_{c},t_{c}^{\prime})=-i\langle{\mathcal{T}}_{c}\big{(}f(t_{c})f^{\dagger}(t_{c}^{\prime})\big{)}\rangle$$

$$G_{f}^{>}(t)=-iw_{0}e^{-i(E_{f}-\mu)t}\text{Det}_{[0,t]}\big{|}\delta(t_{1}-t_{2})-UG_{0}(t_{1}-t_{2})\big{|}$$

$${\mathrm{Det}}_{[0,t]}$$

$$\bar{C}(\omega )=\frac{1}{\bar{P}(\omega )\bar{Q}(-\omega )}.$$

$$A_f(\omega )=-\frac{2}{\pi }{\int }_0^{\infty }dt\cos (\omega t)\text{Im }{G}^{>}(t).$$

Full text

1 Introduction

The Falicov-Kimball model originalFK is perhaps the simplest solid-state model for describing strongly correlated electron systems. The model possesses a Mott-insulator transition, a charge-density-wave ordered phase and can be solved exactly in the limit of infinite dimensions Brandt . Initially the model was employed to understand phase transitions in transition-metal oxides originalFK , and has since been used to study a number of strongly correlated systems. Its solution via dynamical mean-field theory has also been reviewed jimFK .

Recently, the Falicov-Kimball model has been used to investigate core-level X-ray photoemission spectroscopy jimXPS1 jimXPS2 . It has been known for some time that the localized $f$-electron Green’s function in the metallic phase of this model is related to the x-ray edge singularity problem XRAYEDGE1 ; XRAYEDGE2 (we only reference papers immediately relevant to this work here; a more complete history of X-ray edge and the Falicov-Kimball model appears elsewhere jimXPS2 ). In particular, at low temperature, the density of states in the metal displays a power-law divergence with an interaction-dependent power law.

The $f$-electrons, unlike the conduction electrons, possess a non-trivial and highly temperature dependent spectral function. This has been studied previously using the numerical renormalization group NRG illustrating the power law divergence of the $f$-electron spectrum at $\omega=0$ in the metallic phase. Here we seek to understand the gapped behavior of the $f$-electron spectral function in the Mott-insulating regime, where NRG approaches are known to have accuracy issues. Particularly, we ask if the $f$-electron gap approaches the (temperature-independent) conduction-electron gap in the limit $T\to 0$ in the insulating phase. This question may seem like it should have an obvious answer, because the $f$-electron dynamics are inherited through their interaction with the conduction electrons. But at nonzero temperature, we clearly see $f$-electron spectral weight within the gap, so the answer to the question is far from obvious.

2 Formalism

The spinless Falicov-Kimball model describes two electron species, conduction ($c$) and localized ($f$) electrons, which interact via a local Couloumb repulsion $U$ when occupying the same lattice site $i$. We represent the creation (destruction) of a conduction electron at the site $i$ by the second quantized operator $c^{\dagger}_{i}$ ($c_{i}^{\phantom{\dagger}}$) and of a localized electron by $f^{\dagger}_{i}$ ($f_{i}^{\phantom{\dagger}}$). Assuming a common chemical potential $\mu$ for the two species, the spinless Falicov-Kimball Hamiltonian is given by

[TABLE]

where $E_{f}$ is the $f$-electron site energy, and $t_{ij}$ is the hopping matrix, with $t_{ij}=t=t_{ji}$ for nearest-neighbor sites $i$ and $j$ with the hopping integral $t$ a real constant. We work at half filling ($\mu=U/2$ and $E_{f}=-U/2$) on the infinite-coordination Bethe lattice in equilibrium.

Scaling the hopping energy $t$ with coordination number $Z$ as $t=t^{*}/\sqrt{Z}$ then in the limit $Z\to\infty$ we obtain the noninteracting density of states on the Bethe lattice, given by

[TABLE]

with $t^{*}=1$ our energy unit and we restrict to $|\epsilon|\leq 2t^{*}$.

The dynamical mean-field theory (DMFT) approach solves the many-body problem by mapping the lattice problem onto an impurity problem because the self-energy has no momentum dependence. A self-consistent iterative approach is employed to determine the Green’s functions, with details given in jimFK . We do not need the full conduction-electron solution here, instead, we require only the bare time-ordered Green’s function for the conduction electrons $G_{0}(t)$ given by

[TABLE]

with $f(\omega)$ the Fermi-Dirac distribution,

[TABLE]

the effective medium, and the dynamical mean-field $\lambda(\omega)$ obtained from the DMFT algorithm jimFK . Here $\theta(t)$ is the Heaviside unit-step function.

The conduction electrons behave like noninteracting particles, in that they have a temperature-independent density of states, but also like interacting electrons, since they have a Mott metal-insulator transition at $U=2$.

We define the greater and lesser $f$-electron Green’s functions by

[TABLE]

where the angle brackets denote the thermal average, $\langle\ldots\rangle={\text{Tr}(\exp(-\beta\mathcal{H}_{FK})(...))}/{\text{Tr}\exp(-\beta{\mathcal{H}}_{FK})}$ The creation and annihilation operators are in the Heisenberg representation. These Green’s functions can be determined by selecting $t$, $t^{\prime}$ on certain branches of the Kadanoff-Baym-Keldysh contour (Fig. 1) and using the contour-ordered Green’s function,

[TABLE]

where $t_{c}$, $t_{c}^{\prime}$ are two times on the contour (Fig. 1). ${\mathcal{T}}_{c}$ is the contour time-ordering operator, arranging the operators such that $t_{c}^{\prime}$ lies before $t_{c}$ on the contour. For example, if we pick $t_{c}^{\prime}$ on the upper real-time branch of the contour and $t_{c}$ on the lower real-time branch, we recover $G_{f}^{>}(t,t^{\prime})$ from the contour-ordered Green’s function.

It has been shown in Ref. main1 , that in equilibrium the greater Green’s functions for the $f$-electrons take the form of a Toeplitz determinant of a continuous matrix operator over only the positive time branch of the contour

[TABLE]

where $w_{0}$ is the average density of sites without an $f$-electron ($w_{0}=\frac{1}{2}$ at half filling), and $G_{0}(t)$ is the bare time-ordered Green’s function determined from the dynamical mean-field $\lambda(\omega)$. The symbols $t_{1}$ and $t_{2}$ denote the matrix indices of the continuous matrix operator for which we evaluate the determinant; note that both times must fall within the interval $[0,t]$. (There is a similar expression for the lesser Green’s function.) To approximate this continuous matrix operator, we discretize it to a conventional matrix and calculate the determinant for three different discretization time steps $\Delta t$. We then perform a second-order Lagrange interpolation to extrapolate to the limit $\Delta t\to 0$. These numerical results are checked for accuracy against known spectral moments of the Green’s functions SpecMom .

In the limit of large times, an exact analytic formula for the $f$-electron Green’s function in equilibrium can be obtained using a factorization technique from complex analysis described by McCoy and Wu Mccoy and called the Wiener-Hopf sum approach. It relies on a result for infinite-sized determinants of Toeplitz matrices called Szego’s theorem. Further finite-time approximations can be made to improve the short-time agreement of this asymptotically exact result with the determinant calculation given by Eq. (8) main1 , with slightly different formulas for the case when the interaction energy $U$ lies above or below $U_{c}$, a critical interaction strength $U_{c}=\sqrt{2}$ on the Bethe lattice where the complex function $C(\omega)=1-UG_{0}(\omega)$ goes from no winding around the origin for $U<U_{c}$ to winding once in the clockwise direction (Ind$\,C$=-1) for $U\geq U_{c}$; note that this critical $U_{c}$ is smaller than the critical $U$ for the Mott transition. In the regime with $U>U_{c}$, the analytic result for the Toeplitz determinant in Eq. (8) is

[TABLE]

where $\bar{C}(\omega)=\exp[i\omega\Delta t][1-UG_{0}(\omega)]$, $\,\bar{g}(t)=\int_{-\infty}^{\infty}\,d\omega$ 

$\times\exp[-i\omega t]\ln(\bar{C}(\omega))/2\pi$, and $\bar{P}$ $(\bar{Q})$ are integrals over all positive (negative) time of $\bar{g}(t)$ and satisfy

[TABLE]

This expression represents a significant reduction in computational complexity, and allows us to probe a much wider parameter space when directly calculating the discretized determinant is not possible. Our approach will be to calculate the determinant directly for short times and use the analytic expression for long times, allowing us to obtain the spectral function of the $f$-electrons down to temperatures significantly lower than previous calculations. We need to patch the two solutions together smoothly, as described below. We examine the behavior of the $f$-electrons as they approach their $T=0$ limit in both the metallic phase near the Mott-insulator transition $U=2$ and in the Mott insulating regime.

Finally, we define the local density of states of the $f$-electrons, $A_{f}(\omega)$. At half-filling, there exists a particle-hole symmetry in our system [$A_{f}(\omega)=A_{f}(-\omega)$], and consequently the full $f$-electron density of states can be expressed as a Fourier transform of Im $G^{>}(t)$ alone Brandt1992 ,

[TABLE]

3 Numerics

We examine the temperature-dependent dynamics of the $f$-electron spectral function in the metallic regime near the Mott transition as well as in weakly and strongly correlated insulating phases above the Mott transition, which occurs at $U=2$. The spectral function is known to have a power-law like divergence in the metallic phase that disappears as we move into the Mott phase NRG . In the Mott insulator, the $f$-electron density of states develops a gap with decreasing temperature main2 . We explore the character of this gap to see if it is the same as in the temperature-independent conduction electron density of states.

To utilize the Weiner-Hopf technique in Eq. (9), we perform the discretized matrix calculation of Eq. (8) out to the longest time computationally feasible for three different time steps $\Delta t$ in the ratio $1:2:4$. Next we employ a second-order Lagrange extrapolation to take the limit $\Delta t\to 0$. We use Eq. (9) to calculate the determinant out to even longer times, and use a weighted blending of the two functions over a time range where the analytic approximation is roughly parallel to the discretized determinant results. This procedure works exceptionally well for high temperature, as shown in Fig. 4, but as we decrease the temperature, the determinant calculation requires significantly smaller time steps and takes longer to reach the ”long-time” regime where our analytic result holds (see Fig. 5). This effect limits our low-temperature calculations, and along with truncating the calculation at a finite time, leads to numerical artifacts near $\omega=0$ (see Figs. 7, 8, 9).

4 Results

In the metallic regime, for $U=1.5$ (Fig. 3) we see evidence of a power-law divergence as the temperature approaches zero in agreement with NRG . In the time domain, we observe a delayed decay towards zero with decreasing temperature. Note there is no zero crossing in Im $G^{>}(t)$ (Fig. 2), indicating we are still in the metallic phase.

For $U=1.9$, closer to the Mott transition, we still see the power-law like increase in the spectral function (Fig. 6), but with the development of a kink in the center of our density of states, representing a precursor to the insulating Mott phase.

At the Mott transition, $U=2$ (Fig. 7), the $f$-electrons partially fill the insulating gap at high temperature, with a kinked density of states that does not reach zero at the minimum temperature we were able to reach ($T=0.001$). We expect that the density of states should touch zero precisely at $T=0$. This is in contrast with the conduction electrons, whose density of states touches zero at $U=2$ for all temperatures.

For $U=2.5$ (Fig. 7), we see a similar transfer of weight from high to low frequencies as described in main2 . Here we see the spectral function becomes gapped around $T=0.1$, but continues to change shape down to $T=0.001$, below which the gap is frozen into place. Notice in the right inset of Fig. 8 the small region where the density of states becomes slightly negative. This is an artifact of not properly capturing the long-time behavior, likely due to the finite time truncation of some long period oscillations in $G^{>}(t)$ and possibly the blending of the numerical and asymptotic results.

Finally we examine the strongly correlated insulator with $U=3$ (Fig. 9). We see a similar temperature evolution to $U=2$ and $U=2.5$. In this case, the system becomes gapped around $T=0.25$, higher than $U=2.5$. Below this temperature $A_{f}(\omega)$ changes more slowly as the low-temperature behavior is frozen in. We observe a general trend of the gap appearing and freezing into place at a higher temperature for more strongly correlated materials. Again we see the irregular behavior near $\omega=0$ from the finite time truncation. More interesting, however, is the comparison with the conduction-electron density of states. We see from the left inset of Fig. 9 that at low temperature the width of the gap is nearly indistinguishable between the two different electrons—this strongly suggests that the $T=0$ gap of the $f$-electrons is the same as the conduction-electron gap.

5 Conclusion

We extensively studied the properties of the $f$-electron spectra of the Falicov-Kimball model in a number of interesting cases: just below the Mott transition, in the strongly correlated regime, and at temperatures approaching the $T=0$ limit. By using the Weiner-Hopf technique, we were able to examine these more complex situations by obtaining an analytic expression for the long-time behavior which is computationally much more efficient than the direct determinant calculation, because the matrices grow in size with increasing time. Due to the rich temperature-dependent dynamics of the $f$-electrons, we asked if the spectral gap in the $f$-spectrum approaches that of the (temperature-independent) conduction electrons at $T=0$. We found strong evidence that this is the case. This technique could be pushed further by making higher-order finite-time corrections to the asymptotic formulas given in main1 and by carrying the calculations out to longer times.

Acknowledgements.

We would like to thank A. M. Shvaika for valuable conversations and suggestions. This work was supported by the Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Contract No. DE-SC0019126. J. K. F. was also supported by the McDevitt bequest at Georgetown.

Conflict of interest

The authors declare that they have no conflict of interest.