Abstract

The structural, electronic, vibrational, and optical properties of 2D- 2H-WTe2 monolayer are investigated using density functional theory with respect to a plane wave ultrasoft pseudopotentials (PW-USPPs) in a generalized gradient approximation (GGA) and with the Hubbard potential (GGA + U). The equilibrium state properties such as lattice parameters, unit cell volume, bulk modulus, and its derivative are determined. The band gap values of monolayer 2H-WTe2 are investigated for unstrained, 2% biaxial compression, and biaxial tensile stress using GGA, respectively. The obtained band gap values of 2H-WTe2 with respect to GGA are 1.043, 1.1487, and 0.9439 eV for unstrained, biaxial compression, and tensile strain, respectively. Moreover, the band gap values determined using Hubbard correction (GGA + U) are 1.1089 eV (unstrained), 1.2332 eV (2% biaxial compression), and 0.9945 eV (2% biaxial tensile stress), respectively. The band gap value obtained using Hubbard correction predicts the experimental value more precisely. The projected density of state shows W (3d) orbital is more dominant both in the valence band maximum and conduction band minimum. Moreover, a small amount of tensile or compressive strain is used to tune the band gap of the monolayer without affecting its direct band gap nature. In addition to this, the phonon dispersion and optical properties are discussed for tensile strain, unstrained, and compressive strain, respectively.

1. Introduction

Due to their tunable direct band gaps, monolayer transition metal dichalcogenides (TMDs) are a preferred family of two-dimensional (2D) materials. The range of band gap values of 2D transition metal dichalcogenides made them an interesting material for technological applications in optoelectronic and photovoltaic devices [1, 2]. Due to their significant chemical and physical features molybdenum and tungsten dichalcogenides (MoS2, MoSe2, WS2, and WSe2) with a band gap of 1.88, 1.57, 2.03, and 1.67, respectively, have received high attention as a 2D layered material [3]. These band gap values fall within the band gap value of ideal photovoltaic materials (1–1.8 eV) [4]. Moreover, 2D layered structures enable doping with the additional atoms in substitution of one or more atoms to adjust and tune its electronic and optical properties [5, 6]

The hexagonal tungsten ditelluride is one among the family of materials composed of WX2 (W-tungsten and chalcogens X = S, Se, Te), which have been made from a 2D slab of layers with strong covalent bonds which are stacked together with the weak van-der Waals force with the next layer. The 2D-WTe2 has the highest mobility among this family because of its lowest effective mass [7]. On the other hand, strain is vigorous in enhancing device performance because it alters the band structure, electron mobility, and other important properties of 2D-TMDs in general and it is also essential in the case of 2D-WTe2 [8, 9]. One of the most obvious ways to acquire a broad range of tunability in electronic and optical properties is to manipulate the band gap of two-dimensional materials by strain in addition to doping and other techniques. Strain engineering is one of the powerful methods for tuning the band gap and significant performance enhancement of electronic and optoelectronic devices [10].

DFT simulations were used to determine the optical absorption characteristics of WS2, and these characteristics make WS2 a promising photovoltaic material [11]. MoS2, MoSe2, WS2, and WSe2 monolayers complex dielectric functions for photon energies between 1.5 and 3.0 eV have been experimentally calculated. These four materials have a peak absorbance of more than 15% and the dielectric functions predict extensive light-matter interactions in the monolayers. Exciton splitting is extrapolated from the measured dielectric function in conformity with density functional models [12]. For 2D-MoSe2 and 2D-WSe2, the biaxial strain has been used to identify the direct-to-indirect band gap transition in the compressive strain. It has also been reported that the peaks for real and imaginary values of the dielectric constant, optical absorption, and electron energy loss spectra, may allow for appealing applications for photovoltaic and electroluminescent devices. They were strong prospects for digital electronics and optoelectronics device applications due to their optical performance and light absorption under tensile strain [13]. However, most of the studies on WTe2 did not consider the consequences of the residual self-interaction. One of the corrective schemes widely used to alleviate the challenges of residual self-interaction is the GGA + U approach [14, 15]. Moreover, the density functional perturbation theory (DFPT) with Hubbard correction (DFPT + U functional) [16] gives more precise results for phonon dispersion as compared to the DFPT [17].

