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

    類富勒烯團(tuán)簇發(fā)光性能的理論研究*

    2022-07-19 07:45:34楊小偉佘潔周思趙紀(jì)軍
    物理學(xué)報(bào) 2022年12期
    關(guān)鍵詞:能隙富勒烯氮化硼

    楊小偉 佘潔 周思 趙紀(jì)軍

    (大連理工大學(xué)三束材料改性教育部重點(diǎn)實(shí)驗(yàn)室,大連 116024)

    熱激活延遲熒光(TADF)作為一種特殊的分子熒光機(jī)制,對(duì)于提高發(fā)光效率有著重要意義.以C60 和C70 為代表的碳富勒烯具有高對(duì)稱結(jié)構(gòu)和離域π 電子,被廣泛證明具有顯著的TADF 效應(yīng);相比之下,其他類富勒烯團(tuán)簇的光物理性質(zhì)尚不清楚.本文利用含時(shí)密度泛函理論探索了一系列類富勒烯團(tuán)簇的激發(fā)態(tài)性質(zhì),包括實(shí)驗(yàn)合成的具有不同尺寸的氮化硼籠型團(tuán)簇B12N12,B24N24 和B36N36,以及與B12N12 結(jié)構(gòu)相同、元素組成不同的B12P12,Al12N12 和Ga12N12.計(jì)算結(jié)果表明,這些類富勒烯化合物團(tuán)簇具有2.83—6.54 eV 的能隙,主要吸收紫外光,熒光發(fā)射波長在可見光區(qū)間,包括紅光、橙光、藍(lán)光和紫光.它們的第一激發(fā)單重態(tài)和三重態(tài)的能量差較小(ΔEST <0.29 eV),因此可能通過系間竄越和反系間竄越發(fā)生TADF 過程.導(dǎo)致ΔEST 較低的原因是這些化合物團(tuán)簇的最高占據(jù)分子軌道和最低未占據(jù)分子軌道分布在不同的元素上,使得二者重疊度較低.這些理論結(jié)果揭示了類富勒烯團(tuán)簇的發(fā)光性質(zhì)和可能的熒光機(jī)理,為尋找和設(shè)計(jì)穩(wěn)定高效的發(fā)光材料提供了重要指導(dǎo).

    1 引言

    熱激活延遲熒光(thermally activated delayed fluorescence,TADF)是一種特殊的分子發(fā)光機(jī)制,對(duì)于提高發(fā)光效率有著重要意義.常規(guī)的熒光機(jī)制,即瞬時(shí)熒光(prompt fluorescence),是從第一激發(fā)單重態(tài)S1回到基態(tài)S0而發(fā)生的輻射躍遷,由于第一激發(fā)三重態(tài)T1的能量比S1態(tài)更低,理論上只能形成1/4 的單重態(tài)激子,導(dǎo)致熒光內(nèi)量子效率最高只能達(dá)到25%.而當(dāng)S1態(tài)與T1態(tài)能量相近時(shí),在熱激發(fā)下通過系間竄越和反系間竄越可以實(shí)現(xiàn)S1→T1→S1之間的躍遷,最終由S1態(tài)退激發(fā)至基態(tài)而發(fā)生延遲熒光,三重態(tài)激子得以被利用,使得熒光內(nèi)量子效率可達(dá)100%[1,2].通常,TADF相比瞬時(shí)熒光的發(fā)光強(qiáng)度要弱,而壽命顯著延長.1961 年,Parker 和Hatchard[3]發(fā)現(xiàn)伊紅染料(四溴熒光素)能夠發(fā)出微弱的長壽命熒光,首次觀察到了TADF 現(xiàn)象.隨后,呫噸染料[4]、芳香醛和芳香酮類化合物[5,6]、金屬卟啉化合物[7]及芳烴化合物[8,9]等分子體系中的TADF 過程均被實(shí)驗(yàn)證實(shí).2012 年,Adachi 等[10]開發(fā)了由供體和受體結(jié)構(gòu)組成的有機(jī)分子TADF 材料,由于最高占據(jù)軌道(highest occupied molecular orbitals,HOMO)和最低未占軌道(lowest unoccupied molecular orbitals,LUMO)分別局域在供體和受體單元上,軌道重疊較低導(dǎo)致較小的S1態(tài)與T1態(tài)的能量差,這種利用HOMO-LUMO 空間分離的策略已被廣泛證實(shí)可以有效提高TADF 效率[11?14].

    以C60和C70為代表的碳富勒烯因具有獨(dú)特的籠狀結(jié)構(gòu)和離域的π 電子而展現(xiàn)出獨(dú)特的光物理性質(zhì),特別是顯著的TADF 效應(yīng).它們的S1態(tài)與T1態(tài)的能量差較小,S1→T1系間竄越非常迅速,T1態(tài)的量子產(chǎn)率可接近于1 且壽命較長(約24 ms),是發(fā)生TADF 的必要條件[15].C70自身的熒光峰在650—725 nm 范圍[16],并且具有磷光特性,波長在810 nm,磷光壽命約53 ms[17,18].C70的TADF 過程首先被Berberan-Santos 和Garcia[19]在實(shí)驗(yàn)上證實(shí),其TADF 發(fā)光強(qiáng)度遠(yuǎn)高于瞬時(shí)熒光強(qiáng)度,不同溫度下可達(dá)50—167 倍.C60自身的熒光峰在650—760 nm 范圍[20],其延遲熒光強(qiáng)度比C70低,與瞬時(shí)熒光強(qiáng)度相當(dāng)[15,21?23].它們的TADF 特性可應(yīng)用于溫度傳感器和氧探測(cè)器,在較大的溫度區(qū)間都有很高的靈敏度[24].

    繼碳富勒烯之后,許多類富勒烯團(tuán)簇也被預(yù)言和發(fā)現(xiàn),例如過渡金屬內(nèi)嵌摻雜的籠型團(tuán)簇,M@Si16和M@Ge16(M=Ti,Zn,Hf),M@Au12(M=Mo,W),M@Al12(M=Sc-Zn),M@Ag12(M=Sc-Cd)和M@Cu12(M=Sc-Ag)等,與碳富勒烯相似,它們具有高對(duì)稱性、良好的穩(wěn)定性及較大的能隙[25,26].特別地,氮化硼可以像碳一樣形成穩(wěn)定的中空籠型結(jié)構(gòu),Oku 等[27?31]利用電弧熔煉法成功獲得了一系列具有不同尺寸的氮化硼類富勒烯團(tuán)簇,包括B12N12[27],B24N24[28],B28N28[29]和B36N36[30,31],并利用高分辨電子顯微鏡和飛行時(shí)間質(zhì)譜等手段確定了其籠狀結(jié)構(gòu),它們由四元環(huán)、六元環(huán)或八元環(huán)構(gòu)成,滿足孤立四方穩(wěn)定性規(guī)則(isolated tetrahedral rule)[32].這些氮化硼團(tuán)簇具有高穩(wěn)定性,HOMO-LUMO 能隙接近6 eV.此外,實(shí)驗(yàn)還發(fā)現(xiàn)硼氮納米膠囊包裹Ag 納米顆粒具有光致發(fā)光特性,發(fā)光波長約336 nm[33].然而,氮化硼類富勒烯團(tuán)簇自身的光物理性質(zhì)如何、是否與碳富勒烯一樣具有TADF 特性,尚未有人研究.

    本文通過第一性原理計(jì)算探索了一系列類富勒烯團(tuán)簇的發(fā)光性質(zhì),包括B12N12,B24N24,B36N36,B12P12,Al12N12和Ga12N12,計(jì)算了它 們?cè)诨鶓B(tài)、第一激發(fā)單重態(tài)和三重態(tài)下的幾何和電子結(jié)構(gòu),模擬了吸收光譜和熒光發(fā)射光譜,并考察了S1態(tài)與T1態(tài)的能量差、自旋-軌道耦合常數(shù)、重組能等關(guān)鍵物理量,以評(píng)估其TADF 性能,為理解類富勒烯團(tuán)簇的光物理性質(zhì)提供了重要理論參考.

    2 計(jì)算方法

    本文考慮了實(shí)驗(yàn)已報(bào)道的B12N12,B24N24和B36N36籠型團(tuán)簇,以及與B12N12具有相同結(jié)構(gòu)的B12P12,Al12N12和Ga12N12團(tuán)簇.采用密度泛函理論(density functional theory,DFT)[34]方法計(jì)算這6 種團(tuán)簇的基態(tài)結(jié)構(gòu)和電子結(jié)構(gòu),采用含時(shí)密度泛函理論(time-dependent density functional theory,TD-DFT)[35],優(yōu)化得到團(tuán)簇單重激發(fā)態(tài)和三重激發(fā)態(tài)的結(jié)構(gòu).所有結(jié)構(gòu)優(yōu)化、分子軌道、吸收光譜和發(fā)射光譜的計(jì)算均采用Gaussian16 軟件包[36],HSE06 雜化泛函和6-311G*基組,結(jié)構(gòu)建模采用Device Studio[37],發(fā)光譜分析采用MOMAP[38]軟件.我們過去的系列性工作證明HSE06 雜化泛函可以較準(zhǔn)確地模擬多種籠型團(tuán)簇的基態(tài)結(jié)構(gòu)和光電子譜,包括過渡金屬內(nèi)嵌摻雜的硅籠團(tuán)簇[39,40]及硼團(tuán)簇[41,42].這里我們用六方氮化硼和閃鋅礦氮化硼晶體作為參照,測(cè)試了不同泛函(PBE,PBE0,PW91,B3lYP,HSE06)對(duì)能帶結(jié)構(gòu)的影響,其中HSE06 雜化泛函計(jì)算預(yù)言兩種晶體的能隙分別是5.30 和5.70 eV,與實(shí)驗(yàn)值(5.97 和6.10 eV)[43,44]最為接近,詳細(xì)結(jié)果見補(bǔ)充材料表 S1 (online).

    我們還測(cè)試了泛函對(duì)TD-DFT 計(jì)算團(tuán)簇激發(fā)態(tài)性質(zhì)的影響,具體考察了單重激發(fā)態(tài)S1和三重激發(fā)態(tài)T1的結(jié)構(gòu)、能量和對(duì)應(yīng)的發(fā)射光波長.定義S1態(tài)和T1態(tài)的絕熱激發(fā)能分別為2 個(gè)激發(fā)態(tài)與基態(tài)S0的能量差[45]:

    式中,ES1/T1分別是TD-DFT 優(yōu)化得到的2 種激發(fā)態(tài)結(jié)構(gòu)的能量,ES0是基態(tài)S0的能量.?ES1即對(duì)應(yīng)發(fā)射光的能量(以及波長λ).通常,S1態(tài)和T1態(tài)之間的絕熱能量差(?EST)用于評(píng)估熱激活延遲熒光是否容易發(fā)生,定義為

    ES1和ET1能量越接近,即?EST的數(shù)值越小,激發(fā)電子通過反系間竄越從T1態(tài)至S1態(tài)的概率越大,TADF 越容易發(fā)生.

    眾所周知,采用TD-DFT 描述體系激發(fā)態(tài)時(shí),激發(fā)態(tài)的性質(zhì)對(duì)計(jì)算選擇的交換關(guān)聯(lián)泛函非常敏感,會(huì)直接影響發(fā)射光能量、S1態(tài)和T1態(tài)之間的能量差等[46].于是我們考慮了5 種泛函,含有不同Hartree-Fock (HF)交換的百分比,包括B3LYP(20%HF),PBE0 (25%HF),BMK (42%HF),M06-2X (54%HF),M06-HF (100%HF),以及短程HF交換的雜化泛函HSE06 和長程校正泛函CAMB3LYP.以B12N12為例,補(bǔ)充材料表S2 (online)列出了不同泛函計(jì)算的?ES1,λ和?EST.可以看出:?EST的范圍在0.20—0.25 eV,受泛函影響不大;?ES1隨著HF 百分比的增加而增大,范圍在2.40—2.99 eV,相應(yīng)發(fā)射光波長變短,分布在413.73—516.47 nm;總的來說,PBE0 和B3LYP (HF% <25%)以及HSE06 雜化泛函預(yù)言的?ES1和λ值較為接近.過去的研究表明,B3LYP 雜化泛函計(jì)算碳富勒烯的光譜性質(zhì)與實(shí)驗(yàn)測(cè)量結(jié)果[42,47,48]符合得較好.因此,本文采用HSE06 雜化泛函研究氮化硼籠型團(tuán)簇的基態(tài)和激發(fā)態(tài)性質(zhì).

    3 結(jié)果與討論

    3.1 類富勒烯團(tuán)簇的基態(tài)和激發(fā)態(tài)的幾何結(jié)構(gòu)

    首先考察了實(shí)驗(yàn)報(bào)道的不同尺寸的氮化硼籠型團(tuán)簇B12N12,B24N24和B36N36的結(jié)構(gòu),如圖1 所示,其關(guān)鍵結(jié)構(gòu)參數(shù)總結(jié)在表1 中.其中,B12N12團(tuán)簇由6 個(gè)四元環(huán)和8 個(gè)六元環(huán)組成,具有Th對(duì)稱性,直徑約4.80 ?;B24N24團(tuán)簇由12 個(gè)四元環(huán)、8 個(gè)六元環(huán)和6 個(gè)八元環(huán)組成,具有S8對(duì)稱性,直徑約6.68 ?;而B36N36團(tuán)簇是由6 個(gè)四元環(huán)和32 個(gè)六元環(huán)組成的具有Oh對(duì)稱性的八面體籠形團(tuán)簇,直徑約8.44 ?;這些氮化硼籠形團(tuán)簇的B—N 鍵長為1.44—1.46 ?,平均每個(gè)N 原子失去的電荷數(shù)為1.14e—1.45e.為了探究團(tuán)簇的元素組分效應(yīng),還考察了與B12N12具有相同結(jié)構(gòu)的B12P12,Al12N12和Ga12N12團(tuán)簇,籠子直徑在5.96—6.35 ?,它們具有比氮化硼團(tuán)簇更大的鍵長:B—P,Al—N和Ga—N 平均鍵長分別為1.91 ?,1.82 ?和1.88 ?,平均每個(gè)P,N 原子失去的電荷數(shù)為0.16e,1.93e和1.76e.另外,6 種團(tuán)簇的鍵級(jí)在0.62—1.91,表明原子間形成較強(qiáng)的化學(xué)鍵.頻率計(jì)算得到這些化合物團(tuán)簇的簡正模式,其最小頻率范圍為104.80—323.19 cm–1,沒有虛頻,進(jìn)一步證明了這些籠型結(jié)構(gòu)的穩(wěn)定性.

    圖16 種類富勒烯團(tuán)簇的基態(tài)結(jié)構(gòu),其中B,N,P,Al,Ga 原子的顏色分別用淺粉色、藍(lán)色、淺紫色、紫粉色和棕色表示Fig.1.Ground-state structures of six fullerene-like clusters.The B,N,P,Al,and Ga atoms are shown in pink,blue,purple,magenta,and brown colors,respectively.

    光激發(fā)導(dǎo)致的結(jié)構(gòu)變化對(duì)單重態(tài)和三重態(tài)之間的轉(zhuǎn)換有重要影響[49?51],因此,我們仔細(xì)表征了團(tuán)簇在S1態(tài)和T1態(tài)下的結(jié)構(gòu).如補(bǔ)充材料圖S1(online)所示,S1態(tài)相對(duì)基態(tài)S0結(jié)構(gòu)的最小均方根偏差(RMSD)范圍為0.04—0.11 ?,即籠型團(tuán)簇的結(jié)構(gòu)在單重激發(fā)態(tài)下只經(jīng)歷了微小的變化,證明了這些類富勒烯的高穩(wěn)定性.另外,T1態(tài)與S1態(tài)及基態(tài)的結(jié)構(gòu)非常相似,T1態(tài)相對(duì)基態(tài)S1結(jié)構(gòu)的RMSD 范圍為0.06—0.13 ?,這種激發(fā)態(tài)的微小形變可能導(dǎo)致S1態(tài)和T1態(tài)之間的重組能較小,使得T1態(tài)到S1態(tài)的反系間竄越速率較低,不利于TADF 過程,后續(xù)會(huì)詳細(xì)討論這個(gè)問題.

    表16 種類富勒烯團(tuán)簇的基態(tài)性質(zhì),包括對(duì)稱性(Sym.)、直徑(D)、平均鍵長(d)、最小頻率(υmin)、原子間電荷轉(zhuǎn)移(CT)、鍵級(jí)(BO)、HOMO-LUMO 能隙()以及吸收波長(λabs)Table 1. Ground-state properties of six fullerene-like clusters,including the symmetry (Sym.),diameter (D),average bond lengths (d),minimum frequency (υmin),Mulliken charge transfer (CT),average wiberg bond order (BO),HOMO-LUMO gap(),and absorption wavelength (λabs).

    表16 種類富勒烯團(tuán)簇的基態(tài)性質(zhì),包括對(duì)稱性(Sym.)、直徑(D)、平均鍵長(d)、最小頻率(υmin)、原子間電荷轉(zhuǎn)移(CT)、鍵級(jí)(BO)、HOMO-LUMO 能隙()以及吸收波長(λabs)Table 1. Ground-state properties of six fullerene-like clusters,including the symmetry (Sym.),diameter (D),average bond lengths (d),minimum frequency (υmin),Mulliken charge transfer (CT),average wiberg bond order (BO),HOMO-LUMO gap(),and absorption wavelength (λabs).

    3.2 類富勒烯團(tuán)簇的電子結(jié)構(gòu)、吸收光譜和發(fā)射光譜

    類富勒烯團(tuán)簇的電子態(tài)密度如圖2 所示,由HSE06 雜化泛函計(jì)算的B12N12,B24N24和B36N36的HOMO-LUMO 能隙分別是6.54 eV,6.11 eV 和6.20 eV,隨著籠子尺寸的增加能隙變化不大,與固體六方氮化硼和閃鋅礦氮化硼的能隙相近(5.97 eV和6.10 eV)[43,44].因此,這些團(tuán)簇可能與固體六方氮化硼和閃鋅礦氮化硼相似,主要吸收紫外光.而替換元素后,B12P12,Al12N12和Ga12N12的HOMOLUMO 能隙顯著變小,分別為3.30 eV,3.62 eV 和2.83 eV,能夠吸收可見光.這些團(tuán)簇的HOMO 能級(jí)主要分布在–6.42 eV 至–7.89 eV,LUMO 能級(jí)主要分布在–1.34 eV 至–3.72 eV.電子態(tài)密度隨元素的分布在補(bǔ)充材料圖S2 (online)中給出,結(jié)果顯示除了B12P12的其余5 種團(tuán)簇的占據(jù)能級(jí)主要由N 元素的電子態(tài)所貢獻(xiàn),而非占據(jù)能級(jí)主要由B,Al 和Ga 元素的電子態(tài)貢獻(xiàn);而對(duì)于B12P12團(tuán)簇,P 元素對(duì)占據(jù)能級(jí)和非占據(jù)能級(jí)都有較大貢獻(xiàn).6 種團(tuán)簇的吸收光譜如圖3(a)所示,3 個(gè)氮化硼籠型團(tuán)簇的吸收峰在189.60—202.97 nm 范圍,吸收紫外光;B12P12,Al12N12和Ga12N12的吸收光譜的波長范圍為375.75—438.16 nm,主要吸收紫色光.由此可見,氮化硼類富勒烯團(tuán)簇的能隙隨團(tuán)簇尺寸變化不明顯,而元素替換可以有效調(diào)控HOMO-LUMO 能隙,從而改變光吸收性質(zhì).

    圖2 類富勒烯團(tuán)簇的電子態(tài)密度.填充和未填充區(qū)域分別代表占據(jù)態(tài)和未占據(jù)態(tài),數(shù)字給出HOMO-LUMO 能隙,插圖為前線軌道分布,等值面設(shè)為0.04 a.u.Fig.2.Density of states (DOS) of fullerene-like clusters in their ground states.The occupied states are shadowed.The HOMO-LUMO gap is given for each cluster.The insets display the HOMO and LUMO wavefunctions,with an isosurface value of 0.04 a.u.

    圖3 類富勒烯團(tuán)簇的(a)吸收光譜和(b)發(fā)射光譜Fig.3.(a) Absorption spectra and (b) emission spectra of fullerene-like clusters.

    團(tuán)簇激發(fā)態(tài)下的電子態(tài)密度如補(bǔ)充材料圖S3(online)所示,單重激發(fā)態(tài)的HOMO-LUMO 能隙相對(duì)基態(tài)的能隙顯著減小,在2.08—3.69 eV 范圍,其中HOMO 能級(jí)主要分布在–5.83 eV 至–7.07 eV,LUMO 能級(jí)主要分布在–2.83 eV 至–3.85 eV.三重激發(fā)態(tài)的軌道分布與S1態(tài)相似,HOMO-LUMO能隙為1.63—3.70 eV,其中HOMO 能級(jí)主要分布在–5.85 eV 至–7.07 eV,LUMO 能級(jí)主要分布在–3.42 eV 至–4.22 eV.圖3(b)給出6 種團(tuán)簇的熒光發(fā)射光譜,根據(jù)卡莎規(guī)則(Kasha’s rule)[52],對(duì)于多重態(tài)的分子,光子僅能由最低激發(fā)態(tài)以熒光或者磷光的形式發(fā)射,發(fā)射光的波長和激發(fā)光的波長無關(guān).因此,這里我們只考慮S1態(tài)至S0態(tài)對(duì)應(yīng)的熒光發(fā)射光譜.根據(jù)計(jì)算的發(fā)射波長,由高斯展寬得到熒光發(fā)射譜,不同于實(shí)驗(yàn)測(cè)量的熒光譜,其半峰寬通常與樣品純度相關(guān)[53].總的來說,這些類富勒烯的熒光發(fā)射峰都在可見光范圍,其中B24N24發(fā)紫色熒光(405.36 nm),B12N12、B36N36和Al12N12發(fā)藍(lán)色熒光(488.38 nm,462.40 nm 和455.83 nm),Ga12N12發(fā)橙色熒光(644.47 nm),B12P12發(fā)紅色熒光(706.93 nm).相比之下,實(shí)驗(yàn)觀測(cè)C60的熒光發(fā)射峰在655—760 nm[20],C70的熒光發(fā)射峰在650—725 nm 范圍[16],它們主要發(fā)紅色熒光.考慮到實(shí)驗(yàn)上通常將富勒烯溶于正己烷溶劑來測(cè)量熒光譜[54,55],于是我們考察了溶劑效應(yīng)對(duì)類富勒烯團(tuán)簇發(fā)光性質(zhì)的影響,選擇正己烷作為溶劑,6 種團(tuán)簇的熒光譜在補(bǔ)充材料圖S4 (online)中給出,計(jì)算結(jié)果顯示溶劑效應(yīng)對(duì)發(fā)光峰位影響較小,真空和溶劑中發(fā)光波長的差別在0.15—24.24 nm 范圍,各團(tuán)簇的發(fā)光顏色基本保持不變.

    進(jìn)一步地,圖4 給出了類富勒烯團(tuán)簇的發(fā)光波長、基態(tài)能隙、S1和T1態(tài)能隙的分布.由圖4(a)可以看出,激發(fā)態(tài)S1和T1的HOMO-LUMO 能隙相近,且均比S0態(tài)的HOMO-LUMO 能隙值小,是由于最高占據(jù)態(tài)上的電子被激發(fā)到未占據(jù)態(tài)之后,使得原來最高占據(jù)態(tài)的能級(jí)位置升高、未占據(jù)態(tài)能級(jí)降低.圖4(b)給出了這些類富勒烯團(tuán)簇的S1態(tài)的HOMO-LUMO 能隙值與熒光的發(fā)射波長的關(guān)聯(lián).總的來說,3 種氮化硼團(tuán)簇(B12N12,B24N24和B36N36)隨尺寸增加,S1態(tài)能隙增大、熒光波長發(fā)生藍(lán)移;對(duì)于具有相同尺寸而不同元素組成的B12N12,Al12N12,Ga12N12和B12P12團(tuán)簇,S1態(tài)能隙隨著同族元素半徑的增大而減小,熒光波長發(fā)生紅移.由此可見,這些化合物團(tuán)簇具有更大的空間(尺寸和元素組成)來調(diào)控能隙和發(fā)光性質(zhì),相比于碳富勒烯,我們考慮的6 種類富勒烯團(tuán)簇在紫外或可見光激發(fā)下,能夠發(fā)射紅色至紫色區(qū)間的熒光,展現(xiàn)出更豐富的熒光顏色,是潛在的發(fā)光材料.鑒于它們良好的穩(wěn)定性、無毒性,有可能用于顯示、溫度或分子傳感器、醫(yī)學(xué)上的生物成像和標(biāo)記等領(lǐng)域[56].

    圖4 (a)類富勒烯團(tuán)簇的基態(tài)()、激發(fā)單重態(tài)()和三重態(tài)的HOMO-LUMO 能隙();(b)發(fā)射波長與 的關(guān)系圖Fig.4.(a) Histogram of HOMO-LUMO gaps of the ground state (),excited singlet state () and triplet state for the fullerene-like clusters ();(b) relationship between emission wavelength and .

    3.3 類富勒烯團(tuán)簇的TADF 性能

    在熱激活條件下,S1→T1態(tài)的系間竄越(intersystem crossing,ISC)和T1→S1態(tài)的反系間竄越(reverse intersystem crossing,RISC)是兩個(gè)重要的非輻射過程,能夠?qū)⑷貞B(tài)激子轉(zhuǎn)化為單重態(tài)激子,從而提高S1→S0輻射躍遷的內(nèi)量子效率[46].反系間竄越速率K主要由S1態(tài)和T1態(tài)的能量差?EST、自旋軌道耦合及S1態(tài)的重組能λS決定,可以通過半經(jīng)驗(yàn)的Marcus 速率公式來計(jì)算[57?60]:

    式中:kB是玻爾茲曼常數(shù);T是溫度(這里設(shè)為300 K);V為自旋軌道耦合常數(shù),可通過使用具有有效電荷近似的Breit-Pauli 自旋-軌道哈密頓量的PySOC 程序計(jì)算得到[61];λS為重組能,定義為單重激發(fā)態(tài)下分子結(jié)構(gòu)由S1態(tài)構(gòu)型變化至T1態(tài)構(gòu)型所引起的能量弛豫,

    系間竄越速率K′的計(jì)算公式與(3)式形式相似,其中S1態(tài)和T1態(tài)的能量差定義為?EST=ET1–ES1[45],重組能(λT)定義為三重激發(fā)態(tài)下分子結(jié)構(gòu)由T1態(tài)構(gòu)型變?yōu)镾1態(tài)構(gòu)型所引起的能量弛豫,即λT=ET1(S1-geometry)–ET1(T1-geometry).根據(jù)上述公式,系間竄越和反系間竄越速率隨?EST的增大而減小,隨V的增大而增大,重組能較小時(shí)隨λS或λT的增大而增大.

    我們計(jì)算了6 種類富勒烯團(tuán)簇的系間竄越和反系間竄越過程的相關(guān)物理參數(shù),總結(jié)在表2 中.可以看出,這些籠型團(tuán)簇均具有較小的?EST值,在0.14—0.29 eV 范圍,低于C60和C70的?EST(計(jì)算值分別為0.37 eV 和0.36 eV[62],實(shí)驗(yàn)測(cè)量值分別為0.36 eV 和0.30 eV)[19,21],有利于反系間竄越和TADF 過程.這些團(tuán)簇的自旋-軌道耦合常數(shù)在2.30×10–3—9.86 cm–1范圍,其中Ga12N12的V值最大.另外,這些類富勒烯團(tuán)簇具有很高的穩(wěn)定性,其激發(fā)態(tài)結(jié)構(gòu)相對(duì)基態(tài)變化非常微小,且不隨自旋多重度有明顯變化,導(dǎo)致S1態(tài)和T1態(tài)之間的重組能較小,在3.00×10–4—0.15 eV 范圍.基于此,由(3)式計(jì)算所得這些團(tuán)簇體系中Al12N12的系間竄越速率最大為2.09×107s–1,與C60(1.34×107s–1)相近,低于C70的系間竄越速率(2.04×108s–1)[62],Ga12N12的反系間竄越速率最大為5.9×105s–1,表明Al12N12和Ga12N12可能具有TADF 效應(yīng).

    3.4 類富勒烯團(tuán)簇激發(fā)特性分析

    氮化硼籠型團(tuán)簇的低?EST值可以通過它們的前線軌道分布來理解.圖2 為6 種類富勒烯團(tuán)簇的基態(tài)HOMO 和LUMO 的電子云分布,明顯地,HOMO 和LUMO 的電子云分別分布于不同元素上.對(duì)于B12N12,B24N24和B36N36團(tuán)簇,它們基態(tài)的HOMO 軌道均勻地分布在N 原子上,而LUMO則均勻地分布在B 原子上;類似地,對(duì)于B12P12,Al12N12和Ga12N12團(tuán)簇,它們HOMO 軌道主要分布在P 或N 原子上,而LUMO 主要分布在B,Al,Ga 原子上.這種獨(dú)特的化合物籠型團(tuán)簇結(jié)構(gòu)將HOMO 和LUMO 從空間上分離,能夠有效地降低HOMO-LUMO 波函數(shù)的重疊. 根據(jù)?EST和HOMO-LUMO 波函數(shù)的重疊積分公式[63]:

    式中J為電子之間的交換能,由HOMO-LUMO的重疊積分決定.HOMO 和LUMO 在空間上的分離使得重疊積分較小、?EST降低,從而使反系間竄越和TADF 過程更容易發(fā)生.這個(gè)策略在設(shè)計(jì)有機(jī)分子TADF 材料中已被廣泛采用,通過引入電子供體和受體單元來分離HOMO 和LUMO 軌道,供體和受體單元的相對(duì)扭轉(zhuǎn)角度對(duì)軌道重疊程度和?EST有一定的影響[11,14].而本文考慮的類富勒烯團(tuán)簇因具有特殊的化合物結(jié)構(gòu)也能使HOMO和LUMO 在空間上分離,并且由于籠狀結(jié)構(gòu)的穩(wěn)定性,軌道重疊程度受結(jié)構(gòu)擾動(dòng)的影響不大.

    為了進(jìn)一步表征類富勒烯團(tuán)簇的激發(fā)態(tài)特性,我們還利用Multiwfn 程序包分析了激發(fā)態(tài)下空穴和電子的密度分布[64,65].如圖5 所示,對(duì)于B12N12,B24N24,B36N36,Al12N12和Ga12N12團(tuán) 簇,這5 個(gè)團(tuán)簇的S1態(tài)的空穴和電子分別出現(xiàn)在不同的元素上,即在B12N12,B24N24和B36N36團(tuán)簇中,電子主要分布在B 原子上,空穴則主要分布在N 原子上;對(duì)于Al12N12和Ga12N12團(tuán)簇,電子主要分布在Al和Ga 原子上,而空穴則主要分布在N 原子上.激發(fā)態(tài)S1上電子和空穴的分布在空間上分離,說明這些團(tuán)簇具有明顯的電子轉(zhuǎn)移(charge-transfer)激發(fā)特性[66].而在B12P12團(tuán)簇,電子和空穴均勻地分布于整個(gè)團(tuán)簇,沒有明顯的空間分離特征,即表現(xiàn)出局域激發(fā)特征(locally excited-state).為了進(jìn)一步證實(shí)這類團(tuán)簇的激發(fā)特性,我們對(duì)S1態(tài)的空穴和電子在全空間中的分布特征進(jìn)行了定量描述,計(jì)算了S1和T1態(tài)下空穴和電子波函數(shù)的重疊程度(Sr),空穴和電子質(zhì)心之間的距離指數(shù)(Δr),以及空穴和電子的分離程度指數(shù)(t),具體定義見補(bǔ)充材料 (online).與圖5 一致,B12P12激發(fā)態(tài)的空穴和電子波函數(shù)的重疊程度較大(Sr >0.5 a.u.),空穴和電子質(zhì)心之間的距離較小(Δr<0.1 ?),且空穴和電子的分離程度指數(shù)為負(fù)(t<0 表示空穴和電子沒有顯著分離),屬于局域激發(fā)特征;而其他5 種團(tuán)簇激發(fā)態(tài)的空穴和電子波函數(shù)的重疊程度只有0.35—0.39 a.u.,空穴和電子質(zhì)心間距較大,在1.51—2.33 ?范圍,空穴和電子的分離程度指數(shù)t>0,顯示出電荷轉(zhuǎn)移特征.

    圖5 類富勒烯團(tuán)簇S1 態(tài)的空穴(黃色)和電子(粉色)密度分布.數(shù)字給出空穴和電子波函數(shù)的重疊值,等值面設(shè)為0.04 a.u.Fig.5.Electron (magenta) and hole (yellow) density distributions for the excited singlet state of fullerene-like clusters.The overlap between hole and electron wavefunctions are given for each cluster,with an isosurface value of 0.04 a.u.

    表26 種類富勒烯團(tuán)簇的激發(fā)態(tài)性質(zhì),包括S1 態(tài)的發(fā)射波長(λex)、HOMO-LUMO 能隙()、激發(fā)能(ES1)、S1 態(tài)和T1 態(tài)的能量差(?EST)、自旋-軌道耦合常數(shù)(V)、以及重組能(λS,λT)Table 2. Excited-state properties of six fullerene-like clusters,including emission wavelength (λex),HOMO-LUMO gap(),and emission energy of S1 state (ES1),energy difference between the S1 and T1 state (?EST),spin-orbit coupling constant (V),and recombination energy of S1 state and T1 state (λS,λT).

    表26 種類富勒烯團(tuán)簇的激發(fā)態(tài)性質(zhì),包括S1 態(tài)的發(fā)射波長(λex)、HOMO-LUMO 能隙()、激發(fā)能(ES1)、S1 態(tài)和T1 態(tài)的能量差(?EST)、自旋-軌道耦合常數(shù)(V)、以及重組能(λS,λT)Table 2. Excited-state properties of six fullerene-like clusters,including emission wavelength (λex),HOMO-LUMO gap(),and emission energy of S1 state (ES1),energy difference between the S1 and T1 state (?EST),spin-orbit coupling constant (V),and recombination energy of S1 state and T1 state (λS,λT).

    4 結(jié)論

    本文利用第一性原理計(jì)算研究了6 種類富勒烯團(tuán)簇的發(fā)光性質(zhì),包括實(shí)驗(yàn)報(bào)道的具有不同尺寸的氮化硼團(tuán)簇B12N12,B24N24,B36N36,以及對(duì)B12N12進(jìn)行元素替換所得的B12P12,Al12N12和Ga12N12團(tuán)簇.在HSE06 雜化泛函描述下,這些團(tuán)簇的HOMO-LUMO 能隙在2.83—6.54 eV 范圍,吸收峰位置主要在紫外光區(qū)間,而熒光發(fā)射波長均在可見光范圍405.36—706.93 nm,包括紅光(B12P12)、橙光(Ga12N12)、藍(lán)光(B12N12,B36N36和Al12N12)、紫光(B24N24).這些籠型團(tuán)簇的發(fā)光峰位置與激發(fā)態(tài)HOMO-LUMO 能隙有一定的關(guān)聯(lián):3 種氮化硼團(tuán)簇的能隙隨尺寸增加而增大、發(fā)射波長發(fā)生藍(lán)移;對(duì)于具有相同尺寸而不同元素組成的B12P12,Al12N12,Ga12N12和B12N12團(tuán)簇,能隙隨同族元素半徑的增大而減小,發(fā)射波長發(fā)生紅移.另外,這些化合物團(tuán)簇的特殊結(jié)構(gòu)使得HOMO 和LUMO軌道分布在不同元素上,導(dǎo)致團(tuán)簇的單重激發(fā)態(tài)和三重激發(fā)態(tài)之間的能量差較小,有利于反系間竄越和熱激活延遲熒光過程.這些理論結(jié)果揭示了氮化硼等類富勒烯團(tuán)簇的激發(fā)態(tài)特性,為精準(zhǔn)調(diào)控其發(fā)光性質(zhì)提供了重要指導(dǎo).鑒于它們的可制備性、高穩(wěn)定性、無毒性、較大的能隙和可見光發(fā)射熒光光譜,這些類富勒烯團(tuán)簇是潛在的發(fā)光材料,可能用于顯示、生物檢測(cè)、醫(yī)學(xué)等領(lǐng)域.

    感謝大連理工大學(xué)超級(jí)計(jì)算中心和鴻之微提供的計(jì)算支持.

    猜你喜歡
    能隙富勒烯氮化硼
    高效硫硒化銻薄膜太陽電池中的漸變能隙結(jié)構(gòu)*
    Bogoliubov-Tolmachev-Shirkov模型臨界溫度和能隙解的數(shù)值方法
    鉆石級(jí)抗衰 諾貝爾獎(jiǎng)光環(huán)揭開“富勒烯”的神秘面紗
    中國化妝品(2019年4期)2019-11-20 01:47:53
    富勒烯基聚合物制備研究的進(jìn)展
    聚晶立方氮化硼復(fù)合片電火花線切割高效切割研究
    新型富勒烯的合成
    三甲基硅雙苯乙炔取代噻咯的合成及聚集態(tài)誘導(dǎo)發(fā)光
    5d過渡金屬原子吸附氮化硼納米管的第一性原理計(jì)算
    基于曲率和電子結(jié)構(gòu)的摻雜C50和C70富勒烯的穩(wěn)定性研究
    鈷摻雜氮化硼納米管吸附氯酚類污染物的理論研究
    99久久综合免费| 国产片特级美女逼逼视频| 两个人看的免费小视频| av又黄又爽大尺度在线免费看| 97在线视频观看| 午夜老司机福利剧场| 免费看av在线观看网站| 国产淫语在线视频| 国产av精品麻豆| 如日韩欧美国产精品一区二区三区| 欧美中文综合在线视频| 狂野欧美激情性bbbbbb| 国产又色又爽无遮挡免| freevideosex欧美| 青春草亚洲视频在线观看| 女性被躁到高潮视频| 精品人妻一区二区三区麻豆| 老汉色av国产亚洲站长工具| 国产精品一国产av| 看免费av毛片| 国产精品久久久久久精品古装| 赤兔流量卡办理| 亚洲av日韩在线播放| 日日撸夜夜添| 丰满少妇做爰视频| 天天躁狠狠躁夜夜躁狠狠躁| 熟女电影av网| 亚洲一级一片aⅴ在线观看| 热re99久久国产66热| 老汉色∧v一级毛片| 久久亚洲国产成人精品v| 精品人妻偷拍中文字幕| 国产片内射在线| 亚洲欧美中文字幕日韩二区| av视频免费观看在线观看| 亚洲一码二码三码区别大吗| 大香蕉久久成人网| 伦理电影免费视频| 麻豆精品久久久久久蜜桃| 精品人妻在线不人妻| 成年美女黄网站色视频大全免费| 亚洲综合色网址| 我要看黄色一级片免费的| 久久久久久人妻| 亚洲男人天堂网一区| 交换朋友夫妻互换小说| 搡老乐熟女国产| 亚洲欧美一区二区三区国产| a级片在线免费高清观看视频| 成人毛片a级毛片在线播放| 精品午夜福利在线看| 9191精品国产免费久久| 巨乳人妻的诱惑在线观看| 制服丝袜香蕉在线| 精品午夜福利在线看| 不卡av一区二区三区| 国产野战对白在线观看| 满18在线观看网站| av一本久久久久| 一级a爱视频在线免费观看| 秋霞伦理黄片| 大香蕉久久成人网| 男男h啪啪无遮挡| 亚洲国产成人一精品久久久| 久久国产精品大桥未久av| 亚洲av成人精品一二三区| 黄色毛片三级朝国网站| 久久精品国产鲁丝片午夜精品| 女性生殖器流出的白浆| 色吧在线观看| 超碰成人久久| 80岁老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 久久ye,这里只有精品| 少妇人妻精品综合一区二区| 精品少妇黑人巨大在线播放| 国产精品秋霞免费鲁丝片| 91精品伊人久久大香线蕉| 精品99又大又爽又粗少妇毛片| 看免费成人av毛片| 亚洲av男天堂| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 麻豆av在线久日| 中文字幕最新亚洲高清| 亚洲国产成人一精品久久久| 午夜福利在线观看免费完整高清在| 亚洲成国产人片在线观看| 婷婷色综合www| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品天堂在线| 久久 成人 亚洲| 高清在线视频一区二区三区| 边亲边吃奶的免费视频| 青青草视频在线视频观看| 国产无遮挡羞羞视频在线观看| 欧美日本中文国产一区发布| 色吧在线观看| 97人妻天天添夜夜摸| 考比视频在线观看| 99香蕉大伊视频| 精品少妇内射三级| 97在线人人人人妻| 99国产精品免费福利视频| 肉色欧美久久久久久久蜜桃| 亚洲欧美一区二区三区黑人 | 国产人伦9x9x在线观看 | 久久亚洲国产成人精品v| 国产精品av久久久久免费| 黄频高清免费视频| 日本午夜av视频| 搡老乐熟女国产| 在线 av 中文字幕| 18禁动态无遮挡网站| 宅男免费午夜| av网站免费在线观看视频| 亚洲伊人久久精品综合| 久久热在线av| 少妇人妻精品综合一区二区| 日本91视频免费播放| 一区二区日韩欧美中文字幕| videosex国产| 男女边摸边吃奶| 日韩中文字幕视频在线看片| 国产成人a∨麻豆精品| 久久久久久久大尺度免费视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日本色播在线视频| 国产精品欧美亚洲77777| www日本在线高清视频| 精品亚洲成国产av| 免费高清在线观看日韩| 又大又黄又爽视频免费| 欧美精品av麻豆av| 最近2019中文字幕mv第一页| 精品亚洲成a人片在线观看| av有码第一页| 99热国产这里只有精品6| 日本av手机在线免费观看| 日韩 亚洲 欧美在线| 亚洲av在线观看美女高潮| 午夜福利乱码中文字幕| 免费女性裸体啪啪无遮挡网站| 久久久精品国产亚洲av高清涩受| 精品第一国产精品| 日本av手机在线免费观看| 久久国内精品自在自线图片| 久久久久久人妻| 精品一区二区三卡| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品999| 性高湖久久久久久久久免费观看| 夫妻性生交免费视频一级片| 欧美日韩一区二区视频在线观看视频在线| 男女边摸边吃奶| 26uuu在线亚洲综合色| 久久人人爽av亚洲精品天堂| av不卡在线播放| 丝瓜视频免费看黄片| 国产探花极品一区二区| 一级,二级,三级黄色视频| videos熟女内射| 18禁国产床啪视频网站| 男的添女的下面高潮视频| 女人精品久久久久毛片| 啦啦啦在线免费观看视频4| 天堂俺去俺来也www色官网| 热99久久久久精品小说推荐| 国产成人精品久久二区二区91 | av有码第一页| 国产亚洲最大av| 欧美 日韩 精品 国产| 中文字幕人妻丝袜制服| 日本爱情动作片www.在线观看| 亚洲av在线观看美女高潮| av网站免费在线观看视频| 午夜免费鲁丝| 成年动漫av网址| 观看av在线不卡| 又黄又粗又硬又大视频| 国产精品二区激情视频| 亚洲内射少妇av| 日本欧美国产在线视频| 巨乳人妻的诱惑在线观看| 免费人妻精品一区二区三区视频| 国产片内射在线| 亚洲三区欧美一区| 成年美女黄网站色视频大全免费| 日日爽夜夜爽网站| 欧美97在线视频| 在线观看免费日韩欧美大片| 免费在线观看视频国产中文字幕亚洲 | 精品人妻熟女毛片av久久网站| 十分钟在线观看高清视频www| 少妇 在线观看| 国产成人精品福利久久| 中文字幕最新亚洲高清| 午夜影院在线不卡| 制服诱惑二区| 精品少妇久久久久久888优播| 我的亚洲天堂| 亚洲av.av天堂| av国产久精品久网站免费入址| 欧美精品亚洲一区二区| 午夜免费观看性视频| 国产男女超爽视频在线观看| 免费女性裸体啪啪无遮挡网站| 1024视频免费在线观看| av网站在线播放免费| a 毛片基地| 一区二区日韩欧美中文字幕| 高清黄色对白视频在线免费看| 日韩伦理黄色片| 亚洲精品国产一区二区精华液| 亚洲,一卡二卡三卡| 热re99久久国产66热| 午夜精品国产一区二区电影| 久久亚洲国产成人精品v| 国产黄色免费在线视频| av一本久久久久| 欧美日韩av久久| 欧美在线黄色| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av成人精品一二三区| av在线播放精品| 亚洲成人一二三区av| 蜜桃在线观看..| 老鸭窝网址在线观看| 国产成人精品婷婷| 国产精品嫩草影院av在线观看| 亚洲国产精品国产精品| 国产福利在线免费观看视频| 校园人妻丝袜中文字幕| 一区福利在线观看| 少妇人妻久久综合中文| 在线免费观看不下载黄p国产| 国产免费现黄频在线看| 精品第一国产精品| 精品人妻一区二区三区麻豆| 午夜av观看不卡| 99久久综合免费| 热99久久久久精品小说推荐| 最近最新中文字幕免费大全7| 亚洲精品成人av观看孕妇| 亚洲中文av在线| 国产成人精品婷婷| 精品人妻一区二区三区麻豆| 最黄视频免费看| 少妇 在线观看| 国产精品一区二区在线观看99| 国产亚洲午夜精品一区二区久久| 国产乱人偷精品视频| 亚洲精品视频女| 国产在视频线精品| 免费大片黄手机在线观看| 成人国语在线视频| 国产精品久久久久久久久免| 亚洲激情五月婷婷啪啪| 久久久久网色| 人体艺术视频欧美日本| 日本av手机在线免费观看| 久久狼人影院| 国产成人免费观看mmmm| 欧美成人午夜免费资源| 一区在线观看完整版| 在线 av 中文字幕| 成人亚洲欧美一区二区av| 久久久久人妻精品一区果冻| 80岁老熟妇乱子伦牲交| 国产一区亚洲一区在线观看| 三上悠亚av全集在线观看| 欧美精品人与动牲交sv欧美| 一本久久精品| 最近中文字幕2019免费版| 亚洲欧洲精品一区二区精品久久久 | 青春草亚洲视频在线观看| 欧美av亚洲av综合av国产av | 亚洲国产欧美网| 久久毛片免费看一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产日韩欧美视频二区| 天天操日日干夜夜撸| 性色avwww在线观看| 欧美激情 高清一区二区三区| 一本—道久久a久久精品蜜桃钙片| 婷婷色av中文字幕| 18在线观看网站| a级毛片黄视频| 国产白丝娇喘喷水9色精品| 国产亚洲午夜精品一区二区久久| 一个人免费看片子| 亚洲一级一片aⅴ在线观看| 久久久久久人人人人人| 成人手机av| 色婷婷久久久亚洲欧美| 中文字幕另类日韩欧美亚洲嫩草| 嫩草影院入口| 国产探花极品一区二区| 天堂8中文在线网| 午夜福利乱码中文字幕| 久久久久国产一级毛片高清牌| 女的被弄到高潮叫床怎么办| 日产精品乱码卡一卡2卡三| 午夜91福利影院| 精品人妻熟女毛片av久久网站| 亚洲精品一区蜜桃| av在线观看视频网站免费| 欧美日韩av久久| 日本免费在线观看一区| 天天躁日日躁夜夜躁夜夜| 五月开心婷婷网| 成人影院久久| 80岁老熟妇乱子伦牲交| 男女边吃奶边做爰视频| 久久精品aⅴ一区二区三区四区 | 一边摸一边做爽爽视频免费| 老鸭窝网址在线观看| 老司机亚洲免费影院| 久久久久国产网址| 我的亚洲天堂| 国产不卡av网站在线观看| 中文精品一卡2卡3卡4更新| 熟妇人妻不卡中文字幕| 欧美日韩视频精品一区| 丝袜脚勾引网站| 久久精品国产亚洲av高清一级| 91精品国产国语对白视频| 人妻 亚洲 视频| 一本—道久久a久久精品蜜桃钙片| 伊人久久大香线蕉亚洲五| 黄片无遮挡物在线观看| 啦啦啦啦在线视频资源| 国产精品国产三级专区第一集| 不卡av一区二区三区| 精品国产一区二区三区久久久樱花| xxxhd国产人妻xxx| a级片在线免费高清观看视频| 午夜激情久久久久久久| 青春草视频在线免费观看| 色婷婷av一区二区三区视频| 中文字幕亚洲精品专区| 国产欧美日韩综合在线一区二区| 丝袜在线中文字幕| 制服诱惑二区| 国产精品久久久久久久久免| 久久精品久久久久久噜噜老黄| 久久久精品区二区三区| 久久免费观看电影| 下体分泌物呈黄色| 欧美日韩一区二区视频在线观看视频在线| 午夜福利影视在线免费观看| 久久久久久久久久人人人人人人| 王馨瑶露胸无遮挡在线观看| 黑丝袜美女国产一区| 国产免费视频播放在线视频| 亚洲,欧美,日韩| 久久影院123| 男女无遮挡免费网站观看| 免费不卡的大黄色大毛片视频在线观看| 涩涩av久久男人的天堂| 亚洲美女视频黄频| 日日啪夜夜爽| 久久久久国产一级毛片高清牌| 久久久a久久爽久久v久久| 国产成人91sexporn| 蜜桃国产av成人99| 欧美日韩视频精品一区| 黄色 视频免费看| 制服丝袜香蕉在线| 激情五月婷婷亚洲| 纯流量卡能插随身wifi吗| 大陆偷拍与自拍| 爱豆传媒免费全集在线观看| 97人妻天天添夜夜摸| 90打野战视频偷拍视频| 多毛熟女@视频| 精品人妻一区二区三区麻豆| 制服丝袜香蕉在线| 啦啦啦视频在线资源免费观看| 国产野战对白在线观看| 两性夫妻黄色片| 日本欧美国产在线视频| 国产成人精品在线电影| 国产人伦9x9x在线观看 | 满18在线观看网站| 青春草亚洲视频在线观看| 久久人人97超碰香蕉20202| av.在线天堂| www.熟女人妻精品国产| 在线亚洲精品国产二区图片欧美| 久久人妻熟女aⅴ| 精品一区在线观看国产| 韩国精品一区二区三区| 亚洲激情五月婷婷啪啪| 在线观看免费高清a一片| 亚洲国产看品久久| 日本午夜av视频| 老司机影院成人| 在线精品无人区一区二区三| 十分钟在线观看高清视频www| 韩国高清视频一区二区三区| 国产免费福利视频在线观看| 伦理电影免费视频| 丁香六月天网| 国产精品偷伦视频观看了| 欧美精品亚洲一区二区| 天天躁夜夜躁狠狠久久av| 久久久久网色| 欧美亚洲 丝袜 人妻 在线| 亚洲精品成人av观看孕妇| 亚洲激情五月婷婷啪啪| 秋霞在线观看毛片| 国产 一区精品| 飞空精品影院首页| 人人妻人人澡人人看| 精品99又大又爽又粗少妇毛片| 免费播放大片免费观看视频在线观看| 国产乱人偷精品视频| 七月丁香在线播放| 亚洲精品久久成人aⅴ小说| 大片电影免费在线观看免费| 亚洲精品久久成人aⅴ小说| 男女国产视频网站| 久久av网站| 啦啦啦中文免费视频观看日本| 国产欧美亚洲国产| 国产精品久久久久久av不卡| 黑丝袜美女国产一区| 精品久久久精品久久久| 男女免费视频国产| 水蜜桃什么品种好| 亚洲精品久久久久久婷婷小说| 韩国高清视频一区二区三区| 26uuu在线亚洲综合色| 欧美成人午夜免费资源| 少妇的丰满在线观看| 欧美av亚洲av综合av国产av | 五月天丁香电影| 天天操日日干夜夜撸| 97精品久久久久久久久久精品| 只有这里有精品99| 中文字幕制服av| 日韩一卡2卡3卡4卡2021年| 午夜福利视频在线观看免费| 爱豆传媒免费全集在线观看| 午夜免费观看性视频| 久久午夜综合久久蜜桃| 亚洲成av片中文字幕在线观看 | 亚洲,欧美精品.| 校园人妻丝袜中文字幕| 成人亚洲欧美一区二区av| 91精品国产国语对白视频| 美女xxoo啪啪120秒动态图| 黄片小视频在线播放| 午夜免费鲁丝| 日韩熟女老妇一区二区性免费视频| 夜夜骑夜夜射夜夜干| 亚洲成人av在线免费| 久久久久网色| 一级,二级,三级黄色视频| 性色av一级| 免费观看a级毛片全部| 久久久久久久大尺度免费视频| 精品久久久久久电影网| 午夜影院在线不卡| av天堂久久9| 久久鲁丝午夜福利片| 嫩草影院入口| 亚洲av日韩在线播放| 三级国产精品片| 日韩伦理黄色片| 国产一区有黄有色的免费视频| 成年女人毛片免费观看观看9 | 丝袜脚勾引网站| 国产精品国产三级专区第一集| 久久久久久久国产电影| 极品人妻少妇av视频| 亚洲国产精品国产精品| 夫妻午夜视频| 国产一区二区激情短视频 | 免费高清在线观看视频在线观看| 久久久久久久国产电影| 大码成人一级视频| 黄片小视频在线播放| 毛片一级片免费看久久久久| 黄频高清免费视频| 麻豆乱淫一区二区| 婷婷色综合大香蕉| 一个人免费看片子| 一级毛片电影观看| 久久久久国产网址| 不卡av一区二区三区| 这个男人来自地球电影免费观看 | 国产在线视频一区二区| 人妻人人澡人人爽人人| 老鸭窝网址在线观看| 精品国产超薄肉色丝袜足j| 亚洲伊人久久精品综合| 精品人妻在线不人妻| 亚洲精品av麻豆狂野| 叶爱在线成人免费视频播放| 国产黄色视频一区二区在线观看| 观看av在线不卡| 一级爰片在线观看| 中文字幕最新亚洲高清| 久久国产精品大桥未久av| a 毛片基地| 香蕉国产在线看| 久久久久久久国产电影| 精品午夜福利在线看| 久久久欧美国产精品| 黄色毛片三级朝国网站| 国产成人精品无人区| 妹子高潮喷水视频| 日韩大片免费观看网站| 两个人免费观看高清视频| 久久久久久久久免费视频了| 最近中文字幕2019免费版| 欧美精品亚洲一区二区| 搡女人真爽免费视频火全软件| 国产免费又黄又爽又色| 日本av手机在线免费观看| 久久国产精品大桥未久av| 少妇人妻精品综合一区二区| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 国产视频首页在线观看| 99精国产麻豆久久婷婷| 少妇熟女欧美另类| 亚洲国产av新网站| 亚洲人成网站在线观看播放| 久久午夜福利片| 欧美精品一区二区大全| 亚洲精品中文字幕在线视频| 爱豆传媒免费全集在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 秋霞伦理黄片| 啦啦啦视频在线资源免费观看| 边亲边吃奶的免费视频| 黄色 视频免费看| 又粗又硬又长又爽又黄的视频| 亚洲欧美成人精品一区二区| 天天躁日日躁夜夜躁夜夜| 这个男人来自地球电影免费观看 | 夫妻性生交免费视频一级片| 久久久久国产精品人妻一区二区| 最新中文字幕久久久久| 久久这里有精品视频免费| 久久精品久久精品一区二区三区| 肉色欧美久久久久久久蜜桃| 国产精品一二三区在线看| 女人精品久久久久毛片| 两个人免费观看高清视频| 亚洲国产欧美网| 电影成人av| 免费黄网站久久成人精品| 亚洲综合色网址| av有码第一页| 国产精品99久久99久久久不卡 | 一区二区三区乱码不卡18| 国产极品天堂在线| 黑人巨大精品欧美一区二区蜜桃| 丰满乱子伦码专区| 黄色配什么色好看| 国产精品嫩草影院av在线观看| 如日韩欧美国产精品一区二区三区| 欧美av亚洲av综合av国产av | 在线 av 中文字幕| 亚洲精品日韩在线中文字幕| 在线免费观看不下载黄p国产| 边亲边吃奶的免费视频| 好男人视频免费观看在线| 精品国产乱码久久久久久男人| 国产精品久久久久久精品电影小说| 视频在线观看一区二区三区| 在线观看www视频免费| 视频在线观看一区二区三区| 日韩一区二区三区影片| 免费高清在线观看视频在线观看| 精品少妇一区二区三区视频日本电影 | 侵犯人妻中文字幕一二三四区| 一区二区三区精品91| 亚洲成人一二三区av| 看非洲黑人一级黄片| 老女人水多毛片| 亚洲综合色网址| 亚洲精华国产精华液的使用体验| 人人妻人人添人人爽欧美一区卜| 国产在线视频一区二区| 99九九在线精品视频| 永久网站在线| 亚洲精品日本国产第一区| 亚洲激情五月婷婷啪啪| 韩国高清视频一区二区三区| 成人黄色视频免费在线看| 欧美日韩亚洲高清精品| 国产精品秋霞免费鲁丝片| 国产野战对白在线观看| av网站在线播放免费| 亚洲精华国产精华液的使用体验| 黄片播放在线免费| 午夜福利在线免费观看网站| 建设人人有责人人尽责人人享有的| 一本—道久久a久久精品蜜桃钙片| 不卡视频在线观看欧美| 欧美日韩av久久| 日韩一卡2卡3卡4卡2021年| 精品少妇黑人巨大在线播放| 下体分泌物呈黄色| 桃花免费在线播放| 熟女av电影| 国产精品嫩草影院av在线观看| 欧美人与性动交α欧美精品济南到 | 中文字幕人妻丝袜制服| 欧美成人精品欧美一级黄|