Localized-itinerant dichotomy and unconventional magnetism in SrRu2O6.

Electron correlations tend to generate local magnetic moments that usually order if the lattices are not too frustrated. The hexagonal compound SrRu2O6 has a relatively high Neel temperature but small local moments, which seem to be at odds with the nominal valence of Ru5+ in the [Formula: see text] configuration. Here, we investigate the electronic property of SrRu2O6 using density functional theory (DFT) combined with dynamical-mean-field theory (DMFT). We find that the strong hybridization between Ru d and O p states results in a Ru valence that is closer to +4, leading to the small ordered moment ~1.2 μ B . While this is consistent with a DFT prediction, correlation effects are found to play a significant role. The local moment per Ru site remains finite ~2.3 μ B in the whole temperature range investigated. Due to the lower symmetry, the t 2g manifold is split and the quasiparticle weight is renormalized significantly in the a 1g state, while the renormalization in [Formula: see text] states is about a factor of 2-3 weaker. Our theoretical Neel temperature ~700 K is in reasonable agreement with experimental observations. SrRu2O6 is a unique system in which localized and itinerant electrons coexist with the proximity to an orbitally-selective Mott transition within the t 2g sector.

Introduction

For systems with an odd number of electrons per unit cell, correlation effects can lead to insulating magnetic ground states. Depending on the relative strength between the interactions U and the electron bandwidth or kinetic energy W, there are two classes of insulators: A Slater insulator[1] in the weak coupling regime in which the normal state is a non-magnetic metal. Below a critical temperature a gap in the single-particle opens up because of magnetic ordering and consequently the unit cell gets doubled and the Brillouin zone folds up. In this regime, the magnetic transition temperature increases with increasing Coulomb interaction strength. In the opposite strong coupling regime, the system is a Mott insulator[2,3]. In this limit, the normal state is a gapped insulator due to strong Coulomb repulsion. Magnetic ordering sets in below a critical temperature T  ~ W 2/U related to the superexchange scale[4]. Interpolating between these two limits, the maximum T is expected to occur in the crossover regime where U and W are comparable.

The behavior of T in a series of perovskite transition-metal oxides (TMOs) with the formal electron configuration, including SrMnO3, SrTcO3 and NaOsO3, lend support to the above picture. 3d TMO SrMnO3 has a relatively low Néel temperature T  = 260 K[5] for antiferromagnetic (AF) ordering but the insulating behavior persists above T , indicating that it is in the strong coupling limit. On the other hand, 5d TMO NaOsO3 with an AF transition T  = 410 K and a metal-insulator transition that is coincident[6] appears to be in the weak coupling limit. Finally, 4d TMO SrTcO3 shows an extremely high AF T  ~ 1000 K[7], which indicates it is located in the crossover regime[8]. It is also worth to note that the crossover behavior associated with high magnetic transition temperatures is realized in double perovskite TMOs A 2 BB′O6 [9], where A is an alkali or alkaline earth metal, and B and B′ are two transition metals with different Coulomb interaction strengths.

Given this backdrop, an important question that arises is can there be systems that are incompatible with this classification? And more broadly, what are the criteria for weak coupling Slater and strong coupling Mott regimes in multi-orbital systems?

One such systems is pyrochlore Cd2Os2O7. As in NaOsO3, the metal-insulator transition of Cd2Os2O7 accompanies a magnetic ordering[10,11]. However, unlike NaOsO3 the unit cell of Cd2Os2O7 contains four Os sites and therefore can become an insulator without changing the size of the unit cell. In fact, the observed all-in-all-out magnetic structure is compatible with the structural unit cell[12,13]. It turns out that the metal-insulator transition in Cd2Os2O7 is a Lifshitz transition[14] arising from a reconstruction of the Fermi surface, rather than a Slater transition.