In this study, we have investigated the structural, electronic, vibrational, and optical properties of 2H-WTe2 (space group p63/mmc) for photovoltaic applications and digital electronics using density functional theory (DFT) and time-dependent density functional theory (TD-DFT) with the help of a plane wave ultrasoft pseudopotentials (PW-USPPs). The structural, electronic, vibrational, and optical properties are investigated in detail. The structural and electronic properties are investigated using the approximations GGA [18] and GGA + U [14, 15] for exchange-correlation potential. Moreover, the effect of tensile or compressive strain on the electronic structure of WTe2 is determined. The phonon dispersion is calculated using DFPT and DFPT + U. The optical properties are determined using TD-DFPT. Finally, the obtained results are compared with the existing experimental and theoretical results.

2. Computational Details

In this study we have used DFT with the implementation of Quantum Espresso (QE) package [19, 20], applying the generalized gradient approximation (GGA) with the Perdew–Burke–Ernzerhof (PBE) functional and including the Hubbard correction GGA + U for the electron exchange-correlation interactions. Here we have used GGA + U since our material contains a transition heavy metal W which has high-electron–electron correlations. The pseudopotential developed for the QE package (PBE QE UPF) from the standard solid-state pseudopotential (SSSP) library is used [21]. The corresponding valence electrons considered for the calculations are W-[Xe] 4 f14 5d4 6 s2, and Te-[Kr] 4d10 5 s2 5 p4. For structural and electronic parameter calculations, a 2 × 2 × 1 supercell of 2D 2H-WTe2 unit cell was designed. For phonon dispersion investigations, a 4 × 4 × 1 supercell of 2D 2H-WTe2 unit cell was developed. The energy cutoff is attained at 55 Ry (750 eV) with a tolerance of less than  eV and the Brillion zone with the Monkhorst–Pack arrangement for k-point grids is considered for the monolayer. The equilibrium atomic structure is optimized by relaxing lattice atomic coordinates with a tolerance of less than for each atom in the structure and these convergence tests are acceptable according to the standards for the next calculations. Under these conditions, the lattice parameters and the vacuum space, normal to the plane of the layer are calculated and is used for this work. Tensile and compressive straining circumstances are taken into consideration in two dimensions parallel to the layer plane following optimization of the equilibrium lattice parameters. The structural, electrical, lattice vibration, and optical properties are computed under both stressed and unstrained circumstances. Biaxial strain is applied to the monolayer to the desired properties from 2% tension to 2% compression ratios to analyze the electronic property, lattice vibration, and dielectric function.

The effective Hubbard parameter (Ueff) was determined iteratively for W-d orbitals using a linear response formalism utilizing DFPT with an ortho-atomic projection method [22, 23]. The Hubbard parameters are calculated independently for three cases using theQE package and are described in Table 1.

3. Result and Discussion

3.1. Crystal Structure

The crystal structure of 2D 2H-WTe2 is described in Figure 1. In monolayer of 2D 2H-WTe2, c is a constant which determines the inter layer distance and the monolayer is characterized by in plane lattice constant After relaxing the structure in a unit cell to minimum total energy, the lattice constant for a unit cell of 2D 2H-WTe2 is found to be 3.565 Å (GGA) and 3.539 Å (GGA + U). The calculated lattice constants are compared and are in good agreement with the previous theoretical results as shown in Table 2.

We have seen the structural properties computing important parameters to check the structural stability of WTe2 monolayer. The computed structural parameters for the layer are in a good agreement with the Birch Murnaghan equation of state (EOS) [24] given bywhere is for minimum energy, is volume of the unit cell, is bulk modulus and is derivative of the bulk modulus.

The calculated volume versus energy and volume versus pressure diagrams are indicated for 2H-WTe2 monolayer of supercell is described in Figure 2. Moreover, the obtained results are explained in Table 2.

3.2. Electronic Properties: Band Structures and Density of States (DOS)

