W. P.
Xu
and
H.
Xu
*
Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China. E-mail: xu.h@sustc.edu.cn
First published on 30th November 2017
H2O and O2 are believed to be key factors that influence the structural stability of black phosphorus (BP) in ambient conditions. In this work, the interactions of H2O and/or O2 with BP are investigated using first-principles calculations. The results indicate that water molecules prefer to adsorb on the BP surface and form a six-member water ring. The dissociation barrier of O2 is significantly reduced in the presence of H2O, which dramatically promotes the degradation of BP. Moreover, the introduction of O2 also facilitates the adsorption of water on the surface. The effects of H2O and/or O2 on the quasiparticle band gap and exciton binding energy of BP are also investigated. The results suggest that water adsorption has only a slight influence on the electronic properties and exciton binding energy, while O2 adsorption causes obvious changes in the properties of BP, which results in a direct-to-indirect band gap transition in BP.
The interactions of O2 and/or water with BP have recently received considerable academic attention. Theoretically, Ziletti et al.20 found that O2 interacts strongly with the BP surface in the dissociative form with an energy release of 2.08 eV per O atom. The reaction barrier in the process of O2 dissociation is only 0.15 eV, and the dissociated O2 serves as anchors to improve the adsorption ability of water by forming hydrogen bonds.20 In addition, the calculations performed by Wang et al.21 indicated that O2 can spontaneously dissociate on a BP surface at room temperature, while water only weakly binds to pristine BP. Experimentally, the synergetic effect of O2 and H2O in the oxidation of BP has been confirmed by X-ray photoelectron spectroscopy (XPS) measurements.22 Moreover, the site-selective oxidation of BP after exposure to high-purity O2 or H2O has been systematically studied using XPS and transmission electron microscopy (TEM) observations,23 and an oxidation mechanism involving a multistep reaction has been proposed. Despite intensive efforts, the interplay between water molecules on the BP surface and the synergetic effect of O2 and H2O during BP oxidation at the atomic level remain unclear. More importantly, it is desirable to investigate the influence of H2O and/or O2 adsorption on the electronic properties of BP, especially on the quasiparticle (QP) band gap and EBEs.
Fig. 1 The adsorption of (a) a single H2O molecule, (b) three H2O molecules and (c) six H2O molecules on the p(6 × 6) BP surface. |
When the number of H2O molecules further increases, the H2O molecules still tend to form a ring structure with more intermolecular hydrogen bonds. In the H2O ring, most non-hydrogen-bonded H atoms point toward the BP surface. However, several H atoms point upwards, which could reduce the repulsive forces between H atoms and further reduce the total energy. The six-member ring is an example: six hydrogen atoms in the H2O molecules form a ring through hydrogen bonds, four hydrogen atoms in the H2O molecules point toward the BP surface, and the remaining two H atoms of diagonal H2O molecules point upwards (see Fig. 1(c)). The average adsorption energy per H2O of a six-member H2O ring is −0.498 eV. As shown in Table 1, our results indicate that the adsorption energy reaches its minimum (−0.498 eV) when the six-member ring forms, which suggests that the optimal number for H2O adsorption on a BP surface is six. Furthermore, we also investigated the stable configuration at temperatures of 60 and 100 K using molecular dynamics (MD) calculations. The results also show that the average adsorption energy is the lowest when H2O molecules adopt a ring-like structure, which is consistent with the above findings.
nH2O | 1 | 2 | 3 | 455 | 5 | 6 | 7 |
---|---|---|---|---|---|---|---|
E 1 | −0.080 | −0.196 | −0.283 | −0.356 | −0.389 | −0.391 | −0.388 |
E 2(vdw) | −0.185 | −0.306 | −0.392 | −0.470 | −0.490 | −0.498 | −0.496 |
ΔE | 0.095 | 0.110 | 0.111 | 0.114 | 0.101 | 0.107 | 0.108 |
In ambient conditions, H2O molecules always coexist with O2 molecules on the BP surface. In the following, we investigated the possible adsorption structure of the water molecules on the BP surface in the presence of O2. The dissociation of O2 on the BP surface has been intensively studied in the literature.20 To verify the dissociation process, we also calculated the O2 dissociation on the BP surface. The calculated O–P bond length is 1.505 Å, and the adsorption energy per O atom is −2.10 eV. The calculated adsorption energy per O atom for O2 dissociation at the ortho-site (see Fig. 2(a)) is 80 meV lower than at the meta-site (see Fig. 2(b)), which is consistent with prior results.20 The two O atoms at the ortho-site with large separation experience a weak repulsive force, which results in a lower adsorption energy.
Fig. 2 A single O2 molecule that is dissociated at (a) the ortho-site and (b) the meta-site of the p(6 × 6) BP surface. |
During O2 dissociation on the BP surface, the dissociated O2 initially bonds to two meta-site P atoms.20 Although it is not the energetically favorable configuration, it can form on the BP surface because of the larger diffusion barrier of the oxygen atom (∼2 eV) to form the ortho-site configuration. Moreover, the calculated average adsorption energy per H2O is lower for the meta-site configuration (see Fig. 3). Thus, we only considered the meta-site configuration below. The results show that each H atom in H2O points to its neighboring O atom, forming two hydrogen bonds, as shown in Fig. 4(a). Compared to the H2O adsorption on a pristine BP surface, the introduction of O2 is predicted to increase the hydrogen bonding and further reduces the average adsorption energy of the H2O molecules. The calculated adsorption energy per H2O with O2 dissociation (−0.389 eV) is 0.204 eV lower than on the pristine BP surface, which indicates the hydrophilic properties of the oxidized BP. Similarly, as shown in Fig. 4(b), three H2O molecules form a ring structure that is bonded by three pairs of hydrogen bonds. Two of the unpaired H atoms between the neighboring H2O molecules point to the adsorbed O atoms and the rest of the H atoms point upwards. The ring structure again forms when more H2O molecules are adsorbed on the BP surface, such as the five-member ring that is shown in Fig. 4(c). The calculated average adsorption energy per H2O with different numbers of H2O molecules that is shown in Fig. 3 indicates that the five-member water ring has the lowest value of −0.532 eV. Therefore, the optimal number for H2O adsorption on a BP surface with dissociated O2 is five.
Fig. 3 The average adsorption energy of different water molecules with the dissociation of an O2 molecule at the meta-site and ortho-site on the p(6 × 6) BP surface. |
Fig. 4 The adsorption of (a) a single H2O molecule, (b) three H2O molecules, and (c) five H2O molecules on the p(6 × 6) BP surface with the dissociation of O2. |
Next, we performed calculations to reveal the influence of H2O on the dissociation barrier of O2. In the absence of water, the dissociation barrier of O2 on a p(3 × 5) BP surface is approximately 0.2 eV, as shown in Fig. 5(a), which is consistent with prior results.20 In the presence of one H2O molecule, both H atoms of the H2O molecule point toward the dissociated O2 to form two hydrogen bonds, which promote the decomposition of O2. The dissociation barrier of O2 is found to be almost 0 when H2O is introduced as shown in Fig. 5(b), indicating a spontaneous reaction. In other words, the presence of H2O molecules can dramatically accelerate the decomposition of O2 molecules on the BP surface. Simultaneously, the presence of O2 gives rise to an enhancement of water adsorption on a BP surface. Thus, O2 and H2O molecules have a strong synergetic effect in the oxidation of BP. Therefore, the degradation of BP not only starts from edges, but also occurs on the BP surface in wet conditions.
Fig. 5 (a) The O2 dissociation process on the p(3 × 5) BP surface. (b) The O2 dissociation process on the p(3 × 5) BP surface in the presence of one adsorbed H2O molecule. |
Because BP is easily damaged by H2O and O2 molecules, it is important to know whether the electronic properties suffer from oxidation. Because of the expensive computational cost, a p(2 × 2) supercell was used, and the models that were used are shown in Fig. 6. The calculated QP band gap and GW-BSE optical band gap of the pristine BP are 1.958 eV (see Fig. 7(a)) and 1.252 eV (see Fig. 7(e)), respectively. Therefore, the calculated EBE is 0.706 eV, which is consistent with the literature.34–41 As illustrated in Fig. 7(a), BP shows direct band gap semiconductor features, where both the conduction-band minimum (CBM) and valence-band maximum (VBM) are located at the Γ point. The projected density of states (PDOS) indicates that both VBM and CBM are mainly contributed by pz orbitals, which is in good agreement with prior results.42
Because the influence of O2 and/or H2O adsorption on the QP band structure and optical spectra within the GW-BSE have not yet been reported, it is particularly important to perform such calculations. When a single H2O molecule is adsorbed on a p(2 × 2) supercell, the band structure is almost unchanged, and CBM and VBM are still located at the Γ point. The orbitals from the adsorbed H2O are located at deep levels. The calculated QP band gap and GW-BSE optical band gap are 1.972 eV and 1.282 eV, as reported in Fig. 7(b) and (f), respectively. The calculated EBE is 0.690 eV, which is almost the same as the EBE of the pristine BP. Different from water, the adsorption of O2 on a BP surface has a significant influence on the electronic properties as shown in Fig. 7(c). The VBM is still at the Γ point, and it is also originated from the pz orbitals of P atoms. Oxygen adsorption has a significant impact on the CBM. The CBM deviates from the Γ point, which locates along Γ–M direction. The CBM is dominated by the px and py orbitals of P atoms bonding to oxygen, while 2p orbitals of O atoms also contribute to CBM. The formation of O–P bond is accompanied by a strong orbital hybridization, and it causes a decrement of the underlying P–P bond lengths, leading to an enhancement of the splitting between the pz bonding and antibonding orbitals of P atoms. As the pz antibonding orbitals of P atoms are raised in energy and the contributions from the px and py orbitals of P atoms increase, it results in a direct-to-indirect band gap transition. The QP and GW-BSE band gaps at the Γ point are 2.482 eV and 1.956 eV, respectively. The calculated EBE becomes 0.526 eV upon the adsorption of O2. In addition, we find that the co-adsorption of H2O and O2 on the BP surface is similar to O2 dissociation, which also leads to a direct–indirect band gap transition. As shown in Fig. 7(d) and (h), the QP and GW-BSE band gaps at the Γ point are 2.474 eV and 1.926 eV, respectively, resulting in an EBE of 0.548 eV. The above analysis explains why H2O adsorption has scarcely any influence on BP and why O2 adsorption and their co-adsorption change the electronic properties of BP.
This journal is © the Owner Societies 2018 |