Hybrid MXene-Graphene/Hexagonal Boron Nitride Structures: Electronic and Molecular Adsorption Properties.

Recent advances in experimental techniques allow for the fabrication of hybrid structures. Here, we study the electronic and molecular adsorption properties of the graphene (G)/hexagonal boron nitride (h-BN)-MXenes (Mo2C) hybrid nanosheets. We use first-principles calculations to explore the structure and electronic properties of the hybrid structures of G-2H-Mo2C and h-BN-2H-Mo2C with two different oxygen terminations of the Mo2C surface. The embedding of G or h-BN patches creates structural defects at the patch-Mo2C border and adds new states in the vicinity of the Fermi energy. Since this can be utilized for molecular adsorption and/or sensing, we investigate the ability of the G-M-O1 and BN-M-O1 hybrid structures to adsorb twelve molecules. Generally, the adsorption on the hybrid systems is significantly higher than on the pristine systems, except for N2 and H2, which are weakly adsorbed on all systems. We find that OH, NO, NO2, and SO2 are chemisorbed on the hybrid systems. COOH may be chemisorbed, or it may dissociate depending on its location at the edge between the G/h-BN and the MXene. NH3 is chemisorbed/physisorbed on the BN/G-M-O1 systems. CO, H2S, CO2, and CH4 are physisorbed on the hybrid systems. Our results indicate that the studied hybrid systems can be used for molecular filtration/sensing and catalysis.

1. Introduction

Research on two-dimensional (2D) materials has been growing rapidly since the successful exfoliation of graphene (G) [1]. The reduced dimensionality changes a material’s behavior compared to its bulk properties, as we see in the transition from the metallic behavior of graphite (3D) to the semimetallic behavior of graphene (2D), or to the semiconducting or metallic behavior of carbon nanotubes (1D). Such low-dimensional systems generate a lot of interest from both fundamental research and application points of view [2]. 2D materials have attracted considerable attention during the past few years because of their unique physical and chemical properties and various potential applications [2,3]. For example, these materials may be used in ultra-thin integrated circuitry required for modern electronics, optoelectronics, sensors, catalysis, and energy storage devices [4]. 2D world materials can be divided into single-atom nanosheets (e.g., silicene [5], and germanene [6]), two-atom nanosheets (e.g., hexagonal boron nitride (h-BN) [7], CN [8], BeN [9], BeO [10], and transition metal chalcogenides [11,12,13]), and three or more-atom nanosheets (e.g., BCN [14,15], ZnSnN [16], MoSiN [17,18,19], and transition metal carbides and nitrides [20]).

MXenes are 2D transition metal carbides and nitrides [20,21,22]. Their chemical compositions are given by the formula MA (n = 1,2,3), where M is transition metal, and A can be carbon or nitrogen [23]. MXenes have been proposed for energy storage [24], sensors [25], and many other applications [23,26,27]. The surfaces of MXene layers may be chemically terminated/functionalized with certain groups, such as O, F, or OH to avoid dangling bonds (which changes their general formula to MAX, where X represents groups 16 and 17 elements). The tunable electronic properties, excellent ionic conductivity, high Young’s moduli, and the flexible layer spacing of MXenes make them ideal for many applications [28,29].

A heterostructure comprising two or more dissimilar 2D materials have very different properties from its component materials. Heterostructures are thus emerging as new exciting candidates in materials research [2]. They can exist in layered form, where layers of different 2D materials are held together by van der Waals (vdW) forces. They can also be formed when different 2D materials are combined in the same monolayer, forming a lateral 2D heterostructure [30], or an embedded 2D heterostructure [31]. Heterostructures may exhibit exciting electronic properties and a high potential for use in nanodevices [32,33,34]. The MXene TiCX (X = F, O, and OH)/graphene (G) heterostructures have been successfully fabricated and can be used as electrodes for Li- and Mg-ion batteries [35,36,37]. TiCT/h-BN/TiCT heterostructures have been studied as a supercapacitor material [38]. Porous MXene-Polyimide/nanocellulose hybrids have been used for electromagnetic interference shielding [39,40].

