Superconducting properties of doped blue phosphorene: Effects of non-adiabatic approach 2
1. Introduction
Recent advances in the fabrication of two-dimensional (2D) materials reveal many
unprecedented phenomena for further investigation of the various properties of these
structures due to excellent electronic and phononic band energies [1–8]. One of
the distinctive features is the superconducting properties of atomically thin 2D
systems [9–18]. According to the BCS theory [19] the higher superconducting transition
temperature (Tc) can be achieved by having a remarkable electronic density of states
at the Fermi level, N(εF), with a large proportion of the Debye frequency, ωD. On the
other hand, to get a significant Tc, one has to reach a tangible electron-phonon coupling
(EPC) parameter λ= 2N(εF)hg2i/ωD. Nevertheless, a significant λarising from a large
value of N(εF) and hg2i, an average of EPC matrix elements over the first Brillouin zone
(1BZ), is indeed more favorable than a value of λarising from a small ωD. The main
reason for this is that a small ωDgives rise to a smaller ωlog, characteristic phonon
energy, and subsequently suppresses the Tc.
Based on these analyses, one can look for light materials with a large λto achieve
reasonable Tc. In the context of 2D superconductivity, there are several materials with
a sharp peak in their N(ε) spectrum near the valence band maximum. Examples of
such systems are 2D buckled structures based on nitrogen group elements such as blue
phosphorene (BLP), which is the lightest in this material group [20,21], and 2D M2X2
with (M= In, Ga and X= S and Se) [22–27]. For systems with large N(εF), there is
often the possibility that a large nesting, ξ, occurs at some specific phonon wave vectors
(qs) leading to a considerable reduction in the bare charge susceptibility, χ0, when the
temperature reduces. This character leads to an intense softening at specific qs in some
branches of the phonon spectrum, known as the Kohn anomaly [28,29], and in some
cases, a charge density wave (CDW) instability will appear in the form of the imaginary
modes for some particular phonon wave vectors (qCDW ) when the temperature is below
TCDW . Accordingly, access to such a high Tccan be prohibited by this instability [22,30].
Therefore, a more uniform distribution of λin terms of phonon wave vectors over the
1BZ is more favorable to gain large Tc, while, the formation of the CDW phase becomes
less likely.
For systems where the Kohn anomaly exists, one can further include the effects
of phonon spectrum normalization on temperature variations which are applied to the
Eliashberg function (α2F). This can lead to a temperature-dependent λ, so the presence
of the Kohn anomaly can alter λand ωlog simultaneously. Precisely, it is also important
to consider the non-adiabatic renormalization of the phonon spectrum, leading to a more
accurate determination of the formation of either the superconducting or the CDW
phase at low temperatures. The non-adiabatic phonon dispersions can be naturally
obtained by diagonalizing the phonon dynamical matrix related to non-adiabatic non-
self-consistent force constants at a physical temperature. Therefore, a self-consistent
solution should be employed to achieve a better Tc[31].
In this paper, based on first-principles calculations, we investigate the