Tunable Electronic Transport in Pd$_3$O$_2$Cl$_2$ Kagome Bilayers: Interplay of Stacking Configuration and Strain

Kagome lattice bilayers offer unique opportunities for engineering electronic properties through interlayer stacking and strain. We report a comprehensive first-principles study of Pd$_3$O$_2$Cl$_2$ kagome bilayers, examining four stacking configurations (AA, AA$’$, AB, AB$’$). Our calculations reveal dramatic stacking-dependent band gap modulation from 0.08 to 0.76eV, with the AB$’$ configuration being the most thermodynamically stable. All stackings exhibit robust mechanical stability with Young’s moduli of 54.82-61.97N/m and ductile behavior suitable for flexible electronics. Carrier effective masses show significant stacking dependence, ranging from 2.39-6.35~$m_0$ for electrons and 0.67-1.55~$m_0$ for holes. Strain engineering of the AB$’$ bilayer demonstrates non-monotonic band gap tuning and asymmetric modulation of carrier masses, with hole effective masses showing stronger strain sensitivity. These results establish Pd$_3$O$_2$Cl$_2$ bilayers as a promising platform for strain-engineered kagome-based quantum devices, where stacking order and mechanical deformation provide complementary control over electronic transport.

💡 Research Summary

In this work the authors present a systematic first‑principles investigation of Pd₃O₂Cl₂ kagome bilayers, focusing on how interlayer stacking order and biaxial strain jointly influence mechanical stability, elastic response, and electronic transport properties. Four distinct stacking configurations—AA, AA′, AB, and AB′—were constructed by vertically stacking monolayers while preserving the hexagonal symmetry. Geometry optimizations were performed with density‑functional theory (PBE functional) using a 5 × 5 × 1 k‑point mesh, 500 eV plane‑wave cutoff, and stringent convergence criteria (10⁻⁸ eV energy, 0.02 eV Å⁻¹ forces). The relaxed lattice constants range from 5.67 Å to 5.73 Å, and interlayer separations vary from 3.51 Å (AA) down to 2.05 Å (AB′). Formation‑energy analysis shows that AA′ (−0.088 eV/atom) and AB′ (−0.104 eV/atom) are thermodynamically favorable, with AB′ being the most stable configuration.

Mechanical stability was assessed through 2D elastic constants computed with the ElasTool. All four bilayers satisfy the Born criteria: eigenvalues of the elastic tensor are positive (λ₁≈18–21 N/m, λ₂≈35–42 N/m, λ₃≈115–122 N/m). Young’s moduli lie between 54.8 and 61.9 N/m, shear moduli between 18.7 and 20.8 N/m, and Poisson ratios cluster around 0.48–0.52, indicating conventional elastic behavior. The G/B ratio remains below the ductility threshold of 0.57, classifying every stacking as ductile and thus suitable for flexible device applications. Acoustic phonon analysis yields longitudinal sound velocities of 3.9–4.1 km s⁻¹ and transverse velocities of 1.9–2.1 km s⁻¹, with Debye temperatures spanning 461–493 K, confirming robust dynamical stability.

Electronic structure calculations reveal that the characteristic kagome features—flat bands and Dirac‑like crossings at the K point—persist across all stackings. However, the band gap is highly sensitive to stacking registry. AA and AB exhibit narrow direct gaps of 0.08 eV and 0.17 eV, respectively, due to strong Pd‑4d/O‑2p interlayer hybridization that compresses the valence‑conduction separation. In contrast, AA′ and AB′ display much larger gaps (0.71–0.76 eV) because the relative lateral shift and layer flipping reduce Pd‑Pd overlap, enhancing bonding–antibonding splitting. Projected density of states shows valence‑band maxima dominated by Pd‑4d/O‑2p hybridization, while conduction‑band minima involve Pd‑4d with notable Cl‑p contributions. Effective masses are markedly stacking‑dependent: electron masses range from 2.39 to 6.35 m₀, hole masses from 0.67 to 1.55 m₀, with holes generally lighter but more variable.

Strain engineering was performed on the most stable AB′ bilayer by applying biaxial strains from –6 % (compression) to +6 % (tension). The band gap exhibits a non‑monotonic response, decreasing to ≈0.65 eV under compression and increasing to ≈0.78 eV under tension. Electron effective masses change modestly with strain, whereas hole masses display a pronounced asymmetric behavior: they shrink under compression and swell dramatically under tension. This asymmetry originates from the stronger dependence of the hole‑band dispersion on the Pd‑O‑Pd hopping pathways, which are directly altered by lattice deformation.

Overall, the study demonstrates that stacking order and mechanical strain provide complementary, highly tunable knobs for controlling the electronic landscape of Pd₃O₂Cl₂ kagome bilayers. The ability to modulate band gaps over a factor of ten, to engineer carrier effective masses, and to retain ductile mechanical response makes these bilayers promising candidates for flexible electronics, strain‑controlled optoelectronic devices, and platforms where strong electronic correlations or topological phases could be induced by external perturbations. The authors suggest that experimental synthesis of the AB′ stacking, followed by strain‑application techniques (e.g., substrate engineering or flexible substrates), could unlock a new class of kagome‑based quantum materials with designer transport properties.