A First-Principles Study of Gas Molecule Adsorption on Carbon-, Nitrogen-, and Oxygen-Doped Two-Dimensional Borophene

A first-principles study was performed to investigate the adsorption properties of gas molecules (CO, CO2, NO, and NO2) on carbon(C-), nitrogen(N-), and oxygen-doped (O) borophene. /e adsorption energies, adsorption configurations, Mulliken charge population, surface work functions, and density of states (DOS) of the most stable doped borophene/gas-molecule configurations were calculated, and the interaction mechanisms between the gas molecules and the doped borophene were further analyzed. /e results indicated that most of the gas molecules exhibited strong chemisorption at the VB site (the center of valley bottom B–B bond) of the doped borophene (compared to pristine borophene). Electronic property analysis of the C-doped borophene/CO2 and the NO2 adsorption system revealed that there were numerous charge transfers from the C-doped borophene to the CO2 and NO2 molecules. /is indicated that C-doped borophene was an electron donor, and the CO2 and NO2 molecules served as electron acceptors. In contrast to variations in the adsorption energies, electronic properties, and surface work functions of the different gas, C-, N-, and O-doped borophene adsorption systems, we concluded that the C-, N-, and O-doped borophene materials will improve the sensitivity of CO, CO2, and NO2 molecule; this improvement of adsorption properties indicated that C-, N-, and O-doped borophene materials are excellent candidates for surface work functions transistor to detect gas molecules.


Introduction
Graphene is a new two-dimensional nanomaterial that was first prepared in 2004 by Novoselov et al. [1] using mechanical exfoliation. Since then, the structures, properties, and applications of new two-dimensional monolayer materials have become some of the most interesting areas of study. Borophene, a new type of two-dimensional singleatom layer nanomaterial, has been successfully synthesized on single crystal silver (crystallographic plane of (111)) substrates under ultrahigh-vacuum conditions [2]. Its unique geometrical and electronic structure endows borophene with excellent physical, chemical, and mechanical properties, such as high thermal conductivity [3,4], anisotropic electronic [5] and optical properties [6], superconducting properties [7], and ultrahigh mechanical modulus [8]. ese properties make borophene suitable for various applications, such as alkali metal ion batteries [9,10], hydrogen storage [11,12], supercapacitors [13], and gas molecule sensors [14][15][16][17]. In addition, harmful gases such as CO, NO, NO 2 , NH 3 , SO 2 , and HCOH cause environmental pollution, which is becoming increasingly serious.
us, studies on the adsorption and detection of harmful gases are important, and the development of new harmful gas sensors with high sensitivity and fast response is gaining significant attention. Studies have shown that 2-Pmmn two-dimensional borophene possesses an extraordinary wrinkled geometry [2], excellent surface-to-volume ratio, and anisotropic electronic properties, suggesting that the charge states and work functions of borophene could be easily modified by gas molecule adsorption, allowing it to be used in sensors and gas detection. Numerous studies have reported the first-principles calculations, for the adsorption of borophene for gas molecules. Liu et al. [14] calculated the adsorption configurations, adsorption energies, and electronic properties of gas molecules absorbed on borophene and found that CO, CO 2 , NH 3 , NO, and NO 2 were chemisorbed on borophene, while CH 4 physisorbed on borophene. eir results revealed that borophene could be used as a highly sensitive CO and CO 2 gas sensor. Cui [16] used the density functional theory (DFT) method to study the adsorption behavior of borophene for SO 2 molecules and showed that borophene had good adsorption capacity toward SO 2 molecules by chemisorption. e calculations also indicated meaningful applications of the materials as SO 2 adsorbents in the near future. Tan et al. [17] examined the possibility of conductive borophene nanosheets as promising sorbent materials for charge-modulated switchable CO 2 capture, and the results indicated that the binding strength of CO 2 molecules on negatively charged borophene could be significantly enhanced by injecting extra electrons into the adsorbent. Nagarajan and Chandiramouli [15] also calculated the electronic and ethanol adsorption properties of a hydrogenated 2D-borophene nanosheet device using nonequilibrium Green's function (NEGF) and DFT methods. eir results also confirmed that borophene nanosheets could be used as ethanol detectors. Furthermore, the adsorption energy, adsorption sites, and charge transfer of CO and NO on borophene were analyzed at the atomistic level [18], and their finding showed that borophene could be used for the development of chemical nanosensors to detect CO and NO molecules. Shen et al. [19] first designed a borophene/MoS 2 heterostructure and investigated the adsorption behavior of small gas molecules (CO, CO 2 , NO, NO 2 , and NH 3 ) on borophene surfaces, and their results showed that, except for CO 2 , the other four molecules adsorbed on borophene through chemisorption. However, all of the above studies focused on the calculations of pristine borophene on gas molecules. e edge functionalization, the nonedge defects, and chemical decorations of the twodimensional materials can provide additional modes of functionalization to tailor the electronic, chemical, optical, and magnetic properties of materials. e modified borophene materials and devices with excellent adsorption properties are still expected for next-generation gas molecules sensors, with the inadvertent introduction of various topological defects (vacancies, substitutional atoms, or adatoms) and other large defects (nanopores and nanoholes), which were unavoidable during the preparation and manipulation of borophene based devices. Previous studies have shown that structural defects in borophene sheets greatly affected their properties [20,21]. For example, using DFT calculations, Li et al. [12], Chen et al. [22], and Wang et al. [23] found that Li-, Na-, Ca-, and K-decorated borophene could be a promising H 2 storage material for H 2 fuel battery. In addition, Kumar [24] explored the structural parameters and gas adsorption properties of CO, NO, CO 2 , NO 2 , and NH 3 gas molecules adsorbed over transition metal atom embedded borophene sheets and found that gas adsorption could be enhanced by embedding transition metal atoms into borophene sheets. Furthermore, the effects of sodium-decorated borophene on the adsorption of CO and CO 2 molecules were determined [25], and they found that Na-decorated borophene monolayers could be used as sensors to detect CO and CO 2 molecules in the atmosphere. erefore, it is important to investigate the effects of dopants on the gas molecule adsorption properties of borophene sheets.
In this work, we used first-principles calculations based on DFT to study the adsorption behavior of CO, CO 2 , NO, and NO 2 small gas molecules on C-, N-, and O-doped borophene sheets. By analyzing the adsorption energy, adsorption configuration, Mulliken charge population, work functions, and density of states (DOS) of gas molecules on doped borophene, we concluded that the C-, N-, and O-doped borophene materials may be excellent candidates for gas detection applications.

Methods and Calculations
All calculations were performed using the Cambridge Sequential Total Energy Package (CASTEP) based on the DFT [26]. In these calculations, the Perdew-Burke-Ernzerhof (PBE) function of generalized gradient approximation (GGA) was used to assess the exchange correlation energy between the electrons [27], and the ultrasoft pseudopotential was used to describe the interactions between the ion core and electrons [28]. Structural optimization and energy calculations were carried out using a plane wave basis set to a cut-off energy of 400 eV and a 4 × 4 × 1 k-point mesh in the Brillouin zone. e configurations for the gas molecule-borophene systems were optimized by fully relaxing the atomic structures until the remaining forces were smaller than 0.03 eV/Å. Energy convergence was 10 −5 eV/atom between the two steps, and a vacuum of 20Å was employed along the Z-direction of the doped borophene sheets to eliminate interactions between the borophene sheets.  Figure 1). e C-, N-, and O-doped atoms were substituted at the T-site or V-site with the boron atom in borophene. e small gas molecules were then positioned vertically in the four highly symmetrical adsorption sites of doped borophene, and the best adsorption sites and adsorption energies were obtained after structural optimization. e adsorption energy of the small gas molecules on C-, N-, and O-doped borophene was defined by where E molecule/doped−borophene , E doped−borophene , and E molecule are the total energy of the optimized adsorption configuration of the small molecule on the doped borophene system, the total energy of the doped borophene sheet, and the total energy of the isolated small gas molecule, respectively. In the calculations of adsorption energies, the zero-point energy (ZPE) contributions were neglected, owing to Gholizadeh and Yu [29] who found that the zero-point energy (ZPE) contributions can be omitted in most cases. To calculate the DOS, the k-point was set to 15 × 15 × 1 for accuracy.  Table 1. From Figure 2(a) and Table 1, for the pristine borophene sheet, the bond length of B20-B21 and B21-B22 bond is 1.880Å and 1.626Å, respectively (consistent with the results of 1.879Å and 1.613Å from Peng [6]). e bond population of both B20-B21 and B21-B22 bond is 0.20, and the charge populations of boron atom are 0.0 e. After the carbon atom, nitrogen atom, and oxygen atom substituting the boron atom, the perfect geometric structures of borophene were disturbed (the top view and side view of Figures 2

Geometries and Electronic
e C-doped borophene lead to the B20-C bond (length of 1.541Å) and B23-C bond (length of 1.538Å) formation, and our calculated C-B bond length of borophene sheet is in accord with Gholizadeh's reports [30] of B-C bond length of 1.493Å for boron atom-doped graphene sheet. e bond populations of both B20-C bond and B23-C bond are 1.16, which are much larger than the bond populations (0.20) of B-B bonds for pristine borophene. From the point of view of the atomic charge populations, it is can be seen that the carbon atom gains 0.3 e, of which 0.26 e was contributed by the two boron atoms (both B20 and B23 lost 0.13 e). e nitrogen atom doping the borophene results in the geometric structures of perfect borophene sheet extensive distortion (Figure 2(c)).
e nitrogen atom and their nearest three boron atoms (B20, B21, and B23) form a triangular pyramid structure. And the bond populations of B20-N (bond length of 1.499Å), B21-N (1.500Å), and B23-N (1.500Å) are 0.91, 0.58, and 0.76, respectively. e nitrogen atom obtains 0.6 e from the nearest three boron atoms (B20, B21, and B23). When the oxygen atom is doping borophene, the B20, O, and B23 atoms form a B-O-B triangular structure, and a six-membered ring vacancy defect is left in the borophene sheet. e bond populations of B20-O (1.412Å) and B23-O (1.507Å) are 0.71 and 0.53, respectively. e oxygen atom obtains 0.62 e from the nearest two boron atoms (B20 and B23). Based on the analysis of the geometric structures and charge populations of C-doped, N-doped, and O-doped borophene, we can conclude that the doped heteroatoms (carbon nitrogen and oxygen) interact strongly with the boron atoms of borophene by a large amount of charge transfer from boron atoms to doped heteroatom.

Analysis of Doped Borophene Adsorption.
e adsorption energy variations of the CO, CO 2 , NO, and NO 2 molecules on C-, N-, and O-doped borophene at different adsorption sites are shown in Figure 3, and the optimized adsorption structures are shown in Figure 4.
For most adsorption configurations, with variations in adsorption site from the T-site, TB-site, V-site, to the VBsite, the adsorption energies of CO, CO 2 , NO, and NO 2 on the C-, N-, and O-doped borophene decreased. Except for the adsorption configuration of CO on N-doped borophene (the most stable adsorption site was the T-site), the most stable adsorption site of the other adsorption configurations was the VB-site (the center of the valley bottom B-B bond). In addition, the adsorption energies of the most stable adsorption configurations are shown in Figure 3(d), where the adsorption energies of CO on C-, N-, and O-doped borophene were −1.45, −1.70, and −1.86 eV, respectively, and these values were larger than those for the adsorption of CO on pristine borophene (−0.761 eV) [14] and the adsorption of CO on graphene (−0.12 eV). Besides, we also systematically analyzed the variations of bond length and bond angle before and after gas small molecule adsorption, as shown in Table 2. e bond length of C-O is 1.220Å, 1.191Å, and 1.206Å after CO molecule adsorbed on the C-, N-, and O-doped borophene, respectively. Compared with the free CO molecule (bond length of 1.154Å), its bond length increases slightly after adsorption. ese results showed that CO was chemisorbed on the C-, N-, and O-doped borophene in a perpendicular direction to the borophene sheet, and the adsorption structure is shown in Figure 4. In addition, similar adsorption structures were observed, with six-membered vacancy rings of the doped borophene and CO molecule adsorbing vertically on the borophene (Figure 4 (a1)-(a6)). e adsorption energies of CO 2 on C-, N-, and O-doped borophene were −0.64, −1.34, and −0.71 eV, respectively. Compared with the adsorption energy of CO 2 on pristine borophene (−0.898 eV), the carbon atom and oxygen atom doping resulted in a decrease in the adsorption energy of borophene for CO 2 , while nitrogen atom doping led to an increase in adsorption energy. As shown in Figure 4   . While the adsorption mechanism of the NO 2 molecule on doped borophene differed from that of the CO, CO 2 , and NO 2 molecules, as shown in Figure 4 (d1)-(d6), the NO 2 molecule dissociated into N and NO, which then adsorbed on the B atoms and formed O2-N2 bond (bond length of 1.247Å, 1.262Å, and 1.264Å for C-, N-, and O-doped borophene, respectively) and O3-B1 bond (bond length of 1.221Å, 1.219Å, and 1.223Å for C-, N-, and O-doped borophene, respectively). is indicated that NO 2 adsorption on C-, N-, and O-doped borophene was more stable. e calculated adsorption energies and configurations of the small gas molecules and doped borophene revealed that CO 2 chemisorbed on C-and O-doped borophene with a more moderate adsorption energy value compared to pristine borophene. us, the C and O atom-doped borophene materials would be more suitable as CO 2 molecule sensors. Although the C, N, and O atom-doped borophene could partly reduce the adsorption energy of the borophene sheet for NO, the adsorption energy was still approximately −2.0 eV. Additionally, the adsorption energies of CO and NO 2 molecules on C-, N-, and O-doped borophene increased. erefore, CO, NO, and NO 2 chemisorbed on C-, N-, and O-doped borophene with stronger adsorption energies, and desorption did not easily occur.

Work Functions for Different Gas-Doped Borophene
Adsorption Configurations. According to the different detecting principles, gas molecule sensors can be divided into chemical resistors, field effect transistors, and surface   [31]. e chemical resistors are one of the most popular gas sensors, whose working mechanism is to detect the changes in electrical conductance (band gap value) caused by adsorption of small gas molecules. However, we found that the band gap of heteroatom-doped (C-, N-, and O-) borophene is 0 eV before and after gas small molecules adsorption based on the DFT calculations. is phenomenon limits the application of heteroatom-doped borophene as a chemical resistor to detecting gas small molecules. Other widely studied sensing devices are surface work function changed transistors which depend on the work function change caused by the gas small molecules adsorption. Work function is defined as the minimum energy required to move an electron from the Fermi level to the vacuum level. e variation of surface work function is a measurable quantity after gas small molecules adsorption on the adsorbent. erefore, through analysis, the change of surface work function caused by the gas molecules adsorption, the heteroatom-doped borophene can be used as a surface work function transistor to detecting gas small molecules. Figure 5 shows the changes of surface work function caused by the heteroatom-doped borophene and gas small molecules adsorption. e work function of pristine borophene is 5.157 eV, which is consistent with the 5.310 eV from Kistanov's calculations [32]. After the carbon, nitrogen, and oxygen atom doping the borophene, the work function of C-doped, N-doped, and O-doped borophene is 4.928 eV, 4.735 eV, and 5.26 eV, respectively. Due to the charge transfer between the borophene and heteroatom (C, N, and O), the carbon and nitrogen atom doping leads to the work function decreases, but the oxygen doping results in the work function increases. From Figure 5(b) and the red column inset, the CO 2 , NO, and NO 2 molecules adsorbed on the C-doped borophene result in the work function increases of 0.043 eV, 0.083 eV, and 0.473 eV, respectively, but the CO molecule adsorbed on the C-doped borophene leads to the work function decrease of 0.235 eV. e variation of work function of heteroatomdoped borophene larger than 0.40 eV was caused by gas small molecules indicating that the heteroatom-doped borophenes are very sensitive to these gas small molecules [31]. e CO, CO 2 , NO, and NO 2 molecules adsorbed on the N-doped borophene lead to the work function increases, and the increment is 0.548 eV, 0.662 eV, 0.165 eV, and 0.723 eV, respectively ( Figure 5   and the N-doped borophene can be used as surface work function gas sensor to detect the CO, CO 2 , and NO 2 molecule, and O-doped borophene based gas sensor can detect the CO 2 and NO 2 molecule.

Electronic Properties of C-Doped Borophene/CO 2 and NO 2
Systems. Based on the Mulliken [33] populations algorithm, the population distribution of charges for the C-doped borophene/CO 2 and C-doped borophene/NO 2 systems before and after adsorption was calculated, and the results are presented in Table 3. e population changes for the most stable C-doped borophene/CO 2 adsorption system are presented in the first section of Table 3, which shows that the C2, O2, and O3 atoms in CO 2 gained 0.71, 0.02, and 0.01 e, respectively. In addition, the B1, B2, and B3 atoms in C-doped borophene in contact with the adsorbed CO 2 molecule lost 0.09, 0.55, and 0.43 e, respectively (C1 atom of C-doped borophene gained 0.1 e). For the most stable adsorption configuration of the C-doped borophene/NO 2 system, the N2, O2, and O3 atoms in the NO 2 molecule gained 0.51, 0.01, and 0.47 e, respectively. e B1, B2, and B3 atoms lost 0.61, 0.40, and 0.35 e, respectively. erefore, we found significant differences in gained electrons for the O2 and O3 atoms, mainly because the O3 atom first separated from the NO 2 molecule and then bonded with the B1 atom in doped borophene (Figure 4 (d1) and (d2)). us, charge transfer occurred mainly in the 2s and 2p orbitals of the atom. ese results showed that both the CO 2 and NO 2 gas molecules acted as electron acceptors, and the C-doped borophene acted as an electron donor.  To further understand the interaction mechanisms between borophene and the gas molecules after adsorption, the changes in the DOS of the C-doped borophene and gas molecules before and after adsorption were analyzed. We calculated the DOS (including the total density of states (TDOS) and partial DOS (PDOS)) of the C-doped borophene, C-doped borophene/CO 2 adsorption, and C-doped borophene/NO 2 adsorption systems, as shown in Figure 6.
Comparing Figures 6(a), 6(b), and 6(c) and the TDOS and PDOS before and after adsorption of the small gas molecule on C-doped borophene, we concluded that gas molecule adsorption on C-doped borophene caused charge transfer, which then broke and formed covalent bonds. When the CO 2 molecule adsorbed on C-doped borophene, the DOS became deeper (from −16.5∼5 eV to −25∼5 eV).
is DOS range, from −25 to −20 eV, was mainly due to the 2s and 2p orbitals of the C2 atom and the 2s orbitals of the O2, O3 atoms ( Figure 6(b)). In addition, the 2s and 2p orbitals of the B, C1, C2, O2, and O3 atoms overlapped, and the obvious resonance C-doped borophene/CO 2 adsorption structure indicated covalent bonds between the C2-B1, O2-B2, and B3-O3 atoms (Figure 4 (b1) and (b2)). e analysis results of the C-doped borophene/NO 2 configuration of DOS (Figure 6(c)) were similar to the results of the C-doped borophene/CO 2 system. After the NO 2 molecule adsorbed on the C-doped borophene sheet, the DOS also became deeper, and the deeper DOS level was mainly due to the 2s and 2p orbitals of the N2 and O2 atoms. As shown in Figure 6(c), we also found obvious orbital hybridization of the 2s orbitals in the B, N2, and O2 atoms, from approximately −25 to −22.5 eV, as shown by the interactions between the B, N2, and O2 atoms. However, there was no obvious orbital hybridization between the N2 and O3 atoms from −25 to 5 eV, as the O3 atom had already separated from the NO 2 molecule. Between −20 and −17.5 eV, both the 2s and 2p orbitals of the B and O3 atoms overlapped, indicating that B1-O3 formed covalent bonds (Figure 4 (d1) and (d2)). erefore, the DOS analysis results were consistent with the atomic configuration analysis and Mulliken charge population results.

Conclusions
Using first-principles calculations based on DFT, we studied the adsorption behavior of CO, CO 2 , NO, and NO 2 gas molecules on C-, N-, and O-doped borophene sheets. By comparing the adsorption energies and adsorption configurations of the gas molecules on the four adsorption sites in borophene, we determined that most of the gas molecules formed relatively stable chemisorption bonds at the VB site (the center of the valley bottom B-B bond) on C-, N-, and O-doped borophene (only the CO molecule adsorbed on the T-site in N-doped borophene). Compared to pristine borophene, C, N, and O atom doping increased the adsorption energies of the gas small molecules on borophene, making the adsorption configuration more stable. Furthermore, by analyzing the variations in DOS before and after adsorption, the effects of gas molecules on the electronic states of the C-, N-, and O-doped borophene were investigated to provide a microscopic understanding of the interactions between the gas molecules and the borophene sheet. In contrast to variations in the adsorption energies, electronic properties, and surface work functions of the different gas, C-, N-, and O-doped borophene adsorption systems, we concluded that the C-, N-, and O-doped borophene materials will improve the sensitivity of CO, CO 2 , and NO 2 molecule. is improvement of adsorption properties indicated that C-, N-, and O-doped borophene materials are excellent candidates for surface work functions transistor to detect gas molecules.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.