Rational design of transparent p-type conducting non-oxide materials from high-throughput calculations

In this work, high-throughput ab initio calculations are employed to identify the most promising chalcogenide-based semiconductors for p-type transparent conducting materials (TCMs). A large computational data set is investigated by data mining. Binary semiconductors with large band gaps (Eg) and anions that are less electronegative than oxygen are considered. The roles of intrinsic defects and extrinsic dopants are investigated to probe the p-type performance of these semiconductors. Nine novel p-type non-oxide TCMs that have a low hole effective mass, good optical transparency, and hole dopability are proposed (ZnS, ZnSe, ZnTe, MgS, MgTe, GaSe, GaTe, Al2Se3, and BeTe). This study also focuses on a material engineering approach to modulate the electronic properties as a function of the layer thickness and external stress.


Introduction
2][3][4] Potential materials for TCMs should offer two contradictory properties: a wide band gap and a low carrier effective mass. 5Among various TCMs, oxides have been widely studied owing to their stability, low electron effective mass, and optical transparency. 6Previous research on transparent conducting oxides (TCOs) has shown their outstanding properties as high-performance n-type conductors in transparent transistors, UV light-emitting diodes, detectors, gas sensors, and solar cells. 7However, in many optoelectronic applications, such as interlayer electrodes in tandem cells, high-performance p-type TCMs are required.The performance of p-type TCOs is hindered by the localized p states of oxygen in the valence band, which results in heavy holes. 8,9Demands for p-type TCMs and the deficient performance of this class of material have promoted intensive research efforts from academia and industry to rationalize and design potential p-type TCMs. 4,10he discovery of CuAlO 2 as a promising p-type TCO led to the development of rational design rules for p-type TCOs containing Cu 1+ and prompted an extensive investigation of Cu-based materials. 11Experimental and theoretical efforts also revealed that mixing O with other anions, such as S and Se, influences p-type dopability.For example, p-type dopability has been observed in LaCuOS and BaCuOSe based compounds. 12urther, mixing of other anions with oxygen leads to higher valence band energies, making ternary oxides promising p-type conductors. 8,13,14In addition, p-type dopability has been tested for a series of compounds in which Cu is replaced by other cations (X 2 SeO 2 (X = La, Pr, Gd, and Nd)). 15However, p-type TCOs are commonly hampered by the carrier mobility of candidate materials. 16Hence, there is a necessity for enhanced research efforts towards the rational design of suitable materials beyond TCOs.
Experimental studies and theoretical calculations have also been performed on non-oxide TCMs to determine the underlying physics of p-type TCMs.8][19] Further, it has been observed that the electronegativity and size of the anion influence the bonding characteristics and structures of such materials, which determine the p-type performance. 20Studies of the intrinsic and substitutional defects in BaCuChF (Ch = S, Se, Te) based compounds have revealed unintentional p-type conductivity owing to the presence of Cu vacancies. 21,22The optoelectronic properties of BAMChF (A = Sr, Ba; M = Cu, Ag; and Ch = S, Se, Te) based chalcogenide fluorides have been found to vary with atomic substitution or with vacancy formation. 23Chalcogenide-based compounds have been found to exhibit a larger carrier mobility than oxides. 24,25he majority of experimental studies and ab initio calculations performed so far have analyzed the p-type performance of Ba-and Cu-based ternary and quaternary chalcogenides and oxides.Until now, the question remains open as to whether alternative chemistries and design rules could lead to the discovery of promising transparent materials with low hole effective masses.Addressing this question plays a dominant role in the TCM field, with the aim of rationally designing materials sought by the optoelectronic industry.Recently, high-throughput density functional theory (DFT) approach has been employed for non-oxide binary semiconductors and found boron phosphide (BP) to be a promising p-type TCM candidate. 26Their work has opened up a new direction of investigation into binary non-oxide TCMs.However, an in-depth understanding of the factors that influence the p-type performance of these materials is still required.
In this work, we performed high-throughput DFT calculations for chalcogenide-based binary semiconductors to identify the most promising p-type transparent conductors.Our material screening strategy led to the identification of nine potential p-type transparent conductors (ZnS, ZnSe, ZnTe, MgS, MgTe, GaSe, GaTe, Al 2 Se 3 , and BeTe).In addition, we proposed strategies to alter the electronic structure of the materials through material engineering approach.As an example, we demonstrated how to tune the electronic properties of GeS so that it can be employed as a p-type TCM.

