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

    電場(chǎng)作用下邊緣修飾石墨烯納米帶的電磁性質(zhì)

    2020-04-29 10:47:00吳秀解憂宋寶寧曹松張衛(wèi)濤陳立勇

    吳秀 解憂 宋寶寧 曹松 張衛(wèi)濤 陳立勇

    摘要:采用基于密度泛函理論的第一性原理,研究了非金屬原子(B,N)和金屬原子(Na,Mn)邊緣修飾不同寬度zigzag型石墨烯納米帶(ZGNR)的結(jié)構(gòu)穩(wěn)定性和電磁性質(zhì)。并進(jìn)一步研究了外加不同強(qiáng)度垂直電場(chǎng)對(duì)原子邊緣修飾ZGNR電磁性質(zhì)的調(diào)控規(guī)律,所有原子邊緣修飾ZGNR體系均能形成穩(wěn)定結(jié)構(gòu)。在施加垂直于納米帶平面的電場(chǎng)強(qiáng)度后,隨著電場(chǎng)強(qiáng)度的增大,[N-ZG]4體系穩(wěn)定性最好且基本不變,[Mn-ZG]4體系穩(wěn)定性增加,而[B-ZG]4和[Na-ZG]4體系穩(wěn)定性降低;電場(chǎng)促進(jìn)了邊緣修飾原子(B,Na,Mn)的電荷進(jìn)一步向石墨烯C原子轉(zhuǎn)移,阻礙了電荷從石墨烯C原子向邊緣N原子轉(zhuǎn)移;N和Na修飾體系的磁性逐漸減小,而B(niǎo)和Mn修飾體系磁性逐漸增大。研究結(jié)果表明,外加不同強(qiáng)度電場(chǎng)可以有效調(diào)控原子邊緣修飾ZGNR體系的結(jié)構(gòu)穩(wěn)定性、電子性質(zhì)和磁性,能夠?yàn)槭┘{米帶在納米電子領(lǐng)域的應(yīng)用提供理論基礎(chǔ)。

    關(guān)鍵詞:石墨烯納米帶;邊緣修飾;電子結(jié)構(gòu);第一性原理

    中圖分類(lèi)號(hào):O469

    DOI:10.16152/j.cnki.xdxbzr.2020-02-015

    The electronic and magnetic properties of edge-modifiedgraphene nanoribbons in external electric field

    WU Xiu, XIE You, SONG Baoning, CAO Song, ZHANG Weitao, CHEN Liyong

    Abstract: The structure stability and electronic and magnetic properties have been calculated for the non-metal (B, N) and metal (Na, Mn) atoms edge-modified zigzag graphene nanoribbons (ZGNR) with different widths by using the first principle method based on density functional theory. Then the modified regulations of electronic structures are investigated for the atoms edge-modified ZGNR systems in electric fields with different intensities. The stable systems can be formed for all the atoms edge-modified ZGNR structures. With the increasing intensity of electric fields perpendicular to the nanoribbon surface, the stability is almost unchanged for the most stable [N-ZG]4 system, and the stability is increases for the [Mn-ZG]4 system but decreases for the [B-ZG]4 and [Na-ZG]4 systems. The electric field promotes the further transfer of charge from atoms (B, Na, Mn) to the C atoms of graphene in [M-ZG] system, and hinders the transfer of charge from C atom to N atom in [N-ZG] system. The magnetisms? decrease gradually for the [N-ZG] and [Na-ZG] systems while increase gradually for [B-ZG] and [Mn-ZG] systems. The results show that the structure stability, electronic and magnetic properties can be regulated effectively by the electric fields with different intensities for the atoms edge modified ZGNR systems, which can provide a theoretical basis for the application of graphene nanoribbons in the field of nanoeletronics.

    Key words: graphene nanoribbon; edge modification; electronic structures; first-principle

    石墨烯(graphene)是一種由碳原子以sp2雜化軌道組成六角型呈蜂巢晶格的二維碳納米材料。石墨烯具有優(yōu)異的物理化學(xué)性質(zhì),在納米電子器件、新能源材料等方面具有重要的應(yīng)用前景,引起了科研人員的廣泛關(guān)注。對(duì)石墨烯沿不同方向裁剪可以得到準(zhǔn)一維的石墨烯納米帶(GNR)[1-5]。根據(jù)GNR邊緣的不同形狀,可分為zigzag型GNR(ZGNR)和armchair型GNR(AGNR)。不同結(jié)構(gòu)類(lèi)型和不同寬度的GNR可分別表現(xiàn)出半導(dǎo)體性或者金屬性[2-3],GNR的電磁性質(zhì)受其寬度和邊緣形狀的影響[6-7]。同時(shí),通過(guò)摻雜、吸附、缺陷、邊緣修飾、外加電場(chǎng)或磁場(chǎng)等方法能夠進(jìn)一步豐富GNR的電磁性質(zhì)[8-14]。調(diào)控GNR的電磁性質(zhì)對(duì)于擴(kuò)展GNR在納米電子器件領(lǐng)域的應(yīng)用研究具有重要意義。

    前期對(duì)于GNR的研究中,主要是用氫原子進(jìn)行邊緣飽和懸掛鍵。但是,GNR具有較強(qiáng)的邊緣活性,與元素、化學(xué)基團(tuán)等結(jié)合,易于形成局域化的邊緣電子態(tài)。所以,對(duì)GNR進(jìn)行邊緣修飾是調(diào)控GNR性質(zhì)研究的重要手段[15-20]。Gunlycke等人用不同的原子和官能團(tuán)(H,O,OH,NH)對(duì)zigzag型GNR進(jìn)行邊緣飽和[15],發(fā)現(xiàn)ZGNR費(fèi)米能級(jí)附近的電子結(jié)構(gòu)產(chǎn)生明顯改變,可以產(chǎn)生半導(dǎo)體和金屬性。對(duì)armchair型GNR的邊緣用鉀原子修飾后[16],AGNR的電子性質(zhì)由半導(dǎo)體轉(zhuǎn)變?yōu)榻饘傩再|(zhì)。Wi等人研究了單軸拉伸應(yīng)力作用下的邊緣BN對(duì)稱修飾armchair型GNR的力學(xué)和電子性質(zhì)[19],應(yīng)力和帶隙的關(guān)系表明電子結(jié)構(gòu)性質(zhì)對(duì)于單軸拉伸應(yīng)力比較敏感,并主要受納米帶帶寬和BN位置的影響。邊緣BN對(duì)稱修飾,在不改變AGNR力學(xué)性質(zhì)的情況下,可以有效地改變AGNR的電子性質(zhì)。這些研究都表明邊緣修飾對(duì)GNR的穩(wěn)定性及電磁特性有極其重要的影響。在邊緣修飾的基礎(chǔ)上,如何進(jìn)一步調(diào)控石墨烯的性質(zhì),是值得探索的科學(xué)問(wèn)題。研究表明,電場(chǎng)也是調(diào)控石墨烯電磁性質(zhì)的重要手段[12-13]。因此,本文以zigzag型GNR(ZGNR)為研究對(duì)象,把邊緣修飾和電場(chǎng)相結(jié)合,在不同原子M(M=B,N,Na,Mn)邊緣修飾ZGNR的基礎(chǔ)上,對(duì)ZGNR施加垂直電場(chǎng),探索外加不同強(qiáng)度電場(chǎng)對(duì)ZGNR的結(jié)構(gòu)穩(wěn)定性和電磁性質(zhì)的影響規(guī)律。

    1 計(jì)算方法與結(jié)構(gòu)模型

    本文計(jì)算采用基于密度泛函理論的第一性原理軟件VASP[21-22]。 離子與電子間的相互作用選用投影綴加波方法(PAW)來(lái)描述[23], 電子間的交換關(guān)聯(lián)能采用廣義梯度近似(GGA)的PBE(Perdew-Burke-Ernzerhof)泛函理論進(jìn)行處理[24], 用平面波函數(shù)展開(kāi)處理電子波函數(shù),平面波的截?cái)嗄茉O(shè)置為450 eV。 為消除石墨烯納米帶周?chē)溺R像影響, 在垂直于體系方向添加了厚度為0.2×10-8 m的真空層。 通過(guò)Monkorst-Park自動(dòng)生成方法在簡(jiǎn)約布里淵區(qū)中產(chǎn)生1×11×1個(gè)k點(diǎn)[25]。 采用共軛梯度算法弛豫離子到基態(tài), 且離子的弛豫能量收斂標(biāo)準(zhǔn)為1.0×10-4 eV/atom, 作用到每個(gè)原子上的力的收斂標(biāo)準(zhǔn)為0.2×109 V/m。 對(duì)于外加電場(chǎng)的計(jì)算,通過(guò)偶極子方法LDIPOL引入外電場(chǎng),將偶極子置于結(jié)構(gòu)模型的質(zhì)心,即設(shè)置參數(shù)DIPOL=0.5 0.5 0.5。電場(chǎng)值大小由參數(shù)EFIELD控制,本文計(jì)算分別設(shè)置為EFIELD=0.5×1010 V/m,0.1×1011 V/m和0.15×1011 V/m。電場(chǎng)方向由參數(shù)IDIPOL控制,IDIPOL為1,2,3時(shí)分別代表x,y,z方向,本文電場(chǎng)方向垂直于納米帶平面(z方向),即設(shè)置IDIPOL=3。同時(shí),對(duì)于本文研究的邊緣修飾單層石墨烯納米帶來(lái)說(shuō),外加電場(chǎng)垂直于石墨烯納米帶平面(即電場(chǎng)方向與z軸方向平行),而石墨烯納米帶平面在x,y方向上,根據(jù)石墨烯納米帶的對(duì)稱性,電場(chǎng)沿z軸向上和向下的作用效果是一樣的,不會(huì)對(duì)計(jì)算結(jié)果產(chǎn)生影響,所以,本文只計(jì)算了電場(chǎng)方向垂直向上的情況。

    為研究不同原子種類(lèi)對(duì)不同寬度ZGNR的影響,以n(n=3,4,…,8)個(gè)碳原子寬度的ZGNR為研究對(duì)象,分別用非金屬和金屬原子M(M=B,N,Na,Mn)邊緣修飾ZGNR。圖1為不同寬度(3~8個(gè)碳原子)的ZGNR邊緣修飾M原子的復(fù)合體系結(jié)構(gòu),分別記為[M-ZG]3,[M-ZG]4,[M-ZG]5,[M-ZG]6,[M-ZG]7,[M-ZG]8。

    2 結(jié)果與討論

    2.1 結(jié)構(gòu)特性

    對(duì)所有的邊緣修飾[M-ZG]n體系進(jìn)行結(jié)構(gòu)優(yōu)化和馳豫,得到最穩(wěn)定的基態(tài)結(jié)構(gòu)。復(fù)合體系結(jié)構(gòu)的穩(wěn)定性取決于其形成能的大小,邊緣修飾[M-ZG]n體系的形成能定義為:Ef=ET-EZG-EM。其中,Ef為[M-ZG]n體系的形成能;ET為[M-ZG]n體系的總能量;EZG為ZGNR的能量;EM為修飾M原子的能量。表1給出了4種原子邊緣修飾不同寬度ZGNR的[M-ZG]n體系的形成能。負(fù)的形成能Ef表示原子邊緣修飾ZGNR的過(guò)程屬于放熱過(guò)程,能夠自發(fā)形成穩(wěn)定結(jié)構(gòu),且形成能越小,[M-ZG]n體系越穩(wěn)定。由表1可以得出以下幾點(diǎn)結(jié)論:第一,非金屬原子B和N修飾的ZGNR體系([B-ZG]和[N-ZG]),其形成能小于金屬原子Na和Mn修飾的ZGNR體系([Na-ZG]和[Mn-ZG])的形成能,說(shuō)明非金屬原子(B和N)邊緣修飾ZGNR比金屬原子(Na和Mn)邊緣修飾ZGNR更為穩(wěn)定。第二,對(duì)于不同的修飾原子(B,N,Na,Mn),N原子邊緣修飾ZGNR的穩(wěn)定性最好,Na原子邊緣修飾ZGNR的穩(wěn)定性最差。其中,N原子邊緣修飾3個(gè)原子寬度的ZGNR體系([N-ZG]3)形成能最小,Na原子邊緣修飾6個(gè)原子寬度的ZGNR體系([Na-ZG]6)形成能最大。第三,隨著納米帶寬度的增加(3~8),[B-ZG]體系的穩(wěn)定性稍微增加,但變化不大,而[N-ZG]體系的穩(wěn)定性逐漸降低。[Na-ZG]和[Mn-ZG]體系的形成能呈周期性振蕩特性。說(shuō)明納米帶寬度對(duì)金屬原子(Na,Mn)邊緣修飾ZGNR體系穩(wěn)定性的影響較大,而非金屬原子(B,N)邊緣修飾ZGNR體系的穩(wěn)定性受納米帶寬度影響較小。

    對(duì)于原子邊緣修飾ZGNR的穩(wěn)定結(jié)構(gòu),施加垂直于石墨烯納米帶平面的電場(chǎng),進(jìn)一步研究電場(chǎng)對(duì)結(jié)構(gòu)穩(wěn)定性的影響。下面以4個(gè)碳原子寬度的納米帶([M-ZG]4)為例, 計(jì)算了不同強(qiáng)度(0.5×1010 V/m,0.1×1011 V/m,0.15×1011 V/m)垂直電場(chǎng)作用下[M-ZG]4體系(分別記為[M-ZG-0.5]4,[M-ZG-1.0]4,[M-ZG-1.5]4)的形成能,結(jié)果如表2所示。由表2可知,隨著電場(chǎng)強(qiáng)度的增大,[B-ZG]4和[Na-ZG]4體系形成能越來(lái)越大,說(shuō)明施加電場(chǎng)不利于該體系的結(jié)構(gòu)穩(wěn)定性。對(duì)于[N-ZG]4體系,電場(chǎng)強(qiáng)度對(duì)其結(jié)構(gòu)穩(wěn)定性影響較小。而[Mn-ZG]4體系的形成能隨著電場(chǎng)值的增加而減小,穩(wěn)定性增加。此外,通過(guò)對(duì)[M-ZG]n體系中的鍵長(zhǎng)分析表明,施加電場(chǎng)前后的鍵長(zhǎng)變化也證明了上述結(jié)構(gòu)穩(wěn)定性的變化情況,具體鍵長(zhǎng)數(shù)據(jù)不再一一列舉。

    2.2 電子性質(zhì)

    下面通過(guò)電荷轉(zhuǎn)移、差分電荷密度、能帶結(jié)構(gòu)和態(tài)密度進(jìn)一步研究分析電場(chǎng)對(duì)邊緣修飾ZGNR電子性質(zhì)的影響。表3給出了[M-ZG]體系的電荷轉(zhuǎn)移情況,表4為施加不同電場(chǎng)強(qiáng)度下[M-ZG]4體系的電荷分析。由表3可知,首先,[B-ZG],[Na-ZG]和[Mn-ZG]體系的電荷轉(zhuǎn)移量大于零,表明電荷從B,Na和Mn原子轉(zhuǎn)移到近鄰的石墨烯C原子;[N-ZG]體系電荷轉(zhuǎn)移量小于零,說(shuō)明電荷從石墨烯C原子轉(zhuǎn)移到邊緣N原子。其次,對(duì)于4種不同的原子,電荷的轉(zhuǎn)移量基本上不受石墨烯納米帶寬度的影響。由表4可知,對(duì)[M-ZG]4體系施加垂直電場(chǎng)作用后,隨電場(chǎng)強(qiáng)度的增加,[B-ZG]4和[Mn-ZG]4體系電荷轉(zhuǎn)移量增大,[Na-ZG]4體系的電荷轉(zhuǎn)移量先減小后增大,而[N-ZG]4電荷轉(zhuǎn)移量絕對(duì)值減?。◤氖〤原子轉(zhuǎn)移到邊緣N原子的電荷減小)。綜合來(lái)說(shuō),電場(chǎng)促進(jìn)了邊緣修飾原子(B,Na,Mn)的電荷進(jìn)一步向石墨烯中的C原子轉(zhuǎn)移,阻礙了電荷從石墨烯C原子向邊緣N原子轉(zhuǎn)移。這些結(jié)果說(shuō)明電場(chǎng)有利于增加原子邊緣修飾ZGNR體系的導(dǎo)電性,對(duì)于調(diào)控石墨烯的電子結(jié)構(gòu)具有重要作用。

    為了直觀顯示[M-ZG]體系的電荷分布和轉(zhuǎn)移情況,圖2和圖3計(jì)算了[M-ZG]差分電荷密度分布情況。圖2為不同原子邊緣修飾5個(gè)原子寬度的ZGNR差分電荷密度分布,圖3為Mn原子邊緣修飾不同寬度ZGNR體系及不同電場(chǎng)下[Mn-ZG]4體系的差分電荷密度圖。由差分電荷密度分布可以看出,在電荷轉(zhuǎn)移的過(guò)程中,電荷在邊緣修飾M原子和石墨烯邊緣C原子之間聚集,形成了明顯的共價(jià)鍵,這種共價(jià)鍵有利于增強(qiáng)體系的穩(wěn)定性。相對(duì)于B,Na,Mn原子修飾的體系,在N原子修飾體系中,石墨烯中與N原子結(jié)合的C原子虧損更多的電荷。從圖3中可以看出電子從Mn原子附近轉(zhuǎn)移到ZGNR的C原子附近,并且隨著納米帶寬度的增加,電荷轉(zhuǎn)移量逐漸增多。施加垂直電場(chǎng)后,隨著電場(chǎng)值的增大,Mn原子附近虧損更多的電子,說(shuō)明電場(chǎng)值越大,Mn-C鍵的穩(wěn)定性越強(qiáng)。這些結(jié)論進(jìn)一步證明了表2和表3中電荷轉(zhuǎn)移的正確性。

    圖4~6分別給出了[Mn-ZG]體系及其在不同電場(chǎng)作用下的能帶圖、態(tài)密度(DOS)圖和分波態(tài)密度(PDOS)圖,圖5和圖6中正值(向上箭頭表示)和負(fù)值(向下箭頭表示)分別代表自旋向上和自旋向下態(tài)密度。從圖4和圖5對(duì)比分析可知,自旋向上和自旋向下的能帶是不重合的,表明[Mn-ZG]體系具有磁性。Mn原子邊緣修飾ZGNR在費(fèi)米面附近引入了新的能帶,并穿過(guò)費(fèi)米能級(jí),體系表現(xiàn)為典型的金屬能帶結(jié)構(gòu)。這種金屬性質(zhì)同樣也發(fā)生在B,N和Na原子修飾ZGNR的體系中。這與前人的相關(guān)研究是一致的,在以前的研究中,H原子邊緣修飾的ZGNR為金屬性[2];不同原子和官能團(tuán)(H,O,OH,NH,F(xiàn),CH3,NO2)邊緣修飾的ZGNR[15,17,26-29],主要表現(xiàn)為金屬性,但也會(huì)產(chǎn)生半導(dǎo)體性或者半金屬性。但是,對(duì)于Mn原子邊緣修飾AGNR[30],在反鐵磁性基態(tài)下表現(xiàn)出半導(dǎo)體性,在非磁性和鐵磁性基態(tài)下表現(xiàn)出金屬性,這與本文Mn原子邊緣修飾ZGNR始終表現(xiàn)為金屬性的研究結(jié)果是有所不同的。[Mn-ZG]體系的帶隙始終為0.0 eV,說(shuō)明其金屬特性不受修飾原子種類(lèi)、納米帶寬度以及外加電場(chǎng)的影響。在無(wú)電場(chǎng)作用情況下,電子結(jié)構(gòu)受納米帶寬度變化的影響較小。在外電場(chǎng)作用下,對(duì)于Mn邊緣修飾同一寬度的ZGNR,自旋向上的能帶逐漸向低能級(jí)處移動(dòng),而自旋向下的能帶逐漸向高能級(jí)處移動(dòng),最后都越過(guò)費(fèi)米能級(jí),增強(qiáng)了體系的金屬性和導(dǎo)電性。而且,隨著電場(chǎng)強(qiáng)度的增加,自旋簡(jiǎn)并程度降低。從圖6可進(jìn)一步看出,過(guò)渡金屬M(fèi)n原子的3d軌道與C原子的2p軌道發(fā)生雜化,其中,dxy和dyz軌道與ZGNR的σ鍵雜化,dxz和dz2軌道與ZGNR的π鍵雜化,在費(fèi)米能級(jí)附近形成穩(wěn)定化學(xué)鍵。這種軌道間的相互作用有利于邊緣修飾體系形成更穩(wěn)定的結(jié)構(gòu)。

    2.3 磁學(xué)性質(zhì)

    對(duì)于非金屬(B,N)和金屬原子(Na,Mn)邊緣修飾的ZGNR,表5、表6分別列出了[M-ZG]體系及其電場(chǎng)作用下的總磁矩。本征ZGNR是無(wú)磁性的,在原子(B,N,Na,Mn)邊緣修飾后,所有體系[M-ZG]均具有磁性,說(shuō)明邊緣修飾引起了ZGNR的磁化。其中,Mn原子修飾的納米帶[Mn-ZG]體系的磁矩最大,而Na原子修飾的納米帶[Na-ZG]體系的磁矩最小。隨著ZGNR寬度的增加,[N-ZG]和[Na-ZG]體系的磁性逐漸增強(qiáng),而[B-ZG]和[Mn-ZG]體系磁性變化較小。由表6可知隨著電場(chǎng)強(qiáng)度的增加,N和Na修飾體系的磁矩逐漸在減小,而B(niǎo)和Mn修飾體系磁矩逐漸增大。圖7給出了[B-ZG]體系和不同電場(chǎng)下[B-ZG]4體系的自旋電荷密度圖??梢钥闯觯跙-ZG]體系的磁性主要來(lái)源于非金屬B原子,在未加電場(chǎng)時(shí),石墨烯中的C原子沒(méi)有磁性,外加垂直電場(chǎng)使得GNR邊緣C原子出現(xiàn)自旋電荷,說(shuō)明ZGNR邊緣C原子被磁化。且隨電場(chǎng)強(qiáng)度增加,體系磁性增強(qiáng),說(shuō)明電場(chǎng)可以對(duì)邊緣修飾的體系磁性進(jìn)行調(diào)控,這為磁性儲(chǔ)存介質(zhì)及自旋電子學(xué)材料的應(yīng)用提供了理論基礎(chǔ)。

    3 結(jié) 論

    本文研究了原子(B,N,Na,Mn)邊緣修飾zigzag型石墨烯納米帶(ZGNR)的結(jié)構(gòu)穩(wěn)定性、電子性質(zhì)及磁學(xué)性質(zhì),以及外加垂直電場(chǎng)對(duì)原子邊緣修飾ZGNR體系電磁性質(zhì)的調(diào)控規(guī)律,得出以下結(jié)論:

    1)所有原子邊緣修飾石墨烯納米帶[M-ZG]n體系都能夠形成穩(wěn)定結(jié)構(gòu)。非金屬原子(B和N)邊緣修飾ZGNR比金屬原子(Na和Mn)邊緣修飾ZGNR更為穩(wěn)定。納米帶寬度對(duì)金屬原子邊緣修飾ZGNR體系穩(wěn)定性的影響較大,而對(duì)非金屬原子修飾體系的穩(wěn)定性影響較小。隨著電場(chǎng)強(qiáng)度的增大,[B-ZG]4和[Na-ZG]4體系穩(wěn)定性降低,[N-ZG]4體系結(jié)構(gòu)穩(wěn)定性基本不變,而[Mn-ZG]4體系穩(wěn)定性增加。

    2)[N-ZG]體系中電荷從石墨烯C原子轉(zhuǎn)移到邊緣N原子,而B(niǎo),Na,Mn修飾石墨烯體系的電荷轉(zhuǎn)移則反之,電荷的轉(zhuǎn)移量與納米帶寬度無(wú)關(guān)。電場(chǎng)阻礙了電荷從石墨烯C原子向邊緣N原子轉(zhuǎn)移,促進(jìn)了邊緣修飾原子(B,Na,Mn)的電荷進(jìn)一步向石墨烯中的C原子轉(zhuǎn)移。電荷在邊緣修飾原子M和石墨烯邊緣C原子之間形成了共價(jià)鍵。體系具有金屬特性,且不受修飾原子種類(lèi)、納米帶寬度的影響,外加電場(chǎng)增強(qiáng)了金屬性。

    3)在原子邊緣修飾后,所有[M-ZG]體系均具有磁性,其中,Mn原子修飾的納米帶[Mn-ZG]體系的磁矩最大,而Na原子修飾的納米帶[Na-ZG]體系的磁矩最小。[N-ZG]和[Na-ZG]體系的磁矩隨著納米帶寬度的增加而逐漸增強(qiáng)。隨著電場(chǎng)強(qiáng)度的增加,N和Na修飾體系的磁矩逐漸減小,而B(niǎo)和Mn修飾體系磁矩逐漸增大。電場(chǎng)可以對(duì)邊緣修飾石墨烯體系的磁性進(jìn)行有效調(diào)控。

    參考文獻(xiàn):

    [1]CAO L, LI X, JIA C, et al. Spin-charge transport properties for graphene/graphyne zigzag-edged nanoribbon heterojunctions: A first-principles study[J]. Carbon, 2018, 127:519-526.

    [2]SON Y W, COHEN M L, LOUIE S G.Energy gaps in graphene nanoribbons[J]. Physical Review Letters, 2006, 97(21):216803.

    [3]王曉偉,胡慧芳,張照錦,等.邊緣裁剪石墨烯納米帶的電子輸運(yùn)性質(zhì)研究[J]. 材料導(dǎo)報(bào), 2014, 28(5): 145-148.

    WANG X W, HU H F, ZHANG Z J, et al. Study on electronic and transport properties of the edge cutting graphene nanoribbons[J]. Materials Review, 2014, 28(5): 145-148.

    [4]周晨露, 鄭新亮, 鄭繼明, 等. 一種基于石墨烯納米帶的半金屬材料設(shè)計(jì)[J].西北大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 44(2): 201-205.

    ZHOU C L, ZHENG X L, ZHENG J M, et al. Designing a kind of half-metal material based on the sawtooth edge zigzag graphene nanoribbon[J].Journal of Northwest University(Natural Science Edition),2014,44(2):201-205.

    [5]孫凱剛, 解憂, 周安寧,等. 扶手椅型石墨烯納米帶吸附鈦原子鏈的電子結(jié)構(gòu)和磁性[J]. 陜西師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2016, 44(2):27-32.

    SUN K G, XIE Y, ZHOU A N, et al. Electronic structure and magnetic properties of zigzag Ti atomic chains adsorbed armchair graphene nanoribbons[J]. Journal of Shaanxi Normal University (Natural Science Edition), 2016, 44(2):27-32.

    [6]RUFFIEUX P, WANG S, YANG B, et al. On-surface synthesis of graphene nanoribbons with zigzag edge topology[J]. Nature, 2016, 531(7595):489-492.

    [7]KIMOUCHE A, ERVASTI M M, DROST R, et al. Ultra-narrow metallic armchair graphene nanoribbons[J]. Nature Communications, 2015, 6:10177.

    [8]CRUZ-SILVA E, BARNETT Z M, SUMPTER B G, et al. Structural, magnetic, and transport properties of substitutionally doped graphene nanoribbons from first principles[J]. Physical Review B, 2011, 83(15): 155445.

    [9]劉小月, 李林峰, 葛桂賢. 二茂錸分子吸附zigzag型石墨烯納米帶自旋輸運(yùn)性質(zhì)的理論研究 [J]. 固體電子學(xué)研究與進(jìn)展, 2019, 39(1):10-16.

    LIU X Y, LI L F, GE G X. The theoretical study of spin transport properties of zigzag graphene nanoribbons decorated with Re(cp)2 molecules[J]. Research & Progress of SSE, 2019, 39(1):10-16.

    [10] 張華林, 孫琳, 王鼎. 含單排線缺陷鋸齒型石墨烯納米帶的電磁性質(zhì)[J]. 物理學(xué)報(bào), 2016, 65(1): 016101.

    ZHANG H L, SUN L, WANG D. Electromagnetic properties of zigzag graphene nanoribbons with single-row line defect[J]. Acta Physica Sinica, 2016, 65 (1): 016101.

    [11]解憂, 張衛(wèi)濤, 曹松, 等. 過(guò)渡金屬原子鏈對(duì)雙層石墨烯納米帶的電磁性質(zhì)的調(diào)控[J]. 陜西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 46(6):54-60.

    XIE Y, ZHANG W T, CAO S, et al. The modulation of electronic and magnetic properties ofbilayer graphene nanoribbons by zigzag transition metal chains[J].Journal of Shaanxi Normal University (Natural Science Edition), 2018, 46(6):54-60.

    [12]YU X H, XU L, ZHANG J. Topological state engineering by in-plane electric field in graphene nanoribbon[J]. Physics Letters A, 2017, 381(34):2841-2844.

    [13]ZHANG W X, LIU Y X, TIAN H, et al. Electric field effect in ultrathin zigzag graphene nanoribbons[J]. Chinese Physics B, 2015, 24(7):076104.

    [14]WEI W, LIU R, DAN L, et al. Monte Carlo simulation of magnetic properties of a nano-graphene bilayer in a longitudinal magnetic field[J]. Superlattices and Microstructures, 2016, 98:458-472.

    [15]GUNLYCKE D, LI J, MINTMIRE J W, et al. Altering low-bias transport in zigzag-edge graphene nanostrips with edge chemistry[J]. Applied Physics Letters, 2007, 91(11):19912.

    [16]MAO Y, HAO W, WEI X, et al. Edge-adsorption of potassium adatoms on graphene nanoribbon: A first principle study[J]. Applied Surface Science, 2013, 280:698-704.

    [17]LU D B, SONG Y L, YANG Z X, et al. Energy gap modulation of graphene nanoribbons by F termination[J]. Applied Surface Science, 2011, 257(15):6440-6444.

    [18]孫宇, 宋新祥, 汪潔生, 等. 羥基修飾石墨烯納米帶吸附鈉和氯原子的影響[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 47(S1):143-149.

    SUN Y, SONG X X, WANG J S, et al. Effects of sodium and chlorine adatoms on graphene nanoribbons with hydroxyl terminated edge[J]. Periodical of Ocean University of China(Natural Sciences Edition), 2017,47(S1):143-149.

    [19]WI J H, RIA N C, RIA S I, et al. Mechanical and electronic properties of armchair graphene nanoribbons with symmetrically double-lines-doped BN under uniaxial tensile strain: Ab initio study[J]. Solid State Communications, 2019, 300: 113644.

    [20]朱朕, 李春先, 張振華. 功能化扶手椅型石墨烯納米帶異質(zhì)結(jié)的磁器件特性[J].物理學(xué)報(bào), 2016, 65(11): 118501.

    ZHU Z, LI C X, ZHANG Z H. Magnetic device properties for a heterojunction based on functionalized armchair-edged graphene nanoribbons[J]. Acta Physica Sinica, 2016, 65(11): 118501.

    [21]KRESSE G, FURTHMLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set[J]. Physical Review B, 1996, 54(16): 11169.

    [22]KRESSE G, FURTHMLLER J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set[J]. Computational Materials Science, 1996, 6(1): 15.

    [23]KRESSE G, JOUBERT D. From ultrasoft pseudopotentials to the projector augmented-wave method[J]. Physical Review B, 1999, 59(3): 1758.

    [24]PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple[J]. Physical Review Letters, 1996, 77(18): 3865.

    [25]MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations[J]. Physical Review B, 1976, 13(12): 5188-5192.

    [26]WU M H, WU X J, ZENG X C. Exploration of half metallicity in edge-modified graphene nanoribbons[J]. Journal of Physical Chemistry C, 2010, 114(9), 3937-3944.

    [27]蓋永杰. 邊緣化學(xué)修飾鋸齒型石墨烯納米帶的電子結(jié)構(gòu)[J].長(zhǎng)春大學(xué)學(xué)報(bào),2011, 21(12):70-72.

    GAI Y J. Electronic structure of edge-chemically modified zigzag graphene nanoribbons[J]. Journal of Changchun University, 2011, 21(12): 70-72.

    [28]江帆, 何珊.邊緣修飾的石墨烯納米帶電子特性的第一性原理研究[J].原子與分子物理學(xué)報(bào), 2011, 28(6): 1137-1142.

    JIANG F, HE S. First-principle calculations of electronic properties for edge-modified graphene nanoribbons[J]. Journal of Atomic and Molecular Physics, 2011, 28(6): 1137-1142.

    [29]歐陽(yáng)方平, 徐慧, 李明君, 等. Armchair型石墨納米帶的電子結(jié)構(gòu)和輸運(yùn)性質(zhì)[J]. 物理化學(xué)學(xué)報(bào), 2008, 24(2): 328-332.

    OUYANG F P, XU H, LI M J, et al. Electronic structure and transport properties of armchair graphene nanoribbons[J]. Acta Physico-Chimica Sinica, 2018, 24(2):328-332.

    [30]ZHU Z, WANG D, ZHANG Z H, et al. Magnetic structures and magnetic device properties of edge modified armchair-edged graphene nanoribbons[J]. Carbon, 2016, 106: 252-259.

    (編 輯 李 靜)

    收稿日期:2019-08-25

    基金項(xiàng)目:中國(guó)博士后科學(xué)基金資助項(xiàng)目(2014M560798);陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃(2013JM8004)

    作者簡(jiǎn)介:吳秀,女,陜西榆林人,從事納米電子學(xué)研究。

    通信作者:解憂,男,安徽宿州人,博士,教授,從事凝聚態(tài)物理研究。

    美女高潮喷水抽搐中文字幕| 好男人在线观看高清免费视频 | 国产亚洲欧美在线一区二区| 欧美精品啪啪一区二区三区| 亚洲一区高清亚洲精品| АⅤ资源中文在线天堂| 午夜福利高清视频| 高潮久久久久久久久久久不卡| 亚洲成av人片免费观看| 哪里可以看免费的av片| 亚洲 国产 在线| 久久久久国产精品人妻aⅴ院| 亚洲第一欧美日韩一区二区三区| 亚洲片人在线观看| 精品久久久久久成人av| 久久香蕉国产精品| 亚洲精品av麻豆狂野| 国产三级黄色录像| 国产激情欧美一区二区| 亚洲,欧美精品.| 久久人妻av系列| 久久中文看片网| 国产成人av激情在线播放| 国产精品一区二区三区四区久久 | 欧美人与性动交α欧美精品济南到| 久久久久久久午夜电影| 国产精品久久久久久亚洲av鲁大| 欧美久久黑人一区二区| 18禁裸乳无遮挡免费网站照片 | 一区二区三区高清视频在线| 日本一本二区三区精品| 久久精品国产亚洲av高清一级| 国产91精品成人一区二区三区| 精品国产亚洲在线| 国产亚洲精品综合一区在线观看 | 好男人在线观看高清免费视频 | 亚洲男人的天堂狠狠| 青草久久国产| 久久久久久大精品| 88av欧美| 2021天堂中文幕一二区在线观 | 免费在线观看完整版高清| 2021天堂中文幕一二区在线观 | 亚洲精品色激情综合| 欧美不卡视频在线免费观看 | 亚洲中文日韩欧美视频| 啦啦啦免费观看视频1| 白带黄色成豆腐渣| 天堂影院成人在线观看| 无人区码免费观看不卡| 国产高清有码在线观看视频 | 桃红色精品国产亚洲av| 日本免费一区二区三区高清不卡| 国产aⅴ精品一区二区三区波| 制服诱惑二区| 亚洲狠狠婷婷综合久久图片| 亚洲电影在线观看av| 日本一本二区三区精品| 国产成人欧美在线观看| 午夜激情福利司机影院| 桃红色精品国产亚洲av| 欧美性猛交╳xxx乱大交人| 久久精品国产清高在天天线| 女人高潮潮喷娇喘18禁视频| 黑人操中国人逼视频| 后天国语完整版免费观看| www.999成人在线观看| 嫩草影视91久久| 少妇 在线观看| 精品日产1卡2卡| 亚洲欧美一区二区三区黑人| 国产av不卡久久| 在线永久观看黄色视频| 日韩欧美 国产精品| 久久香蕉精品热| 久久国产亚洲av麻豆专区| 99久久综合精品五月天人人| 大型av网站在线播放| 黄色视频不卡| 久久天堂一区二区三区四区| 日本五十路高清| 亚洲成人久久性| 老汉色∧v一级毛片| 男人舔女人下体高潮全视频| 老汉色av国产亚洲站长工具| 搡老妇女老女人老熟妇| 国产精品av久久久久免费| 亚洲狠狠婷婷综合久久图片| 国产一卡二卡三卡精品| 国产精品香港三级国产av潘金莲| 精品电影一区二区在线| 亚洲专区中文字幕在线| 12—13女人毛片做爰片一| 久久国产精品影院| 伦理电影免费视频| 女性生殖器流出的白浆| 日本 av在线| 香蕉av资源在线| 日本a在线网址| 一进一出好大好爽视频| 成人午夜高清在线视频 | 窝窝影院91人妻| 欧美三级亚洲精品| 真人做人爱边吃奶动态| 女人爽到高潮嗷嗷叫在线视频| 99久久久亚洲精品蜜臀av| 亚洲成人精品中文字幕电影| 不卡一级毛片| 日韩精品中文字幕看吧| 久久精品国产综合久久久| 国产精品一区二区免费欧美| 在线观看免费日韩欧美大片| 18美女黄网站色大片免费观看| 日韩视频一区二区在线观看| 极品教师在线免费播放| 亚洲精品久久成人aⅴ小说| 亚洲全国av大片| 女人高潮潮喷娇喘18禁视频| 成人特级黄色片久久久久久久| 成人永久免费在线观看视频| 熟妇人妻久久中文字幕3abv| 亚洲午夜理论影院| 一区二区三区国产精品乱码| 两个人看的免费小视频| 欧美一级a爱片免费观看看 | 丝袜人妻中文字幕| 免费人成视频x8x8入口观看| 99精品欧美一区二区三区四区| 男女下面进入的视频免费午夜 | 国产精品二区激情视频| 久久精品aⅴ一区二区三区四区| 日韩欧美 国产精品| 国产午夜福利久久久久久| 侵犯人妻中文字幕一二三四区| 国产亚洲精品久久久久久毛片| 精品国产亚洲在线| 黄色 视频免费看| 亚洲精品在线美女| 亚洲精品中文字幕在线视频| 999久久久国产精品视频| 欧美日韩精品网址| 色老头精品视频在线观看| 亚洲熟女毛片儿| 精品国产亚洲在线| 熟女电影av网| 久久久久久九九精品二区国产 | 国产伦在线观看视频一区| 夜夜躁狠狠躁天天躁| 国产精品永久免费网站| 欧美中文综合在线视频| 老司机午夜十八禁免费视频| 国产精品美女特级片免费视频播放器 | 日韩欧美免费精品| 99在线人妻在线中文字幕| 国产av不卡久久| 黄片小视频在线播放| 丁香欧美五月| 国产不卡一卡二| 此物有八面人人有两片| 一区二区三区激情视频| 欧美 亚洲 国产 日韩一| 中文字幕高清在线视频| 一级毛片精品| 国产成人精品久久二区二区免费| 可以在线观看毛片的网站| 国产高清videossex| 在线看三级毛片| 亚洲精品久久国产高清桃花| 欧美成狂野欧美在线观看| a级毛片在线看网站| 啪啪无遮挡十八禁网站| 免费搜索国产男女视频| 国产欧美日韩一区二区精品| 在线观看日韩欧美| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧洲精品一区二区精品久久久| 黑人欧美特级aaaaaa片| 一本精品99久久精品77| 成人国语在线视频| 午夜激情福利司机影院| 久久久国产精品麻豆| 久久久久久久久久黄片| 国产成年人精品一区二区| 露出奶头的视频| 91av网站免费观看| 一a级毛片在线观看| 欧美绝顶高潮抽搐喷水| 国产伦一二天堂av在线观看| 美女国产高潮福利片在线看| 两性午夜刺激爽爽歪歪视频在线观看 | 无遮挡黄片免费观看| 黄色女人牲交| 国产欧美日韩精品亚洲av| 黄色视频,在线免费观看| 成人免费观看视频高清| 成在线人永久免费视频| 国产一区二区三区在线臀色熟女| 国产精品自产拍在线观看55亚洲| 中文字幕另类日韩欧美亚洲嫩草| 久久中文字幕人妻熟女| 99久久无色码亚洲精品果冻| 俺也久久电影网| 久久精品国产亚洲av高清一级| 国产av又大| 在线观看免费日韩欧美大片| 99久久国产精品久久久| aaaaa片日本免费| xxxwww97欧美| 久久伊人香网站| 日韩精品中文字幕看吧| 精品少妇一区二区三区视频日本电影| 久久欧美精品欧美久久欧美| 在线免费观看的www视频| 久久中文看片网| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲成人久久爱视频| 99国产综合亚洲精品| 变态另类丝袜制服| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美一区二区三区黑人| 国产精品精品国产色婷婷| 巨乳人妻的诱惑在线观看| 国产欧美日韩一区二区三| 亚洲av日韩精品久久久久久密| 中文字幕另类日韩欧美亚洲嫩草| 51午夜福利影视在线观看| 欧美乱妇无乱码| 香蕉国产在线看| 国产成人一区二区三区免费视频网站| 国产真实乱freesex| 又大又爽又粗| 免费在线观看成人毛片| 久久精品国产99精品国产亚洲性色| 日韩欧美在线二视频| 一卡2卡三卡四卡精品乱码亚洲| 免费在线观看日本一区| 人妻丰满熟妇av一区二区三区| 男男h啪啪无遮挡| 国产激情偷乱视频一区二区| 成人18禁在线播放| 99国产精品99久久久久| 精品国内亚洲2022精品成人| 欧美中文综合在线视频| 亚洲国产中文字幕在线视频| 亚洲五月天丁香| 午夜免费激情av| 久久性视频一级片| 99精品久久久久人妻精品| 亚洲成av人片免费观看| 午夜激情福利司机影院| 国产aⅴ精品一区二区三区波| 久久久久久久久免费视频了| 色播在线永久视频| 久久久久久免费高清国产稀缺| 亚洲精品国产区一区二| 黄色 视频免费看| 日韩成人在线观看一区二区三区| 男女午夜视频在线观看| 国产主播在线观看一区二区| 亚洲精品粉嫩美女一区| 午夜福利免费观看在线| 成在线人永久免费视频| 成人亚洲精品一区在线观看| 午夜亚洲福利在线播放| 成在线人永久免费视频| 麻豆av在线久日| 亚洲人成电影免费在线| 亚洲五月婷婷丁香| 欧美另类亚洲清纯唯美| 99精品在免费线老司机午夜| 中文字幕高清在线视频| 9191精品国产免费久久| 久久香蕉精品热| 老司机靠b影院| 日韩欧美国产一区二区入口| 在线永久观看黄色视频| 亚洲成av片中文字幕在线观看| 亚洲男人天堂网一区| 女生性感内裤真人,穿戴方法视频| 精品国产亚洲在线| 老司机靠b影院| 好男人电影高清在线观看| 99re在线观看精品视频| 99热只有精品国产| 亚洲精品美女久久av网站| 动漫黄色视频在线观看| 国产亚洲精品一区二区www| 欧美日韩精品网址| 国产野战对白在线观看| 1024香蕉在线观看| 亚洲男人天堂网一区| 啪啪无遮挡十八禁网站| tocl精华| 婷婷丁香在线五月| 高清在线国产一区| 黄片播放在线免费| 成人18禁高潮啪啪吃奶动态图| 一区二区三区国产精品乱码| 一本大道久久a久久精品| 亚洲专区国产一区二区| 久久欧美精品欧美久久欧美| 搡老岳熟女国产| 19禁男女啪啪无遮挡网站| cao死你这个sao货| tocl精华| 三级毛片av免费| 成年女人毛片免费观看观看9| 黄色毛片三级朝国网站| 久久国产精品影院| 午夜视频精品福利| 国产男靠女视频免费网站| 中文字幕人成人乱码亚洲影| 国产亚洲欧美在线一区二区| 中文字幕另类日韩欧美亚洲嫩草| 18美女黄网站色大片免费观看| 夜夜夜夜夜久久久久| 色尼玛亚洲综合影院| 757午夜福利合集在线观看| 亚洲国产欧美日韩在线播放| 久久久久亚洲av毛片大全| 亚洲精品久久成人aⅴ小说| 亚洲狠狠婷婷综合久久图片| av超薄肉色丝袜交足视频| 精品欧美国产一区二区三| 欧美色欧美亚洲另类二区| 精品第一国产精品| 久久国产精品影院| 国产男靠女视频免费网站| 19禁男女啪啪无遮挡网站| 亚洲成av人片免费观看| 欧美中文日本在线观看视频| 日本在线视频免费播放| 99国产精品99久久久久| 国语自产精品视频在线第100页| 2021天堂中文幕一二区在线观 | 日韩 欧美 亚洲 中文字幕| 午夜免费观看网址| 免费高清在线观看日韩| 黄片大片在线免费观看| 久久精品国产99精品国产亚洲性色| 老司机午夜福利在线观看视频| 免费观看精品视频网站| www日本在线高清视频| 男女午夜视频在线观看| 国产蜜桃级精品一区二区三区| 久久久久久国产a免费观看| 亚洲成av人片免费观看| 亚洲精华国产精华精| 欧美日韩一级在线毛片| 久久草成人影院| 久久久久久久久免费视频了| 一级毛片女人18水好多| 岛国在线观看网站| 久久 成人 亚洲| 天天一区二区日本电影三级| 97碰自拍视频| 国产一区二区三区在线臀色熟女| 俺也久久电影网| 中文字幕精品亚洲无线码一区 | 老司机午夜十八禁免费视频| 欧美大码av| 午夜两性在线视频| 在线永久观看黄色视频| 18美女黄网站色大片免费观看| 欧美成人免费av一区二区三区| 日韩国内少妇激情av| 欧美又色又爽又黄视频| 91麻豆精品激情在线观看国产| 欧美激情高清一区二区三区| 国产真实乱freesex| www.精华液| 精品国产国语对白av| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 免费女性裸体啪啪无遮挡网站| 免费看日本二区| 国产精品野战在线观看| 黑人巨大精品欧美一区二区mp4| av在线播放免费不卡| 免费搜索国产男女视频| 午夜精品久久久久久毛片777| 国产人伦9x9x在线观看| 一本综合久久免费| 黑人巨大精品欧美一区二区mp4| 波多野结衣av一区二区av| 性色av乱码一区二区三区2| 国产成人av教育| 国产在线精品亚洲第一网站| 白带黄色成豆腐渣| 成人av一区二区三区在线看| 午夜福利欧美成人| 天天添夜夜摸| 三级毛片av免费| 国产精品野战在线观看| 婷婷精品国产亚洲av在线| 老汉色av国产亚洲站长工具| 精品欧美国产一区二区三| 久久人妻av系列| 少妇粗大呻吟视频| 国产av不卡久久| 美女国产高潮福利片在线看| 国产精品久久久av美女十八| videosex国产| 成年人黄色毛片网站| 久久热在线av| 999精品在线视频| 美女高潮喷水抽搐中文字幕| 精品欧美国产一区二区三| 看片在线看免费视频| 老司机靠b影院| 久久精品国产亚洲av香蕉五月| 香蕉丝袜av| 成年女人毛片免费观看观看9| 亚洲av五月六月丁香网| 国产成人系列免费观看| 亚洲欧美一区二区三区黑人| 很黄的视频免费| 淫妇啪啪啪对白视频| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 久久精品91无色码中文字幕| 久久久国产成人精品二区| 老司机在亚洲福利影院| 美女大奶头视频| 听说在线观看完整版免费高清| 欧美色视频一区免费| 男男h啪啪无遮挡| 每晚都被弄得嗷嗷叫到高潮| 国产成年人精品一区二区| 真人一进一出gif抽搐免费| 可以在线观看毛片的网站| 91成人精品电影| 免费观看人在逋| 中出人妻视频一区二区| 亚洲成人久久爱视频| 母亲3免费完整高清在线观看| 老司机午夜福利在线观看视频| 久久精品亚洲精品国产色婷小说| 美女免费视频网站| 欧美在线一区亚洲| 久久人妻av系列| 中文在线观看免费www的网站 | √禁漫天堂资源中文www| 精品人妻1区二区| 十八禁人妻一区二区| 怎么达到女性高潮| 麻豆久久精品国产亚洲av| 99国产综合亚洲精品| 国产视频内射| 欧美 亚洲 国产 日韩一| 国产91精品成人一区二区三区| 久久亚洲真实| 国内久久婷婷六月综合欲色啪| av免费在线观看网站| 久久久水蜜桃国产精品网| 两性夫妻黄色片| 国内精品久久久久精免费| 在线观看www视频免费| 久久久精品欧美日韩精品| 国产97色在线日韩免费| 两性午夜刺激爽爽歪歪视频在线观看 | 国产欧美日韩精品亚洲av| 一本一本综合久久| 久久 成人 亚洲| 久久精品国产亚洲av香蕉五月| 一本久久中文字幕| 免费一级毛片在线播放高清视频| 久久精品夜夜夜夜夜久久蜜豆 | 日韩欧美在线二视频| 精品国产超薄肉色丝袜足j| 亚洲熟妇熟女久久| 无人区码免费观看不卡| 精品高清国产在线一区| 国产高清激情床上av| 久久久久九九精品影院| 脱女人内裤的视频| 国产三级在线视频| 欧美zozozo另类| 久久热在线av| 一区二区三区高清视频在线| cao死你这个sao货| 欧美大码av| 精品电影一区二区在线| 中亚洲国语对白在线视频| 看免费av毛片| 神马国产精品三级电影在线观看 | www.自偷自拍.com| 91av网站免费观看| 欧美久久黑人一区二区| 日韩欧美 国产精品| a级毛片a级免费在线| 亚洲成国产人片在线观看| 精品少妇一区二区三区视频日本电影| 久久中文看片网| 国产成人av教育| 国产高清视频在线播放一区| 亚洲七黄色美女视频| 免费一级毛片在线播放高清视频| 国产av不卡久久| 成人18禁高潮啪啪吃奶动态图| 国产又色又爽无遮挡免费看| 国产真实乱freesex| 欧美不卡视频在线免费观看 | 成在线人永久免费视频| 99久久无色码亚洲精品果冻| 看免费av毛片| av欧美777| 国产成人欧美在线观看| 妹子高潮喷水视频| 大型av网站在线播放| 精品久久久久久久人妻蜜臀av| 18禁国产床啪视频网站| 欧美中文日本在线观看视频| 一区二区三区高清视频在线| 99久久无色码亚洲精品果冻| www.熟女人妻精品国产| 成人av一区二区三区在线看| av在线天堂中文字幕| 亚洲第一欧美日韩一区二区三区| 久久香蕉精品热| 欧美黑人巨大hd| 欧美乱色亚洲激情| 在线天堂中文资源库| 黄色丝袜av网址大全| 精华霜和精华液先用哪个| 精品一区二区三区四区五区乱码| 黄网站色视频无遮挡免费观看| 国产成人精品久久二区二区91| 午夜日韩欧美国产| 亚洲成av片中文字幕在线观看| 两人在一起打扑克的视频| 亚洲精品中文字幕在线视频| 色精品久久人妻99蜜桃| 日韩精品中文字幕看吧| av在线天堂中文字幕| 可以在线观看毛片的网站| 久久中文看片网| 国产久久久一区二区三区| 一级a爱片免费观看的视频| 非洲黑人性xxxx精品又粗又长| 欧美日韩黄片免| 在线观看午夜福利视频| 岛国在线观看网站| 久久人妻福利社区极品人妻图片| 99久久99久久久精品蜜桃| 老汉色av国产亚洲站长工具| 色尼玛亚洲综合影院| 一级a爱片免费观看的视频| xxxwww97欧美| 神马国产精品三级电影在线观看 | 欧美性猛交黑人性爽| 人妻丰满熟妇av一区二区三区| 两个人免费观看高清视频| 久久久久久亚洲精品国产蜜桃av| 亚洲专区中文字幕在线| tocl精华| 亚洲电影在线观看av| 色精品久久人妻99蜜桃| 亚洲色图av天堂| 夜夜爽天天搞| 99国产精品一区二区三区| 免费高清在线观看日韩| 波多野结衣巨乳人妻| 美女高潮到喷水免费观看| 两个人看的免费小视频| 他把我摸到了高潮在线观看| 亚洲av电影不卡..在线观看| 精品久久久久久,| 欧美黄色淫秽网站| 国产精品影院久久| 亚洲人成伊人成综合网2020| 久久婷婷成人综合色麻豆| 两性夫妻黄色片| 搡老妇女老女人老熟妇| 亚洲电影在线观看av| 亚洲精品中文字幕一二三四区| 精品国产乱子伦一区二区三区| 在线十欧美十亚洲十日本专区| 天堂动漫精品| 久久久久国产一级毛片高清牌| 好男人电影高清在线观看| 国内毛片毛片毛片毛片毛片| 久久久久国产一级毛片高清牌| 亚洲自拍偷在线| 在线观看午夜福利视频| 午夜久久久久精精品| 亚洲美女黄片视频| 美女 人体艺术 gogo| 国产视频一区二区在线看| 69av精品久久久久久| 色综合婷婷激情| 88av欧美| 亚洲在线自拍视频| 国产精品香港三级国产av潘金莲| 国产真实乱freesex| 久久婷婷成人综合色麻豆| 满18在线观看网站| 国产区一区二久久| 亚洲成人免费电影在线观看| 丝袜人妻中文字幕| 国产激情欧美一区二区| 久久香蕉国产精品| 国产1区2区3区精品| 久久久久久国产a免费观看| 最新美女视频免费是黄的| 久久精品国产清高在天天线| 日韩欧美国产在线观看| 亚洲国产精品sss在线观看| 国产视频内射| 午夜影院日韩av| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 成人手机av| 老司机靠b影院| 亚洲成人久久爱视频|