Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5363068

Formats

Article sections

- Abstract
- Overview of the SCAN Methodology
- Computational Details
- Results and Discussion
- Conclusions
- Additional Information
- References

Authors

Related links

Sci Rep. 2017; 7: 44766.

Published online 2017 March 23. doi: 10.1038/srep44766

PMCID: PMC5363068

I. G. Buda,^{a,}^{1} C. Lane,^{1} B. Barbiellini,^{1} A. Ruzsinszky,^{2} J. Sun,^{3} and A. Bansil^{1}

Received 2016 October 19; Accepted 2017 February 14.

Copyright © 2017, The Author(s)

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

We discuss self-consistently obtained ground-state electronic properties of monolayers of graphene and a number of ’beyond graphene’ compounds, including films of transition-metal dichalcogenides (TMDs), using the recently proposed strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) to the density functional theory. The SCAN meta-GGA results are compared with those based on the local density approximation (LDA) as well as the generalized gradient approximation (GGA). As expected, the GGA yields expanded lattices and softened bonds in relation to the LDA, but the SCAN meta-GGA systematically improves the agreement with experiment. Our study suggests the efficacy of the SCAN functional for accurate modeling of electronic structures of layered materials in high-throughput calculations more generally.

Discovery of graphene, a one-atom-thick crystal of carbon, has spurred an intense interest in the electronic properties of 2D materials more generally^{1}^{,2}. Recent research has turned to ‘beyond graphene’ materials, which exhibit novel spin and charge transport properties, including the possibility of harboring quantum spin Hall and other topological phases^{3} relevant for next generation electronics applications and as materials platforms for replacing the current Si-based technologies. For example, unlike the flat structure of graphene, silicene, germanene and stanene, which are Si, Ge and Sn based cousins of graphene, assume a crystal structure that is naturally buckled^{4}^{,5}^{,6}. As a result, these materials exhibit spin-split states, which can be controlled via external electric fields^{7}. Phosphorene displays remarkable mechanical flexibility and sensitive tuning of electronic properties by mechanical strain^{8}. Ultra-thin films of transition metal dichalcogenides (TMDs) undergo a transition from an indirect to a direct band gap semiconductor in the monolayer limit, and have become attractive candidates for nanoelectronics^{9}, water-splitting^{10}, photocatalysis^{11} and other applications.

The need for theoretical methods capable of accurate and efficient prediction of structural and electronic properties of atomically thin films and layered materials is clear. In this connection, improvements in density functional theory (DFT)^{12} based first-principles computations, which have been the workhorse in the field for over five decades^{13}^{,14}, have centered around the development of new classes of exchange-correlation functionals. One of the latest advances in this direction is the SCAN meta-GGA scheme, which has been proposed recently^{15}. Our purpose in this study is to assess the efficacy of the SCAN meta-GGA for addressing the ground state properties of 2D materials. SCAN meta-GGA has been tested in diversely bonded systems^{16}, where it has been shown to capture a wide range of physical structures without being fitted to any specific type of bonding. These SCAN-based existing studies include: MnO_{2} polymorphs^{17}; Cu-intercalated birnessite^{18}; and, band gaps of semiconductors and insulators^{19}. Here we show that SCAN meta-GGA yields a systematic improvement over the LDA and GGA (at a comparable cost) in modeling ground state properties of 2D materials. For this purpose, we consider the application of SCAN to monolayers of graphene and a number of ‘beyond graphene’ compounds, including films of transition-metal dichalcogenides (TMDs) as exemplar 2D systems.

An outline of this article is as follows. The introductory remarks above are followed by an overview of the SCAN functional and its construction. We then describe the relevant computational details, followed by a presentation and discussion of our results, and a summary of our conclusions.

Within the framework of the DFT, the total energy of the many-body electron system, *E*_{total}[*n*], can be written in terms of the electron density, *n*, as

where *K* is the independent-electron kinetic energy, *E*_{ie} is the Coulomb energy between the electrons and ions, *E*_{ee} describes the classical electron-electron Coulomb interaction, and *E*_{xc} the exchange-correlation energy. Approximation schemes for *E*_{xc} can be arranged conceptually on the rungs of the so-called DFT *Jacob’s Ladder*^{20} in the sense that this ladder leads to the ‘heaven’ of chemical accuracy. Various rungs of this ladder, beginning with the lowest rung, are: LDA^{21}^{,22}; GGA^{23}^{,24}; meta-GGA^{25}^{,26}; hybrid functionals^{27}; and, finally the random phase approximation (RPA)^{28}. Computational demands, along with the accuracy of the schemes increase as we go up the rungs of the ladder.

