Next Article in Journal
An Inhalable Powder Formulation Based on Micro- and Nanoparticles Containing 5-Fluorouracil for the Treatment of Metastatic Melanoma
Next Article in Special Issue
Synthesis of Large-Scale Single-Crystalline Monolayer WS2 Using a Semi-Sealed Method
Previous Article in Journal
Preparation of Cuprous Oxide Mesoporous Spheres with Different Pore Sizes for Non-Enzymatic Glucose Detection
Previous Article in Special Issue
Exploring Reaction Conditions to Improve the Magnetic Response of Cobalt-Doped Ferrite Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First-Principles Study on the Structural and Electronic Properties of Monolayer MoS2 with S-Vacancy under Uniaxial Tensile Strain

1
School of Mechano-Electronic Engineering, Xidian University, Xi’an 710071, China
2
Department of Mechanical Engineering, Northwestern University, Evanston, IL 60208, USA
3
School of Mechanical Engineering and Automation, Fuzhou University, Fuzhou 350108, China
4
ZNDY of Ministerial Key Laboratory, Nanjing University of Science and Technology, Nanjing 210094, China
5
McCormick School of Engineering and Applied Science, Northwestern University, Evanston, IL 60208, USA
*
Authors to whom correspondence should be addressed.
Nanomaterials 2018, 8(2), 74; https://doi.org/10.3390/nano8020074
Submission received: 24 November 2017 / Revised: 19 January 2018 / Accepted: 25 January 2018 / Published: 29 January 2018
(This article belongs to the Special Issue Experimental Nanosciences, Computational Chemistry, and Data Analysis)

Abstract

:
Monolayer molybdenum disulfide (MoS2) has obtained much attention recently and is expected to be widely used in flexible electronic devices. Due to inevitable bending in flexible electronic devices, the structural and electronic properties would be influenced by tensile strains. Based on the density functional theory (DFT), the structural and electronic properties of monolayer MoS2 with a sulfur (S)-vacancy is investigated by using first-principles calculations under uniaxial tensile strain loading. According to the calculations of vacancy formation energy, two types of S-vacancies, including one-sulfur and two-sulfur vacancies, are discussed in this paper. Structural analysis results indicate that the existence of S-vacancies will lead to a slightly inward relaxation of the structure, which is also verified by exploring the change of charge density of the Mo layer and the decrease of Young’s modulus, as well as the ultimate strength of monolayer MoS2. Through uniaxial tensile strain loading, the simulation results show that the band gap of monolayer MoS2 decreases with increased strain despite the sulfur vacancy type and the uniaxial tensile orientation. Based on the electronic analysis, the band gap change can be attributed to the π bond-like interaction between the interlayers, which is very sensitive to the tensile strain. In addition, the strain-induced density of states (DOS) of the Mo-d orbital and the S-p orbital are analyzed to explain the strain effect on the band gap.

1. Introduction

