M. Kamaraj,
J. Vijaya Sundar and
V. Subramanian*
Chemical Laboratory, CSIR – Central Leather Research Institute, Adyar, Chennai 600020, India. E-mail: subuchem@hotmail.com; subbu@clri.res.in; Fax: +91 44 24911589; Tel: +91 44 24411630
First published on 31st October 2016
The changes in the electronic properties of single and bilayers of graphene (G) and hexagonal boron nitride (h-BN) two dimensional (2D) sheets have been investigated upon interaction with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) by employing density functional theory (DFT) based calculations. The calculated interaction energy, band gap and charge transfer reveal that bilayer h-BN (BLBN) may serve as a better potential candidate for sensing TCDD than the other systems. To gain further insight into the sensing properties of these materials, the transmission spectra and current–voltage (I–V) characteristics have been calculated by using Non-Equilibrium Green's Function (NEGF) combined with DFT approach. It is interesting to mention that a Negative Differential Resistance (NDR) effect has been observed in a single layer BN (SLBN) nanoribbon upon interaction with TCDD. The calculated I–V characteristics of the BLBN nanoribbon–TCDD complex reveal that the interaction between the two systems enhances the current flow through the BLBN nanoribbon, which is appreciably higher than that of pristine BLBN. These findings may open up new avenues for the application of bilayers as possible sensors for pollutants.
Dioxin families are the most toxic organic pollutants. They are stable pollutants in water, soil, and air.40 It is lipophilic in nature.41 Dioxins are released during many industrial processes as an unwanted by-product. In some natural process such as volcanic eruptions and forest fire, dioxins are released.42 Any source of organic materials in the presence of chlorine or other halogens can generate dioxins and furans during combustion.41 The half-life of dioxin in soil and sediments is estimated to be 50–100 years.43 Among the members of dioxin family, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is highly toxic due to its half-life (10 years) in the human body.44,45 TCDD is a well-known carcinogen. The exposure of TCDD leads to skin problem called chloracne. Also, it affects the immune system and causes severe reproductive and developmental problems. Therefore, various experimental and theoretical studies have been carried out to sense TCDD with the help of new functional materials. Various nano materials are widely used in sensing of TCDD. Some of the significant studies on this topic are covered in the following section.
Single wall carbon nanotube (SWCNT) has been used as an efficient material for the removal of dioxin from waste incinerators.46 By using first principle calculation Fagan et al. showed that creation of lattice defect on SWCNT can increase the interaction of TCDD when compared to pristine SWCNT. The calculated interaction energy of the complex is −1.21 eV.47 The interaction of CNT and BNNT with TCDD has been investigated and it is found that TCDD interacts strongly with BNNT than CNT.48 The calculated interaction energy of TCDD with BNNT is −0.55 eV. Metal-doped graphene also been used as an effective adsorbent for TCDD.49,50 Ti doped graphene sheet shows interaction energy of −1.94 eV for TCDD and Ca decorated graphene has interaction energy of −3.20 eV. The dissociative interaction of dioxin on various metal oxide surfaces has also been studied and results suggested that the surface oxygen atoms of metal oxides with more nucleophilic character and less tightly bonded atoms might be promising sites for the decomposition of the dioxins and dibenzofurans.51
Experimentally, the sensing of environmental pollutants is a challenging task because of its high toxicity. Hence, computational chemistry approaches have been utilized to predict and design possible materials for the development of sensors. The specific surface to volume ratio of the nanomaterial is an important parameter in determining the sensitivity of the sensors. This value is much higher for nanomaterials based sensors than conventional micro sensors. Particularly, the 2D material provides the large adsorption area for sensing molecules as all of their atoms are exposed to the adsorbing molecules. Thus, highest sensing efficiency per unit volume can be achieved.21 The simple mechanism of the nanomaterial-based sensor is to detect the changes in the conductivity of material upon adsorption of analyte molecules. The changes in the conductivity could be attributed to the variations in the charge carrier concentration in the material as induced by adsorbed target molecules.
Since, G is a semimetal with zero band gap, and h-BN is wide band gap insulator, it is interesting to unravel the changes in the electronic properties of these materials upon surface interaction with TCDD. In this investigation, the following points have been addressed by using DFT based calculations combined with NEGF formalism.
(1) To understand the nature of interaction and binding modes of TCDD with G and h-BN single and their bilayers.
(2) To unravel the changes in the electronic structure, transport properties and I–V characteristics of G and h-BN single and their bilayer after interaction with TCDD.
(3) To predict a suitable material for sensing of TCDD based on the binding energy as well as the transport properties of single and bilayers of G and h-BN.
The transport properties were calculated using TransSIESTA code57 as implemented in SIESTA package.58 The computational details are presented in the ESI 1.† The LDA-DFT was previously used to obtain reliable results on lattice constants and binding energy of bilayers of G and h-BN systems.59,60 The same method was also yielded other electronic properties of these 2D systems in close agreement with the experimental results.61–64 The LDA-DFT method mimics features of the electronic band structure obtained by many-electron Green's function approach within the GW approximation reasonably well.65 LDA functional shows promising performance for noncovalent interactions between nanomaterials with various molecules.30,61–68 The choice of the LDA is not fortuitous and is more suitable than the GGA to study weakly interacting systems like the π-stacking interaction on sp2-like materials.69 It has been widely applied in the computational study of adsorption of the molecules on graphene and h-BN. Previous reports61–68 has also highlighted the accuracy of LDA functional. It should be worth to point out that the LDA-DFT method has been shown to provide reasonable good descriptions of the physics and chemistry of graphitic and hBN systems.61,62 Based on these solid reports, we have selected this functional in the present investigation.
As a benchmark, the geometrical and electronic properties of single layer graphene (SLG) (98C) and single layer of h-BN (SLBN) (49 B and 49 N) was calculated using 7 × 7 × 1 hexagonal supercell. The size of the supercell was chosen to keep the distance between TCDD molecules large enough to eliminate their interactions with neighbouring supercells. The supercell was constructed from the optimized unit cell with lattice constant of 2.44 Å for SLG and 2.48 Å for SLBN. The calculated C–C and B–N bond lengths are 1.41 and 1.45 Å, respectively. The calculated lattice constants (SLG – 2.44 Å, SLBN – 2.48 Å) and band gaps (SLG – 0 eV and SLBN – 4.39 eV) are consistent with previous results.70,71 Similarly, the bilayer graphene (BLG), bilayer h-BN (BLBN) and their hetero bilayers (GBN and BNG) were optimized. The calculated interlayer distances and band gap of these pristine homo and hetero bilayers are presented in ESI Table S1.† These values are also in good agreements with the previous reports.70,72,73
It has been found from previous studies13,73–76 that AB stacking arrangement in both homo and hetero bilayers was found to be the most stable one. The homo bilayer was constructed by placing the one G (h-BN) layer on the top of another G (h-BN) in the AB stacking arrangement. Similarly, the hetero bilayer system was constructed by stacking of G (h-BN) monolayer with the monolayer of h-BN (G). In this hetero bilayer arrangement, the carbon (boron) atom is on top of the boron (carbon) atom and the nitrogen atom is in the centre of hexagon of graphene. Thus, the AB stacking configuration was selected in this study for the construction of bilayers. The calculated mismatch in lattice constants between G and h-BN is less than 2%. The one-dimensional (1D) armchair nanoribbons were used for the calculation of transport properties (further details are included in ESI 2†).
The interaction energy (Eint) of TCDD adsorbed single and bilayers of G and h-BN were calculated by conventional definition Eint = Esub+ads − (Esub + Eads), where Eint, Esub+ads, Esub and Eads are interaction energy, total energy of the complexes, total energy of substrate and the total energy of adsorbate, respectively. The energetics of various binding modes of TCDD with single layer was calculated and it is found that AB stacking like pattern of TCDD is be the most stable arrangement. Hence, the same was chosen for the exploration of interaction of TCDD with bilayers.
Atom | Bond lengths (Å) | |
---|---|---|
This work | Previous experimental work77 | |
C2–H2 | 1.09 | — |
C1–Cl2 | 1.70 | 1.72 |
C1–C6 | 1.39 | 1.38 |
C6–C5 | 1.37 | 1.38 |
C5–C4 | 1.37 | 1.37 |
C4–C3 | 1.39 | 1.38 |
O2–C3 | 1.35 | 1.37 |
Both top and side view of the optimized geometries of all the complexes are shown in Fig. 2. It can be seen from Fig. 2(a) that the two aromatic rings of TCDD orient in AB stacking arrangement with reference to the basal plane of the graphene sheet. A similar kind of arrangement has also been observed to be the most stable one in the previous computational studies, where the aromatic molecules are stacked with graphene.78–80 Similarly, in the case of h-BN sheet, the carbon atoms (C1, C3, C5 and C8, C10, C12) in the aromatic ring of TCDD lie on the top of the boron atoms of sheet and the remaining three carbons lie at the centre of the hexagonal ring of the sheet as shown in Fig. 2(b).
![]() | ||
Fig. 2 The top and side view of optimized geometry of TCDD on (a) SLG (b) SLBN (c) BLG (d) BLBN (e) GBN and (f) BNG. |
The calculated interaction energy (Eint), band gap (Eg), interlayer distance (dI) and optimized average vertical distance between the TCDD and the top layer of the sheet (dTCDD–sheet) are provided in Table 2. The calculated interaction energy values imply that TCDD forms stable complex with single and bilayers of G and h-BN. It can be seen from Table 2 that the optimized average vertical distance between the TCDD and the SLG is 3.27 Å. In the case of BLG, the average vertical distance between TCDD and the top layer of BLG is 3.24 Å. TCDD interacts strongly with the surface of BLG when compared to the interaction of TCDD with SLG. The calculated interaction energy of TCDD with SLG at LDA (vdW-DF2) level of theory is −0.75 (−1.25) eV. Upon inclusion of additional layer, the interaction energy increases. The interaction energy of TCDD with BLG at LDA (vdW-DF2) level of calculation is found to be −1.17 (−2.69) eV. A similar trend of increasing the interaction energy upon addition of another layer has already been observed in the interaction of benzene with multilayer graphene.81
System | Eint (eV) | Eg (eV) | dI (Å) | dTCDD–sheet (Å) |
---|---|---|---|---|
a The values within the brackets was calculated with vdW-DF2 functional. | ||||
SLG–TCDD | −0.75 (−1.25) | 0 | — | 3.27 (3.36) |
BLG–TCDD | −1.17 (−2.69) | 0 | 3.33 (3.35) | 3.24 (3.26) |
SLBN–TCDD | −0.91 (−1.89) | 3.12 | — | 3.08 (3.30) |
BLBN–TCDD | −1.23 (−2.78) | 2.99 | 3.20 (3.21) | 3.06 (3.12) |
GBN–TCDD | −0.84 (−1.83) | 0.06 | 3.19 (3.22) | 3.20 (3.18) |
BNG–TCDD | −1.14 (−2.56) | 0.05 | 3.17 (3.19) | 3.13 (3.16) |
As expected, the weak vdW interaction plays the predominant role in the stabilization of the complexes. Hence, there is only a marginal change in the interlayer distance between two graphene sheets upon interaction with TCDD. The calculated band structures of SLG, SLG-TCDD and BLG-TCDD complexes are depicted in Fig. 3. The interaction of TCDD on both SLG and BLG does not significantly affect the electronic properties of graphene. The linear band near the Dirac point of pristine graphene survives after the interaction of TCDD on graphene. The presence of Dirac point and associated band structure clearly reveals that the stacking interaction does not alter the D6h symmetry of the graphene.
![]() | ||
Fig. 3 Band structures of (a) pristine SLG, (b) SLG–TCDD and (c) BLG–TCDD. The Fermi level was set to be zero. |
The partial density of states (PDOS) of both SLG and BLG after interaction with TCDD is shown in Fig. 4. The highest occupied and lowest unoccupied states of TCDD hybridize with the G orbitals. As a result, some hybrid peaks are observed in PDOS of both SLG and BLG sheets. Overall findings reveal that the interaction of TCDD does not alter the electronic structure of graphene. However, significant increase in the interaction energy can be observed in the BLG when compared to SLG. Lowdin population analysis reveals that, significant amount of charge (electron) is transferred from TCDD to graphene systems. The calculated charge transfer from TCDD to BLG is 0.22e and for SLG–TCDD system is found to be is 0.09e. The increase in charge transfer may be responsible for the enhanced interaction of TCDD with BLG.
![]() | ||
Fig. 4 Partial density of states (PDOS) of (a) pristine SLG, (b) SLG–TCDD and (c) BLG–TCDD. The Fermi level was set to be zero. |
The optimized geometries of SLBN–TCDD and BLBN–TCDD are presented in Fig. 2(b) and (d). The calculated interaction energies at LDA (vdw-DF2) level of theory for SLBN–TCDD and BLBN–TCDD complexes are −0.91 (−1.89) eV and −1.23 (−2.78) eV, respectively. The trend in the interaction energy for these systems is similar to that of SLG–TCDD and BLG–TCDD complexes. However, the calculated interaction energy for BLBN–TCDD complex is relatively higher than that of BLG–TCDD complexes by about ∼1.38 kcal mol−1. The optimized average distances between TCDD and SLBN as well as BLBN are 3.08 and 3.06 Å, respectively. It can be seen from the Table 2 that, the interaction distance between TCDD with BLBN is less than all other systems implying the strong interaction between TCDD and BN layers.
The band structures of SLBN, SLBN–TCDD and BLBN–TCDD are given in Fig. 5. It is well known that SLBN is a wide band gap semiconductor with a value of 4.39 eV. It is possible to observe from Table 2 that band gap of SLBN reduces to 3.12 eV after interaction with TCDD. Similarly, in the case of BLBN, the band gap decreases after interaction with TCDD from 4.32 eV to 2.99 eV. The strong interaction of TCDD introduces new flat bands in the band gap region of both SLBN and BLBN. The reduction of bandgap may be due to the recombination of the energy levels of BN sheet with those of TCDD. This is also reinforced by the PDOS of SLBN and BLBN as shown in Fig. 6. The PDOS indicates that both valence band and conduction bands originate from C2p orbitals of TCDD. From the Lowdin charge analysis, it is found that 0.20e is transferred from TCDD to SLBN. However, the same from TCDD to BLBN is 0.42e. This may be the reason for the highest interaction energy in case of BLBN systems. Thus, the interaction of TCDD with these layers affect the conductivity of the both SLBN and BLBN and thus concomitant enhancement in the electronic response can be exploited for the development of sensors.
![]() | ||
Fig. 5 Band structures of (a) pristine SLBN, (b) SLBN–TCDD and (c) BLBN–TCDD. The Fermi level was set to be zero. |
![]() | ||
Fig. 6 PDOS of (a) pristine SLBN, (b) SLBN–TCDD and (c) BLBN–TCDD. The Fermi level was set to be zero. |
The optimized structures of graphene–h-BN (GBN) and h-BN–graphene (BNG) with TCDD are presented in Fig. 2(e) and (f). The optimised distance between TCDD and GBN is found to be 3.20 Å. The same distance between TCDD and BNG is 3.13 Å. Among these hetero bilayers, the stronger interaction between TCDD with BNG is evident from these distances. The difference in the interaction energy of BNG–TCDD and GBN–TCDD complexes is about ∼6.92 kcal mol−1. The calculated band structures of GBN and BNG after interaction of TCDD are depicted in ESI Fig. S1.† The predicted bandgaps of both GBN and BNG bilayers are 0.063 and 0.053 eV, respectively. The interaction of TCDD does not influence band gap of hetero bilayers. The calculated charge transfers from TCDD to BNG as well as GBN are 0.31e and 0.25e, respectively. The calculated PDOS (ESI Fig. S2†) indicates that there are no changes near Fermi levels of both BNG and GBN after interaction with TCDD. Thus, electronic properties of these hetero bilayers remain unchanged after interaction of TCDD.
Results obtained from the calculation highlighted that BLBN–TCDD complex has higher interaction energy than rest of the complexes. Interaction energy is one of the important criteria to quantify the sensing ability of host materials. Among the six model (SLBN, SLG, BLG, BLBN, GBN and BNG) systems, BLBN nano sheet effectively sense TCDD. This stronger interaction allows adsorption of TCDD over BLBN for long time. It is well known that the parameter that determines the employment of these materials as a suitable substrate for sensing adsorbates is the recovery time.31 Hence a discussion on the recovery time has been qualitatively included in the following section.
The recovery time (τ) can be expressed by the conventional transition state theory31 as given below
τ = ν0−1e(−Eint/κT) |
In general, the selectivity of guest molecules (H2, No, NH2, CH3, TTF, TCNQ, etc.) interacting over 2D materials can be easily monitored by measuring the output current of pristine and complex materials. The interaction of various target molecules on 2D h-BN system shows different band gap. For example, the interaction of TCNQ on BN sheet decreases the band gap of BN sheet to 0.43 eV, and TTF to 2.07 eV. Similarly, the interaction of O atom on BN sheet decreases the band gap of the BN systems to 3.8 eV. From molecular level calculations, the interaction of para-nitrophenol on the BN nanosheet decreases the band gap of BN nanosheet from 5.93 eV to 3.85 eV.32 So interaction of different target molecules shows different values of band gap, so obviously the conductivity of these systems are different from TCDD interacting systems. In the present investigation, we highlighted that the adsorption of TCDD decreases the band gaps of BLBN and SLBN to 2.9 eV and 3.1 eV, these changes in the band gap will increase the electrical conductivity. Results reveal that, the interaction of TCDD over BLBN reduces the band gap of the sheet from 4.3 eV to 2.9 eV. Therefore, the electrical conductivity of the BLBN sheet increases. We have measured the output current of BLBN sheet before and after interaction of TCDD complex. It was pictorially represented in Section 5. So selectivity of TCDD can be monitored by measuring the output currents before and after interaction. It is found from previous report36 that, the small changes in the band gap of Al-doped BN sheet, influences the large different in electrical conductivity of the material.
Both transmission spectra and DOS for SLG nanoribbon before and after interaction with TCDD at zero bias are represented in Fig. 8. The step-like behaviour in transmission spectra is observed due to the conductance channels in the band structure. The interaction of TCDD with SLG nanoribbon leads to the appearance of dips in the transmission spectra. These are also in accordance with the DOS that contains a localized peak at −0.8 eV below the Fermi level. It can be noticed from Fig. 8 that the interaction of TCDD causes a decrease in the conductance of SLG nanoribbon. The decrease in the conduction can be attributed to the electronic perturbation in the nanoribbon after interaction of TCDD. This perturbation may lead to intrinsic backscattering along the sheet, which alters conduction manifolds and reduces the transmission capability of the system.84 The computed transmission spectra and DOS of BLG nanoribbon before and after interaction with TCDD are shown in Fig. 9. The step-like behaviour is maintained after interaction of TCDD in BLG nanoribbon with a small localized peak at −0.85 eV. The zero bias transmission spectra and DOS of SLBN nanoribbon before and after interaction of TCDD are depicted in Fig. 10. The disappearance of step-like behaviour can be noticed after interaction with TCDD. As observed from ESI Table S2,† pristine SLBN nanoribbon is a wide gap semiconductor with a band gap of 4.55 eV. After interaction of TCDD the band gap decreases to 3.14 eV. The reduction in band gap is reflected in transmission spectra and DOS. As evident from DOS, peak is observed at 1.82 eV above the Fermi level. Similarly, in case of BLBN nanoribbon, the band gap decreases from 4.34 eV to 3.02 eV after interaction with TCDD. The computed zero bias transmission spectra and DOS of BLBN nanoribbon are displayed in Fig. 11. It is clearly evident from Fig. 11(a) that, after interaction of TCDD on BLBN nanoribbon a localized peak can be found at −1.30 eV below the Fermi level. Above Fermi level, another peak can be seen at 2.43 eV. The interaction of TCDD is responsible for the increase in the conductance of bilayer BN system. The zero bias transmission spectra of hetero bilayers (both GBN and BNG) nanoribbons before and after interaction with TCDD are given in ESI Fig. S4(a) and S5(a).† The step-like behaviour is observed in both cases. However, in the case of GBN system, there are no appreciable changes in the transmission spectra near Fermi level before and after interaction of TCDD. In the case of BNG systems, considerable variations in the transmission spectra can be seen at −0.75 eV below the Fermi level which arises due to the interaction of TCDD. The zero bias DOS of both the hetero bilayers are presented in ESI Fig. S4(b) and S5(b).†
![]() | ||
Fig. 8 Zero bias transmission spectra (a) and DOS (b) of SLG nanoribbon before and after interaction with TCDD. |
![]() | ||
Fig. 9 Zero bias transmission spectra (a) and DOS (b) of BLG nanoribbon before and after interaction with TCDD. |
![]() | ||
Fig. 10 Zero bias transmission spectra (a) and DOS (b) of SLBN nanoribbon before and after interaction with TCDD. |
![]() | ||
Fig. 11 Zero bias transmission spectra (a) and DOS (b) of BLBN nanoribbon before and after interaction with TCDD. |
The current as a function of applied voltage of SLG–TCDD nanoribbon complex is shown in Fig. 12(a). The current flow through SLG–TCDD nanoribbon is slightly less than that of SLG. The interaction of TCDD decreases the conductivity of SLG nanoribbon which is in good agreement with the calculated transmission spectra and DOS. The same trend is observed after interaction of TCDD with BLG as presented in Fig. 12(b). Interestingly, the conductivity of SLBN nanoribbon after interaction with TCDD increases the flow of current when the input bias voltage is 3.0 V. The Negative Differential Resistance (NDR) effect has been observed in the SLBN after interaction with TCDD as shown in Fig. 13(a). The NDR effect is an important quantum transport phenomenon. It is very useful in many device applications like analog to digital conversions, amplifications, noise rectifications, etc. The NDR effect is a nonlinear changes in the current with respect to an applied bias voltage. The origin of the observed NDR effect may be related to the changes in the transmission spectra and the DOS with different input bias voltage.
The DOS for the SLBN–TCDD (at Vb = 3.8 V) nanoribbon is shown in Fig. 14(a). The peaks are observed in the bias window due to the interaction of TCDD with SLBN nanoribbon which facilitates the tunneling of electron that leads to sharp increase in the current. When Vb is equal to 4.0 V, the reduction in the degree of tunneling occurs due to the absence of such peaks in the bias window. Thus, the decrease in the current with increase in Vb is observed (Fig. 14(b)). This NDR effect has already been observed in fluorinated BNNR, CNT, GNR and porous graphene.85–87
![]() | ||
Fig. 14 The DOS of SLBN before and after interaction with TCDD at input bias voltage of (a) 3.8 V and (b) 4.0 V. The dashed line represents the bias window. |
The I–V characteristics of BLBN nanoribbon is shown in Fig. 13(b). In the bias voltage range from 0 to 2.8 V, the output current is zero. When Vb reaches to 2.8 V, an increase in the flow of current can be observed for the BLBN–TCDD nanoribbon systems. Further increasing the Vb, the flow of current also increases. The interaction of TCDD increases the flow of current. Hence, the interaction of TCDD with BLBN nanoribbon enhances the tunneling of electrons and as a result, significant increase in the current is observed. Thus, these variations are useful to detect or monitor the presence of TCDD. The I–V characteristics of hetero bilayers of GBN and BNG before and after interaction with TCDD are depicted in ESI Fig. S6(a) and (b).† It is evident from these figures that there are no noticeable changes in the I–V characteristics of these systems before and after interaction of TCDD.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6ra18976h |
This journal is © The Royal Society of Chemistry 2016 |