10.8
CiteScore
 
5.3
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Letter to Editor
Original Article
Research article
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
10.8
CiteScore
5.3
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Letter to Editor
Original Article
Research article
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
View/Download PDF

Translate this page into:

Original Article
8 (
2
); 168-173
doi:
10.1016/j.arabjc.2010.10.035

Effects of influence of carbon ring-doping on NMR parameters of boron phosphide nanotubes: A DFT study

Department of Applied Chemistry, Faculty of Science, Malayer University, Malayer 65174, Iran

*Tel./fax: +98 851 3339843 mrsameti@gmail.com (M. Rezaei-Sameti)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University.

Available online 15 February 2011

Abstract

The electronic structure of boron phosphide nanotubes (BPNTs) and influence of carbon ring doping (C-doping) in the horizontal region (model A) and vertical region (model B) of BPNTs is studied by density functional theory (DFT). At first, each form was optimized at B3LYP level of theory using 6-31G bases set. After, the computed chemical shielding (CS) tensors at the sites of 11B and 31P nuclei were converted to isotropic chemical shielding (CSI) and anisotropic chemical shielding (CSA). The calculated results reveal that the CS parameters of B and P nuclei in C-ring doped on vertical region (model B) undergo more significant changes than horizontal region (model A).

Keywords

BPNTs
Chemical shielding
DFT
Carbon ring doping
1

1 Introduction

In recent years carbon nanotubes (CNTs) have became one of significant topics in the research of nanomaterial science for their unique geometry, band structures, as well as the outstanding mechanical, electrical, thermal, optical properties, molecular device and biosensor (Shen, 2009; Mirzaei and Giahi, 2010; Li et al., 2006; Mirzaei, 2009a,b). Nanotubes consisting of combinations of atoms of groups three and five of the periodic table of elements have been viewed as proper candidates instead of the CNTs (Seif et al., 2007; Mirzaei et al., 2008a,b; Zhang and Zhang, 2003). The properties of group three nitrides have been studied more than those of group III phosphide. Whereas BNNTs, which have much wider energy, band gaps weakly depending on the diameter, chirality, and the number of the walls of the tube, are semiconductors, BNNTs are very attractive materials for application in nanoscale devices (Mpourmakis and Froudakis, 2007; Zhang et al., 2006; Zhou et al., 2002; Baughman et al., 1999; Harris, 1999). The considerable efforts dedicated to the synthesis of non-carbon nanotubes with properties independent of tubular diameter and chirality led to the conclusion that the group III nitrides, especially boron and aluminum nitrides (BN and AlN), are proper materials (Loiseau et al., 1996; Chen et al., 2005; Blase et al., 1994). The stable tubular structures of the nanotubes of BN (Guo et al., 2006; Chopra et al., 1995), were initially characterized by the calculations and then were successfully synthesized by the experiments (Loiseau et al., 1996; Sieh et al., 1995; Oku, 2002; Arenal et al., 2006; Nirmala and Kolandaivel, 2007; Mirzaei, 2009a,b). The similarity between electronic properties of boron phosphide (BP) and silicon carbide (SiC) has made BP nanotube (BPNT) as an important subject of numerous studies (Mirzaei, 2010; Schroten et al., 1998; Ferreira and Leite Alves, 2008). In this computational research, at first time, the properties of the electronic structure and chemical shielding parameters of the armchair (4, 4) BPNTs and influence of carbon ring doping (C-doping) on the vertical and horizontal region of this nanotubes (Figs. 1 and 2) have been investigated by performing density functional theory (DFT) calculations of the NMR parameters. At first, the pristine and the C-doped ring structures have been allowed to relax by all atomic geometrical optimization. Subsequently, the isotropic and anisotropic chemical shielding (CSI and CSA) parameters have been calculated for the 31P and 11B atoms present in both of the pristine and the C-ring doped structures.

(a) 2D and (b) 3D views of the pristine model of BPNTs.
Figure 1
(a) 2D and (b) 3D views of the pristine model of BPNTs.
(a) 2D and (b) 3D views of model (A) and model (B) carbon ring doping in BPNTs.
Figure 2
(a) 2D and (b) 3D views of model (A) and model (B) carbon ring doping in BPNTs.

2

2 Computational methods

In this computational work, the (4, 4) armchair BPNTs and C-ring doped model of BPNTs on vertical and horizontal region with 1 nm length (Figs. 1 and 2) have been investigated.

In both models, the mouths of nanotube are saturated by H atoms. All models are individually optimized by using density functional theory (DFT) at B3LYP (Kohn and Sham, 1965; Lee et al., 1988) level of theory using the Gaussian 98 set of programs (Frisch et al., 2001). The standard 6-31G basis set was used for all models. Calculations were performed on a Pentium-PC computer with a 3000 MHz processor. The chemical shielding (CS) tensors at the sites of 11B and 31P nuclei are calculated based on the gauge included atomic orbital (GIAO) approach (Ditchfield et al., 1971).

