Cyclosarin nerve agent interaction with the pristine, Stone Wales defected, and Si-doped BN nanosheets: Theoretical study

Cyclosarin nerve agent interaction with the pristine, Stone Wales defected, and Si-doped BN nanosheets: Theoretical study

Physica E 90 (2017) 143–148 Contents lists available at ScienceDirect Physica E journal homepage: www.elsevier.com/locate/physe Cyclosarin nerve ag...

882KB Sizes 0 Downloads 48 Views

Physica E 90 (2017) 143–148

Contents lists available at ScienceDirect

Physica E journal homepage: www.elsevier.com/locate/physe

Cyclosarin nerve agent interaction with the pristine, Stone Wales defected, and Si-doped BN nanosheets: Theoretical study ⁎

MARK



K. Nejatia, S. Arshadia, E. Vessallya, , A. Bekhradniab, , A. Hosseinianc a b c

Department of Chemistry, Payame Noor University, Tehran, Iran Pharmaceutical Sciences Research Center, Department of Medicinal Chemistry, Mazandaran University of Medical Sciences, Sari, Iran Department of Engineering Science, College of Engineering, University of Tehran, P.O. Box 11365-4563, Tehran, Iran

A R T I C L E I N F O

A BS T RAC T

Keywords: BN hexagonal sheet Cyclosarin White graphene DFT Sensor

Never agent identification and disposal is vital for both civilian and military defense resources. Herein, using density functional theory calculations, the reactivity and electronic sensitivity of pristine, Stone Wales (SW) defected, and Si-doped BN (Si-BN) nanosheets toward cyclosarin nerve agent were investigated. It was found that the interaction of cyclosarin with the pristine BN sheet is very weak and also that is not energetically favorable with SW defected one. Unlike the SW defect, replacing a B atom by Si atom significantly makes the cyclosarin adsorption energetically favorable. Calculations show that the carbonyl and etheric oxygen atoms of cyclosarin attack the Si atom of Si-BN with the adsorption energies of −73.5 and −136.9 kJ/mol, respectively. The cyclosarin nerve agent can be decomposed by the Si-BN sheet which is thermodynamically highly favorable. Upon this process, the HOMO and LUMO levels are significantly unstabilized and the HOMO-LUMO gap significantly changed by about 24.2%. The cyclosarin presence and its decomposition by Si-BN sheet can be recognized because of the electrical conductivity change of the sheet.

1. Introduction Decomposition of nerve agents is an important chemical process for hazardous chemical wastes disposal [1]. Also, their identification is vital for both civilian and military defense resources [2]. To date, many methods have been introduced for nerve agent detection and decomposition [3–5]. By advent of nanotechnology, nanostructures have attracted an extensive attention as gas sensors, and surfaces for gas adsorption and decomposition because of their surface/volume ratio [6–17]. Single-walled carbon nanotubes have been used to recognize nerve agents with reproducible resistance change [18]. However, practically, working with carbon nanotubes is a hard task due to their electronic dependency on the chirality [19]. Recently, a great attention has been devoted to the inorganic BN nanostructures whose electronic properties are less depended on the chirality which makes them more suitable for electronic devices [20–23]. An important isostructural arrangement of graphene is hexagonal BN nanosheet ‘white graphene’ which has been synthesized using various methods [24]. Many researchers investigated the gas adsorption and sensing properties of the withe graphene [25–29]. Also, doping, creating defects, and functionalization approaches have been investigated as strategies to tailor the reactivity and sensitivity of



nanomaterials [30–39]. It has been indicated that BN nanosheet can be employed in the pristine form for detection of NO2 gas and in the Aldoped form for sensing the nitrophenol [30]. Computational methods significantly help the experimentalist to understand the adsorption and sensing behavior of different surfaces to chemicals. In this letter, we explore the adsorption behavior and sensitivity the pristine, StoneWales detected (SW-BN) and Si-doped BN (Si-BN) nanosheets to cyclosarin nerve agent by means of density functional theory (DFT) calculations. Experimentally, several works have been published on the synthesis of different Si-BN nanomaterials [40–42]. Fan et al. have reported the growth of silicon-doped BN nanotubes via catalyst-assisted pyrolysis of a boron-containing polymeric precursor, characterizing their morphologies and structures using electron microscopy and Raman spectroscopy [40]. Also, Si-doped BN films with various Si concentrations were obtained by in situ cosputtering during ion beam assisted deposition by Ying et al. [41]. Si-doped multiwalled BN nanotubes were synthesized via thermal chemical vapor deposition by Cho et al. [42]. Electron energy-loss spectroscopy showed that 5% of Si atoms were homogeneously doped into the nanotubes. X-ray absorption and photoelectron spectroscopy measurements indicated that the Si−B and Si−N bonding structures are produced.