In this paper, we focus on another compound, the hexagonal SrRu2O6 [15] [Fig. 1(a)]. Given a half filled t 2 configuration for Ru sites from its formal valence +5, like Os5+ in NaOsO3 or Cd2Os2O7, and also because Ru is a 4d element, one may expect that it would show behavior similar to SrTcO3. Indeed, a fairly high T  ~ 565 K is observed in SrRu2O6 [16]. However, the ordered moment ~1.3 μ is found to be much smaller than that expected for Ru5+, 3 μ [16,17]. To solve this puzzle, there have appeared a number of theoretical studies using density functional theory (DFT). Singh[18] and Tian et al.[16] argued that the strong hybridization between Ru and O ions is responsible for the reduction of the ordered moment. Streltsov and coworkers proposed the formation of molecular orbitals within a Ru hexagonal plane owing to such a strong hybridization[19]. They also performed the DFT calculations combined with the dynamical mean field theory[20] (DFT + DMFT[21-23]). It was reported that the ordered moment on a Ru site becomes 3 μ as anticipated for Ru5+ and that the transition temperature becomes ~2000 K, which is very similar to a theoretical result for SrTcO3 [8] and in return supports the importance of the molecular orbital picture that is not considered in DMFT.

Figure 1

(a) Top view of SrRu2O6. The cell is doubled along the a and b directions. (b) Ordered 〈M 〉 and equal-time spin-spin correlation as a function of temperature T. From the high temperature value, , the local moment is estimated to be M ~ 2.3 μ . (c) Schematic view of the localized vs. itinerant dichotomy between a 1 electrons and electrons. Correlation effects are stronger for a 1 electrons resulting in localized moments, while electrons maintain a rather itinerant character and show spin-density wave (SDW) like behavior when magnetically ordered.

In this paper, we investigate the electronic and magnetic properties of SrRu2O6 using DFT + DMFT. Our main finding is summarized in Fig. 1(b). In contrast to the previous report of Ref.[19], we find that the ordered moment is indeed ~1.2 μ , which is nearly identical to the DFT result 1.3 μ [18] and also consistent with the experimental value 1.4 μ [16,17]. These observations might indicate that SrRu2O6 is in the weak coupling limit, similar to NaOsO3, not in the crossover regime. However, the local moment estimated by the equal-time spin-spin correlation remains relatively large throughout the temperature range analyzed. We also find that the transition temperature estimated from DMFT is T  ~ 700 K, reasonably close to the experimental T  ~ 565 K [16,17]. Fluctuations not included in DMFT are expected to reduce the ordering temperature bringing it closer to the experimental estimate. Even more significantly, we find that one of t 2 states shows strong mass enhancement. Thus, our results indicate SrRu2O6 is in close proximity to an orbitally-selective Mott insulating state and localized and itinerant electrons coexist as schematically shown in Fig. 1(c).

Results

DFT analysis

We utilize DFT + DMFT techniques[21-23] to better account for correlation effects associated with the localized Ru d electrons. We first perform the DFT calculations for the non-magnetic state with the Elk package[24]. The experimental lattice parameters are considered with atomic configurations determined at room temperature[15]. It turns out that these calculations themselves provide important insights into the novel behavior of SrRu2O6. The spin-orbit coupling (SOC) λ is neglected in our analyses because λ ~ 200 meV in ruthenium oxides[25] is smaller than our estimation for the Hund coupling parameters. We anticipate that the main effect of the SOC is to introduce the uniaxial spin anisotropy.

Within the DFT, a small band gap ~0.05 eV exists at the Fermi level in a paramagnetic state[16], indicating that SrRu2O6 is a band insulator.

We then construct a model Hamiltonian consisting of Ru t 2 states and O p states using the Wannier functions[26-29] and transfer integrals between them on different sites extracted from the DFT band structure. We estimate the interaction parameters from the constrained random phase approximation (cRPA)[30]. Details of our model construction are provided in the Methods section. The interacting model thus obtained is solved by means of DMFT with the continuous-time quantum Monte-Carlo impurity solver[31-33].

