Fermi surface geometry and momentum dependent electron-phonon coupling drive the charge density wave in quasi-1D ZrTe$3$

ZrTe$_3$ is a prototypical quasi-one-dimensional compound undergoing a charge density wave transition via a very sharp Kohn anomaly in phonon momentum space. While Fermi surface geometry has long been considered the primary driver of the instability, a full understanding of the lattice dynamics and electron-phonon role has remained elusive. Our first principles calculations in the high-symmetry phase show that the Fermi surface is correctly reproduced only when the Hubbard interaction on the Te $5p$ orbitals is included, which in turn is essential for the appearance of a soft harmonic phonon mode at the CDW wavevector. Analyzing the mode and momentum dependence of the electron-phonon coupling, we find that its variations with phonon momentum dominate over electronic effects. These results identify unambiguously the CDW origin in ZrTe$_3$ as a cooperative effect of Fermi surface geometry and momentum-dependent electron-phonon coupling, with the latter playing the leading role. The mechanisms revealed in our work are directly relevant to other quasi-1D systems, including trichalcogenides and compounds hosting Peierls-like chains.

💡 Research Summary

ZrTe₃ is a prototypical quasi‑one‑dimensional (quasi‑1D) material that undergoes a charge‑density‑wave (CDW) transition at 63 K and becomes superconducting below 2 K. While earlier studies have emphasized Fermi‑surface (FS) nesting as the primary driver of the CDW, the detailed role of lattice dynamics and electron‑phonon interaction (EPI) remained unclear. In this work the authors present a comprehensive first‑principles investigation of the high‑temperature (HT) phase, combining density‑functional theory (DFT), density‑functional perturbation theory (DFPT), and Wannier‑function techniques.

Standard PBE calculations, even when spin‑orbit coupling (SOC) is included, fail to reproduce the ARPES‑observed quasi‑1D FS sheets: the Te 5pₓ bands are overly delocalized and the resulting FS topology does not match experiment. By computing on‑site Hubbard U parameters for the Te 5p orbitals via a linear‑response approach (U≈4.3 eV) and performing DFT+U calculations, the authors obtain a markedly improved electronic structure. The corrected FS consists of two nearly parallel, flat sheets derived from the Te(2)–Te(3) 5pₓ bands, in excellent agreement with ARPES data and with the experimentally determined CDW wave vector q_CDW = (0.07, 0, 0.33) r.l.u.

Having established a realistic FS, the nesting function ζ(q) is evaluated along the (h, 0, 1/3) direction. A sharp peak appears at h ≈ 0.07 r.l.u., confirming that the FS geometry provides a nesting condition that selects the CDW wave vector. However, the real part of the non‑interacting electronic susceptibility χ₀(q) remains essentially flat across the same momentum range, showing only a modest softening when the full Brillouin zone is considered. This indicates that the electronic instability associated with pure FS nesting is weak and not sufficient to drive the transition.

The lattice dynamics are then examined with DFPT+U. A harmonic phonon mode softens dramatically at q_CDW, reproducing the experimentally observed narrow‑in‑q Kohn anomaly. The electron‑phonon linewidth γ(q) is decomposed to reveal the momentum dependence of the electron‑phonon matrix elements gₙₘ(k,q). The analysis shows that variations of gₙₘ(k,q) with q dominate over the nesting‑derived contributions (ζ and χ₀) by a factor of two or more. In other words, the strength of the EPI is strongly q‑dependent, and this enhancement at q_CDW is the primary driver of the lattice distortion and CDW formation. Inclusion of SOC does not materially alter these conclusions, confirming that spin‑orbit effects are secondary in this material.

The study thus arrives at three central conclusions: (i) an accurate description of the Te 5p electronic correlations via Hubbard U is essential for reproducing the correct FS; (ii) FS nesting supplies the correct CDW wave vector but does not generate a pronounced electronic susceptibility peak; (iii) the momentum‑dependent electron‑phonon coupling is the decisive factor that actually triggers the CDW transition. The authors argue that this cooperative mechanism—FS geometry setting the wave vector and a sharply q‑dependent EPI providing the driving force—is likely relevant to other quasi‑1D systems such as transition‑metal trichalcogenides, blue bronzes, and platinum‑chain compounds. Their work provides a clear framework for disentangling nesting and electron‑phonon effects in low‑dimensional CDW materials.