Corresponding authors. E-mail addresses: [email protected] (E. Vessally), [email protected] (A. Bekhradnia).

http://dx.doi.org/10.1016/j.physe.2017.03.023 Received 30 November 2016; Received in revised form 13 February 2017; Accepted 27 March 2017 Available online 27 March 2017 1386-9477/ © 2017 Elsevier B.V. All rights reserved.

Physica E 90 (2017) 143–148

K. Nejati et al.

(~5.50 eV) [53]. To obtain the stable cyclosarin/BN complexes, several initial structures entailing single atom (F, O, and H), two, three, and more atoms near the surface of the BN were explored. It was predicted that the molecule tends to escape from the surface of BN sheet upon the relax optimization and a weak interaction is predicted. Finally, in the most stable form, the cyclosarin located on the surface of the BN as shown in Fig. 1 with adsorption energy of about −1.32 kJ/mol. Because of this weak interaction, the electronic properties of the BN sheet are not significantly affected by the adsorption process.

2. Computational methods The all calculations were performed using the B3LYP functional augmented with an empirical dispersion term (B3LYP-D) and the 6– 31G* basis set as implemented in the GAMESS code [43]. B3LYP delivers an efficient and robust basis for III–V semiconductor calculation and has been frequently used for nanomaterials [44–51]. We used a model of BN nanosheet which is consisted from 36 B and 36 N atoms that its dangling bonds is saturated with hydrogen atoms to reduce the boundary effects. Previously, this model has been frequently used for different purposes and it has been indicated that it is a good representative for larger BN nanosheets [25–27,29,30]. GaussSum code was hired to attain DOS results [52]. The adsorption energy (Ead) is predicted as follows:

3.2. The adsorption of cyclosarin on the SW-BN sheet The optimized structure of SW-BN sheet is shown in Fig. 2. As it was schematically indicated (Fig. 2), for creating SW defect, four adjacent hexagonal rings at the center of BN sheet are changed into two pentagonal and two heptagonal rings when the bond uniting two of the adjacent rings rotates. Length of the rotated bond is change from 1.44 to 1.39 Å due to a structural distortion. Upon the SW-defect formation an N-N and a B-B bond are formed with length of about 1.48 and 1.67 Å, respectively. The energy is needed for SW defect formation is calculated to be about −133.6 kcal/mol. In the other words, the pristine BN sheet is more stable than the SW-BN by about 133.6 kcal/ mol. The HOMO is unstabilized by about 0.48 eV by shifting from −6.30 eV in the pristine BN sheet to −5.82 eV in the SW-BN sheet. While the LUMO is stabilized by changing from −0.42 to −0.81 eV. Also, their shapes are significantly changed by localizing on the newly formed B-B (LUMO) and N-N (HOMO) bonds as shown in Fig. 2. Finally, the Eg is somewhat decreased from 5.88 to 5.01 eV, indicating that existence of SW defects decreases the electrical resistance of the BN sheet. As it was mentioned before the calculated Eg of the pristine BN is larger than that of the experimental value which may be due to the naturally existence of defects in the real BN sheet. It is expected that the high electron heads (O-head and F-head) of cyclosarin molecule attack the electron deficient B-B site of the SW defect. Our calculations indicate that when the molecule is located from its F head on the B-B site, it reoriented and tends to attach from its carbonyl O-head to a B atom of B-B bond as shown in Fig. 3. The length of newly formed O-B is about 1.75 Å and a local deformation is occurred in the adsorption site. The adsorbing B atom is projected out of plan and the lengths of its surrounding bonds are perturbed. These findings show that the interaction between the cyclosarin molecule and SW-BN sheet is much stronger than that with the pristine BN sheet. But the calculated adsorption energy is positive for this interaction (Table 1) because of the large structural deformation that compensates the released energy. To show the strength of the adsorption we used the following equation to find the binding energy:

