Orbital-Selective Mott and Antiferromagnetic Phases in Diagonally Compressed Kagome Lattice

Abstract: We perform determinant quantum Monte Carlo simulations of the half-filled Hubbard model on a diagonally compressed kagome lattice, introducing exponential decay long-range hopping $t(r) = t_0 \exp\bigl(-r / r_0\bigr)$ to account for the evolving bond length. By varying the lattice angle $θ$ and the on-site interaction $U$, double occupancy, charge compressibility, and spin-spin correlation functions of the whole system and each sub-lattice are measured. We find that geometric compression induces a clear sublattice differentiation: for $θ\gtrsim52<sup>\circ$, the A sublattice establishes long-range hoppings, which in turn suppresses the metallic behavior of the $B/C$ sublattice and drives a selective Mott transition; for $θ\lesssim52<sup>\circ$, the $B$-$C$ chains develop long-range antiferromagnetic correlations within the finite-size simulations, which in turn suppresses the metallic behavior of the $A$ sublattice and drives a selective Mott transition. The critical interaction $U<sup>c_A$ for the $A$ sites decreases sharply near the onset of $B$-$C$ antiferromagnetic correlations, while $U<sup>c_{B/C}$ increases. These competing orders give rise to an orbital-selective Mott phase and a rich $U$-$θ$ phase diagram featuring paramagnetic-metal, paramagnetic-Mott, antiferromagnetic-metal, and antiferromagnetic-Mott states. Our results highlight the complex interplay between lattice geometry, magnetic frustration, and strong correlations in frustrated two-dimensional systems.

• The paper demonstrates a two-stage orbital-selective Mott transition driven by lattice compression and sublattice differentiation.
• It employs determinant quantum Monte Carlo to uncover the interplay of geometric frustration, electronic compressibility, and antiferromagnetic correlations.
• Findings reveal anisotropic charge responses and a comprehensive phase diagram, offering insights for engineered quantum simulators.

Orbital-Selective Mott and Antiferromagnetic Phases in the Diagonally Compressed Kagome Lattice

Introduction

This paper, "Orbital-Selective Mott and Antiferromagnetic Phases in Diagonally Compressed Kagome Lattice" (2607.04621), uses determinant quantum Monte Carlo (DQMC) to study the half-filled Hubbard model on the kagome lattice under diagonal compression. By introducing an exponential decay in long-range hopping as a function of bond length, the model accurately captures the effect of geometric deformations, departing from previous works that rely primarily on bond-anisotropic hopping or artificially constructed interpolations between Lieb and kagome lattices.

The study systematically explores the interplay among geometric frustration, strong correlations, and magnetic ordering. The continuous compression parameter, represented by the angle $\theta$, controls both the degree of frustration and the relative strengths of different hopping channels. The consequences of this subtle network restructuring are evaluated in terms of sublattice-resolved double occupancy, electronic compressibility, and real/momentum-space spin correlation functions, providing a microscopic understanding of orbital-selective Mott transitions and the emergence of antiferromagnetic order along specific lattice directions.

Figure 1: (a) Schematic of the kagome lattice before and after compression; (b) exponential decay of hopping strength with distance; (c) change in geometry with $\theta$; (d, e) sublattice-resolved hopping profiles.

Model, Methodology, and Numerical Formulation

The model consists of the Hubbard Hamiltonian on the kagome lattice with onsite repulsion $U$ and exponentially decaying hopping $t(r) = t_0 e^{-r/r_0}$. The key geometric parameter, $\theta$, denotes the angle between the two primitive vectors. Reducing $\theta$ from $60^\circ$ compresses the lattice diagonally, shortening the $BC$ bonds while the $AB$ and $AC$ bonds are held fixed. As a result, the model interpolates between an isotropic kagome lattice and a regime where $\theta$0-chain character dominates, generating strong sublattice differentiation and an anisotropic network of kinetic and interaction terms.

DQMC provides finite-temperature access to thermodynamic, charge, and spin observables. The simulations are performed at half-filling by numerically imposing the chemical potential to satisfy $\theta$1, compensating for the loss of strict particle-hole symmetry due to complex hopping patterns. System sizes up to $\theta$2 are explored, and a multi-parameter sweep over $\theta$3 and $\theta$4 is conducted for sublattice-resolved analysis.

Figure 2: Average sign $\theta$5 at $\theta$6 as a function of compression angle $\theta$7 for different $\theta$8, showing that compression generally alleviates the sign problem.

Electronic Structure Evolution and Sublattice Differentiation

The non-interacting electronic spectrum exhibits a remarkable transformation under compression. At $\theta$9, the band structure features the characteristic kagome flat band and Dirac point. As $U$0 decreases and $U$1 hopping is enhanced, the flat band acquires significant dispersion and the bandwidth increases, leading to strong differentiation between sublattice $U$2 (with increased long-range hopping) and sublattices $U$3 (experiencing dominant short-range kinetic processes).

Figure 3: Band structures for $U$4, and $U$5; decreasing $U$6 drives the flat band dispersive and enhances bandwidth/differentiation.

This geometry-induced inequivalence is already visible at the level of sublattice occupation numbers, and becomes more pronounced at low temperatures and intermediate interactions, with $U$7 maximized for small $U$8 and large compression.

Figure 4: Temperature dependence of the chemical potential $U$9 required to achieve half-filling, showing strong $t(r) = t_0 e^{-r/r_0}$0 and $t(r) = t_0 e^{-r/r_0}$1 dependencies.