Since the discovery of graphene, 2D materials such as boron nitride (BN), molybdenum disulfide (MoS2), and tungsten disulfide (WS2) have become one of the hot topics in the scientific community because of the fantastic physical properties and promising applications of these materials in flexible electronic devices [1]. Among the 2D materials mentioned above, MoS2 is one candidate for flexible electronic devices because of its good mechanical and electrical properties. The bulk of MoS2 belongs to the space group P63/mmc. The monolayer MoS2 can be viewed as a cleaved form of the (001) surface of bulk MoS2 [2]. It is different from the structure of graphene in that it has three atomic layers composed of one Mo layer plus two S layers on both sides, which are held together by van der Waals interactions [3]. During its application, flexible electronic material will inevitably suffer bending deformations. Up to now, the preparation of ultra-thin MoS2 films, even a monolayer MoS2 structure on a substrate, can be successfully realized by many methods such as electron irradiation and chemical vapor deposition (CVD) [4]. However, vacancy defects can be produced with high density on the MoS2 surface, and this can influence the physical properties and band structures. Therefore, it is very important to explore the influence of these factors on its mechanical and electrical properties.
Recently, much interest has been focused on the strain or vacancy effect on the properties of MoS2 from both experimental [5,6,7,8] and theoretical [9,10,11,12] aspects. Cooper et al. [5] explored the nonlinear elastic properties by indenting suspended circular MoS2 membranes with an atomic force microscope. Bertolazzi et al. [6] found that the in-plane stiffness of monolayer MoS2 is 180 ± 60 N·m−1, which corresponds to an effective Young’s modulus of 270 ± 100 GPa. Wang et al. [7] observed a strong shrinkage of the band gap with a rate of 48 meV/% of strain for monolayer MoS2. Santosh et al. [8] examined by scanning tunneling microscopy the possible defect structures and their impact on the electronic properties of the MoS2 monolayer. In addition to the experimental work, there has been a lot of work on MoS2 such as molecular dynamics (MD) simulations [9,10] and first-principles studies [2,5,11,12]. Wang and his co-authors [9] used nanoindentation simulations to find Young’s modulus of 280 ± 21 GPa. Li et al. [10] investigated the effect of VMoS3 point defects on the elastic properties of monolayer MoS2 sheets by using the classical MD simulation with reactive empirical bond-order (REBO) potential. Cooper et al. [5] applied density functional calculations based on local density approximation (LDA) and estimated Young’s modulus of MoS2 nanosheets equal to 210 GPa. Lu et al. [2] explored through first-principles calculations the tensile or compressive axial strain in different charities and found different effects on the band gap. Recently, much work has been done to calculate the properties of strained MoS2 using the DFT method. For instance, Spirko et al. [11] investigated the defect structure of monolayer MoS2, Defo et al. [12] proved the effect of strain on band gap regulation by applying compressive and tensile strain to single layer MoS2, and Maniadaki et al. [13] calculated the band gap of transition metal dichalcogenide monolayers and found it to decrease while dielectric constants increased for heavier chalcogens X. In addition, Coehoorn et al. [14] investigated the defect structure of monolayer MoS2 using DFT calculations and Lebègue et al. [15] calculated the electronic structure of MoS2 using the ab-initio method. All of these studies have many important conclusions; however, there is little research about the influence of vacancy on the structural and electronic properties of MoS2 when its energy gap is controlled by uniaxial tensile strain. Therefore, it is necessary to investigate the influence of these two factors comprehensively.
In the present work, we investigate the structural and electronic properties of monolayer MoS2 with S-vacancy defects under tensile strain loading by using first-principles calculations. The effect of tensile strain on the density of electronic states is revealed. The mechanical properties of monolayer MoS2 are weakened due to the existence of defects, and the band gap of monolayer MoS2 undergoes a descending trend with increasing strain, which means a semiconductor with a lower band or metal transition can be achieved with the application of strain.

2. Method

First-principle calculations were carried out using density functional theory (DFT) calculations in the projector augmented-wave (PAW) pseudopotential as implemented within the Vienna Ab-initio Simulation Package (VASP) [16,17]. The electronic exchange-correlation potential was used by the Perdew–Burke–Ernzerhof (PBE) flavor [18]. Electronic kinetic-energy cutoff was set to 450 eV [19]. The Brillouin zone integration was carried out by using a 12 × 12 × 1 k-mesh according to the Monkhorst–Pack scheme. For the physical modeling of monolayer MoS2, the supercell kept periodical boundary conditions in the x- and y-directions, as well as during DFT simulations, and a large enough vacuum region in the z-direction (20 Å) was set to prevent inter-layer interactions.
The structure was fully relaxed with an energy convergence of 1.0 × 10−6 eV. The atomic Hellman–Feynman force was minimized to less than 0.01 eV/Å. Especially, a rectangle supercell was chosen in our calculations, and the red dashed frame indicates the original cell that involves four Mo atoms and eight S atoms. The 3 × 3 × 1 supercell structure of a perfect monolayer MoS2 from top view in x-and y-directions is shown in Figure 1a and its side view is given in Figure 1b. The original cell is delineated with a red dashed line frame including 2 × 1 primitive cells. The first Brillouin zone is described in Figure 1c for this model.
In the present study, the monolayer MoS2 with a sulfur-vacancy defect under different tensile strains along the x- and y-directions was investigated, respectively. As shown in Figure 2a,c, we considered three MoS2 lattice structures: a perfect lattice (PL), a one-sulfur vacancy (V1S) lattice, and a two-sulfur vacancy (V2S) lattice. To further understand the vacancy defect, we calculated the vacancy formation energy for the lattice with a vacancy in Figure 2b,c. The vacancy formation energy Ef is defined as Ef = EdEp + nEi where Ed is the total energy of the vacancy lattice, Ep is the corresponding perfect lattice, n is the number of atoms, and Ei is the chemical potential energy of the S atoms. We calculated the vacancy formation energies of V1S and V2S to be 6.59 eV and 13.08 eV, respectively. In addition, we calculated the Mo vacancy formation energy, which was much higher than the S vacancy, indicating that the S vacancy is more easily produced. The uniaxial tensile strain was simulated by enlarging the x and y cell lattice parameters, respectively, to a fixed larger value; then, the atomic positions were fully optimized. It is worth noting that in the original version of VASP, it cannot be achieved that a lattice parameter is fixed while other lattice constants are optimized. For this reason, we complied with the original VASP source code to enable the realization of this function.

