Translate this page into:
Effects of influence of carbon ring-doping on NMR parameters of boron phosphide nanotubes: A DFT study
*Tel./fax: +98 851 3339843 mrsameti@gmail.com (M. Rezaei-Sameti)
-
Received: ,
Accepted: ,
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 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.
(a) 2D and (b) 3D views of model (A) and model (B) carbon ring doping in BPNTs.
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:
3 Results and discussion
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. (a) and (b) are the values of bond length of the C-doping of models A and B. (a) and (b) are the values of bond angles of C-doping of models A and B.
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
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
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 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. 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.
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
–
–
–
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 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
- Nano Lett.. 2006;6:1812.
- Science. 1999;284:1340.
- Europhys. Lett.. 1994;28:335.
- J. Am. Chem. Soc.. 2005;127:17144.
- Science. 1995;269:966.
- J. Chem. Phys.. 1971;54:724.
- J. Cryst. Growth. 2008;310:3973.
- Gaussian 98. Pittsburgh, PA: Gaussian Inc.; 2001.
- Solid State Commun.. 2006;137:549.
- Carbon Nanotubes and Related Structures. Cambridge: Cambridge University Press; 1999.
- Phys. Rev.. 1965;140:1133.
- Phys. Rev. B. 1988;37:785.
- Phys. Lett. A. 2006;357:369.
- Phys. Rev. Lett.. 1996;76:4737.
- J. Phys. Chem.. 2009;223:815.
- Physica E. 2009;41:883.
- J. Mol. Model.. 2010;17:89.
- Physica E. 2010;42:1667.
- Physica E. 2008;40:3060.
- Chem. Phys. Lett.. 2008;461:246.
- Catal. Today. 2007;120:341.
- J. Mol. Struct. (THEOCHEM). 2007;817:137.
- Physica B. 2002;323:357.
- J. Appl. Phys.. 1998;83:1660.
- Z. Naturforsch.. 2007;62:711.
- Comput. Mater. Sci.. 2009;47:220.
- Phys. Rev. B. 1995;51:11229.
- Chem. Phys. Lett.. 2003;371:426.
- J. Mol. Struct. (THEOCHEM). 2006;765:1.
- Acc. Chem. Res.. 2002;35:1045.