Maja
Feierabend
*,
Samuel
Brem
and
Ermin
Malic
Chalmers University of Technology, Department of Physics, 412 96 Gothenburg, Sweden. E-mail: maja.feierabend@chalmers.se
First published on 13th November 2019
Point defects, local strain or impurities can crucially impact the optical response of atomically thin two-dimensional materials as they offer trapping potentials for excitons. These trapped excitons appear in photoluminescence spectra as new resonances below the bright exciton that can even be exploited for single photon emission. While large progress has been made in deterministically introducing defects, only little is known about their impact on the optical fingerprint of 2D materials. Here, based on a microscopic approach we reveal direct signatures of localized bright excitonic states as well as indirect phonon-assisted side bands of localized momentum-dark excitons. The visibility of localized excitons strongly depends on temperature and disorder potential width. This results in different regimes, where either the bright or dark localized states are dominant in optical spectra. We trace back this behavior to an interplay between disorder-induced exciton capture and intervalley exciton–phonon scattering processes.
In this work, we shed light on the temperature-dependent optical fingerprint of TMDs spectrally below the bright exciton in presence of defects. We include direct photoluminescence (PL) from the bright exciton X and its localized states XLoc as well as phonon-assisted PL from momentum-dark excitons XD and the corresponding localized states XDLoc, cf.Fig. 1. We model both free and localized excitons on the same microscopic footing. While free excitons are formed by a Coulomb potential leading to the quantization of the relative motion of electron and hole, localized excitons are formed due to trapping into a disorder potential giving rise to a quantization of the exciton center-of-mass motion. We show that depending on the width of the trapping potential, the existence of localized states as well as exciton capture and recombination rates can be controlled. The calculated optical spectra agree well with recent experiments, observing localized states due to local strain18,22 or in presence of nanopillars.21,28 Moreover, depending on temperature and disorder potential characteristics we predict phonon- or localization-dominated regimes in PL spectra of TMDs.
To obtain access to the optical response of these materials, the knowledge of the number of emitted photons is crucial as it determines the steady-state photoluminescence . The dynamics of the photon number nq on the other side depends on the photon-assisted polarization17,34 resulting in . This microscopic quantity is a measure for optically induced transitions from the state (v,k1) to the state (c,k2) under annihilation (creation) of a photon . The states are characterized by the electronic momentum ki, and the band index λ = v, c denoting valence or conduction band, respectively. Note that we take into account the conduction band minima at the K, Λ, and K′ valley, which are crucial for the formation of momentum-dark exciton states.12
To account for excitonic effects, which are dominant in these materials,1,4,35,36 we project the photon-assisted polarization into an excitonic basis. We use the relation , where our original observable SvcqQ is projected to a new excitonic quantity SμQ that is weighted by the excitonic wave function φμq. Here, we have introduced the center-of-mass momentum Q = k2 − k1 and the relative momentum q = αk1 + βk2 with and with the electron (hole) mass me(h). The free excitonic eigenfunctions φq and eigenenergies εμ are obtained by solving the Wannier equation, which presents an eigenvalue problem for excitons.4,12,17,29,30 The corresponding Coulomb matrix elements are calculated using a Keldysh potential.4,37,38
To obtain the temporal evolution of SμQ(t), we exploit the Heisenberg equation of motion iħṠμQ(t) = [H,SμQ(t)],29,31 which requires the knowledge of the many-particle Hamilton operator H. The latter reads in this work H = H0 + Hc-l + Hc-phon + Hc-dis including the free carrier contribution H0, the carrier-light interaction Hc-l, the carrier-phonon coupling Hc-phon and the carrier-disorder interaction Hc-dis. To calculate the matrix elements, we apply the nearest-neighbor tight-binding approach39,40 including fixed (not adjustable) input parameters from DFT calculations of the electronic bandstructure.41
The crucial part of the Hamilton operator in this work is the carrier-disorder matrix element with a Gaussian disorder potential with σ denoting the full-width half maximum in real space, V0 the depth of the potential and R0 the position of the disorder. Assuming disorder centers far away from each other, R0 only leads to a phase in the momentum space. Although the Gaussian disorder potential is a simple approximation,42–44 in certain situations experimentally measured potential traps show a Gaussian-like behavior.19,22 Note that our theory is not limited to Gaussian potentials but could be extended to other potential forms, such as elliptical45 or nanobubbles.46 The investigated Gaussian potential in its general form can have its origin in a variety of physical phenomena, such as local strain gradients stemming e.g. from nanopillars on TMD monolayers21,22,28 but also in defect- or impurity-induced potentials on atomic level.47 In our work, we do not take into account variations in the trapping potential due to the atomic structure nor defect-induced changes in phonon characteristics. The main focus lies on providing qualitative microscopic insights into the impact of localized bright and dark exciton states on the optical response of monolayer TMDs.
By transforming the disorder potential into momentum space we find that disorder leads to a quantization of the exciton center-of-mass momentum, i.e. excitons are trapped into the disorder potential and cannot move freely anymore.43 Then, we project the phonon-assisted polarization into a localized excitonic basis with , where χμnQ are the localized excitonic wavefunctions. Note that all types of free excitons including bright KK or momentum-dark intervalley KΛ and KK′ excitons48 can be localized in a disorder potential, each with the localization quantum number n = 1s, 2s,…, cf.Fig. 1. The center-of-mass localization approach leading to new states below the free excitons is consistent with the observed mid-bandgap states in the electron–hole picture in density functional theory studies on defects.49–51
To obtain access to the wavefunctions and eigenenergies of localized excitons, we solve the eigenvalue problem
(1) |
(2) |
Here, we use the index b for the bright KK excitons within the light cone and n for the quantum number of the localized state. The equation resembles the well-known Elliott formula.53 The position of excitonic peaks is determined by the energy εbn, while the peak width is given by radiative (γbn) and non-radiative dephasing (Γb,n) due to phonons. The oscillator strength is determined by the optical matrix element, which reads in the new basis with localized and free exciton wave functions χ and φ, respectively, as well as the original momentum-dependent optical matrix element Mq.4 The delta distribution assures that only bright excitons within the light cone with Q = 0 are visible. In the limit of no localized excitons, we find , which results in the free-exciton photoluminescence formula.17
Moreover, our calculations show that the wave functions of localized excitons are an order of magnitude weaker compared to free excitons (Fig. 2(b)), and hence the oscillator strength of localized excitons |Mbn|2 becomes rather small. On the other side, the exciton occupation Nbn of the energetically lower localized states is high at low temperatures. The occupation is determined by the capture rate of excitons into the considered disorder potential. The capture process has been treated microscopically taking into account exciton–phonon scattering processes between free and localized excitonic states within the second-order Born–Markov approximation54 and applying the basis of orthogonal plane waves for localized states.55,56 More details can be found in the ESI.†
Interestingly, the reports in literature show that localized excitons are most pronounced in tungsten-based TMDs.16,57 The latter are characterized by the presence of momentum-dark KΛ and KK′ excitons that are located below the bright KK excitons.12,48 Their role in presence of disorder has not been investigated yet. Depending on their energetic position and the ratio between capture and phonon-assisted intervalley scattering rates, the free KK excitons can either first be captured in a disorder potential or first scatter into a dark intervalley exciton state (KΛ, KK′) and then be captured, cf.Fig. 1. To account for these additional dark exciton scattering channels, we extend eqn (2) including now phonon-assisted PL IdPL from the states ν = KK, KΛ, KK′:17
(3) |
Here, we have introduced the abbreviation . The position of excitonic resonances in the PL spectrum is now determined by the energy of the exciton εμm and the energy of the involved phonon ±ΩαQ. The sign describes either the absorption (+) or emission (−) of phonons. We take into account all in-plane optical and acoustic phonon-modes. Moreover, the appearing phonon occupation is assumed to correspond to the Bose equilibrium distribution according to a bath approximation.58 Since dark states can not decay radiatively, the peak width is only determined by non-radiative dephasing processes Γμn. The oscillator strength of phonon-assisted peaks scales with the exciton–phonon scattering element , where gνμQ is the exciton–phonon coupling element including free exciton wavefunctions.17,32
In this study, we focus on processes involving the bright KK excitons as initial states, as those states are expected to be the most dominant in PL. Note that in eqn (3), the sum over n describes the additive contribution of excited localized excitonic states. However, our microscopic calculations show that the intraexcitonic scattering of localized excitonic states (μn → μn′) appears on a much faster timescale than the capture processes itself. Hence, it is sufficient to take into account only the ground 1s state for the calculation of the optical response. This reduces the complexity to determine the exciton occupations Nνn appearing in eqn (2) and (3). Taking into account all scattering processes on a microscopic footing, we calculate the exciton densities within the states X, XLoc, XD, XDLoc, cf.Fig. 1. We would like to emphasize that we calculate both the exciton decay via intervalley scattering X → XD and the capture process X(D) → X(D)Loc on a microscopic level driven by the exciton–phonon coupling strength , cf. the ESI,† for more details.
Now we have all ingredients to investigate the photoluminescence of free and localized bright and dark excitonic states, i.e. IPL(ω) = IbPL(ω) + IdPL(ω). Fig. 3 shows the time-integrated and temperature dependent PL spectra for the exemplary material tungsten diselenide (WSe2) choosing a fixed disorder width of 30 nm and a disorder depth of V0 = 40 meV. These values are chosen by taking into account recent experimental measurements on nanopillars and deterministic local strain in TMD samples.21 The temperature dependence enters through exciton and phonon occupations in eqn (3). We observe a variety of peaks at low temperatures, which we can microscopically ascribe to: (i) the bright KK exciton X, (ii) the localized KK exciton XLoc (located about 40 meV below the X resonance), (iii) a series of peaks around 100–130 meV below the X resonance reflecting phonon-assisted PL from dark localized states (D1,2 = KΛ, KK′). In particular, we predict pronounced phonon side bands from the energetically lowest KK′ exciton involving TO/LO phonons (around 130 meV) and TA/LA phonons (around 115 meV). Our findings support the hypothesis of recent experimental findings, predicting photoluminescence resonances stemming from dark localized excitons.59,60
We find that with the increasing temperature the linewidth of all peaks increases (Fig. 3(b)) due to the enhanced exciton–phonon scattering. For free excitons, this is intuitive and well established in literature.61 However, for localized excitons the temperature behavior is still not well understood. Our calculations show (i) a temperature-independent contribution due to radiative decay and (ii) an linear increase in the linewidth due to the enhanced exciton–phonon scattering. The radiative part is determined by the wavefunction overlap and we find γ0 ≈0.5–1.0 meV, which is in good agreement with the experimental obtained value of 0.9 meV.60 The increase with temperature stems from the temperature dependence of the phonon-driven capture rate, cf. the inset in Fig. 3. Even at 0 K we find a contribution to the capture rate due to phonon emission processes. This agrees well with the overall temperature behavior of localized excitons observed in experiments, revealing stable localized excitons up to 100 K.60,62
Furthermore, we find that the localized bright exciton XLoc shows a fast intensity decrease with temperature. This can be explained by the behavior of the capture rate (cf. the inset in Fig. 3(b)): For temperatures below 20 K the capture of excitons is faster than the intervalley scattering with phonons, and hence the intensity of XLoc is high. The calculated capture rates are in good agreement with other theoretical studies.63 For temperatures above 20 K, the intervalley exciton–phonon scattering becomes faster, which means that excitons are more likely to scatter to the KΛ and KK′ states first, before they are captured in the respective localized states. This presents a crucial difference between bright and dark localized states, and their temperature behavior could be used to trace back their origin. Moreover, as XLoc signature disappears with temperature, the probability to occupy the dark states XD increases and from about 50 K we can see traces of phonon side17,32 bands from these actually dark states. We find additional peaks between 100 and 50 meV reflecting phonon emission and absorption processes from KK′ and KΛ excitons. To sum up, at low temperatures (<100 K) the PL is dominated by disorder-driven signatures, while at higher temperatures phonon-dominated signatures determine the optical response.
As expected from the calculated eigenenergies in Fig. 2(a), we observe clear spectral shifts of localized excitons towards lower energies with increasing disorder width, cf.Fig. 4. However, at about 60 nm a saturation is reached and the position of the peaks remains constant. This behavior can be understood as follows: for very narrow disorder potentials (in the range of 1 nm and below), i.e. σ < rB with rB as the exciton Bohr radius, excitons do not get trapped and hence we only observe phonon side bands (50–100 meV below X) from the free dark excitons (Fig. 4(b) upper panel). Here, we are in the phonon-dominated regime. As soon as the disorder potential is broad enough, excitons can be captured and localized bright excitons XLoc as well as phonon side bands of localized dark excitons appear (Fig. 4(b) lower panels). Here, we enter a disorder-dominated regime in the PL. When increasing the disorder width in real space, first excitons become stronger localized, however at some width a saturation is reached. This is due to the fact that with very wide potentials, excitons are not tightly confined anymore and tend to behave like free excitons again. They are characterized by a broad wavefunction in real space (and hence a narrow wavefunction in momentum space, cf.Fig. 2(b)).
Beside the spectral shift discussed above, we also observe a clear disorder-induced increase of the resonance linewidth of localized excitons. This is due to narrower exciton wavefunctions (cf.Fig. 2(b)) and hence stronger radiative dephasing, as the latter directly scales with |χμn|2. Interestingly, we also find that the intensity of localized dark and bright states clearly changes with the disorder width. This can be ascribed to the width dependence of capture processes, cf. the inset in Fig. 4(b). We find an overall enhanced capture efficiency at larger disorder widths due to enhanced overlap of the free and localized wavefunctions. As we increase the disorder width, the localized wavefunctions become narrower and hence more similar to the free exciton wavefunction resulting in a larger overlap. However, the capture rate is not only determined by the wavefunctions, but it also crucially depends on the energetic position of the highest localized state. At low temperatures, free excitons can scatter only by emission of a phonon into the energetically lower localized state. Therefore, the energetic difference between free and localized state has to be in the range of the phonon energy (≈20–30 meV). For larger disorder widths, the spectral distance between free and localized state increases and hence phonon-driven scattering becomes less probable. However, at some disorder width, an additional localized state emerges, which offers new scattering channels for phonons. Whenever this occurs, we observe a new peak in the capture rate (cf. different colors in the inset in Fig. 4(b)).
In conclusion, we have developed a microscopic approach to describe trapping processes of excitons in atomically thin 2D materials. We have calculated wavefunctions and energies of localized excitons depending on characteristics of the disorder potential. Using this knowledge, we have determined the optical fingerprint of the exemplary tungsten-based TMDs. We find a variety of pronounced peaks below the bright exciton stemming from both bright and momentum-dark localized states. The predicted signatures are strongly sensitive to temperature. We find a disorder-dominated regime at low temperatures and a free exciton phonon-dominated regime at higher temperatures. The gained insights can be extended to a broader class of 2D materials and might help to design tailored trapping potentials.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9cp05763c |
This journal is © the Owner Societies 2019 |