It is instructive to start from the DFT results as summarized in Fig. 2. Figure 2(a) shows the DFT dispersion relations near the Fermi level. The valence-band maximum is set to E = 0. The first Brillouin zone for the hexagonal lattice is shown as an inset of Fig. 2(b). In this calculation, magnetic ordering is suppressed. RuO6 octahedra are slightly compressed along the c direction, and O-Ru-O angle is about 94.4° (two oxygen sites are on the same plane). Reflecting this distortion and the hexagonal symmetry, i.e., a 1 and are intrinsically inequivalent, the a 1 level is about 0.3 eV higher than the level in the Wannier basis.

Figure 2

DFT results on SrRu2O6. (a) Band dispersion (solid lines) compared with Wannier dispersions (broken lines). (b) Total and orbital-resolved (Ru , a 1 and O p) DOS. The valence-band maximum is set to E = 0. DOS for involves contributions from two degenerate bands. The inset shows the first Brillouin zone for the hexagonal lattice.

These bands are well separated from the other bands located below −15 eV or above +3 eV. As shown in Fig. 2(b), these bands primarily come from Ru t 2 states (split to twofold degenerate states and non-degenerate a 1 states) and O p states, justifying our choice of the basis set of the model Hamiltonian. It is clearly seen that Ru t 2 and O p states are strongly hybridized in the energy window [−7.0: +1.0] eV, which is used to construct Wannier functions in this study. Because of this large energy window, Wannier orbitals with the t 2 symmetry centered on a Ru site are well localized with the typical spatial spread 1.0–1.2 Å. Interestingly, a substantial amount of O p states are above the Fermi level, corresponding to 2.596 holes in the O p states per unit cell or 0.433 holes per oxygen ion. This point will become crucial in the discussion below. It should also be mentioned that Ru a 1 states have the largest weight near the Fermi level (Note that Ru states have the twofold degeneracy, therefore each band has the smaller spectral weight than the a 1 band near the Fermi level). This indicates that the effect of correlations is different between and a 1 states.

Figure 2 (a) also shows the Wannier dispersion relations, which completely overlap with the DFT dispersion relations. There are 48 bands including the spin degeneracy, i.e., 6 per Ru (t 2) and 6 per O (p , p , p ). A Sr does not contribute to these low-energy bands because Sr 5s(4p) level is so high (low) and its valence state is +2 by donating 2 electrons to these bands. Considering the nominal valence, Sr2+Ru2 5+O6 2−, and the charge counting, 3 per Ru5+ and 6 per O2−, these bands are filled by 42 electrons. In this study, the orbital occupancy is determined in the Wannier basis. This is also used to determine the valence state of a Ru ion and the magnetic moment. Within DFT with the Wannier basis, the filling is about 1.42 per a 1 orbital and about 1.45 per orbital, so these orbitals are nearly equally filled. Since the total t 2 occupation 4.3 is larger than 3, the Ru valence state is deviated from the nominal value +5 to +3.7.

Double-counting correction

For methods using DFT supplemented by many-body calculations, it is critical that we first correct for the Hartree contribution arising from correlated orbitals that is already included in the DFT calculation; we call this the double-counting correction (DC). This is introduced as a uniform level shift of correlated orbitals ε . Two methods for accounting for such a DC correction have been investigated during the development of local (spin) density approximation [L(S)DA] + U methods[34,35]. Since the combination of DFT + DMFT techniques started to appear, a number of DC correction methods have been proposed[36-41]. However, it remains unclear which of these methods is best suited for the current system. Here, instead of examining each DC method, we take a different approach: we treat ε as an adjustable parameter and, based on the experimental input, i.e., the ordered moment, determine the parameter range suitable for the material and then use it for the more detailed calculations of the quasi-particle weight.

The inset of Fig. 3(a [b]) shows the total (orbital-resolved) Ru d electron densities N(N ) as a function of ε . As one can see, the total density N remains larger than 3, the nominal value for Ru5+, for the whole ε range examined. To realize N = 3, i.e., the d 3 electron configuration, ε must be unrealistically small. N and N are not sensitive to the magnetic state, compare paramagnetic (PM) solutions and G-type antiferromagnetic (G-AFM) solutions. In all ε , a 1 orbitals have electron density closer to 1 than orbitals. From this, one could anticipate that a 1 electrons contribute to magnetism more strongly than electrons.

