Coulomb Correlations and Electronic Structure of CuCo2S4: a DFT + DMFT Study

By employing a combined method of density functional theory and dynamical mean-field theory (DFT + DMFT) we study the spectral properties and local moment formation in the normal state of a Co-based superconductor CuCo2S4. The obtained quasiparticle mass enhancement \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m{\text{*/}}m \sim 1.2$$\end{document} is smaller than that in iron-based parent compounds indicating weak correlations in the Co \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3d$$\end{document} states. We show that the mass enhancement is accompanied by an insignificant renormalization and shift of the bands in the vicinity of the Fermi level. This subtle modification is sufficient to induce a substantial transformation of one of the Fermi surface sheets and occurs due to the presence of shallow Fermi surface pockets in the electronic structure of CuCo2S4. The calculated local and fluctuating moments are smaller than those in iron-based superconductors questioning the importance of spin fluctuations for understanding the pairing mechanism of CuCo2S4.

INTRODUCTION Superconductivity (SC) of Co-based compounds occurs extremely rare compared to other transition metal counterparts. For a long time, the only known example of a Co-based SC was the oxide hydrate Na x CoO 2 ⋅yH 2 O ( K) [1,2]. Later, SC was found in LaCo 2 B 2 ( K) [3], Lu 3 Co 4 Ge 13 ( K) [4], and very recently in CuCoSi ( K) [5]. Also, it is worth noting that there is another exotic Co-containing system PuCoGa 5 exhibiting a rather high 18 K [6,7] however its SC stems from Pu electrons.
The thiospinel compound CuCo 2 S 4 shows many similarities with Na x CoO 2 ⋅yH 2 O including Co coordination, Co formal valence, and geometric frustration. First report on SC of Cu 1 + x Co 2 -x S 4 multiphase powder samples with an estimate K is dated back to 1993 [8]. Later, the SC-like behavior was confirmed by nuclear magnetic resonance (NMR) measurements of a multiphasic Cu 1.5 Co 1.5 S 4 revealing the growth of antiferromagnetic spin correlations at low temperatures [9]. Recent magnetic and transport measurements of the new generation of nearly stoichiometric CuCo 2 S 4 samples detected Pauli paramagnetism and a SC transition at K [10]. Investigations of the symmetry of the SC gap of CuCo 2 S 4 arrived at different conclusions. Early NMR measurements suggested a gapless SC state developing in line with an enhancement of antiferromagnetic fluctuations [9]. By contrast, later NMR experiments indicated a full isotropic SC gap yet confirmed the presence of weak antiferromagnetic spin correlations [11,12]. In this context, it is interesting that the order parameter of the related Co-based system Na x CoO 2 ⋅yH 2 O is presumably anisotropic [13][14][15]. Thus, it is not clear whether CuCo 2 S 4 should be classified an unconventional SC [16].
Until now, there are only few research papers reporting microscopic analysis of the electronic properties of CuCo 2 S 4 using standard density functionalbased methods neglecting the effect of electron-electron correlations in the partially filled Co shell [17,18]. At the same time importance of electron correlations in CuCo 2 S 4 is traced from resistivity measurements as indicated by the authors of [19]. Therefore, a careful theoretical study of the interplay between the Coulomb interaction and the electronic structure and magnetism in the normal state of CuCo 2 S 4 would be very helpful for understanding the anomalous properties of this material and can be used in further works to clarify the nature of its SC.
In this Letter, we explore the effect of Coulomb correlations on the electronic structure and magnetic properties of CuCo 2 S 4 employing the DFT + DMFT method [20,21]. We analyze the impact of electronelectron interaction on the spectral functions, Fermi surface (FS) shape, quasiparticle effective mass renormalization , and formation of the local

CONDENSED MATTER
moments. Our results demonstrate that CuCo 2 S 4 is a weakly correlated compound with a high degree of electronic itinerancy. Nevertheless, it is found that correlation effects are strong enough to cause a significant modification of one the FS sheets.
COMPUTATIONAL DETAILS In our work we use the DFT + DMFT computational scheme fully self-consistent in charge density [22,23]. This approach merges the many-body dynamical mean-field theory of correlated electrons (DMFT) [24,25] and density functional theory (DFT) and has proven to be a powerful theoretical tool for investigation of the spectral, structural and magnetic properties of real materials including the normal state of unconventional Fe-based superconductors [26]. We adopt the experimental crystal structure of CuCo 2 S 4 [27] (space group ) and employ the pseudopotential based generalized gradient approximation within DFT [28,29] to calculate the noninteracting band structure. The basis set of atomiccentered symmetry-constrained Wannier functions is used to construct the effective Hamiltonian in DMFT. The Wannier functions are computed by the projection procedure [30,31] for the energy window including Co , S , and Cu states. The DMFT realistic many-body problem is solved by the continuoustime quantum Monte Carlo (CT-QMC) method [32,33] using the on-site Hubbard parameter eV and Hund's rule coupling eV for the partially filled Co states of CuCo 2 S 4 . In DFT + DMFT we employ the fully localized limit for the double-counting correction and use the Padé approximant method [34] for analytical continuation of the self-energy.

