Theory of unconventional magnetism in a Cu-based kagome metal

Kagome metals have established a new arena for correlated electron physics. To date, the predominant experimental evidence centers around unconventional charge order, nematicity, and superconductivity, while magnetic fluctuations due to electronic interactions, i.e., beyond local atomic magnetism, have largely been elusive. We find the challenge of locating the appropriate parameter regime for such exotic order to center around two aspects. First, the correlations implied by low-energy orbitals have to be sufficiently large to yield a dominance of magnetic fluctuations and weak to retain an itinerant parent state. Second, the kinematic kagome profile at the Fermi level demands an efficient mitigation of sublattice interference causing the suppression of magnetic fluctuations descending from electronic on-site repulsion. We elucidate our methodology by analyzing the potential copper-based kagome compound CsCu$_3$Cl$_5$: From ab initio design and many-body analysis, we develop a model framework of realistic Cu-based kagome materials the simulations of which reveal unconventional magnetic order in a kagome metal.

💡 Research Summary

The manuscript investigates the emergence of unconventional magnetic order in a copper‑based kagome metal, using CsCu₃Cl₅ as a concrete example. The authors begin by highlighting that most kagome metals studied to date (e.g., the AV₃Sb₅ family) display charge‑order, nematicity, and superconductivity, while magnetic fluctuations driven by electronic correlations have remained largely absent. They argue that to realize correlation‑driven magnetism one needs (i) sufficiently strong low‑energy correlations to dominate magnetic fluctuations yet remain weak enough to preserve an itinerant metallic state, and (ii) a mitigation of the sublattice‑interference (SI) mechanism that suppresses on‑site Hubbard scattering at the kagome van Hove singularities.

First‑principles density‑functional calculations reveal that CsCu₃Cl₅ crystallizes in the P6/mmm space group, forming a layered Cu₃Cl₅ kagome sheet separated by Cs layers. The Cu atoms adopt a mixed valence (+1.33 on average) with both Cu(I) (d¹⁰) and Cu(II) (d⁹) character. The band structure shows three isolated kagome‑derived bands near the Fermi level: a flat band at the bottom, a Dirac point at K, and, crucially, an m‑type van Hove singularity (vHS) at the three inequivalent M points. Unlike the p‑type vHS of AV₃Sb₅, where each vHS is localized on a single sublattice, the m‑type vHS distributes electronic weight over two sublattices while the third remains empty. This “mixed” vHS weakens SI, allowing on‑site interactions to act more efficiently.

To quantify electronic correlations the authors perform constrained random‑phase approximation (cRPA) calculations on the low‑energy Wannier subspace. They obtain a screened on‑site Hubbard U ≈ 3.6 eV (bare U ≈ 14.7 eV) and a nearest‑neighbor repulsion V of comparable magnitude. With a bandwidth t extracted from the tight‑binding fit, the ratio U/t ≈ 4.7 places CsCu₃Cl₅ in an intermediate regime—stronger than the weakly correlated AV₃Sb₅ (U/t < 1) but far below the Mott‑localized Herbertsmithite (U/t ≈ 20). This intermediate coupling is identified as the “sweet spot” where magnetic fluctuations can become dominant without destroying metallicity.

Maximally localized Wannier functions are constructed from a linear combination of Cu dₓz/d_yz and Cl p orbitals, yielding three identical orbitals centered on the kagome sites and related by 60° rotations. The resulting three‑orbital tight‑binding model reproduces the DFT bands with high fidelity, providing a minimal Hamiltonian for many‑body analysis.

The interacting Hamiltonian includes both on‑site Hubbard U and nearest‑neighbor density‑density V: H_int = U ∑i n{i↑} n_{i↓} + V ∑{⟨i,j⟩,σ,σ′} n{iσ} n_{jσ′}. Because the m‑type vHS reduces SI, both U and V contribute on similar energy scales, unlike the p‑type case where V can often be neglected.

Functional renormalization group (FRG) calculations are performed using the truncated‑unity implementation (divERGe code). The FRG flow integrates out high‑energy modes, renormalizing the effective two‑particle interaction. The leading instability emerges in the spin‑bond channel: a non‑magnetic but spin‑dependent ordering pattern that can be interpreted as a bond‑centered spin density wave or a loop‑current state. This order is driven primarily by the enhanced nesting at the m‑type vHS and the comparable strength of U and V, which together overcome the residual SI suppression. The calculated phase diagram shows that modest variations of U/t or V/t can switch the system between a paramagnetic metal, the spin‑bond ordered phase, and, at stronger coupling, a more conventional spin‑density‑wave state.

The authors compare their findings with AV₃Sb₅, where phonon‑mediated charge order dominates due to strong SI, and with Herbertsmithite, where strong Mott physics leads to spin‑liquid behavior. CsCu₃Cl₅ occupies a distinct niche: the m‑type vHS provides a momentum‑space localization that reduces electron‑phonon coupling, making electronic correlations the primary driver of ordering. This suggests that Cu‑Cl based 135‑type kagome metals can host a new class of correlation‑driven magnetic phenomena absent in previously studied kagome compounds.

In conclusion, the paper presents a comprehensive workflow—from crystal‑structure design, through ab‑initio electronic structure and interaction extraction, to many‑body FRG analysis—that identifies CsCu₃Cl₅ as a realistic material platform for unconventional magnetism in a metallic kagome system. It opens avenues for experimental verification (e.g., neutron scattering, muon spin rotation) and for extending the methodology to other transition‑metal kagome compounds where mixed valence and tailored van Hove singularities could be engineered to realize exotic magnetic orders.