The calculated CS tensors in principal axes system (PAS) (δ33 > δ22 > δ11) are converted to measurable NMR parameters, chemical shielding isotropic (CSI) and chemical shielding anisotropic (CSA) by using Eqs. (1) and (2), respectively (Mirzaei, 2010; Schroten et al., 1998; Ferreira and Leite Alves, 2008). The evaluated NMR parameters at the sites of 11B and 31P nuclei in the pristine and C-ring doped models (A and B) are presented in Table 3:

(1)
CSI ( ppm ) = 1 3 ( σ 11 + σ 22 + σ 33 )
(2)
CSA ( ppm ) = δ 33 - ( δ 11 + δ 22 ) / 2

3

3 Results and discussion

3.1

3.1 The optimized geometries of the BPNTs and C-ring doped

The electronic structural properties of the (4, 4) armchair BPNTs and two models of carbon ring doping on horizontal region (model A) and vertical region (model B) (see Figs. 1 and 2) consisting of the B–P bond lengths and B–P–B/P–B–P bond angels are presented in Tables 1 and 2. The results show that the bond lengths for equivalent positions in the three forms of BPNTs are similar. It was shown that when C atoms are substituted in BPNTs, in models (A) and (B), the bond length of B–P is decreased about 1.89 up to 1.46 Å (Table 1). The electronegativity of C (eC = 2.55) is considerably larger than that of B (eB = 2.04) and P (eP = 2.19) leading a charge transfer from B to C yielding asymmetric electronic charge density distribution along CB bond, therefore, the C–P length is shorter than B–P length. This trend would mean that the influence of C-doping instead of the B and P atoms on the properties of the electronic structure of the BPNTs is significant.

Table 1 The optimized geometries of (4, 4) BPNT and C-doped.A.
B–P bonding Undoped C-doped B–P bonding Undoped C-doped
(a) (b) (a) (b)
Bond length (Å)
1–2 1.89 1.89 1.96 9–13 1.90 1.44 1.88
2–3 1.91 1.91 1.56 13–17 1.90 1.44 1.90
3–4 1.90 1.91 1.43 17–18 1.90 1.41 1.90
4–5 1.89 1.89 1.90 17–21 1.90 1.89 1.91
5–6 1.89 1.89 1.93 18–22 1.90 1.53 1.40
6–7 1.91 1.91 1.85 19–23 1.90 1.84 1.45
1–8 1.89 1.89 1.89 20–24 1.90 1.53 1.88
7–12 1.91 1.53 1.97 22–23 1.90 1.92 1.46
11–12 1.90 1.42 1.54 21–25 1.90 1.89 1.93
11–15 1.90 1.44 1.48 24–28 1.90 1.91 1.88
14–15 1.90 1.45 1.46 25–26 1.90 1.89 1.91
15–19 1.90 1.45 1.46 26–30 1.90 1.90 1.38
19–20 1.90 1.41 1.53 25–29 1.90 1.91 1.93
16–20 1.90 1.46 1.90 27–28 1.90 1.89 1.54
14–18 1.90 1.46 1.43 27–31 1.90 1.91 1.40
10–14 1.90 1.46 1.46 28–32 1.90 1.90 1.94
9–10 1.90 1.42 1.90 30–31 1.90 1.89 1.40

(a) and (b) are the values of bond length of the C-doping of models A and B.

See Figs. 1 and 2 for details.
Table 2 The optimized geometries of (4, 4) BPNT and C-doped.A
P–B–P/B–P–P angle Undoped C-doped P–B–P/B–P–P angle Undoped C-doped
(a) (b) (a) (b)
Bond angles
2–1–8 122 116 98 16–20–19 111 115 111
1–2–3 108 104 122 17–21–25 117 113 108
2–3–4 122 122 105 18–17–21 122 128 102
2–3–10 115 121 124 19–20–24 111 123 123
3–4–11 112 99 122 20–24–28 121 125 101
4–5–6 122 116 104 21–25–26 121 116 108
5–6–7 108 104 111 22–26–25 111 106 112
6–7–12 115 122 109 22–26–30 115 115 118
7–12–11 113 123 123 22–23–27 113 112 117
9–10–14 122 115 125 23–22–26 121 118 118
9–13–17 111 114 105 24–28–27 121 106 118
10–14–15 122 121 117 25–26–30 112 111 112
11–15–14 111 122 117 26–25–29 124 128 103
12–11–15 122 119 121 26–30–31 123 123 118
13–17–18 122 119 100 27–28–32 112 111 124
14–18–17 111 115 121 30–31–27 108 108 122
15–19–20 122 119 121 31–27–28 115 129 117
15–14–18 122 122 118