3. Results and Discussion

3.1. Structural Properties

In order to find the equilibrium lattice constant, full relaxations were conducted on the three models in Figure 2a–c. The lattice constant is 3.166 Å, the bond length of Mo–S and S–S are 2.415 Å and 3.131 Å, respectively, which agrees well with the available numerical results [20]. Figure 3 shows the atomic structures of the one-sulfur vacancy (V1S) and the two-sulfur vacancy (V2S) in monolayer MoS2. The bonds in the black dotted circle are colored according to a decrease (blue) or increase (red) in the bond length. The effects of a dark color are stronger than those of the light color.
Compared to the perfect supercell model, there is little change in the lattice constant in the supercell containing vacancies after relaxation, indicating that a small amount of vacancy does not produce a significant impact on the crystal structure; however, the bond length around the vacancy is slightly different. Table 1 lists the lengths of the bonds of V1S and V2S of monolayer MoS2 before and after structure relaxation compared with perfect MoS2 in which a, b, c, d, and e represent the distance between the two adjacent Mo–S atoms in the sandwich layer around the vacancy. It is intuitively clear that ions surrounding S vacancies appear with slight inward relaxation because of the absence of sulfur atoms. For S vacancies in Figure 3b,c, the inward relaxations of the nearest Mo ions are larger than those of the other Mo ions [21]. The phenomenon can be explained by the redistribution of electrons. Compared with perfect MoS2, the adjacent S–Mo bond change is slightly larger than the sub-adjacent one, which conforms to the lattice distortion of a point defect.
To calculate Young’s modulus (E) and ultimate strength (σmax), we applied the uniaxial tensile strains on three models in the x- and y-directions, respectively. The applied strains in the x- and y-directions can be defined by Equation (1) [22]:
ϵ x = a a 0 a 0   and   ϵ y = b b 0 b 0 ,
where a and b are the lengths of the supercell in the x- and y-directions after equilibrium and the lattice constants are a0 = 16.5672 Å and b0 = 19.1307 Å. The strains (ϵ) range from 0.00 to 0.30 with an increment of 0.01 in each step for three deformation cases. There are 90 ab initio DFT calculations in total. The per atom strain energy (Es) with tensile strain for PL, V1S, and V2S sheet in the x- and y-directions are shown in Figure 4. For small strains, Es goes up with increased uniaxial tensile strain, while the influence of both uniaxial tensile orientation and sulfur vacancy types on Es are very small. Once the strain goes beyond ~23%, Es goes up with an increase in strain in the x-direction but slows in the y-direction. These strain–energy curves suggest the anisotropic properties of monolayer MoS2 with the tensile strain direction [23].
Young’s modulus can be derived through the second derivative of the strain energy with respect to the strain and can be expressed by Equation (2) [24]:
E = 1 V 0 2 U ϵ 2 | ϵ = 0
where V0 is the volume of monolayer MoS2, the thickness of it is defined as 6.145 Å, and U is the strain energy calculated by subtracting the total energy of the strained system from the equilibrium state. As shown in Figure 5, the strain–stress curves are plotted for the PL, V1S, and V2S monolayer MoS2 sheets, respectively. The tensile orientation effect on Young’s modulus is so slight that it can be neglected at small strains; however, the effect becomes apparent when the strain increases.
From Figure 5 it can be found that along the x-direction both the maximum stress and its corresponding strain are higher than those along the y-direction. The ultimate strength for a perfect lattice in the x-direction is larger than that in the y-direction [6]. In addition, the S-vacancy defects reduce the energy of the system. Compared with the PL MoS2 sheet, the S-vacancy MoS2 sheets reach the yielding point more quickly and, moreover, the yielding point of x is higher than the y-direction. In addition, for the V1S case, the symmetry is destroyed due to the sulfur vacancy and the oscillation occurs at the yielding point in the x-direction until the failure point is reached. In contrast, V2S reaches steady growth after reaching the yielding point until it reaches a point of failure. As shown in Table 2, the existence of a vacancy lowers the strength and stiffness of the material. The effects of orientation and sulfur defects on Young’s modulus and ultimate stress are not obvious. Moreover, Bertolazzi et al. have approximately obtained Young’s modulus of monolayer MoS2, which is equal to 270 ± 100 GPa [4]; Castellanos–Gomez et al. calculated the average of monolayer molybdenum disulfide to be around 330 ± 70 GPa using atomic force microscopy (AFM) [25]. Therefore, the simulation results of Young’s modulus in the present study agree well with the experimental data.