Different methods have been applied when analyzing the electronic properties of the system including the Hubbard model. The band structure of monolayer WTe2 is shown in Figure 3 (a)–3(c) and it has a direct band gap. When the system is imposed to small strain (0%–2%) both in tension and compression, the band gap increases with biaxial compression and decreases with biaxial tensile strain.

The band structure of monolayer 2D 2H-WTe2 using GGA and GGA + U approximations for exchange correlation potentials has been indicated in Figure 3 with the black (solid) line and magenta (broken) line, respectively. The band gap values obtained for 2% biaxial compression of the system are 1.1487 eV (GGA) and 1.2332 eV (GGA + U), respectively, as shown in Figure 3(a). For unstrained condition, the band gap values are 1.0439 eV (GGA) and 1.1089 eV (GGA + U) as described in Figure 3(b). In case of 2% biaxial tensile strain the band gap values are 0.9439 eV (GGA) and 0.9945 eV (GGA + U) as demonstrated in Figure 3(c). Moreover, the obtained band gap values are compared with the previous theoretical results and shown in Table 3. We determined some basic and important properties of the 2D 2H-WTe2 from these computations. The Hubbard on site correction improves the band gap value of the monolayer as expected. Moreover, the band gap value is amplified when the compressive biaxial strain is applied to the monolayer. However, the band gap of the monolayer decreases in the presence of tensile biaxial strain. This result is in contrast to the trend observed in black SiSe monolayer in which the band gap increases with the tensile biaxial strain and reduces in the presence of compressive biaxial strain as investigated in [27]. The direct band gap cannot be altered by small strain and this may result in materials preferable for some electronic and optical properties. Moreover, the effect of electron–electron interaction in the band diagram is noticeable with the same manner in both strained and unstrained conditions and hence the Hubbard correction is significant for the system. The range of band gap values of monolayer 2D 2H-WTe2 both in unstrained and strained conditions are interesting for the photovoltaic applications.

In addition to the band structure, the total density of states (TDOS) was calculated for 2% biaxial compression, unstrained, and 2% biaxial tensile strain using GGA and GGA + U approximation for exchange correlation potentials as shown in Figure 4(a)4(c), respectively. The gap from the valence band maximum to the conduction band minimum decreases when the system is imposed to biaxial tensile strain than the biaxial compression. Moreover, the partial density of states (PDOS) was also calculated using GGA and GGA + U approximation for exchange correlation potential and described in Figures 5(a) and 5(b). The PDOS result describes that for the low-lying states, the W (3d) orbitals contributed most to the maximum valence band and minimum conduction band.

3.3. Phonon Vibration

Concerning the lattice vibration, it is important to know the dynamic property of the material for its stability. We have calculated the phonons dispersion for 2% biaxial compression, unstrained, and 2% biaxial tensile strain in accord with DFPT and DFPT + U for exchange-correlation potential, respectively, as shown in Figure 6(a)6(c). The obtained result for DFPT shows that the unstrained condition is dynamically stable and has no imaginary part as shown in Figure 6(b). However, when the monolayer of 2H-WTe2 is subjected to biaxial strain the dispersion relation has negative frequency and it is dynamically unstable with respect to DFPT as described in Figures 6(a) and 6(c), respectively. However, the phonon dispersions obtained for biaxial compression, unstrained, and biaxial tensile strain using the Hubbard correction (DFPT + U) are dynamically stable as described in Figure 6 (magenta, broken lines). That is, DFPT + U is used to capture the effect of electronic localization of vibrational frequencies and give full access to the all quantities requiring well-converged sums over the entire vibration spectrum.

3.4. Optical Properties

To calculate the optical properties, the complex dielectric function ε(ω) within an independent particle formalism is used [28];

The imaginary part is determined from the momentum dipole transition matrix elements between the valence band and the conduction band electronic states with long wavelength approximations [29],where and are the electron states. The dispersion of the real part of the dielectric function was determined using the imaginary part by Kramers–Kronig relations [30],

The symbol P represents the principal value of the integral. The complex index of refraction is written as follows:

Here, is refractive index and is extinction coefficient that can be evaluated from complex dielectric function. The refractive index (), the reflectivity (), absorption coefficient (), and energy loss function () can be determined from the dielectric functions by,

The optical properties including dielectric function, optical absorption, and energy loss function, reflectivity, and the refractive index are investigated in parallel and perpendicular direction to the plane of the monolayer 2H-WTe2 in the photon energy range from 0 to 10 eV, and the details of the calculations are presented.

The real part of the dielectric function () when the monolayer of 2H-WTe2 is exposed to the polarization parallel and perpendicular direction to the monolayer plane is shown in Figures 7(a) and 7(b), respectively. From the figure it is possible to deduce that the static dielectric function values for is 10.6 (unstrained), 7.6 (biaxial strain), and 7.3 (biaxial compression) arbitrary unit for the polarization parallel to the plane of monolayer. However, the static dielectric constant decreases when the polarization is perpendicular to the plane of monolayer and the values are 4.38 (unstrained) and 1.9 for biaxial compression or strain. The magnitude of the real part of the dielectric function for the polarization parallel to the plane is magnified as compared to the perpendicular direction. The peak of the real part of the dielectric constant is 14.9 at 0.9 eV for biaxial tensile strain, 17.7 at 1 eV for unstrained, and 14.9 at 1.1 eV for biaxial compression (parallel to the plane) and 4.2 at 4.1 eV for biaxial tensile strain, 12.5 at 4.5 eV for unstrained, and 4.2 at 4.8 eV for biaxial compression (perpendicular to the plane). Moreover, the peak of linear dielectric function for the polarization perpendicular to the plane of monolayer shifts toward high energy for all cases. The biaxial compression and tensile stress also reduce the magnitude of the real part of the dielectric function as compared to the unstrained one. On the other hand, the imaginary part of the dielectric function is described in Figure 7(c) for the polarization parallel to the plane and Figure 7(d) perpendicular to the plane with respect to the incoming photon. The imaginary part of the dielectric function has a peak value 16.25 at 0.9 eV for biaxial tensile strain, 18.5 at 3.9 eV for unstrained, and 16.3 at 1.1 eV for biaxial compression, respectively, when the polarization is parallel to the plane of monolayer. The peak value is 7.5 at 4.1 eV for biaxial tensile strain, 17.5 at 4.5 eV for unstrained, and 7.1 at 4.8 eV for biaxial compression when the polarization is perpendicular to the plane of monolayer. The magnitude of the imaginary part of the dielectric function is quenched as the monolayer of 2H-WTe2 forced to tensile stress or compression. The negative value in the spectra at 4.5, 4.8, and 5 eV corresponds to the peak in the reflectivity spectra.

The other important parameters of optical properties are the refractive index and extinction coefficient as shown in Figures 8(a) and 8(b) for polarization perpendicular to the plane of monolayer, respectively. The value of static refractive index is 2.2 (unstrained) and 1.35 for biaxial compression or strain. The peak values of the refractive index are 2.5 at 4.25 eV for tensile strain, 3.85 at 4.48 eV for unstrained, and 2.5 at 4.6 eV for biaxial compression. The magnitude of the refractive index decreases when the monolayer of 2H-WTe2 is imposed to tensile stress or compression. Moreover, the peak of the refractive index shifts toward higher energy when the monolayer of 2H-WTe2 is imposed to biaxial compression than tensile strain. The peak value of the extinction coefficient is 1.9 at 4.3 eV for tensile stress, 2.75 at 4.5 eV for unstrained, and 2.2 at 5 eV for biaxial compression. The peak value of the extinction coefficient for monolayer of WTe2 shows blue shift as compared to the tensile stress.