Figure 3

(a) DMFT results of the total and the orbital-resolved ordered magnetic moments, 〈M 〉 and 〈M 〉, respectively, as a function of ε . (b) Quasiparticle weight Z as a function of ε . The insets of (a) and (b) respectively show the total d electron density N and the orbital-resolved density N as a function of ε .

In the inset of Fig. 3(a), we also plot the analytical linewhich is the so-called fully localized limit (FLL)[34], with U and J being averaged Coulomb interactions and exchange interactions computed from matrix elements presented in the Method section, Eqs (4) and (5), respectively. We notice that the analytic curve and the numerical data practically overlap at . Thus, if one takes the FLL DC correction and imposes the condition that impurity occupation N and the one in the DC correction Eq. (1) coincide, any point at fulfills this condition. This means that the full convergence for this scheme is very difficult. Instead, we take the Ru d occupancy from the DFT calculation, N = 4.3, and use this value in Eq. (1) to obtain ε  = 15.92 eV. This ε gives impurity occupation N = 4.28(4.31) in the PM (G-AFM) solution at T = 348 K. It is reported that the FLL formula gives a slightly larger impurity occupancy than the “nominal” DC (in this case N = 3), which is close to the “exact” DC[41]. However, the true nominal occupancy in Ru d states is not known yet. While it remains to be justified, our results using the FLL DC with N = 4.3, as discussed below, suggest that our choice of ε provides a reasonable physical picture of SrRu2O6. It is very interesting to apply the improved DC correction to SrRu2O6, but this remains a future research.

Given the vital information on the Ru d occupation number N vs. ε in the inset of Fig. 3(a), we examine physical quantities as a function of ε . Figure 3(a) shows the total ordered moment 〈M 〉 and orbital-resolved ordered moment 〈M 〉 for G-AFM as a function of ε . Here, 〈M 〉 = ∑〈M 〉 with 〈M 〉 = N  − N . One notices that the ordered moment on the a 1 orbital is larger compared to the orbitals. This is because the a 1 occupation is closer to one than the occupations [see the inset of Fig. 3(b)] and, therefore, the spin polarization is easily induced.

It is remarkable that the ordered moment varies drastically with ε . This comes from the smaller occupation in a 1 orbitals than orbitals, i.e., < [see Fig. 3(b) inset]. This is partly explained by DFT results; it is found that the bare a 1 level is about 0.3 eV higher than the level, and = 1.42 < = 1.45. However, the orbital polarization   is found to be significantly enhanced in our DMFT calculations; from 0.03 in DFT to at most 0.25. This is because the orbital polarization is enhanced by correlations[42,43]. Furthermore, as shown later, spectral weight right below the Fermi level is dominated by a 1 orbital [Fig. 4(a)]. Since reducing ε roughly corresponds to lowering E , a 1 orbital changes the occupation number more sensitively than orbitals, leading to the large suppression in the ordered moment with decreasing ε .

Figure 4

DMFT results of the spectral function A (ω) for a PM phase (a) and a G-AFM phase (b) at T = 348 K. (c) Orbital-resolved quasiparticle wright Z as a function of T. The a 1 band shows a stronger mass enhancement than the band by a factor of 2–3. (d) Temperature dependent spectral function A (ω).

With ε  = 15.92 eV (the FLL DC correction with N = 4.3), we found the ordered moment . This value is consistent with the experimental observation 1.4 μ [16,17], giving us confidence in our choice of ε . This is quite different from the nominal occupation of Ru5+, N = 3. Such a large difference of N from its nominal value comes from the strong mixture between Ru d states and O p states in low-energy bands. It should be noted that in order to achieve N = 3, one needs to adopt unphysically small ε .