DFT calculations
All DFT calculations were performed using the Vienna Ab initio Simulation Package (VASP) 27 with PAW pseudopotentials 28,29 and a plane-wave cutoff of 500 eV.The exchange correlation is described with the Perdew-Burke-Ernzerhof (PBE) 30 form of the generalized gradient approximation for structural optimization.All the atoms in the system were allowed to relax until the force on each atom was less than 0.01 eV Å À1 .Brillouin zone integration was performed on a Monkhorst k-point mesh.We adapted the k-mesh of each supercell to maintain a constant k-point density for each system.Since PBE functional tends to underestimate the band gaps of semiconductors, we employed the HSE06 screened hybrid functional (HSE) 31 to calculate the band structures and hole effective masses of selected compounds.Effective masses were calculated using the EMC code. 32We first determined the top of the valence band for each compound.Then, the EMC code has been employed to evaluate q 2 E/qk i k j (i, j = x, y, and z) numerically using the finite difference method.The effective mass tensor is calculated from this derivative.For non-cubic systems we calculated the band structure along the high-symmetry directions as well.Assuming the valence band has a parabolic shape in the small area close to the VBM, we calculated the aforementioned derivative by curve fitting.We note in passing that we compared our calculated effective masses with those calculated by BoltzTraP code 33 and found a reasonable agreement although the methods are different.
The supercell calculation is the most common method to calculate the formation energy of defects. 34In this method a defect a with the charge state q is introduced into the supercell.The periodic boundary is held by repeating the supercell periodically in x, y, and z directions.The dimension of the supercell has to be rather large to reduce the spurious interactions between a defect and its periodic replicas. 35In this work defect formation energies and the corresponding charge transition states were calculated on large supercells (more than 100 atoms) in order to minimize finite-size effects.The formation energy of a defect is defined as where E tot.[a q ]is the total energy of the supercell with the defect a in the charge state q and E bulk tot. is the total energy of the equivalent defect-free bulk supercell.n i represents the number of atoms of type i that is added to/removed from the supercell to create defects.m i + Dm i are associated chemical potential of elements in elemental phase (m i ) subjected to the thermodynamic limits (Dm i ).The thermodynamic limits to the chemical potentials were computed by determining the stability region for all competing structures with respect to the reference structure (For more details see the Tables S1  and S2 of the ESI †).E VBM is the valence band maximum of the defect-free bulk and m e is the chemical potential of electrons.The thermodynamic charge transition level is defined as the chemical potential of electrons at which the charge state of defect a changes from q to q 0 .The defect formation energies were corrected according to the average electrostatic potential as well as the interaction of charged defects. 36,37To do so, first the VBMs were aligned by q(D v0/b ), where D v0/b is the difference between the average electrostatic potential of the bulk and the supercell containing the associated neutral defect.Then, the formation energies were corrected by E q corr.that accounts for the interaction of a charged defect with its periodic replicas. 37n order to ensure the validity of our approach and results, we performed calculations for our final compounds with the recently developed PyCDT script. 38This script gets the atomic structure and formation energies of the compounds from the Materials Project (MP) database. 39The correction terms are calculated based on the Freysoldt approach. 40As it is expected, 36 there is a very good agreement between the results achieved by two approaches.

Material screening and engineering strategy
High-throughput computational screening was performed to identify potential transparent conductors.Further, the electronic and optical properties of selected compounds were tuned by material engineering.Our approach utilized data mining of chalcogenide-based non-oxide binary semiconductors that have already been synthesized and documented in the literature. 41he crystalline structures of the binary compounds were obtained from the Inorganic Crystal Structure Database (ICSD) or Materials Project database. 39Our screening strategy is based on the flowchart depicted in Fig. 1.The following steps were followed for the rational design of materials that are potential p-type candidates: Steps 1 and 2. Database screening of chalcogenide-based binary non-oxide semiconductors was performed for more than 500 compounds.Our aim was to find nontoxic compounds with E g 4 1.7 eV.Some layered chalcogenides that have E g values larger than 0.8 eV were further selected for band engineering studies.We identified less than 40 binary chalcogenides that satisfy E g 4 1.7 eV.Our compounds consist of series of cations (Zn, Mg, Sn, In, Ga, Al, Na, K, Ge, Ca, Be, and Ba) and anions (S, Se, and Te).
Step 3. The atomic structures of the selected compounds were optimized and their stabilities were calculated using DFT.The chemical potential diagram corresponding to each compound was also computed.
Step 4. Materials with a hole effective mass o 1 and E g 4 1.7 eV were further selected to understand the defect physics.
Role of intrinsic defects.Intrinsic defects such as vacancies and antisites were considered to determine the p-type dopability.
Role of extrinsic defects.Materials that upheld intrinsic defects were further selected to understand foreign atom chemistry.Extrinsic dopants, such as various cationic and anionic dopants, were substituted to determine their effects on p-type performance.
Band engineering.Materials with E g o 1.7 eV or a hole effective mass 41 were selected for band engineering.3][44] If the band gap size increases or the size of the effective mass decreases by band engineering then the role of intrinsic and extrinsic defects on p-type performance was explored.