The interaction between G and MXenes has been investigated using first-principles calculations for the cases TiCX with different terminal groups, X=O, OH, and F [41]. G with MXene monolayers (MCX where M = Sc, Ti, V and X = OH, O) heterostructures have been studied for Li-battery applications [42]. Furthermore, the G-CN/MCO (M = Hf, Zr, and Sc) heterostructures not only exhibit favorable structure stability but also show potential for photocatalytic degradation of volatile organic compounds [43]. The TiCX MXene/G hybrid structure has been synthesized and used for NH gas sensing [44] and for highly efficient electrocatalytic hydrogen evolution [45]. The hybrids of MoC anchored on G sheets have been studied as anode materials for lithium-ion lithium storage [46].

In this work, we use first-principles calculations to investigate the properties of hybrid structures made from G and -BN patches inside monolayer 2H-MoC. We study their structural and electronic properties with two different oxygen terminations of the MoC layer [3,23]. We also explore the adsorption of 12 molecules (H, N, CO, OH, NO, CO, NO, HS, SO, NH, COOH and CH) on the hybrid structures.

2. Computational Method

First-principle calculations based on density functional theory (DFT) are used within the projector-augmented wave method [47], which is implemented in the Quantum Espresso package [48]. The exchange-correlation interaction is described by the generalized gradient approximation (GGA) via the Perdew-Burke-Ernzerhof (PBE) functional [49]. A 45 Ry energy cut-off is utilized, and a 6 × 6 2H-MoC supercell is used as the base MXene system, with a vacuum space of 20 Å to avoid any interaction between periodic images. All systems are fully relaxed (volume and ionic position) until the forces on the atoms become less than 0.001 Ry/Bohr. Van der Waals interactions (vdW-DF) [50,51] are included in our study. Brillouin zone integrations are performed using the tetrahedron method with Blöchl corrections and a 9 × 9 × 1 k-point grid to calculate the density of states. Charge transfer from/to the MXene systems are determined by calculating the Löwdin charges. Initial spin-polarized calculations show that all considered structures are non-magnetic.

3. Results and Discussions

To properly set the scene for our investigation of the hybrid structures, we first calculate the properties of single-layer MXenes (Figure 1). The single-layer MXene structure can be viewed as a 2D triangular lattice of X atoms sandwiched between two triangular lattices of M atoms, similar to single-layer transition metal dichalcogenides [52,53]. We investigate multiple MoC structures: pristine (M), G-MoC hybrid (G-M), and h-BN-MoC hybrid (BN-M), each with two different oxygen terminations (Figure 1a). The systems with the first termination (O1) have an oxygen atom atop the Mo layer between 3 Mo atoms (above a carbon site), while systems with the second termination (O2) have the oxygen atom atop a Mo atom. According to prior work, functional groups seem to prefer hollow sites over metal-atom sites [26]. To assess the stability of the G-M-O1, G-M-O2, BN-M-O1, and BN-M-O2, we calculate their formation energies () using the formula: where and () are the total energies of the hybrid sheet and its isolated constituent atoms, respectively, with a total number of atoms . for the G-patch and for the h-BN-patch. and (z = Mo, C, O, B, and N) are the energy of the isolated atom and the number of the z atoms in the hybrid or patch. The calculated values are summarized in Table 1. As we see, the values of indicate that these structures are stable. Furthermore, we notice that the O1 structures are more stable than their O2 counterparts, as reported in previous work [54].

Figure 1

The top and side views of the unterminated (a) M, (b) G-M, and (c) BN-M systems. The projections of the terminations O1 and O2 are shown in (a). Mo, C, B, and N atoms are shown in gray, yellow, green, and pink, respectively.

Table 1

The calculated formation energy ( eV), lattice constant (Å), average bond length (Å) of Mo-Mo (), Mo-C (), Mo-O (), Mo-patch () and C-patch (), and the number of O-termination atoms ().