3.2. Electronic Properties

Strain modulation has been commonly used in low-dimensional systems to tune the electronic structures and band gaps. We further calculated the band gap of monolayer MoS2 sheets with two different types of sulfur-vacancy defects influenced by strain in the x- and y-directions, as shown in Figure 6. It can be found that the band gaps of the three cases—including PL, V1S, and V2S—decrease gradually as the tensile strain increases. The tendencies in the x- and y-directions are basically the same; the only difference is that the band gap of the monolayer MoS2 with sulfur vacancy defects will be reduced to zero. As shown in Figure 6, the calculated direct band gap for perfect monolayer MoS2 is about 1.77 eV at zero strain, which is consistent with the experimental value (1.90 eV [26]) and theoretical values (1.80 eV [27], 1.70 eV [28]). Both the V1S and V2S monolayer MoS2 sheets behave as semiconductors with band gaps of 0.72 eV and 0.73 eV, respectively. The reason why they changed a lot compared to the perfect one (PL) is mainly that the sulfur vacancy defect makes the original band gap impure.
Considering the tensile strain, the electronic property of the perfect and S-vacancy MoS2 are revealed through investigating their band structures, which are shown in Figure 7. From the band structure, it is found that the presence of sulfur vacancies leads to the introduction of new energy bands near the Fermi level and the hydrogen adsorption on the S-vacancies may be the reason for it. Because of the introduction of the sulfur vacancy, the direct band gap of monolayer MoS2 has been changed into an indirect band gap, which is due to the contraction of the neighboring atoms of S-vacancy; a local state or transition state occurs between the VBM and CBM [10,19]. A defect energy band is introduced (red line). The direct band gap of monolayer MoS2 has been changed into the indirect band gap, which is due to the contraction of the neighboring atoms of the S-vacancy. A local state or transition state occurs between the VBM and CBM [10,19].
With an increase in the tensile strain, the conduction band minimum (CBM) moves down but the valence band maximum (VBM) does not visibly change; the CBM point moves from the Γ point, which leads to the band gap changing from direct into indirect. At the same time, the influence of orientation is not obvious, so the x-direction is taken as an example.
Figure 8 shows the density of states (DOS) of the PL, V1S, and V2S monolayer MoS2 and their 10% strain under the uniaxial tension in the x-direction. For the perfect monolayer MoS2, as well as one with defects, CBM and VBM are mainly composed of S-p orbitals and Mo-d orbitals. At the same time for V1S and V2S defects of MoS2, as shown in Figure 8a–c, taking V1S as an example, the existence of a transition state that reduces the band gap is mainly due to the interaction above the Fermi level of the S-p and Mo-d orbitals. The main reason for the transition state is that the sulfur vacancy causes contraction of the Mo–Mo bond around the defect. At the same time, the highly symmetric bond structure of S–Mo–S forms the role of a weak π bond-like interaction, which is extremely sensitive to strain, resulting in a significant change in the energy band structure once the strain is applied. As shown in Figure 8d–f, after the strain is applied, the DOS around the Fermi level changes significantly and the variation of the S-p orbit is more obvious than that of the Mo atomic orbital. The PL MoS2 sheet presents a direct band gap structure and there is no other energy level in the band gap. Compared with the PL electronic structure, the presence of S-vacancies makes the three Mo atoms around the vacancies have six Mo-d orbital electrons that are not bonded with the S atoms, which are acting as electron donors in the defect system. The excess Mo atoms’ Mo-d orbit electrons fill into these defect levels. The analysis holds that the two levels are more likely to belong to the acceptor’s energy level.
We also examined the electron densities of the Mo layer, which can reflect changes in the chemical bonds around the vacancies. Figure 9 compares the electron density of the Mo layer between intrinsic and sulfur vacancies of the monolayer MoS2. It can be seen from the figure that in the absence of sulfur in a symmetrical position in the MoS2 structure, the electrons near the vacancy are concentrated mainly in three adjacent Mo atoms. The high density of electrons gathered makes the three Mo atoms almost together. Also, a strong coulomb attraction is produced to the nearest Mo atom or S atom and that leads to the S–Mo bonds in the foreground becoming smaller. Compared with the perfect state, the electron distribution near the vacancy area is obviously localized and should be related to the formation of the defect level.

