3.1. Crystal structure and stability
The BiOX has a unique layered structure [
40] in which the [Bi
2O
2]
2+ layer is interleaved by double layers of halogen ions [X]
−. In the monolayer [X-Bi-O-Bi-X], each Bi atom is coordinated by four oxygen atoms and four halogen atoms, as shown in
Figure 1 (a) and (b). The optimized lattice constants are calculated to be
a =
b = 3.907 Å and
c = 7.492 Å for bulk BiOCl,
a =
b = 3.943 Å and
c = 8.258 Å for bulk BiOBr, and
a =
b = 4.02 Å and
c = 9.271 Å for bulk BiOI, which are in good agreement with the previous studies in both experiment and theory [
20,
22,
41]. The calculated lattice constants of BiOX monolayer (3.876, 3.926, and 4.015 Å for BiOCl, BiOBr, and BiOI monolayer, respectively) are also consistent with previous theoretical results [
27,
41], verifying the validity in the optimization calculations. We have found that the optimized lattice constants (
a) increase slightly with the number of layers, as shown in
Table 1. Meanwhile, we design a new structure of bi-halogen BiOX (X = Cl
0.5Br
0.5, Cl
0.5I
0.5, Br
0.5I
0.5) monolayer as shown in
Figure 1 (c) - (h), replacing one [X]
− layer of BiOX (X = Cl, Br, I) with another halogen atoms and maintaining the tetragonal structure. The lattice constants of the BiOX (X = Cl, Cl
0.5Br
0.5, Br, Cl
0.5I
0.5, Br
0.5I
0.5, I) monolayer upon full relaxation by PBE-based DFT calculation are listed in Table S1, which increase monotonically with the increase of halogen atomic number.
BiOX possesses a layered crystal structure to offer flexibility and the possibility of exfoliation into ultrathin flakes. This flexibility enables the integration of single-layer BiOX into various optoelectronic device architectures. We first try stressing the cleavage energy
Ecl to access the feasibility of the mechanical exfoliation, as is shown in
Figure 2 as a function of the distance for a BiOX (X = Cl, Br, I) monolayer to be separated from a 5-layered structure. The calculated cleavage energy of monolayer BiOCl, BiOBr and BiOI is 0.455, 0.346 and 0.302 J/m
2, respectively. Referring to those of graphene (0.37 J/m
2), single-layer SnP
3 (0.71 J/m
2) and δ-InP
3 monolayer (0.827 J/m
2) [42-44], BiOX monolayer is a fairly typical 2D material with moderate van der Waals interactions, implying that they can be easily prepared by mechanical exfoliation.
There is no imaginary vibration frequency observed in the phonon spectra of the BiOX monolayer, as shown inFigure S1, indicating the dynamic stability of BiOX. Meanwhile, no structure disruption occurs after 3000 fs AIMD simulations at 300 K (Figure S2), demonstrating that BiOX monolayer are thermodynamically stable. In addition, the cohesive energy was calculated to be 4.92, 4.77, 4.58, 4.84, 4.74, and 4.67 eV/atom for BiOCl, BiOBr, BiOI, BiOCl
0.5Br
0.5, BiOCl
0.5I
0.5, and BiOBr
0.5I
0.5 monolayer, respectively. This implies that BiOCl is the most stable among the above materials. The cohesive energy of the BiOX is greater than that of Bi
2O
2S (3.85 eV/atom) [
21], p-GeN
2 (3.96 eV/atom) [
45] and Silicene (3.71 eV/atom) [
46], indicating their higher stability.
3.2. Thickness-dependent behavior
Because of the Van der Waals interlayer interactions, 2D materials usually possess layer-number-dependent electronic properties [
47,
48]. We investigate the band structure and optical properties of few-layer (two to five layers) BiOX (X = Cl, Br, I) by HSE06 functional [
33]. It is first confirmed that the BiOX (X = Cl, Br, I) monolayer exhibit indirect semiconductor band structure, where the CBM is located at the Г point while the VBM is located between the Г and X points in the irreducible Brillouin zone. As shown in
Figure 3 (a), BiOBr is taken as a representor to display the relationship between the band structure and the thickness (layer number). With increasing layer number, BiOX band structure does not change much, i.e., the CBM always stays at the Г point, while the VBM locates along the Г-X direction. Specifically, the VBM of BiOI moves to and then stay at the X point from 2 single layers on. It is thus clear that BiOX maintains an indirect bandgap character as the layer number increases, which may be helpful for the material design of BiOX. In the meantime, the calculated bandgaps of BiOCl, BiOBr, and BiOI monolayers are 3.745, 3.354, and 2.278 eV, respectively (Table S2), as normally expected. With increasing layer number, the bandgaps do vary as shown in
Figure 3 (b). They decrease monotonically and gradually approach those of their bulks. Those of BiOCl and BiOBr monolayers are reduced obviously with the layer number due to the weak quantum size effect [
41,
49], but that of BiOI shows less sensitivity to the layer number, which may be due to the relatively large charge transfer and orbital hybridization between Bi and I atoms. According to the projected density of states in
Figure 3 (c), the CBM is dominated by Bi 6p states, while the VBM is mainly comprised of O 2p and X
np (n = 3, 4, 5 for X = Cl, Br, I, respectively) states. Thus the reduced electronegativity (Cl > Br > I) can explain the decreasing band gaps. As the halogen atom X gets heavier, the contribution of X
np to the density of states becomes more and more obvious, and the interatomic covalent characteristics from strong to weak are in this order: Bi-O, Bi-I, Bi-Br, and Bi-Cl. The influence of I on Bi is more obvious than that of Cl and Br, so the deviation of Bi at the CBM in BiOI is more significant than those in BiOCl and BiOBr.
The thickness-dependent optical properties of BiOX (X = Cl, Br, I) is also investigated. At the same layer number, the calculated static dielectric constant
ɛ(0) and static refractive index
n(0) become larger with the increase of atomic number of halogen X (Table S3), as is normally expected.
Figure 4 takes the optical properties of BiOCl as an example to show that the thickness dependence is similar to those of multi-layer Bi
2O
2X (X = S, Se, Te) [
21] and MoS
2 [
50]. As the layer thickness increases from monolayer to bulk,
ɛ(0),
n(0), and the optical absorption coefficient become larger, and the peaks of the imaginary part of the dielectric function and the optical absorption spectrum shift to low-energy region (redshift).
Holding the above data showing that monolayer number is effective in tuning the bandgap, dielectric constant, light absorption spectrum of 2D BiOX, one may be more achievable to design BiOX-related optoelectronic devices in extended application ranges.
3.3. Strain engineering and composition tuning
Strain engineering is often considered an effective technique for tuning the electronic structure of 2D semiconductors [
21,
51]. Thus we here investigate the tunable properties of 2D BiOX (X = Cl, Br, I) with strain.
Figure 5 (a) depicts the band gap variation of BiOX monolayer with uniaxial and biaxial strain, and Figure S3 shows the band structures of BiOX monolayer under biaxial strain -8%, 0, and 8%. Under uniaxial strain, the band gaps of BiOCl and BiOBr monolayers decrease monotonically as the strain varies from -8% compression to 8% tension. In contrast, the band gap of BiOI decreases linearly with the increase of tensile strain, while it increases first and then decreases as the compressive strain becomes larger with a critical strain at -4%. BiOX monolayers are all indirect band gap semiconductors with or without uniaxial strain, and their CBM always remains at the Г point. When biaxial strain is applied, the band gaps of BiOCl, BiOBr, and BiOI show similar variation trend. With the increase of compressive strain, the band gaps increase linearly first and then decrease with a critical strain at -6% and -4% for BiOCl/BiOBr and BiOI, respectively. When the compressive strain increases to -8%, the CBM of BiOX shifts from the Г point to the M points, and in the range of compressive strain (-8%, -6%), the CBM and VBM of BiOI are located at the M point and the X point, respectively. With the increase of tensile strain, the band gaps of BiOX decrease monotonically and the CBM of BiOX remains at the Г point. All of monotonic changes versus strain sounds almost linear, suggesting high controllability of the 2D BiOX under strain.
In Figure S3, the CBM shows a sharper band curve, implying BiOX (X = Cl, Br, I) monolayer possesses smaller electron effective mass. Consistently, the calculated electron (hole) effective masses are 0.28 m0 (1.33 m0), 0.24 m0 (1.12 m0), and 0.19 m0 (1.96 m0) for BiOCl, BiOBr, and BiOI monolayer, respectively. Under tensile strain, the CBM of BiOX monolayer exhibits greater curvature than under compressive strain. It is thus expected that tensile strain will bring smaller electron effective mass and better carrier transport property than compressive strain. This is rather helpful to practical performance since tensile strain is more widely used in application of 2D materials. Under biaxial strain, both BiOBr and BiOI monolayers always exhibit indirect band gap regardless of the compressive or tensile strain applied. However, when the tensile strain increases to 8%, we have found that BiOCl monolayer has a transition from indirect band gap to direct band gap, and the CBM and VBM are both located at the Г point. For both uniaxial and biaxial strain, the band gaps of BiOX monolayers tuned under tensile strain exhibit better linearity, while a critical strain occurs under compressive strain. The critical strain of BiOI is smaller because the smaller electronegativity of I makes it more susceptible to strain. In the strain range of (-4%, 8%), uniaxial and biaxial strain have similar influence on BiOI. BiOCl and BiOBr exhibit better tunability under biaxial strain in the strain range from -6% to 8%. Generally, the bandgap of BiOX is found tunable in a wide strain range, suggesting application potential in the fields such as flexible optoelectronics.
Band edge position is an important factor in applications such as photocatalytic water splitting [
35,
52,
53]. As depicted in
Figure 5 (b), the band edge positions of BiOX (X = Cl, Br, I) monolayer (referred to the vacuum level) shift downwards under tensile strain, and the downward trend tends to be gentler when the tensile strain goes larger. With the increase of compressive strain, the VBM position goes upwards while the CBM position shifts first upwards and then downwards with a critical strain at -6%. It is thus a feasible method to adjust the band edge positions by adopting appropriate strain. In the range of -2% to -8%, the CBM and VBM positions of strained BiOX shift upwards to higher energy positions with respect to the redox potential levels of water [
27], meeting the band edges alignment requirements [
53] for photocatalytic water splitting at pH = 0, that is, the CBM is higher than the reduction potential of H
+/H
2 (-4.44 eV) and the VBM is lower than the oxidation potential of O
2/H
2O (-5.67 eV). Obviously, the strain engineering has the potential to improve the photocatalytic performance of BiOX at different pH conditions. For instance, water splitting can be switched on/off simply by applying/relaxing a compressive strain, and the efficiency of photocatalytic water splitting at pH = 0 might be controlled by applying compressive strain from -2% to -8%.
The tunability of 2D BiOX is further studied by adjusting the halogen composition
x in BiOCl
xBr
1-x, BiOCl
xI
1-x, and BiOBr
xI
1-x. Similar to the case x = 0.5 as in
Figure 1, the models for
x = 0.75, 0.875 are constructed by replacing one of the 4, 8 halogen atoms with another element in a 2 × 1 × 1, 2 × 2 × 1 supercell, respectively. In unstrained cases, the band gaps decrease gradually as the
x goes smaller, making these materials suitable for absorbing wider range of visible light compared with the original BiOX (X = Cl, Br, I). Combined with strain engineering, tensile strain monotonically reduces the band gap of composition-tuned BiOX monolayer, as shown in
Figure 6 (a) - (c). When the proportion of I in X is relatively large, the bandgap of BiOX is more affected by I, which may be because I is heavier and its bonding with Bi is different from Cl and Br. According to the projected density of states of BiOBr
0.5I
0.5 (
Figure 6 (d) - (f)), the unstrained CBM is dominated by Bi 6p orbitals, while its VBM is mainly contributed by I 5p states with small amount of O 2p states. We find that there is no Br 4p states below the Fermi level, which might be due to the lower energy of Br 4p orbital than that of I 5p. With the increase of compressive strain, the contributions of I 5p states to VBM and Bi 6p states to CBM are more prominent. On the contrary, as the tensile strain increases, the contribution of O 2p and Br 4p states to VBM becomes more and more obvious. Therefore, the band gap difference between BiOBr
0.5I
0.5 and BiOI becomes greater with the increase of the tensile strain. These results indicate that both composition tuning and strain engineering are effective ways to tune the band gap of BiOX (X = Cl, Br, I). It is predicted that composition tuning and/or tensile strain can decrease the band gap, enhance the optical absorption coefficients and the photocatalyst activity of BiOCl and BiOBr in the visible light region. This flexibility in bandgap tuning will benefit the material design of 2D BiOX and its applications in optoelectronic devices. As a whole, although those with I is weakly changeable when x ≤ 0.5, these BiOX exhibit effective controllability under composition tuning, which is beneficial to many applications.
Next we try studying the optical properties, as are important to optoelectronic materials [
18,
52], of BiOX (X = Cl, Br, I, Cl
0.5Br
0.5, Cl
0.5I
0.5, Br
0.5I
0.5). For unstrained BiOX, the static dielectric constant
ɛ(0) and the static refractive index
n(0) decrease in the order of BiOI, BiOBr
0.5I
0.5, BiOCl
0.5I
0.5, BiOBr, BiOCl
0.5Br
0.5, and BiOCl. This order follows the decreasing average atomic number of halogen X, so it naturally brings about the blueshifting absorption spectrum as shown in
Figure 7 (a). This shift suggests that the absorption spectrum can be tuned on demand in the range from visible to ultraviolet. Specifically, BiOI, BiOBr
0.5I
0.5 and BiOCl
0.5I
0.5 maintain absorption coefficient as high as 10
4 - 10
5 cm
-1 within the energy range of ~ 3.1 - 4.1 eV. These materials may thus be used to compensate the absorption drawback of silicon solar cells in this range to improve the utilization rate of solar energy. Hence the halogen composition tuning enhances the application potential of BiOX in the visible to violet region. With applied strains, the
ɛ(0) and
n(0) change slightly, suggesting a certain stability as the material is tuned by strain. As represented by BiOBr,
Figure 7 (b) shows the absorption spectra under strains. It is observed that the optical absorption edges and the absorption peaks redshift under tensile strain, which can improve the ability of absorbing visible light. As the compressive strain increases from 0% to -6%, the optical absorption edges blueshift. Together with the fact that the absorption coefficients of strained BiOX in the ultraviolet region is not degraded by compressive or tensile strain, the above properties imply that BiOX has application potential in flexible optoelectronic devices.