Long living carriers in a strong electron-phonon interacting two-dimensional doped semiconductor

Carrier doping by the electric field effect has emerged recently as an ideal route for monitoring many-body physics in two-dimensional (2D) materials where the Fermi level is tuned in a way that -- indirectly -- the strength of the interactions can also be scanned. The possibility of systematic doping in combination with high resolution photoemission has allowed to uncover a genuinely many-body electron spectrum in single-layer MoS2 transition metal dichalcogenide, resolving three clear quasi-particle states, where only one state should be expected if the electron-phonon interaction vanished. Our analysis combines first-principles and consistent complex plane analytic approaches and brings into light the presence and the physical origin of two gaps and the three quasi-particle bands which are unambiguously present in the photoemission spectrum. One of these states, though being strongly interacting with the accompanying virtual phonon cloud, presents a notably long lifetime which is an appealing property when trying to understand and take advantage of many-body interactions to modulate the transport properties.

trons do not have sufficient energy for emitting any phonon, and therefore they appear long lived.
However, they are allowed to emit and reabsorb phonons in virtual processes, which produces an entire phonon cloud around electrons having the effect of augmenting their effective mass, similar to polaron states in insulators 11 . On the other hand, for electrons above ω 0 the emission of phonons is now allowed, the effect being that virtual processes are less probable in favour of real emission of phonons leading to a decreasing of their lifetime and effective mass. This idealized picture will be useful for understanding the more intricate situation in MoS 2 and questions the simple view about quasi-particle properties themselves, since the interaction with phonons produces two different electron states with radically different properties. A bit more technically, as the energies of quasi-particle states are determined by the poles of the electron propagation amplitude or Green's function (G(k, z)) 12 , the presence of several poles in this function evidences the existence of as many quasi-particle states (Fig. 1b), with very different physical properties in terms of lifetime and dispersion ( Fig. 1c-d). Similarly, in a simple system consisting of two optical phonon modes with different energies, the electron band would split twice when reaching the energy corresponding to each phonon mode.
There have been many good examples of ARPES measurements on metal surfaces 13-15 , high-temperature superconductors 16 , and doped polar insulators 2,17 , where strong deviations from the single quasi-particle picture have been observed and attributed to the electron-phonon interaction. However, the recent measurements of Kang et al. 3 show that MoS 2 is probably the first system where a double splitting in the ARPES spectrum has been observed unambiguously, therefore the importance of understanding the physical nature of these broken bands. These results go also hand in hand with the spin-valley locking protected superconductivity found in this system 4 , and while some experimental works illustrated the superconducting state as due to the filling of K(K ) valleys 5-7 , this interpretation clashes with the picture which holds that superconductivity emerges as soon as the Fermi energy (E F ) crosses the bottom of the Q(Q ) valleys [18][19][20] . In fact, the onset of superconductivity coincides with an outstanding enhancement of the electron-phonon coupling strength 18 , stemming from phonon-mediated intervalley interactions on the Fermi surface (FS) [18][19][20] , also in agreement with the strong softening of several phonon modes observed in the gate-induced superconducting monolayer MoS 2 19 . We analyse these points and more specifically the case of the observed multiple band splitting 3 in MoS 2 , where the experiments seemed to point out a direct interpretation in terms of multiple-phonon excitations.
We made an in-depth theoretical analysis of the electron-phonon interaction in electron-doped monolayer MoS 2 by means of first-principles calculations in order to shed light on several aspects of the electron-phonon coupling in this system (see also Methods). We chose a doping carrier density of n 2D = 9 × 10 13 cm −2 , for which the conduction-band minima at K(K ) are filled and almost spin-degenerate (Fig. 2a)  We have computed the ab initio electron spectral function A(k, ω) including the electronphonon effects without any adjustable parameter (Methods). This function displays two sharp band-splittings for electron momenta close to K at ω M A and ω M O binding energies and width 22 meV ( Fig. 3a-b), each of these band-splittings evoking the one sketched in Fig. 1a for the simplified Einstein model. The calculated spectral function is in close agreement with the measurements by Kang et al. 3 and describes the three spectral bands observed in experiment. We find that from the two concentric bands around the K valley, it is the outer one (spin-up) which shows the most intense signatures of the electron-phonon coupling, as easily recognized when considering spin conservation arguments (Fig. 2 d). This is also the reason why the inner spin-locked state shows much weaker spectral features. Let us now focus on the strongly interacting outer spin-split band.
For this state, the frequency dependent imaginary part of the electron self-energy (ImΣ kj (ω)) at momentum k = k A close to K (Fig. 2a) shows a rather uncommon rectangular shaped double structure of width ∼ |E Q | = 22 meV separated by a narrow window at ω ∼ 42 meV with almost vanishing value (Fig. 3c right). These two rectangular shapes have onsets precisely at ω M A and ω M O binding energies and they are easily rationalized in terms of energy conservation for the corresponding phonon emission processes connecting states close to Q with those near K (Fig. 2b).
It is actually the large DOS at the occupied Q pockets (Fig. 3c left) which enhances the phase-space of these strongly-interacting scattering processes, yielding the maximum values of ImΣ k A ↑ (ω) (yellow shaded areas in Fig. 3c). The dip in the imaginary part of the self-energy appears as long as the condition ω M O − ω M A > |E Q | holds, since the quasi-particles near k = k A with energies between the two maxima can not scatter to Q valleys. This means that the most uncommon spectral feature found in ARPES close to ∼ 42 meV, also found in our calculations (Fig. 3b), if demonstrated to be a well-defined quasi-particle, it would correspond to a very long-lived state even if it is strongly interacting. This is a unique characteristic of MoS 2 , arising from the fact that the relevant doped valleys have different binding energies and that there are basically two relevant phonon modes connecting precisely those valleys. Recall that in ordinary metals, where the DOS is practically constant for typical phonon energies, the imaginary part of the electron-phonon selfenergy is a monotonically increasing function, which is radically different to what is observed in MoS 2 . Note that thermal and electron-electron interactions are expected to slightly soften the spectral features obtained considering only the electron-phonon coupling. An estimation of the electron-electron interaction considering RPA 21 adds a slowly monotonically increasing function yielding a broadening of Γ ee ∼ 2.8 meV even for electron energies close to ω ∼ 40 meV.
It's worth looking a bit more closely at the general properties of the quasi-particles even if it's succinctly, in order to see if the observed spectral features can be defined as actual quasiparticles. These find their appropriate mathematical definition as poles of the electron Green's , that is, the solutions of the so-called Dyson equation The key point here is that due to the nonlinear character of this equation it may lead to several solutions, even for a single unperturbed state with energy k j . The consistent use of the complex plane allows to treat the Green's function in the entire plane and the energy (E qp ) and the lifetime broadening (Γ qp ) of possible quasi-particles are compactly expressed as an ordinary complex number z qp n = E qp − iΓ qp . Then the Dyson equation leads to a pair of coupled non-linear equations 22 , Capturing the influence of the quasi-particle lifetime broadening on its shifted energy and vice versa, as above, is not the most standard procedure but appears absolutely crucial when trying to understand the spectral features of many strongly interacting systems in terms of elementary excitations, as shall be the case also in MoS 2 .
We have solved Eq.(1) for monolayer MoS 2 , and the results are shown in Fig.3d. For the sake of clarity, we focused on the same momentum region as in Fig.3b, and only on the strongly interacting outer spin-split band. Close to the Fermi momentum k F , an Engelsberg-Schriefferlike 10 state appears with a strongly renormalized dispersion, and it is denoted as the n = 1 solution.
Far enough from k F , we find a dispersive and damped state identified as the n = 3 solution. An important conclusion is that for some intermediate values of the momentum we find an additional solution (n = 2) with an important spectral weight which is practically flat, and it is therefore a strongly interacting state which tends to localization. However, this state appears long lived as it lies just in the energy window where the imaginary part of the self-energy has almost a gap (Fig.3c). More specifically, the electron-phonon limited lifetime broadening at this energy range is almost negligible Γ qp n=2 ∼ 0.35 meV. Considering the Dyson equation in the complex plane allows also to associate a precise spectral weight to each quasi-particle state, in a systematic way, and is simply given by the residue of the poles (Z qp n = 1/ 1 − Σ (z qp n ) ) evaluated at complex quasi-particle energies. It is not therefore necessary to visually analyse the spectral function A(ω) in order to check for the energies of possible quasi-particle states and/or their relative importance, as these are given directly by z qp and Z qp n . This allows to define the coherent part of the spectralfunction systematically ( ), and explicitly broken down into separate contributions from each quasi-particle pole. Note for a moment the severe consistency condition n Re (Z n ) 1 imposed here as the integral of A qp (k, ω) must be of the order but less than unity. Said in pass, this result is not obtained by any means when only the real part of the self-energy is considered for calculating the dispersion of quasi-particles (see Supplementary   Information). When the contribution of each separate pole A qp n=1,3 (k, ω) is plotted as in Fig.3f, the comparison with the full spectral function (Eq.7 in Methods) in Fig.3e is strikingly good. In this figure it is also shown how strongly the spectral weight from one quasi-particle into others is transferred as a function of k, even when for some values of momentum all the three manybody solutions are present. Altogether, the obtained quasi-particle band structure -in its complex version-almost perfectly resembles the three-peak and double-gap structure observed in ARPES. Therefore, the threefold band structure observed in experiment corresponds certainly to quasiparticle states and not to multiple-phonon (high order) processes nor to side-bands or satellites without a clear physical meaning.
These findings provide a theoretical explanation for the singular spectral features observed on monolayer MoS 2 3 in terms of three elementary many-body quasi-particle states. Close to the K valley, the available scattering phase-space appears severely restricted for some narrow energy windows, with the result that one quasi-particle branch (n = 2) appears exceptionally long lived, even when the strong coupling induces a practically flat dispersion for this band. Flat dispersion indicates a sort of real space localization property of these states and, therefore, the effective interactions between them should be expected to be profoundly modified. As the transport and superconducting properties 4-7 depend so dramatically on the lifetime, dispersion and effective interactions between elementary excitations, we believe that these results may serve as a guidance to understand, explore and eventually take adavantage of many-body interactions.

