Robust Orbital-Selective Flat Bands in Transition-Metal Oxychlorides

Abstract

Flat electronic bands, which amplify electron correlations by quenching kinetic energy, provide an ideal foundation for exotic quantum phases. However, prevailing strategies -- including geometrically frustrated lattices, moire superlattices and heavy-fermion physics -- suffer from inherent trade-offs among robustness, tunability and orbital selectivity, limiting their broad applicability. Here, we unveil an intrinsic orbital-selective flat-band mechanism in the van der Waals materials NbOCl2 and TaOCl2, directly observed by angle-resolved photoemission spectroscopy (ARPES) and understood through density functional theory (DFT) and Wannier analysis. Crucially, we experimentally demonstrate that this momentum-independent flat band exhibits remarkable robustness, surviving from the bulk crystal down to the few-layer limit at room temperature. Our theoretical analysis traces its origin to the hybridization between Nb-dz2 orbital chains and the Lieb-like dx2-y2 sublattice, which is further reinforced by Peierls dimerization. Our findings not only establish transition-metal oxychlorides as a robust and tunable platform for flat-band-driven correlated phases under ambient conditions, but also uncover a new orbital-selective design principle for realizing flat bands in quantum materials.

Figures (5)

• Figure 1: Crystal structure, cooperative Peierls distortions and emergence of an orbital-selective flat band in NbOCl$_{2}$. (a-c) Crystal structure of NbOCl$_{2}$ viewed along various crystallographic directions. Nb (blue), O (red) and Cl (yellow) atoms form quasi-2D layers stacked along the a axis. Within each layer, alternating Nb–Nb bond lengths $d_{1}$ and $d_{2}$ highlight pronounced dimerization. (d) Schematic lattices and tight-binding band structures of kagome and Lieb systems, where symmetry-protected flat bands arise from destructive interference. Owing to distinct lattice geometries, the flat band positions differ between these two systems. NbOCl$_{2}$ adopts a related Lieb-like lattice configuration. (e) Orbital energy levels of Nb$^{4+}$ ($4d^{1}$) atoms and bonding configurations in the dimerized state. (f) Density functional theory (DFT) band structure of the monolayer dimerized NbOCl$_{2}$. The rectangular BZ is defined with $k_{x}$ along the $\Gamma$–X direction (parallel to the b axis) and $k_y$ along the $\Gamma$–Y direction (parallel to the c axis). The narrow, well-isolated flat band (FB) near the Fermi level remains robust in the ultrathin limit, with the Fermi level inside the gap between the FB and dispersive bands. (g) Schematic of the cooperative mechanism: dimerization driven by Peierls distortion, Lieb-like lattice interference and orbital localization collectively stabilize the orbital-selective flat band.
• Figure 2: Momentum-independent flat band and its orbital character in NbOCl$_{2}$ characterized by ARPES and DFT. (a) Constant-energy contours of NbOCl$_{2}$ measured at room temperature (T= 300 K) using photon energy h$\nu$ = 66 eV, shown at representative binding energies $E_{B}$ = 2.2, 4.5, 5.5, 6.5 and 7.5 eV. High-symmetry points of the BZ are marked by orange circles. Momentum cuts (Cut 1 and Cut 2) indicated by black lines are chosen for detailed analysis. (b) Band structures along the $\Gamma$–X direction (Cut 1) recorded at h$\nu$ = 100 eV. The middle panel shows the corresponding EDCs integrated at $\Gamma$ and X points. The right panel presents the DFT band structure along the same path for comparison. (c) Same as (b), but along the $\Gamma$–Y direction (Cut 2). (d) Orbital-resolved projected density of states (PDOS) for monolayer NbOCl$_{2}$. Inset shows the real-space charge density of the FB's wavefunction from DFT. (e) Photon-energy-dependent ARPES spectra along the $\Gamma$–X direction, taken at h$\nu$= 40, 70, 100, 120 and 140 eV, revealing the $k_{z}$ dependence of the valence bands.
• Figure 3: Universality and tunability of orbital-selective flat bands across NbOCl$_{2}$ and TaOCl$_{2}$. (a) Experimental band dispersions of NbOCl$_{2}$ (left) and TaOCl$_{2}$ (right) obtained by ARPES at h$\nu$= 70 eV. Momentum cuts measured with LH and LV polarized light were combined, and are shown overlaid with DFT calculations for the corresponding monolayer structures. High-symmetry points of the BZ are defined in Fig. 1d. (b) Photoemission spectra of NbOCl$_{2}$ (red) and TaOCl$_{2}$ (blue), integrated along $\Gamma$–X direction. Core-level peaks corresponding to Nb/Ta, O and Cl elements are labeled. (c) Band dispersions along $\Gamma$–X direction for NbOCl$_{2}$ and TaOCl$_{2}$, measured separately with LH and LV light polarizations at h$\nu$ = 70 eV. (d) Same as in (c) but along the $\Gamma$–Y direction. (e) EDCs from (c), extracted at $\Gamma$ ($k_{x}$ = 0) and X point. Peaks associated with NbOCl$_{2}$ are marked by red triangles, while those of TaOCl$_{2}$ are marked by black triangles. (f) Same as in (e) but for the data in (d).
• Figure 4: Persistence of the flat band in few-layer NbOCl$_{2}$ revealed by ARPES. (a) Schematic of the micro-ARPES setup for NbOCl$_{2}$/graphene heterostructures, with few-layer NbOCl$_{2}$ encapsulated by graphene and supported on hBN and Au contacts on a SiO$_{2}$/Si substrate. (b) Optical image of the device, showing patterned Au electrodes, hBN support, graphene encapsulation and NbOCl$_{2}$ flakes. (c) Atomic force microscopy (AFM) image of the NbOCl$_{2}$ flake with height profile (inset), confirming a thickness of $\sim$2 nm (three layers). The position of the profile is marked by the grey dashed line. (d–f) ARPES dispersions along $\Gamma$–X for graphene capped few-layer NbOCl$_{2}$ (d), bare graphene (e) and bulk NbOCl$_{2}$ (f), measured with photon energy h$\nu$= 80 eV under LH polarization. The flat band (FB) in the few-layer NbOCl$_{2}$ sample is marked by an orange arrow in (d). (g) EDCs at $\Gamma$ point comparing few-layer (purple) and bulk (red) NbOCl$_{2}$, showing that the flat band persists in reduced thickness.
• Figure 5: Microscopic origin of the flat band formation from orbital hybridization and Peierls dimerization. (a) Orbital-projected band structure of monolayer NbOCl$_{2}$ in the non-CDW phase, with the color scale indicating the relative contributions from Nb $d$ orbitals and Cl/O $p$ orbitals. The inset highlights the nearly dispersionless Nb $d_{z^{2}}$ and $d_{yz}$ bands, while retaining finite dispersion in certain momentum regions. (b) Same as in (a) but for the CDW phase. Due to the Peierls distortion, the $d_{z^{2}}$ and $d_{yz}$ -derived bands exhibit enhanced flatness relative to the non-CDW phase and open up a CDW gap, indicating an origin predominantly rooted in orbital and charge character rather than structural effects alone. (c-e) Schematic orbital–lattice building blocks and corresponding tight-binding band structures for representative real-space configurations in the CDW phase. The Lieb-lattice type configuration (c) is primarily contributed by Nb $d_{xy}$, Cl $p_{x}$ and O $p_{y}$ orbitals. The SSH chain (d) originates from Nb $d_{yz}$, Cl $p_{y}$ and $p_{z}$ orbitals. The hybridization configuration (e) arises from strong hybridization between Nb $d_{x^{2}-y^{2}}$ and $d_{z^{2}}$ orbitals that are further coupled to Cl/O $p$ states, enhancing the flatness of the flat band.