The absorption coefficient of monolayer 2H-WTe2 for unstrained, biaxial compression, and biaxial tensile strain is also described in Figure 9(a)–9(c), respectively. The peaks of the absorption coefficients are at 8.0636 eV (2% biaxial compression), at 7.54 eV (unstrained), and at 7.15 eV (2% biaxial tension). The magnitude of the absorption coefficient of the monolayer of 2H-WTe2 has a little variation in three cases. However, red shift is observed when the monolayer is imposed to biaxial tension as compared to the biaxial compression. This result is a consequence of the variation of the band gap value from 1.2332 to 0.9945 eV as the monolayer is exposed to a biaxial compression and biaxial tensile stress. The large value of absorption coefficient in 2H-WTe2 monolayer coupled with its optimum band gap enables it for application of solar cell with high-power conversion efficiency.

The electron energy loss function of monolayer 2H-WTe2 for unstrained, biaxial compression, and biaxial tensile strain is described in Figure 9(d)–9(f), respectively. The energy loss function (EELS) which shows the energy loss of fast electron moving across the medium is a complicated mixture of interband transition and plasmon. The electron energy loss spectra are related to the plasma frequency. Peaks in the measured spectra can thus derive zeros of (ω) (when (ω) is not too large), which correspond to the collective plasmon excitations, or from structures (ω) which are due to valence-conduction interband transitions. The electron energy loss function for a polarization perpendicular to the plane of monolayer has a peak 9.1 at 8.5 eV (unstrained), 9.6 at 8.1 eV (biaxial compression), and 9.2 at 7.5 eV (biaxial strain).

Finally, the reflectivity is computed for unstrained (a), biaxial compression (b), and biaxial stress (c) as shown in Figure 10 for polarization perpendicular to the plane of monolayer 2H-WTe2. The value of static reflectivity are 12.5%, 4.2%, and 12.5% for biaxial compression, unstrained, and biaxial, respectively. The peak values of optical reflectivity are 52% at 4.8 eV (biaxial compression), 65% at 5 eV (unstrained), and 52% at 4.9 eV (biaxial stress).

4. Conclusions

In this study we have used the first principle calculations implemented in QE software package to study the structural, electronic, dynamic, and optical properties of 2D 2H-WTe2 monolayer. The lattice constant, bulk modulus, and the derivative of bulk modulus are calculated using GGA and GGA + U approximations for exchange correlation potentials in PBE functional. The obtained results are in good agreement with the previous theoretical results. The band gap values calculated with respect to GGA are 1.0439, 1.1487, and 0.9439 eV for unstrained, 2% biaxial compression, and 2% biaxial tensile strain, respectively. However, the band gap values determined using GGA + U approximations for exchange correlation potential are 1.1089, 1.2332, and 0,9945 eV for unstrained, biaxial compression, and biaxial tensile strain, respectively. The onsite Hubbard corrections are used to remove self-interaction in the band gap of the monolayer 2H-WTe2. Moreover, the discontinuity of TDOS from the valence band maximum to conduction band minimum reproduces the obtained band gap values. In addition to this imposing the monolayer to a small amount biaxial compression or tensile stress tunes the band gap value without altering the direct band gap nature of the system. The phonon dispersion is calculated for 2% biaxial compression, unstrained, and 2% biaxial tensile stress using DFPT and DFPT + U. The DFPT result shows that the phonon dispersion calculated for unstrained condition is stable. In contrast to this the phonon dispersion obtained employing DFPT for the biaxial compression and tensile stress are not dynamically stable. On the other hand, the phonon dispersions determined using DFPT + U are dynamically stable for unstrained and strained conditions. This is due to the consequence that DFPT + U is used to capture the effect of electronic localization of vibrational frequencies over the entire vibrational spectrum. Finally, the optical properties of 2D 2H-WTe2 monolayer such as real and imaginary part of the dielectric function, refractive index, absorption coefficient, electron energy loss function, and the reflectivity are investigated. The high-absorption power in the visible spectral region complemented with the direct band gap make 2D 2H-WTe2 the best candidate for photovoltaic devices. After all the experimental verification of the electronic, phonon dispersion, and optical properties of the monolayer 2H-WTe2 are recommended in the future works.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors express their thanks and appreciation to the Department of Physics, College of Natural Sciences, Jimma University (grant number: CNSPHYS-04-2021/2022), for financially supporting the Ph.D. student.