Edge-Dependent Superconductivity in Twisted Bismuth Bilayers

Abstract: Twisted bilayers offer a compelling and, at times, confounding platform for the engineering of new twistronic materials. Whereas standard studies almost exclusively focus on the explicit enigma that is presented by twist-angles, perhaps better epitomized by the related phenomena that have been observed in twisted bilayer graphene, functional devices necessarily face a fundamental concern: boundary heterogeneity in their structures. In this study, we address this concern by strictly investigating the electronic properties of twisted bismuth bilayers at the flake's edges and the vibrational properties of the flake. Twisted flakes exhibit continuous variations of these properties, away from the bulk, as we herein report using ab initio density functional theory, by systematically mapping the drastic evolution of band topology, electronic density of states, and possible superconductivity. Our work reveals a dramatic, non-fortuitous consequence of the structural disorder at the edges of the flakes: an enhanced electronic density of states at the Fermi level. This enhancement reaches a maximum of 10 times that of perfect-crystalline bismuth. Given that the superconducting critical temperature, Tc, is exponentially dependent on the electronic density of states at the Fermi level, this substantial structural variation immediately suggests a powerful mechanism for vastly increasing Tc. We also identify the twist-angle as a new critical parameter in designing novel engineering devices with topologically enhanced properties. Our results provide a necessary theoretical framework for interpreting new data for the upcoming generation of twistronic heterogeneous materials, and pave the way to search for atomic disordered metastable structures that could lead to enhanced superconducting transition temperatures.

• The paper reveals that edge disorder and twist-angle heterogeneity dramatically enhance superconductivity in twisted bismuth bilayers, with a 7000-fold increase in Tc at specific angles.
• It employs large-scale DFT calculations and detailed atomistic modeling to delineate a core-edge dichotomy, featuring a crystalline core and a disordered edge.
• Findings suggest that deliberate edge engineering can optimize superconducting devices in 2D materials beyond conventional infinite-layer models.

Edge-Dependent Superconductivity in Twisted Bismuth Bilayers: An Expert Analysis

Introduction and Motivation

The study of twistronic materials—stacked monolayers with precisely tuned twist-angles—has shifted paradigms in quantum materials research, particularly after the discovery of magic-angle superconductivity in twisted bilayer graphene. However, most prior investigations have focused on pristine, infinite systems, largely neglecting the roles of physical boundaries and structural disorder that are intrinsic to real materials and devices. "Edge-Dependent Superconductivity in Twisted Bismuth Bilayers" (2511.21657) systematically addresses this gap by focusing on how edge disorder and twist-angle heterogeneity in finite bismuth bilayer flakes fundamentally reshape electronic and superconducting properties, with direct design implications for future 2D quantum devices.

Structural Dichotomy: Crystalline Core versus Disordered Edge

Using large-scale DFT calculations with explicit atomistic modeling of finite flakes, the authors reveal a distinct core-edge dichotomy in twisted bismuth bilayers (TBBs). The optimized geometries show the flake's core remains a compressed, crystalline-like phase closely resembling compressed Bi-I, while the edge region exhibits significant amorphization and structural disorder. Quantitative analysis through pair distribution functions (PDFs) and plane-angle distributions robustly confirms this dichotomy. The edge is not merely a zone of geometric defect, but constitutes a unique, highly-strained, non-crystalline phase with no evidence for a liquid-like character.

This boundary disorder arises unavoidably from lattice relaxations and local strain gradients near the edges, which are exacerbated by the finite flake size and rotation-induced moiré distortions—central real-world effects omitted by periodic, infinite-layer models. The transition region separating core and edge further highlights a gradual spatial modulation of structural motifs, including triangles, squares, and pentagonal rings.

Electronic Structure: Edge-Driven Enhancement of Density of States