The strong sensitivity of 〈M 〉 on ε indicates the magnetism and the band structure are strongly coupled. This might suggest that SrRu2O6 is in the weak coupling regime despite relatively large values of the matrix elements of . But is this true? Thus, we now turn to the quasiparticle weight, which is the direct measure of the correlation strength. Figure 3(b) shows orbital dependent quasiparticle weight Z as a function of ε . It is seen that Z for α = a 1 becomes extremely small at ε  < 10 eV and ε  ~ 14 eV approaching integer fillings N = 3 and 4, respectively, while Z for α =  remains larger than 0.02. This indicates that the correlation effect is quite strong and the system is in the vicinity of orbital-selective Mott insulating regimes[44], in which α = a 1 electrons tend to be localized but electrons maintain an itinerant character. This dichotomy between localized and itinerant characters will be discussed in detail later. We also note that the coexistence of localized and itinerant electrons would lead to the bad-metal behavior that has been suggested in an extended range of coupling and carrier densities in strongly-correlated multiorbital systems in the presence of Hund’s coupling[45,46].

In many magnetic materials, the ordered moment increases with increasing correlation strength. Thus, the quasiparticle weight in a PM solution and the ordered moment in a magnetic solution tend to have a negative correlation. As shown in Fig. 3(a) and (b), on the contrary, the ordered moments become large where the quasiparticle weights become large, i.e., they have a positive correlation. All these observations suggest that SrRu2O6 is in the strong coupling regime, but Mott physics and itinerant band physics coexist.

Temperature dependence

Based on the detailed comparison with experiments, now we focus on ε = 15.92 eV and examine the temperature dependence of the electronic property of SrRu2O6.

We start from the low-temperature electronic property. Figure 4(a,b) show the Ru d density of states computed by the maximum entropy analytic continuation for the impurity Green’s functions on the Matsubara axis[47]. states have a clear gap in both PM and G-AFM solutions with the peak-to-peak distance ~0.6 eV and ~0.9 eV, respectively. For a 1 states, the gap structure is rather vague especially in the PM solution. The minimum peak-to-peak distance is ~0.2 eV for both PM and G-AFM. These values are comparable to the experimental estimate on the activation energy ~0.4 eV[16].

Magnetic properties are summarized in Fig. 1(b). Local ordered moment 〈M 〉 is computed for G-type AFM and plotted as a function of temperature. At a glance, one notices a close resemblance between the theoretical 〈M 〉 vs. T curve and experimental ones reported in refs[16] and [17]. Theoretical results for the ordered moment at low temperatures is and the Néel temperature T  ~ 700 K are reasonably close to the experimental results, and T  ~ 565 K. Theoretical T is overestimated by about 20%, which could be attributed to the mean field nature of the current “single-site” DMFT and the uniaxial anisotropy existing in real material due to the finite spin-orbit coupling but absent in our calculations.

In Slater-type systems in the weak coupling limit, local moments disappear when magnetic ordering disappears above T . On the other hand, in Mott insulators in the strong coupling limit, local moments remain unchanged above T . In the intermediate crossover regime, local moments decrease with increasing temperature but survive at high temperatures. We have also computed the equal-time spin-spin correlation in an impurity model by suppressing magnetic ordering and plotted in Fig. 1(b). It is remarkable that remains constant in the whole temperature range. From and using the spin rotational symmetry, the size of the local moment is deduced as , which is roughly 2 times larger than the size of the ordered moment at low temperatures. Clearly SrRu2O6 is not in the weak-coupling limit and most likely in the strong-coupling limit. Similar dichotomy between large local moments and small ordered magnetic moments was theoretically suggested for iron-based superconductors[48,49] and has indeed been experimentally reported[50].

We also found other evidence of strong correlation effects from the temperature dependent quasiparticle weight. As shown in Fig. 4(c), Z in PM solutions are decreased with decreasing T. This behavior contradicts with what one expects for a good metal where the electron coherence is enhanced with decreasing temperature. It is also interesting to point out that the renormalization in is moderate for high temperatures. This is because bands maintain the band-insulator-like character and have larger spectral weights away from ω = 0. On the other hand, the renormalization in Z is significant. Thus, SrRu2O6 is in the vicinity of the correlation-induced orbital-selective Mott transition.