Sys.	E_form	Lattice Constant	d_1	d_2	d_3	d_4	d_5	N_O
M-O1	−9.16	6.07	2.84	2.17	2.05	-	-	72
G-M-O1	−9.11	6.07	2.83	2.17	2.08	2.26	2.35	66
BN-M-O1	−9.09	6.07	2.85	2.20	2.04	2.26	1.73	70
M-O2	−8.30	6.07	2.93	2.23	1.72	-	-	72
G-M-O2	−8.82	6.07	2.86	2.18	1.95	2.20	1.44	66
BN-M-O2	−8.30	6.07	2.86	2.24	1.72	2.40	1.65	70

In the M-O1 system, we find the lattice constant of the supercell, and the bond lengths Mo-Mo, Mo-C, and Mo-O are 6.07 Å, 2.84 Å, 2.17 Å and 2.05 Å, respectively, (see Table 1). The calculated values of the lattice constant and bond lengths agree well with those in previous reports (lattice constant = 5.64–6.00 Å, Mo-Mo = 2.73 Å, Mo-C = 2.14 Å, Mo-0 = 2.06 Å) [3,23,55,56]. The electronic properties of the pristine M-O1 show that this mono-layered 2D material is metallic. The total and projected densities of states (DOS and PDOS, respectively) are shown in Figure 2a. At the Fermi energy and above it, the Mo-4d states dominate the spectrum. The O-2p states contribute from −0.25 eV below the Fermi energy to 2 eV above it. The C-2p states contribute far from the Fermi energy in the range (−2.0,−3.0) eV and above 2.0 eV. The DOS/PDOS spectra agree with recent literature [55]. The MXenes provide more 4d electrons because of the two metal atoms in the primitive cell. Therefore, pristine MoC are metallic structures with states near the Fermi energy, making them superior candidates for molecular adsorption.

Figure 2

The optimized structures with O termination, top and side views, and DOS/PDOS of (a) M-O1, (b) G-M-O1, (c) BN-M-O1, (d) M-O2, (e) G-M-O2, and (f) BN-M-O2. Mo, C, O, B, and N atoms are shown in gray, yellow, red, green, and pink, respectively.

