Next Article in Journal
Fabrication of Scaffold-Based 3D Magnetic Nanowires for Domain Wall Applications
Previous Article in Journal
Boron-Doping Proximity Effects on Dislocation Generation during Non-Planar MPCVD Homoepitaxial Diamond Growth
Previous Article in Special Issue
Properties, Preparation and Applications of Low Dimensional Transition Metal Dichalcogenides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nature of Excitons in Bidimensional WSe2 by Hybrid Density Functional Theory Calculations

Dipartimento di Scienza dei Materiali, Università di Milano Bicocca, via R. Cozzi 55, 20125 Milano, Italy
*
Author to whom correspondence should be addressed.
Nanomaterials 2018, 8(7), 481; https://doi.org/10.3390/nano8070481
Submission received: 7 June 2018 / Revised: 27 June 2018 / Accepted: 27 June 2018 / Published: 29 June 2018
(This article belongs to the Special Issue Synthesis, Structure and Applications of 2D Nanomaterials)

Abstract

:
2D tungsten diselenide (2D-WSe2) is one of the most successful bidimensional materials for optoelectronic and photonic applications, thanks to its strong photoluminescence properties and to a characteristic large exciton binding energy. Although these optical properties are widely recognized by the scientific community, there is no general understanding of the atomistic details of the excitonic species giving rise to them. In this work, we present a density functional theory investigation of excitons in 2D-WSe2, where we compare results obtained by standard generalized gradient approximation (GGA) methods (including spin-orbit coupling) with those by hybrid density functionals. Our study provides information on the size of the self-trapped exciton, the number and type of atoms involved, the structural reorganization, the self-trapping energy, and the photoluminescence energy, whose computed value is in good agreement with experimental measurements in the literature. Moreover, based on the comparative analysis of the self-trapping energy for the exciton with that for isolated charge carriers (unbound electrons and holes), we also suggest a simplified approach for the theoretical estimation of the excitonic binding energy, which can be compared with previous estimates from different approaches or from experimental data.

Graphical Abstract

1. Introduction

Transition metal dichalcogenides (TMDCs) have been studied for a long time [1,2]. The structure of the progenitor of this class of materials, i.e., molybdenum disulfide (MoS2), was determined in 1923 [3]. However, only after the successful discovery of graphene, which has paved the way for the incredibly fast development of the research field on 2D materials, have the bidimensional forms of TMDCs become the object of intense study and activity, especially thanks to their interesting electronic and optical properties [4,5,6]. Indeed, they present strong photoluminescence (PL) and a large exciton binding energy, which make them ideal materials for optoelectronic and photonic applications [7,8,9,10,11,12]. TMDCs are different from graphene and hexagonal boron nitride (h-BN) which have semi-metallic and insulating characters, respectively, as they are semiconductors with a direct band gap in the visible energy range [4,5]. The electronic structure of TMDCs is far more complex than that of graphene or h-BN because of the metal d state’s interaction with the chalcogen p states. This feature can explain the excellent quantum yield of photoluminescence for the single layers that are, thus, particularly suitable materials for LEDs (light-emitting diodes) [13,14]. Due to the low dielectric screening, there is a large Coulomb interaction between the charge carriers (electron and hole), resulting in an unusually high exciton binding energy [8,9,15,16,17]. Understanding the nature of excitons in two-dimensional (2D) TMDCs is of key importance to make use of their optical and charge transport properties in optoelectronic and photonic applications.
The 2D tungsten diselenide (WSe2), which is the object of the present work, belongs to the 2D space group P 6 ¯ m 2 and comprises three atomic layers: The two external layers are triangular sublattices of selenium atoms, whereas the middle layer is a triangular sublattice of tungsten atoms (Figure 1a). The unit cell is hexagonal and is defined by one lattice parameter. Each W atom presents a trigonal prismatic coordination, leading to d states splitting shown in Figure 1b and to the semiconducting nature of this material. Although the electronic structure of 2D-WSe2 has been extensively studied from both an experimental and theoretical point of view, there is still some debate on the true nature of the band gap as we will discuss in detail below [18,19]. Calculations were performed at both density functional theory (DFT) and GW levels of theory [5,9,20,21,22]. The main and most peculiar features are related to the presence of two conduction band (CB) minima (at Q and K points in the Brillouin zone) and two valence band (VB) maxima (at K and Γ points in the Brillouin zone). Progressively reducing the number of layers, Zhao et al. observed a band gap inversion from indirect to direct when going from bi- to the mono- (triatomic) layer [13]. Correspondingly, a very large increase in the experimental photoluminescence intensity is registered, which could be assigned to the shift from an indirect to a direct band gap. However, a scanning tunneling spectroscopic (STS) study suggests a minimum indirect band gap (2.20 eV) and the presence of a direct band gap at a very close energy (2.12 eV) [19]. Further complication comes from the sensible difference between optical and electronic gaps [4,5] due to a large exciton binding energy. A typical intense photoluminescence peak deriving from the exciton radiative recombination can be observed at about 1.6–1.7 eV [7,13,23]. In general, when there is a strong coupling between the just formed free exciton and the atomic lattice, self-trapping phenomena takes place, leading to large lattice (polaronic) distortions around the exciton. This atomic relaxation is mainly responsible for the Stokes-shift between the absorption and emission peaks.
The accurate theoretical description of the electronic properties of the ground and excited states of 2D-WSe2 is not a simple task. First, local density approximation (LDA) or generalized gradient approximation (GGA) functionals are not well suited to correctly reproducing a semiconductor band gap. One should resort to hybrid density functionals or DFT+U [24,25,26] or GW [27] approaches. The excitation or optical absorption is an even more demanding task. Available theoretical approaches are time-dependent DFT (TD-DFT) [28,29] or the Bethe-Salpeter equation (BSE) [30,31]. The latter includes excitonic effects which are crucial for these types of materials. However, electron-phonon coupling and local relaxation effects cannot be easily included. In this work, we apply a simplified but effective approach that we have used before to model free excitons in semiconducting materials, the exciton trapping process, and finally, the exciton radiative sections. The exciton wave function is computed as the ground triplet recombination [32]. The energy quantities associated with those processes are estimated quite accurately by means of hybrid density functional calculations (HSE06) [33] as will be discussed in the next state of the systems. Then, atomic positions are relaxed to allow for the polaronic distortion. Finally, the spin constraint is released, so that the system can go back to the electronic ground state through a vertical radiative transition. We compare the HSE06 with Perdew-Burke-Ernzerhof (PBE) results, as implemented in the CRYSTAL14 (CRY14) code [34,35]. Moreover, we compare the PBE results based on atomic localized basis sets (CRY14) with those based on plane wave basis sets (Quantum ESPRESSO, QE) [36]. For the latter, we can even introduce the spin-orbit coupling effect, which could be relevant in some respects. To gain a deeper understanding of the physical processes involved, we also investigate self-trapping of separated charge carriers (electron and hole). This will allow, for the first time, a theoretical estimation of the binding energy in a self-trapped exciton in 2D-WSe2 to be provided and to be compared to experimentally obtained values in the range between −0.4/−0.8 eV [15,16].

