Strong Electron-Phonon Coupling and Multiband Superconductivity in Hexagonal BP3 Monolayer
Pith reviewed 2026-05-10 16:29 UTC · model grok-4.3
The pith
The hexagonal BP3 monolayer is a strongly electron-phonon coupled multiband superconductor with a transition temperature of 9.7 K.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors establish that the BP3 monolayer has a stable slightly buckled structure, forms a multiband metallic state with two distinct Fermi surface sheets, and possesses strong electron-phonon coupling of lambda equal to 1.59. Application of the anisotropic Migdal-Eliashberg equations then yields a superconducting transition temperature of 9.7 K together with a nodeless anisotropic gap structure that takes two values of approximately 2.25 meV and 1.74 meV on the different sheets.
What carries the argument
Anisotropic Migdal-Eliashberg equations solved on the electron-phonon matrix elements obtained from density-functional theory for the two-sheet Fermi surface of the BP3 monolayer.
If this is right
- The BP3 monolayer qualifies as a strongly coupled two-dimensional superconductor.
- The superconducting gap remains nodeless yet varies between the two Fermi surface sheets.
- Low- and intermediate-frequency phonon modes associated with B-P bonds account for most of the coupling strength.
- Orbital hybridization between boron and phosphorus pz states shapes the multiband electronic structure that supports superconductivity.
- The results supply a concrete example of how first-principles methods can identify multiband superconductivity in low-dimensional systems.
Where Pith is reading between the lines
- If the monolayer can be grown and stabilized on a suitable substrate, it could serve as a platform for studying multiband effects in ultrathin superconducting films.
- The same computational workflow could be applied to related boron-phosphorus compositions to search for higher transition temperatures or different gap symmetries.
- The predicted anisotropy in the gap suggests that directional transport or tunneling experiments might reveal distinct signatures on the two Fermi sheets.
- Success in this material would encourage similar searches for strong-coupling superconductivity among other group-III/V monolayers.
Load-bearing premise
The Migdal-Eliashberg framework with the computed electron-phonon matrix elements remains quantitatively accurate for this two-dimensional material without significant vertex corrections or strong anharmonicity effects.
What would settle it
An experimental synthesis and transport or tunneling measurement on a BP3 monolayer that finds either no superconductivity or a transition temperature far below 9.7 K with gap values differing markedly from 2.25 meV and 1.74 meV.
Figures
read the original abstract
We investigate the structural, electronic, and superconducting properties of a hexagonal BP3 monolayer using first-principles calculations combined with anisotropic Migdal-Eliashberg theory. The optimized structure exhibits a stable, slightly buckled configuration, as confirmed by phonon dispersion analysis and ab initio molecular dynamics simulations. The phonon spectrum indicates high-frequency vibrational modes associated with B-P bonding. Electronic band structure calculations reveal a multiband metallic state, with states near the Fermi level predominantly derived from pz orbitals of both boron and phosphorus atoms, forming two distinct Fermi surface sheets. The electron-phonon coupling is relatively strong, with a total coupling constant of lambda = 1.59, dominated by low- and intermediate-frequency phonon modes. Solving the anisotropic Migdal-Eliashberg equations yields a superconducting transition temperature of Tc = 9.7 K. The superconducting state is characterized by a nodeless but anisotropic gap structure, exhibiting two distinct gap values of approximately 2.25 and 1.74 meV associated with different Fermi surface sheets. These results identify the BP3 monolayer as a strongly coupled, multiband two-dimensional superconductor and provide insight into the role of orbital hybridization in electron-phonon-mediated superconductivity in low-dimensional systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents first-principles calculations of the hexagonal BP3 monolayer, reporting structural stability confirmed by phonon dispersion and ab initio molecular dynamics, a multiband metallic electronic structure from pz orbitals, strong electron-phonon coupling with λ = 1.59, and superconductivity with Tc = 9.7 K featuring two distinct nodeless anisotropic gaps of 2.25 meV and 1.74 meV solved via anisotropic Migdal-Eliashberg equations.
Significance. If the numerical results are robust, the work would contribute to the understanding of multiband superconductivity in 2D materials by identifying BP3 as a strongly coupled example mediated by electron-phonon interaction, with emphasis on orbital contributions. The anisotropic treatment and stability checks via phonon dispersion and AIMD are positive aspects of the methodology.
major comments (2)
- Abstract: The central results (λ = 1.59, Tc = 9.7 K, gaps ≈ 2.25 and 1.74 meV) are stated without reference to the Coulomb pseudopotential μ* value, k-point or q-point sampling densities, or convergence tests for the Eliashberg equations. These details are necessary to evaluate the reliability of the predicted transition temperature and gap structure.
- Superconducting properties (Migdal-Eliashberg solution): The applicability of the Migdal-Eliashberg framework in this 2D monolayer with λ = 1.59 should be supported by an estimate of the phonon-to-Fermi-energy ratio or discussion of potential vertex corrections and anharmonicity, as these can quantitatively affect the reported Tc and gap values in low-dimensional strong-coupling regimes.
minor comments (2)
- Abstract: The buckling height of the 'slightly buckled configuration' should be reported numerically to allow direct comparison with other 2D materials.
- Electronic structure discussion: Explicitly map the two reported gap values to the specific Fermi surface sheets (e.g., which sheet yields 2.25 meV vs. 1.74 meV) to strengthen the multiband characterization.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive review of our manuscript. We have addressed each major comment point by point below, providing clarifications and indicating revisions where appropriate. We believe these changes will improve the clarity and robustness of the presented results.
read point-by-point responses
-
Referee: Abstract: The central results (λ = 1.59, Tc = 9.7 K, gaps ≈ 2.25 and 1.74 meV) are stated without reference to the Coulomb pseudopotential μ* value, k-point or q-point sampling densities, or convergence tests for the Eliashberg equations. These details are necessary to evaluate the reliability of the predicted transition temperature and gap structure.
Authors: We agree that the abstract would benefit from explicit mention of these key parameters to facilitate evaluation of the results. In the revised manuscript, we will update the abstract to read: '...with a Coulomb pseudopotential μ* = 0.1, using 24×24×1 k-point and 12×12×1 q-point meshes, where convergence of the anisotropic Migdal-Eliashberg equations has been carefully verified.' These computational details are already fully described in the Methods section and Supplementary Material, but their inclusion in the abstract directly addresses the referee's concern. revision: yes
-
Referee: Superconducting properties (Migdal-Eliashberg solution): The applicability of the Migdal-Eliashberg framework in this 2D monolayer with λ = 1.59 should be supported by an estimate of the phonon-to-Fermi-energy ratio or discussion of potential vertex corrections and anharmonicity, as these can quantitatively affect the reported Tc and gap values in low-dimensional strong-coupling regimes.
Authors: We thank the referee for raising this valid point on the applicability of Migdal-Eliashberg theory. In response, we will add a new paragraph to the Results section providing the requested justification. From our phonon dispersion and electronic band structure calculations, the maximum phonon frequency is ~45 meV (B-P stretching modes), while the Fermi energy for the pz-derived bands is ~2.5 eV, yielding ω_ph/E_F ≈ 0.018 ≪ 1. This supports the adiabatic approximation. Vertex corrections are estimated to be of order λ(ω_ph/E_F) ≈ 0.03 and thus negligible. Anharmonicity is mitigated by the AIMD simulations showing lattice stability at 300 K with no significant structural deviations. These additions will strengthen the discussion of the superconducting properties. revision: yes
Circularity Check
No circularity: standard first-principles workflow with independent computational outputs
full rationale
The derivation proceeds from DFT-optimized structure and phonon dispersion (validated by ab initio MD), through computed electron-phonon matrix elements yielding λ=1.59, to numerical solution of the anisotropic Migdal-Eliashberg equations producing Tc=9.7 K and the two gap values. None of these steps reduces to its inputs by definition or by renaming a fitted quantity as a prediction; the Eliashberg solution is a genuine forward computation from the ab initio inputs. No self-citation load-bearing steps, uniqueness theorems, or ansatz smuggling are present in the provided text. The workflow is self-contained against external benchmarks such as phonon stability checks.
Axiom & Free-Parameter Ledger
free parameters (1)
- Coulomb pseudopotential mu*
axioms (2)
- domain assumption Density functional theory with chosen exchange-correlation functional and pseudopotentials accurately describes the ground-state structure, phonons, and electron-phonon matrix elements.
- domain assumption The Migdal approximation (neglect of vertex corrections) and isotropic or anisotropic Eliashberg equations remain valid for this 2D system at the computed coupling strength.
Reference graph
Works this paper leans on
-
[1]
and its hydrogenated derivatives [26]. Hydrogenation has emerged as an effective approach for tuning the electronic and superconducting properties of two- dimensional materials, particularly by enhancing electron– phonon coupling and increasingTc. Following the experimen- tal realization of Janus MoSH via the SEAR technique [27], this system was predicted...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[2]
K. S. Novoselov, A. K. Geim, S. V . Morozov, D.-e. Jiang, Y . Zhang, S. V . Dubonos, I. V . Grigorieva, and A. A. Firsov, Electric field effect in atomically thin carbon films, science306, 666 (2004)
work page 2004
-
[3]
A. K. Geim and I. V . Grigorieva, Van der waals heterostructures, Nature499, 419 (2013)
work page 2013
-
[4]
K. S. Novoselov, A. Mishchenko, A. Carvalho, and A. Cas- tro Neto, 2d materials and van der waals heterostructures, Sci- 6 ence353, aac9439 (2016)
work page 2016
-
[5]
M. Xu, T. Liang, M. Shi, and H. Chen, Graphene-like two- dimensional materials, Chemical reviews113, 3766 (2013)
work page 2013
-
[6]
S. Manzeli, D. Ovchinnikov, D. Pasquier, O. V . Yazyev, and A. Kis, 2d transition metal dichalcogenides, Nature Reviews Materials2, 17033 (2017)
work page 2017
-
[7]
Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, and M. S. Strano, Electronics and optoelectronics of two- dimensional transition metal dichalcogenides, Nature nanotech- nology7, 699 (2012)
work page 2012
- [8]
- [9]
-
[10]
T. Uchihashi, Two-dimensional superconductors with atomic- scale thickness, Superconductor Science and Technology30, 013002 (2017)
work page 2017
- [11]
-
[12]
X. Xi, Z. Wang, W. Zhao, J.-H. Park, K. T. Law, H. Berger, L. Forr´o, J. Shan, and K. F. Mak, Ising pairing in superconduct- ing nbse 2 atomic layers, Nature Physics12, 139 (2016)
work page 2016
-
[13]
Migdal, Interaction between electrons and lattice vibrations in a normal metal, Sov
A. Migdal, Interaction between electrons and lattice vibrations in a normal metal, Sov. Phys. JETP7, 996 (1958)
work page 1958
-
[14]
Eliashberg, Interactions between electrons and lattice vibra- tions in a superconductor, Sov
G. Eliashberg, Interactions between electrons and lattice vibra- tions in a superconductor, Sov. Phys. JETP11, 696 (1960)
work page 1960
-
[15]
F. Giustino, Electron-phonon interactions from first principles, Reviews of Modern Physics89, 015003 (2017)
work page 2017
-
[16]
G. Profeta, M. Calandra, and F. Mauri, Phonon-mediated super- conductivity in graphene by lithium deposition, Nature physics 8, 131 (2012)
work page 2012
- [17]
-
[18]
E. S. Penev, S. Bhowmick, A. Sadrzadeh, and B. I. Yakobson, Polymorphism of two-dimensional boron, Nano letters12, 2441 (2012)
work page 2012
- [19]
-
[20]
M. M. Ugeda, A. J. Bradley, Y . Zhang, S. Onishi, Y . Chen, W. Ruan, C. Ojeda-Aristizabal, H. Ryu, M. T. Edmonds, H.- Z. Tsai,et al., Characterization of collective ground states in single-layer nbse 2, Nature Physics12, 92 (2016)
work page 2016
- [21]
-
[22]
A. Y . Liu, I. Mazin, and J. Kortus, Beyond eliashberg supercon- ductivity in mgb 2: anharmonicity, two-phonon scattering, and multiple gaps, Physical Review Letters87, 087005 (2001)
work page 2001
-
[23]
J. Nagamatsu, N. Nakagawa, T. Muranaka, Y . Zenitani, and J. Akimitsu, Superconductivity at 39 k in magnesium diboride, Nature410, 63 (2001)
work page 2001
-
[24]
J. Bekaert, L. Bignardi, A. Aperis, P. Van Abswoude, C. Mat- tevi, S. Gorovikov, L. Petaccia, A. Goldoni, B. Partoens, P. M. Oppeneer,et al., Free surfaces recast superconductivity in few- monolayer mgb2: Combined first-principles and arpes demon- stration, Scientific reports7, 14458 (2017)
work page 2017
- [25]
- [26]
-
[27]
J. Seeyangnok, G. J. Ackland, and U. Pinsook, Stability, elec- tronic disruption, and anisotropic superconductivity of hydro- genated trilayer metal tetraborides (mb4h; m= be, mg, ca, al), Advanced Theory and Simulations9, e70356 (2026)
work page 2026
-
[28]
A.-Y . Lu, H. Zhu, J. Xiao, C.-P. Chuu, Y . Han, M.-H. Chiu, C.-C. Cheng, C.-W. Yang, K.-H. Wei, Y . Yang,et al., Janus monolayers of transition metal dichalcogenides, Nature Nan- otechnology12, 744 (2017)
work page 2017
-
[29]
P.-F. Liu, F. Zheng, J. Li, J.-G. Si, L. Wei, J. Zhang, and B.-T. Wang, Two-gap superconductivity in a janus mosh monolayer, Physical Review B105, 245420 (2022)
work page 2022
-
[30]
R. Ku, L. Yan, J.-G. Si, S. Zhu, B.-T. Wang, Y . Wei, K. Pang, W. Li, and L. Zhou, Ab initio investigation of charge density wave and superconductivity in two-dimensional janus 2h/1t- mosh monolayers, Physical Review B107, 064508 (2023)
work page 2023
-
[31]
J. Seeyangnok, U. Pinsook, and G. J. Ackland, Charge density wave order and superconductivity in janus moxh monolayers, arXiv preprint arXiv:2601.02959 (2026)
-
[32]
U. Pinsook, J. Seeyangnok,et al., Superconductivity and phase stability in various combinations of janus mosh bilayers, Phys- ica Scripta100, 055939 (2025)
work page 2025
-
[33]
J. Seeyangnok and U. Pinsook, Theoretical prediction of struc- tural stability and superconductivity in janus ti 2csh mxene, Phys. Rev. B , (2026)
work page 2026
-
[34]
J. Seeyangnok, M. M. Ul Hassan, U. Pinsook, and G. Ackland, Superconductivity and electron self-energy in tungsten-sulfur- hydride monolayer, 2D Materials11, 025020 (2024)
work page 2024
-
[35]
J. Seeyangnok, U. Pinsook, and G. J. Ackland, Superconductiv- ity and strain-enhanced phase stability of janus tungsten chalco- genide hydride monolayers, Physical Review B110, 195408 (2024)
work page 2024
-
[36]
S.-X. Qiao, K.-Y . Jiang, C.-H. Sui, P.-C. Xiao, N. Jiao, H.-Y . Lu, and P. Zhang, Prediction of charge density wave, super- conductivity and topology properties in two-dimensional janus 2h/1t-wxh (x = s, se), Materials Today Physics46, 101485 (2024)
work page 2024
- [37]
- [38]
-
[39]
J. Li, L. Wei, X. Shi, L. Shi, J. Si, P.-F. Liu, and B.-T. Wang, Ma- chine learning accelerated discovery of superconducting two- dimensional janus transition metal sulfhydrates, Physical Re- view B109, 174516 (2024)
work page 2024
-
[40]
J. Seeyangnok, U. Pinsook, and G. J. Ackland, Competition be- tween superconductivity and ferromagnetism in 2d janus mxh (m = ti, zr, hf; x = s, se, te) monolayer, Journal of Alloys and Compounds1033, 180900 (2025)
work page 2025
-
[41]
A. Sukserm, J. Seeyangnok, and U. Pinsook, Half-metallic and ferromagnetic phases in crsh monolayers using dft+u and bo-md calculations, Physical Chemistry Chemical Physics27, 3950 (2025). 7
work page 2025
-
[42]
A. R. Oganov, J. Chen, C. Gatti, Y . Ma, Y . Ma, C. W. Glass, Z. Liu, T. Yu, O. O. Kurakevych, and V . L. Solozhenko, Ionic high-pressure form of elemental boron, Nature457, 863 (2009)
work page 2009
- [43]
-
[44]
A. Carvalho, M. Wang, X. Zhu, A. S. Rodin, H. Su, and A. H. Castro Neto, Phosphorene: from theory to applications, Nature Reviews Materials1, 16061 (2016)
work page 2016
-
[45]
P. D. Pancharatna, S. H. Dar, U. D. Chowdhury, and M. M. Balakrishnarajan, Anatomy of classical boron–boron bonding: overlap and sp dissonance, The Journal of Physical Chemistry A126, 3219 (2022)
work page 2022
-
[46]
H. Liu, A. T. Neal, Z. Zhu, Z. Luo, X. Xu, D. Tom ´anek, and P. D. Ye, Phosphorene: an unexplored 2d semiconductor with a high hole mobility, ACS nano8, 4033 (2014)
work page 2014
-
[47]
F. Xia, H. Wang, and Y . Jia, Rediscovering black phosphorus as an anisotropic layered material for optoelectronics and elec- tronics, Nature communications5, 4458 (2014)
work page 2014
-
[48]
T. V . Vu, D.-Q. Hoang, T. H. Ho, H. Van Chi, and K. D. Pham, The bp 3 monolayer as a high-capacity and rapid-diffusion an- ode for sodium-ion batteries: a first-principles study, Nanoscale Advances8, 673 (2026)
work page 2026
-
[49]
P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo,et al., Quantum espresso: A modular and open-source software project for quantum simulations of materials, Journal of Physics: Condensed Matter21, 395502 (2009)
work page 2009
-
[50]
H. J. Monkhorst and J. D. Pack, Special points for brillouin- zone integrations, Physical Review B13, 5188 (1976)
work page 1976
-
[51]
M. Methfessel and A. Paxton, High-precision sampling for brillouin-zone integration in metals, physical review B40, 3616 (1989)
work page 1989
-
[52]
Hamann, Optimized norm-conserving vanderbilt pseudopo- tentials, Physical Review B88, 085117 (2013)
D. Hamann, Optimized norm-conserving vanderbilt pseudopo- tentials, Physical Review B88, 085117 (2013)
work page 2013
-
[53]
M. Schlipf and F. Gygi, Optimization algorithm for the genera- tion of oncv pseudopotentials, Computer Physics Communica- tions196, 36 (2015)
work page 2015
-
[54]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Physical Review Letters77, 3865 (1996)
work page 1996
-
[55]
D. C. Liu and J. Nocedal, On the limited memory bfgs method for large scale optimization, Mathematical Programming45, 503 (1989)
work page 1989
-
[56]
Y . Nambu, Quasi-particles and gauge invariance in the theory of superconductivity, Physical Review117, 648 (1960)
work page 1960
-
[57]
U. Pinsook, N. Natkunlaphat, K. Rientong, P. Tasee, and J. Seeyangnok, Analytic solutions of eliashberg gap equations at superconducting critical temperature, Physica Scripta99, 065211 (2024)
work page 2024
-
[58]
J. Noffsinger, F. Giustino, B. D. Malone, C.-H. Park, S. G. Louie, and M. L. Cohen, Epw: A program for calculating the electron–phonon coupling using maximally localized wan- nier functions, Computer Physics Communications181, 2140 (2010)
work page 2010
-
[59]
S. Ponc ´e, E. R. Margine, C. Verdi, and F. Giustino, Epw: Elec- tron–phonon coupling, transport and superconducting proper- ties using maximally localized wannier functions, Computer Physics Communications209, 116 (2016)
work page 2016
-
[60]
F. Giustino, M. L. Cohen, and S. G. Louie, Electron–phonon interaction using wannier functions, Physical Review B76, 165108 (2007)
work page 2007
-
[61]
F. Mouhat and F.-X. Coudert, Necessary and sufficient elastic stability conditions in various crystal systems, Physical review B90, 224104 (2014)
work page 2014
-
[62]
E. Margine and F. Giustino, Two-gap superconductivity in heavily n-doped graphene: Ab initio migdal–eliashberg theory, Physical Review B90, 014518 (2014)
work page 2014
-
[63]
Y . Zhao, C. Lian, S. Zeng, Z. Dai, S. Meng, and J. Ni, Two-gap and three-gap superconductivity in alb 2-based films, Physical Review B100, 094516 (2019)
work page 2019
- [64]
- [65]
-
[66]
J. Seeyangnok and U. Pinsook, Phase stability and supercon- ductivity in hydrogenated and lithiated janus gaxs2 (x = ga, in) monolayers, Journal of Applied Physics138(2025)
work page 2025
-
[67]
H. Xie, Z. Huang, Y . Zhao, H. Huang, G. Li, Z. Gu, and S. Zeng, Strong electron–phonon coupling and multigap superconduc- tivity in 2h/1t janus mosli monolayer, The Journal of Chemical Physics160(2024)
work page 2024
-
[68]
J. Seeyangnok, U. Pinsook, and G. J. Ackland, Tunable super- conductivity in functionalized janus mosea (a= h, li) monolay- ers: Competition between lattice instability and electron pair- ing, ACS Applied Energy Materials (2026)
work page 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.