Results and discussion
Our study focused on material screening and material engineering for optoelectronic applications.Our criteria for material screening were as follows: materials with E g 4 1.7 eV (the minimum value needed for Cu(In, Ga)Se 2 -based tandem cells 45 ) and a hole effective mass o1.Among 500 compounds, we found less than 40 chalcogenides that satisfied our screening criterion for E g .The calculated hole effective mass versus the experimentally measured band gap for some compounds with hole effective masses of o3 are shown in Fig. 2. We divided the compounds into three groups: group A compounds satisfy our screening criteria for E g and hole effective mass.The compounds to the left of the vertical green line (group B) have low hole effective masses but E g is too small whereas the compounds above the horizontal green line (group C) have large E g value but their hole effective masses are larger than 1.Material engineering was carried out to overcome the deficiencies in the group B and C compounds.
We further examined the dopability of the group A compounds by introducing intrinsic defects such as vacancies and antisites.Anion vacancies are usually donor defects that compensate for the holes provided by acceptor-like defects.The p-type dopability of a material depends predominantly on the formation of anion vacancies.It is possible to hinder the formation of anion vacancies during synthesis by controlling the chemical potential of the elements.Anion-rich conditions impede the formation of anion vacancies, but can also lead to the formation of antisites.Hence, we calculated the intrinsic defect formation energies under anion-rich conditions.Among the various chalcogenides, only a few compounds were found to uphold the intrinsic defects.For example, cubic CaTe (Fm% 3m) has a wide E g (B4 eV) 46 and a relatively small hole effective mass (0.72).Despite its promising electronic properties, the formation of Te Ca antisites in CaTe compensated for the holes generated by V Ca .We neglected compounds from group A whose hole-compensating defects have low formation energies (red circles in Fig. 2).
The materials we considered here demonstrated propensity for increased anion vacancy formation energy and decreased cation vacancy formation energy on moving from sulfides to tellurides.For example, in ZnTe, the formation of anion vacancies is about 1.2 eV larger than that in ZnS.Further, cation vacancies are more favorable in ZnTe than in ZnS or ZnSe.Similar trends were observed for antisite defects, although the change in formation energies is much smaller.From these trends, we cannot say a priori whether tellurides are better candidates for p-type TCMs than sulfides.While lower anion vacancies in tellurides are beneficial for hole conductivity, the formation of donor-like antisites can be detrimental for p-type conductivity.
Another vital prerequisite for high-performance TCMs is plausibility of doping with dopants that enhance the freecarrier density in the system.Regardless of the possibility that the intrinsic defects will increase the hole concentration, a TCM must be doped to improve performance.Therefore, we doped our compounds that showed good p-type dopability.We also studied the p-type doping properties of compounds for which the formation energy of intrinsic defects is high, even if they did not show the desired properties.A potential p-type TCM candidate may also exhibit considerably lower formation energy for some p-type dopants compared with intrinsic defects.
When dopants are incorporated into the host semiconductor, depending on experimental conditions, the chemical potentials of the dopant are subject to specific bounds. 47We considered two extreme limits in our calculations, namely when the dopants are in the elemental phases, and when the chemical potentials of dopants are bounded by avoiding the formation of the secondary phases.The upper limits for the chemical potential of dopants to avoid the formation of the secondary phases are listed in the ESI.† In practice, the chemical potentials of dopants vary between these two limits depending on the experimental conditions.In Fig. 3, the formation energies of extrinsic defects under the later condition are shown.The formation energies under the former condition are shown in the ESI.† We found nine chalcogenides that satisfied the abovementioned criteria and could be employed as p-type TCMs (blue squares in Fig. 2).The compounds, their symmetries, calculated/measured band gaps, and calculated hole effective masses are listed in Table 1.We note that these compounds have different atomic and electronic structures.The required application determines whether a layered-structure compound (GaSe) or a symmetric structure (ZnTe) is desirable.The following sections outline the properties of the most promising p-type TCM candidates.