2. Computational Details

All the calculations were carried out using the CRYSTAL14 (CRY14) package [34,35] based on density functional theory (DFT) where the Kohn−Sham orbitals are expanded in Gaussian-type orbitals. For selenium and tungsten, we adopted the well-assessed effective core pseudopotential (ECP) techniques [37,38]. The small-core ECP derived by Hay and Wadt [37] was chosen for W, leaving only the 5p, 5d, and 6sp valence electrons to be treated explicitly. For the tungsten valence shell electrons, we used a modified Hay Wadt double-ζ basis set [37], as adopted before for WO3 [39]. The ECP and valence shell electrons for Se are treated as reported previously, where the 3s, 3p, 3d, 4s and 4p electrons are viewed as valence electrons [40,41]. Hybrid functional HSE06 [42] was used for both geometry optimization and electronic properties calculations. For the WSe2 monolayer, the hexagonal unit cell with the symmetry group of P 6 ¯ m 2 was fully optimized with respect to both atomic coordinates and cell parameters. K point mesh with uniform spacing of 0.015 Å−1 was used to get the converged total energy and band gap. The convergence criterion of 0.023 eV/Å for force was used during geometry optimization, and the convergence criterion for total energy was set at 10−7 Hartree for all the calculations.
Different sizes of supercell for the monolayer 2D-WSe2, including 2 × 2, 4 × 4, 5 × 5, 6 × 6, 8 × 8 supercells, were first tested for the exciton species, which suggested that the 6 × 6 supercell is big enough not to observe spurious effects due to interaction with repeated images. Therefore, we used this throughout the work. To model exciton and charge carriers, triplet and doublet state calculations were performed, respectively. Energy of absorption (ΔEABS) is calculated as the energy difference between the triplet and singlet states of the neutral system in its ground state geometry. Self-trapping energies (ΔEST) are calculated as the energy difference between the isolated exciton (or charge: e or h+) in the trapping geometry and in the ground state geometry. The vertical emission energy (photoluminescence) (ΔEPL) is obtained as the difference between the energy of a neutral system in the triplet and singlet states in trapping geometry. The distortion energy (ΔEDIST) is defined as the difference between the energy of singlet neutral system in trapping geometry and in ground state geometry. Symmetry was removed for these calculations to allow the distortion in geometry induced by exciton or charge trapping.
To compare the results, including spin-orbital coupling, from different basis sets, some calculations are also performed using the plane-wave-based QE package [36]. The projector augmented wave (PAW) potentials [43] were adopted to describe the electron-ion interactions with W (4f, 5s, 5p, 5d, 6s) and Se (3d, 4s, 4p) treated as valence electrons. The exchange and correlation interaction were described by the PBE [44] functional within the GGA. Energy cutoffs of 55 Ryd and 220 Ryd (for kinetic energy and charge density expansion, respectively) were adopted for all calculations. The convergence criterion of 0.026 eV/Å for force was used during geometry optimization, and the convergence criterion for total energy was set at 10−6 Ryd for all the calculations.

3. Results and Discussion

3.1. Atomic and Electronic Structure of 2D-WSe2 by Hybrid Density Functional Calculations

We have first assessed the performance of hybrid density functionals in the description of 2D-WSe2 since those data are missing in the literature. For this preliminary assessment, we have compared different hybrid functional (HSE06, B3LYP [45,46] and B3PW [45,47]) with PBE results, using the CRY14 code. 2D-WSe2 structure belongs to the P 6 ¯ m 2 hexagonal spatial group. The lattice parameters are summarized in Table S1 and are rather close to the experimental value from X-ray diffraction (XRD) of the bulk phase (3.282 Å). In particular, HSE06 outperforms all the others (3.278 Å) with an error of only 0.13%. The correct reproduction of the lattice constant is particularly relevant in the present case since it has been shown that the electronic structure of 2D-WSe2 is very sensitive to this parameter. Some artificial compressive/tensile stress leads to an increase/decrease of band gap or to contrasting results on its direct or indirect nature [18].
Although the W d-states splitting caused by the trigonal prismatic coordination can well explain the basic electronic properties, most of the interesting properties of 2D-TMDCs come from the peculiar dispersion of the band states. In Figure 2, we report the band structure of 2D-WSe2 as obtained with the different density functionals. All the hybrid density functionals, in contrast with PBE (Kv–Qc), present a minimum direct (Kv–Kc) gap. Moreover, PBE largely underestimated the band gap value (1.71 vs. 2.12 eV [19]) as expected. Although the presence of a direct gap is a hallmark of 2D-TMDCs, for the specific case of 2D-WSe2, the crossover between direct/indirect band gap is still under debate. This could be easily rationalized considering the very tiny difference in energy between the two conduction band minima in Q and in K. In HSE and PBE they are practically degenerate since the respective differences of +0.03 and −0.01 eV are within the error of the methods. These differences are very sensitive to the lattice parameter: Note that B3LYP, whose lattice parameter is the largest, provides the largest ΔE (Qc–Kc). From the experimental side, Zhang et al. tried to shed some light on the issue by performing some electronic transport measurements with STS [19]. This work has revealed the existence of both an indirect (Kv–Qc) gap at 2.12 eV and of a direct (Kv–Kc) gap at 2.20 eV, with an error range of ~0.1 eV. These data are in good agreement with the HSE06 values as we have just discussed above. HSE06 also provides the best estimate of the ΔE (Kv–Γv) of 0.49 eV versus the experimental measurement of 0.64 eV [19]. A direct band gap for 2D-WSe2 is quite surprising, given the extraordinary PL properties [7]. However, the existence of an almost degenerate indirect band gap can justify the observed phenomenon and is in line with accurate DFT calculations in this work. These results fully agree also with the analysis by Hsu et al. of heterojunction-induced strain effects on PL energy and intensity [18].
It is noteworthy now to compare the PBE band structure by CRY14 to the corresponding one by QE (Figure 2 vs. Figure 3a). They are essentially identical, which suggests that there is no basis set effect on the band structure. More importantly, it is then relevant to analyze what happens upon the introduction of spin–orbit coupling (SOC) in the PBE/QE calculations (Figure 3b) which cannot be done for CRY14 calculations. We clearly observe a splitting of the α and β states, especially in the valence band at the K point (by 0.46 eV), but not in Γ. In the conduction band, the largest splitting is at the Q point (0.23 eV) and very feebly in K. Our observations are in perfect agreement with a previous work where they have been largely discussed [13,15,16,31,48,49,50,51,52,53,54,55,56,57].
In the last part of this section, we will discuss the details of 2D-WSe2 electronic structure in terms of total and projected density of states (TDOS and PDOS, respectively), which are reported in Figure 4a. The states of the valence and conduction bands are essentially made up by W 5d and Se 4p. The contribution from both elements to both VB and CB suggests a strong covalent character of the W–Se bond in WSe2 and, consequently, a strong d-p hybridization. Moreover, we can clearly distinguish a first large peak in the VB, from the Fermi level (set to the zero energy) to −1.5 eV, where the W contribution is predominant, whereas, from −1.5 eV down in energy, the Se contribution becomes the largest, indicating some polarity in the bond, too. The d-p mixing is crucial for the absorption properties of this material since it makes the transition permitted, in contrast to what would happen if top of the VB and bottom of the CB had both a d character since d-d transition are forbidden for the electric dipole selection rules. In Figure 4b, we report the projected density of states on the different d states of the W atoms following the trigonal prismatic d splitting. The dxz,yz contribution is negligible in the energy range of the VB close to the Fermi level, where the dz2 and the dxy,x2−y2 are more relevant, in agreement with the crystal field splitting described in Figure 1. At a deeper level of analysis, we can even provide the main contribution of the d states to the highest occupied band and to the lowest unoccupied band, in particular, where at the maximum and minimum points, respectively. The data are presented in Table 1. We can clearly observe that the bands in Kc and Kv are primarily made up of the W dz2 and the dxy,x2−y2 states. A small contribution comes from the Se px and py states but not from pz ones. On the contrary, pz states contribute to the states at the top of the valence and at the bottom of the conduction band at the Qc and Γv points, causing the different gap character (direct) for the monolayer system with respect to the bulk.