(a) and (b) are the values of bond angles of C-doping of models A and B.

See Figs. 1 and 2 for details.

In contrast with the same value of bond lengths for different B–P bonds, the bond angles significantly deviate for different bonds. The results of Table 2 show that the value of the bond angle (3–4–11) placed at sides of carbon doped in model (A) is decreased largely from 112° to 99°. Similar results are also observed in C-doped on model (B) for the bond angle (2–1–18) from 122° to 98° and other bond angle such as ([2–3–4],[4–5–6],[13–17–18],[18–17–21],[21–25–26] and[26–25–29]). The bond angle (31–27–28) in model (A) increased from 115° to 129° and this similar trend is shown in model (B) for bond angels (27–28–32) and (30–31–27) from (112° to 124°) and (108° to 122°), respectively. Comparison of the results of models (A) and (B) reveal that the deviation in bond angle for vertical ring carbon doped (model B) is larger than horizontal ring carbon doped (model A) in BPNTs. This means that the properties of electronic structures by C-doping in vertical region undergo more changes than C-doping in the horizontal region.

3.2

3.2 The NMR parameters of the BPNTs

In this work, we investigated the influences of carbon ring doping in vertical and horizontal region of the BPNTs on the electronic structure properties, the chemical shielding (CS) tensors at the sites of various 11B and 31P atoms are calculated. In order to directly relate the calculated tensors to the experimentally measurable parameters, the tensors are converted to the isotropic CS (CSI) and the anisotropic CS (CSA) parameters. The CSI is the average value of the eigenvalues of the CS tensors, Eq. (1), and the orientation of the eigenvalues of the CS tensors into the z-axis plays a dominant role in determining the value of the CSA parameter, Eq. (2). Table 3 shows the evaluated CSI and CSA values for 11B and 31P in the pristine model of (4, 4) armchair BPNTs and two models of C-ring doping (models A and B, Figs. 1 and 2). It shows that the 32 B and 32 P atoms in pristine nanotube of BPNTs is divided into eight layers with equivalent 31P and 11B, NMR parameters for each layer, The results show that the CSI values of layers (1, 8), (2, 7), (3, 6) and (4, 5) for B sites are (40, 38, 42 and 44 ppm), respectively, and for P sites are (433, 381, 380 and 384 ppm), respectively; which means that the CSI values at the sites of the nuclei of each layer feel equivalent chemical environment.

Table 3 The NMR parameters of 11B and 31P sites.
B–11 atom CSI (ppm) CSA (ppm) P–31 atom CSI (ppm) CSA (ppm)
Undoped C-doped Undoped C-doped Undoped C-doped Undoped C-doped
(a) (b) (a) (b) (a) (b) (a) (b)
B1.1 40 40 37 88 133 94 P1.1 433 462 446 147 130 68
B1.2 40 40 43 88 133 64 P1.2 433 462 395 147 130 100
B1.3 40 37 88 94 P1.3 433 446 147 68
B1.4 40 43 88 64 P1.4 433 395 147 100
B2.1 38 58 23 83 43 100 P2.1 381 385 152 164
B2.2 38 58 83 43 P2.2 381 385 385 152 164 87
B3.1 42 70 P3.1 380 364 296 176
B3.2 42 45 70 70 P3.2 380 296
B4.1 44 16 88 128 P4.1 383 163
B4.2 44 88 P4.2 383 349 163 75
B5.1 44 79 P5.1 384 391 273 218
B5.2 44 37 79 105 P5.2 384 273
B6.1 42 52 14 87 27 117 P6.1 380 395 151 171
B6.2 42 52 87 27 P6.2 380 395 368 151 171 67
B7.1 38 52 95 47 P7.1 382 414 343 276 148 172
B7.2 38 52 33 95 47 82 P7.2 382 414 276 148
B8.1 40 52 24 109 109 48 P8.1 433 442 117 117 54
B8.2 40 39 109 109 P8.2 433 442 378

In each row, (a) and (b) are the CS values of the C-doping of models A and B. See Figs. 1 and 2 for details.

The CSI values for 11B sites in the C-ring doping of BPNTs in two models (A) and (B) have different trend. So that in model (A) CSI values at the sites of B2.1, B2.2, B6.1, B6.2 and B7.1, B7.2 increased, and in model (B), CSI values at sites B1.2, B1.4 and B3.2 increased and at the other sites of two models decreased.