Formally, the *E*_{xc}[*n*] term can be cast as a double integral over space, which involves half of the Coulomb interaction between electrons and the associated exchange-correlation holes^{26}^{,29}, but it is computationally expensive to evaluate. In the semilocal approximation, this term is reduced to a single integral of the general form

where is the electron density, *n* its gradient, the positive orbital kinetic energy density, and Ψ_{i,σ} are the Kohn-Sham orbitals. Nonempirical functionals are generally built to satisfy exact constraints as far as possible. It is here that the SCAN meta-GGA^{15} makes a substantial advance as it is the only semilocal exchange-correlation functional which satisfies the complete set of 17 known exact constraints that can be satisfied by semilocal functionals. Moreover, SCAN is ‘appropriately normed’ in that it accurately captures interactions in rare-gas atoms and unbonded systems (see Supplementary Material of Sun *et al*.^{15}). The earlier nonempirical meta-GGAs such as the Tao-Perdew-Staroverov-Scuseria (TPSS)^{25} and revTPSS^{30} meta-GGA have been shown to be less accurate than the Perdew-Burke-Ernzerhof (PBE) GGA for the critical pressures of structural phase transitions of solids^{31}^{,32}. SCAN meta-GGA eliminates this problem by introducing the dimensionless parameter

where is the single-orbital limit of *τ*, and is the uniform density limit. The case of *α*=0 corresponds to covalent single bonds while *α*≈1 to metallic, and the *α*1 limit describes weak bonds. The rare-gas-atom norm contains information about 0<*α*<∞, and some information about *α*1, while the non-bonded-interaction norm (the compressed Ar_{2}) provides more information about *α*1.

Recently, SCAN meta-GGA has been tested in diversely bonded systems^{16}, where it was shown to be sophisticated enough to model a wide range of physical structures without being fitted to any bonded system. In the present work, we apply it further to the class of thin film materials and we show a similar trend of successful predictions of ground-state structural and electronic properties. In particular, SCAN improves the overall agreement with experiment compared to LDA and GGA, at a comparable computational cost.

We have performed first-principles calculations using the pseudopotential projector augmented-wave method^{33} as implemented in the Vienna Ab-Initio Simulation Package (VASP)^{34}^{,35}, with a kinetic energy cutoff of 400eV (TMD monolayers and Bi_{2}Se_{3} quintuple layer) and 800eV (graphene, silicene, germanene, and phosphorene) for the plane-wave basis set. The exchange-correlation functional was treated using LDA^{15}^{,36}, GGA-PBE^{23}^{,24} and SCAN meta-GGA^{15}. A 12×12×1 Γ-centered *k*-point mesh was used to sample the Brillouin zone. Spin-orbit coupling effects were included in the case of TMD monolayers and Bi_{2}Se_{3} quintuple layer in a self-consistent manner. We used a vacuum layer of at least 15Å thickness in the z-direction to simulate the films. The equilibrium positions of the ions were calculated via structural optimization, where the internal degrees of freedom, along with the shape and volume of the unit cell, were allowed to vary until the residual forces per atom were less than 0.005eV/Å. The resulting equilibrium unit cell was subsequently expanded and compressed uniformly around the equilibrium volume, while keeping the shape of the unit cell fixed. The equilibrium lattice constants were calculated by fitting the total energy per cell as a function of volume using the Birch-Murnaghan^{37}^{,38} equation of state:

