Charge fluctuations, phonons and superconductivity in multilayer graphene

Ziyan Li,Xueheng Kuang,Alejandro Jimeno-Pozo,Héctor Sainz-Cruz,Zhen Zhan,Shengjun Yuan,Francisco Guinea

DOI: https://doi.org/10.1103/PhysRevB.108.045404

2024-01-31

Abstract:Motivated by the recent experimental detection of superconductivity in Bernal bilayer (AB) and rhombohedral trilayer (ABC) graphene, we study the emergence of superconductivity in multilayer graphene based on a Kohn-Luttinger (KL)-like mechanism in which the pairing glue is the screened Coulomb interaction. We find that electronic interactions alone can drive superconductivity in AB bilayer graphene and ABC trilayer graphene with the critical temperatures in good agreement with the experimentally observed ones, allowing us to further predict superconductivity from electronic interactions in Bernal ABA trilayer and ABAB tetralayer and rhombohedral ABCA tetralayer graphene. By comparing the critical temperatures ($T_c$) of these five non-twisted graphene stacks, we find that the ABC trilayer graphene possesses the highest $T_c\sim100$ mK. After considering the enhancement of superconductivity due to Ising spin-orbit coupling, we observe that the AB bilayer graphene has the largest enhancement in the critical temperature, increasing from 23 mK to 143 mK. The superconducting behaviors in these non-twisted graphene stacks could be explained by the order parameters (OPs). The OPs of Bernal stacks preserve intravalley $C_3$ symmetry, whereas rhombohedral stacks break it. In all stacks, the OPs have zeroes and change signs between valleys, which means that these multilayers of graphene are nodal spin-triplet superconductors. Moreover, dressing the purely electronic interaction with acoustic phonons, we observe minor changes of the critical temperatures in these five stacks. We adopt the KL-like mechanism to investigate the tendency of superconductivity in multilayer graphene without fitting parameters, which could provide guidance to future experiments exploring superconductivity in non-twisted graphene.

Mesoscale and Nanoscale Physics,Superconductivity

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the superconductivity mechanism and its characteristics in multilayer graphene. Specifically, the research aims to:

1. **Explain experimental phenomena**: In recent years, superconductivity has been detected in Bernal bilayer (AB) and rhombohedral trilayer (ABC) graphene, but the mechanism behind it is still unclear. This paper explains these experimental phenomena through a Kohn - Luttinger (KL) - like mechanism, especially how electron - electron interactions drive superconductivity.

2. **Predict superconductivity in new structures**: Based on the KL mechanism, predict whether other multilayer graphene structures that have not been experimentally verified (such as ABA trilayer, ABAB four - layer and ABCA four - layer) also have superconductivity, and provide theoretical basis.

3. **Explore the enhancing factors of superconductivity**: Study the influence of Ising spin - orbit coupling (SOC) on the superconducting critical temperature (\(T_c\)), especially the significant enhancement phenomenon observed in AB bilayer graphene.

4. **Understand the influence of different stacking methods**: Compare the differences in superconductivity in multilayer graphene with different stacking methods (Bernal vs. rhombohedral), and reveal their underlying physical mechanisms, such as the symmetry and nodal structure of the order parameters (OPs).

### Research background

In recent years, superconductivity in non - twisted multilayer graphene (such as AB bilayer and ABC trilayer) has been discovered experimentally, which provides a new platform for studying correlation effects in two - dimensional materials. However, the superconductivity mechanism has not been fully clarified, especially regarding the issues of pairing symmetry and interaction mechanism.

### Main findings

- **Superconductivity mechanism**: The research shows that electron - electron interactions alone can drive superconductivity in AB bilayer and ABC trilayer graphene, and the calculated critical temperature is consistent with the experimental observations.

- **Predict superconductivity in new structures**: It is predicted that ABA trilayer, ABAB four - layer and ABCA four - layer graphene should also have superconductivity, among which ABC trilayer graphene has the highest critical temperature (about 100 mK).

- **Influence of Ising SOC**: Ising SOC significantly enhances the critical temperature of AB bilayer graphene, from 23 mK to 143 mK, while expanding the energy range where superconductivity occurs.

- **Order parameter analysis**: All the studied multilayer graphene are valley - singlet, spin - triplet superconductors, and the order parameter has sign changes and nodes within each valley, indicating that these materials are nodal spin - triplet superconductors.

### Conclusion

This research not only explains the existing experimental phenomena, but also provides theoretical guidance for future exploration of superconductivity in non - twisted multilayer graphene. The research results show that superconductivity does not simply depend on the number of layers, but is related to more subtle factors such as Fermi surface characteristics. In addition, Ising SOC plays an important role in enhancing superconductivity, especially in AB bilayer and ABCA four - layer graphene.