Zn-Based chalcogenides
Zinc chalcogenides have zincblende atomic structure (F% 43m).The electronic structures of ZnS, ZnSe, and ZnTe are very similar, namely, the top of the valence band is at the G point and has mainly anion p character with a small contribution from Zn.The band gaps and hole effective masses of Zn chalcogenides decrease with decreasing anion electronegativity (Fig. 2).
ZnS is a direct band gap semiconductor with a wide measured E g of 3.7 eV. 48When the Fermi level was close to the valence band maximum (VBM), hole-compensating defects formed more easily than cation vacancies (Fig. 3a).The formation energy of cation vacancies, however, decreased rapidly as the electron chemical potential m e increased.Considering the intrinsic defects, ZnS is not a promising p-type conductor. 55owever, the role of extrinsic dopants for ZnS was explored owing to its promising electronic properties.We observed that extrinsic dopants such as Na and K for Zn, as well as N for S, increase the hole concentration.The formation energies for these defects were found to be smaller than those for intrinsic defects.It has been reported that the transparency and conductivity of Cu-doped ZnS are satisfactory for a p-type TCM. 56,57n agreement with the previous results, we also found that Cu Zn acts as an acceptor and increases the hole concentration.Doping with the suggested alkali metals will give rise to better p-type conductivity by preventing the formation of intrinsic antisite defects.N S defects, on the other hand, will prevent the formation of hole-compensating anion vacancies.
In ZnTe, cation vacancies formed more easily than anion vacancies and antisites (Fig. 3b).Cation vacancies introduced holes into the system and improved conductivity in ZnTe.Antisites and anion vacancies impacted the transport properties of ZnTe.However, the formation energies of these intrinsic defects were too high.Similar to ZnS, ZnTe could also be doped with Na, Cu, and N. Na Zn , Cu Zn , and N Te substitutional defects had low formation energies and formed easily in ZnTe.Cu-doped ZnTe films and nanowires have been synthesized previously and show good p-type conductivity and transparency. 58,59In this respect, our results are in agreement with experimental observations.The behavior of ZnSe falls between that of ZnS and ZnTe (ESI, † Fig. S2).There are several reports on the good conductivity of p-doped ZnSe 60,61 and our results agree well with the reported observations.We found Zn-based chalcogenides to be promising p-type candidates.

Mg-Based chalcogenides
Mg-Based chalcogenides have wide band gaps, which is the primary requirement for optoelectronic applications.Cubic MgS has an indirect E g of 6.4 eV, 51 whereas hexagonal MgTe has a direct band gap of 4.7 eV 52 (Fig. S5 of the ESI †).In order to understand defect physics in these compounds, we examined the influence of point defects on the carrier concentration.The formation of intrinsic defects was less likely in MgS than in MgTe (Fig. 3c and d).Although both materials had small hole densities without doping, Na-or K-doped MgS and MgTe are good candidates for p-type conductors.Adding these alkali metals to these systems led to the formation of Na Mg and K Mg substitutional acceptor defects.Anion substitutional defects (N S and N Te ) also resulted in acceptor defects in these compounds.
The VBM of MgS consists of p states of the anion with a hole effective mass of 0.96.By substituting S with Te, the characteristics of the VBM did not change, but the effective mass and E g decreased.We also studied MgS with the zinc-blende structure (F% 43m).In contrast to MgS with the Fm% 3m symmetry, zinc-blende MgS is not a good candidate for p-type TCMs owing to easy formation of donor-like defects.
Our results indicate that with intentional doping, MgS and MgTe could be used as a high-performance TCMs.To the best of our knowledge, Mg-based chalcogenides have not been considered so far as p-type TCMs.In contrast to MgS and similar to other layered-structure materials, the effective mass of MgTe is anisotropic.While the hole effective mass is as small as 0.18 parallel to the layers, it is 8.06 perpendicular to the layers.It has been shown that depending on the growth conditions, the orientation of the layers in layered-structure materials can be controlled. 62This property gives a degree of freedom in designing devices to control the conductivity in the preferred direction.