4. Conclusions

In summary, through DFT calculations, the structural and electronic properties of monolayer MoS2 with one or two sulfur vacancies have been investigated systematically under uniaxial tensile strain in the x- and y-directions. Regarding structural properties, it was found that the existence of vacancies can weaken the stiffness and strength of the MoS2 sheet. Young’s modulus and strength of perfect MoS2 have been determined to be 315 GPa and 33.2 GPa, respectively, in the x-direction; the one-sulfur vacancy MoS2 reduces to 272 GPa and 21.9 GPa, respectively. Differences exist in the weakening effect in different uniaxial tensile loading directions (i.e., the x- and y-directions).
Besides, an analysis of the electronic properties shows that both of the two vacancy types behave as semiconductors and that their band gap is much smaller than perfect monolayer MoS2 because of the impurity. The band gap of perfect monolayer MoS2 is about 1.77 eV while the one and two sulfur-vacancy MoS2 are 0.72 eV and 0.73 eV. The strains also have a huge impact on the band gap of defected MoS2 sheets because of the relaxation near the vacancy. By comparing the electronic structures with that of perfect MoS2, the existence of vacancy defects has obvious effects on the electronic structures of monolayer MoS2, especially for the high-energy region of the conduction band density of state. These effects may be related to the defect energy levels introduced by the vacancy defects. Above all, our study explored MoS2 subjected to external uniaxial tensile strain, which can present band gap changes to tune its resistance, which makes it applicable in MoS2 piezo-resistive sensors. At the same time, the effects of sulfur vacancies on the mechanical properties and the band gaps of the material have been studied to determine the measuring range of the sensor and the interference from the vacancy.

Acknowledgments

We acknowledge the financial support provided by the National Natural Foundation of China (grant NOs. 51205302 and 50903017) and the Natural Science Basic Research Plan in Shaanxi Province of China (grant NO. 2017JM5003).

Author Contributions

