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

    探討多原子分子的Schr?dinger 方程組的嚴(yán)格解析解法——以血紅蛋白電離能的計(jì)算為例

    2012-07-02 00:50:36任愛娣王文信
    兵器裝備工程學(xué)報(bào) 2012年4期
    關(guān)鍵詞:電離能原子核原子

    任愛娣,王文信

    ( 海軍工程大學(xué)( 天津校區(qū)),天津 300450)

    眾所周知,許許多多科學(xué)家都研究過多原子分子Schr?dinger 方程組[1]

    其中:方程(1)是固定核的位置時(shí)電子運(yùn)動(dòng)方程;E 為所需要的能量;方程(2)是核運(yùn)動(dòng)方程;ε 為總能量;M 和N 分別代表原子核和電子的總數(shù);V 系靜電勢(shì)能。

    可見,方程組(1)和(2)的解是十分重要的。Bohr.N,Hartree.D.R ,F(xiàn)ock.B.A,Oppenheimer.J.R 和Born.M 等[1-5]的研究成果一直沿用到今天,極大地豐富了人類的知識(shí)寶庫。唐敖慶先生在文獻(xiàn)[1]中只是給出雙原子分子Schr?dinger 的近似解,雖然已經(jīng)給出雙原子分子、三原子分子和四原子分子的Schr?dinger 方程組的嚴(yán)格解析解[6],但是,這并不意味著能夠給出所有多原子分子Schr?dinger 方程組的嚴(yán)格解析解。那么對(duì)于任意一個(gè)多原子分子Schr?dinger 方程組究竟有沒有嚴(yán)格解析解呢? 至今沒見答案。為了能夠找到這個(gè)答案,筆者發(fā)展了我國(guó)著名科學(xué)家唐敖慶先生在文獻(xiàn)[1]中的方法,在超球坐標(biāo)系中,首先獲得了電子運(yùn)動(dòng)方程的嚴(yán)格解析解。從所得解中找到了原子核運(yùn)動(dòng)的所需的能量,自然就能夠求得原子核運(yùn)動(dòng)方程的嚴(yán)格解析解。又從所得解中找到總能量和電離能的計(jì)算公式,就能夠計(jì)算血紅蛋白分子的電離能,從而驗(yàn)證該探討是成功的。答案就自動(dòng)產(chǎn)生了。

    1 電子運(yùn)動(dòng)的Schr?dinger 方程的嚴(yán)格解析解

    對(duì)于電子運(yùn)動(dòng)的Schr?dinger 方程(1),用超球諧法[7-9]可以得到它的解析解

    其中:rij、r0ij和Rij分別表示原子核i 到電子j,電子i 到電子j 以及原子核i 到原子核j 之間的距離;Zj系原子序數(shù);N0=1,Nj=N0+Z1+…Zj;Ω 是球極角和超球角的集。

    為了能夠方便地應(yīng)用所得到的嚴(yán)格解析解(3),需要進(jìn)一步剖析它。首先看它的靜電勢(shì)能V。從式(5)可以看出,它是各個(gè)電子之間、各個(gè)原子核之間和原子核與電子之間所形成的靜電勢(shì)能之和,其中包括各個(gè)電子繞自己的原子核運(yùn)動(dòng)所產(chǎn)生的靜電勢(shì)能之和V*,即V*?V。顯然V*與核間距Rij的變化無關(guān),而V-V*則相反,它們的變化就意味著核間距Rij的變化。由此就可以把靜電勢(shì)能V 分解成

    其中V*與核間距Rij的變化無關(guān),而V-V*則相反。

    既然靜電勢(shì)能V 可以寫成2 項(xiàng)之和,它們?cè)诔蛑C中的展開式自然也能夠?qū)懗? 項(xiàng)之和,即

    由于V*與核間距Rij的變化無關(guān),所以它在超球諧中的展開式與核間距Rij的變化也無關(guān),而則相反,它的變化就意味著核間距Rij的變化。類似地,也能把式(6)一分為二:其一與核間距Rij無關(guān);其二則相反,即

    其中:

    能量公式改寫成式(10)后,其意義更加明。先看第2 項(xiàng),由于α*與V*有關(guān),與核間距Rij變化無關(guān),V*是各個(gè)電子繞自己的核運(yùn)動(dòng)所產(chǎn)生的靜電勢(shì)能之和,所以-α2*/2 所表示的能量與V*有關(guān),與核間距Rij變化無關(guān)。第1 項(xiàng)則相反,E*變化意味著核間距Rij變化,也就是說,E*變化決定著核間距Rij變化,可見,E*應(yīng)該是核間距Rij變化所需要的能量。

    E*是核間距Rij的變化所需要的能量,而核間距Rij的變化是由原子核運(yùn)動(dòng)決定的。那么原子核是如何運(yùn)動(dòng)呢? 根據(jù)熟知的Born-Oppenheimer 理論[1],由于原子核的質(zhì)量比電子大103~105倍,電子速度比原子核快得多,這就使得當(dāng)核間距任意微小運(yùn)動(dòng)時(shí),迅速運(yùn)動(dòng)的電子都能立即進(jìn)行調(diào)整,建立起與變幻后的核力場(chǎng)相應(yīng)的運(yùn)動(dòng)狀態(tài);原子核相對(duì)于電子來說,速度慢得多,好像是不運(yùn)動(dòng),其實(shí)也可能在運(yùn)動(dòng),不可能絕對(duì)不運(yùn)動(dòng)。既然原子核也在運(yùn)動(dòng)。就需要去解原子核運(yùn)動(dòng)的Schr?dinger 方程(2)。

    2 核運(yùn)動(dòng)Schr?dinger 方程的一種解法

    為了能夠解Schr?dinger 方程(2),需要求其折合質(zhì)量。

    圖1 四核坐標(biāo)示意圖

    1)折合質(zhì)量的求法

    為了求其折合質(zhì)量取坐標(biāo)原點(diǎn)為O 點(diǎn),M 個(gè)原子,核A,B,C,D….的坐標(biāo)向量分別是OA=R1,OB=R2,…,OC=RM( 見圖1),向量R*1的終點(diǎn)為質(zhì)心的坐標(biāo),它們質(zhì)量分別為m1,m2,…,mM,并且滿足

    其中m*=m1+m2+…+ mM。把上面的式子寫成矩陣形式

    再由式(11)和復(fù)合函數(shù)求全導(dǎo)數(shù)法,解得速度之間的關(guān)系為

    其中:v*和vi分別為質(zhì)心和第i 原子核的速度; vij為第i 與第j 個(gè)原子核的相對(duì)速度。式( 12)的兩端同乘一個(gè)三角矩陣可得

    和其轉(zhuǎn)置矩陣

    從式(12)~(14)和Laplace 算子▽2易知,動(dòng)能T 可以表達(dá)為

    它們的折合質(zhì)量或約化質(zhì)量分別為

    其中:R*2=R1M;R*j= |R*j|;θij是向量Ri+1i與向量Rj+1j之間的夾角;tj為向量模間的比例系數(shù),即

    為了解多原子分子核運(yùn)動(dòng)的Schr?dinger 方程(2),需要分解它。

    2)多原子分子核運(yùn)動(dòng)的Schr?dinger 方程的分解

    用分離變量法式(17)可化成下面等價(jià)的方程組[1]:

    其中:方程(19)是質(zhì)心平動(dòng)的方程,ε1表示其所的需能量;方程(20)是核相對(duì)運(yùn)動(dòng)的方程,ε-E-ε1表示核相對(duì)運(yùn)動(dòng)所的需能量,也就是核間距變化所的需能量。有趣的是,在分析式(10)時(shí)已知核間距變化所需要的能量為E*,所以二者應(yīng)該相等,即

    把式(10)代人式(21)消去E 得到

    把式(22)代入方程(20),方程(20)就轉(zhuǎn)化為

    這樣以來,只須解方程(19)和(23)就可以了。

    3)多原子分子核運(yùn)動(dòng)的Schr?dinger 方程的解析解

    ①對(duì)于質(zhì)心平動(dòng)方程(19),ε1是未知能量,應(yīng)用式(42),方程(23)的解是

    當(dāng)l1=0 時(shí)為常數(shù),質(zhì)心只能作平移運(yùn)動(dòng)[1],無旋轉(zhuǎn)運(yùn)動(dòng),于是質(zhì)心平動(dòng)方程(19)的嚴(yán)格解析解為

    其中C 是常數(shù)。再由式(43)和(44)得

    再把式(24)和式(26)代入式(18),得到方程(17)通解

    其中Cj為歸一化系數(shù)。至此,就找到了多原子分子Schr?dinger 方程組(1)和(2)的嚴(yán)格解析解,其完整表達(dá)式就是式(3)和式(29)。

    3 多原子分子解析解的重要性質(zhì)

    為了了解嚴(yán)格解析解(3)和(29)的重要性質(zhì),首先需要剖析它的各個(gè)因子。

    3)Jacobi 多項(xiàng)式F( -nj,nj+λj-1+3j/2-1,Lj+3/2|y)和Kummer 函數(shù)F(1,2,ρ1)及F( l2+1,2l+2,ρ2)顯然都是連續(xù)、單值、可導(dǎo)和平方可積的函數(shù),并且該Jacobi 多項(xiàng)式還能夠形成完全正交基函數(shù)的集[8,11]。于是,ω( nj,Li,λj-1|y)和YλNu( Ω)都是連續(xù)、單值、可導(dǎo)和平方可積的函數(shù)[11]。

    從1)~3)的論證可知,φ( r,Ω)和ψ( ρ1,ρ2,l2)都符合波函數(shù)的標(biāo)準(zhǔn)化條件[2-3],所以它們都是多原子分子的品優(yōu)波函數(shù)[2-3],并且利用它們還能夠求得其核間距。

    4 多原子分子核間距的計(jì)算方法

    由于核間距|Rij|的變化是由原子核運(yùn)動(dòng)決定的,所以確切的說,核間距是隨機(jī)變量,是不能夠十分準(zhǔn)確確定它的值的。但是,根據(jù)需要,人們還得去求它,測(cè)它。人們之所以能夠求到它,測(cè)到它,說明它出現(xiàn)的概率是很大的。因此,可以利用原子核相對(duì)運(yùn)動(dòng)的Schr?dinger 方程的解析解(29)來求。

    解析解ψ( ρ1,ρ2,l2)既然是個(gè)品優(yōu)波函數(shù),根據(jù)波函數(shù)統(tǒng)計(jì)規(guī)律,ψ( ρ1,ρ2,l2)應(yīng)該具有表示微粒運(yùn)動(dòng)幾率的功能。具體說來,以原子核A 為球心,以為ρ2半徑,另外一個(gè)原子核C 為動(dòng)點(diǎn),分布函數(shù)為[2]

    它表征發(fā)現(xiàn)該動(dòng)點(diǎn)處于單位厚度的球殼的幾率[2]。當(dāng)D 極大時(shí),dD/dρ2=0,即

    其中,Cj0是相應(yīng)的常數(shù)。這樣以來,就可以用式(30)求出極值點(diǎn)ρ2=ρ*,這說明ρ2=ρ*的幾率極大,也就是說,ρ2=ρ*的可能性極大。再由式(25)可知

    這就表明,R1M=ρ*/(2β2)的可能性也極大。正因?yàn)槿绱耍瑹o論是測(cè)量核間距|R1M|,或是計(jì)算核間距|R1M|,得到ρ*/(2β2)的可能性極大。這就是把ρ*/(2β2)當(dāng)作核間距|R1M|的值的理由。再用式(11)和(16)可以求出其它的核間距|Rij|的值。在這里自然需要假設(shè)各個(gè)核間距|Rij|的變化是同步的,即各個(gè)核間距|Rij|變化時(shí),它們的比值保持不變。從式(30)又可以看出,這些核間距|Rij|的值不僅與β2有關(guān),而且與Kummer 函數(shù)也有關(guān)。理所當(dāng)然,在某一個(gè)核間距Rij的值很大時(shí),就意味著分子的瓦解。那么瓦解時(shí)所需要的能量是多少呢?

    5 多原子分子的分解所需要的能量和總能量

    為了求總能量,先求出多原子分子核相對(duì)運(yùn)動(dòng)所需要的能量。

    1)多原子分子核相對(duì)運(yùn)動(dòng)所需要的能量

    由式(22)知

    再由式(31)知

    于是

    利用核相對(duì)運(yùn)動(dòng)所需要的能量就可以求得多原子分子的分解時(shí)所需的能量。

    2)多原子分子的分解所需要的能量

    在式(32)中,β2可由式(28)得到,即

    3)質(zhì)心平動(dòng)所的需要的能量

    從式(25)的第1 式能夠求得質(zhì)心平動(dòng)所的需要的能量

    特別地,當(dāng)質(zhì)心作勻速平動(dòng)或靜止時(shí),由方程(2)可以得到ε1=0。

    4)核運(yùn)動(dòng)所需要的能量

    由式(32)和(33)得

    式(34)表示質(zhì)心平動(dòng)所的需要的能量ε1和核相對(duì)運(yùn)動(dòng)所需要的能量E*之和,也就是核運(yùn)動(dòng)所需要的能量。下面就該求出多原子分子電子運(yùn)動(dòng)需要的能量。

    5)多原子分子電子運(yùn)動(dòng)需要的能量

    利用能量公式( 式(3)第2 式)就可以來計(jì)算電子運(yùn)動(dòng)需要的能量:

    也可以用式(10)來求電子運(yùn)動(dòng)需要的能量E,即

    再由式(32),上式變成

    用式(35)來計(jì)算電子運(yùn)動(dòng)所需的能量,可以省去積分運(yùn)算,比較方便,由此還可以去求總能量。

    6)多原子分子的總能量

    從式(21)可知總能量

    利用式(32)~(35)得

    式(36)就是總能量ε 的計(jì)算公式。利用式(35)和式(3)第2 式還可以得到

    由多原子分子總能量ε 的計(jì)算式(37),就可以導(dǎo)出其電離能的計(jì)算公式。

    7)多原子分子電離能的計(jì)算公式

    對(duì)于多原子分子,共有N 個(gè)電子。不妨選中一個(gè)原子的一個(gè)外層電子,可用式(36)或(37)來計(jì)算總能量ε。當(dāng)這個(gè)原子的一個(gè)外層電子丟失后,剩下N-1 個(gè)電子。也可以用式(36)或(37)重新計(jì)算其總能量:

    它的電離能應(yīng)該等于丟失電子前、后的總能量之差Δε,利用式(37)和(38)得

    當(dāng)多原子分子在丟失電子時(shí),在那一剎那間,分子來不及調(diào)整結(jié)構(gòu),也就是說,結(jié)構(gòu)是不變的,于是

    把式(40)代人(39)式中得

    式(41)就是多原子分子的一個(gè)外層電子丟失前、后的總能量差,即多原子分子電離能的計(jì)算公式。還可以看出,它也是被指定的原子的一個(gè)外層電子丟失前、后能量的差,即被指定原子的電離能的計(jì)算公式。也就是說,它既是分子電離能的計(jì)算公式,也是原子的電離能的計(jì)算公式。對(duì)于它的應(yīng)用,舉例說明如下。

    金屬鋁的晶胞是由4 個(gè)鋁原子構(gòu)成,Z=13,N=52,m=26.915 4。計(jì)算鋁的晶胞的電離能,其結(jié)果見表1。

    表1 鋁的晶胞電離能Δε(ev)和Lagueree 系數(shù)n 的關(guān)系

    從表1 的第1 列能夠看出,金屬鋁的晶胞電離能與其實(shí)驗(yàn)值8.986 ev 很接近[2],說明所得的解析解具有一定的可靠性。另外,從表1 還能夠看出,隨著Lagueree 函數(shù)的次數(shù)n 增大,電離能而減小,說明電離能與Lagueree 函數(shù)的次數(shù)n 等環(huán)境條件有密切關(guān)。

    血紅蛋白分子[5]Fe4C3032H4816O872N780S8有9 512 個(gè)原子核,原子序數(shù)分別是:28,6,1,8,7,16,共有38 076 個(gè)電子。現(xiàn)在不妨選中血紅蛋白分子的碳原子的一個(gè)外層電子,應(yīng)用式(41)計(jì)算,結(jié)果見表2。

    表2 血紅蛋白分子電離能Δε( -kj/mol)和Lagueree 系數(shù)

    從表2 能夠看出,血紅蛋白分子的電離能與Lagueree 函數(shù)的次數(shù)n 等環(huán)境條件有關(guān)。隨著Lagueree 函數(shù)的次數(shù)n 增大而減小。從表2 還可以看出,血紅蛋白分子的電離能很小很小,也正因?yàn)槿绱?,它很用容易失去電子,這就是很多病的病根和病原[2]。

    值得注意的是,在尋求解析解(3)和(29)的過程中,引用了式(42)~(44),它們來自下面的定理。

    定理已知Schr?dinger 方程

    其中E 表示未知能量,則其解析解為

    其中l(wèi) 是正整數(shù),F(xiàn)( l+1,2l+2,ρ)為Kummer 函數(shù),

    這個(gè)定理筆者已經(jīng)證明過,它是正確的[6]。這就說明筆者的探討是成功的。

    [1]唐敖慶.量子化學(xué)[M].北京:科學(xué)出版社,1982: 2-11.

    [2]徐光憲,王祥云.物質(zhì)結(jié)構(gòu)[M].北京:高等教育出版社,1959:3-55,589-591.

    [3]潘道曀,趙成大,鄭載興.物質(zhì)結(jié)構(gòu)[M].北京:等教育出版社,1982:5-101.

    [4]周公度,段連運(yùn).結(jié)構(gòu)化學(xué)基礎(chǔ)[M].北京:北京大學(xué)出版社,1997:392-396,429-434.

    [5]維達(dá)多夫.生物學(xué)與量子力學(xué)[M].李培廉,李景惠,張維華,譯.北京:科學(xué)出版社,1990,40-49,36-37.

    [6]任愛娣,王文信. 探討雙原子分子Schr?dinger 方程組的嚴(yán)格解析解[J].四川兵工學(xué)報(bào),2011,32(11):144-148.

    [7]任愛娣,王文信.探討多電子原子和分子Schr?dinger 方程的嚴(yán)格解析解[J].分子科學(xué)學(xué)報(bào),2011,27:291-296.

    [8]DENG C,ZHANG R Q,F(xiàn)ENG D C.Solution of Atom and MolecularSchr?dinger Equation Described Hyperspherical Coordinates[J].J Quant Chem Approach to the Description of Using Hyperspherical Coordinates,1993,45:385.

    [9]KINRK D L.Approach to the Description of Using Hyperspherical Coordinates[J].Chem Phys,1974,60:67.

    [10]劉式適,劉式達(dá).特殊函數(shù)[M].北京:Meteorologic Press,2002:597.

    [11]江澤堅(jiān),吳智泉.實(shí)變函數(shù)論[M].北京:高等教育出版社,1961:137-138,157-159.

    猜你喜歡
    電離能原子核原子
    學(xué)科融合視野下的探究教學(xué)實(shí)踐*
    ——以物質(zhì)結(jié)構(gòu)與性質(zhì)模塊“元素周期律”教學(xué)為例
    原子究竟有多???
    原子可以結(jié)合嗎?
    帶你認(rèn)識(shí)原子
    淺談A Level 化學(xué)中電離能的影響因素及變化趨勢(shì)與普高教學(xué)的差異
    關(guān)于原子核結(jié)構(gòu)的討論
    物質(zhì)構(gòu)成中的“一定”與“不一定”
    基于圖表分析的課堂知識(shí)建構(gòu)與能力提升
    卷煙煙氣化學(xué)成分的電離能計(jì)算及其在定性分析中的應(yīng)用
    煙草科技(2011年12期)2011-01-15 10:09:04
    走出半衰期的認(rèn)識(shí)誤區(qū)
    激情视频va一区二区三区| 精品国产乱子伦一区二区三区| 一级毛片女人18水好多| 不卡av一区二区三区| 亚洲成人免费电影在线观看| 精品亚洲成国产av| 十八禁网站免费在线| 在线观看舔阴道视频| 女性被躁到高潮视频| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦 在线观看视频| 午夜福利在线观看吧| 国产成人免费观看mmmm| 在线观看舔阴道视频| 纯流量卡能插随身wifi吗| 欧美老熟妇乱子伦牲交| 青草久久国产| 黄片播放在线免费| 国产一区二区三区视频了| 男女午夜视频在线观看| 午夜福利免费观看在线| 一进一出抽搐动态| 亚洲中文日韩欧美视频| 一级a爱视频在线免费观看| 99久久人妻综合| 亚洲精品美女久久久久99蜜臀| 一级,二级,三级黄色视频| 久久久久视频综合| 五月天丁香电影| 99久久人妻综合| 成年人免费黄色播放视频| 午夜91福利影院| 国产高清激情床上av| 国产精品自产拍在线观看55亚洲 | 青草久久国产| 看免费av毛片| 伊人久久大香线蕉亚洲五| 欧美变态另类bdsm刘玥| 欧美亚洲 丝袜 人妻 在线| 法律面前人人平等表现在哪些方面| 精品高清国产在线一区| 女人被躁到高潮嗷嗷叫费观| 美女高潮到喷水免费观看| 如日韩欧美国产精品一区二区三区| www.熟女人妻精品国产| 一区二区三区激情视频| 91精品三级在线观看| 成人手机av| 国产人伦9x9x在线观看| 999久久久精品免费观看国产| 国产不卡av网站在线观看| 在线看a的网站| 日韩中文字幕视频在线看片| 9热在线视频观看99| 亚洲第一av免费看| av网站免费在线观看视频| 欧美成人免费av一区二区三区 | 久热爱精品视频在线9| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人| 国产精品欧美亚洲77777| 国产精品久久久久久精品电影小说| 在线观看一区二区三区激情| cao死你这个sao货| 久久ye,这里只有精品| 国产精品一区二区精品视频观看| 国产精品免费一区二区三区在线 | 国产国语露脸激情在线看| 亚洲 欧美一区二区三区| 蜜桃国产av成人99| 国产成人精品在线电影| 欧美一级毛片孕妇| 又紧又爽又黄一区二区| 日本黄色日本黄色录像| 精品久久久久久电影网| 欧美精品啪啪一区二区三区| 18禁国产床啪视频网站| 男人操女人黄网站| 国产免费视频播放在线视频| 男女午夜视频在线观看| 啪啪无遮挡十八禁网站| 国产男女内射视频| 亚洲精品在线美女| 欧美日韩中文字幕国产精品一区二区三区 | 久久久国产精品麻豆| 99热国产这里只有精品6| 久久 成人 亚洲| 亚洲专区中文字幕在线| 免费观看人在逋| 在线av久久热| 午夜激情av网站| 80岁老熟妇乱子伦牲交| 后天国语完整版免费观看| 青青草视频在线视频观看| 婷婷成人精品国产| 黄色视频,在线免费观看| 精品国产一区二区久久| 99国产精品一区二区三区| 国产日韩欧美在线精品| 久久久久国内视频| 精品国产乱子伦一区二区三区| 12—13女人毛片做爰片一| 亚洲精品粉嫩美女一区| 老司机午夜福利在线观看视频 | 亚洲性夜色夜夜综合| 亚洲人成电影免费在线| 欧美精品啪啪一区二区三区| 一本一本久久a久久精品综合妖精| 亚洲色图 男人天堂 中文字幕| 99精品在免费线老司机午夜| 免费少妇av软件| 亚洲一码二码三码区别大吗| 天堂8中文在线网| 欧美 日韩 精品 国产| 国产免费现黄频在线看| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女 | 夜夜骑夜夜射夜夜干| 日本欧美视频一区| 一区二区三区精品91| 亚洲 欧美一区二区三区| 日韩制服丝袜自拍偷拍| 精品福利永久在线观看| av又黄又爽大尺度在线免费看| 我要看黄色一级片免费的| 男女无遮挡免费网站观看| av超薄肉色丝袜交足视频| 国产精品一区二区在线不卡| 亚洲精品国产区一区二| 久久人妻福利社区极品人妻图片| av一本久久久久| 999久久久精品免费观看国产| av有码第一页| 国产男女内射视频| 日韩中文字幕欧美一区二区| 大陆偷拍与自拍| 少妇的丰满在线观看| 无限看片的www在线观看| 日本欧美视频一区| 视频在线观看一区二区三区| 亚洲中文日韩欧美视频| 精品人妻在线不人妻| 国产精品成人在线| 天天躁日日躁夜夜躁夜夜| 免费少妇av软件| 18在线观看网站| 精品少妇久久久久久888优播| 男人舔女人的私密视频| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图| 国产成人欧美| 亚洲av日韩在线播放| av天堂久久9| 12—13女人毛片做爰片一| 成在线人永久免费视频| 久久久精品94久久精品| 国产深夜福利视频在线观看| av网站免费在线观看视频| 人人澡人人妻人| 亚洲精品成人av观看孕妇| 极品少妇高潮喷水抽搐| 亚洲人成电影观看| av天堂久久9| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利,免费看| 黄色片一级片一级黄色片| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久人人做人人爽| 日本av免费视频播放| 亚洲午夜理论影院| a级片在线免费高清观看视频| 在线观看免费午夜福利视频| 欧美精品亚洲一区二区| 国产精品一区二区在线观看99| 欧美日韩亚洲国产一区二区在线观看 | 久久精品国产综合久久久| 亚洲精品久久成人aⅴ小说| 亚洲色图 男人天堂 中文字幕| 美国免费a级毛片| av视频免费观看在线观看| 精品熟女少妇八av免费久了| 在线 av 中文字幕| 视频区欧美日本亚洲| 久久久国产一区二区| 欧美日韩黄片免| 久久中文字幕人妻熟女| 精品国产乱子伦一区二区三区| 午夜福利视频在线观看免费| 国产熟女午夜一区二区三区| 国产精品久久久久久精品电影小说| 捣出白浆h1v1| 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 中亚洲国语对白在线视频| 亚洲精品av麻豆狂野| 日韩视频在线欧美| 欧美日韩福利视频一区二区| 黑人巨大精品欧美一区二区mp4| 国产亚洲欧美在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 韩国精品一区二区三区| 欧美日韩成人在线一区二区| 女性被躁到高潮视频| 国产在线视频一区二区| 色综合欧美亚洲国产小说| 丝袜在线中文字幕| cao死你这个sao货| 天天躁日日躁夜夜躁夜夜| 1024香蕉在线观看| 欧美日韩亚洲国产一区二区在线观看 | 黄色视频在线播放观看不卡| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看 | 老司机在亚洲福利影院| 叶爱在线成人免费视频播放| 精品少妇久久久久久888优播| 日本欧美视频一区| 性少妇av在线| 精品国产亚洲在线| 国产精品 国内视频| 免费一级毛片在线播放高清视频 | 99国产精品99久久久久| 在线永久观看黄色视频| 国产成人免费观看mmmm| 亚洲专区中文字幕在线| 免费看十八禁软件| 黑人猛操日本美女一级片| 淫妇啪啪啪对白视频| 女人被躁到高潮嗷嗷叫费观| 99精品欧美一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 老熟妇仑乱视频hdxx| 久久亚洲真实| 亚洲成人手机| 天堂动漫精品| 高清av免费在线| 久久婷婷成人综合色麻豆| 狠狠狠狠99中文字幕| 精品少妇久久久久久888优播| 久久狼人影院| 18禁裸乳无遮挡动漫免费视频| 亚洲精品在线美女| 久久久久国内视频| 97人妻天天添夜夜摸| 久久久国产精品麻豆| 久久ye,这里只有精品| 波多野结衣av一区二区av| 男人操女人黄网站| 天堂俺去俺来也www色官网| 91av网站免费观看| 99九九在线精品视频| 免费在线观看完整版高清| 日韩中文字幕视频在线看片| videos熟女内射| 亚洲色图 男人天堂 中文字幕| 亚洲午夜精品一区,二区,三区| 国产免费现黄频在线看| 国产av国产精品国产| 午夜福利一区二区在线看| e午夜精品久久久久久久| 亚洲精品国产区一区二| 亚洲人成电影观看| 国产精品98久久久久久宅男小说| 亚洲av第一区精品v没综合| 国产1区2区3区精品| 亚洲精品一二三| 亚洲第一欧美日韩一区二区三区 | 国产精品偷伦视频观看了| 在线十欧美十亚洲十日本专区| 日韩三级视频一区二区三区| 久久久精品94久久精品| av欧美777| 操出白浆在线播放| 亚洲中文字幕日韩| 国产精品久久久久久人妻精品电影 | 国产在线精品亚洲第一网站| 欧美在线黄色| 成年女人毛片免费观看观看9 | tube8黄色片| 最黄视频免费看| 亚洲成人免费电影在线观看| 国产黄色免费在线视频| 99国产精品一区二区三区| 国产免费现黄频在线看| 日日摸夜夜添夜夜添小说| 宅男免费午夜| 2018国产大陆天天弄谢| 中文字幕制服av| 一边摸一边抽搐一进一小说 | 精品国产一区二区三区四区第35| 视频区图区小说| 热re99久久精品国产66热6| 男女下面插进去视频免费观看| 黑人操中国人逼视频| 亚洲精品粉嫩美女一区| 99精品在免费线老司机午夜| 亚洲精品美女久久av网站| 欧美精品av麻豆av| 在线永久观看黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 91字幕亚洲| 香蕉国产在线看| 久久国产精品影院| 国产精品久久久久久人妻精品电影 | 亚洲精品美女久久av网站| 欧美大码av| 在线观看免费视频网站a站| 国产成人欧美在线观看 | 黑丝袜美女国产一区| 色视频在线一区二区三区| 另类亚洲欧美激情| svipshipincom国产片| 国产成人系列免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲国产成人一精品久久久| 久久香蕉激情| 激情视频va一区二区三区| 亚洲一区二区三区欧美精品| 亚洲五月色婷婷综合| 亚洲七黄色美女视频| 美国免费a级毛片| 露出奶头的视频| 日日夜夜操网爽| 人人妻人人爽人人添夜夜欢视频| www.熟女人妻精品国产| 国产精品国产av在线观看| 婷婷丁香在线五月| 日韩成人在线观看一区二区三区| 国产麻豆69| 亚洲国产欧美日韩在线播放| 美女主播在线视频| 日韩免费av在线播放| 波多野结衣一区麻豆| 妹子高潮喷水视频| 香蕉国产在线看| 黄网站色视频无遮挡免费观看| av超薄肉色丝袜交足视频| a级毛片黄视频| 肉色欧美久久久久久久蜜桃| 成人永久免费在线观看视频 | 黄色成人免费大全| 一级,二级,三级黄色视频| 亚洲午夜精品一区,二区,三区| 亚洲欧美一区二区三区黑人| a级毛片黄视频| 在线看a的网站| 久久精品亚洲熟妇少妇任你| 欧美在线一区亚洲| 亚洲欧美精品综合一区二区三区| 一边摸一边抽搐一进一小说 | 久久av网站| 中文字幕最新亚洲高清| 免费在线观看完整版高清| 最新美女视频免费是黄的| 宅男免费午夜| 黄色怎么调成土黄色| 天堂动漫精品| 午夜免费鲁丝| 国产一区二区激情短视频| 亚洲色图综合在线观看| 女性生殖器流出的白浆| 亚洲欧美一区二区三区久久| 免费女性裸体啪啪无遮挡网站| 成人永久免费在线观看视频 | 日韩中文字幕欧美一区二区| 如日韩欧美国产精品一区二区三区| 另类精品久久| 久久香蕉激情| av线在线观看网站| 国产成人欧美在线观看 | 精品福利永久在线观看| 麻豆成人av在线观看| a在线观看视频网站| 亚洲中文字幕日韩| 久9热在线精品视频| 国产一区二区三区视频了| 亚洲天堂av无毛| 美女视频免费永久观看网站| 免费在线观看黄色视频的| 狠狠精品人妻久久久久久综合| 99香蕉大伊视频| 国产高清视频在线播放一区| 18禁黄网站禁片午夜丰满| 欧美+亚洲+日韩+国产| 亚洲人成伊人成综合网2020| 午夜福利视频在线观看免费| 大片电影免费在线观看免费| 久久天堂一区二区三区四区| 黄色怎么调成土黄色| 在线天堂中文资源库| 一级毛片女人18水好多| 国产精品国产av在线观看| 精品国产一区二区久久| 亚洲专区国产一区二区| 成人免费观看视频高清| 在线av久久热| 九色亚洲精品在线播放| 国产男女超爽视频在线观看| 人人澡人人妻人| 亚洲精品久久成人aⅴ小说| 欧美人与性动交α欧美精品济南到| 欧美精品啪啪一区二区三区| 国产欧美亚洲国产| e午夜精品久久久久久久| 午夜视频精品福利| a在线观看视频网站| 天天躁日日躁夜夜躁夜夜| 日本黄色日本黄色录像| 国产精品欧美亚洲77777| 午夜精品国产一区二区电影| 久久久久久久久久久久大奶| 精品国产国语对白av| 国产一卡二卡三卡精品| 操美女的视频在线观看| 麻豆乱淫一区二区| 国产男女超爽视频在线观看| 色视频在线一区二区三区| av在线播放免费不卡| 亚洲人成电影免费在线| 欧美精品亚洲一区二区| 国产午夜精品久久久久久| 欧美国产精品一级二级三级| 新久久久久国产一级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看66精品国产| 国产伦人伦偷精品视频| 亚洲 欧美一区二区三区| 日韩欧美一区视频在线观看| 男人操女人黄网站| 一区二区三区乱码不卡18| 亚洲国产看品久久| 女警被强在线播放| 国产亚洲av高清不卡| 精品久久久久久久毛片微露脸| 人人妻人人添人人爽欧美一区卜| 午夜福利免费观看在线| 五月天丁香电影| 亚洲av国产av综合av卡| 在线亚洲精品国产二区图片欧美| 久久精品国产99精品国产亚洲性色 | 女人爽到高潮嗷嗷叫在线视频| 中文字幕人妻丝袜一区二区| 正在播放国产对白刺激| 日韩精品免费视频一区二区三区| 69精品国产乱码久久久| 亚洲色图av天堂| 国产有黄有色有爽视频| 在线观看免费高清a一片| 人人妻人人爽人人添夜夜欢视频| 最新美女视频免费是黄的| 亚洲视频免费观看视频| 久久久久视频综合| 日韩成人在线观看一区二区三区| 国产精品亚洲av一区麻豆| 纯流量卡能插随身wifi吗| 青草久久国产| 亚洲熟女毛片儿| 亚洲黑人精品在线| 亚洲av成人一区二区三| av视频免费观看在线观看| 亚洲avbb在线观看| 国产野战对白在线观看| 成人手机av| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一小说 | 久久国产亚洲av麻豆专区| 亚洲成人国产一区在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产熟女午夜一区二区三区| 久久精品亚洲精品国产色婷小说| 久久精品人人爽人人爽视色| 免费在线观看日本一区| 亚洲男人天堂网一区| 久久久精品94久久精品| 国产成人系列免费观看| av片东京热男人的天堂| 少妇粗大呻吟视频| 制服诱惑二区| 亚洲五月色婷婷综合| 久久人妻熟女aⅴ| 久久国产亚洲av麻豆专区| 国产免费av片在线观看野外av| 热99re8久久精品国产| 999久久久精品免费观看国产| 亚洲av日韩在线播放| 午夜日韩欧美国产| 一本一本久久a久久精品综合妖精| 人妻久久中文字幕网| 91成人精品电影| 无人区码免费观看不卡 | 成人免费观看视频高清| 亚洲精品国产一区二区精华液| 在线亚洲精品国产二区图片欧美| 最新的欧美精品一区二区| 国产精品一区二区精品视频观看| 天堂中文最新版在线下载| 日日爽夜夜爽网站| 亚洲精品国产精品久久久不卡| 最新美女视频免费是黄的| 两个人免费观看高清视频| 一本—道久久a久久精品蜜桃钙片| 国产高清视频在线播放一区| 免费在线观看黄色视频的| 欧美+亚洲+日韩+国产| 精品国产一区二区久久| 男人舔女人的私密视频| 99香蕉大伊视频| 飞空精品影院首页| 欧美精品一区二区免费开放| 9热在线视频观看99| 人成视频在线观看免费观看| 欧美性长视频在线观看| 香蕉丝袜av| 热99国产精品久久久久久7| 男女之事视频高清在线观看| 免费在线观看视频国产中文字幕亚洲| 黄色视频,在线免费观看| 一区二区三区国产精品乱码| 性高湖久久久久久久久免费观看| 两个人看的免费小视频| 国产精品免费一区二区三区在线 | 亚洲黑人精品在线| 99国产精品一区二区三区| 黄频高清免费视频| 91精品国产国语对白视频| 日韩熟女老妇一区二区性免费视频| 国产在线精品亚洲第一网站| 黄片小视频在线播放| 亚洲欧洲精品一区二区精品久久久| av欧美777| 在线播放国产精品三级| h视频一区二区三区| 亚洲国产欧美在线一区| 精品一区二区三卡| 色综合欧美亚洲国产小说| 国产精品香港三级国产av潘金莲| 欧美一级毛片孕妇| 99re在线观看精品视频| 熟女少妇亚洲综合色aaa.| 丁香六月天网| 亚洲国产av影院在线观看| 飞空精品影院首页| bbb黄色大片| 好男人电影高清在线观看| 后天国语完整版免费观看| 国产成人欧美| 一进一出抽搐动态| 国产在线一区二区三区精| 色婷婷av一区二区三区视频| www.自偷自拍.com| 精品久久久久久电影网| 97在线人人人人妻| 午夜老司机福利片| 欧美精品高潮呻吟av久久| 婷婷丁香在线五月| 人妻 亚洲 视频| 一级毛片电影观看| 一进一出好大好爽视频| 男女高潮啪啪啪动态图| 午夜福利一区二区在线看| 国产在线精品亚洲第一网站| 99re6热这里在线精品视频| 脱女人内裤的视频| 视频区欧美日本亚洲| 日本av手机在线免费观看| 国产精品免费一区二区三区在线 | 香蕉国产在线看| 成人特级黄色片久久久久久久 | 精品一区二区三区av网在线观看 | 日本av免费视频播放| 在线观看www视频免费| 欧美日韩视频精品一区| 成人av一区二区三区在线看| 人人妻人人澡人人看| 女人高潮潮喷娇喘18禁视频| 亚洲成国产人片在线观看| 国产又爽黄色视频| 高清在线国产一区| 深夜精品福利| 日日摸夜夜添夜夜添小说| 大陆偷拍与自拍| 欧美国产精品va在线观看不卡| 男女免费视频国产| 国产精品98久久久久久宅男小说| 亚洲人成电影观看| 欧美 日韩 精品 国产| 久久这里只有精品19| 一个人免费看片子| av网站在线播放免费| 午夜福利一区二区在线看| 天堂8中文在线网| 亚洲中文av在线| 一区二区三区乱码不卡18| 亚洲精品中文字幕一二三四区 | 狠狠婷婷综合久久久久久88av| 国产精品久久久久久精品电影小说| 亚洲精品自拍成人| 精品国产超薄肉色丝袜足j| 久久午夜亚洲精品久久| 亚洲精品一卡2卡三卡4卡5卡| 国产黄频视频在线观看| 狠狠精品人妻久久久久久综合| 国产成人一区二区三区免费视频网站| 人人妻,人人澡人人爽秒播| 大码成人一级视频| 国产精品香港三级国产av潘金莲| 中文字幕最新亚洲高清| 老汉色av国产亚洲站长工具| 成人精品一区二区免费| 亚洲精品中文字幕一二三四区 | 最近最新中文字幕大全免费视频| 侵犯人妻中文字幕一二三四区| 国产欧美日韩一区二区精品| 老熟女久久久| 国产区一区二久久|