Al-Based chalcogenides
Al 2 Se 3 has a wide direct band gap (E g = 3.1 eV) 41 (Fig. S6 of the ESI †).The top of the valence band has mainly Se p character and is very close to the G point.This band structure results in a rather low hole effective mass of 0.30.The formation of anion vacancies generated transition states at about 0.60 eV above the VBM (Fig. 3e).Antisite and cation vacancies had low formation energies compared with anion vacancies.However, the intrinsic defects in Al 2 Se 3 have very high formation energies and consequently their concentrations are low.In contrast to intrinsic defects, N Se and Mg Al substitutional defects had low formation energies and both defects were beneficial for p-type conductivity.P Se was also an acceptor dopant, but the formation energy of this defect was comparable to those of intrinsic defects (not shown here).Al 2 Se 3 has a wurtzite-type structure with four ordered Al vacancies per unit cell.Owing to this symmetry, the effective mass of this compound is anisotropic (0.30, 2.71, and 3.33).This anisotropy is less pronounced in devices where TCMs are grown as polycrystalline films.
Replacing Se with Te resulted in a smaller E g and hole effective mass in Al 2 Te 3 .The formation of antisites was more favorable in Al 2 Te 3 than in Al 2 Se 3 .In contrast to Al 2 Se 3 , cation vacancies created charge transition levels in the band gap of Al 2 Te 3 .Although these defects hindered p dopability, substitutional dopants such as Mg Al and N Te could be employed to improve the hole conductivity of Al 2 Te 3 .Al 2 S 3 has a larger band gap and better p-dopability than Al 2 Se 3 .Nevertheless, we could not find an appropriate dopant to increase the hole concentration in Al 2 S 3 .Hence, among aluminum chalcogenides, Al 2 Se 3 was found to be the most promising candidate.

Be-Based chalcogenides
BeTe is an indirect band gap semiconductor with a large band gap (E g = 2.9 eV) 41 and a small hole effective mass of 0.35.Similar to other good candidates, the VBM consists of p states.Without doping, BeTe is a potential p-type conductor because the formation energies for cation vacancies are low compared with other intrinsic defects (Fig. 3f).To enhance the hole conductivity, BeTe was intentionally doped with Li, Na, and N atoms.Li Be has low formation energies when the Fermi level is at the VBM.This defect is an acceptor and increases the hole concentration in the valence band.BeS has a larger band gap and hole effective mass than BeTe.Antisite defects were detrimental for p-type conductivity in BeS.Doping with alkali metals was also not beneficial since the formation energies were comparable to those of antisites.Hence, BeTe was found to be the most promising candidate among the various Be-based chalcogenides.

Ga-Based chalcogenides
GaSe is another chalcogenide that has a low hole effective mass and relatively large band gap (2.1 eV). 53The VBM of GaSe has p character, which is a mixture of the Ga and Se p states.Anion vacancies and antisites acted as donors at small electron chemical potentials and created transition levels in the band gap (Fig. 3g).Cation vacancies created a transition state in the band gap.All interstitial defects had relatively high formation energies and consequently low concentrations in GaSe.GaSe could be doped with P and Zn to enhance the p-type conductivity.Zn Ga and P Se defects created shallow acceptor levels.It has been shown experimentally that Zn enhances the hole conductivity of GaSe, 63 which is also evident from our calculations.
Replacing Se with S resulted in a wider band gap and larger effective mass.The dopability of GaS is worse than that of GaSe owing to the formation of anion vacancies rather than cation vacancies.Zn Ga and N S substitutional defects in GaS remained neutral and were not beneficial for p-type conductivity.Replacing Se with Te, on the other hand, increased the probability of cation vacancy formation in GaTe.Although cation vacancies are acceptor defects, they are harmful for conductivity since they create charge transition levels in the band gap.To enhance the hole conductivity, GaTe can be doped with Zn and N (Fig. 3h).The size of the band gap and hole effective mass decreases in Ga chalcogenides from GaS to GaTe.Despite the low effective mass and good dopability of GaTe, it is not suitable for many optoelectronic applications owing to its small band gap.For this reason, we investigated GaTe and similar compounds from the band engineering perspective as discussed in the next section.We note in passing that the effective masses of GaSe and GaTe are anisotropic owing to the crystal symmetry.GaSe and GaTe are layered compounds and the effective masses are small parallel to the layers; 0.30 and 0.14 for GaSe and GaTe, respectively.The values of effective mass perpendicular to the layers are 1.78 and 18.92 for GaSe, and 0.91 and 10.01 for GaTe.