RESULTS AND DISCUSSION
We start with an analysis of the effect of Coulomb correlations on the spectral properties. To this end, we first compare the spectral functions of CuCo 2 S 4 calculated by DFT with those obtained within DFT + DMFT (Fig. 1). Our DFT results (Fig. 1, shaded  areas) show that the total spectral function corresponds to a metal and features a broad peak at the Fermi energy ( ). This peak originates from the almost fully occupied Co orbitals (filling per spin-orbital) strongly hybridized with the S states. By contrast, the Co states are close to half-filling ( ), however their spectral function exhibits a dip in the vicinity of . We also note that the Cu orbitals show a minor contribution to the spectral weight close to forming a broad band from -5 to -1.5 eV. These results imply that the effect of local Coulomb correlations on the spectral properties of CuCo 2 S 4 is presumably not significant due to Indeed, we find that the spectral functions of CuCo 2 S 4 computed by DFT + DMFT (Fig. 1, solid  lines) show little difference with those obtained by DFT. In particular, we observe that correlation effects only slightly shift the peaks of the spectral functions below the Fermi level leaving their shape in the vicinity of almost unchanged. To give a quantified characterization of the effect of Coulomb correlations we calculate the orbitally resolved enhancement of the quasiparticle mass, = where is the diagonal element of the self-energy. We obtain for the orbitals and for the states. These values accord to the observed small renormalization of the DFT spectral functions by correlation effects accompanied by an insignificant reduction of the quasiparticle spectral weight in the vicinity of . We also note that these renormalization factors are approximately twice smaller than characteristic values of obtained in DFT + DMFT calculations for different Fe-based unconventional superconductors [35,36] and resemble those in the classical itinerant magnet ZrZn 2 [37].
To investigate the effect of Coulomb correlations on the spectral properties of CuCo 2 S 4 in more detail we compute the Fermi surface and the momentumresolved spectral function . We found that within DFT the FS exhibits three sheets (Fig. 2, lower panel). The first (inner) and the second hole-like sheets represent closed geometrical shapes centered at and enveloping the point of the fcc Brillouin zone (BZ). The first FS to a good approximation can be considered a sphere while the second FS features "outgrowths" directed to the X points of the BZ. The third (outer) FS exhibits a complex shape that intersects the BZ border.
Unlike the integral spectral functions our DFT + DMFT calculations show significant correlation induced transformations of one of the FS sheets (Fig. 2, upper panel). In particular, we observe an increase of the radius of the inner sheet and a stretch of the middle sheet along the Γ-L direction. These parts of the FS however exhibit a weak sensitivity to correlation effects and mostly preserve their shape in DFT + DMFT. On the contrary, the outer sheet now represents a closed shape that fully fits the BZ.
The above mentioned modifications of the FS can be understood from a comparison of the DFT + with the energy bands from DFT (Fig. 3). We note that due to the small quasiparticle damping (both Im and Im eV) the spectral function describes a coherent distribution of the spectral weight close to that allows an unambiguous interpretation of the DFT + DMFT results. One can see that the effect of Coulomb correlations on the electronic structure primarily corresponds to a shift and renormalization of the bands in close similarity to the situation in parent compounds of Fe-based superconductors [38,39]. In CuCo 2 S 4 however due to the presence of shallow pockets along the W-L-K path these transformations induce a significant change of the FS shape. Specifically, we observe that the DFT bands crossing the Fermi level along the W-L-K path and one band in the vicinity of the K point on the Γ-K direction shift below in DFT + DMFT. As a result, the outer FS sheet loses its intersection with the W-L-K plane of the BZ and transforms into a closed shape. The observed transformations of the FS presumably indicate the need of taking into account the electronic correlations for interpretation of the temperature behavior of the resistivity of Cu-Co-S superconducting samples as noted in [19] and require a thorough study.
As the final step of our study of the electronic properties of CuCo 2 S 4 we analyze the local moment formation. Namely, our DFT + DMFT calculations give the fluctuating magnetic moment . This value is significantly (by a factor ~2.5) smaller than that reported for parent compounds of pnictide and chalcogenide Fe-based superconductors [40,41] which show well-defined local moments and therefore implies a weaker localization of electrons in = 290 T CuCo 2 S 4 . To illustrate the higher degree of electronic itinerancy we compute the local spin correlator (where is the instantaneous Co spin at imaginary time computed by CT-QMC) and estimate the local moment M loc = . Our calculations give 0. 18 which is at least two times smaller than that obtained by DFT + DMFT for the Fe-based pnictides and chalcogenides [40,41]. The small magnitude of both and (compared to those in Fe-based superconductors) is indicative of weak spin fluctuations. This result questions the spin-fluctuationmediated origin of the superconducting state in CuCo 2 S 4 . We note however that absence of strong spin fluctuations due to correlation effects does not necessarily exclude importance of the latter for understanding the pairing mechanism. For example, it is known that not only Coulomb correlations affect the phonon anomalies [42] but also significantly enhance the electron-phonon coupling [43] leaving the possibility to explain the superconductivity of CuCo 2 S 4 within the standard Bardeen-Cooper-Schrieffer theory.

CONCLUSIONS
In conclusion, using the DFT + DMFT method we have investigated the effect of Coulomb correlations on the spectral properties and local moments of CuCo 2 S 4 . Our calculations have shown that CuCo 2 S 4 is a weakly correlated material characterized by the quasiparticle effective mass enhanced by a factor of . We have demonstrated that the mass enhancement is accompanied by a weak renormalization and shift of the DFT spectral functions and energy bands. Surprisingly, we have found that these shifts induce a significant transformation of one of the FS sheet due to the presence of shallow FS pockets along the W−L−K path indicating the possible importance of correlation effects for transport properties of CuCo 2 S 4 . Finally, our results for the local and fluctuating magnetic moments allow us to speculate that the spin-fluctuation scenario of unconventional pairing is unlikely applicable to explain the pairing mechanism in CuCo 2 S 4 . FUNDING This work was supported by the Russian Science Foundation (project no. 19-12-00012).