Our findings are remarkably different from the previous DFT + DMFT study on SrRu2O6. Reference[19] found a large Ru moment ~3 μ and an extremely high transition temperature ~2000 K. While sufficient details are not provided in ref.[19], we suspect from our calculations that the main difference arises from the inclusion of ligand O p states in our calculations that are crucial to describe the electronic properties of the late transition-metal oxides with deeper d levels.

As a closely related system, SrTcO3 was studied in ref.[7]. The experimental ordered moment is found to be 2.13(1.87) μ at low (room) temperature. The ordered moment deduced by DFT is 1.3 μ . All these values are substantially smaller than S = 3/2 moment 3.87 μ . SrTcO3 as well as SrMnO3, both of which are d 3 systems, were examined using DFT + DMFT in ref.[8]. The ordered moment is found to be 2.5 μ for SrTcO3 and 3 μ for SrMnO3 at low temperatures. These values are close to the experimental values. The local moment for the two materials are estimated from above T to be 2.7 μ for the Tc and 3.8 μ for Mn compound. The fluctuating magnetic moment of Mn is close to the maximal value 3.87 μ , while that for Tc is suppressed because of the more-itinerant nature of SrTcO3. The reported dependence of ordered and local moments on the “B” site elements of perovskite SrBO3 is consistent with the Slater-Mott crossover for half-filled Hubbard models. Our results on SrRu2O6 are closer to SrTcO3 than SrMnO3. However, the reduction from the local moment to the ordered moment is much stronger 1.3/23 ~ 0.57 for our SrRu2O6 as compared with 2.5/2.7 ~ 0.93 for SrTcO3. This clearly indicates the importance of the coexistence of localized and itinerant electrons in SrRu2O6. It would be also very instructive to see the difference between SrRu2O6 and other ruthenium compounds with the nominal Ru valence of +5. Double perovskite ruthenates Ba2LaRuO6 (T  = 29.5 K), Ca2LaRuO6 (T  = 12 K)[51,52] and Sr2YRuO6 (T  = 26 K)[51,53] are among such compounds. The high temperature fluctuating (low temperature ordered) moment of Ba2LaRuO6, Ca2LaRuO6 and Sr2YRuO6 was deduced as 4.00, 4.27 and 3.13 μ (1.96, 1.92 and 1.85 μ ), respectively. These values are factor 1.4–1.8 larger than the corresponding values obtained in this study for SrRu2O6. The difference between these double perovskite ruthenates and SrRu2O6 can be ascribed to the more localized nature of the double perovskite ruthenates, in which Ru sites are separated by rare earth ions, and the Ru valence expected to be much closer to +5.

Discussion

In the previous sections, we discussed magnetic and electronic properties of SrRu2O6. Here we discuss some implications of our results to the transport property. Because of the large imaginary part of the self-energies, it remains challenging to identify the low-energy feature of the spectral function. Nevertheless, the temperature dependence of A (ω) might give a hint. As shown in Fig. 4(d), bands show the gapped feature and, therefore, do not contribute to the transport much. On the other hand, the spectral function for the a 1 band is gradually diminished at ω ~ 0 with decreasing temperature. Thus, we expect the electric resistivity ρ(T) to increase with decreasing temperature in the whole temperature range, including T  ≫ T .

In the weak-coupling picture, SrRu2O6 is a band insulator at both above and below T , the gap amplitude is fixed at T > T and grows rapidly according to the magnetic order parameter at T < T . Therefore, the ρ(T) curve is expected to have a kink at T . On the other hand, in our strong coupling picture, the spectral functions evolve dynamically even in the PM regime. Due to the large self-energy or scattering rate, details of low-energy features in A(ω) are smeared out. Thus, ρ(T) is expected to change smoothly across T .