All of the authors made important contributions on this work. Specifically, W.W. and W.L. proposed the project and analyzed the simulation results. C.Y. and L.B. contributed to the DFT simulations. W.W., M.L. and W.L. composed the manuscript. All the authors participated in the discussions of the results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Topsakal, M.; Cahangirov, S.; Ciraci, S. The response of mechanical and electronic properties of graphane to the elastic strain. Appl. Phys. Lett. 2010, 96, 091912. [Google Scholar] [CrossRef] [Green Version]
  2. Lu, P.; Wu, X.; Guo, W.; Zeng, X.C. Strain-dependent electronic and magnetic properties of MoS2 monolayer, bilayer, nanoribbons and nanotubes. Phys. Chem. Chem. Phys. 2012, 14, 13035–13040. [Google Scholar] [CrossRef] [PubMed]
  3. Scalise, E.; Houssa, M.; Pourtois, G.; Afanas’ev, V.V.; Stesmans, A. First-principles study of strained 2D MoS2. Physica E 2014, 56, 416–421. [Google Scholar] [CrossRef]
  4. Tao, P.; Guo, H.; Yang, T.; Zhang, Z. Strain-induced magnetism in MoS2 monolayer with defects. J. Appl. Phys. 2014, 115, 054305. [Google Scholar] [CrossRef]
  5. Cooper, R.C.; Lee, C.; Marianetti, C.A.; Wei, X.; Hone, J.; Kysar, J.W. Nonlinear elastic behavior of two-dimensional molybdenum disulfide. Phys. Rev. B 2013, 87, 035423. [Google Scholar]
  6. Bertolazzi, S.; Brivio, J.; Kis, A. Stretching and breaking of ultrathin MoS2. ACS Nano 2011, 5, 9703–9709. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, M.; Edmonds, K.W.; Gallagher, B.L.; Rushforth, A.W.; Makarovsky, O.; Patanè, A.; Campion, R.P.; Foxon, C.T.; Novak, V.; Jungwirth, T. High Curie temperatures at low compensation in the ferromagnetic semiconductor (Ga,Mn)As. Phys. Rev. B 2013, 87, 121301. [Google Scholar] [CrossRef]
  8. Santosh, K.C.; Longo, R.C.; Addou, R.; Wallace, R.M.; Cho, K. Impact of intrinsic atomic defects on the electronic structure of MoS2 monolayers. Nanotechnology 2014, 25, 375703. [Google Scholar]
  9. Wang, W.; Li, L.; Yang, C.; Soler-Crespo, R.A.; Meng, Z.; Li, M.; Zhang, X.; Keten, S.; Espinosa, H.D. Plasticity resulted from phase transformation for monolayer molybdenum disulfide film during nanoindentation simulations. Nanotechnology 2017, 28, 164005. [Google Scholar] [CrossRef] [PubMed]
  10. Li, M.; Wan, Y.; Tu, L.; Yang, Y.; Lou, J. The effect of VMoS3 point defect on the elastic properties of monolayer MoS2 with REBO potentials. Nanoscale Res. Lett. 2016, 11, 155. [Google Scholar] [CrossRef] [PubMed]
  11. Spirko, J.A.; Neiman, M.L.; Oelker, A.M.; Klier, K. Electronic structure and reactivity of defect MoS2 II. Bonding and activation of hydrogen on surface defect sites and clusters. Surf. Sci. 2004, 572, 191–205. [Google Scholar] [CrossRef]
  12. Defo, R.K.; Fang, S.; Shirodkar, S.N.; Tritsaris, G.A.; Dimoulas, A.; Kaxiras, E. Strain dependence of band gaps and exciton energies in pure and mixed transition-metal dichalcogenides. Phy. Rev. B 2016, 94, 155310. [Google Scholar] [CrossRef]
  13. Maniadaki, A.E.; Kopidakis, G.; Remediakis, I.N. Strain engineering of electronic properties of transition metal dichalcogenide monolayers. Solid State Commun. 2016, 227, 33–39. [Google Scholar] [CrossRef]
  14. Coehoorn, R.; Haas, C.; de Groot, R.A. Electronic structure of MoSe2, MoS2, and WSe2. II. The nature of the optical band gaps. Phys. Rev. B. 1987, 35, 6203. [Google Scholar] [CrossRef]
  15. Lebègue, S.; Eriksson, O. Electronic structure of two-dimensional crystals from ab initio theory. Phy. Rev. B 2009, 79, 115409. [Google Scholar] [CrossRef]
  16. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. [Google Scholar] [CrossRef]
  17. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758. [Google Scholar] [CrossRef]
  18. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef] [PubMed]
  19. Yue, Q.; Kang, J.; Shao, Z.; Zhang, X.; Chang, S.; Wang, G.; Qin, S.; Li, J. Mechanical and electronic properties of monolayer MoS2 under elastic strain. Phys. Lett. A 2012, 376, 1166–1170. [Google Scholar] [CrossRef]
  20. Lu, N.; Guo, H.; Li, L.; Dai, J.; Wang, L.; Mei, W.-N.; Wu, X.; Zeng, X.C. MoS2/MX2 heterobilayers: Bandgap engineering via tensile strain or external electrical field. Nanoscale 2014, 6, 2879–2886. [Google Scholar] [CrossRef] [PubMed]
  21. Krasheninnikov, A.V. Strain fields and electronic structure of vacancy-type defects in graphene from first-principles simulation. In Nanodevices and Nanomaterials for Ecological Security; Spriger Nature: New York, NY, USA, 2012. [Google Scholar]
  22. Van der Zande, A.M.; Huang, P.Y.; Chenet, D.A.; Berkelbach, T.C.; You, Y.; Lee, G.-H.; Heinz, T.F.; Reichman, D.R.; Muller, D.A.; Hone, J.C. Grains and grain boundaries in highly crystalline monolayer molybdenum disulphide. Nat. Mater. 2013, 12, 554–561. [Google Scholar] [CrossRef] [PubMed]
  23. Ansari1, R.; Malakpour, S.; Faghihnasiri, M.; Ajori, S. Characterization of the mechanical properties of monolayer molybdenum disulfide nanosheets using first principles. J. Nanotechnol. Eng. Med. 2014, 4, 034501. [Google Scholar] [CrossRef]
  24. Gan, Y.; Zhao, H. Chirality effect of mechanical and electronic properties of monolayer MoS2 with vacancies. Phys. Lett. A 2014, 378, 2910–2914. [Google Scholar] [CrossRef]
  25. Priezjev, N. Interfacial friction between semiflexible polymers and crystalline surfaces. J. Chem. Phys. 2012, 136, 224702. [Google Scholar] [CrossRef] [PubMed]
  26. Böker, T.; Severin, R.; Müller, A.; Janowitz, C.; Manzke, R.; Voß, D.; Krüger, P.; Mazur, A.; Pollmann, J. Band structure of MoS2, MoSe2, and α-MoTe2: Angle-resolved photoelectron spectroscopy and ab initio calculations. Phys. Rev. B 2001, 64, 235305. [Google Scholar] [CrossRef]
  27. Remskar, M.; Mrzel, A.; Virsek, M.; Godec, M.; Krause, M.; Kolitsch, A.; Singh, A.; Seabaugh, A. The MoS2 nanotubes with defect-controlled electric properties. Nanoscale Res. Lett. 2011, 6, 26. [Google Scholar] [CrossRef] [PubMed]
  28. Miao, Y.-P.; Ma, F.; Huang, Y.-H.; Xu, K.-W. Strain effects on electronic states and lattice vibration of monolayer MoS2. Physica E 2015, 71, 1–6. [Google Scholar] [CrossRef]