3.2. Self-Trapped Charge Carriers in WSe2

3.2.1. Self-Trapped Electron

An excess electron was added to the neutral and stoichiometric monolayer (here for the 6 × 6 supercell model). It was injected into the conduction band. It was found to be non-uniformly delocalized on several W atoms, as can be clearly observed in Figure 5a. However, when using a hybrid functional (HSE06/CRY14) and upon atomic relaxation, the electron undergoes a self-trapping process with an energy gain of −0.20 eV ( Δ E S T , e ), as shown below:
  W S e 2 + 1 e [ W S e 2 ] R e l a x a t i o n [ W S e 2 ] t r a p
  Δ E S T , e = E [ W S e 2 ] t r a p E [ W S e 2 ]
During the relaxation, some symmetry is lost, even though we can still observe a perpendicular plane of symmetry, which is represented by a dotted line in Figure 5b. The electron is then found to be localized on three W atoms of one hexagonal ring of the monolayer (65% on one W and 11% on each of the other two, summing to 87% of the unpaired electron density (Figure 5b). There is a further 8% on the W atom lying in the plane of symmetry of the system (Figure 5b). As a point of difference, any attempt to localize the electron by the standard PBE functional was unsuccessful, due to the intrinsic electron self-interaction error, which causes an excessive delocalization. The polaronic distortions leading to the self-trapping with the hybrid functional involve some structural rearrangement reported in Figure 5e to be compared with Figure 5d. We observe an elongation of the W–Se bonds, for those W that host the extra electron, an elongation of the distance between the W with 65% of the extra electron and the Se atom on the opposite side of the hexagon as shown in Figure 5e, and a shortening of the distance between the two Se atoms bound to the W with the 65% of the extra electron.
A second solution for the self-trapped electron exists, which was obtained starting from a different initial atomic configuration and is only 5 meV higher in energy than the one we have just described. For this second solution, some spin polarization is observed, as it is evident in Figure 5c. Again, one plane of symmetry perpendicular to the system layer is kept, however, the electron is mainly localized on the two W atoms at the two sides of the plane (up a component of 50–55%) with some down component (22%) on the third W of the hexagon. The polaronic distortion is similar but with some differences. When we compare Figure 5e with 5f, we observe an elongation of the W–Se bonds again for those W that host the extra electron, an elongation of the distance between the two W atoms with 50–55% of the extra electron, and a shortening of the distance between the W and the Se atoms in the plane of symmetry.
The band structure of the two solutions is very similar (Figure 6) and presents a clear flat state in the gap of the spin up component that is occupied by the extra electron. This state features a strong W d character (see Figures S1 and S2 in the Supplementary Literature) and is −0.57 eV below the bottom of the conduction band for the first solution and −0.50 eV for the second one.

3.2.2. Self-Trapped Hole

The removal of one electron from the stoichiometric neutral systems leads to the formation of an electron hole, as shown below:
  W S e 2 [ W S e 2 ] + + 1 e R e l a x a t i o n [ W S e 2 ] t r a p +
  Δ E S T , h + = E [ W S e 2 ] t r a p + E [ W S e 2 ] +
Any attempt to localize this hole, with supercells models of different size and with both the standard PBE and the hybrid HSE functionals, was unsuccessful (see Figure 7a for the 6 × 6 supercell and HSE). The spin density is found to be fully delocalized on all the W atoms of the model. The relaxed geometry is almost identical to that of the neutral system and the self-trapping energy ( Δ E S T , e ) is only −0.05 eV. The band structure (Figure 7b) confirms delocalized states in the spin down component crossing the Fermi level.

3.3. Self-Trapped Exciton in WSe2

We now consider the formation, self-trapping and recombination processes of an exciton in a monolayer of 2D-WSe2. We define the energy quantities involved in those processes in the scheme of Figure 8.
We will first present the results obtained with the hybrid functional HSE with CRY14 code on the 6 × 6 supercell model. Density functional theory is a ground state theory. However, we can bypass this critical issue by computing the first triplet state. This is a ground state for the triplet spin configuration and can, therefore, be correctly described by DFT [32]. Therefore, the vertical excitation will be simulated as a S0→T1 transition ( Δ E a b s = 2.44   eV ). Although, this is forbidden, it highly resembles the S0→S1 transition, since, in both cases, one electron has been promoted to the conduction band, where it is largely delocalized and, therefore, its spin character (up or down) is not expected to influence the energy of the system very much when in the same S0 ground state atomic configuration. We do not expect that this approximated absorption energy reproduces the experimental value. However, we can use this to obtain the exciton self-trapping energy, as the difference between the energy of the T1 state in the S0 atomic configuration and the energy of the T1 state in its minimum atomic configuration:
  Δ E S T , exciton = E T 1 ( S 0 ) E T 1 ( T 1 ) =   0.74   eV
In the ground state atomic configuration, the T1 state is confirmed to be fully delocalized as shown by the spin plot in Figure 9a. Upon geometry optimization, however, the exciton is found to almost fully localize on three W atoms of one hexagon of the system. It is noteworthy to mention that, although we have removed all symmetry constraints, the structure maintains some symmetry, in particular, it keeps a perpendicular plane of symmetry going through a W−Se bond (see dotted line in Figure 9b).
The presence of an exciton leads to two unpaired electrons, one excited in the conduction and one left in the valence band (see the band structure of the T1 states in its minimum energy configuration as reported in Figure 10 with the highest occupied band state in the spin up component in red in the middle of the band gap, 1.25 eV below the bottom of the conduction band, and the lowest unoccupied band state in the spin down component in blue about 0.74 eV above the top of the valence band). Therefore, the total spin density must sum to 2. Indeed, we observe a spin localization of 0.93 on one W atom and of 0.53 on the other two W atoms in one hexagonal ring of the lattice (see Figure 9b). We have made several attempts to remove even this residual symmetry operation; however, the system would always go back to this optimized structure, which we consider the atomic configuration of minimum energy for the T1 state or triplet self-trapped exciton.
The spin density plot describes the overall exciton species. However, one would like to distinguish the contribution from the electron and from the hole. This can be obtained following the procedure proposed in Reference [58] by Van Ginhove et al., where the contribution of e and h+ are separated by subtracting the spin up and spin down density of the S0 solution from the T1 one when they are in the same atomic configuration. The resulting plots for the e and the h+ in the triplet exciton are shown in Figure 11. They actually look quite similar. A fully analogous result is obtained by plotting the highest occupied spin up state and the lowest unoccupied beta state (from the band structure in Figure 10) which are the corresponding one-electron states for the e and the h+ in the triplet exciton.
For the 6 × 6 supercell model, the distance between two repeating excitons is of ~16 Å. We have shown that both the electron and hole are localized on few atoms (three W atoms), which means that the hole is very much attracted by the electron and gets highly localized, in contrast with what was observed for an excess isolated hole in the previous Section 3.2.2. The lattice distortions around the relaxed excitons (Figure 9c) involve both the W−Se bond distance, which are all systematically enlarged, and the W---W or Se---Se or W---Se distances in the hexagon. The hexagon is deformed so that the W---Se distance on the opposite sides becomes 4.38 vs. 4.14 Å in the ground state S0 structure (Figure 5d). While the Se---Se distance is reduced to 3.20 Å, the W---W one is enlarged to 3.42 Å, both with respect to 3.28 Å in the ground state S0 structure.
Using the data discussed above in this and previous sections, we may now provide a theoretical estimation of the exciton binding energy (EB), based on electronic structure calculations and density functional theory, according to the following equation:
  E B = Δ E S T ,   e x c i t o n ( Δ E S T , e + Δ E S T ,   h + )
The exciton binding energy is an extremely delicate quantity, which is very difficult to measure and to theoretically compute. From a theoretical point of view, one should resort to the solution of the Bethe-Salpeter equation, but still with the limitation that this would not allow for atomic relaxation and electron-phonon coupling phenomena. On the contrary, here we compute the self-trapping energy of an extra electron and an extra hole in the system as two separated and independent objects, and then we compute the self-trapping energy of the trapped exciton, where electron and hole interact and are bound. This provides a simple route to the exciton binding energy. We obtain a value of EB = −0.49 eV in excellent agreement with previous experimental and theoretical evaluations in the range between −0.4/−0.8 eV. [9,15,16,17] The large exciton binding energy in this bidimensional material is well explained by the low dielectric screening.
We now suppose that the two charge carriers recombine through a radiative emission. We may compute the photoluminescence energy ( Δ E P L ) as the vertical transition from the T1 state to the S0 state in the T1 minimum energy configuration:
Δ E P L = E S 0 ( T 1 ) E T 1 ( T 1 ) =   1.04   eV
This energy quantity can more easily be compared to an experimental measurement of 1.6–1.7 eV [7,13,23]. Given the approximations used in this single-particle approach, purely based on DFT total energy calculations, with no electron–phonon coupling treatment, the agreement is satisfactory.
We now wish to comment on the use of supercells of different size and then of different functionals with respect to the exciton description. We have performed the same series of calculations for supercells models of different size (N × N) going from a 2 × 2, which was proven to be too small for describing the exciton trapping phenomenon, to 4 × 4, 5 × 5, 6 × 6 and 8 × 8. The results are summarized in Table 2 and indicate that they are rather independent on the supercell size, with only minor variations. Not only are the energy parameters similar but also, the spin density values on the relevant trapping W atoms (in Table 2) and the structural distortions (see Figure S3 in the Supplementary Materials). We tested one model with N = odd number, to verify that a different periodicity would not lead to different exciton features, but that was proven not to be the case.
We will now comment on the opportunity of using a different hybrid functional (B3LYP) and then of using a standard GGA functional (PBE). The hybrid HSE06 functional contains 25% of exact exchange in the exchange functional, and it is a long-range screened functional. On the contrary, B3LYP contains only 20% of exact exchange, and it is not screened. Therefore, it is relevant to compare the results of B3LYP with those of HSE to highlight any difference between the two approaches. We observe that B3LYP tend to localize the electron density on the trapping W atoms (0.94, 0.54, 0.54) slightly more with respect to HSE (0.93, 0.53, 0.53). However, the differences are very tiny. Similar considerations can be made on the structural distortions that are reported for B3LYP in Figure S4 in the Supplementary Materials and on the energetics of absorption, self-trapping and photoluminescence that are reported in Table 3.
If we now look at the results with a standard GGA functional such as PBE, the conclusions are rather different. In this case, we tested both an atomic localized basis set with the CRY14 code and a planewave basis set with the QE code. Results with PBE, independent of the type of basis set are characterized by a lower spin density localization for the trapped exciton with respect to what was observed with the hybrid functionals HSE and B3LYP: (0.70, 0.38, 0.38) with PBE/CRY14 and (0.53, 0.30, 0.27) with PBE/QE. However, one can notice that the structural distortion is rather similar to those with hybrid functionals. Energy changes for the cycle in Figure 8, which are reported in Table 3, are rather underestimated, except for Δ E P L , because of the too small band gap by PBE affecting Δ E A B S , and because of the excessive delocalization reducing the trapping effect ( Δ E S T , e x c i t o n ) and enhancing the Δ E P L . Therefore, we may conclude that even the standard PBE functional, which is affected by a severe electron self-interaction error, may qualitatively correctly describe the exciton in 2D-WSe2. Note that in the case of bulk excitons in TiO2, PBE was found to be totally unsuccessful in the self-trapping process, in contrast with hybrid functionals [32]. We rationalize the unexpected PBE performance on 2D-WSe2 as the consequence of an extremely strong and peculiar exciton binding energy. It is indeed an interesting observation and should be compared with the results for separated charge carriers (electrons and holes) for which any attempt of localization and self-trapping process was totally unsuccessful with PBE. The band structure for the relaxed triplet exciton is shown for comparison with the HSE one (see Figure 12). The defect states are clearly less deep in the band gap than for HSE due to the lower degree of trapping. The plot of those defect states is shown in Figure S5 in the Supplementary Materials, as obtained with the PBE/QE calculations. We also show the corresponding density of states (DOS) for these defects, as obtained with PBE/QE and HSE/CRY14 in Figure S6. Again, we observe that the excitonic levels are deeper in the gap for HSE than for PBE, with the Fermi level above the defect state in the spin up component and the Fermi level below the defect state in the spin down component.

4. Conclusions

To summarize, in this work we have presented a density functional theory investigation of excitons in 2D-WSe2, where we have compared results obtained by standard GGA methods (including spin-orbit coupling) with those by hybrid density functionals. Our study has provided information on the size of the self-trapped exciton, the number and type of atoms involved, the structural reorganization, the self-trapping energy, and the photoluminescence energy, whose computed value (−1.04 eV) is found to be in good agreement with experimental measurements in the literature (1.6–1.7 eV [7,13,23]). In addition, based on the comparative analysis of the self-trapping energy for the exciton with that for isolated charge carriers (unbound electrons and holes), we have also suggested a simplified approach for the theoretical estimation of the excitonic binding energy (−0.49 eV), which can be compared with previous estimates from different approaches or from experimental data (−0.4/−0.8 eV) [9,15,16,17]. To conclude, through this work we have provided the community working with this interesting and useful bidimensional material with an atomistic model of self-trapped exciton in 2D-WSe2, whose reliability and accuracy has been corroborated by comparison of computed relevant quantities, such as photoluminescence and exciton binding energy, with available experimental measurements.

Supplementary Materials

Supplementary Materials are available online at https://www.mdpi.com/2079-4991/8/7/481/s1. Table reporting computed lattice parameters of WSE2 with different density functionals. TDOS and PDOS for the electron doped WSe2 monolayer. Electron charge density plot of the states hosting the extra electron in the electron doped WSe2 monolayer. Ball-and-stick representation reporting atomic distances of the exciton models as calculated with HSE/CRY14 and B3LYP/CRY14. Exciton density of states (DOS). PBE/QE charge density plots of defect states for the self-trapped exciton model.

Author Contributions

Conceptualization, C.D.V.; Data curation, P.L.; Formal analysis, H.L. and P.L.; Investigation, H.L. and P.L.; Methodology, C.D.V.; Project administration, C.D.V.; Resources, C.D.V.; Supervision, C.D.V.; Validation, H.L.; Visualization, H.L.; Writing—original draft, C.D.V.; Writing—review & editing, H.L. and C.D.V.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wilson, A.D.; Yoffe, J.A. The transition metal dichalcogenides discussion and interpretation of the observed optical, electrical and structural properties. Adv. Phys. 1969, 18, 193–335. [Google Scholar] [CrossRef]
  2. Bromley, R.A.; Murray, R.B.; Yoffe, A.D. The band structures of some transition metal dichalcogenides. III. Group VIA: Trigonal prism materials. J. Phys. C Solid Stat. Phys. 1972, 5, 759. [Google Scholar] [CrossRef]
  3. Dickinson, R.G.; Pauling, L. The crystal structure of molybdenite. J. Am. Chem. Soc. 1923, 45, 1466–1471. [Google Scholar] [CrossRef]
  4. Kolobov, A.V.; Tominaga, J. Two-Dimensional Transition Metal Dichalcogenides; Springer: Basel, Switzerland, 2016. [Google Scholar]
  5. Liu, G.B.; Xiao, D.; Yao, Y.; Xu, X.; Yao, W. Electronic structure and theoretical modelling of two-dimensional group-VIB transition metal dichalcogenides. Chem. Soc. Rev. 2015, 44, 2643–2663. [Google Scholar] [CrossRef] [PubMed]
  6. Chhowalla, M.; Shin, H.S.; Eda, G.; Li, L.; Loh, K.P.; Zhang, H. The chemistry of two-dimensional layered transition metal dichalcogenide nanosheets. Nat. Chem. 2013, 5, 263–275. [Google Scholar] [CrossRef] [PubMed]
  7. Yan, T.; Qiao, X.; Liu, X.; Tan, P.; Zhang, X. Photoluminescence properties and excitation dynamics in monolayer WSe2. Appl. Phys. Lett. 2014, 105, 101901. [Google Scholar] [CrossRef]
  8. Chernikov, A.; Berkelbach, T.C.; Hill, H.M.; Rigosi, A.; Li, Y.; Aslan, O.B.; Reichman, D.R.; Hybertsen, M.S.; Heinz, T.F. Exciton binding energy and nonhydrogenic rydberg series in monolayer WS2. Phys. Rev. Lett. 2014, 113, 076802. [Google Scholar] [CrossRef] [PubMed]
  9. Ramasubramaniam, A. Large excitonic effects in monolayers of molybdenum and tungsten dichalcogenides. Phys. Rev. B 2012, 86, 115409. [Google Scholar] [CrossRef]
  10. Choi, W.; Choudhary, N.; Han, G.H.; Park, J.; Akinwande, D.; Lee, Y.H. Recent development of two-dimensional transition metal dichalcogenides and their applications. Mater. Today 2017, 20, 116–130. [Google Scholar] [CrossRef]
  11. Srivastava, A.; Sidler, M.; Allain, A.V.; Lembke, D.S.; Kis, A.; Imamoglu, A. Optically active quantum dots in monolayer WSe2. Nat. Nanotechnol. 2015, 10, 491–496. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Q.H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J.N.; Strano, M.S. Electronics and optoelectronics of two-dimensional transition metal dichalcogenides. Nat. Nanotechnol. 2012, 7, 699–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhao, W.; Zhang, H.; Zhu, G.; Lee, S.; Lin, Z.; Wang, Z.L. Evolution of electronic structure in atomically thin sheets of WS2 and WSe2. ACS Nano 2013, 7, 791–797. [Google Scholar] [CrossRef] [PubMed]
  14. Ross, J.S.; Klement, P.; Jones, A.M.; Ghimire, N.J.; Yan, J.; Mandrus, D.G.; Taniguchi, T.; Watanabe, K.; Kitamura, K.; Yao, W.; et al. Electrically tunable excitonic light-emitting diodes based on monolayer WSe2 p-n junctions. Nat. Nanotechnol. 2014, 9, 268–272. [Google Scholar] [CrossRef] [PubMed]
  15. He, K.; Kumar, N.; Zhao, L.; Wang, Z.; Mak, K.F.; Zhao, H.; Shan, J. Tightly Bound Excitons in Monolayer WSe2. Phys. Rev. Lett. 2014, 113, 026803. [Google Scholar] [CrossRef] [PubMed]
  16. Hanbicki, A.T.; Currie, M.; Kioseoglu, G.; Friedman, A.L.; Jonker, B.T. Measurement of high exciton binding energy in the monolayer transition-metal dichalcogenides WS2 and WSe2. Solid State Commun. 2015, 203, 16–20. [Google Scholar] [CrossRef]
  17. Kylänpää, I.; Komsa, H.P. Binding energies of exciton complexes in transition metal dichalcogenide monolayers and effect of dielectric environment. Phys. Rev. B 2015, 92, 205418. [Google Scholar] [CrossRef]
  18. Hsu, W.; Lu, L.; Wang, D.; Huang, J.; Li, M.; Chang, T.; Chou, Y.; Juang, Z.; Jeng, H.; Li, L.; et al. Evidence of indirect gap in monolayer WSe2. Nat. Commun. 2017, 8, 929. [Google Scholar] [CrossRef] [PubMed]
  19. Zhang, C.; Chen, Y.; Johnson, A.; Li, M.; Li, L.; Mende, P.C.; Feenstra, R.M.; Shih, C. Probing critical point energies of transition metal dichalcogenides: Surprising indirect gap of single layer WSe2. Nano Lett. 2015, 15, 6494–6500. [Google Scholar] [CrossRef] [PubMed]
  20. Ding, Y.; Wang, Y.; Ni, J.; Shi, L.; Shi, S.; Tang, W. First principles study of structural, vibrational and electronic properties of graphene-like MX2 (M = Mo, Nb, W, Ta; X = S, Se, Te) monolayers. Phys. B 2011, 406, 2254–2260. [Google Scholar] [CrossRef]
  21. Rasmussen, F.A.; Thygesen, K.S. Computational 2D Materials Database: Electronic Structure of Transition-Metal Dichalcogenides and Oxides. J. Phys. Chem. C 2015, 119, 13169–13183. [Google Scholar] [CrossRef] [Green Version]
  22. Zhuang, H.L.; Hennig, R.G. Computational Search for Single-Layer Transition-Metal Dichalcogenide Photocatalysts. J. Phys. Chem. C 2013, 117, 20440–20445. [Google Scholar] [CrossRef]
  23. Tonndorf, P.; Schmidt, R.; Böttger, P.; Zhang, X.; Börner, J.; Liebig, A.; Albrecht, M.; Kloc, C.; Gordan, O.; Zahn, D.R.T.; et al. Photoluminescence emission and raman response of monolayer MoS2, MoSe2 and WSe2. Opt. Express 2013, 21, 4908–4916. [Google Scholar] [CrossRef] [PubMed]
  24. Anisimov, V.I.; Zaanen, J.; Andersen, O.K. Band theory and Mott insulators: Hubbard U instead of Stoner I. Phys. Rev. B 1991, 44, 943. [Google Scholar] [CrossRef]
  25. Anisimov, V.I.; Solovyev, I.V.; Korotin, M.A.; Czyżyk, M.T.; Sawatzky, G.A. Density-functional theory and NiO photoemission spectra. Phys. Rev. B 1993, 48, 16929. [Google Scholar] [CrossRef]
  26. Liechtenstein, A.I.; Anisimov, V.I.; Zaanen, J. Density-functional theory and strong interactions: Orbital ordering in Mott-Hubbard insulators. Phys. Rev. B 1995, 52, R5467. [Google Scholar] [CrossRef]
  27. Hybertsen, M.S.; Louie, S.G. Electron Correlation in Semiconductors and Insulators: Band Gaps and Quasiparticle Energies. Phys. Rev. B Condens. Matter Mater. Phys. 1986, 34, 5390–5413. [Google Scholar] [CrossRef]
  28. Runge, E.; Gross, E.K.U. Density-functional theory for time-dependent systems. Phys. Rev. Lett. 1984, 52, 997. [Google Scholar] [CrossRef]
  29. Turkowski, V.; Ud Din, N.; Rahman, T.S. Time-dependent density-functional theory and excitons in bulk and two-dimensional semiconductors. Computation 2017, 5, 39. [Google Scholar] [CrossRef]
  30. Onida, G.; Reining, L.; Rubio, A. Electronic excitations: Density-functional versus many-body Green’s-function approaches. Rev. Mod. Phys. 2002, 74, 601. [Google Scholar] [CrossRef]
  31. Botti, S.; Schindlmayr, A.; Del Sole, R.; Reining, L. Time-dependent density-functional theory for extended systems. Rep. Prog. Phys. 2007, 70, 357. [Google Scholar] [CrossRef]
  32. Di Valentin, C.; Selloni, A. Bulk and Surface Polarons in Photoexcited Anatase TiO2. J. Phys. Chem. Lett. 2011, 2, 2223–2228. [Google Scholar] [CrossRef]
  33. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207. [Google Scholar] [CrossRef]
  34. Dovesi, R.; Orlando, R.; Erba, A.; Zicovich-Wilson, C.M.; Civalleri, B.; Casassa, S.; Maschio, L.; Ferrabone, M.; De La Pierre, M.; D’Arco, P.; et al. CRYSTAL14: A Program for the Ab Initio Investigation of Crystalline Solids. Int. J. Quantum Chem. 2014, 114, 1287–1317. [Google Scholar] [CrossRef]
  35. Dovesi, R.; Saunders, V.R.; Roetti, C.; Orlando, R.; Zicovich-Wilson, C.M.; Pascale, F.; Civalleri, B.; Doll, K.; Harrison, N.M.; Bush, I.J.; et al. CRYSTAL14 User’s Manual; University of Torino: Torino, Italy, 2014. [Google Scholar]
  36. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A Modular and Open-source Software Project for Quantum Simulations of Materials. J. Phys. Condens. Matter 2009, 21, 395502. [Google Scholar] [CrossRef] [PubMed]
  37. Hay, P.J.; Wadt, W.R. Ab Initio Effective Core Potentials for Molecular Calculations. Potentials for K to Au including the Outermost Core Orbitals. J. Chem. Phys. 1985, 82, 299. [Google Scholar] [CrossRef]
  38. Durand, P.H.; Barthelat, J.C. A Theoretical Method to Determine Atomic Pseudopotentials for Electronic Structure Calculations of Molecules and Solids. Theor. Chim. Acta 1975, 38, 283. [Google Scholar] [CrossRef]
  39. Wang, F.; Di Valentin, C.; Pacchioni, G. Electronic and Structural Properties of WO3: A Systematic Hybrid DFT Study. J. Phys. Chem. C 2011, 115, 8345–8353. [Google Scholar] [CrossRef]
  40. Heyd, J.; Peralta, J.E.; Scuseriaa, G.E. Energy Band Gaps and Lattice Parameters Evaluated with the Heyd-Scuseria-Ernzerh of Screened Hybrid Functional. J. Chem. Phys. 2005, 123, 174101. [Google Scholar] [CrossRef] [PubMed]
  41. Pernot, P.; Civalleri, B.; Presti, D.; Savin, A. Prediction Uncertainty of Density Functional Approximations for Properties of Crystals with Cubic Symmetry. J. Phys. Chem. A 2015, 119, 5288–5304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Krukau, A.V.; Vydrov, O.A.; Izmaylov, A.F.; Scuseria, G.E. Influence of the Exchange Screening Parameter on the Performance of Screened Hybrid Functionals. J. Chem. Phys. 2006, 125, 224106. [Google Scholar] [CrossRef] [PubMed]
  43. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758. [Google Scholar] [CrossRef]
  44. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  45. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 8, 5648–5652. [Google Scholar] [CrossRef]
  46. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti Correlation Energy Formula into a Functional of the Electron-Density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  47. Becke, A.D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic-Behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  48. Tanabe, I.; Gomez, M.; Coley, W.C.; Le, D.; Echeverria, E.M.; Stecklein, G.; Kandyba, V.; Balijepalli, S.K.; Klee, V.; Nguyen, A.E.; et al. Band structure characterization of WS2 grown by chemical vapor deposition. Appl. Phys. Lett. 2016, 108, 252103. [Google Scholar] [CrossRef] [Green Version]
  49. Kozawa, D.; Kumar, R.; Carvalho, A.; Amara, K.K.; Zhao, W.; Wang, S.; Toh, M.; Ribeiro, R.M.; Castro Neto, A.H.; Matsuda, K.; et al. Photocarrier relaxation pathway in two-dimensional semiconducting transition metal dichalcogenides. Nat. Commun. 2014, 5, 4543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Baugher, B.W.H.; Churchill, H.O.H.; Yang, Y.; Jarillo-Herrero, P. Optoelectronic devices based on electrically tunable p–n diodes in a monolayer dichalcogenide. Nat. Nanotechnol. 2014, 9, 262. [Google Scholar] [CrossRef] [PubMed]
  51. Klots, A.R.; Newaz, A.K.M.; Wang, B.; Prasai, D.; Krzyzanowska, H.; Lin, J.; Caudel, D.; Ghimire, N.J.; Yan, J.; Ivanov, B.L.; et al. Probing excitonic states in suspended two-dimensional semiconductors by photocurrent spectroscopy. Sci. Rep. 2014, 4, 6608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Kosmider, K.; Fernandez-Rossier, J. Electronic properties of the MoS2-WS2 heterojunction. Phys. Rev. B 2013, 87, 075451. [Google Scholar] [CrossRef]
  53. Kormányos, A.; Burkard, G.; Gmitra, M.; Fabian, J.; Zolyomi, V.; Drummond, N.D.; Fal’ko, V. k•p theory for two-dimensional transition metal dichalcogenide semiconductors. 2D Mater. 2015, 2, 022001. [Google Scholar] [CrossRef]
  54. Zhu, Z.Y.; Cheng, Y.C.; Schwingenschlogl, U. Giant spin-orbit-induced spin splitting in two-dimensional transition-metal dichalcogenide semiconductors. Phys. Rev. B 2011, 84, 153402. [Google Scholar] [CrossRef]
  55. Wickramaratne, D.; Zahid, F.; Lake, R.K. Electronic and thermoelectric properties of few-layer transition metal dichalcogenides. J. Chem. Phys. 2014, 140, 124710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Yuan, H.; Bahramy, M.S.; Morimoto, K.; Wu, S.; Nomura, K.; Yang, B.J.; Shimotani, H.; Suzuki, R.; Toh, M.; Kloc, C.; et al. Zeeman-type spin splitting controlled by an electric field. Nat. Phys. 2013, 9, 563–569. [Google Scholar] [CrossRef]
  57. Kang, J.; Tongay, S.; Zhou, J.; Li, J.; Wu, J. Band offsets and heterostructures of two-dimensional semiconductors. Appl. Phys. Lett. 2013, 102, 012111. [Google Scholar] [CrossRef]
  58. Van Ginhoven, R.M.; Jónsson, H.; René Corrales, L. Characterization of Exciton Self-trapping in Amorphous Silica. J. Non-Cryst. Solids 2006, 352, 2589–2595. [Google Scholar] [CrossRef]
Figure 1. (a) Top and side views of a repeated unit cell of 2D tungsten diselenide (2D-WSe2); (b) A view of the trigonal prismatic coordination and a scheme of the d states splitting in WSe2 monolayer deriving from it.
Figure 1. (a) Top and side views of a repeated unit cell of 2D tungsten diselenide (2D-WSe2); (b) A view of the trigonal prismatic coordination and a scheme of the d states splitting in WSe2 monolayer deriving from it.
Nanomaterials 08 00481 g001
Figure 2. Band structures of monolayer tungsten diselenide (WSe2) calculated by (a) HSE; (b) PBE; (c) B3LYP and (d) B3PW with CRYSTAL14 code. The Fermi level is scaled to zero (the top of valence band).
Figure 2. Band structures of monolayer tungsten diselenide (WSe2) calculated by (a) HSE; (b) PBE; (c) B3LYP and (d) B3PW with CRYSTAL14 code. The Fermi level is scaled to zero (the top of valence band).
Nanomaterials 08 00481 g002
Figure 3. Band structure of monolayer WSe2 unit cell obtained by Perdew-Burke-Ernzerhof/Quantum ESPRESSO (PBE/QE) with (a) and without (b) spin–orbit coupling. The Fermi level is scaled to zero (the top of valence band).
Figure 3. Band structure of monolayer WSe2 unit cell obtained by Perdew-Burke-Ernzerhof/Quantum ESPRESSO (PBE/QE) with (a) and without (b) spin–orbit coupling. The Fermi level is scaled to zero (the top of valence band).
Nanomaterials 08 00481 g003
Figure 4. Total and projected density of states (TDOS and PDOS) of monolayer WSe2 calculated by hybrid density functional calculations/CRYSTAL14 (HSE06/CRY14). The grey area represents the TDOS. In (a), blue and red lines represent PDOS on 4p and 5d states of Se and W; In (b), black, blue, green and red lines represent PDOS on 5d, d z 2 , dxz + dyz and dxy +   d x 2 y 2 , respectively. The Fermi level is scaled to zero (the top of valence band).
Figure 4. Total and projected density of states (TDOS and PDOS) of monolayer WSe2 calculated by hybrid density functional calculations/CRYSTAL14 (HSE06/CRY14). The grey area represents the TDOS. In (a), blue and red lines represent PDOS on 4p and 5d states of Se and W; In (b), black, blue, green and red lines represent PDOS on 5d, d z 2 , dxz + dyz and dxy +   d x 2 y 2 , respectively. The Fermi level is scaled to zero (the top of valence band).
Nanomaterials 08 00481 g004
Figure 5. (a,d) are spin density and geometry of WSe2 6 × 6 supercell with one more electron before geometry relaxation by HSE/CRY14; (b,c) are spin densities of two solutions for WSe2 6 × 6 supercell with a self-trapped electron; (e,f) are the local geometries of solution (b) and (c), respectively. Distances are in Å. The isosurface level is 0.006 electron/bohr3 for the spin density plots. Yellow and cyan are positive and negative spin densities, respectively. Mulliken electron spin density values on W atoms are given in (b) and (c).
Figure 5. (a,d) are spin density and geometry of WSe2 6 × 6 supercell with one more electron before geometry relaxation by HSE/CRY14; (b,c) are spin densities of two solutions for WSe2 6 × 6 supercell with a self-trapped electron; (e,f) are the local geometries of solution (b) and (c), respectively. Distances are in Å. The isosurface level is 0.006 electron/bohr3 for the spin density plots. Yellow and cyan are positive and negative spin densities, respectively. Mulliken electron spin density values on W atoms are given in (b) and (c).
Nanomaterials 08 00481 g005
Figure 6. Band structure of the 6 × 6 supercell WSe2 monolayer with one self-trapped electron obtained by HSE06/CRY14. (a) is for the most stable solution (in Figure 5b) and (b) is for the less stable solution (in Figure 5c). Red and blue lines represent the spin up (majority) and spin down (minority) bands, respectively. The Fermi level is scaled to zero. One energy difference at Γ is reported in eV.
Figure 6. Band structure of the 6 × 6 supercell WSe2 monolayer with one self-trapped electron obtained by HSE06/CRY14. (a) is for the most stable solution (in Figure 5b) and (b) is for the less stable solution (in Figure 5c). Red and blue lines represent the spin up (majority) and spin down (minority) bands, respectively. The Fermi level is scaled to zero. One energy difference at Γ is reported in eV.
Nanomaterials 08 00481 g006
Figure 7. (a) Spin density and band structure of the 6 × 6 supercell WSe2 monolayer with one hole obtained by HSE06/CRY14. In (a), the isosurface level is 0.006 electron/bohr3; In (b), red and blue lines represent the spin up (majority) and spin down (minority) bands. The Fermi level is scaled to zero.
Figure 7. (a) Spin density and band structure of the 6 × 6 supercell WSe2 monolayer with one hole obtained by HSE06/CRY14. In (a), the isosurface level is 0.006 electron/bohr3; In (b), red and blue lines represent the spin up (majority) and spin down (minority) bands. The Fermi level is scaled to zero.
Nanomaterials 08 00481 g007
Figure 8. Schematic representation of the S0-T1 excitation (ΔEABS), exciton self-trapping (ΔEST) and emission or photoluminescence (ΔEPL) energies for 2D-WSe2.
Figure 8. Schematic representation of the S0-T1 excitation (ΔEABS), exciton self-trapping (ΔEST) and emission or photoluminescence (ΔEPL) energies for 2D-WSe2.
Nanomaterials 08 00481 g008
Figure 9. Spin density plots for the triplet state (a) before and (b) after geometry relaxation calculated by HSE06/CRY14 with a 6 × 6 supercell of WSe2; (c) Geometry parameters (in Å) after the relaxation due to the exciton. The isosurface level is 0.0015 and 0.006 electron/bohr3 for (a) and (b), respectively. Mulliken electron spin density values on W atoms are given in (b).
Figure 9. Spin density plots for the triplet state (a) before and (b) after geometry relaxation calculated by HSE06/CRY14 with a 6 × 6 supercell of WSe2; (c) Geometry parameters (in Å) after the relaxation due to the exciton. The isosurface level is 0.0015 and 0.006 electron/bohr3 for (a) and (b), respectively. Mulliken electron spin density values on W atoms are given in (b).
Nanomaterials 08 00481 g009
Figure 10. Band structure for the triplet state calculated by HSE06/CRY14 with a 6 × 6 supercell of WSe2. The red and blue lines represent the spin up and spin down bands, respectively. The black dashed line indicates the Fermi level and the Fermi level for spin up states is scaled to zero. Some energy difference is reported in eV.
Figure 10. Band structure for the triplet state calculated by HSE06/CRY14 with a 6 × 6 supercell of WSe2. The red and blue lines represent the spin up and spin down bands, respectively. The black dashed line indicates the Fermi level and the Fermi level for spin up states is scaled to zero. Some energy difference is reported in eV.
Nanomaterials 08 00481 g010
Figure 11. Charge plots for the e (a) and the h+ (b) in the triplet exciton by HSE/CRY14. The isosurface level is 0.006 electron/bohr3.
Figure 11. Charge plots for the e (a) and the h+ (b) in the triplet exciton by HSE/CRY14. The isosurface level is 0.006 electron/bohr3.
Nanomaterials 08 00481 g011
Figure 12. The band structures for the relaxed triplet exciton from PBE calculation by CRY14 (a) and QE (b). The red and blue lines represent the spin up and spin down bands, respectively. The black dashed line indicates the Fermi level and the Fermi level for spin up states is scaled to zero. Some energy difference is reported in eV.
Figure 12. The band structures for the relaxed triplet exciton from PBE calculation by CRY14 (a) and QE (b). The red and blue lines represent the spin up and spin down bands, respectively. The black dashed line indicates the Fermi level and the Fermi level for spin up states is scaled to zero. Some energy difference is reported in eV.
Nanomaterials 08 00481 g012
Table 1. Contribution of atomic orbitals to the valence band maximum (VBM) and to the conduction band minimum (CBM) states projected on different k-points calculated with HSE06/CRY14.
Table 1. Contribution of atomic orbitals to the valence band maximum (VBM) and to the conduction band minimum (CBM) states projected on different k-points calculated with HSE06/CRY14.
K PointMain ContributionMinor Contribution
KCW- d z 2 Se-px, py
KVW- d x 2 y 2 ,   d x y Se-px, py
QCW- d x 2 y 2 ,   d x y W- d z 2 ; Se-px, py, pz
ΓVW- d z 2 Se-pz
Table 2. Comparison of the energy parameters and spin density for the intermediates involved in the processes schematically described in Figure 8 for different supercells using the hybrid density functional calculations (HSE06/CRY14). The definition of the energy parameters is given in Section 2 of Computational details.
Table 2. Comparison of the energy parameters and spin density for the intermediates involved in the processes schematically described in Figure 8 for different supercells using the hybrid density functional calculations (HSE06/CRY14). The definition of the energy parameters is given in Section 2 of Computational details.
SupercellEnergy Parameters (eV)Spin Density
ΔEABSΔESTΔEPLΔEDIST
4 × 42.47−0.75−1.040.680.922
0.534
0.534
5 × 52.38−0.66−1.040.680.930
0.530
0.530
6 × 62.44−0.74−1.040.660.934
0.533
0.533
8 × 82.35−0.69−1.030.630.935
0.535
0.535
Table 3. Comparison of the energy parameters for the intermediates involved in the processes schematically described in Figure 8 as obtained with different functionals and a 6 × 6 supercell. The definition of the energy parameters is given in Section 2 of Computational details.
Table 3. Comparison of the energy parameters for the intermediates involved in the processes schematically described in Figure 8 as obtained with different functionals and a 6 × 6 supercell. The definition of the energy parameters is given in Section 2 of Computational details.
FunctionalEnergy Parameters (eV)
ΔEABSΔESTΔEPLΔEDIST
HSE06/CRY142.44−0.74−1.040.66
B3LYP/CRY142.38−0.81−0.900.67
PBE/CRY141.91−0.17−1.170.57
PBE/QE1.71−0.19−1.300.22

Share and Cite

MDPI and ACS Style

Liu, H.; Lazzaroni, P.; Di Valentin, C. Nature of Excitons in Bidimensional WSe2 by Hybrid Density Functional Theory Calculations. Nanomaterials 2018, 8, 481. https://doi.org/10.3390/nano8070481

AMA Style

Liu H, Lazzaroni P, Di Valentin C. Nature of Excitons in Bidimensional WSe2 by Hybrid Density Functional Theory Calculations. Nanomaterials. 2018; 8(7):481. https://doi.org/10.3390/nano8070481

Chicago/Turabian Style

Liu, Hongsheng, Paolo Lazzaroni, and Cristiana Di Valentin. 2018. "Nature of Excitons in Bidimensional WSe2 by Hybrid Density Functional Theory Calculations" Nanomaterials 8, no. 7: 481. https://doi.org/10.3390/nano8070481

APA Style

Liu, H., Lazzaroni, P., & Di Valentin, C. (2018). Nature of Excitons in Bidimensional WSe2 by Hybrid Density Functional Theory Calculations. Nanomaterials, 8(7), 481. https://doi.org/10.3390/nano8070481

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop