Gas sensors have received considerable attentions because of their important applications in detecting toxic and deleterious gases, supervising air quality, and monitoring human health [
1,
2]. Compared to other gas detection techniques, electrochemical sensing has various advantages, such as, low energy linear output with high resolution, good selectivity and repeatability, high detection accuracy, and being more inexpensive than other techniques [
2]. As motivated by the driving force of enhancing the performance of gas sensors, great efforts have been exhausted on their design with tiny architectures or configurations [
1]. At present, numerous distinct gas sensors have been prepared, for instance, Zhang
et al. [
3] have built a gas sensor based on anatase TiO
2 that is used to detect decomposed gases in SF
6, Cui
et al. [
4] have constructed a field effect transistor sensor device which exhibits an ultrahigh sensitivity to NO
2 in dry air, Cho
et al. [
5] have demonstrated atomically thin heterostructured gas sensors which comprising a combination of exfoliated MoS
2 and CVD graphene, Zhao
et al. [
6] have made a MoS
2-based flexible and wearable gas sensor which was fabricated based on the direct synthesis of MoS
2 on a polyimide substrate.
A clear trend is, with the emergence of two-dimensional (2D) materials [
7–
9], gas sensors based on 2D materials have drawn more and more attentions because of their superior properties, including large surface to volume ratio, associated charge transfer between gas and surface, and high electron mobility [
10–
13]. Recently, graphene and other 2D materials-based gas sensors have been widely investigated both experimentally and theoretically [
14–
17]. Among these 2D materials, borophene, which has various excellent properties, such as anisotropic mechanical [
18–
20], optical [
21,
22] and electronic properties [
23–
25], plays a crucial role in promoting the development of 2D gas sensors [
8,
26,
27]. Shukla
et al. [
10] reported the realization of 2D monolayer borophene-based gas sensor using first principle density functional theory methods, and concluded that 2D monolayer borophene can be a potentially important substrate for gas sensing of CO, NO, NO
2 and NH
3 gas molecules. Furthermore, Shen
et al. [
28] studied the transport properties of gas adsorbed on 2D monolayer borophene with a MoS
2 substrate, i.e., on borophene–MoS
2 heterostructure, and results show that the heterostructure is also very sensitive to CO, NO, NO
2 and NH
3 gas molecules, especially to NO moleclue, and therefore can be used as a 2D gas sensor. However, these studies ignored the effect of metal electrodes, which is necessary for not only the collecting of electrical signal during the operation of the sensor, but also adjusting the Fermi energy and/or electrostatic potential of the 2D materials [
29–
31]. In addition, a comprehensive understanding on the effects of MoS
2 substrate, gas molecules, as well as metal electrodes, to the transport properties of 2D monolayer borophene is still lack.
In this paper, we systematically investigated the transport properties of borophene-based nano gas sensors with or without considering the influence of gas molecules, MoS
2 substrate and/or metal electrodes. Results show that the borophene and borophene–MoS
2 heterostructure with gold electrodes is sensitive to various gas molecules, therefore, they can be used as a gas sensor, which qualitatively agrees with previous studies [
10,
28]. The calculated current–voltage (
I–
V) characteristic through pristine borophene without considering the effect of MoS
2 substrate and gold electrodes under 0.5 V has a linear behavior, however, with considering the MoS
2 substrate and/or the metal electrodes, the current is reduced and a non-linear
I–
V behavior can be observed. Further band structure, charge and potential analyses well explained the obtained results. Our studies not only provided comprehensive information regarding the underlying physics of borophene-based gas sensors, but also can be used to understand more transport behaviors of 2D monolayer materials and heterostructures with metal electrodes, which are usually not included when theoretically studying the transport properties of 2D metallic materials in the literature.
We started by considering the fabrication and operation process of a real borophene-based nano gas sensor with gold electrodes, as shown in Fig. 1(a). Monolayer borophene is placed on MoS
2 substrate, and gold electrodes are deposited onto borophene. The distance between two gold electrodes is big enough to contain a gas molecule, so that the molecule can be regarded as no interaction with the electrodes. When a gas molecule is approaching and adsorbing on the borophene, the current is collected through the electrodes. Figures 1(b, c) show the side view and top view of the optimized atomic structure of the entire gas sensor. The distance between each two parts (borophene, MoS
2 or two gold electrodes) of the structure is optimized by vasp software [
32], the optimization details can be found in the supplementary material. To calculate the transport properties without MoS
2 substrate and/or gold electrodes, we just need to simply remove the corresponding part from the entire structure in Fig. 1. According to the previous studies [
10,
28,
33], in this work we considered four gas molecules, namely, CO, NO, NO
2 and NH
3, which are all sensitive to borophene-based gas sensors due to chemisorption between the molecules and borophene. These molecule–borophene structures are fully relaxed (see the supplementary material), while pristine borophene is used for all the systems without molecules. After the structure was determined, we carried the quantum transport calculations by density functional theory combined with the nonequilibrium Green’s function (NEGF-DFT) as implemented in nanodcal software package [
34,
35]. In our NEGF-DFT numerical calculations, double-zeta polarized atomic orbital basis set was used to expand the physical quantities of B and molecule elements, single-zeta polarized atomic orbital basis set was used for Mo, S and Au [
36]. The exchange and correlation were treated at the level of the Perdew–Burke–Ernzerhof (PBE) level [
37], and atomic cores are defined by the standard norm conserving nonlocal pseudopotentials [
38]. In the self-consistent calculations, 3 k-points are used to sample the Brillouin zone in the periodic direction, and a cutoff energy of 80 Hartree were used. To obtain the transmission and
I–
V curve, 50 k-points are used. Spin orbital coupling is not considered in this work.
Figure 2 shows the I–V curves of borophene-based nano gas sensors under 0.5 V from our NEGF-DFT numerical calculations on eleven different systems with or without considering the influence of gas molecules, MoS2 substrate and/or gold electrodes, from which the following results can be clearly observed:
(i) The
I–
V characteristics of pristine borophene (black curve) show a nearly perfect linear behavior, which maintains after a molecule (CO) is adsorbed (red curve), and the scattering of the CO molecule leads to the decreasing of the current. Our calculated
I–
V values quantitatively agree with the previous reports [
10]. We should note that for the pristine borophene system (black) and borophene–MoS
2 heterostructure system (orange) that will be discussed in the following, the atomic structures are actually perfect periodic along the transport direction. For a periodic system, it is experimentally unfeasible to direct apply finite bias on a selected region, but numerical calculations of these periodic systems are very useful for one to better understand the physics of transport properties through the structure.
(ii) After MoS
2 is introduced, the current values of heterostructure systems (orange and gray curves) decrease compared to the current values of the corresponding borophene systems. More interestingly, under 0.15 V the
I–
V characteristic show non-linear behavior, and then become nearly linear after 0.2 V, which is more obvious for the orange curve of borophene–MoS
2 system in Fig. 2. This non-linear behavior was not mentioned in previous studies, and should be further understood. Quantitatively, our calculated current values under the same voltage is smaller than that in the previous reports on borophene–MoS
2 heterostructure [
28], which is because we applied different strain on MoS
2 or borophene layer and used different stacking way to make the bilayer structure.
(iii) The rest seven
I–
V curves are all for the systems with gold electrodes. We can see from the supplementary material that all these curves lose the perfect linear behavior, which must stem from the existence of metal electrodes. In addition, after gold electrodes are included, the calculated current values under the same voltage is much smaller than the current values of the corresponding system without gold electrodes, although gold is a good conductor from a macro point of view. Similar experimental observations are reported by Ref. [
39] in the Al: Graphene system, where Al electrodes reduced the current of pure graphene.
(iv) It is obvious that for all the sensors, the current values of systems with molecules are smaller than that without molecules, which agrees with all the previous reports [
10,
28]. Moreover, one can see that the four purple curves quantitatively have nearly the same
I–
V characteristic, which qualitatively agrees with the previous results on monolayer borophene gas sensor for CO, NO, NO
2 and NH
3 gas molecules [
10]. Therefore, to further judge if the borophene-based nano-gas sensors can well distinguish these molecules, more studies on the sensitivities of the sensors to various molecules are necessary.
From the I–V curves in Fig. 2, we can preliminarily conclude that borophene-based nano heterostructure gas sensor with gold electrodes are sensitive enough to detect some small molecules. Till now, two behaviors observed in the I–V curves (Fig. 2) still need to be further understood. First, why MoS2 substrate lead to a non-linear I–V characteristic under small bias? Second, why gold electrode, which is a good conductor, lead to the decrease of current values compared to that in the corresponding systems without gold electrodes? In the following we will solve these two problems (See Fig. 5 and related discussions).
Figure 3(a) shows the band structure and density of states (DOS) of pristine borophene calculated from 8-atom unit cell that matches MoS2 lattice (see supplementary materials), one can observe that each energy band in the band structure cross the entire energy range from −0.5 to 0.5 eV in a nearly linear fashion, and the corresponding DOS nearly keep a constant value (∼1.5). This means, if an external bias under 0.5 V is applied onto the borophene, each incoming electron state in the left electrode at a certain moment k and energy E can find a way out through and outgoing electron state in the right electrode at the same moment k with a different energy E′. In this way, the current will be proportional to the width of the bias widow, namely, the applied external bias voltage, leading to a linear behavior on the I–V characteristic. After MoS2, whose band structure and DOS is shown in Fig. 3(b), and borophene forming a heterostructure, charge transfer will take place between these two materials. Figure 3(c) shows the band structure and DOS of borophene–MoS2 heterostructure, and the contribution of each element is decomposed in the band structure. Clearly, in G–X and G–Y region above the Fermi energy, Mo element contributes to the energy band, due to the charge transfer from Mo to B (See later for charge discussion). Moreover, this energy band is in a non-linear fashion and ends at 0.143 eV, leading to an abrupt increase of the DOS (∼4) above the Fermi energy. In this way, this energy band in both electrodes will contribute to the transport in a resonance way when the external bias is below 0.143 V. After the applied bias is over 0.143 V, this energy band in two electrodes will lose resonance, namely, the incoming electron state corresponding to the DOS peak in the left electrode miss the outgoing electron state corresponding to the DOS peak in the right electrode. We believe that the resonance of the DOS peak in both electrodes is the main reason that the I–V characteristic of borophene–MoS2 system has a non-linear behavior under smaller external bias. Figure 3(d) shows the I–V curve of borophene–MoS2 system with a smaller bias interval of 0.02 V, and a piecewise function is used to fit the curve as: I1 =−935.83 V 2 + 629.89 V (V≤0.14 V), and I2 = 423.46 V + 10.71 (V≥0.16 V). The crossover point of these two fitted curves is at V = 0.139 V, which is close to 0.143 V.
We now turn to understand why currents of the systems with gold electrodes are drastically reduced compared to the current values of the systems without gold electrodes. Ref. [
40] studied the ballistic conduction across the intermetallic interface between gold and tungsten, and concluded that the conduction is reduced compared to the conduction of pure gold or tungsten because of dissimilar electronic structure that gold has
s wave-like modes but tungsten has
d wave-like modes. Refs. [
41] and [
42] also think that one factor to influence the spin injection efficiency from magnetic metal to graphene or transition-metal dichalcogenide semiconductors is orbital overlaps. In our case, borophene is actually extended along the transport direction and gold acts as top contacts [
42], that is, the
p orbitals of borophene can contribute to the electron transport through the entire system, therefore, we still seek another understanding on the current reducing with gold electrodes.
Figure 4(a) shows the averaged electrostatic potential of four systems without molecules over the transverse direction. Clearly, for the pristine borophene and heterostructure systems, the potential is periodic along the transport direction, which is determined by the periodicity of the 2D crystal. For the systems with gold electrodes, a potential barrier is formed in the scattering region due to charge transfer. In this way, the incoming electron should pass through this potential barrier before it goes out through the right electrode. We also made a charge analysis as shown in Fig. 4(b), in which the valence charge values of several B atoms in a selected row as indicted in Fig. 4(c) is plotted. For pristine borophene the charge is +3 on all the B atoms. After MoS2 substrate is added, electron transfer takes place from Mo atoms to B atoms, leading to slightly increase of the charge on B atoms, these charges contribute to the energy bands of the heterostructure around the Fermi energy as shown in Fig. 3(c). After gold electrodes are introduced, more electrons transfer from the electrodes to B atoms, leading to potential barrier forming in the scattering region. Although gold is a good conductor and borophene also shows metallic state, this potential barrier hinders the electron transport and results in the decreasing of the current values.
We finally calculated the sensitivity
S of the borophene-based nano gas sensor to various molecules (Fig. 5), which is defined as
S = |
I0−
I|/
I, where
I and
I0 are the currents with and without molecules, respectively [
28,
43]. The sensitivity of the pristine borophene gas sensor keeps a high and nearly constant value (0.22), the sensitivity of the borophene–MoS
2 heterostructure gas sensor quickly and smoothly reduced with increasing of the bias voltage. After gold electrodes are added, the sensitivities keep below 0.2. We should point that under 0.1 V bias voltage the sensitivities of heterostructure gas sensor with gold electrodes (purple curves) can be clearly distinguished, which means that the four molecules we studied in this paper can be well detected and differed by the gas sensor.
In summary, we studied the quantum transport properties of borophene-based nano gas sensors with or without considering the influence of gas molecules, MoS2 substrate and/or metal electrodes using density functional theory combined with nonequilibrium Green’s function method. Numerical results show that the gas sensors with gold electrodes can well detected and distinguished CO, NO, NO2 and NH3 gas molecules. The effects of MoS2 substrate and gold electrodes to the transport properties of borophene are calculated and analyzed by electronic structure, charge, and electrostatic potential. Our studies provide useful information to understand the transport behaviors between 2D metallic materials with metal electrodes, which is usually not considered when theoretically investigating the transport properties of 2D metallic materials from first-principle in the literature.
See the supplementary material for (i) the optimized atomic details of borophene and MoS2 heterostructure; (ii) Relaxed structures of borophene monolayer after adsorption of four different gas molecules; (iii) The I–V curve of ten different systems; (iv) transmission for eleven different systems with or without considering the influence of gas molecules, MoS2 substrate or gold electrodes at equilibrium state and 0.5 V bias voltage and the transmission versus moment k at Fermi energy.