The CSI values at the site of 31P nuclei in model (A) at all sites increased with respect to the pristine model; the presence of lone pair of electrons in valence shell of 31P caused the significant change more than that at the 11B nuclei. On the other hand, the CSI values for 31P nuclei in model (B) at the sites P1.1, P1.3, P2.2, increased and for other sites decreased with respect to the pristine model. An important point to note is that the CSI values for 31P nuclei in models (A) and (B) at the sites of P7.1 and P1.2, respectively, undergo high significant changes.

Since the natures of the CSI and the CSA values are different, their behaviors due to the C-ring doping in the studied BPNTs are also different. The comparison CSA and CSI values show that, CSA values for B nuclei in model (A) at layers B2, B6, B7 decreased with respect to the pristine model. Also the CSA tensors for P nuclei in model (A) at layers P2, P6 increased and at layers P1, P7 decreased. Due to high number of carbons doping in structures of BPNTs in model (B), the CSA values for B nuclei at sites B1.1, B1.3, B2.1, B4.1, B5.2 and B6.1 increased with respect to pristine model, and at the other sites decreased.

The CSA values in model (B) for P atom for all layers decreased with respect to the pristine model.

The CSA values in model (A) and model (B) at the site of B7.1 also undergo notable changes due to the chemical bonding to the B atom and C atom, respectively; therefore, the value of CSA for this nucleus significantly reduces in comparison to that of the pristine model.

The comparison of CSI and CSA values between C-ring doped in model (A) and model (B) reveal that the changes of electrical structures in model (B) are more than model (A), because of more number of carbons doped in the vertical region. It is worth noting that just the electrostatic properties of those nuclei are mainly dependent on electronic density at the sites of nuclei, therefore, due to influence of the C-doping the electronic densities of all nuclei and the CS parameters undergo changes, so the CSI values at all sites that neighbor of C doped undergo notable decrease.

4

4 Conclusions

We performed density functional theory (DFT) to calculate the CS parameters to investigate the properties and electronic structure of C-ring doping on vertical and horizontal region of the armchair models of (4, 4) BPNTs. The optimized bond lengths reveal that when C atoms are substituted in BPNTs, in models (A) and (B), the bond lengths of B–P is decreased from about 1.89 to 1.46 Å, and the bond angel (3–4–11) placed at sides of carbon doped in model (A) and the bond angle (2–1–18) in model (B) are significantly decreased. The calculated CSI and CSA values for two models reveal that at the site of P7.1 nucleus the electrical structures undergo significant changes because of contributing to the chemical bonding with the C-doping atoms. The results show that the effect of C-ring doping on CS parameters in vertical region model (B) is more sensible than horizontal region model (A).

References

  1. , , , , , , , , . Nano Lett.. 2006;6:1812.
  2. , , , , , , , , , , , , , . Science. 1999;284:1340.
  3. , , , , . Europhys. Lett.. 1994;28:335.
  4. , , , , , . J. Am. Chem. Soc.. 2005;127:17144.
  5. , , , , , , , . Science. 1995;269:966.
  6. , , , . J. Chem. Phys.. 1971;54:724.
  7. , , . J. Cryst. Growth. 2008;310:3973.
  8. , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , . Gaussian 98. Pittsburgh, PA: Gaussian Inc.; .
  9. , , , , . Solid State Commun.. 2006;137:549.
  10. , . Carbon Nanotubes and Related Structures. Cambridge: Cambridge University Press; .
  11. , , . Phys. Rev.. 1965;140:1133.
  12. , , , . Phys. Rev. B. 1988;37:785.
  13. , , , , , , , . Phys. Lett. A. 2006;357:369.
  14. , , , , , . Phys. Rev. Lett.. 1996;76:4737.
  15. , . J. Phys. Chem.. 2009;223:815.
  16. , . Physica E. 2009;41:883.
  17. , . J. Mol. Model.. 2010;17:89.
  18. , , . Physica E. 2010;42:1667.
  19. , , , , . Physica E. 2008;40:3060.
  20. , , , . Chem. Phys. Lett.. 2008;461:246.
  21. , , . Catal. Today. 2007;120:341.
  22. , , . J. Mol. Struct. (THEOCHEM). 2007;817:137.
  23. , . Physica B. 2002;323:357.
  24. , , , . J. Appl. Phys.. 1998;83:1660.
  25. , , , , . Z. Naturforsch.. 2007;62:711.
  26. , . Comput. Mater. Sci.. 2009;47:220.
  27. , , , , , , , , , , . Phys. Rev. B. 1995;51:11229.
  28. , , . Chem. Phys. Lett.. 2003;371:426.
  29. , , , , , , , , , , , . J. Mol. Struct. (THEOCHEM). 2006;765:1.
  30. , , , , , , . Acc. Chem. Res.. 2002;35:1045.
Show Sections