Effect of Interlayer Stacking on the Electronic Properties of 1$T$-TaS$_2$

Controlled stacking of van der Waals materials is a powerful tool for exploring the physics of quantum condensed matter. Given the small binding between layers, exploitation for engineering will require a breakthrough in stacking methodology, or an ability to take advantage of thicker defective stacks. Here we describe computational groundwork for the latter, using – on account of its promise for cold memory applications – 1$T$-TaS$_2$ as a model system. Comparing recursive Hendricks-Teller calculations and Monte Carlo simulations to published X-ray diffraction data, we obtain the key parameters describing the random stacking in mesoscopic flakes. These then regulate the electronic structures via specification of the random stacks in dynamical mean-field theory simulations. Hubbard repulsion induces strongly correlated metallic, band and Mott insulating layers, providing compelling evidence that electronic properties follow from the coexistence of more than the metallic and insulating planes associated by ordinary band theory.

💡 Research Summary

This paper investigates how random interlayer stacking in the low‑temperature commensurate charge‑density‑wave (CCDW) phase of 1T‑TaS₂ controls its electronic properties. The authors first quantify the stacking disorder by combining recursive Hendricks‑Teller (HT) calculations with Monte Carlo (MC) simulations and compare the resulting structure factors to published X‑ray diffraction (XRD) data. By fitting the out‑of‑plane Bragg peaks, especially the half‑ordered “dimer” peaks around l ≈ 15 r.l.u., they extract the probabilities of the three possible stacking configurations (Tₐ, T_b, T_c). The analysis shows that the CCDW state is dominated by T_c stacking, with a mixture of Tₐ‑stacked dimers and monolayers in an approximate 2 : 1 ratio (≈ 23 % dimers). The T_b configuration is essentially absent because its inclusion would shift the diffraction peaks away from experiment.

Using these statistically realistic stacks, the authors construct a minimal single‑orbital tight‑binding model for each layer, where the orbital corresponds to the half‑filled d_{z²} band localized on the central Ta atom of each polaron star. Interlayer hopping t_v is obtained from density‑functional theory (DFT) for Tₐ and T_c configurations, while a Hubbard repulsion U = 0.2 eV is added to capture strong electron correlations. Dynamical mean‑field theory (DMFT) is then applied layer‑by‑layer (real‑space DMFT) to compute the local self‑energies and spectral functions A_n(ω) at T ≈ 30 K.

The DMFT results reveal three distinct electronic behaviors depending on the local stacking environment: (i) Dimerized Tₐ layers are band‑insulating; the non‑interacting band structure already shows a ~0.4 eV hybridization gap, and correlations have little effect. (ii) Isolated monolayers surrounded by dimers become Mott insulators, as indicated by a divergent low‑frequency self‑energy and a correlation‑driven gap that appears only when U is included. (iii) Sequences of monolayers connected by T_c stacking develop a strongly correlated metallic state, with a finite quasiparticle weight at the Fermi level despite the presence of a gap‑like feature in the non‑interacting DOS.

Averaging over ten representative 20‑layer stacks (each with the 2 : 1 dimer‑to‑monolayer ratio) yields an overall spectral function A(ω) that retains spectral weight at the Fermi level, confirming that the bulk material is not a simple insulator but a heterogeneous mixture of band‑insulating dimers, Mott‑insulating monolayers, and correlated metallic layers. The authors also note that the degree of dimerization correlates with the out‑of‑plane diffraction peak positions reported in different XRD studies, suggesting that cooling rate or sample preparation can tune the stacking disorder and thus the electronic response.

In summary, the paper provides a comprehensive workflow: (1) extract stacking probabilities from XRD using HT and MC methods; (2) generate realistic multilayer configurations; (3) feed these configurations into DMFT to obtain layer‑resolved electronic spectra. The key finding is that interlayer stacking, rather than merely in‑plane charge order, dictates whether a given layer behaves as a band insulator, a Mott insulator, or a correlated metal. This insight resolves longstanding controversies over the nature of the low‑temperature insulating state of 1T‑TaS₂ and points to the possibility of buried Mott‑insulating layers that could host quantum spin‑liquid physics. The methodology is broadly applicable to other van‑der‑Waals materials where weak interlayer coupling coexists with strong electronic correlations, offering a powerful bridge between diffraction experiments and many‑body electronic theory.