Ead = E (cyclosarin/adsorbent) − E (adsorbent) − E (cyclosarin) + E (BSSE) (1) where E(adsorbent) is the total energy of a pristine, SW-BN or Si-BN sheet. E (cyclosarin/ adsorbent) is the total energy of the complex of the adsorbed cyclosarin molecule on the adsorbent. E (BSSE) is the basis set superposition error (BSSE) corrected for all adsorption energies. To predict the electronic sensitivity Eq. (2) was used which relates the HOMO-LUMO gap (Eg) to the electrical conductivity of a semiconductor as follows:

σ = A T3/2 exp (−Eg/2kT)

(2) 3

3/2

where k is the Boltzmann's constant and A (electrons/m K constant.

) is a

3. Results and discussion 3.1. The cyclosarin on the pristine BN sheet We have shown the structure of BN nanosheet in Fig. 1, revealing that the B-N bond lengths are about 1.44–1.47 Å in agreement with previous reports [25]. The HOMO and LUMO of the sheet are lied at −6.30 and −0.42 eV, respectively, representing an Eg of about 5.88 eV which is approximately in accordance with the experimental bad gap

Ebin = E (cyclosarin/SW−BN) − Esp (SW−BN) − Esp (cyclosarin) + E (BSSE)

(3)

where E (cyclosarin/SW-BN) is total energy of cyclosarin/SW-BN complex, and Esp (SW-BN) is single point energy of the SW-BN after removing the cyclosarin from the complex and also Esp (cyclosarin) is the single point energy of cyclosarin after removing the SW-BN from the complex. Finally the predicted binding energy is predicted to be about −63.2 kJ/mol, indicating a strong interaction. The HOMO and LUMO levels are unstabilized upon the cyclosarin adsorption process. But Table 1 shows that the Eg of the SW-BN is not sensibly affected. Thus, the electrical conductivity of the SW-BN will not vary meaningfully based on the Eq. (2). It can be decided that although the SW defect considerably increases the reactivity of the BN sheet, its effect on the electrical sensitivity is negligible. 3.3. The cyclosarin adsorption on the Si-BN sheet Fig. 1. Optimized structure of BN nanosheet and its complex with cyclosarin nerve agent.

Herein, we will explore the effect of Si-doping on the reaction of the 144

Physica E 90 (2017) 143–148

K. Nejati et al.

Fig. 2. (a) A schematic view for Stone Wales (SW) defect creation, (b) the SW defected BN nanosheet, and its (c) LUMO and, (d) HOMO levels. Distances in Å. Table 1 The adsorption energy (Ead, kJ/mol) for adsorption of cyclosarin on the pristine, Stone Wales (SW) defected, and Si-doped BN nanosheets. The energies of HOMO, LUMO, and HOMO-LUMO gap (Eg) are in eV. The %ΔEg indicates the change of Eg after the adsorption process. Structure

Ead

EHOMO

ELUMO

Eg

%ΔEg

BN BN/Sarin SW-BN SW-BN/Sarin Si-BN Complex A Complex B

– −1.32 – 37.6 – −73.5 −136.6

−6.30 −6.34 −5.82 −5.27 −4.84 −4.91 −4.11

−0.42 −0.39 −0.81 −0.25 −1.95 −1.72 −0.52

5.88 5.96 5.01 5.02 2.89 3.19 3.59

– 1.3 – 0.1 – 10.3 24.2