DFT-based electronic structure calculations demonstrate a strong spatial dependence of the density of states at the Fermi level, $N(E_F)$. The core’s $N(E_F)$ remains semi-metallic and nearly invariant across the entire studied twist-angle range (0°–30°), with peak values of 0.59 states·eV⁻¹·atom⁻¹ at 15°. In contrast, the edge region exhibits a highly non-monotonic, oscillatory $N(E_F)$ with sharp local maxima at twist angles of 3.5°, 10.0°, 20.5°, and 30.0°. At these peaks, the edge $N(E_F)$ can be up to 10 times that of perfect-crystalline Bi-I (∼0.15 states·eV⁻¹·atom⁻¹). This enhancement is attributed to the localization of electronic states, reminiscent of van Hove singularities, due to disorder-induced boundary states.

Importantly, the edge $N(E_F)$ sensitivity to twist-angle is disconnected from the core’s properties, establishing the edges as the key locus for emergent phenomena in finite TBB systems, potentially dominating device performance in realistic flakes.

Superconductivity: Disorder- and Angle-Induced Critical Temperature Enhancement

The superconducting transition temperature $T_c$ is estimated using the Mata-Valladares approach, combining calculated $N(E_F)$ with angle-dependent Debye temperatures, $\Theta_D$, derived from ab initio vibrational density of states. The Debye temperature shows periodic oscillations with twist-angle, but its effect on $T_c$ is secondary to the dramatic variations in $N(E_F)$.

A critical, quantifiable result is the order-of-magnitude difference in $T_c$ between core and edge regions:

• Core: $T_c$ remains bounded, peaking at 4.67 K (15° twist), consistent with phonon-mediated BCS behavior in crystalline bismuth phases.
• Edge: $T_c$ is highly oscillatory, with four prominent maxima that directly track $N(E_F)$, culminating in a maximum $T_c$ of 37.53 K at 30° twist—representing more than a 7000-fold enhancement compared to the superconducting transition in crystalline Bi-I (0.53 mK).

This result underscores that maximal $T_c$ is achieved not at "magic" twist-angles classically associated with flat bands in graphene, but instead arises from a synergy between edge-induced disorder, van Hove singularities, and moiré patterning. The findings reinforce theoretical and experimental claims from other 2D systems that controlled disorder can enhance superconductivity [50-54].

Implications, Limitations, and Future Directions

This work establishes a theoretical framework for engineering high-$T_c$ superconductivity in finite, twisted bilayer bismuth structures by deliberately exploiting edge disorder and twist-angle tuning. The identification of the edge as the site of profound $T_c$ enhancement has significant consequences:

• Device Engineering: Edge-angle engineering becomes a critical design variable for next-generation 2D quantum devices.
• Twistronics Beyond Graphene: These results generalize the twistronics paradigm to materials with intrinsically strong spin-orbit coupling and topological features, such as bismuthene.
• Disorder as a Tool: Rather than treating disorder as a purely detrimental effect, controlled structural disorder at boundaries emerges as a pathway to realize new superconducting phases—directly challenging traditional approaches based solely on bulk or periodic models.

Prospective future work should focus on:

• Experimental Validation: In situ spectroscopic and transport experiments in edge-engineered TBB flakes to observe the predicted $T_c$ oscillations and electronic signatures.
• Beyond Bismuth: Application of this framework to other van der Waals materials and topological systems with strong spin-orbit coupling.
• Multiscale Modeling: Integration of atomistic disorder with mesoscale models and many-body theory to quantify effects beyond mean-field DFT.

Conclusion

"Edge-Dependent Superconductivity in Twisted Bismuth Bilayers" (2511.21657) demonstrates that atomic edge disorder, coupled with twist-angle engineering, fundamentally controls superconducting properties in finite bismuth bilayer flakes. The results challenge the prevailing focus on homogeneous, infinite systems, and provide a rigorous theoretical underpinning for experimentally leveraging heterogeneity and edge effects in future 2D superconducting devices. The paradigm of edge-angle design opens a robust new avenue for probing and optimizing emergent quantum phenomena in twistronic materials.