Role of the layer thickness
5][66] Here, we investigated the electronic and optical properties of the group B compounds (Fig. 2) as a function of the layer thickness.As shown in Table 2, E g shows an increasing trend with a reduction in the layer thickness.The change in the band gap size could be as large as 60% for SnS.As an example, we show the modulation of the band gap of GeS with the layer thickness (Fig. S3 of the ESI †).Reducing the layers of GeS widened the band gap and at the same time induced a direct to indirect band gap transition.
Changes in the hole effective masses were comparatively predominant in Ge-based compounds (Table S2 of the ESI †).Although Te-based compounds normally have smaller hole effective masses than sulfides, the band gap modulation was less effective in Te-based compounds and most promising in S-based compounds.These results show that anions influence the optoelectronic properties significantly.We doped two S-based compounds (GeS and SnS) that showed large band gap modulation with the layer thickness to understand the defect physics.In GeS, cation vacancies had lower formation energies than anion vacancies and antisites (Fig. S4 of the ESI †).In contrast, for SnS, the anion vacancies and antisites had lower formation energies than cation vacancies, which hindered the performance.Further, extrinsic dopants, such as Na Ge and N S , were found to increase the hole concentration in GeS.These results indicate that by controlling the layer thickness and by substitutional doping, GeS can be considered as a promising p-type TCM.

External stress
In devices, TCMs are normally used as thin films.The lattice parameters of epitaxially grown films are either compressed or expanded owing to lattice mismatch.Normally, for materials under pressure, when the lattice parameters are compressed, the valence band becomes more disperse, which results in a smaller hole effective mass. 44,67,68This change in the electronic structure by stress can be used to enhance hole conductivity.
If TCM films are grown coherently on substrates in such a way that the lattice parameters of the film are compressed, then one expects to observe a smaller hole effective mass.We studied the evolution of the optoelectronic properties of a few compounds in groups A and C (Fig. 2) when the lattice parameters were compressed.The lattice parameters of MgS, CaSe, Na 2 Te, and BaTe were compressed by about 5% and the effective mass was calculated.The changes in the effective masses were about 9% for Na 2 Te and about 15% for MgS, CaSe, and BaTe.We suggest that this strategy can be employed for reducing the hole effective mass further to improve the performance of p-type TCMs.

Concluding remarks
Data mining and high-throughput calculations were performed to identify high-performance p-type TCMs.Our criteria for material selection are a band gap larger than 1.7 eV and a hole effective mass smaller than 1.Our calculations suggested nine promising p-type TCM candidates (ZnS, ZnSe, ZnTe, MgS, MgTe, GaSe, GaTe, Al 2 Se 3 , and BeTe).The VBM for the promising candidates consists of anion p states, which demonstrates the importance of the anions in p-type TCMs.For a particular cation, the anion vacancy formation energies were found to increase from sulfides to tellurides, whereas the formation energies of the cation vacancies and antisites decreased.We showed that with intentional external doping, some of the compounds with high-energy intrinsic defects could be promising p-type TCMs.In addition, material engineering was carried out for materials with E g o 1.7 eV or a hole effective mass larger than 1.It was observed that the optoelectronic properties of a material could be modulated by tuning the layer thickness and external stress.

Fig. 1
Fig. 1 Flowchart for data screening and material engineering.

Fig. 2
Fig. 2 Experimentally measured band gap 41 vs. calculated hole effective mass.Blue squares indicate the materials we suggest for p-type TCM applications.Red circles show materials that are not suitable for highperformance TCMs either owing to a small band gap (group B), or a large hole effective mass (group C), or the lack of dopability.

Table 1
Selected compounds, their symmetries and their calculated band gaps and effective masses.Experimentally measured band gaps are given in parentheses This journal is © The Royal Society of Chemistry 2018

Table 2
Band gaps modulation (HSE calculations) with the thickness for layered compounds