Methods
First-principles calculations. The first-principles calculations were performed within the noncollinear density functional theory (DFT) 23,24 and density functional perturbation theory (DFPT) 25 with fully relativistic norm-conserving pseudopotentials as implemented in the QUANTUM ESPRESSO package 26 and using the Perdew-Zunger local density approximation parametrization for the exchange-correlation functional 27 . In order to correctly describe the ground-state electronic and phononic structures of the single-layer MoS 2 , we used the experimentally measured in-plane bulk lattice parameter a = 3.16Å 28 and a vacuum layer of five times the lattice parameter, which is large enough for avoiding any interplay between adjacent layers. Carrier doping effects were simulated by the addition of excess electronic charge into the unit cell, which is compensated by a uniform positive jellium background. All atomic forces were relaxed up to at least 10 −6 Ry/a.u..
We used a 36 × 36 Monkhorst-Pack grid for the self-consistent electronic calculations, while the lattice vibrational properties were evaluated on a coarse 9 × 9 q-point mesh.
Electron-phonon interaction calculations. The self-consistent electron-phonon magnitudes were computed considering doping-sensitive full-spinor electron states, as well as doping-sensitive phonon states and deformation potentials. In a first step we considered a coarse mesh of 9 × 9 for phonons (q) and a 18 × 18 mesh for electron states (k), but for fine integrals, we used the Wannier interpolation scheme [29][30][31] , which allowed us to consider 10 7 points in the Brillouin zone for electrons and 10 6 for phonons. The Fermi surface integrated squared matrix elements (Fig. 2c), equivalent to a sort of electron-phonon weighted nesting function defined for the phonon branch (ν) at q, is defined by the following integral: where k j is the single-particle bare eigenvalue of the electron state with momentum k in the band j, ω q ν is the frequency of the lattice vibrational normal mode with momentum q and mode ν, N (E F ) the DOS at E F , and g ν ij (k, q) the electron-phonon matrix element connecting the |k, j and |k + q, i electron states by a phonon |q, ν . This quantity allows us to identify the relevant phonon modes interplaying with electrons at E F and offers a measure of the strength of the coupling for each phonon mode |q, ν .
The mass enhancement parameter or Fermi surface averaged electron-phonon coupling strength is defined as 9 : This factor reveals how much the effective mass is enhanced m = m 0 (1 + λ) at the Fermi level due to the electron-phonon interaction. The state |k, j resolved electron-phonon coupling strength ( Fig. 2d) was evaluated as: The semi-empirical McMillan-Allen-Dynes formula for estimating the superconducting critical temperature T c is defined as, where ω log is a logarithmic average of the phonon frequencies and µ * is a parameter describing the Coulomb repulsion. We estimated the superconducting critical temperature of the doped monolayer MoS 2 for some values of µ * within the range of 0.1 − 0.3, obtaining values in the range of 4 − 8 K, with a calculated electron-phonon coupling of λ ∼ 0.63. This estimation is in agreement with the experimentally measured critical temperature of 8.5 K at a 2D-carrier density n 2D = 9 × 10 13 cm −2 4-7 .
Electron self-energy and spectral function. The ab inito expression for the Fan-Migdal electron self-energy valid at T = 0 K used in our calculations was 29 : where ω represents the quasi-particle energy in the real-energy axis, f denotes the Fermi-Dirac occupation factor, and η is a real positive infinitesimal. The electron self-energy as calculated in Eq. 6 allows to obtain the spectral function for real frequencies 8,9,29 : Analytic continuation of the electron self-energy. In order to find solutions of the complex quasi-particle equation Eq. 1, we need to perform an analytical continuation of the electron selfenergy Σ(ω) from the upper half complex plane into the lower half. First, we calculated the selfenergy from first principles only for real frequencies (Eq. 6), and next we generalized the method outlined in Ref. 22 to handle self-energies without assuming particle-hole symmetry. For doing so, we used the Kramers-Kronig relation integrated by parts, As it is known, a direct substitution of ω by z in Eq. 8 is not valid, as the branch-cuts introduced by the log term makes the direct numerical integration inappropriate. However, a piecewise polynomial interpolation of dIm(Σ(ω)) dω and a subsequent analytical integration is a possible numerical strategy. We used a cubic spline interpolation, ensuring the continuity of the first and second derivatives all over the ω axis. represented by a dotted line (ω 0 ). The electron decay processes by phonon emission are energetically allowed only for electrons above ω 0 (red background). Therefore, the situation is completely different for electrons below this energy range, though even there, quantum field laws allow for virtual excitations of phonons. The result is that the coupling produces two excitation branches E qp 1 and E qp 2 , whose dispersions are represented by blue and red solid lines, respectively. Fixing the electron momentum at momentum at k 0 (dashed line in a), one obtains a spectral function with two peaks at real frequencies (b) corresponding to the above pair of bands. Looking at the complex frequency plane, these two excitations are traced back to poles of the electron Green's function and their position in the complex plane determines the renormalized energy (E qp ) and lifetime broadening (Γ qp ). When most of the spectral function (coherent part) is recovered from the superposition of the separate contributions of the poles, one may say that a multiple quasi-particle picture is valid. c, Electrons below ω 0 are strongly renormalized and tend to localize as they are followed by a dense virtually emitted phonon cloud. d, On the contrary, the higher energy band is energetically allowed to emit phonons and acquires a more extended character and a lighter effective mass.  to K (Fig. 2a). The onsets of the rectangular maxima are at ω M A and ω M O , the energies of the acoustic and optical phonons at q = M, while their width is related to the enhanced DOS at the occupied Q pockets (yellow shaded area in left panel). d, Dispersion of the three quasi-particle poles found for the outer band. The real part of the poles -quasi-particle energies -with respect to the momentum are shown by the blue (n = 1), green (n = 2) and red (n = 3) dots, respectively. The length of the bars represent the spectral weight of each pole, given by the real part of their residues Z qp n , while the imaginary part of the residues are represented by the rotation of the bars, Im(Z qp n ) = 1 giving a rotation of θ qp n = π radians. e,f Comparison between the full ab initio spectral function (see Eq. 7 in Methods) (e), and the contribution coming from each complex quasi-particle pole (f), shown by different colors following the same convention as in d.