We now discuss our hybrid systems (G and h-BN inside the M-O1, Figure 2b,c. The embedded graphene comprises a patch made up of 6 carbon atoms inside M-O1 after removing 2 Mo, 6 O atoms, and 6 C atoms, then adding a carbon hexagon. Structural relaxation gives us a lattice constant of 6.07 Å similar to that of M-O1. The hybrid G-M-O1 structure has a minor symmetric distortion at the Mo sites surrounding the patch where the bond lengths slightly change (see Table 1), and the distance between Mo-patch ( and C-patch () is 2.26 Å and 2.35 Å, respectively. The atoms moved away from their original positions to release the stress in the system, leading to the slight increase or decrease of the bond length as shown in Figure 2b. The DOS of G-M-O1 slightly shifts towards low energy as compared to the pristine structure. The contribution of C states is still tiny as compared to other states, such as in the energy range from −1.9 eV to 3.0 eV. The DOS in the energy range from 0.2 eV to 2.7 eV has many peaks because of the effect of the carbon ring and Mo atoms around the patch. The DOS at the Fermi energy is higher than the corresponding value for the pristine structure, which indicates that the G patch increases the conductivity of M-O1.

The BN-M-O1 system is constructed by removing 2 Mo, 6 O, and 3 C atoms, then adding a BN hexagon (Figure 2c). After structural relaxation, the lattice constant does not change. The bond lengths slightly change as compared to the previous structures to release the stress in the system (Table 1), and and are 2.26 Å and 1.73 Å, respectively. The DOS shifts towards the low energy, and the DOS in the energy range from 0.22 eV to 2.7 eV has many peaks due to the contribution of the B and N state effect of Mo atoms around the patch.

Turning to the pristine M-O2 system, Figure 2d, the lattice constant of a 6 × 6 supercell is similar to that of M-O1. The system has 72 O termination atoms. The Mo-Mo and Mo-C bond lengths increase to 2.93 Å, and 2.23 Å, respectively, while the Mo-O bond length decreases to 1.72 Å. All calculated values of the lattice constants and bond lengths agree well with those in previous reports (Mo-C = 2–2.5 Å, Mo-O = 1.9–2.3 Å) [3,23,56]. The DOS shows the metallic behavior of the M-O2 sheet. Mo states dominate from −2.8 eV to 3.0 eV, with some contributions from the O states. The C states have large contributions from −1.8 eV to −0.8 eV. The contribution below −2.0 eV is small as compared to the corresponding states in pristine M-O1. The DOS in the energy range from 0.2 eV to 2.7 eV is not constant as compared to the DOS of M-O1.

The relaxed G-M-O2 structure, Figure 2e, exhibits more distortion than the G-M-O1 structure, especially at the interface between the G patch and the M-O2. The lattice constant of hybrid-O2 structures does not change as compared to the lattice constant of M-O2 system. The bond lengths Mo-Mo, and Mo-C decrease to 2.86 Å, 2.18 Å, but that of Mo-O increases to 1.95 Å as compared to M-O2 system. The and are 2.20 Å and 1.44 Å, respectively, less than the corresponding values in G-M-O1 structures. The graphene patch increases the density of states near the Fermi energy. The contributions of O states and C states are small as compared to Mo states and the corresponding states in the pristine structure. The BN-M-O2 structure does not suffer from any significant deformation after structural optimization (Figure 2f). The bond length of Mo–Mo (around the h-BN patch) decreases to 2.86 Å, and Mo-C increases to 2.24 Å, but Mo-O does not change as compared to M-O2 structure. The and are 2.40 Å and 1.65 Å, respectively, (Table 1). The DOS of BN-M-O2 is very similar to that of the pristine structure. It has been reported that the M-O2 structure is less stable than M-O1 [54]. Furthermore, the large distortions observed in the G-M-O2 structure may render the molecular adsorption on it surface highly dependent on the adsorption site. Therefore, we limit our molecular adsorption studies to the G/BN-M-O1 structures [54].

Gas Adsorption on MXenes

We started the structural relaxation of the adsorbent-sheet systems on three distinct adsorption sites. On the pristine structure, we placed the adsorbent close to an oxygen atom (position O1). In the hybrid structures, we have two distinct starting points: above a C/N atom of the (G/h-BN hexagon (points C1/N1), and above the Mo atom at the edge of the patch (point Mo1, Figure 3c,e).

Figure 3

Adsorption locations of various molecules on (a,b) pristine M-O1, (c,d) G-M-O1, and (e,f) BN-M-O1. C and (B,N) refer to the center of the patch.Mo, C, O, B, and N atoms are shown in gray, yellow, red, green, and pink, respectively.

We consider 12 small molecules, and divide them into four groups, diatomic molecules (OH, NO, CO, N, H), triatomic molecules (NO, HS, SO, CO), tetratomic molecules (COOH, NH), and polyatomic molecules (CH). The adsorption energy is calculated using the formula: where is the total energy of the sheet with the adsorbate, is the total energy of the pristine or hybrid sheet, and is the total energy of the isolated molecule.

In order to find the minimum configuration for the adsorbent, the diatomic molecules are initially placed at (O1) (Figure 3a). After geometrical optimization, the molecules do not change their positions. By analyzing the structural data for diatomic molecules, We found that the OH bond length of the isolated OH molecule is 0.98 Å, which is close to the literature value of (0.97 Å [57]). It slightly increases to 0.99 Å after adsorption on all the pristine and hybrid sheets. The bond lengths of the isolated NO, CO are 1.16 Å, and 1.14 Å, respectively, which match previous reports (1.15 Å,and 1.13 Å, respectively, [58,59]). The NO bond slightly decreases to 1.12 Å (Figure 4a) on M-O1 and increases to 1.20 Å on G/BN-M-O1 (Figure 4a). The CO bond is the same as for the isolated CO (1.13 Å), and increases to 1.17 Å on G/BN-M-O1. The bond lengths of the isolated N and H are 1.10 Å and 0.75 Å, respectively, which are close to those in previous reports (1.09 Å and 0.74 Å, respectively, [58]). These become 1.11 Å and 0.75 Å, respectively, on M-O1, and 1.11 Å and 0.76 Å, respectively, on G/BN-M-O1.

Figure 4

The optimized structures of (a) NO, (b) HS, (c) COOH, and (d) NH adsorbed molecules on M-O1, G-M-O1, and BN-M-O1 nanosheets in different side views. Mo, C, O, B, and N, S, and H atoms are shown in gray, yellow, red, green, pink, black, and blue, respectively.

The adsorption energy of the most stable configurations was calculated by Equation (2). The detailed results of all the studied cases are presented in Table 2. Energies up to ∼2 eV indicate physisorption, while higher values indicate chemisorption [60] Calculations of and Q (e) of diatomic molecules on M-O1 show that OH, NO, and CO have the largest energies, which can be explained by the polarity of oxide molecules. Thus, OH and NO are chemisorbed on M-O1 while CO, N, and H are physisorbed. The OH molecule stands at a distance of 1.45 Å above the M-O1 surface, with its O atom facing the sheet. For NO and CO, N and C are the atoms of contact with the sheet, respectively. As for the bond distance of NO (O-N), it is 2.05 Å, while for CO the distance (O-C) is 4.17 Å.

Table 2

(X-X) refers to the adsorbed site (before-after) relaxation, E is the adsorption energy (eV), the closest distance between the sheet and molecule is d(Y-Z) (Å), where Y and Z indicate the sheet atom and adsorbed-molecule atom, respectively. Q (e) is the charge transfer, which is negative when a charge is transferred from the structure to a molecule. The references in the adsorption energy columns refer to the corresponding values of the gas adsorption on G, and h-BN monolayers.

Gas	M-O1				G-M-O1				BN-M-O1
	X _i -X _f	E _ad	Δ Q	d(Y_s-Z_m)	X _i -X _f	E _ad	Δ Q	d(Y _s -Z _m )	X _i -X _f	E _ad	Δ Q	d(Y_s-Z_m)
OH	O1-O1	−1.15	−0.21	1.45 (O-O_m)	C1-C1	−5.48/−0.5 [61]	0.07	2.11 (Mo-O_m)	N1-N1	−5.45	0.07/−2.3 [61]	2.11 (Mo-O_m)
	–	–	–	–	Mo1-C1	−5.48	0.06	2.11 (Mo-O_m)	Mo1-N1	−5.45	0.07	2.11 (Mo-O_m)
NO	O1-O1	−1.10	−0.56	2.05 (O-N_m)	C1-C1	−2.58	−0.06	2.07 (Mo-N_m)	N1-N1	−3.15/−0.026 [62]	−0.06	2.09 (Mo-N_m)
	–	–	–	–	Mo1-C1	−2.58	−0.07	2.07 (Mo-N_m)	Mo1-N1	−3.15	−0.07	2.09 (Mo-N_m)
CO	O1-O1	−0.10	−0.07	4.17 (O-C_m)	C1-C1	−1.26/−1.13 [63]	−0.21	2.21 (Mo-C_m)	N1-N1	−1.81/−0.02 [62]	−0.21	2.21 (Mo-C_m)
	–	–	–	–	Mo1-C1	−1.27	−0.19	2.21 (Mo-C_m)	Mo1-N1	−1.81	−0.19	2.22 (Mo-C_m)
N_2	O1-O1	−0.09	−0.07	3.08 (O-N_m)	C1-C1	−0.23	−0.20	3.12 (Mo-N_m)	N1-N1	−0.64	−0.21	3.28 (O-N_m)
	–	–	–	–	Mo1-C1	−0.24	−0.17	3.03 (Mo-N_m)	Mo1-N1	−0.64	−0.18	3.26 (Mo-N_m)
H_2	O1-O1	−0.05	−0.014	2.46 (O-H_m)	C1-C1	−0.07/−0.08 [64]	−0.02	3.58 (O-H_m)	N1-N1	−0.09/−0.21 [65]	−0.02	2.90 (Mo-H_m)
	–	–	–	–	Mo1-Mo1	−0.08	−0.06	2.69 (O-H_m)	Mo1-Mo1	−0.06	−0.02	2.74 (O-H_m)

For adsorption on G-M-O1 hybrid, all the diatomic molecules starting from the C1 position end up in the same place after structural optimization (C1 → C1). When we start the relaxation from the Mo1 position, the diatomic molecules move to C1 (Mo1 → C1) except the H molecule, which stays in Mo1 (Mo1 → Mo1), (Figure 3c,d).

For the BN-M-O1 system, all adsorbed diatomic molecules starting from the N1 position do not change their positions after structural relaxation (N1 → N1). These molecules at the Mo1 position move to the N1 position after optimization (Mo1 → N1), except for H, which stays at Mo1 (Mo1 → Mo1) (Figure 3e,f). The available theoretical and experimental studies show that the adsorption of diatomic molecules on pristine (G) [61,63,64] and pristine h-BN)[61,65,66,67] are smaller than the corresponding value on the hybrid sheets, except for H molecule, (see Table 2). The E in the case of the G or h-BN hybrid is much higher than in that of M-O1, pristine G, and pristine h-BN for all molecules. This is because of electrostatic irregularities at the border between the patch and the original sheet. The molecules are adsorbed at the same atom as in the case of M-O1, and the distance between the molecules and the sheet becomes larger as compared to the M-O1 sheet except for CO (see Table 2).

Regarding triatomic molecules, the initial and final positions of the adsorbed molecules on the different sheets are shown in (Figure 3) and listed in Table 3. For the M-O1 system, and for all triatomic molecules, we find that the final position is the same as the initial position (Figure 3a,b). For the hybrid systems, an NO molecule starting at the C1 or Mo1 positions ends up at the midpoint between C1 and C6 (C1 → C1-C6). The N–O bond length and the O–N–O angle of the isolated molecule are 1.21 Å and 133.5, respectively, which agree with values reported in the literature (1.20 Å and 134.3[68]). After adsorption, these change to 1.17 Å and 145.62 for the pristine sheet, and 1.26 Å and 118.25 for G-M-O1 and BN-M-O1 sheets. For SO, the S–O bond length and the O–S–O angle are 1.45 Å and 119.5, respectively, which also agree with the corresponding values in the literature (1.43 Å and 119[69]). These become 1.45 Å and 118.57 on the pristine sheet, and 1.62 Å and 108.97 on the hybrid systems G-M-O1 and BN-M-O1.

Table 3

X-X refers to the adsorbed site before-after relaxation, E is the adsorption energy (eV). The closest distance between the sheet and an adsorbent is d(Y-Z) (Å), where Y and Z indicate the sheet atom and adsorbent atom, respectively, and Q (e) refers the charge transfer where the negative sign means charge is transferred from the structure to an adsorbent. C and (B,N) refer to the final position located at the hexagonal ring of G and h-BN, respectively. The references in the adsorption energy columns refer to the corresponding values of the adsorbent gases on the pristine G and h-BN monolayers.

Gas	M-O1				G-M-O1				BN-M-O1
	X _i -X _f	E _ad	Δ Q	d (Y_s-Z_m)	X _i -X _f	E _ad	Δ Q	d (Y _s -Z _m )	X _i -X _f	E _ad	Δ Q	d (Y _s -Z _m )
NO_2	O1-O1	−0.58	−0.59	1.97 (O-N_m)	C1-(C1,C6)	−2.46/−2.17 [71]	−0.03	2.21 (Mo-N_m)	N1-(B1,N1)	−2.91/−0.03 [62]	−0.03	2.21 (Mo-N_m)
	–	–	–	–	Mo1-(C1,C6)	−2.46	−0.04	2.21 (Mo-N_m)	Mo1-(B1,N1)	−2.91	−0.04	2.20 (Mo-N_m)
SO_2	O1-O1	−0.49	−0.06	3.66 (O-O_m)	C1-C_H	−2.31/−1.97 [72]	−0.12	2.13 (Mo-O_m)	N1-(B,N)_H	−2.78/−0.08 [67]	−0.13	2.14 (Mo-O_m)
	–	–	–	–	Mo1-C_H	−2.28	−0.13	2.15 (Mo-O_m)	Mo1-N1	−1.93	−0.24	2.15 (Mo-O_m)
H_2S	O1-O1	−0.26	−0.27	3.46 (O-H_m)	C1-C1	−1.19/−0.45 [73]	−0.61	2.78 (Mo-S_m)	N1-N1	−1.65/−0.11 [74]	−0.59	2.69 (Mo-S_m)
	–	–	–	–	Mo1-C1	−1.22	−0.62	2.72 (Mo-S_m)	Mo1-N1	−1.65	−0.61	2.71 (Mo-S_m)
CO_2	O1-O1	−0.19	−0.06	2.80 (O-C_m)	C1-C3	−0.56/-0.32 [75]	−0.32	2.54 (Mo-O_m)	N1-B2	−0.78/−0.02 [62]	−0.05	2.95 (Mo-O_m)
	–	–	–	–	Mo1-Mo1	−0.12	−0.22	3.07 (O-O_m)	Mo1-Mo1	−0.55	−0.05	3.05 (O-O_m)

The changes occurring in the HS and CO molecules are noticeably less. The calculated H–S bond length and the H–S–H angle are 1.35 Å and 61.61 on the M-O1 sheet (Figure 4b), and 1.37 Å and 94.52 for the hybrid sheets (Figure 4b), compared to 1.35 Å and 91.7 for the isolated molecule (literature values are 1.34 Å and 92.1[70]). For CO, we find that C–O bond length and the O–C–O angle to be 1.17 Å and 179.86 for the pristine sheet, and 1.45 Å and 178.47 for the hybrid sheets, compared to 1.17 Å and 180 for the free molecule (literature values are 1.16 Å and 180 [68]).

We notice that in the 4 cases, the bond length increases upon adsorption on the hybrid systems due to the redistribution of the electronic charge. For NO and SO, we see large adsorption energies and small charge transfers, indicating that the bonds formed are more covalent in nature, which can be described as chemisorption. For HS and CO, the energies are small, but we see a moderate charge transfer, with longer bonds with the hybrid Mo atom, indicating physisorption. We also notice that the adsorption occurs through the N and O atoms for the case of NO and SO, respectively, while for HS (Figure 4b) and CO, it occurs through the H and C atoms, respectively. This can be explained by the difference in electronegativity between the adsorbent atom and the O atom of the sheet. Table 3 shows that NO has the highest adsorption energy for all sheets. Compared to the G, h-BN, and M-O1 sheets, the hybrid sheets are superior in adsorbing the studied triatomic molecules, with BN-M-O1 offering the highest adsorption. The adsorption energies and the distances between various adsorbents and the underlying systems are shown in Table 3.

Now we consider the adsorption of the tetratomic molecules COOH, and NH, and the polyatomic molecule CH on the pristine and hybrid systems. The COOH molecule is physisorbed on the M-O1 sheet at the C1 location, while it is chemisorbed with a larger energy on the hybrid systems at the same location. If the molecule is positioned at the Mo1 location on either hybrid sheet, it dissociates into a CO molecule and a H atom that gets attached to a sheet O atom (Figure 4c). The bond lengths and the angles of the isolated NH and CH are 1.02Å and 106.67, and 1.09 Å and 110.41, respectively, which match those of previous reports (1.00 Å and 109.47, and 1.09 Å and 109.5, respectively, [76,77]). The bond lengths and the angles of both molecules do not change significantly after adsorption. The optimized location for NH is at C1/N1 on G-M-O1/BN-M-O1. For CH, the molecule stays at its initial locations: C1/Mo1 on G-M-O1, and Mo1 on BN-M-O1.

As expected, the adsorption energy of COOH is larger than those of NH and CH due to the structure of the molecule which allows for an easier redistribution of charge (Figure 4c). This translates to a closer distance d(Y-Z) between the molecule and the hybrid systems (Table 4). As before, and for the three molecules, the adsorption on the hybrid systems is higher than on the pristine system.

Table 4

(X-X) refers to the adsorbed site (before-after) relaxation, E is the adsorption energy (eV). The closest distance between the sheet and molecule is d(Y-Z) (Å) where Y and Z indicate the sheet atom and the adsorbed molecule atom, respectively, and Q (e) refers the charge transfer where a negative sign means the charge is transferred from the structure to a molecule. C and (B,N) mean the final position is at the hexagonal ring of G and h-BN, respectively. The “D” indicates the dissociation of the molecule. The references in the adsorption energy columns refer to the corresponding values of the adsorbents on the pristine G and h-BN monolayers.

Gas	M-O1				G-M-O1				BN-M-O1
	X _i -X _f	E _ad	Δ Q	d (Y_s-Z_m)	X _i -X _f	E _ad	Δ Q	d (Y_s-Z_m)	X _i -X _f	E _ad	Δ Q	d (Y _s -Z _m )
COOH	O1-O1	−2.95	−0.45	1.54 (O-C_m)	C1-C1	−3.32	−0.24	2.26 (Mo-O_m)	N1-N1	−3.11	−0.53	2.41 (Mo-C_m)
	–	–	–	–	Mo1-d	D	D	D	Mo1-d	D	D	D
NH_3	O1-O1	−0.40	−0.30	2.66 (O-N_m)	C1-C1	−1.78/−1.09 [71]	−0.45	2.95 (O-N_m)	N1-N1	−2.18/−0.17 [78]	−0.45	2.80 (O-N_m)
	–	–	–	–	Mo1-C1	−1.78	−0.46	2.76 (Mo-N_m)	Mo1-N1	−2.66	−0.46	2.75 (Mo-N_m)
CH_4	O1-O1	−0.21	−0.03	2.81 (O-H_m)	C1-C1	−0.41/−0.33 [75]	−0.13	2.36 (O-H_m)	N1-N1	−0.80/−0.12 [79]	−0.03	2.40 (O-H_m)
	–	–	–	–	Mo1-Mo1	−0.24	−0.13	2.60 (O-H_m)	Mo1-Mo1	−0.66	−0.04	2.60 (O-H_m)

4. Conclusions

In this paper, we use density functional theory to investigate the structure and electronic properties of the hybrid structures of graphene (G)/2H-MoC (G-M), and hexagonal boron nitride ( h-BN)/2H-MoC (BN-M), with two different O-functionalizations (O1 and O2). The G-M-O1/O2 and BN-M-O1/O2 hybrid sheets are metallic as 2H-MXenes (M-O1 and M-O2). The G and h-BN patches create localized states near the Fermi energy, which can be utilized for the adsorption/sensing of some molecules. We thus investigate the ability of G-M-O1 and BN-M-O1 hybrid structures to adsorb twelve molecules: diatomic (OH, NO, CO, N, H), triatomic (NO, SO, HS, CO), and polyatomic (COOH, NH, and CH). Generally, the change in the bond lengths and angles of the molecules upon adsorption are larger on the hybrid systems compared to the pristine system (except for H, which is weakly adsorbed on all systems). This is primarily due to the structural deformations occurring at the borders between the G/ h-BN and the MXene regions, which act as electrostatic/covalent traps. We found that OH, NO, NO, and SO are chemisorbed on G-M-O1 and BN-M-O1. COOH may be chemisorbed, or it may dissociate depending on its location at the edge between the G/ h-BN and the MXene. NH is chemisorbed/physisorbed on the BN/G-M-O1 systems. The other studied molecules are all physisorbed on the hybrid systems. The results of our work indicate that the hybrid nanosheets can be used in catalysis and molecular sensing.