atom (~1.17 Å) than the N one (~0.75 Å) [56]. The optimized structure of the BN in which a B atom is replaced by Si atom (Si–BN sheet) is shown in Fig. 4. The geometry of BN sheet is noticeably distorted around the Si-doped site and the dopant Si atom projects out of the sheet surface to decrease the compressive stress. The length of newly formed Si-N bond is approximately 1.74 Å, and the Si NBO (natural bond orbitals) hybridization is almost sp3. The HOMO level of Si-BN sheet is singly occupied molecular orbital which significantly tends to accept electron to be fulfilled. As shown in Fig. 4, this level is largely localized on the dopant Si atom which makes it suitable for a nucleophilic attack. The energy of HOMO level of Si-BN sheet is higher than that of the pristine BN sheet by about 1.46 eV (Table 1) which makes its more accessible for the nucleophilic attack. The DOS plot (Fig. 4) indicates that two new states are appeared within the Eg of the BN sheet at −4.84 and −1.95 eV after the Si-doping process. The LUMO level of BN sheet is shifted from −0.42 to −1.95 eV and the Eg is narrowed by about 2.99 eV. It is well kwon that the smaller Eg causes lower kinetic stability and thereby more reactivity. Our calculations indicate that the cyclosarin molecule energetically

Fig. 3. The complex between Stone Wales defected BN sheet and cyclosarin molecule. Distance in Å.

cyclosarin with the BN sheet and also its effect on the electronic properties. Our desired purpose is increasing the sensitivity and reactivity of this sheet toward the nerve agent. Several Si-doped nanostructures have been previously inspected theoretically and also their synthesis has been reported [54,55]. It can be seen that a B or N atom can be replaced by Si atom. It has been previously shown [56] that replacing N atom by Si atom is energetically unfavorable and the structure of doped sheet is a saddle point with large imaginary frequencies. The imaginary frequency orientation has indicated that the dopant Si tends to escape from the sheet surface in this complex. This is due to very larger size of the Si atom compared to the N atom. However, the radius of Si atom (~1.46 Å) is more close to that of the B

145

Physica E 90 (2017) 143–148

K. Nejati et al.

tends to attach to the Si atom via its carbonyl or etheric oxygen atom. When cyclosarin is located from its F atom on the Si atom, it reoriented and attached from oxygen atom to the Si dopant. When it is adsorbed from carbonyl oxygen, the complex A (Fig. 5) is obtained in which the newly formed O-Si bond is about 1.65 Å and an energy of about 73.5 kJ/mol is released. It should be noted that this large energy and structural deformation indicate that the recovery of the Si-BN will be a hard task. However, the main potential application of this compound may be decomposition and removal of the nerve agent in the state of emergency and war. Although at this condition the recovery of Si-BN may not be a vital concern, this may be a disadvantage of this compound. We used the Synchronous Transit-Guided Quasi-Newton (STQN) method [57] to predict the transition state (TS) structure for the adsorption process on the Si-BN. The predicted TS structure is shown in Fig. 1S (Supplementary material), indicating that the formed Si-O bond is about 1.89 Å which is shorter about 0.22 Å form that in the final structure. Also, the O-C bond is weakened and increased to 2.32 Å which is smaller than that in the final structure. The calculated energy barrier is about 61.2 kJ/mol which is not too large to be overcome at room temperature. The HOMO and LUMO are slightly changed and the Eg is changed by about 10.3% (Table 1). When the molecule attaches from its etheric oxygen atom, the cyclohexane ring is separated and an intermediate complex (B, Fig. 5) is created which tends to transform to the complex A. The adsorption energy about −136.6 kJ/mol indicates that the interaction is energetically highly favorable. This decomposition process is useful for disposal and removing the toxicity of this nerve agent which can be recognized because of a large change in the electronic properties of the Si-BN sheet. Upon the decomposition process, the HOMO and LUMO levels are significantly unstabilized. The HOMO is shifted from −4.84 to −4.11 eV and the LUMO from −1.95 to −0.52 eV. In accordance with the energy change, the shapes of the HOMO and LUMO are meaningfully changed. As shown in Fig. 6, the HOMO level of complex B is shifted on the hexagonal ring of the cyclosarin molecule and the LUMO is mainly focused on the edge of the sheet which is ended by the electron deficient B atoms. The change of the energy of the LUMO level is much more than that of the HOMO level, and thus, the Eg is significantly changed by about 24.2%. Based on the Eq. (2), this large change of the Eg will exponentially alter the electrical conductivity of the Si-BN sheet which can be converted to an electrical signal. Therefore, by appearing this signal the existence and decomposition of the cyclosarin can be recognized.

Fig. 4. Optimized structure of Si-doped BN sheet and its DOS plot and HOMO level profile. Distance in Å.