where *E* is the total energy per cell, *E*_{0} the equilibrium total energy per cell, *B*_{0} the equilibrium bulk modulus, *V* the unit cell volume, *V*_{0} the equilibrium unit cell volume and the first derivative of the bulk modulus with respect to *V*. In this way, we determine *V*_{0} (from which the equilibrium lattice constant *a* was extracted), *B*_{0} and . It should be noted that we are extending the Murnaghan fit to 2D materials, and quantities such as the bulk modulus should be regarded as fitting parameters rather than physical quantities as discussed by Behera and Mokhopadhyay [BM]^{39}. BM simulated the 2D-hexagonal structure of graphene and silicene using 3D-hexagonal supercells with large values of the lattice parameter c to keep the interlayer interaction negligibly small. They calculated for fixed values of *a* the values of *c* and the ground state energy *E*_{0} for various cell volumes *V*, corresponding to different in-plane lattice constants *a*. Then, by fitting *E*_{0} as a function of *V* with the Birch-Murnaghan equation of state, they extracted *a* from the value of *V* at the minimum of *E*_{0}. Finally, the value of the lattice constant *a* corresponding to *c* going to infinity, was obtained by a linear fit of the data set (*a*, 1/*c*). Here, we have followed a similar procedure.

We present ground-state structural and electronic properties of a series of free-standing monolayer (ML) materials, which are: graphene, silicene, germanene and phosphorene, TMD monolayers MX_{2} in the semiconducting 2H phase^{40} (M=Mo, W; X=S, Se, Te), and one quintuple layer (QL) film of Bi_{2}Se_{3}. The crystal structures are depicted in Fig. 1. We tested how SCAN performs compared to the LDA and PBE-GGA by calculating the lattice constants *a*, the nearest-atom bond lengths *d* for graphene, silicene, germanene and phosphorene, the buckling heights Δ for silicene, germanene and phosphorene, and X-M distances *d*_{M−X} for the TMD monolayers (Fig. 2). These parameters are defined in Fig. 1, and their values are given in Tables 1 and and2.2. For the TMD monolayers, we also extracted the band gaps *E*_{g}, as well as the spin-splittings at the *K* point of the conduction band Δ*E*_{CB}, and the valence band Δ*E*_{VB}, as defined in Fig. 3.

Comparison between the calculated structural parameters using different exchange-correlation functionals: (**a**) and (**b**) Lattice constants, *a*. (**c**) and (**d**) Nearest-atom bond lengths, *d*, and transition metal-chalcogen distances, *d*_{M−X}. The experimental **...**

(**a**) WTe_{2} monolayer band structure along the high-symmetry lines *M-K*′*-*Γ*-K-M* in the Brillouin zone. The colored dots denote spin polarization: blue is for spin *up*, and red is for spin *down*. (**b**) Energy band gap values for the MX_{2} monolayers, **...**

Table 1 lists values of lattice constants *a*, bulk moduli *B*_{0} and their first derivatives , the last two being fitting parameters as we discussed in the Computational Details section above. Trends in the lattice constants are visualized in frames (a) and (b) of Fig. 2. The LDA is seen to underestimate *a*, in agreement with the expectation that it leads to overbinding in solids^{41}. On the other hand, the GGA overcorrects *a*, especially for heavier elements as seen by comparing germanene with silicene and graphene in Fig. 1(a), a behavior observed in 3D metals more generally^{24}. Figure 2(a) and (b) show that the SCAN meta-GGA values lie between the LDA and GGA predictions, suggesting that SCAN meta-GGA cures the overcorrection of the GGA, and generally yields better agreement with experiment, within about 0.5%, although experimental data on freestanding silicene, germanene and phosphorene are not currently available. Remarkably, for the QL Bi_{2}Se_{3}, the SCAN-based lattice parameter is also in excellent accord with the experimental value reported by Kou *et al*.^{42}.

The role of electron correlations in graphene remains an open problem. Accurate Quantum Monte Carlo (QMC) simulations suggest that the ground state of graphene is highly nontrivial, with significant contributions from resonating valence bond (RVB) type states^{43}. [RVB effects appear to be important in systems of low dimensionality more generally, such as the Li clusters]^{44}. The fact that SCAN reproduces the experimental lattice constant of graphene quite well thus indicates that SCAN can reasonably capture features of complex ground states in 2D systems.

It is interesting to consider the QMC result for the lattice constant along the armchair direction in phosphorene^{45}. Surprisingly, we find that LDA already overestimates the QMC armchair lattice constant, even though one normally expects overbinding from the LDA. [Note, experimental lattice constants for phosphorene are not currently available]. Furthermore, we find that SCAN also overestimates the QMC result, as seen in Table 1, and performs at the level of the optB88-vdW^{46} functional, see Fig. 2(a). It is not clear to what extent relaxing the fixed-node approximation in QMC might expand the armchair lattice constant in monolayer black phosphorus, and restore the usual paradigm of LDA underestimating lattice constants more generally. We emphasize that when phosphorene layers are coupled, it becomes crucial to include van der Waals corrections. For example, in bulk black phosphorus, SCAN+rvv10 yields lattice constants in close agreement with both experiment and QMC^{47}, and represents a considerable improvement over PBE+vdW. Concerning the phosphorene lattice constant along the zigzag direction, our results in Table 1 show that is fairly insensitive to corrections beyond the LDA.

Table 2 shows that the equilibrium structures assumed by all 2D films considered (other than graphene) are buckled, i.e. exhibit non-zero values of Δ, and that the buckling is amplified in going from the LDA to the GGA. In sharp contrast, SCAN predicts smaller buckling heights for silicene and germanene compared to the LDA. A possible reason for this flattening trend is that SCAN satisfies the non-uniform coordinate scaling constraint^{15}, while the LDA and GGA do not. In phosphorene, since the buckling height is much larger than that in silicene and germanene, and lies at the scale of a typical chemical bond, SCAN predicts a value comparable to LDA and GGA. Turning to bond lengths, here also we see that, like the lattice constants, SCAN systematically rectifies GGA’s tendency to overcorrect LDA, see Table 2 and Fig. 2(c) and (d). For the TMD films trends in bond lengths between the LDA, SCAN and GGA are similar. Notably, spin-orbit effects, which are included in the calculations, do not seem to influence the trends in bond lengths.

Given the interest in potential applications of TMD films^{11}, Fig. 3(a) shows the band structure of a WTe_{2} monolayer, which is typical of the family of TMD monolayers considered. Table 3 gives the band gaps obtained from the band structures based on different functionals computed *at the equilibrium crystal structures*, see also Fig. 3(b). Note that our band structures arise from a ground-state theory^{41}^{,48}, and thus do not accurately model the band gaps. Nevertheless, the LDA is well-known to reasonably capture optical energy gaps in many materials. GGA expands the lattice, and it generally worsens the band gap. In contrast, consistent with the findings of Yang *et al*.^{19}, SCAN restores an improved agreement with the experimental band gaps, together with improved lattice structures. This good agreement can be understood to be a result of using the generalized Kohn-Sham theory^{48} within SCAN meta-GGA. Incidentally, within the many-body body perturbation theory, Qiu *et al*.^{49} have noticed an interesting compensation between the quasiparticle (QP) and excitonic corrections in the case of transition metal dichalcogenides. For example, in MoS_{2}, the GW approximation yields a direct gap of 2.67eV. The observed optical gap is about 0.8eV smaller, which could be explained as the exciton binding energy.

Returning to Fig. 3(a), note that the conduction band (CB) and the valence band (VB) are split at the *K*-point, which is a consequence of spin-orbit coupling^{50}^{,51}. Furthermore, because TMD monolayers lack inversion symmetry, there is an inversion in the spin-resolved band structures near the Fermi level between the *K* and *K*′ symmetry points (Fig. 3(a)), where blue dots denote spin *up* and red dots spin *down*. These features of band structures of TMD monolayers have been predicted in earlier DFT calculations^{52} and observed in experiments^{53}^{,54}^{,55}^{,56}^{,57}. We define the CB and VB spin-splitting energies as: , and . The values of these splitting energies are listed in Table 3. We see that Δ*E*_{CB}<0 for MoX_{2} monolayers, and Δ*E*_{CB}>0 for the WX_{2} counterparts. This sign change can be explained in terms of the material-dependent spin-orbit coupling effects^{52}. The results of Table 3 indicate that SCAN predicts the spin-splittings in TMD monolayers, at least in some cases (MoS_{2}, MoSe_{2} and WS_{2}), more accurately than the LDA and GGA (see Fig. 3(d)). We thus adduce that SCAN reasonably describes the delicate balance between the exchange, correlation and spin-orbit coupling interactions, which underlie spin-resolved band structures. The exquisite ability of SCAN to capture such subtle effects will allow the study of controlled magnetism in 2D crystals. Interesting proposals have been put forward for monolayer transition metal dichalcogenides^{58}, but magnetic order has not been proven so far in experiments. SCAN meta-GGA could thus accelerate the discovery of these fascinating materials.

In order to test the efficacy of the recently proposed SCAN functional toward capturing improved ground-state properties of layered materials, we have carried out SCAN based computations on monolayers of graphene and a number of ‘beyond graphene’ compounds, including films of transition-metal dichalcogenides (TMDs). The results are compared and contrasted with those based on the commonly used LDA and GGA schemes. SCAN is shown to yield systematic improvements in the equilibrium lattice constants and the nearest-atom bond lengths. We also consider band gaps and spin-splittings in the TMD films, and show that here also the SCAN functional leads to improvements, difficulties of interpreting band gaps in a ground-state computation notwithstanding. We thus conclude that SCAN would provide an improved description of the ground-state electronic and geometric structures of layered materials more generally, at a cost comparable to the LDA and GGA.

**How to cite this article:** Buda, I. G. *et al*. Characterization of Thin Film Materials using SCAN meta-GGA, an Accurate Nonempirical Density Functional. *Sci. Rep.*
**7**, 44766; doi: 10.1038/srep44766 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

It is a pleasure to thank Professor J. P. Perdew for important discussions. This work was supported by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences grant number DE-FG02-07ER46352 (core research), and benefited from Northeastern University’s Advanced Scientific Computation Center (ASCC), the NERSC supercomputing center through DOE grant number DE-AC02-05CH11231 and the support (applications to layered materials) from the DOE EFRC: Center for the Computational Design of Functional Layered Materials (CCDM) under grant number DE-SC0012575.

The authors declare no competing financial interests.

**Author Contributions** B.B., J.S. and C.L. conceived and initiated the study. J.S. and A.R. provided the SCAN meta-GGA functional. I.G.B. performed the calculations and analysis. I.G.B., C.L., B.B. and A.B. prepared the manuscript. All authors contributed to the discussions and reviewed the manuscript.

- Novoselov K. S. et al. . Two-dimensional atomic crystals. Proceedings of the National Academy of Sciences of the United States of America 102, 10451–10453 (2005). [PubMed]
- Castro Neto A. H., Guinea F., Peres N. M. R., Novoselov K. S. & Geim A. K. The electronic properties of graphene. Rev. Mod. Phys. 81, 109–162 (2009).
- Bansil A., Lin H. & Das T. Colloquium: Topological band theory. Rev. Mod. Phys. 88, 021004 (2016).
- Jose D. & Datta A. Structures and chemical properties of silicene: Unlike graphene. Accounts of Chemical Research 47, 593–602 (2014). [PubMed]
- Dávila M. E., Xian L., Cahangirov S., Rubio A. & Lay G. L. Germanene: a novel two-dimensional germanium allotrope akin to graphene and silicene. New Journal of Physics 16, 095002 (2014).
- Zhu F.-f. et al. . Epitaxial growth of two-dimensional stanene. Nat Mater 14, 1020–1025 (2015). [PubMed]
- Tsai W.-F. et al. . Gated silicene as a tunable source of nearly 100% spin-polarized electrons. Nat Commun 4, 1500 (2013). [PubMed]
- Wei Q. & Peng X. Superior mechanical flexibility of phosphorene and few-layer black phosphorus. Applied Physics Letters 104 (2014).
- Wang Q. H., Kalantar-Zadeh K., Kis A., Coleman J. N. & Strano M. S. Electronics and optoelectronics of two-dimensional transition metal dichalcogenides. Nat Nano 7, 699–712 (2012). [PubMed]
- Yu X., Prevot M. S., Guijarro N. & Sivula K.
Self-assembled 2D Wse
_{2}thin films for photoelectrochemical hydrogen production. Nat Commun 6 (2015). [PMC free article] [PubMed] - Singh A. K., Mathew K., Zhuang H. L. & Hennig R. G. Computational screening of 2D materials for photocatalysis. The Journal of Physical Chemistry Letters 6, 1087–1098 (2015). [PubMed]
- Hohenberg P. & Kohn W. Inhomogeneous electron gas. Phys. Rev. 136, B864–B871 (1964).
- Zangwill A. A half century of density functional theory. Physics today 68, 34 (2015).
- Perdew J. P. & Ruzsinszky A. Fourteen easy lessons in density functional theory. International Journal of Quantum Chemistry 110, 2801–2807 (2010).
- Sun J., Ruzsinszky A. & Perdew J. P. Strongly constrained and appropriately normed semilocal density functional. Phys. Rev. Lett. 115, 036402 (2015). [PubMed]
- Sun J. et al. . Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional. Nature Chemistry (2016). [PubMed]
- Kitchaev D. A.
et al. . Energetics of MnO
_{2}polymorphs in density functional theory. Phys. Rev. B 93, 045132 (2016). - Thenuwara A. C. et al. . Copper-intercalated birnessite as a water oxidation catalyst. Langmuir 31, 12807–12813 (2015). [PubMed]
- Yang Z.-h., Peng H., Sun J. & Perdew J. P. More realistic band gaps from meta-generalized gradient approximations: Only in a generalized Kohn-Sham scheme. Phys. Rev. B 93, 205205 (2016).
- Perdew J. P. & Schmidt K. Jacobs ladder of density functional approximations for the exchange correlation energy. AIP Conf. Proc. 577 (2001).
- Kohn W. & Sham L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133–A1138 (1965).
- Vosko S. H., Wilk L. & Nusair M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Canadian Journal of physics 58, 1200–1211 (1980).
- Perdew J. P., Burke K. & Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996). [PubMed]
- Barbiellini B., Moroni E. & Jarlborg T. Effects of gradient corrections on electronic structure in metals. Journal of Physics: Condensed Matter 2, 7597 (1990).
- Tao J., Perdew J. P., Staroverov V. N. & Scuseria G. E. Climbing the density functional ladder: Nonempirical meta–generalized gradient approximation designed for molecules and solids. Physical Review Letters 91, 146401 (2003). [PubMed]
- Sun J., Xiao B. & Ruzsinszky A. Communication: Effect of the orbital-overlap dependence in the meta generalized gradient approximation. The Journal of Chemical Physics 137 (2012). [PubMed]
- Heyd J., Scuseria G. E. & Ernzerhof M. Hybrid functionals based on a screened coulomb potential. The Journal of Chemical Physics 118, 8207–8215 (2003).
- Ren X., Rinke P., Joas C. & Scheffler M. Random-phase approximation and its applications in computational chemistry and materials science. Journal of Materials Science 47, 7447–7471 (2012).
- Gunnarsson O. & Lundqvist B. I. Exchange and correlation in atoms, molecules, and solids by the spin-density-functional formalism. Phys. Rev. B 13, 4274–4298 (1976).
- Perdew J. P., Ruzsinszky A., Csonka G. I., Constantin L. A. & Sun J. Workhorse semilocal density functional for condensed matter physics and quantum chemistry. Phys. Rev. Lett. 103, 026403 (2009). [PubMed]
- Batista E. R. et al. . Comparison of screened hybrid density functional theory to diffusion Monte Carlo in calculations of total energies of silicon phases and defects. Phys. Rev. B 74, 121102 (2006).
- Xiao B.
et al. . Testing density functionals for structural phase transitions of solids under pressure: Si, SiO
_{2}, and Zr. Phys. Rev. B 88, 184103 (2013). - Kresse G. & Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758–1775 (1999).
- Kresse G. & Furthmüller J.
Efficient iterative schemes for
*ab initio*total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996). [PubMed] - Kresse G. & Hafner J.
*Ab initio*molecular dynamics for open-shell transition metals. Phys. Rev. B 48, 13115–13118 (1993). [PubMed] - Perdew J. P. & Zunger A. Self-interaction correction to density-functional approximations for many-electron systems. Phys. Rev. B 23, 5048–5079 (1981).
- Murnaghan F. The compressibility of media under extreme pressures. Proc. Natl. Acad. Sci. USA 30, 244–247 (1944). [PubMed]
- Birch F. Finite elastic strain of cubic crystals. Phys. Rev. 71, 809–824 (1947).
- Behera H. & Mukhopadhyay G. Structural and electronic properties of graphene and silicene: An FP-(L)APW+lo study. AIP Conference Proceedings 1313, 152–155 (2010).
- Kappera R.
et al. . Phase-engineered low-resistance contacts for ultrathin MoS
_{2}transistors. Nature materials 13, 1128–1134 (2014). [PubMed] - Jones R. O. & Gunnarsson O. The density functional formalism, its applications and prospects. Rev. Mod. Phys. 61, 689–746 (1989).
- Kou L. et al. . Graphene-based topological insulator with an intrinsic bulk band gap above room temperature. Nano Letters 13, 6251–6255 (2013). [PubMed]
- Marchi M., Azadi S. & Sorella S. Fate of the resonating valence bond in graphene. Phys. Rev. Lett. 107, 086807 (2011). [PubMed]
- Nissenbaum D., Spanu L., Attaccalite C., Barbiellini B. & Bansil A. Resonating-valence-bond ground state of lithium nanoclusters. Phys. Rev. B 79, 035416 (2009).
- Shulenburger L., Baczewski A., Zhu Z., Guan J. & Tomk D. The nature of the interlayer interaction in bulk and few-layer phosphorus. Nano Letters 15, 8170–8175 (2015). [PubMed]
- Qiao J., Kong X., Hu Z.-X., Yang F. & Ji W. High-mobility transport anisotropy and linear dichroism in few-layer black phosphorus. Nature communications 5 (2014). [PMC free article] [PubMed]
- Peng H., Yang Z.-H., Perdew J. P. & Sun J. Versatile van der Waals density functional based on a meta-generalized gradient approximation. Phys. Rev. X 6, 041005 (2016).
- Perdew J. P.
et al. . Understanding Band Gaps of Solids in Generalized Kohn-Sham Theory.
*arXiv:cond-mat.mtrl-sci 1608.06715v1*(2016). [PubMed] - Qiu D. Y., Cao T. & Louie S. G. Nonanalyticity, valley quantum phases, and lightlike exciton dispersion in monolayer transition metal dichalcogenides: Theory and first-principles calculations. Phys. Rev. Lett. 115, 176801 (2015). [PubMed]
- Chang T.-R., Lin H., Jeng H.-T. & Bansil A.
Thickness dependence of spin polarization and electronic structure of ultra-thin films of MoS
_{2}and related transition-metal dichalcogenides. Scientific reports 4, 6270 (2014). [PMC free article] [PubMed] - Alidoust N.
et al. . Observation of monolayer valence band spin-orbit effect and induced quantum well states in MoX
_{2}. Nat Commun 5 (2014). [PubMed] - Kośmider K., González J. W. & Fernández-Rossier J. Large spin splitting in the conduction band of transition metal dichalcogenide monolayers. Phys. Rev. B 88, 245436 (2013).
- Kormányos A. et al. . k. p theory for two-dimensional transition metal dichalcogenide semiconductors. 2D Materials 2, 022001 (2015).
- Zhang Y.
et al. . Direct observation of the transition from indirect to direct bandgap in atomically thin epitaxial
*MoSe*_{2}. Nat Nano 9, 111–115 (2014). [PubMed] - Böker T.
et al. . Band structure of MoS
_{2}, MoSe_{2}, and*α*-MoTe_{2}angle-resolved photoelectron spectroscopy and*ab initio*calculations. Phys. Rev. B 64, 235305 (2001). - Dendzik M.
et al. . Growth and electronic structure of epitaxial single-layer WS
_{2}on Au(111). Phys. Rev. B 92, 245442 (2015). - Le D.
et al. . Spinorbit coupling in the band structure of monolayer Wse
_{2}. Journal of Physics: Condensed Matter 27, 182201 (2015). [PubMed] - Cheng Y. C., Zhu Z. Y., Mi W. B., Guo Z. B. & Schwingenschlögl U.
Prediction of two-dimensional diluted magnetic semiconductors: Doped monolayer MoS
_{2}systems. Phys. Rev. B 87, 100401 (2013). - Baskin Y. & Meyer L. Lattice constants of graphite at low temperatures. Phys. Rev. 100, 544–544 (1955).
- Brown A. & Rundqvist S. Refinement of the crystal structure of black phosphorus. Acta Crystallographica 19, 684–685 (1965).
- Jiang Y. C., Gao J. & Wang L.
Raman fingerprint for semi-metal Wte
_{2}evolving from bulk to monolayer. Scientific Reports 6 (2016). [PMC free article] [PubMed] - Liu G.-B., Xiao D., Yao Y., Xu X. & Yao W.
Electronic structures and theoretical modelling of two-dimensional group-
*VIB*transition metal dichalcogenides. Chem. Soc. Rev. 44, 2643–2663 (2015). [PubMed] - Froehlicher G., Lorchat E. & Berciaud S.
Direct versus indirect band gap emission and exciton-exciton annihilation in atomically thin molybdenum ditelluride MoTe
_{2}. Phys. Rev. B 94, 085429 (2016). - Miwa J. A.
et al. . Electronic structure of epitaxial single-layer MoS
_{2}. Phys. Rev. Lett. 114, 046802 (2015). [PubMed]

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |