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

    石墨烯條帶的電子結(jié)構(gòu)與性質(zhì):電場及長度效應(yīng)

    2014-10-18 05:27:52梁萬珍
    物理化學(xué)學(xué)報(bào) 2014年3期
    關(guān)鍵詞:電子密度基態(tài)條帶

    孫 進(jìn) 梁萬珍

    (1安徽大學(xué)物理與材料科學(xué)學(xué)院,合肥 230039;2廈門大學(xué)化學(xué)系,福建廈門 361005)

    1 引言

    自從單層的石墨也就是石墨烯在實(shí)驗(yàn)上被制備出來以后,1-3人們不斷發(fā)現(xiàn)這是一種具有有趣的物理特性以及顯著應(yīng)用前景的納米材料.目前,石墨烯以及以石墨烯為基礎(chǔ)的其它結(jié)構(gòu)的材料已引起人們廣泛的關(guān)注和研究,例如,s-p磁性的早期發(fā)現(xiàn),納米尺度石墨烯在技術(shù)上的開發(fā),以及表面修飾、外場、壓力對于石墨烯性質(zhì)的影響等等.4-13石墨烯的可控性為人們制備出高效率、低成本的電磁分子器件提供了一條可能的途徑.

    如果將單層的石墨烯沿著一些特殊的線剪開,就可以得到具有準(zhǔn)一維結(jié)構(gòu)的石墨烯納米帶(GNR).石墨烯納米帶材料與碳納米管有著天然的聯(lián)系.沿著周期性方向,按照帶的兩邊不同的形狀,石墨烯納米帶也可以分為椅子型(armchair)和鋸齒型(zigzag)兩種不同的類型.不同的邊帶形狀或不同的長度都可能會影響石墨烯納米帶的電子和磁性質(zhì).我們可以通過研究納米帶材料來增加對納米管的認(rèn)識,同時(shí),在納米管的研究中,一些已經(jīng)比較完備的實(shí)驗(yàn)手段也可以在納米帶材料上得到充分應(yīng)用.最近,這種在10 nm以下尺寸的石墨烯納米帶也的確可以由化學(xué)方法獲得.14相應(yīng)的在這個(gè)尺度上納米帶的物理化學(xué)性質(zhì)研究必將會有更大的發(fā)展.

    加州大學(xué)伯克利分校的Louie研究組15利用密度泛函理論和緊束縛近似研究了armchair型和zigzag型納米帶材料的電子結(jié)構(gòu).研究顯示,兩種結(jié)構(gòu)的石墨烯納米帶都是半導(dǎo)體性質(zhì).而zigzag型由于邊界態(tài)在費(fèi)米能級附近,會通過自旋極化來降低體系的能量,導(dǎo)致磁性的產(chǎn)生,引起了更廣泛的研究,包括它的自旋性質(zhì)、磁性、半金屬性等等.16-21

    由于石墨烯特殊的電導(dǎo)率,22它的光學(xué)及激發(fā)態(tài)性質(zhì)一直是人們關(guān)注的一個(gè)焦點(diǎn).理論上,通過對激發(fā)態(tài)上光學(xué)性質(zhì)的研究,單、雙層以及其它類型的石墨烯的光學(xué)電導(dǎo)率相繼被計(jì)算出來,23-27從而有效地促進(jìn)了石墨烯材料作為分子器件的應(yīng)用與發(fā)展.我們知道具有n個(gè)重復(fù)單元的zigzag型石墨烯納米帶(n-ZGNR)的基態(tài)具有半導(dǎo)體性質(zhì),如果它的周期性重復(fù)單元的個(gè)數(shù)n<32的話,那么它的反鐵磁態(tài)的能量會總是小于鐵磁態(tài)的能量.因而在這種情況下,基態(tài)上n-ZGNR兩個(gè)邊上的電子自旋方向相反,具有反鐵磁性.17,20但是,目前還很少有文章報(bào)道過石墨烯納米帶在激發(fā)態(tài)上,不同自旋方向電子的不同之處,以及它們對于激發(fā)態(tài)和光學(xué)等性質(zhì)的影響.Lu等28曾經(jīng)對較短的石墨烯納米帶進(jìn)行過計(jì)算,他們利用極化子模型,比較了不同自旋的自旋波函數(shù),并且研究了不同自旋的電荷響應(yīng)動力學(xué)過程.此外,Zhao等29則在平均場近似下,研究了以石墨烯作為分子器件在輸運(yùn)過程中不同的自旋電流.

    本文我們選擇了體系較小的zigzag型石墨烯條帶分子,利用密度泛函理論(DFT)計(jì)算了其基態(tài)的電子結(jié)構(gòu),并在此基礎(chǔ)上采用實(shí)時(shí)的含時(shí)密度泛函理論(TDDFT)方法,對超短激光脈沖影響下實(shí)時(shí)的電子動力學(xué)性質(zhì)進(jìn)行了研究.避免線性展開,求解完整的含時(shí)密度泛函方程或含時(shí)Hartree-Fock方程的方法,可以直接用來研究在不同的外激光脈沖影響下的電子動力學(xué)過程以及多光子激發(fā).因而這種方法目前被廣泛地用于強(qiáng)場過程的研究.30-33對于石墨烯條帶分子,深入研究它們在超短激光脈沖影響下的電子動力學(xué)過程有助于其在分子器件方面的應(yīng)用與發(fā)展,因而是必須和迫切的.

    本文采用密度泛函及含時(shí)密度泛函方法討論了寬度上包含8個(gè)zigzag鏈的石墨烯條帶基態(tài)及激發(fā)態(tài)的性質(zhì),包括條帶長度、靜電場對于石墨烯條帶磁性的影響以及激發(fā)態(tài)上隨著激光脈沖的變化,自旋向上和自旋向下的電子的不同表現(xiàn).

    2 計(jì)算方法

    在含時(shí)密度泛函理論框架下,第i個(gè)分子軌道波函數(shù)Ψi滿足含時(shí)薛定諤方程

    其中υext(r)和υxc[]ρ(r)分別是外勢和交換相關(guān)勢.等式(1)可以寫成密度矩陣形式的含時(shí)密度泛函方程,也就是標(biāo)準(zhǔn)的Liuviolle方程

    其中F(t)和ρ(t)分別為Fock矩陣和密度矩陣.對于等式(3)的求解,可以在微擾理論基礎(chǔ)上分別對Fock矩陣和密度矩陣進(jìn)行展開,從而得到線性和非線性的性質(zhì).

    目前,含時(shí)密度泛函理論已經(jīng)成功地被應(yīng)用于線性光學(xué)性質(zhì)的計(jì)算.但是對于一些超短的強(qiáng)激光脈沖過程,非線性性質(zhì)非常顯著,如果采用線性響應(yīng)的方法則無法體現(xiàn)出來.為了充分考慮系統(tǒng)與外電場之間的非線性相互作用,本文中,我們避免使用微擾展開理論,而是求解完整的含時(shí)密度泛函方程,從而使得在求解的過程中能夠平等地對待外電場和原子或分子之間的庫侖力.含時(shí)密度泛函方程和含時(shí)Hartree-Fock方程一樣是一個(gè)微分方程,因而可以用一些標(biāo)準(zhǔn)的數(shù)值計(jì)算方法進(jìn)行求解.但是低階的Verlet和Runge-Kutta數(shù)值方法34都需要非常小的時(shí)間步長來確保密度矩陣的厄米性.除此以外,在每一個(gè)時(shí)間步長里,Fock矩陣都要被構(gòu)造很多次,這是非常耗時(shí)的.為了提高整個(gè)計(jì)算的效率,我們充分利用厄米變換,并采用改進(jìn)的中間點(diǎn)計(jì)算方法.35考慮到隨著時(shí)間步長的推進(jìn),密度和外場的線性變化,通過對時(shí)間步長中間點(diǎn)的Fock矩陣的計(jì)算來演化密度矩陣.這種方法能夠非常有效地處理求解含時(shí)的密度泛函方程過程中對于時(shí)間演化算符的積分.我們已經(jīng)用此方法成功計(jì)算過碳納米管等納米材料的光學(xué)及動力學(xué)性質(zhì).36-39具體來看,在含時(shí)密度泛函理論中,時(shí)間演化過程可概括如下:

    在此過程中,首先利用ti時(shí)間點(diǎn)的Fock矩陣F(ti)按照公式(4,5)將密度矩陣從ρ(ti-1/2)演化到ρ(ti+1/2),然后用同樣的方法按照公式(6,7)可以得到ti+1時(shí)間點(diǎn)的密度矩陣ρ(ti+1)和Fock矩陣F(ti+1),從而為下一個(gè)時(shí)間步長的計(jì)算做準(zhǔn)備.由于在演化過程中,矩陣U和U'都是厄米矩陣,在每一個(gè)時(shí)間步長中,矩陣的等冪性都能夠得到保證.此外,由于在每一個(gè)時(shí)間步長中,只需要構(gòu)造一次Fock矩陣,因而計(jì)算速度得到了很大的提高.對于演化過程中矩陣的冪指數(shù)計(jì)算,我們用Taylor展開進(jìn)行處理.

    在過去的研究中我們發(fā)現(xiàn),當(dāng)展開階數(shù)k=100時(shí),計(jì)算結(jié)果就能夠保持穩(wěn)定.36在研究過程中,由于希望能夠觀察到石墨烯條帶分子在激發(fā)態(tài)上的動力學(xué)過程,因而加入一個(gè)極化的強(qiáng)激光脈沖,

    其中F、ω、τ和t分別指外激光場的強(qiáng)度、頻率、脈沖持續(xù)寬度及時(shí)間.

    3 計(jì)算結(jié)果及討論

    本文選擇了長度較短的,僅僅包含了8個(gè)zigzag鏈(8-ZGNR)的zigzag型石墨烯納米條帶,其化學(xué)式為C168H36.由于長度較短,C168H36在有些方面表現(xiàn)得更像一個(gè)分子,但是作為zigzag型石墨烯納米條帶的單元分子,也能夠反映出納米條帶的部分性質(zhì).我們在DFT/TDDFT的理論框架下,采用B3LYP雜化泛函,對C168H36的基態(tài)和激發(fā)態(tài)的性質(zhì)進(jìn)行研究.由于要進(jìn)行實(shí)時(shí)的TDDFT計(jì)算,考慮到計(jì)算量的問題,采用小基組STO-3G,對應(yīng)要處理876個(gè)軌道基函數(shù).這一基組在我們之前的研究工作40中也采用過,并且得到了較好的結(jié)果.對于C168H36分子的結(jié)構(gòu)在UB3LYP/STO-3G的量子化學(xué)理論水平上進(jìn)行了優(yōu)化,得到的分子結(jié)構(gòu)如圖1所示.分子在同一平面上,C―C鍵的長度為0.144 nm左右,C―H鍵的長度為0.110 nm.所有的密度泛函及含時(shí)密度泛函計(jì)算都是在Q-Chem商業(yè)軟件包41發(fā)展版本的框架下進(jìn)行的.

    3.1 基態(tài)電子密度

    圖1 8-ZGNR分子結(jié)構(gòu)Fig.1 Geometry of 8-ZGNR

    圖2 (a)C168H36的基態(tài)自旋密度;(b)在基態(tài)模擬中加入了一個(gè)靜電場(沿Y方向)情況下的自旋電子密度Fig.2 (a)Spin density of C168H36in ground state;(b)spin density of C168H36where a static field(along the Y direction)is added in the simulation

    首先,我們研究了有限長石墨烯納米條帶C168H36的基態(tài)性質(zhì).在UB3LYP/STO-3G的計(jì)算水平上,對C168H36分子模型進(jìn)行了自旋非限制的自洽計(jì)算,圖2為計(jì)算出的自旋密度(ρα(r)-ρβ(r)).根據(jù)以往的研究,對于zigzag型石墨烯納米條帶,自旋電子密度在其兩邊的自旋方向是相反的,因而zigzag型石墨烯納米條帶具有反鐵磁性的基態(tài)性質(zhì).但是在我們的計(jì)算中,對于C168H36分子,它的基態(tài)自旋電子密度并沒有顯示出反鐵磁的性質(zhì)(圖2(a)).盡管在相鄰的C原子之間自旋密度交替變換,但是在分子兩邊帶上的電子自旋方向并不是完全相反的,在有一些地方,電子自旋還是平行的.這與Hod等20的計(jì)算結(jié)果不同,我們認(rèn)為這首先是由于B3LYP雜化泛函所起的關(guān)鍵作用,類似于B3LYP泛函可以影響到對于石墨烯納米帶的半金屬性質(zhì)的計(jì)算結(jié)果.19其次,在我們的研究中,當(dāng)石墨烯納米條帶分子的長度較短時(shí),如C168H36,它更像是一個(gè)分子,而不能體現(xiàn)出石墨烯納米帶反鐵磁性的基態(tài)性質(zhì).為了證明長度的影響,我們又計(jì)算了C328H56的基態(tài)自旋電子密度(圖3),它的長度是C168H36的兩倍.圖3的計(jì)算結(jié)果非常清楚地顯示,α自旋和β自旋的電子密度在相鄰的C原子間交替出現(xiàn),而在石墨烯條帶分子的兩邊電子密度的自旋方向是反平行的,這說明隨著長度的增加,石墨烯納米條帶會顯示出反鐵磁的基態(tài)特性.本文中我們也考慮了加入一個(gè)垂直于石墨烯的一維方向,強(qiáng)度為0.01 a.u的靜電場的情況.結(jié)果表明對于C168H36,靜電場會大大改變它基態(tài)的性質(zhì).首先,靜電場會改變電子的自旋方向.圖2(b)是加入靜電場后,自旋電子密度的計(jì)算結(jié)果.對于C168H36,靜電場使得自旋電子密度往中間集中,并且兩邊帶上的電子自旋呈現(xiàn)出反平行的特性.其次,也對最高占據(jù)分子軌道(HOMO)與最低非占據(jù)分子軌道(LUMO)之間的帶隙進(jìn)行了分析.在不加靜電場時(shí),對于α自旋和β自旋電子帶隙是相等的,為1.17 eV.但是加入靜電場后,α自旋和β自旋電子顯示出了不同的性質(zhì),α自旋的帶隙會稍微增大,為1.22 eV,而β自旋的帶隙則會明顯減小,為0.52 eV.對于這一部分的計(jì)算,也用稍大的基組3-21G進(jìn)行了驗(yàn)證.發(fā)現(xiàn)基組的增加雖然可以大大降低基態(tài)的能量,但是對于外場的作用并沒有本質(zhì)的影響.在零外場時(shí)得到的不同自旋能量帶隙是相同的,為0.76 eV,當(dāng)電場強(qiáng)度達(dá)到0.01 a.u時(shí),α自旋的帶隙增大,為1.06 eV,而β自旋的帶隙會減小,為0.46 eV.通過比較足以表明,STO-3G基組對于加入橫向外場后不同自旋態(tài)的計(jì)算結(jié)果是合理的.Louie等16曾經(jīng)報(bào)道過,如果將石墨烯納米帶放入一個(gè)靜電場中,它的α自旋和β自旋電子的HOMO-LUMO帶隙會發(fā)生不同的變化.α自旋態(tài)的能量帶隙會增大,而β自旋態(tài)的則會減小,從而顯示出半金屬的特性.我們的計(jì)算結(jié)果也確實(shí)證明了這一變化過程,與之前的報(bào)道16也是符合的.

    圖3 C328H56的基態(tài)自旋電子密度Fig.3 Spin density of C328H56in ground sate

    在加入橫向電場后,α自旋和β自旋的不同表現(xiàn)主要是由于外電場對于不同的自旋態(tài)有著相反的影響.它會使得α自旋和β自旋電子分別朝著不同的方向重新排布,并且使得β自旋態(tài)的占據(jù)軌道和非占據(jù)軌道在能量上更加靠近,而α自旋態(tài)則分離得更遠(yuǎn).為了更好地觀察α自旋和β自旋電子對于外電場的作用所顯示出的不同表現(xiàn),在下一部分對激發(fā)態(tài)性質(zhì)的計(jì)算中,我們一直沿用了在基態(tài)計(jì)算時(shí)加入這樣一個(gè)靜電場的計(jì)算結(jié)果.

    3.2 激光場誘導(dǎo)的電子密度

    對于8-ZGNR激發(fā)態(tài)性質(zhì)的研究,考慮了在外激光場的作用下,誘導(dǎo)電子密度隨時(shí)間的演化.在對C168H36基態(tài)的計(jì)算時(shí),為了保持基態(tài)的反鐵磁性質(zhì),我們?nèi)匀槐A粢粋€(gè)強(qiáng)度為0.01 a.u.的靜電場.而在基態(tài)計(jì)算的基礎(chǔ)上,對激發(fā)態(tài)的計(jì)算模擬中,則另加入了一個(gè)隨時(shí)間變化的激光脈沖.如圖4所示,在時(shí)間域中,8-ZGNR經(jīng)歷了兩個(gè)激光脈沖,模擬的總時(shí)間為1600 a.u.,時(shí)間步長為0.4 a.u.通過對等式(3)的計(jì)算,可以在時(shí)間域中計(jì)算出誘導(dǎo)電子密度δρ(t)=ρ(t)-ρ(0).由此可以得到誘導(dǎo)電子密度,即為某一時(shí)刻t和初始時(shí)刻的電子密度之差.

    圖4 外激光場脈沖隨時(shí)間的變化Fig.4 External laser field changes with time

    在上式中,φλ代表第λ個(gè)原子軌道.圖5和圖6分別為計(jì)算得到的在一個(gè)激光脈沖中不同的四個(gè)時(shí)間點(diǎn)上α自旋和β自旋誘導(dǎo)電子密度隨時(shí)間的變化.隨著外激光場的不同,首先不論是α自旋還是β自旋誘導(dǎo)電子密度都會隨時(shí)間發(fā)生明顯的變化并且誘導(dǎo)電子會隨著外場的方向而發(fā)生移動;但是對于不同的自旋,有著不同的變化特點(diǎn).如圖5(a)所示,對于α自旋,當(dāng)t=240.8 a.u.時(shí),外激光場的強(qiáng)度接近其正向最大值(圖4中的a點(diǎn)),此時(shí)產(chǎn)生了較強(qiáng)的誘導(dǎo)電子密度,特別是在石墨烯納米帶的上下邊緣上.沿著外激光場的方向,正負(fù)誘導(dǎo)電子密度發(fā)生分離,相鄰C原子上的電子密度極性相反.當(dāng)t=384.8 a.u.時(shí),此時(shí)外激光場強(qiáng)度很小,接近于零(圖4中的b點(diǎn)).受此影響,α自旋誘導(dǎo)電子密度也變得很弱,幾乎消失不見(圖5(b)).當(dāng)t=480.8 a.u.時(shí),外激光場接近于其負(fù)最大值(圖4中的c點(diǎn)),圖5(c)顯示出誘導(dǎo)電子密度又增強(qiáng)很多,且多集中在石墨烯分子的邊緣.但是對比圖5(a)發(fā)現(xiàn),此時(shí)誘導(dǎo)電子密度的正負(fù)正好相反.這說明誘導(dǎo)電子會沿著外激光場的方向移動,當(dāng)外場的極化方向由正向轉(zhuǎn)為反向時(shí),誘導(dǎo)電子的極性也會反向.最后,當(dāng)t=768.8 a.u.時(shí),外激光場又一次接近于零(圖4中的d點(diǎn)),此時(shí)誘導(dǎo)電子密度除了在zigzag邊的中間部分,其它地方也幾乎消失不見(圖5(d)).

    圖5 在外激光場的作用下C168H36的α自旋誘導(dǎo)電子密度Fig.5 Induced α-spin charge density of C168H36interacting with an external laser field

    圖6 在外激光場的作用下C168H36的β自旋誘導(dǎo)電子密度Fig.6 Induced β-spin charge density of C168H36interacting with an external laser field

    對于β自旋,圖6顯示出誘導(dǎo)電子密度隨著外激光場變化與α自旋共同的特點(diǎn),如β自旋誘導(dǎo)電子會沿著外激光場的極化方向而發(fā)生移動,會隨著外激光場強(qiáng)度的變化而增強(qiáng)或減弱.但是對比圖5(a)和圖6(a),發(fā)現(xiàn)在同樣強(qiáng)度的外激光場的影響下,產(chǎn)生的β自旋誘導(dǎo)電子密度要比α自旋弱很多.同樣的情況在t為384.8和480.8 a.u.時(shí)刻也會出現(xiàn).但是當(dāng)t=768.8 a.u.時(shí),雖然外激光場的強(qiáng)度已接近于零,但是β自旋誘導(dǎo)電子密度并沒有出現(xiàn)明顯的減弱.從誘導(dǎo)偶極矩隨時(shí)間的變化(圖7)發(fā)現(xiàn),在t=768.8 a.u.時(shí)刻,誘導(dǎo)偶極矩出現(xiàn)了非常明顯的非絕熱效應(yīng).而圖6(d)更是證明了這一點(diǎn),此時(shí)的β自旋誘導(dǎo)電子密度已經(jīng)并不完全受外激光場的影響,非絕熱效應(yīng)使得即使外場的強(qiáng)度接近于零,誘導(dǎo)電子密度的強(qiáng)度也不會發(fā)生明顯的減弱.而同樣的情況之所以沒有出現(xiàn)在圖4(d)所示的α自旋誘導(dǎo)電子密度中,則說明了β自旋比α自旋電子態(tài)更容易出現(xiàn)非絕熱效應(yīng).

    圖7 誘導(dǎo)偶極矩隨時(shí)間的變化Fig.7 Induced dipole moment changes with time

    4 結(jié)論

    在密度泛函和含時(shí)密度泛函理論的基礎(chǔ)上,對于含有8個(gè)zigzag鏈長度較短的石墨烯納米條帶C168H36的基態(tài)和激發(fā)態(tài)的性質(zhì)進(jìn)行了研究.發(fā)現(xiàn)在基態(tài)上,B3LYP雜化泛函的計(jì)算結(jié)果顯示較短的石墨烯納米條帶并不能表現(xiàn)出反鐵磁性質(zhì),但是隨著長度的增加,納米條帶上兩邊的自旋密度會出現(xiàn)反平行,從而顯示出反鐵磁的基態(tài)性質(zhì).靜電場的加入,會降低石墨烯納米條帶C168H36的反鐵磁態(tài)的能量,成為能量最穩(wěn)定的基態(tài).同時(shí),靜電場也會使得α自旋和β自旋的HOMO-LUMO的帶隙大小出現(xiàn)不同的變化,其中α自旋帶隙的大小基本不變,而β自旋帶隙則明顯變小.這一計(jì)算結(jié)果也證明了zigzag型石墨烯納米帶的半金屬特性.

    在激發(fā)態(tài)上,加入了一個(gè)隨時(shí)間變化的激光脈沖,通過在時(shí)間域里求解含時(shí)密度泛函方程,從而觀察石墨烯納米條帶C168H36在激發(fā)態(tài)上的動力學(xué)過程.研究發(fā)現(xiàn),隨著外激光脈沖的變化,誘導(dǎo)電子密度會隨著外場強(qiáng)度的增強(qiáng)或減弱而相應(yīng)地增大或減少,并且會隨著外場的極化方向而發(fā)生明顯移動.此外,α自旋和β自旋電子在激發(fā)態(tài)上的動力學(xué)過程仍然會有不同的表現(xiàn).相比而言,在同樣強(qiáng)度的外激光場的影響下,α自旋電子更容易被激發(fā)而產(chǎn)生較明顯的誘導(dǎo)電子密度,而β自旋電子則更容易脫離外激光場的控制而產(chǎn)生非絕熱現(xiàn)象.

    (1)Novoselov,K.S.;Geim,A.K.;Morozov,S.V.;Jiang,D.;Zhang,Y.;Dubonos,S.V.;Grigorieva,I.V.;Firsov,A.A.Science 2004,306,666.doi:10.1126/science.1102896

    (2)Novoselov,K.S.;Geim,A.K.;Morozov,S.V.;Jiang,D.;Katsnelson,M.I.;Grigorieva,I.V.;Dubonos,S.V.;Firsov,A.A.Nature 2005,438,197.doi:10.1038/nature04233

    (3)Zhang,Y.;Tan,Y.W.;Stormer,H.L.;Kim,P.Nature 2005,438,201.doi:10.1038/nature04235

    (4)Novoselov,K.S.;Jiang,D.;Schedin,F.;Booth,T.J.;Khotkevich,V.V.;Morozov,S.V.;Geim,A.K.Proc.Natl.A cad.Sci.U.S.A.2005,102,10451.doi:10.1073/pnas.0502848102

    (5)Novoselov,K.S.;Jiang,Z.;Zhang,Y.;Morozov,S.V.;Stormer,H.L.;Zeitler,U.;Maan,J.C.;Boebinger,G.S.;Kim,P.;Geim,A.K.Science 2007,315,1379.doi:10.1126/science.1137201

    (6)Zhang,Y.;Small,J.P.;Pontius,W.V.;Kim,P.Appl.P hys.L ett.2005,86,073104.doi:10.1063/1.1862334

    (7)Berger,C.;Song,Z.;Li,T.;Li,X.;Ogbazghi,A.Y.;Feng,R.;Dai,Z,;Marchenkov,A.N.;Conrad,E.H.;First,P.N.;de Heer,W.A.J.Phys.Ch em.B 2004,108,19912.doi:10.1021/jp040650f

    (8)Peres,N.M.R.;Guinea,F.;Castro Neto,A.H.Phys.R ev.B 2006,73,125411.doi:10.1103/PhysRevB.73.125411

    (9)Kane,C.L.;Mele,E.J.Phys.Rev.Lett.2005,95,226801.doi:10.1103/PhysRevLett.95.226801

    (10)Pereira,V.M.;Guinea,F.;Lopes dos Santos,J.M.B.;Peres,N.M.R.;Castro Neto,A.H.Phys.Rev.Lett.2006,96,036801.doi:10.1103/PhysRevLett.96.036801

    (11)Sun,D.L.;Peng,S.L.;Ouyang,J.;Ouyang,F.P.Acta Phys.-Chim.Sin.2011,27,1103.[孫大立,彭盛霖,歐陽俊,歐陽方平.物理化學(xué)學(xué)報(bào),2011,27,1103.]doi:10.3866/PKU.WHXB20110345

    (12)Li,Y.X.;Wei,Z.D.;Zhao,Q.L.;Ding,W.;Zhang,Q.;Chen,S.G.Acta Phys.-Chim.Sin.2010,26,1570.[李云霞,魏子棟,趙巧玲,丁 煒,張 騫,陳四國.物理化學(xué)學(xué)報(bào),2010,26,1570.]doi:10.3866/PKU.WHXB20100620

    (13)Hu,Y.J.;Jin,J.;Zhang,H.;Wu,P.;Cai,C.X.A cta Phys.-Chim.Sin.2010,26,2073.[胡耀娟,金 娟,張 卉,吳 萍,蔡稱心.物理化學(xué)學(xué)報(bào),2010,26,2073.]doi:10.3866/PKU.WHXB20100812

    (14)Yang,X.;Dou,X.;Rouhanipour,A.;Zhi,L.;R?der,H.J.;Müllen,K.J.Am.Ch em.Soc.2008,130,4216.doi:10.1021/ja710234t

    (15)Son,Y.;Cohen,M.L.;Louie,S.G.P hys.R ev.L ett.2006,97,216803.doi:10.1103/PhysRevLett.97.216803

    (16)Son,Y.W.;Cohen,M.L.;Louie,S.G.Nature 2006,444,347.doi:10.1038/nature05180

    (17)Jiang,D.;Sumpter,B.G.;Dai,S.J.Chem.P hys.2007,127,124703.doi:10.1063/1.2770722

    (18)Hod,O.;Barone,V.;Peralta,J.E.;Scuseria,G.E.Nano Lett.2007,7,2295.doi:10.1021/nl0708922

    (19)Rudberg,E.;Selek,P.;Luo,Y.Nano L ett.2007,7,2211.doi:10.1021/nl070593c

    (20)Hod,O.;Barone,V.;Scuseria,G.E.Phys.Rev.B 2008,77,035411.doi:10.1103/PhysRevB.77.035411

    (21)Xu,B.;Yin,J.;Xia,Y.D.;Wan,X.G.;Jiang,K.;Liu,Z.G.Appl.P hys.Lett.2010,96,163102.doi:10.1063/1.3402762

    (23)Zhang,C.;Chen,L.;Ma,Z.Phys.Rev.B 2008,77,241402.doi:10.1103/PhysRevB.77.241402

    (24)Stauber,T.;Peres,N.M.R.;Geim,A.K.Ph ys.Rev.B 2008,78,085432.doi:10.1103/PhysRevB.78.085432

    (25)Cserti,J.;Csordas,A.;David,G.Phys.Rev.Lett.2007,99,066802.doi:10.1103/PhysRevLett.99.066802

    (26)Wright,A.R.;Liu,F.;Zhang,C.Nanotechnology 2009,20,405203.

    (27)Liu,J.;Wright,A.R.;Zhang,C.;Ma,Z.A ppl.Ph ys.Lett.2008,93,041106.doi:10.1063/1.2964093

    (28)Lu,J.Q.;Zhang,X.G.;Pantelides,S.T.Phys.Rev.B 2009,79,073408.doi:10.1103/PhysRevB.79.073408

    (29)Zhao,Z.Y.;Zhai,X.C.;Jin,G.J.Appl.Phys.L ett.2012,101,083117.doi:10.1063/1.4748110

    (30)Zhang,G.P.P hys.Rev.L ett.2005,95,047401.doi:10.1103/PhysRevLett.95.047401

    (31)Gordon,A.;Kartner,F.X;Rohringer,N.;Santra,R.Ph ys.Rev.Lett.2006,96,223902.doi:10.1103/PhysRevLett.96.223902

    (32)Onida,G.;Reining,L.;Rubio,A.Rev.Mod.P hys.2002,74,601.doi:10.1103/RevModPhys.74.601

    (33)Suzuki,M.;Mukamel,S.J.Chem.Phys.2003,119,4722.doi:10.1063/1.1594721

    (34)Press,W.H.;Teukolsky,S.A.;Vetterling,W.T.;Flannery,B.P.Numerical Recipes in Fortran 90;Cambridge University Press:Cambridge,1996.

    (35)Li,X.;Smith,S.M.;Markevitch,A.N.;Romanov,D.A.;Levis,R.J.;Schlegel,H.B.Phys.Chem.Ch em.Phys.2005,7,233.doi:10.1039/b415849k

    (36)Sun,J.;Song,J.;Zhao,Y.;Liang,W.Z.J.Ch em.Phys.2007,127,234107.doi:10.1063/1.2805396

    (37)Sun,J.;Guo,Z.Y.;Liang,W.Z.Ph ys.Rev.B 2007,75,195438.doi:10.1103/PhysRevB.75.195438

    (38)Sun,J.;Liang,W.Z.;Liu,J.J.T h eo.Comp.Chem.2008,7,579.doi:10.1142/S0219633608003976

    (39)Sun,J.;Liu,J.;Liang,W.Z.;Zhao,Y.J.Ph ys.Chem.A 2008,112,10442.doi:10.1021/jp804408a

    (40)Sun,J.;Liang,W.Z.;Yang,J.L.;Gao,J.S.J.Mol.S truct.-Theochem 2005,755,105.doi:10.1016/j.theochem.2005.07.030

    (41)Shao,Y.;Molnar,L.F.;Jung,Y.;Kussmann,J.;Ochsenfeld,C.;Brown,S.T.;Gilbert,A.T.B.;Slipchenko,L.V.;Levchenko,S.V.;O′Neill,D.P.;DiStasio,R.A.;Lochan,R.C.;Wang,T.;Beran,G.J.O.;Besley,N.A.;Herbert,J.M.;Lin,C.Y.;Van Voorhis,T.;Chien,S.H.;Sodt,A.;Steele,R.P.;Rassolov,V.A.;Maslen,P.E.;Korambath,P.P.;Adamson,R.D.;Austin,B.;Baker,J.;Byrd,E.F.C.;Dachsel,H.;Doerksen,R.J.;Dreuw,A.;Dunietz,B.D.;Dutoi,A.D.;Furlani,T.R.;Gwaltney,S.R.;Heyden,A.;Hirata,S.;Hsu,C.P.;Kedziora,G.;Khalliulin,R.Z.;Klunzinger,P.;Lee,A.M.;Lee,M.S.;Liang,W.;Lotan,I.;Nair,N.;Peters,B.;Proynov,E.I.;Pieniazek,P.A.;Rhee,Y.M.;Ritchie,J.;Rosta,E.;Sherrill,C.D.;Simmonett,A.C.;Subotnik,J.E.;Woodcock,H.L.;Zhang,W.;Bell,A.T.;Chakraborty,A.K.,Chipman,D.M.;Keil,F.J.;Warshel,A.;Hehre,W.J.;Schaefer,H.F.;Kong,J.;Krylov,A.I.;Gill,P.M.W.;Head-Gordon,M.Phys.Ch em.Chem.Phys.2006,8,3172.doi:10.1039/b517914a

    猜你喜歡
    電子密度基態(tài)條帶
    一類非線性Choquard方程基態(tài)解的存在性
    擬相對論薛定諤方程基態(tài)解的存在性與爆破行為
    一類反應(yīng)擴(kuò)散方程的Nehari-Pankov型基態(tài)解
    非線性臨界Kirchhoff型問題的正基態(tài)解
    顧及地磁影響的GNSS電離層層析不等像素間距算法*
    不同GPS掩星電離層剖面產(chǎn)品相關(guān)性分析
    等離子體電子密度分布信息提取方法研究
    一種適用于電離層電子密度重構(gòu)的AMART算法
    基于條帶模式GEOSAR-TOPS模式UAVSAR的雙基成像算法
    基于 Savitzky-Golay 加權(quán)擬合的紅外圖像非均勻性條帶校正方法
    十八禁人妻一区二区| videosex国产| 国产成人一区二区在线| 国产精品一区二区在线观看99| 亚洲中文av在线| 国产男女超爽视频在线观看| 国产精品av久久久久免费| 午夜福利免费观看在线| 日韩人妻精品一区2区三区| 男人操女人黄网站| 亚洲免费av在线视频| 欧美国产精品一级二级三级| 中国国产av一级| 欧美国产精品一级二级三级| av国产精品久久久久影院| 国语对白做爰xxxⅹ性视频网站| 少妇 在线观看| 国产高清视频在线播放一区 | 久久久久久久国产电影| 国产精品麻豆人妻色哟哟久久| 日本五十路高清| 国产日韩欧美在线精品| 尾随美女入室| 一级片免费观看大全| 18禁黄网站禁片午夜丰满| xxx大片免费视频| 人人澡人人妻人| 国语对白做爰xxxⅹ性视频网站| 久久精品国产综合久久久| 国产精品99久久99久久久不卡| 好男人电影高清在线观看| 美女午夜性视频免费| 国产在线视频一区二区| 亚洲av美国av| 交换朋友夫妻互换小说| 国产一区二区在线观看av| 国产精品亚洲av一区麻豆| 久久影院123| 一级毛片我不卡| 一级黄色大片毛片| 极品少妇高潮喷水抽搐| 久久av网站| 亚洲九九香蕉| 如日韩欧美国产精品一区二区三区| 国产成人一区二区三区免费视频网站 | 99re6热这里在线精品视频| 久9热在线精品视频| 国产免费又黄又爽又色| 亚洲男人天堂网一区| 精品人妻熟女毛片av久久网站| 女人久久www免费人成看片| 亚洲精品自拍成人| 91麻豆av在线| 王馨瑶露胸无遮挡在线观看| 黄色视频在线播放观看不卡| 亚洲午夜精品一区,二区,三区| 亚洲精品在线美女| 在线观看免费视频网站a站| 久久天堂一区二区三区四区| videosex国产| 丝袜在线中文字幕| 精品亚洲成国产av| 午夜福利免费观看在线| 国产精品成人在线| 国产亚洲午夜精品一区二区久久| 黄色毛片三级朝国网站| 老司机在亚洲福利影院| 伊人久久大香线蕉亚洲五| 80岁老熟妇乱子伦牲交| 大话2 男鬼变身卡| 日本一区二区免费在线视频| 各种免费的搞黄视频| 欧美日韩成人在线一区二区| 国产精品成人在线| 国产精品欧美亚洲77777| 丰满少妇做爰视频| 18禁观看日本| 宅男免费午夜| 男女边摸边吃奶| 婷婷色综合大香蕉| 美女扒开内裤让男人捅视频| 日日爽夜夜爽网站| 婷婷色av中文字幕| 久久久欧美国产精品| 最近中文字幕2019免费版| 丝袜人妻中文字幕| 午夜久久久在线观看| 男女国产视频网站| 亚洲精品日本国产第一区| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区| 成年人免费黄色播放视频| 满18在线观看网站| 日本一区二区免费在线视频| 久久午夜综合久久蜜桃| 美国免费a级毛片| 国产精品一二三区在线看| 悠悠久久av| 国产99久久九九免费精品| 欧美日韩黄片免| 两个人免费观看高清视频| 精品国产一区二区三区久久久樱花| 男女之事视频高清在线观看 | 人人妻人人澡人人看| 日韩一卡2卡3卡4卡2021年| 中文精品一卡2卡3卡4更新| 欧美日韩视频高清一区二区三区二| 女性被躁到高潮视频| 成年av动漫网址| 王馨瑶露胸无遮挡在线观看| 中文乱码字字幕精品一区二区三区| 丝袜美足系列| 一级毛片电影观看| 成人影院久久| 国产精品成人在线| 亚洲精品美女久久久久99蜜臀 | 免费在线观看视频国产中文字幕亚洲 | 婷婷成人精品国产| 999久久久国产精品视频| 又黄又粗又硬又大视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美色中文字幕在线| 看免费av毛片| 欧美国产精品一级二级三级| 国产极品粉嫩免费观看在线| videosex国产| 人人妻人人澡人人爽人人夜夜| 美女高潮到喷水免费观看| 美女脱内裤让男人舔精品视频| 欧美在线一区亚洲| 亚洲欧美成人综合另类久久久| 久久久国产一区二区| 男女高潮啪啪啪动态图| 在线观看免费高清a一片| 一级黄片播放器| 中文字幕亚洲精品专区| 免费黄频网站在线观看国产| 久久精品aⅴ一区二区三区四区| 国产一区有黄有色的免费视频| 巨乳人妻的诱惑在线观看| 欧美黑人精品巨大| 亚洲久久久国产精品| 视频区欧美日本亚洲| av视频免费观看在线观看| 老熟女久久久| 欧美成人午夜精品| 欧美日韩av久久| 热re99久久国产66热| 好男人视频免费观看在线| 亚洲欧美一区二区三区久久| www日本在线高清视频| 精品福利观看| 每晚都被弄得嗷嗷叫到高潮| 午夜福利视频在线观看免费| 国产午夜精品一二区理论片| 午夜福利,免费看| 久久女婷五月综合色啪小说| 一级毛片黄色毛片免费观看视频| 免费观看人在逋| 精品国产一区二区三区久久久樱花| 国产成人精品久久二区二区91| 欧美亚洲 丝袜 人妻 在线| 欧美精品一区二区免费开放| 少妇被粗大的猛进出69影院| 亚洲,欧美,日韩| 一本一本久久a久久精品综合妖精| 免费在线观看影片大全网站 | 欧美久久黑人一区二区| 青春草亚洲视频在线观看| 国产av国产精品国产| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品一区三区| 亚洲国产最新在线播放| 超碰成人久久| 国产免费一区二区三区四区乱码| 精品少妇一区二区三区视频日本电影| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品一区蜜桃| 成年av动漫网址| 女人高潮潮喷娇喘18禁视频| 国产成人系列免费观看| 国产91精品成人一区二区三区 | 在线看a的网站| 69精品国产乱码久久久| www.精华液| 亚洲欧美一区二区三区久久| 飞空精品影院首页| 99久久人妻综合| 国产亚洲精品第一综合不卡| 亚洲av片天天在线观看| 我的亚洲天堂| 亚洲 欧美一区二区三区| 男女高潮啪啪啪动态图| 人妻人人澡人人爽人人| 9191精品国产免费久久| 好男人视频免费观看在线| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧洲国产日韩| 精品人妻熟女毛片av久久网站| 国产成人精品在线电影| 国产精品av久久久久免费| 久久女婷五月综合色啪小说| 欧美激情高清一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产男女超爽视频在线观看| 十八禁人妻一区二区| 久久热在线av| 国产精品一区二区精品视频观看| 欧美老熟妇乱子伦牲交| 日本午夜av视频| 亚洲中文日韩欧美视频| 免费黄频网站在线观看国产| 亚洲欧美一区二区三区黑人| 婷婷色综合大香蕉| 一本大道久久a久久精品| 新久久久久国产一级毛片| 99精国产麻豆久久婷婷| 一本—道久久a久久精品蜜桃钙片| av不卡在线播放| 日韩 欧美 亚洲 中文字幕| 亚洲国产欧美一区二区综合| av一本久久久久| 亚洲国产最新在线播放| 免费黄频网站在线观看国产| 高清黄色对白视频在线免费看| 国产精品亚洲av一区麻豆| 咕卡用的链子| avwww免费| 久久天躁狠狠躁夜夜2o2o | 一个人免费看片子| 精品免费久久久久久久清纯 | 在现免费观看毛片| 亚洲精品久久成人aⅴ小说| 亚洲精品自拍成人| 精品第一国产精品| 亚洲专区中文字幕在线| 高清视频免费观看一区二区| 69精品国产乱码久久久| 天天躁日日躁夜夜躁夜夜| 免费少妇av软件| 亚洲国产欧美一区二区综合| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 国产免费又黄又爽又色| 51午夜福利影视在线观看| 日本wwww免费看| 又黄又粗又硬又大视频| 赤兔流量卡办理| 久久天躁狠狠躁夜夜2o2o | 亚洲熟女毛片儿| 少妇精品久久久久久久| 色精品久久人妻99蜜桃| av国产精品久久久久影院| 亚洲精品第二区| 日本av手机在线免费观看| 国产成人精品久久二区二区免费| 精品国产一区二区三区四区第35| 观看av在线不卡| 日韩免费高清中文字幕av| 国产又色又爽无遮挡免| 大码成人一级视频| 亚洲国产成人一精品久久久| 欧美在线一区亚洲| 好男人视频免费观看在线| 国产一卡二卡三卡精品| 啦啦啦在线观看免费高清www| 亚洲成人免费电影在线观看 | 亚洲欧美清纯卡通| 免费在线观看完整版高清| 操美女的视频在线观看| 亚洲国产成人一精品久久久| 母亲3免费完整高清在线观看| 脱女人内裤的视频| 久久亚洲国产成人精品v| 亚洲精品乱久久久久久| 十八禁高潮呻吟视频| 免费高清在线观看视频在线观看| 欧美性长视频在线观看| 老司机影院毛片| 久久久久久久久免费视频了| 日韩视频在线欧美| 亚洲av成人精品一二三区| 亚洲成人免费电影在线观看 | 国产亚洲一区二区精品| 人成视频在线观看免费观看| 国产免费视频播放在线视频| 国产伦理片在线播放av一区| 国产男人的电影天堂91| 精品人妻1区二区| 如日韩欧美国产精品一区二区三区| 亚洲午夜精品一区,二区,三区| 国产男女超爽视频在线观看| 免费在线观看日本一区| 亚洲av美国av| 中国国产av一级| 日韩 欧美 亚洲 中文字幕| 日韩制服骚丝袜av| 午夜久久久在线观看| 亚洲国产精品国产精品| 中文字幕人妻熟女乱码| 啦啦啦啦在线视频资源| 欧美日韩成人在线一区二区| 老司机在亚洲福利影院| 国产成人精品无人区| 丰满少妇做爰视频| 大陆偷拍与自拍| 精品少妇久久久久久888优播| 天天影视国产精品| 操美女的视频在线观看| 女人被躁到高潮嗷嗷叫费观| 一区二区三区精品91| 美女国产高潮福利片在线看| 狠狠婷婷综合久久久久久88av| 老司机靠b影院| 日韩一区二区三区影片| 人妻人人澡人人爽人人| 亚洲精品一卡2卡三卡4卡5卡 | 国产欧美日韩一区二区三 | 一本—道久久a久久精品蜜桃钙片| 亚洲成人免费电影在线观看 | 欧美精品亚洲一区二区| 国产精品一区二区在线观看99| 亚洲欧洲日产国产| 波多野结衣av一区二区av| 国产熟女欧美一区二区| 午夜福利视频精品| 国产在线视频一区二区| 亚洲av电影在线观看一区二区三区| av在线app专区| 亚洲熟女精品中文字幕| 一区二区三区激情视频| 最近最新中文字幕大全免费视频 | 看免费av毛片| 日韩视频在线欧美| 亚洲精品在线美女| 人人妻,人人澡人人爽秒播 | 无限看片的www在线观看| 亚洲,欧美,日韩| 高清av免费在线| 国产成人欧美| 国产精品成人在线| 欧美精品亚洲一区二区| 欧美精品人与动牲交sv欧美| 婷婷色综合www| 美女高潮到喷水免费观看| 99热网站在线观看| 新久久久久国产一级毛片| 波多野结衣一区麻豆| 亚洲欧美日韩另类电影网站| 高清黄色对白视频在线免费看| 91精品伊人久久大香线蕉| 国产一卡二卡三卡精品| www.精华液| 天天躁狠狠躁夜夜躁狠狠躁| 久久这里只有精品19| 999精品在线视频| 婷婷色麻豆天堂久久| 国产老妇伦熟女老妇高清| 叶爱在线成人免费视频播放| 啦啦啦在线免费观看视频4| 欧美黑人精品巨大| 午夜老司机福利片| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| 水蜜桃什么品种好| 男的添女的下面高潮视频| 一本大道久久a久久精品| videosex国产| 日日夜夜操网爽| 欧美人与善性xxx| 欧美精品av麻豆av| 国产精品国产av在线观看| 国产精品人妻久久久影院| 十八禁高潮呻吟视频| 99热国产这里只有精品6| 国产高清视频在线播放一区 | 精品第一国产精品| 宅男免费午夜| 日韩伦理黄色片| 黄网站色视频无遮挡免费观看| 一级,二级,三级黄色视频| 欧美性长视频在线观看| 天堂中文最新版在线下载| av电影中文网址| 国产高清不卡午夜福利| 丝袜美足系列| 欧美日韩综合久久久久久| 啦啦啦在线免费观看视频4| www.999成人在线观看| 制服人妻中文乱码| 欧美老熟妇乱子伦牲交| 久久国产精品人妻蜜桃| 久久av网站| 在线观看国产h片| 成人午夜精彩视频在线观看| 熟女少妇亚洲综合色aaa.| 国产又爽黄色视频| 啦啦啦在线观看免费高清www| 国产在线观看jvid| 中文精品一卡2卡3卡4更新| tube8黄色片| 亚洲免费av在线视频| 99久久人妻综合| 午夜老司机福利片| 男女午夜视频在线观看| 国产精品亚洲av一区麻豆| 精品亚洲成a人片在线观看| 国产成人一区二区三区免费视频网站 | 午夜福利影视在线免费观看| 久久人妻福利社区极品人妻图片 | 亚洲一区二区三区欧美精品| 一二三四在线观看免费中文在| 永久免费av网站大全| 色婷婷久久久亚洲欧美| 满18在线观看网站| 狂野欧美激情性bbbbbb| 久久精品成人免费网站| 丰满饥渴人妻一区二区三| 丝袜人妻中文字幕| 一边摸一边做爽爽视频免费| 欧美日韩成人在线一区二区| 亚洲精品国产区一区二| 亚洲人成电影免费在线| 老熟女久久久| 首页视频小说图片口味搜索 | 久久ye,这里只有精品| av视频免费观看在线观看| 国产精品 欧美亚洲| 欧美黑人精品巨大| 丰满迷人的少妇在线观看| 亚洲第一青青草原| 午夜免费鲁丝| www.999成人在线观看| 久久精品久久精品一区二区三区| 国产在线视频一区二区| 国产成人精品久久二区二区免费| 午夜福利在线免费观看网站| 一本大道久久a久久精品| 亚洲欧美一区二区三区黑人| 一边摸一边抽搐一进一出视频| 久久99一区二区三区| 国产片特级美女逼逼视频| 欧美在线一区亚洲| 亚洲中文字幕日韩| 亚洲天堂av无毛| 777久久人妻少妇嫩草av网站| 国产成人精品久久二区二区91| 一边摸一边抽搐一进一出视频| 久久久久久久久久久久大奶| 少妇 在线观看| 免费看十八禁软件| 又紧又爽又黄一区二区| 久久九九热精品免费| 女人被躁到高潮嗷嗷叫费观| 国产国语露脸激情在线看| 成人国产一区最新在线观看 | 久久精品熟女亚洲av麻豆精品| 美女扒开内裤让男人捅视频| 蜜桃在线观看..| 国产一区二区 视频在线| 亚洲,欧美,日韩| 国产成人免费无遮挡视频| 啦啦啦 在线观看视频| 日日摸夜夜添夜夜爱| 成人亚洲欧美一区二区av| 亚洲精品国产av成人精品| 亚洲欧美一区二区三区黑人| 欧美日韩亚洲高清精品| 亚洲av美国av| 高潮久久久久久久久久久不卡| 欧美日本中文国产一区发布| 99久久综合免费| 美女脱内裤让男人舔精品视频| 五月天丁香电影| 日韩av在线免费看完整版不卡| 精品福利观看| 又大又爽又粗| 久久久久久久大尺度免费视频| 一区二区av电影网| 精品一区二区三卡| 丰满人妻熟妇乱又伦精品不卡| 国产伦理片在线播放av一区| 啦啦啦视频在线资源免费观看| 岛国毛片在线播放| 我要看黄色一级片免费的| 精品免费久久久久久久清纯 | 久久久亚洲精品成人影院| 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| av线在线观看网站| 国产成人欧美在线观看 | 青春草亚洲视频在线观看| 日韩,欧美,国产一区二区三区| 亚洲av片天天在线观看| cao死你这个sao货| 亚洲精品在线美女| 七月丁香在线播放| 日韩欧美一区视频在线观看| 久久精品久久精品一区二区三区| 狠狠婷婷综合久久久久久88av| 久久久久国产精品人妻一区二区| 亚洲天堂av无毛| 国产野战对白在线观看| 夫妻午夜视频| 国产一区二区三区综合在线观看| 美国免费a级毛片| 日韩av不卡免费在线播放| 美女视频免费永久观看网站| 国产精品亚洲av一区麻豆| 国产又色又爽无遮挡免| 国产免费视频播放在线视频| 亚洲欧美一区二区三区黑人| 啦啦啦啦在线视频资源| 最新在线观看一区二区三区 | 久久久久久亚洲精品国产蜜桃av| 人妻 亚洲 视频| 色视频在线一区二区三区| 美女脱内裤让男人舔精品视频| 一二三四在线观看免费中文在| 别揉我奶头~嗯~啊~动态视频 | 欧美激情高清一区二区三区| 一级黄片播放器| 日韩伦理黄色片| 男人添女人高潮全过程视频| 色视频在线一区二区三区| 欧美激情高清一区二区三区| 一本综合久久免费| 国产在线免费精品| av一本久久久久| 久久久国产欧美日韩av| 制服诱惑二区| 欧美日韩av久久| 久久精品人人爽人人爽视色| 久久午夜综合久久蜜桃| 国产日韩欧美在线精品| 亚洲一卡2卡3卡4卡5卡精品中文| 最近手机中文字幕大全| videosex国产| 啦啦啦中文免费视频观看日本| 午夜福利一区二区在线看| 国产麻豆69| 三上悠亚av全集在线观看| 国产色视频综合| 亚洲国产精品成人久久小说| 色94色欧美一区二区| 青青草视频在线视频观看| 搡老岳熟女国产| 1024视频免费在线观看| 国产精品 欧美亚洲| 亚洲av综合色区一区| 欧美+亚洲+日韩+国产| 中国美女看黄片| 丰满少妇做爰视频| 欧美黑人欧美精品刺激| 欧美另类一区| 在线亚洲精品国产二区图片欧美| 日韩熟女老妇一区二区性免费视频| 国产男女内射视频| 亚洲国产成人一精品久久久| 亚洲欧美一区二区三区黑人| 日本av手机在线免费观看| 18禁国产床啪视频网站| 久久久久视频综合| 久久久久久久国产电影| 中国美女看黄片| 国语对白做爰xxxⅹ性视频网站| 母亲3免费完整高清在线观看| 天天躁夜夜躁狠狠躁躁| 午夜激情久久久久久久| 久久天堂一区二区三区四区| 亚洲国产最新在线播放| 在线观看免费日韩欧美大片| 久热爱精品视频在线9| 黄色一级大片看看| 成年美女黄网站色视频大全免费| 亚洲 欧美一区二区三区| 免费看不卡的av| 午夜激情久久久久久久| 黄片小视频在线播放| 人人妻人人添人人爽欧美一区卜| 久久亚洲精品不卡| av又黄又爽大尺度在线免费看| 97在线人人人人妻| 一本色道久久久久久精品综合| 日本五十路高清| 国产精品亚洲av一区麻豆| 久久人人97超碰香蕉20202| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 日韩伦理黄色片| 精品亚洲成国产av| 人人妻人人澡人人看| 丝袜脚勾引网站| 嫁个100分男人电影在线观看 | 一级黄片播放器| 婷婷色综合大香蕉| 飞空精品影院首页| 色婷婷av一区二区三区视频| 精品亚洲成国产av| 一区二区三区四区激情视频| 精品国产乱码久久久久久男人| 婷婷成人精品国产| 亚洲av成人精品一二三区| 亚洲视频免费观看视频| 亚洲精品国产av成人精品| 婷婷色麻豆天堂久久| 天天影视国产精品| 国产高清videossex| 国产伦人伦偷精品视频| 97人妻天天添夜夜摸| 国产成人免费观看mmmm| 美女视频免费永久观看网站| 一区二区三区精品91| 国产成人精品久久二区二区免费| av网站免费在线观看视频|