• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    單壁碳納米管改變手性外延生長的密度泛函理論研究

    2016-11-22 07:31:24尹利長
    新型炭材料 2016年5期
    關(guān)鍵詞:生長

    諶 為, 李 峰, 劉 暢, 尹利長

    (中國科學(xué)院金屬研究所 沈陽材料科學(xué)國家(聯(lián)合)實(shí)驗(yàn)室,遼寧 沈陽110016)

    ?

    單壁碳納米管改變手性外延生長的密度泛函理論研究

    諶 為, 李 峰, 劉 暢, 尹利長

    (中國科學(xué)院金屬研究所 沈陽材料科學(xué)國家(聯(lián)合)實(shí)驗(yàn)室,遼寧 沈陽110016)

    采用密度泛函理論計(jì)算系統(tǒng)研究了單壁碳納米管(Single-walled carbon nanotube,SWCNT)改變手性外延生長(手性指數(shù)從(n,m)變化到(n±Δ,m?Δ),其中Δ=1和2)的熱力學(xué)過程。結(jié)果表明,碳管手性變化后外延生長在熱力學(xué)上都需要吸收能量,其所需吸收的能量隨著管徑的減小線性減小。在Δ=1的情況下,由于近扶手椅型碳管改變手性時(shí),所引入的5~7元環(huán)對與管軸的夾角比近鋸齒型碳管更大,導(dǎo)致5~7元環(huán)對的形成能增加,使得管徑相同的近扶手椅型碳管比近鋸齒型碳管在改變手性生長時(shí)需要吸收更多的能量。在Δ=2的情況下,發(fā)現(xiàn)只有當(dāng)兩個(gè)必須引入的5~7元環(huán)對相互毗鄰,手性改變的外延生長所需能量最小,預(yù)測其為實(shí)驗(yàn)上最易于實(shí)現(xiàn)的碳管手性指數(shù)由(n,m)變化到(n±Δ,m?Δ)的外延生長模式。這些理論研究結(jié)果有助于深入理解SWCNTs手性變化后外延生長的熱力學(xué)行為,可為基于外延生長可控制備單一手性SWCNTs提供理論依據(jù)。

    單壁碳納米管; 手性變化; 5~7缺陷; 外延生長; 密度泛函理論

    1 Introduction

    A single-walled carbon nanotube (SWCNT) can be conceptually considered as a cylinder rolled up from a graphene sheet. The geometry structure of a SWCNT is exclusively determined by a pair of chirality indices (n,m) according to the rolling direction of the graphene sheet[1]. The electronic properties of SWCNTs are closely related to their geometry structures. For example, SWCNTs can be either metallic (m-SWCNTs) or semiconducting (s-SWCNTs) depending on their chiralities. To be specific, when (2n+m)/3 is an integer, the SWCNTs are metallic, and otherwise semiconducting[2]. The unique electronic properties of SWCNTs show appealling potential applications in many fields. For instance, s-SWCNTs can be used for high mobility nano transistors[3]and computers[4], while m-SWCNTs are suitable for fabricating transparent electrically conductive films[5], field-emission electron sources[6], nanocables[7], and so on. Unfortunately, as-grown SWCNT samples are always a mixture of m-SWCNTs and s-SWCNTs[8]due to the negligible formation energy difference among different SWCNTs with comparable diameters[9], highly limiting their widespread applications. Basically, the above issue can be possibly solved by a direct synthesis of uniformly metallic (or semiconducting) or even single-chirality nanotubes through controlling the composition, morphology and/or size of catalysts for SWCNT growth by chemical vapor deposition[10]. For example, dominantly semiconducting SWCNTs with a narrow distribution of diameter and chirality were synthesized by using bimetallic CoMn[11], FeRu[12]catalysts, or monometallic Co nanoparticles with a well-defined crystal structure[13]. Moreover, SWCNTs with only a few specific chiralities were realized by a fine regulation and design of catalysts[14-20]. Also, the chirality distribution of as-grown SWCNTs can be altered by varying the composition of NixFe1-xnanocatalysts based on an epitaxial growth model[21]. Very recently, Yang et al. reported a direct growth of the (12,6) tube with an abundance higher than 92% using WCo bimetallic nanocatalysts[22], and Fasel et al. have achieved single-chirality SWCNTs with only one chirality of (6, 6) by using surface-catalyzed cyclodehydrogenation of C96H54precusors on Pt(111) surface[23]. However, despite enormous efforts and remarkable progress, real achievements of uniformly conducting and even chirally pure SWCNTs in a large amount by a direct growth is still highly impractical due to the insufficient understanding of the growth mechanism and the multiple factors (catalyst structure and morphology[24], temperature[25], and so on) relating to the growth of SWCNTs.

    On the other hand, we also noted that the chirality of SWCNTs could be uniformly controlled by an epitaxial growth with open-end short SWCNT seeds, namely “cloning growth”. For example, Yao et al. grew SWCNTs with a controlled chirality by using open-end SWCNTs as seeds[26], Liu et al. synthesized SWCNTs with predefined chiralities using purified single-chirality nanotubes as seeds[27]. It is also found that the diameter or even the chirality of SWCNTs can be finely controlled via a tube junction formation by changing the temperature during SWCNT growth[28], which is a typical case for a chirality-changed. In principle, such a chirality-unchanged (cloning growth) or chirality-changed epitaxial growth would achieve a chirally pure SWCNT growth in a large amount with finely controlled growth conditions. However, the underlying mechanism of this cloning or chirality-changed epitaxial growth of SWCNTs is still unclear, thus highly hindering the experimental progress in this aspect.

    In this work, we systematically studied the thermodynamics for the cloning and chirality-changed growth of SWCNTs by DFT calculations, aiming to obtain guideful hints to achieve single-chirality SWCNTs by an epitaxial growth. Taking the SWCNTs with diameters ranging from 0.4 to 1.3 nm as examples, we have calculated the energy changes (ΔE) for SWCNTs with chirality change from (n,m) to (n±Δ, m?Δ),Δ=1 and 2. The energy changes for the chirality-changed epitaxial growth of SWCNTs increase linearly with increasing the tube diameters, ranging from 1.1 (2.0) to 3.7 (4.4) eV in the case ofΔ=1 (2). It is also found that the energy changes in the case ofΔ=1 are about 0.2- 0.4 eV larger for near armchair (nAC) tubes than those for near zigzag (nZZ) tubes, resulting from the larger angle difference between pentagon-heptagon defects (5,7-defects) and the tube axis for nAC tubes. In the case ofΔ=2, the topological connection with two adjacent 5,7-defects, which are indispensable for the chirality-changed growth of SWCNTs in this case, is found to be the most energetically stable. Similar to the case ofΔ=1, the energy needed to change chirality of a SWCNT forΔ=2 case increases with increasing the tube diameter. Our DFT results provide a deep understanding on the thermodynamics of the chirality-changed epitaxial growth and may guide for the chirality-controlled synthesis of SWCNTs.

    2 Computational methods

    In this work, all calculations were carried out within the framework of DFT implemented in the Vienna Ab Initio Simulation Package (VASP)[29], using the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof[30]. The electron-ion interactions were described by using the frozen-core projector augmented wave approach[31, 32]. Testing calculations show that the total energy of a (5, 1)/(4, 2) tube junction obtained by using a cut-off energy of 450 eV is only 30 meV lower than that by using the cut-off energy of 280 eV, which is negligible in comparison with the energy input (1.1 eV) for changing a (5, 1) tube to a (4, 2) tube. Thus, the energy cut-off for plane waves was set to be 280 eV for all calculations in this work. Considering that different tube junctions were constructed based on the cluster model to simulate the chirality-changed epitaxial growth for different SWCNTs, only Gamma point was used to sample the first Brillouin zone for all calculations in this work. Optimized geometry structures for all tube junctions considered in the present work were determined until the residual forces acted on each atom were less than 0.01 eV/?.

    Different tube junctions consisting of two short segments of different SWCNTs with chirality of (n,m) and (n±Δ, m?Δ) (Δ=1, 2) were constructed to study the atomic structure change and corresponding energy change induced by the chirality-changed epitaxial growth of different SWCNTs. Eight nZZ (n, 1) (n=5-11 and 15) and eight other tubes including (5,2), (5,3), (6,3), (6,4), (7,4), (7,5), (8,5) and (10,7) were considered for the case of chirality change withΔ=1, and twelve (n,m) tubes with m=2, 3, 4 and n ranging from m+4 to m+7 for case of chirality change withΔ=2, in order to achieve a convincing comparison of the energy change among different tubes in this work. It should be noted that, one of the initial purposes of this work is to understand the difference of energy changes between highly symmetrical tubes, like ZZ and AC tubes, during their chirality-changed epitaxial growths. Unfortunately, both ZZ and AC tubes will change their helicities from the left-handed to the right-handed (or from right-handed to left-handed) when their chiralities are changed from (n,m) to (n±Δ, m?Δ) (Δ=1, 2). Alternatively, we choose these nZZ and nAC tubes in this work without considering the possible influence of the tube helicity on the energy changes during the chirality-changed epitaxial growth. We use the cluster models for all the calculations in this work since the periodic models for a tube junction are too large for standard DFT calculations, and a large supercell containing a short tube (or a tube junction) was constructed to model the chirality-changed epitaxial growth of SWCNTs. The inter-tube distance was set to be 10 ? to minimize the periodic image interactions for all calculations. For comparison, a monolayer graphene sheet with a 5,7-defect (or two 5,7-defects) was constructed to model the hypothetically chirality-changed growth of a SWCNT with an infinite diameter.

    3 Results and discussion

    It is intuitively easy to understand that the cloning growth of a SWCNT just grows the tube length along the direction of the tube axis by continuously adding carbon atoms at the open-end without inducing any defects, while the chirality-changed epitaxial growth for the SWCNTs inevitably induces topological defects at the tube junction position to change the tube chirality of (n, m) to (n′, m′). As for the chirality-changed epitaxial growth considered in this work, we only focus on the cases in which the tube diameters do not change too much, for example from (n, m) to (n±Δ, m?Δ) withΔ=1 and 2, since it is very difficult to enlarge or reduce the tube diameters in a typically epitaxial growth without catalyst at the open-end. As well demonstrated, at least one and two pentagon-heptagon topological defects have to be introduced for the (n,m)/(n±1, m?1) and (n,m)/(n±2, m?2) tube junctions, respectively[33]. In energetics, the topological 5,7-defects in sp2carbon materials like graphene and carbon nanotubes contribute the minimum energy gain to the total energy due to the zero net curvature change, which has been proved to be the dominant and thermodynamically favorable defects in the SWCNT junctions[34-36]and graphene grain boundary[37-39]. Following the rule proposed in ref. 31, the tube junctions of different SWCNTs with chirality changed from (n,m) to (n±1, m?1) and from (n, m) to (n±1, m?1) can be easily constructed.

    Fig. 1 presents the fully relaxed atomic structures of four tube junctions showing the chirality change for (10, 1) and (7, 5) tubes from (n, m) to (n±1, m?1). As clearly shown in Fig. 1, by inducing only one 5,7-defect, the tube chirality changes from initial (n, m) to (n±1, m?1) without an obvious diameter change for both (10, 1) and (7, 5) tubes. And such a topological 5,7-defect connects two tube segments with different chiralities, thus forming a tube junction.

    We expect that different SWCNTs need different energies to change their chiralities during the epitaxial growth. In order to check this, the energies needed for changing the chiralities of different SWCNTs are calculated based on the following equation,

    ΔE=(2×E(n1,m1)→(n2,m2)-E(n1,m1)-E(n2,m2))/2

    (1)

    where,E(n1,m1)→(n2,m2) is the total energy of an isolated tube junction of (n1,m1)/(n2,m2),E(n1,m1) andE(n2,m2) are the total energies of the correponding and defect-free SWCNTs with chiralities of (n1,m1) and (n2,m2), respectively. Therefore, the enengy (ΔE) needed for the chirality-unchanged epitaxial growth of a SWCNT should be zero according to Eq. (1). It should be also noted that the number of C and H atoms in an isolated tube junction should be exactly the same as the half number of C and H atoms in two corresponding SWCNTs to obtain the accurate energy change. For example, the number of C and H atoms in junction of (10, 1)/(11, 0) is 176 and 22, respectively, and the number of total C and H atoms of two correspondingly separated tubes (10, 1) and (11, 0) is 88 and 11, respectively.

    Fig. 1 Schematic cluster models of tube junctions showing the tube chirality changed from (a) (10, 1) to (11, 0), (b) (10, 1)to (9, 2), (c) (7, 5) to (8, 4), and (d) (7, 5) to (6, 6). The open-ends of all junctions are terminated by H atoms. The C and H atoms are represented by purple and blue balls, respectively. Those C atoms at the 5,7-defect sites are highlighted in green.

    Fig. 2 presents a plot of the calculated energy needed (ΔE) to change chirality from (n, m) to (n±1, m?1) for different SWCNTs as a function of the inverse of tube diameter. As we can see, the chirality change from (n, m) to (n±1, m?1) for different SWCNTs is an endothermal process due to an indispensable 5,7-defect introduced in this case. Clearly, the calculatedΔE for different SWCNTs due to the chirality change from (n, m) to (n±1, m?1) increases linearly with increasing the tube diameter, ranging from 1.1 eV for (5, 1) tube to 6.3 eV for an infinite large tube (graphene). This value (6.3 eV) of the formation energy of a single 5,7-defect in graphene is comparable to that in the previous work (7.5 eV)[37]. In order to better understand the nearly linear relationship between energies and tube diameters, we fit all the data plotted in Fig. 2 and obtained a linear relationship asE=-19.7/d+5.4. This is reasonable considering the fact that a small SWCNT has a larger curvature than the large ones, resulting in a larger formation energy and lower stability[9], thus it is more energetically feasible to introduce a 5,7-defect in small SWCNTs. Moreover, it is found that the energies needed to change the chirality of nAC tubes from (n, m) to (n±1, m?1) are about 0.2- 0.4 eV larger than those of nZZ ones with comparable diameters (Fig. 2), implying that it is more energetically difficult to change the chirality from (n, m) to (n±1, m?1) for nAC than nZZ tubes. It seems hard to understand why different tubes with the same diameter show different difficulties in changing their chiralities from (n, m) to (n±1, m?1). Specifically, why a nZZ tube needs less energy to change its chirality than a nAC tube with almost the same diameter during the epitaxial growth, considering that only one 5,7-defect is introduced and the tube diameter does not change too much in this case.

    Fig. 2 The calculated energy change (ΔE in eV) of different SWCNTs due to the chirality change from (n, m) to (n±1, m?1)as a function of the inverse of tube diameter. Open triangle (square)and plus sign (cross) symbols indicate the change of tube chirality from (n, m) to (n+1, m-1) and (n-1, m+1), respectively.The data of ΔE (the formation energy of a 5, 7-defect) for graphene is also included in this figure to represent the SWCNT with an infinite diameter.

    For example, the diameter of (10, 1) tube changes from 8.3 to 8.6 and 8.0 ? after its chirality is changed to (11, 0) and (9, 2), respectively, and the (7, 5) tube changes its diameter from 8.2 to 8.3 and 8.1 ? after its chirality is changed to (8, 4) and (6, 6), respectively. However, as we can see from the atomic structures of tube junctions shown in Fig. 3a-d, the 5,7-defects show different orientations with the tube axis for different junctions, resulting in different intersection angles (θ) between the defect orientation and the tube axis. For example, the intersection angles for (9, 1) /(10, 0), (9, 1)/(8, 2), (7, 4)/(8, 3), and (7, 4)/(6, 5) junctions are 66°, 61°, 83°, and 79°, respectively. Fig. 3e shows the intersection angles (θ) of all (n, m) /(n±1, m?1) tube junctions involved in this work as a function of the inverse of tube diameter. Interestingly, it is found that the intersection angles of different (n, m)/(n±1, m?1) tube junctions for the nAC tubes are about 10°- 20° larger than those for the nZZ tubes. Therefore, more energies needed to change the chirality for nAC tubes than those for nZZ tubes (Fig. 2) can be ascribed to the largerθfor nAC tubes. Actually, it was reported that the larger intersection angle between the 5,7-defect and tube axis results in higher formation energies of the 5,7-defects in SWCNTs[40].

    Fig. 3 The fully relaxed atomic structures of four tube junctions showing the tube chirality change from (a) (9, 1) to (10, 0), (b) (9, 1)to (8, 2), (c) (7, 4) to (8, 3), and (d) (7, 4) to (6, 5), the intersection angles (θ) between the orientation of 5, 7-defect and the tube axis are highlighted by red dashed lines. (e) The plot of intersection angles of all (n, m) /(n±1, m?1) tube junctions for nZZ (nAC) tubes as function of the inverse of tube diameter.

    In addition to the cases of chirality change from (n, m) to (n±1, m?1) as discussed above, the chirality-changed epitaxial growth of SWCNTs with other chirality changes might occur if the energy inputs are large enough by severely altering the epitaxial growth condition. For the case of chirality change from (n, m) to (n±2, m?2), in which the diameter change of a (n, m) tube is larger than that in the case of chirality change from (n, m) to (n±1, m?1) and in consequence a larger energy is expected. In order to check this, we further calculate the energies needed to change the chirality from (n, m) to (n±2, m?2) for different SWCNTs. As reported in ref. 32 that at least two 5,7-defects have to be introduced into a SWCNT to change its chirality from (n, m) to (n±2, m?2). Following this speculation, we constructed the corresponding SWCNT junctions with the tube chirality change from (n, m) to (n±2, m?2) for different SWCNTs. Fig. 4 shows fully relaxed atomic structures of four SWCNT tube junctions showing the chirality change from (7, 3) to (9, 1). Obviouly, two 5,7-defects introduced to change tube chirality from (n, m) to (n±2, m?2) can be separated by different numbers of 6-member ring. As we can see in Fig. 4, the required energy for chirality change in this case increases with increasing the number of 6-member ring between two 5,7-defects. Thus, the tube junction with two adjacent 5,7-defects is the one that is most energetically stable. Actually, it was reported that the formation energy of a pair of 5,7-defects in SWCNTs increases with increasing the spacing distance between two 5,7-defects[41],well consistent with our theoretical results. Therefore, tube junctions with two adjacent 5,7-defects were used to study the energy changes of different SWCNTs with the chirality changed from (n, m) to (n±2, m?2) in this work.

    Fig. 4 The energy needed to change a (7, 3) to a (9, 1) tube as a function of the number of 6-member rings between two 5, 7-defects.The fully relaxed atomic structures of different junctions with different numbers of 6-member rings between two 5, 7-defects are also presented. The open-ends of all junctions are terminated by H atoms. The C and H atoms are represented by purple and blue balls, respectively. The C atoms at the 5, 7-defect sites are highlighted in green.

    Similar to the case of chirality change from (n, m) to (n±1, m?1), the chirality change of SWCNTs from (n, m) to (n±2, m?2) are also found to be an endothermal process and the calculated energy changes of different SWCNTs considered in this work decrease linearly with decreasing the tube diameters (Fig. 5), ranging from 2.07 to 4.4 eV for tubes with diameters ranging from 5.7 to 10.5 ?. This linear relainship between the energy changes and tube diameters shown in Fig. 5 can be further described by an analytic linear equation as, Energy=-66.42/d+12.01. Interestingly, it is found that the required energies for chirality change of different SWCNTs from (n, m) to (n±2, m?2) are less than two times of those for SWCNTs with a similar diameter from (n, m) to (n±1, m?1). Considering that one 5,7-defect and two 5,7-defects are indispensably introduced to change the chirality of different SWCNTs from (n, m) to (n±1, m?1) and from (n, m) to (n±2, m?2), respectively, it seems that energy requirement for the former case should be about two times of that for the latter. However, it should be noted that two isolated 5,7-defects cause a larger lattice strain in SWCNTs than two adjacent 5,7-defects do. In fact, we found that the total energy of tube junction with two adjacent 5,7-defects is remarkably lower than that with two isolated ones for the same tube junction. For example, the total energy of the tube junction of (7, 3)/(9, 1) with two 5,7-defects separated by three 6-member rings was calculated to be 3.3 eV higher than that of the most stable one with two adjacent 5,7-defects (Fig. 4).

    Fig. 5 The calculated energy changes (ΔE in eV) of different SWCNTs due to the chirality change from (n, m) to (n±2, m?2)as a function of the inverse of tube diameter. The data of ΔE(the formation energy of a pair of neighboring 5, 7-defects)for graphene is also shown to represent the SWCNT with an infinite diameter.

    4 Conclusions

    In contrast to the chirality-unchanged epitaxial growth, the energy required is found to be several electron volts for the chirality-changed epitaxial growth of SWCNTs based on the DFT calculations. The energy needed for a chirality change of different SWCNTs decreases linearly with decreasing tube diameter. More energy is required for the nAC tubes to change the chirality from (n, m) to (n±1, m?1) than those for the nZZ tubes with the comparable diameters, due to a larger formation energy for 5,7-defects introduced in the nAC tubes. It is also found that the energy requirement to change the tube chirality from (n, m) to (n±2, m?2) is less than two times of those in the case of chirality change from (n, m) to (n±1, m?1) for SWCNTs with the comparable diameters, due to a less strain induced by two adjacent 5,7-defects than that by two isolated 5,7-defects. We believe that our DFT results may provide a better understanding on the chirality changed epitaxial growth of SWCNTs and guide the synthesis of single-chirality SWCNTs by an epitaxial growth in the future.

    Acknowledgements

    We acknowledge the support from Shenyang Supercomputing Center (CAS). Part of this work was carried out at National Supercomputer Center in Tianjin, and the calculations were performed on TianHe-1(A). We also acknowledge Dr. Bi-Lu Liu of University of Southern California for fruitful discussions.

    [1] Jorio A, Dresselhaus G, Dresselhaus M S. Carbon Nanotubes : Synthesis, Structure, Properties, and Applications[M]. Springer, Berlin, New York, 2001.

    [2] Saito R, Fujita M, Dresselhaus G, et al. Electronic structure of chiral graphene tubules[J]. Appl Phys Lett, 1992, 60(18): 2204-2206.

    [3] Tans S J, Verschueren A R M, Dekker C. Room-temperature transistor based on a single carbon nanotube[J]. Nature, 1998, 393(6680): 49-52.

    [4] Shulaker M M, Hills G, Patil N, et al. Carbon nanotube computer[J]. Nature, 2013, 501(7468): 526-530.

    [5] Wu Z C, Chen Z H, Du X, et al. Transparent, conductive carbon nanotube films[J]. Science, 2004, 305(5688): 1273-1276.

    [6] Deheer W A, Chatelain A, Ugarte D. A carbon nanotube field-emission electron source[J]. Science, 1995, 270(5239): 1179-1180.

    [7] Kreupl F, Graham A P, Duesberg G S, et al. Carbon nanotubes in interconnect applications[J]. Microelectron Eng, 2002, 64(1-4): 399-408.

    [8] Cambre S, Wenseleers W, Goovaerts E, et al. Determination of the metallic/semiconducting ratio in bulk single-wall carbon nanotube samples by cobalt porphyrin probe electron paramagnetic resonance spectroscopy[J]. ACS Nano, 2010, 4(11): 6717-6724.

    [9] Kato K, SaitoS, Geometries. electronic structures and energetics of small-diameter single-walled carbon nanotubes[J]. Physica E, 2011, 43(3): 669-672.

    [10] Wang H, Yuan Y, Wei L, et al. Catalysts for chirality selective synthesis of single-walled carbon nanotubes[J]. Carbon, 2015, 81: 1-19.

    [11] Loebick C Z, Podila R, Reppert J, et al. Selective synthesis of subnanometer diameter semiconducting single-walled carbon nanotubes[J]. J Am Chem Soc, 2010, 132(32): 11125-11131.

    [12] Li X L, Tu X M, Zaric S, et al. Selective synthesis combined with chemical separation of single-walled carbon nanotubes for chirality selection[J]. J Am Chem Soc, 2007, 129(51): 15770-15771.

    [13] He M S, Jiang H, Liu B L, et al. Chiral-selective growth of single-walled carbon nanotubes on lattice-mismatched epitaxial cobalt nanoparticles[J]. Sci Rep, 2013, 3: 1460(1)-1460(7).

    [14] Wang H, Wei L, Ren F, et al. Chiral-selective CoSO4/SiO2Catalyst for (9,8) single-walled carbon nanotube growth[J]. Acs Nano, 2013, 7(1): 614-626.

    [15] Dutta D, Chiang W H, Sankaran R M, et al. Epitaxial nucleation model for chiral-selective growth of single-walled carbon nanotubes on bimetallic catalyst surfaces[J]. Carbon, 2012, 50(10): 3766-3773.

    [16] Fouquet M, Bayer B C, Esconjauregui S, et al. Highly chiral-selective growth of single-walled carbon nanotubes with a simple monometallic Co catalyst[J]. Phys Rev B, 2012, 85: 235411(1)-235411(7).

    [17] Liu B L, Ren W C, Li S S, et al. High temperature selective growth of single-walled carbon nanotubes with a narrow chirality distribution from a CoPt bimetallic catalyst[J]. Chem Commu, 2012, 48(18): 2409-2411.

    [18] He M, Chernov A I, Fedotov P V, et al. Predominant (6,5) single-walled carbon nanotube growth on a copper-promoted iron catalyst[J]. J Am Chem Soc, 2010, 132(40): 13994-13996.

    [19] Wei L, Wang B, Goh T H, et al. Selective enrichment of (6,5) and (8,3) single-walled carbon nanotubes via cosurfactant extraction from narrow (n,m) distribution samples[J]. J Phys Chem B, 2008, 112(10): 2771-2774.

    [20] Wang H, Wang B, Quek X Y, et al. Selective synthesis of (9, 8) single-walled carbon nanotubes on cobalt incorporated TUD-1 catalysts[J]. J Am Chem Soc, 2010, 132(47): 16747-16749.

    [21] Chiang W H, SankaranR M. Linking catalyst composition to chirality distributions of as-grown single-walled carbon nanotubes by tuning NixFe1-x nanoparticles[J]. Nat Mater, 2009, 8(11): 882-886.

    [22] Yang F, Wang X, Zhang D Q, et al. Chirality-specific growth of single-walled carbon nanotubes on solid alloy catalysts[J]. Nature, 2014, 510(7506): 522-524.

    [23] Sanchez-Valencia J R, Dienel T, Groning O, et al. Controlled synthesis of single-chirality carbon nanotubes[J]. Nature, 2014, 512(7512): 61-64.

    [24] Lei Z X, Liu J, Wang J B, et al. The effects of catalyst structure and morphology on the growth of carbon nanotubes[J]. New Carbon Materials, 2003, 18(4): 271-276.

    [25] Song J L, Wang L, Feng S A, et al. Growth of carbon nanotubes by the catalytic decomposition of methane over Fe-Mo/Al2O3catalyst: Effect of temperature on tube structure[J]. New Carbon Materials, 2009, 24(4): 307-313.

    [26] Yao Y G, Feng C Q, Zhang J, et al. “Cloning” of single-walled carbon nanotubes via open-end growth mechanism[J]. Nano Lett, 2009, 9(4): 1673-1677.

    [27] Liu J, Wang C, Tu X M, et al. Chirality-controlled synthesis of single-wall carbon nanotubes using vapour-phase epitaxy[J]. Nat Commun, 2012, 3: 1199(1)-1199(7).

    [28] Yao Y G, Li Q W, Zhang J, et al. Temperature-mediated growth of single-walled carbon-nanotube intramolecular junctions[J]. Nat Mater, 2007, 6(4): 283-286.

    [29] Kresse G, Furthmuller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set[J]. Phys Rev B, 1996, 54(16): 11169-11186.

    [30] Perdew J P, Burke K, Ernzerhof M. Generalized gradient approximation made simple[J]. Phys Rev Lett, 1996, 77(18): 3865-3868.

    [31] Blochl P E. Projector augmented-wave method[J]. Phys Rev B, 1994, 50(24): 17953-17979.

    [32] Kresse G, Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method[J]. Phys Rev B, 1999, 59(3): 1758-1775.

    [33] Saito R, Dresselhaus G, Dresselhaus M S. Physical Properties of Carbon Nanotubes [M]. Imperial College Press, London, 1996.

    [34] Chico L, Crespi V H, Benedict L X, et al. Pure carbon nanoscale devices: Nanotube heterojunctions[J]. Phys Rev Lett, 1996, 76(6): 971-974.

    [35] Ouyang M, Huang J L, Cheung C L, et al. Atomically resolved single-walled carbon nanotube intramolecular junctions[J]. Science, 2001, 291(5501): 97-100.

    [36] Wei D C, LiuY Q. The intramolecular junctions of carbon nanotubes[J]. Adv Mater, 2008, 20(15): 2815-2841.

    [37] Yazyev O V, LouieS G. Topological defects in graphene: Dislocations and grain boundaries[J]. Phys Rev B, 2010, 81: 195420(1)-195420(7).

    [38] Grantab R, Shenoy V B, Ruoff R S. Anomalous strength characteristics of tilt grain boundaries in graphene[J]. Science, 2010, 330(6006): 946-948.

    [39] Huang P Y, Ruiz-Vargas C S, van der Zande A M, et al. Grains and grain boundaries in single-layer graphene atomic patchwork quilts[J]. Nature, 2011, 469(7330): 389-392.

    [40] Zhou L G, Shi S Q. Formation energy of stone-wales defects in carbon nanotubes[J]. Appl Phys Lett, 2003, 83(6): 1222-1224.

    [41] Yuan Q H, Xu Z P, Yakobson B I, et al. Efficient defect healing in catalytic carbon nanotube growth[J]. Phys Rev Lett, 2012, 108: 245505(1)-245505(5).

    Changing the chirality of single-wall carbon nanotubes during epitaxial growth: A density functional theory study

    SHEN Wei, LI Feng, LIU Chang, YIN Li-chang

    (ShenyangNationalLaboratoryforMaterialsScience,InstituteofMetalResearch,ChineseAcademyofSciences,Shenyang110016,China)

    The energetic of the change in the chirality of single-wall carbon nanotubes (SWCNTs)during epitaxial growth from (n, m) to (n±Δ, mΔ)(Δ=1, 2) was investigated by density functional theory calculations. The calculated energies for changing the chirality of different SWCNTs show a nearly linear decrease with decreasing tube diameter. In the case ofΔ=1, more energy input is needed for near armchair (nAC) SWCNTs to change their chiralities than those for near zigzag (nZZ) SWCNTs with comparable diameters, due to the larger formation energies of pentagon-heptagon defects (5,7-defects) introduced in the nAC-SWCNTs. These larger formation energies for the nAC-SWCNTs come from the larger angles between the orientation of a 5,7-defect and the tube axis than those for nZZ-SWCNTs. The topological connection oftwo adjacent 5,7-defects, which is indispensable for changing the chirality during growth in the case ofΔ=2, is found to be energetically most stable. The energies needed to change chirality in the case ofΔ=2 are calculated to be less than twice those in the case ofΔ=1 for SWCNTs with comparable diameters. These results may help us understand the change in chirality during the epitaxial growth of SWCNTs and guide the future synthesis of SWCNTs with a single-chirality .

    Single-walled carbon nanotube; Chirality-changed; Pentagon-heptagon defect; Epitaxial growth; Density functional theory

    National Natural Science Foundation of China (51272257, 51202255, 51472249).

    YIN Li-chang, Associate Professor. E-mail: lcyin@imr.ac.cn

    國家自然科學(xué)基金項(xiàng)目(51272257, 51202255, 51472249).

    尹利長,副研究員. E-mail: lcyin@imr.ac.cn

    作者介紹:諶 為,博士研究生. E-mail: weichen@imr.ac.cn

    1007-8827(2016)05-0525-07

    TQ127.1+1

    A

    Authorintroduction: SHEN Wei, Ph. D Candidate. E-mail: weichen@imr.ac.cn

    10.1016/S1872-5805(16)60030-6

    Receiveddate: 2016-07-29;Reviseddate: 2016-10-03

    English edition available online ScienceDirect (http:www.sciencedirect.comsciencejournal18725805).

    猜你喜歡
    生長
    野蠻生長
    碗蓮生長記
    小讀者(2021年2期)2021-03-29 05:03:48
    生長的樹
    自由生長的家
    美是不斷生長的
    快速生長劑
    共享出行不再“野蠻生長”
    生長在哪里的啟示
    野蠻生長
    NBA特刊(2018年21期)2018-11-24 02:48:04
    生長
    文苑(2018年22期)2018-11-19 02:54:14
    久久久久国产一级毛片高清牌| 嫩草影视91久久| 国产爱豆传媒在线观看 | 久久精品91蜜桃| 欧美又色又爽又黄视频| 黄片小视频在线播放| 久久精品国产亚洲av高清一级| 精品欧美国产一区二区三| 免费看日本二区| 夜夜躁狠狠躁天天躁| 曰老女人黄片| 国产不卡一卡二| 黄色视频,在线免费观看| 国产黄a三级三级三级人| 最好的美女福利视频网| 日韩欧美 国产精品| 国产熟女xx| 特大巨黑吊av在线直播 | 91大片在线观看| 这个男人来自地球电影免费观看| 天堂√8在线中文| 亚洲欧美日韩无卡精品| 一区二区三区激情视频| 无限看片的www在线观看| 日韩大尺度精品在线看网址| 特大巨黑吊av在线直播 | 此物有八面人人有两片| 国产精品av久久久久免费| 精品第一国产精品| 国产麻豆成人av免费视频| 日韩av在线大香蕉| 久久久国产精品麻豆| 婷婷六月久久综合丁香| 悠悠久久av| 不卡av一区二区三区| 99精品在免费线老司机午夜| 法律面前人人平等表现在哪些方面| 久久这里只有精品19| 成年人黄色毛片网站| 亚洲精品粉嫩美女一区| 亚洲精品久久国产高清桃花| 国产aⅴ精品一区二区三区波| 欧美精品啪啪一区二区三区| 精品国产一区二区三区四区第35| 国产精品爽爽va在线观看网站 | 国产精品野战在线观看| 自线自在国产av| 丁香欧美五月| 男女之事视频高清在线观看| 精品欧美一区二区三区在线| 亚洲国产日韩欧美精品在线观看 | 免费看日本二区| 国产精品电影一区二区三区| 天天躁夜夜躁狠狠躁躁| or卡值多少钱| 首页视频小说图片口味搜索| 欧美日韩瑟瑟在线播放| 亚洲精品久久成人aⅴ小说| 亚洲七黄色美女视频| 国产视频内射| 天天一区二区日本电影三级| avwww免费| 一夜夜www| 国产成人欧美在线观看| 色老头精品视频在线观看| 波多野结衣av一区二区av| 色综合站精品国产| 亚洲欧美激情综合另类| 国产亚洲精品一区二区www| ponron亚洲| 亚洲熟妇熟女久久| 国产国语露脸激情在线看| 免费在线观看日本一区| 久久精品人妻少妇| 国产精品av久久久久免费| 免费在线观看亚洲国产| av电影中文网址| av中文乱码字幕在线| 淫秽高清视频在线观看| 非洲黑人性xxxx精品又粗又长| 99精品在免费线老司机午夜| av中文乱码字幕在线| 黑丝袜美女国产一区| 国内精品久久久久久久电影| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久久久亚洲av鲁大| 日韩欧美免费精品| 熟女少妇亚洲综合色aaa.| 香蕉久久夜色| 国产高清激情床上av| 老鸭窝网址在线观看| 国产久久久一区二区三区| av免费在线观看网站| 香蕉av资源在线| 在线观看免费日韩欧美大片| 精品高清国产在线一区| 久久国产精品男人的天堂亚洲| 成人午夜高清在线视频 | 久久久水蜜桃国产精品网| 女警被强在线播放| www日本在线高清视频| 国产又色又爽无遮挡免费看| 国产亚洲av嫩草精品影院| 亚洲欧美日韩高清在线视频| 波多野结衣高清无吗| 女性生殖器流出的白浆| 制服丝袜大香蕉在线| 免费在线观看视频国产中文字幕亚洲| 脱女人内裤的视频| 欧美乱色亚洲激情| 热99re8久久精品国产| 人成视频在线观看免费观看| e午夜精品久久久久久久| a级一级毛片免费在线观看| 日韩精品中文字幕看吧| 日韩人妻高清精品专区| 男人舔女人下体高潮全视频| 国产精品一区二区三区四区久久| 亚洲av美国av| 校园人妻丝袜中文字幕| 日本一二三区视频观看| 成熟少妇高潮喷水视频| 亚洲av中文字字幕乱码综合| 精品不卡国产一区二区三区| 成人三级黄色视频| 亚洲精品国产成人久久av| 禁无遮挡网站| 精品99又大又爽又粗少妇毛片| 国产亚洲精品综合一区在线观看| 免费观看在线日韩| 成人国产麻豆网| 18禁在线播放成人免费| 简卡轻食公司| 国产又黄又爽又无遮挡在线| 国产成人freesex在线 | 精品日产1卡2卡| 亚洲专区国产一区二区| 久久精品影院6| 麻豆乱淫一区二区| av福利片在线观看| 狂野欧美激情性xxxx在线观看| 99久国产av精品| 国产精品亚洲一级av第二区| 噜噜噜噜噜久久久久久91| 免费看日本二区| 日本成人三级电影网站| 淫秽高清视频在线观看| 在线播放国产精品三级| 国产精品久久久久久亚洲av鲁大| 日本a在线网址| 日本黄色片子视频| 少妇的逼水好多| 亚洲精品日韩在线中文字幕 | 国产成人a∨麻豆精品| 国产精品久久电影中文字幕| 欧美激情久久久久久爽电影| 国产一区二区在线观看日韩| 欧美在线一区亚洲| 欧美成人免费av一区二区三区| 日韩国内少妇激情av| 久久精品人妻少妇| 美女大奶头视频| 听说在线观看完整版免费高清| 国产精品美女特级片免费视频播放器| 最近最新中文字幕大全电影3| av天堂中文字幕网| 男女下面进入的视频免费午夜| 成人特级av手机在线观看| 美女cb高潮喷水在线观看| 国产av一区在线观看免费| 国产淫片久久久久久久久| 中国美女看黄片| 日韩欧美免费精品| 免费搜索国产男女视频| 免费黄网站久久成人精品| 成人高潮视频无遮挡免费网站| 国产男靠女视频免费网站| 成人漫画全彩无遮挡| 1024手机看黄色片| 国产成人a区在线观看| 色哟哟哟哟哟哟| 免费高清视频大片| 男人的好看免费观看在线视频| 在线播放国产精品三级| 亚洲,欧美,日韩| 国产午夜精品久久久久久一区二区三区 | 1024手机看黄色片| 国产aⅴ精品一区二区三区波| 天天一区二区日本电影三级| 日本与韩国留学比较| a级毛片a级免费在线| 国产黄色视频一区二区在线观看 | 欧美zozozo另类| 国产v大片淫在线免费观看| 一个人看视频在线观看www免费| av在线观看视频网站免费| 嫩草影院精品99| 乱系列少妇在线播放| 久久久精品大字幕| 全区人妻精品视频| 亚洲18禁久久av| 午夜爱爱视频在线播放| 国产欧美日韩精品一区二区| 久久久久久九九精品二区国产| 国产一区二区在线av高清观看| 欧美在线一区亚洲| 一区福利在线观看| 九九热线精品视视频播放| 亚洲四区av| 精品人妻一区二区三区麻豆 | 日韩av不卡免费在线播放| 久久久久久久午夜电影| 能在线免费观看的黄片| 亚洲中文日韩欧美视频| 波多野结衣高清无吗| 成人漫画全彩无遮挡| 美女高潮的动态| 伦精品一区二区三区| 日韩av不卡免费在线播放| 久久久精品94久久精品| 国产日本99.免费观看| 中文资源天堂在线| 成人三级黄色视频| 中国美女看黄片| 又黄又爽又刺激的免费视频.| 日产精品乱码卡一卡2卡三| 亚洲va在线va天堂va国产| 亚洲国产日韩欧美精品在线观看| 哪里可以看免费的av片| 国产在线精品亚洲第一网站| 老熟妇乱子伦视频在线观看| 国产伦一二天堂av在线观看| 日本成人三级电影网站| 少妇的逼好多水| 成人鲁丝片一二三区免费| 99久久精品一区二区三区| 欧美三级亚洲精品| 免费电影在线观看免费观看| 久久午夜亚洲精品久久| 日日摸夜夜添夜夜添小说| 亚洲精品在线观看二区| 观看美女的网站| 一级黄片播放器| 老司机影院成人| 伦精品一区二区三区| 午夜影院日韩av| 激情 狠狠 欧美| 欧美成人一区二区免费高清观看| 国产乱人偷精品视频| 精品一区二区三区视频在线观看免费| 成人国产麻豆网| 欧美zozozo另类| 成年av动漫网址| 欧美区成人在线视频| 国产国拍精品亚洲av在线观看| 国产精品美女特级片免费视频播放器| 日韩欧美在线乱码| eeuss影院久久| 久久午夜福利片| 晚上一个人看的免费电影| 日日摸夜夜添夜夜添小说| 国产免费一级a男人的天堂| 久久久精品欧美日韩精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人av在线免费| 搡老岳熟女国产| 少妇猛男粗大的猛烈进出视频 | or卡值多少钱| 女人被狂操c到高潮| 搞女人的毛片| 国产三级在线视频| 欧美激情在线99| 亚洲精品乱码久久久v下载方式| 熟女电影av网| 变态另类丝袜制服| 97超级碰碰碰精品色视频在线观看| 国产探花极品一区二区| 搡女人真爽免费视频火全软件 | 联通29元200g的流量卡| 午夜福利在线在线| 欧美激情在线99| 亚洲人成网站高清观看| 亚洲精品在线观看二区| 性插视频无遮挡在线免费观看| 亚洲欧美日韩高清在线视频| 一区二区三区高清视频在线| 在线观看美女被高潮喷水网站| 精品人妻视频免费看| 精品乱码久久久久久99久播| 性欧美人与动物交配| 国产亚洲精品av在线| 男女做爰动态图高潮gif福利片| 日韩 亚洲 欧美在线| 成年版毛片免费区| or卡值多少钱| 国产乱人偷精品视频| 国产乱人视频| 亚洲欧美日韩卡通动漫| 人妻丰满熟妇av一区二区三区| 亚洲乱码一区二区免费版| 91久久精品国产一区二区成人| 色尼玛亚洲综合影院| 人妻夜夜爽99麻豆av| 久久久精品欧美日韩精品| 老司机福利观看| 欧美日韩一区二区视频在线观看视频在线 | 久久草成人影院| 国产爱豆传媒在线观看| 亚洲国产色片| 99视频精品全部免费 在线| 99国产精品一区二区蜜桃av| 看免费成人av毛片| 精品免费久久久久久久清纯| 亚洲av中文字字幕乱码综合| 欧美一区二区国产精品久久精品| 午夜视频国产福利| 精品久久国产蜜桃| 亚洲在线自拍视频| 亚州av有码| 久久久久精品国产欧美久久久| 欧美日韩精品成人综合77777| 久久人人爽人人片av| 国产精品一区二区性色av| 亚洲精品国产av成人精品 | 成人av在线播放网站| 国产亚洲精品久久久久久毛片| 岛国在线免费视频观看| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站| 在线国产一区二区在线| 亚洲最大成人av| 国产欧美日韩精品一区二区| 久久久久久国产a免费观看| 大又大粗又爽又黄少妇毛片口| 97超视频在线观看视频| 国产大屁股一区二区在线视频| 日韩精品有码人妻一区| 亚洲第一电影网av| 色吧在线观看| 十八禁国产超污无遮挡网站| 丰满人妻一区二区三区视频av| 精品久久久久久久久久久久久| 国产高清激情床上av| 国产麻豆成人av免费视频| 国产伦一二天堂av在线观看| 久久久久久久久中文| 日韩精品青青久久久久久| 免费大片18禁| 别揉我奶头 嗯啊视频| 国产爱豆传媒在线观看| videossex国产| 欧美日韩国产亚洲二区| 亚洲国产欧洲综合997久久,| 精品99又大又爽又粗少妇毛片| 男人舔奶头视频| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 欧美精品国产亚洲| 亚洲图色成人| 午夜影院日韩av| 国产片特级美女逼逼视频| 久久久久久久久久久丰满| 国产高清不卡午夜福利| 国产精品免费一区二区三区在线| 国产精华一区二区三区| 99精品在免费线老司机午夜| 欧美一区二区精品小视频在线| 国产成年人精品一区二区| 国产麻豆成人av免费视频| 黄色一级大片看看| 国产一区二区在线av高清观看| 精品无人区乱码1区二区| 亚洲最大成人中文| 日日撸夜夜添| 亚洲人与动物交配视频| 精品少妇黑人巨大在线播放 | 免费在线观看成人毛片| 国产精品99久久久久久久久| 国产精品av视频在线免费观看| 亚洲精品国产av成人精品 | 国产一区亚洲一区在线观看| 亚洲成人久久爱视频| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 俺也久久电影网| 亚洲av中文字字幕乱码综合| 美女免费视频网站| 18禁裸乳无遮挡免费网站照片| 天堂av国产一区二区熟女人妻| 国产 一区精品| 18禁裸乳无遮挡免费网站照片| 联通29元200g的流量卡| 又粗又爽又猛毛片免费看| 精品免费久久久久久久清纯| 精品久久久久久久久久久久久| 国产高清视频在线播放一区| 亚洲性夜色夜夜综合| 内地一区二区视频在线| 国内精品美女久久久久久| 欧美性猛交╳xxx乱大交人| 亚洲国产精品久久男人天堂| 人妻夜夜爽99麻豆av| 国产免费一级a男人的天堂| 久久久精品欧美日韩精品| 国产aⅴ精品一区二区三区波| 看非洲黑人一级黄片| 日本a在线网址| 老师上课跳d突然被开到最大视频| 午夜福利高清视频| 成人午夜高清在线视频| 淫妇啪啪啪对白视频| 国产三级在线视频| 97热精品久久久久久| 天堂av国产一区二区熟女人妻| 欧美日韩精品成人综合77777| 久久九九热精品免费| 国产精品不卡视频一区二区| 久久久成人免费电影| 国产美女午夜福利| 日韩,欧美,国产一区二区三区 | 日本一本二区三区精品| 婷婷精品国产亚洲av在线| 欧美在线一区亚洲| 国产av麻豆久久久久久久| 国产精品av视频在线免费观看| 97热精品久久久久久| 十八禁国产超污无遮挡网站| 国产三级中文精品| 亚洲国产精品成人综合色| 精品久久久噜噜| 色在线成人网| ponron亚洲| 噜噜噜噜噜久久久久久91| 亚洲第一电影网av| av天堂在线播放| 午夜福利在线观看吧| 国产伦精品一区二区三区四那| 少妇被粗大猛烈的视频| 亚洲国产色片| 老熟妇乱子伦视频在线观看| 国产一区二区在线av高清观看| 午夜精品国产一区二区电影 | 黄色视频,在线免费观看| 日本三级黄在线观看| 男人舔奶头视频| 成人无遮挡网站| 中文字幕av在线有码专区| 小说图片视频综合网站| 乱人视频在线观看| 久久精品国产鲁丝片午夜精品| 国产老妇女一区| 国产国拍精品亚洲av在线观看| 欧美3d第一页| 国产老妇女一区| 91久久精品国产一区二区三区| 欧美极品一区二区三区四区| 成年女人看的毛片在线观看| 亚洲精品国产成人久久av| 亚洲美女搞黄在线观看 | 欧美xxxx黑人xx丫x性爽| 一级毛片aaaaaa免费看小| 日韩中字成人| 亚洲欧美中文字幕日韩二区| 国产精品乱码一区二三区的特点| 久久久久久久午夜电影| 如何舔出高潮| 亚洲第一区二区三区不卡| 偷拍熟女少妇极品色| 国产色爽女视频免费观看| 精品人妻视频免费看| 亚洲成人久久爱视频| 国产精品亚洲一级av第二区| 一夜夜www| 成人精品一区二区免费| 国产精品久久久久久精品电影| 干丝袜人妻中文字幕| 三级国产精品欧美在线观看| 99国产精品一区二区蜜桃av| 久久久欧美国产精品| 亚州av有码| 成年av动漫网址| av中文乱码字幕在线| 18禁裸乳无遮挡免费网站照片| 日韩成人伦理影院| 亚洲国产欧洲综合997久久,| 人妻久久中文字幕网| 91久久精品电影网| 成人午夜高清在线视频| 青春草视频在线免费观看| 久久精品国产亚洲av香蕉五月| 麻豆乱淫一区二区| 久久精品91蜜桃| 美女xxoo啪啪120秒动态图| 久久久久久久亚洲中文字幕| 日本色播在线视频| 午夜爱爱视频在线播放| 欧美最新免费一区二区三区| 91久久精品国产一区二区三区| 亚洲天堂国产精品一区在线| 国产真实伦视频高清在线观看| 国产精品嫩草影院av在线观看| 久久久久国产网址| 一级毛片电影观看 | 身体一侧抽搐| 欧美激情在线99| 搡老熟女国产l中国老女人| 少妇的逼好多水| 真实男女啪啪啪动态图| 3wmmmm亚洲av在线观看| a级毛色黄片| 成人毛片a级毛片在线播放| 国产精品久久视频播放| 欧美激情久久久久久爽电影| 免费大片18禁| 搡老熟女国产l中国老女人| 日韩强制内射视频| 精品久久久久久久人妻蜜臀av| 波多野结衣高清作品| 又黄又爽又刺激的免费视频.| 99久久中文字幕三级久久日本| 亚洲av电影不卡..在线观看| 久久欧美精品欧美久久欧美| 可以在线观看的亚洲视频| 国产免费一级a男人的天堂| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 欧美zozozo另类| 不卡视频在线观看欧美| 偷拍熟女少妇极品色| 在线播放无遮挡| 久久综合国产亚洲精品| 亚洲无线观看免费| 深爱激情五月婷婷| 少妇的逼水好多| 久久精品综合一区二区三区| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久| 乱码一卡2卡4卡精品| 免费av毛片视频| 直男gayav资源| 久久亚洲精品不卡| 欧美成人精品欧美一级黄| 亚洲美女黄片视频| 在线看三级毛片| 亚洲精品在线观看二区| 国产中年淑女户外野战色| 日韩av不卡免费在线播放| 乱系列少妇在线播放| 91久久精品电影网| 成人性生交大片免费视频hd| 我的老师免费观看完整版| 一级毛片我不卡| 搡老熟女国产l中国老女人| or卡值多少钱| 国产av麻豆久久久久久久| 长腿黑丝高跟| 国产综合懂色| 国产精品,欧美在线| 尤物成人国产欧美一区二区三区| 亚洲精品成人久久久久久| 少妇熟女欧美另类| 悠悠久久av| 国产美女午夜福利| 日韩av不卡免费在线播放| 国产黄片美女视频| 色播亚洲综合网| 亚洲经典国产精华液单| 18禁在线无遮挡免费观看视频 | 男女之事视频高清在线观看| 久久久久国产精品人妻aⅴ院| 18禁在线播放成人免费| 大香蕉久久网| 天天躁日日操中文字幕| 中文亚洲av片在线观看爽| 麻豆成人午夜福利视频| 日本免费一区二区三区高清不卡| 男人的好看免费观看在线视频| 两个人的视频大全免费| 久99久视频精品免费| av天堂在线播放| 欧美成人精品欧美一级黄| 久久亚洲精品不卡| 青春草视频在线免费观看| 亚洲性夜色夜夜综合| 欧美激情久久久久久爽电影| 最近手机中文字幕大全| 久久精品人妻少妇| 久久草成人影院| 一进一出好大好爽视频| 真实男女啪啪啪动态图| a级毛片免费高清观看在线播放| 黄色配什么色好看| 久久久久九九精品影院| 国产高潮美女av| 欧美丝袜亚洲另类| 国内精品美女久久久久久| 日本a在线网址| 国产黄片美女视频| 给我免费播放毛片高清在线观看| 嫩草影视91久久| 老司机福利观看| 亚洲av美国av| 久久久久久久久中文| а√天堂www在线а√下载| 女同久久另类99精品国产91| 国产淫片久久久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲国产精品久久男人天堂| 亚洲国产欧洲综合997久久,| 成人永久免费在线观看视频| 亚洲五月天丁香| 国产精品综合久久久久久久免费| 如何舔出高潮| 免费电影在线观看免费观看| 青春草视频在线免费观看|