Figure 1. The 3 × 3 × 1 supercell structure of a perfect molybdenum disulfide (MoS2) monolayer: (a) top view and (b) side view. The original cell is delineated with the red dashed line frame and the first Brillouin zone is described (c).
Figure 1. The 3 × 3 × 1 supercell structure of a perfect molybdenum disulfide (MoS2) monolayer: (a) top view and (b) side view. The original cell is delineated with the red dashed line frame and the first Brillouin zone is described (c).
Nanomaterials 08 00074 g001
Figure 2. The top view of the atomic configurations of monolayer MoS2: (a) perfect lattice (PL) structure, (b) one-sulfur vacancy (V1S) structure, and (c) two-sulfur vacancy (V2S) structure.
Figure 2. The top view of the atomic configurations of monolayer MoS2: (a) perfect lattice (PL) structure, (b) one-sulfur vacancy (V1S) structure, and (c) two-sulfur vacancy (V2S) structure.
Nanomaterials 08 00074 g002
Figure 3. The atomic structures of PL (a), V1S (b), and V2S (c) in a monolayer MoS2. The bonds in the black dotted circle are colored according to the decrease (blue) or increase (red) in the bond length. The effects of a darker color are stronger than those of the lighter color.
Figure 3. The atomic structures of PL (a), V1S (b), and V2S (c) in a monolayer MoS2. The bonds in the black dotted circle are colored according to the decrease (blue) or increase (red) in the bond length. The effects of a darker color are stronger than those of the lighter color.
Nanomaterials 08 00074 g003
Figure 4. The per-atom strain energy Es of PL, V1S, and V2S monolayer MoS2 sheet versus tensile strain in the x- and y-directions.
Figure 4. The per-atom strain energy Es of PL, V1S, and V2S monolayer MoS2 sheet versus tensile strain in the x- and y-directions.
Nanomaterials 08 00074 g004
Figure 5. The variations of stress for the PL, V1S, and V2S monolayer MoS2 sheets versus the uniaxial tensile strain in the x- and y-directions.
Figure 5. The variations of stress for the PL, V1S, and V2S monolayer MoS2 sheets versus the uniaxial tensile strain in the x- and y-directions.
Nanomaterials 08 00074 g005
Figure 6. The band gaps of the PL (a), V1S (b), and V2S (c) monolayer MoS2 variation with the tensile strain in the x- and y-directions.
Figure 6. The band gaps of the PL (a), V1S (b), and V2S (c) monolayer MoS2 variation with the tensile strain in the x- and y-directions.
Nanomaterials 08 00074 g006aNanomaterials 08 00074 g006b
Figure 7. The band structure of the PL (a), V1S (b), and V2S (c) monolayer MoS2 sheet at natural state (zero strain) and under 10% uniaxial tensile strain in the x-direction (df). The dotted line represents the Fermi energy level.
Figure 7. The band structure of the PL (a), V1S (b), and V2S (c) monolayer MoS2 sheet at natural state (zero strain) and under 10% uniaxial tensile strain in the x-direction (df). The dotted line represents the Fermi energy level.
Nanomaterials 08 00074 g007
Figure 8. The electron density of states (DOS) for the PL (a), V1S (b) and V2S (c) monolayer MoS2 sheet at natural state (zero strain) and under 10% uniaxial tensile strain in the x-direction (df). Zero energy represents the Fermi energy level.
Figure 8. The electron density of states (DOS) for the PL (a), V1S (b) and V2S (c) monolayer MoS2 sheet at natural state (zero strain) and under 10% uniaxial tensile strain in the x-direction (df). Zero energy represents the Fermi energy level.
Nanomaterials 08 00074 g008
Figure 9. The electron density of a Mo layer. (a) The Mo bonding in intrinsic monolayer MoS2. (b) The Mo bonding in two-sulfur vacancy of monolayer MoS2.
Figure 9. The electron density of a Mo layer. (a) The Mo bonding in intrinsic monolayer MoS2. (b) The Mo bonding in two-sulfur vacancy of monolayer MoS2.
Nanomaterials 08 00074 g009
Table 1. The length of bonds of V1S and V2S of the monolayer MoS2 before and after structure relaxation compared with the perfect MoS2 (Unit: Å).
Table 1. The length of bonds of V1S and V2S of the monolayer MoS2 before and after structure relaxation compared with the perfect MoS2 (Unit: Å).
Model Typeabcde
PL2.4152.4152.4152.4152.415
V1S2.4302.4502.3982.4022.371
V2S2.4282.4532.4002.405/
Table 2. Young’s modulus (E) and maximum strength (σmax) of monolayer MoS2 for a perfect film and one-sulfur (V1S) and two-sulfur (V2S) vacancies.
Table 2. Young’s modulus (E) and maximum strength (σmax) of monolayer MoS2 for a perfect film and one-sulfur (V1S) and two-sulfur (V2S) vacancies.
E (GPa) σ max (GPa)
xyxy
Perfect31533533.017.7
V1S27229022.117.1
V2S22724116.115.1