To verify our results, it would be very important to experimentally examine the resistivity in the whole temperature range, from T ≪ T to T ≫ T . Our preliminary results of the resistivity above room temperature appear to support this. Also, (angle resolved) photoemission spectroscopy measurements would provide direct evidence of the dynamical change in the spectral function, the evolution of the pseudogap feature to the full gap with decreasing temperature across T .

The current work suggests that SrRu2O6 has to be considered as a different class of material than a series of TMOs, including SrMnO3, SrTcO3 and NaOsO3, where the low-energy behavior is well captured by multi(3)-orbital Hubbard models. To understand the unique behavior of SrRu2O6, it is essential to consider strong hybridization between Ru d and O p states. This could largely modify the T vs. U phase diagram by reducing both T and the local moment M , while the correlation strength in TM d states remains strong. There might be a large number of materials that could be classified into a similar category as SrRu2O6.

The main effect of the hybridization with O p states is driving the electron configuration of Ru to rather than . In this respect, it would be very interesting to look for systems with the similar electron configurations as SrRu2O6 but stronger SOC. As discussed in ref.[54], the d 4 electron configuration would nominally imply an atomic moment J = 0, but, interatomic exchange can lead to the formation of a local moment and novel spin-orbital entangled magnetism.

To summarize, we have investigated the novel electronic and magnetic properties of the hexagonal compound SrRu2O6 using DFT + DMFT. The small moment on Ru ~1.2 μ and the relatively high Néel temperature ~700 K are consistent with the experimental report, ~1.4 μ and T  ~ 565 K. These seemingly contradictory characters of SrRu2O6 are caused by the strong hybridization between Ru d states and ligand O p states, which increases the Ru d occupancy substantially from the nominal value 3. The strong Ru d-O p hybridization does not imply a weak-coupling nature. In fact, SrRu2O6 is in the strong-coupling regime, on the verge of the orbital-selective Mott localization. In contrast to the DFT results, a band insulator when magnetic ordering is suppressed, we predict that the electron spectral functions evolve dynamically with temperature. In particular above Néel temperature, the spectral functions exhibit pseudogap features, which turn to full gap with decreasing temperature. To verify our predictions, further experimental studies are desirable, including temperature dependent photoemission spectroscopy measurements and resistivity measurements.

Methods

DFT calculations were performed for the non-magnetic state with the Elk package[24] using the exchange-correlation functional proposed by Perdew et al.[55]. We considered the experimental lattice parameters and atomic configurations determined at room temperature[15].

Using the obtained band structure, we constructed the Wannier functions[26-29] containing the Ru t 2 and O p orbitals and calculated the transfer integrals among them. The energy window was set as [−7.0:+1.0] eV. Further, we evaluated the interaction parameters, the Coulomb repulsion and the exchange coupling , by the constrained random phase approximation (cRPA)[30]. In the calculation of the partially screened Coulomb interaction, we took 120 unoccupied bands and used a 4 × 4 × 4 grid. The double Fourier transform of the constrained susceptibility was done with the cutoff of 5 (1/a.u.). We neglected the SOC since it is small compared with the Hund coupling constants.These calculations were performed using a density response code[56] recently developed for the Elk branch of the original EXCITING FP-LAPW code[24].

The obtained effective model consists of two parts as H = H  + ∑ H . The band part H is given byHere, () is a vector consisting of annihilation operators of Ru t 2 (O p) electrons with spin σ in a momentum space. The dispersion consists of Wannier parameters, with α and β running through Ru t 2 and O p states. ε for interacting Ru t 2 states is the DC correction.

The interaction part H is given byHere, the index i for Ru sites is suppressed for simplicity. The matrix elements of and , evaluated by using the cRPA, are given byandin units of eV with the basis of {e′, e′, a 1}.

Recently, we become aware of a preprint by Hariki et al. (Ref.[57]). They performed similar DFT+DMFT calculations with different interaction parameters and a DC scheme and obtained similar conclusion as Ref.[19].