Metal-Insulator and Orbital-Selective Mott Physics

Double occupancy and compressibility are analyzed both globally and per sublattice. The temperature dependence of $t(r) = t_0 e^{-r/r_0}$2 and $t(r) = t_0 e^{-r/r_0}$3 (and their suppression with $t(r) = t_0 e^{-r/r_0}$4) are consistent with Mott physics, but display marked sublattice sensitivity due to the geometry-driven bandwidth and effective-coupling rearrangements.

Figure 5: Average particle occupation per sublattice upon compression, illustrating anisotropic charge response.

Figure 6: Double occupancy $t(r) = t_0 e^{-r/r_0}$5 vs temperature for (a-c) whole lattice, (d-f) $t(r) = t_0 e^{-r/r_0}$6 sublattice, (g-i) $t(r) = t_0 e^{-r/r_0}$7 sublattices, at representative $t(r) = t_0 e^{-r/r_0}$8 and $t(r) = t_0 e^{-r/r_0}$9.

Electronic compressibility $\theta$0 is used to extract an effective charge activation gap via low-temperature fits ($\theta$1). The sign change of $\theta$2 is used to identify the Mott transition point. Notably, the critical interaction strength $\theta$3 differs significantly between sublattices as a function of $\theta$4, yielding a two-stage, orbital-selective Mott transition: for weak compression, $\theta$5; for strong compression, the ordering reverses, $\theta$6.

Figure 7: Temperature dependence of $\theta$7 for whole lattice, $\theta$8, and $\theta$9 sublattices.

Figure 8: Extracted charge activation scale $\theta$0 for the whole lattice, $\theta$1, and $\theta$2 sublattices as a function of $\theta$3 and $\theta$4, showing a dramatic crossover in $\theta$5.

Magnetism: Real- and Momentum-Space Signatures

The study computes both longitudinal ($\theta$6) and transverse ($\theta$7) spin correlation functions. $\theta$8 is preferred due to superior statistical properties in DQMC. Along the $\theta$9 path, correlations are always short-ranged, indicating persistent paramagnetic character. However, for strong compression and large $60^\circ$0, the $60^\circ$1 path develops long-ranged, alternating (antiferromagnetic) spin correlations, as seen both in real space and in enhanced $60^\circ$2 peaks in the momentum-resolved spin structure factor.

Figure 9: (a–c) $60^\circ$3 along $60^\circ$4; (d–f) along $60^\circ$5; (g) difference in sublattice (B vs C) local moments vs $60^\circ$6.

Figure 10: Logarithmic decay of $60^\circ$7 for $60^\circ$8 and $60^\circ$9—only the $BC$0 path in the strongly-compressed, high-$BC$1 regime develops persistent long-range correlations.

Figure 11: (a) $BC$2 at $BC$3 for various $BC$4; (b) correlation ratio $BC$5 vs $BC$6 for different $BC$7, highlighting the AFM transition.

Comprehensive Phase Diagram and Physical Interpretation

The extracted global phase diagram shows four regions: paramagnetic metal, paramagnetic Mott insulator, antiferromagnetic metal, and antiferromagnetic Mott insulator. The metal-insulator boundary is set by the charge gap, whereas the magnetic boundary is operationally determined by the onset of statistically significant long-distance AFM correlations on the $BC$8 path. The boundaries do not generally coincide, leading to intermediate metallic regimes with substantial AFM correlations.

Figure 12: Phase diagram in the ($BC$9) plane, with metal-Mott and AFM boundaries distinguished, revealing orbital-selective transitions.

The phase structure results from the competition between kinetic delocalization (bandwidth-controlled), interaction-induced localization, and frustration-driven suppression of long-range magnetic order. For strong frustration (large $AB$0), the system remains paramagnetic even deep into the Mott regime, while for reduced $AB$1, the suppression of $AB$2 hopping pathways by AFM order on the $AB$3 chains generates a rapid drop in $AB$4 and cooperatively generates orbital-selective and magnetically ordered phases.

Implications and Future Directions

The study demonstrates the potential for geometry-driven engineering of orbital-selective Mott phenomena and anisotropic magnetism in frustrated 2D lattices. The mechanism—based on the manipulation of hopping strengths via lattice compression—directly translates to experimental platforms such as moiré superlattices and cold-atom optical lattices, where tuning of tunneling amplitudes is feasible. The two-stage nature of the Mott transition and the emergence of AFM correlations localized to certain lattice directions provide a microscopic blueprint for emergent orbital/magnetically selective behaviors in materials with tunable frustration.

Further work is needed to address thermodynamic-limit extrapolations, the possibility of true phase transitions as opposed to sharp crossovers, and to explore doping and finite-temperature phase structures. Extensions to multi-orbital real-material models and time-dependent control offer routes to designer quantum simulators and new strongly correlated phases.

Conclusion

This paper provides a comprehensive DQMC study of the half-filled Hubbard model on a diagonally compressed kagome lattice with distance-dependent hopping, revealing rich orbital-selective Mott physics and frustration-driven AFM ordering. The key results include the identification of a two-stage sublattice-resolved Mott transition, the quantification of the AFM boundary as a function of geometry and $AB$5, and a complete ($AB$6, $AB$7) phase diagram capturing the interplay of strong correlations and geometric frustration. The findings have direct relevance for synthetic quantum-matter platforms and lay the groundwork for theoretical and experimental advances in frustrated correlated systems.