4. Conclusions We have studied the interaction between the cyclosarin nerve agent and pristine BN, SW-BN, and Si-BN nanosheets by means of DFT calculation. It was found that the interaction of this gas with the pristine BN sheet is very weak with adsorption energy of about −1.32 kJ/mol, and the electronic properties of the sheet are not affected sensibly. Creating a SW defect, significantly improves the interaction, but it is thermodynamically is unfavorable due to a large structural deformation. When a B atom is replaced by an Si one, the reactivity and electronic sensitivity of the sheet are considerably increased toward cyclosarin agent. It was revealed that cyclosarin is decomposed by the Si-BN, releasing an energy of 136.6 kJ/mol. The decomposition process significantly changes the Eg of the sheet by about 24.2%. This large change exponentially alters the electrical conductivity of the Si-BN sheet which can be converted to an electrical signal which identifies the existence and decomposition of the cyclosarin. Acknowledgements The authors gratefully acknowledge the financial support of this work by the Mazandaran University of Medical Sciences "Professor's Projects Funds".

Fig. 5. Optimized structure of cyclosarin/Si-doped BN sheet complexes. Distances in Å.

146

Physica E 90 (2017) 143–148

K. Nejati et al.

[11] [12] [13] [14] [15]

[16] [17] [18] [19]

[20]

[21] [22]

[23] [24]

[25] [26]

[27] [28]

Fig. 6. The HOMO and LUMO profiles of the complex B (Fig. 5).

[29]

Appendix A. Supplementary material

[30]

Supplementary data associated with this article can be found in the online version at doi:10.1016/j.physe.2017.03.023.

[31]

References

[32]

[1] D. Mandal, B. Mondal, A.K. Das, Isomerization and decomposition of a model nerve agent: a computational analysis of the reaction energetics and kinetics of dimethyl ethylphosphonate, J. Phys. Chem. A 114 (2010) 10717–10725. [2] R. Dawson, Review of oximes available for treatment of nerve agent poisoning, J. Appl. Toxicol. 14 (1994) 317–331. [3] (a) L. Safari, E. Vessally, A. Bekhradnia, A. Hosseinian, L. Edjlali, A DFT study on the sensitivity of two-dimensional BN nanosheet to nerve agents cyclosarin and tabun, Thin Solid Films 623 (2017) 157–163; (b) S.A. Siadati, E. Vessally, A. Hosseinian, L. Edjlali, Possibility of sensing, adsorbing, and destructing the Tabun-2D-skeletal (Tabun nerve agent) by C20 fullerene and its boron and nitrogen doped derivatives, Synth. Met. 220 (2016) 606–611; (c) E. Vessally, S.A. Siadati, A. Hosseinian, L. Edjlali, Selective sensing of ozone and the chemically active gaseous species of the troposphere by using the C20 fullerene and graphene segment, Talanta 162 (2017) 505–510. [4] S. Kim, B. Cho, H. Sohn, Detection of nerve agent stimulants based on photoluminescent porous silicon interferometer, Nanoscale Res. Lett. 7 (2012) 1–5. [5] M.T. Baei, A.A. Peyghan, Z. Bagheri, A computational study of AlN nanotube as an oxygen detector, Chin. Chem. Lett. 23 (2012) 965–968. [6] J. Beheshtian, A.A. Peyghan, Z. Bagheri, M. Kamfiroozi, The alkali metal interactions with MgO nanotubes, Bull. Korean Chem. Soc. 33 (2012) 1925–1928. [7] R.B. dos Santos, F. de Brito Mota, R. Rivelino, A. Kakanakova-Georgieva, G.K. Gueorguiev, Van der Waals stacks of few-layer h-AlN with graphene: an ab initio study of structural, interaction and electronic properties, Nanotechnology 27 (2016) 145601–145609. [8] M.T. Baei, A.A. Peyghan, M. Moghimi, S. Hashemian, Electric field effect on the zigzag (6, 0) single-wall BC2N nanotube for use in nano-electronic circuits, J. Mol. Model. 19 (2013) 97–107. [9] A.A. Peyghan, M.T. Baei, S. Hashemian, P. Torabi, Adsorption of CO molecule on AlN nanotubes by parallel electric field, J. Mol. Model. 19 (2013) 859–870. [10] W. Zhou, J. Zhou, J. Shen, C. Ouyang, S. Shi, First-principles study of high-capacity

[33] [34]

[35] [36] [37] [38]

[39]

[40]

[41]

[42]

[43]

[44]

147

hydrogen storage on graphene with Li atoms, J. Phys. Chem. Solids 73 (2012) 245–251. M. Noei, DFT study on the sensitivity of open edge graphene toward CO2 gas, Vacuum 131 (2016) 194–200. J. Beheshtian, A. Ahmadi Peyghan, Z. Bagheri, A first-principles study of H2S adsorption and dissociation on the AlN nanotube, Physica E 44 (2012) 1963–1968. A.A. Peyghan, M.T. Baei, S. Hashemian, M. Moghimi, Adsorption of nitrous oxide on the (6,0) magnesium oxide nanotube, Chin. Chem. Lett. 23 (2012) 1275–1278. M.T. Baei, A.A. Peyghan, Z. Bagheri, Quantum chemical analysis on hydrogenated Zn12O12 nanoclusters, C.R. Chim. 16 (2012) 122–128. M. Moradi, A.A. Peyghan, Z. Bagheri, Tuning the electronic properties of C30B15N15 fullerene via encapsulation of alkali and alkali earth metals, Synth. Met. 177 (2013) 94–99. M. Pashangpour, A.A. Peyghan, Adsorption of carbon monoxide on the pristine, Band Al-doped C3N nanosheets, J. Mol. Model. 21 (2015) 116. A.A. Peyghan, S.A. Aslanzadeh, M. Noei, A density functional study on the acidity properties of pristine and modified SiC nano-sheets, Physica B 443 (2014) 54–59. C. Kyle, D.K. Rashmi, K. Mikhail, K.M. Sanjeev, Flexible carbon nanotube sensors for nerve agent simulants, Nanotechnology 17 (2006) 4123–4126. D. Sánchez-Portal, E. Artacho, J.M. Soler, A. Rubio, P. Ordejón, Ab initio structural, elastic, and vibrational properties of carbon nanotubes, Phys. Rev. B 59 (1999) 12678–12683. L. Song, L. Ci, H. Lu, P.B. Sorokin, C. Jin, J. Ni, A.G. Kvashnin, D.G. Kvashnin, J. Lou, B.I. Yakobson, Large scale growth and characterization of atomic hexagonal boron nitride layers, Nano Lett. 10 (2010) 3209–3215. A. Soltani, A. Ahmadi Peyghan, Z. Bagheri, H2O2 adsorption on the BN and SiC nanotubes: a DFT study, Physica E 48 (2013) 176–180. R. Freitas, G.K. Gueorguiev, F. de Brito Mota, C. de Castilho, S. Stafström, A. Kakanakova-Georgieva, Reactivity of adducts relevant to the deposition of hexagonal BN from first-principles calculations, Chem. Phys. Lett. 583 (2013) 119–124. M.T. Baei, Z. Bagheri, A.A. Peyghan, Transition metal atom adsorptions on a boron nitride nanocage, Struct. Chem. 24 (2013) 1039–1044. L. Song, L. Ci, H. Lu, P.B. Sorokin, C. Jin, J. Ni, A.G. Kvashnin, D.G. Kvashnin, J. Lou, B.I. Yakobson, Large scale growth and characterization of atomic hexagonal boron nitride layers, Nano Let. 10 (2010) 3209–3215. M. Samadizadeh, A.A. Peyghan, S.F. Rastegar, Sensing behavior of BN nanosheet toward nitrous oxide: a DFT study, Chin. Chem. Lett. 26 (2015) 1042–1045. J. Beheshtian, H. Soleymanabadi, A.A. Peyghan, Z. Bagheri, A DFT study on the functionalization of a BN nanosheet with PC-X,(PC=phenyl carbamate, X=OCH3, CH3, NH2, NO2 and CN), Appl. Surf. Sci. 268 (2012) 436–441. M. Eslami, V. Vahabi, A.A. Peyghan, Sensing properties of BN nanotube toward carcinogenic 4-chloroaniline: a computational study, Physica E 76 (2016) 6–11. J. Beheshtian, A.A. Peyghan, Z. Bagheri, Detection of phosgene by Sc-doped BN nanotubes: a DFT study, Sens. Actuators B: Chem. (2012) 846–852. J. Beheshtian, A.A. Peyghan, Z. Bagheri, Functionalization of BN nanosheet with N2H4 may be feasible in the presence of Stone–Wales defect, Struct. Chem. 24 (2013) 1565–1570. A.A. Peyghan, M. Noei, S. Yourdkhani, Al-doped graphene-like BN nanosheet as a sensor for para-nitrophenol: DFT study, Superlattices Microstruct. 59 (2013) 115–122. M. Moradi, A.A. Peyghan, Role of sodium decoration on the methane storage properties of BC3 nanosheet, Struct. Chem. 25 (2014) 1083–1090. A.A. Peyghan, S.P. Laeen, S.A. Aslanzadeh, M. Moradi, Hydrogen peroxide reduction in the oxygen vacancies of ZnO nanotubes, Thin Solid Films 556 (2014) 566–570. A.A. Peyghan, M. Moradi, Influence of antisite defect upon decomposition of nitrous oxide over graphene-analogue SiC, Thin Solid Films 552 (2014) 111–115. R.P. Gandhiraman, V. Gubala, C.C. O’Mahony, T. Cummins, J. Raj, A. Eltayeb, C. Doyle, B. James, S. Daniels, D.E. Williams, PECVD coatings for functionalization of point-of-care biosensor surfaces, Vacuum 86 (2012) 547–555. A.A. Peyghan, S. Yourdkhani, Exohedral functionalization of C60 by [4+2] cycloaddition of multiple anthracenes, Struct. Chem. 25 (2014) 785–791. J. Beheshtian, M.T. Baei, A.A. Peyghan, Z. Bagheri, Nitrous oxide adsorption on pristine and Si-doped AlN nanotubes, J. Mol. Model. 19 (2013) 943–949. A.A.P. Javad Beheshtian, Theoretical study on the functionalization of BC2N nanotube with amino groups, J. Mol. Model. 19 (2013) 2211–2216. D. Robati, S. Bagheriyan, M. Rajabi, O. Moradi, A.A. Peyghan, Effect of electrostatic interaction on the methylene blue and methyl orange adsorption by the pristine and functionalized carbon nanotubes, Physica E 83 (2016) 1–6. N.L. Hadipour, A. Ahmadi Peyghan, H. Soleymanabadi, Theoretical study on the Al-doped ZnO nanoclusters for CO chemical sensors, J. Phys. Chem. C 119 (2015) 6398–6404. Y. Fan, Y. Wang, J. Lou, S. Xu, L. Zhang, L. An, H. Heinrich, Formation of silicon‐ doped boron nitride bamboo structures via pyrolysis of a polymeric precursor, J. Am. Ceram. Soc. 89 (2006) 740–742. J. Ying, X. Zhang, Z. Yin, H. Tan, S. Zhang, Y. Fan, Electrical transport properties of the Si-doped cubic boron nitride thin films prepared by in situ cosputtering, J. Appl. Phys. 109 (2011) 023716. Y.J. Cho, C.H. Kim, H.S. Kim, J. Park, H.C. Choi, H.-J. Shin, G. Gao, H.S. Kang, Electronic structure of Si-doped BN nanotubes using x-ray photoelectron spectroscopy and first-principles calculation, Chem. Mater. 21 (2009) 136–143. M.W. Schmidt, K.K. Baldridge, J.A. Boatz, S.T. Elbert, M.S. Gordon, J.H. Jensen, S. Koseki, N. Matsunaga, K.A. Nguyen, S. Su, T.L. Windus, M. Dupuis, J.A. Montgomery, J. Comput. Chem. 14 (1993) 1347–1363. J. Beheshtian, A. Ahmadi Peyghan, Z. Bagheri, Hydrogen dissociation on diene-