Share and Cite

MDPI and ACS Style

Wang, W.; Yang, C.; Bai, L.; Li, M.; Li, W. First-Principles Study on the Structural and Electronic Properties of Monolayer MoS2 with S-Vacancy under Uniaxial Tensile Strain. Nanomaterials 2018, 8, 74. https://doi.org/10.3390/nano8020074

AMA Style

Wang W, Yang C, Bai L, Li M, Li W. First-Principles Study on the Structural and Electronic Properties of Monolayer MoS2 with S-Vacancy under Uniaxial Tensile Strain. Nanomaterials. 2018; 8(2):74. https://doi.org/10.3390/nano8020074

Chicago/Turabian Style

Wang, Weidong, Chenguang Yang, Liwen Bai, Minglin Li, and Weibing Li. 2018. "First-Principles Study on the Structural and Electronic Properties of Monolayer MoS2 with S-Vacancy under Uniaxial Tensile Strain" Nanomaterials 8, no. 2: 74. https://doi.org/10.3390/nano8020074

APA Style

Wang, W., Yang, C., Bai, L., Li, M., & Li, W. (2018). First-Principles Study on the Structural and Electronic Properties of Monolayer MoS2 with S-Vacancy under Uniaxial Tensile Strain. Nanomaterials, 8(2), 74. https://doi.org/10.3390/nano8020074

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