Physica E 90 (2017) 143–148

K. Nejati et al. functionalized carbon nanotubes, J. Mol. Model. 19 (2012) 255–261. [45] E. Vessally, S. Soleimani-Amiri, A. Hosseinian, L. Edjlali, A. Bekhradnia, The Hartree-Fock exchange effect on the CO adsorption by the boron nitride nanocage, Physica E 87 (2017) 308–311. [46] J. Beheshtian, A.A. Peyghan, Z. Bagheri, M. Kamfiroozi, Interaction of small molecules (NO, H2, N2, and CH4) with BN nanocluster surface, Struct. Chem. 23 (2012) 1567–1572. [47] S. Tomic, B. Montanari, N.M. Harrison, The group III–V's semiconductor energy gaps predicted using the B3LYP hybrid functional, Physica E 40 (2008) 2125–2127. [48] J. Beheshtian, A.A. Peyghan, Z. Bagheri, Quantum chemical study of fluorinated AlN nano-cage, Appl. Surf. Sci. 259 (2012) 631–636. [49] (a) K. Nejati, A. Hosseinian, L. Edjlali, E. Vessally, The effect of structural curvature on the cell voltage of BN nanotube based Na-ion batteries, J. Mol. Liq. 229 (2017) 167–171; (b) A. Hosseinian, A. Bekhradnia, E. Vessally, L. Edjlali, M.D. Esrafili, A DFT study on the central-ring doped HBC nanographenes, J. Mol. Graph. Model. 73 (2017) 101–107; (c) A. Hosseinian, Z. Asadi, L. Edjlali, A. Bekhradnia, E. Vessally, NO2 sensing properties of a borazine doped nanographene: A DFT study, Comput. Theor. Chem. 1106 (2017) 36–42; (d) K. Nejati, A. Hosseinian, A. Bekhradnia, E. Vessally, L. Edjlali, Na-ion batteries based on the inorganic BN nanocluster anodes: DFT studies, J. Mol. Graph. Model. 74 (2017) 1–7; (e) S. Bashiri, E. Vessally, A. Bekhradnia, A. Hosseinian, L. Edjlali, Utility of

[50] [51] [52] [53]

[54] [55]

[56]

[57]

148

extrinsic [60] fullerenes as work function type sensors for amphetamine drug detection: DFT studies, Vacuum 136 (2017) 156–162; (f) E. Vessally, S. Soleimani-Amiri, A. Hosseinian, L. Edjlali, A. Bekhradnia, A comparative computational study on the BN ring doped nanographenes, Appl. Surf. Sci. 396 (2017) 740–745. J. Beheshtian, A.A. Peyghan, Z. Bagheri, Nitrate adsorption by carbon nanotubes in the vacuum and aqueous phase, Mon. Chem. 143 (2012) 1623–1626. J. Beheshtian, A.A. Peyghan, Z. Bagheri, Theoretical investigation of C60 fullerene functionalization with tetrazine, Comput. Theor. Chem. 992 (2012) 164–167. N. O’Boyle, A. Tenderholt, K. Langner, Cclib: a library for package‐independent computational chemistry algorithms, J. Comput. Chem. 29 (2008) 839–845. L. Song, L. Ci, H. Lu, P.B. Sorokin, C. Jin, J. Ni, A.G. Kvashnin, D.G. Kvashnin, J. Lou, B.I. Yakobson, Large scale growth and characterization of atomic hexagonal boron nitride layers, Nano Let. 10 (2010) 3209–3215. J. Beheshtian, A.A. Peyghan, Z. Bagheri, Detection of phosgene by Sc-doped BN nanotubes: a DFT study, Sens. Actuators B: Chem. (2012) 846–852. Y.C. Kim, C.H. Yoon, J. Park, J. Yoon, N.S. Han, J.K. Song, S.M. Park, J.S. Ha, Effect of Sb doping on the opto-electronic properties of SnO2 nanowires, Thin Solid Films 520 (2012) 6471–6475. F. Behmagham, E. Vessally, B. Massoumi, A. Hosseinian, L. Edjlali, A computational study on the SO2 adsorption by the pristine, Al, and Si doped BN nanosheets, Superlattices Microstruct. 100 (2016) 350–357. C. Peng, H. Bernhard Schlegel, Combining synchronous transit and quasi‐newton methods to find transition states, Isr. J. Chem. 33 (1993) 449–454.