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

    分子動力學(xué)方法中速度標(biāo)度因子的推導(dǎo)

    2022-04-08 06:16:58王圣軍
    大學(xué)物理 2022年4期
    關(guān)鍵詞:方法系統(tǒng)

    齊 楠,王圣軍

    (陜西師范大學(xué) 物理學(xué)與信息技術(shù)學(xué)院,陜西 西安 710119 )

    計算機(jī)模擬現(xiàn)已成為實驗方法和理論研究以外的一種有效的物理研究手段,而分子動力學(xué)方法作為一種確定性模擬方法,其在計算機(jī)模擬中顯得尤為重要.在許多計算物理教材[1-3]中,對于模擬由大量粒子組成的系統(tǒng),分子動力學(xué)方法都是尤為重要的內(nèi)容.

    分子動力學(xué)方法是指利用每個分子的運動方程通過統(tǒng)計方法來計算系統(tǒng)的性質(zhì).其以粒子間的相互作用為出發(fā)點,在計算機(jī)上求解運動方程的數(shù)值解.在這一求解過程中,涉及的物理量包括能量E、力F、質(zhì)量m、溫度T、位置r、時間t、速度v等.然而,在系統(tǒng)中,這些物理量的大小在國際單位制下具有過于小的數(shù)量級,這不利于在計算機(jī)上進(jìn)行運算.于是,在分子動力學(xué)方法中,首先需要對物理量重新標(biāo)度.這樣,對于重新標(biāo)度后的物理量的大小,其數(shù)量級在1附近,并且能夠得到形式簡單的運動方程.

    在分子動力學(xué)模擬中,重新標(biāo)度的方法是理解分子動力學(xué)方法的重要環(huán)節(jié).在重新標(biāo)度的過程中,需要考慮系統(tǒng)的運動特征,這具有一定的物理內(nèi)涵.然而其在多數(shù)教材中顯得過于簡略,例如沒有介紹時間單位的物理內(nèi)涵.另外在模擬過程中標(biāo)度因子是一個尤為關(guān)鍵的量.在微正則系綜中,給定初始位置和初始速度以后,在系統(tǒng)趨衡過程中需要調(diào)節(jié)速度,使系統(tǒng)的溫度達(dá)到給定系統(tǒng)的溫度.對速度的重新標(biāo)度通過標(biāo)度因子實現(xiàn).然而,與之相關(guān)的書籍[1-3]往往直接給出了標(biāo)度因子的數(shù)學(xué)表達(dá)式,沒有明確地指出如何合理地對其進(jìn)行推導(dǎo).這使得分子動力學(xué)模擬的原理難以被學(xué)生充分理解.

    接下來,本文將對標(biāo)度因子的表達(dá)式予以推導(dǎo),并且對物理量重新標(biāo)度的方法給出合理的說明.這樣會得到重新標(biāo)度的運動方程、物理量和標(biāo)度因子,幫助初學(xué)者掌握完整的模擬原理,以便于在面對新的情景時推導(dǎo)相應(yīng)的公式開展計算機(jī)模擬.

    1 重新標(biāo)度前的公式

    在分子動力學(xué)模擬過程中,首先要設(shè)定模擬所采用的模型.這里考察最基本的情況微正則系綜,并且采用單原子粒子系統(tǒng),分子間的相互作用勢采用Lennard-Jones勢.系統(tǒng)的哈密頓量表示為

    (1)

    其中分子間相互作用勢u(rij)為

    (2)

    這里rij是粒子i與粒子j的距離,ε和σ是兩個參數(shù).

    由于統(tǒng)計物理系統(tǒng)的特性,精確地控制初態(tài)是困難的.模擬過程中人為設(shè)置的系統(tǒng)初態(tài)不會具有所要求的能量,并且系統(tǒng)的狀態(tài)不是一個平衡態(tài).為了使該系統(tǒng)的能量達(dá)到一個確定值,在系統(tǒng)趨衡的過程中調(diào)節(jié)系統(tǒng)能量,使其逐步達(dá)到所要求的狀態(tài).熱力學(xué)平衡狀態(tài)下該系統(tǒng)具有給定的溫度T,根據(jù)能量均分定理有

    (3)

    這里kB是玻耳茲曼常量.

    (4)

    (5)

    從而得到標(biāo)度因子為

    (6)

    為了進(jìn)行計算機(jī)模擬,通常對系統(tǒng)進(jìn)行無量綱化處理.教科書中的速度標(biāo)度因子都是在無量綱化的模型中給出的.通常它的形式與式(6)的形式并不相同,并且它依賴于無量綱化的處理方式.

    2 物理量的重新標(biāo)度

    粒子i的牛頓運動方程為

    (7)

    (8)

    于是

    (9)

    其中rij=ri-rj.對于計算機(jī)模擬,需要用到式(6)、(7)、(9).

    這里需要對物理量進(jìn)行重新標(biāo)度,進(jìn)而把分子動力學(xué)模擬中用到的物理量變成無量綱的量,使運動方程具有簡潔的形式.

    另外,如果在模擬中仍然采用國際單位制,那么粒子的質(zhì)量m、位置r和基本時間步長h的數(shù)量級均很小,可能會導(dǎo)致舍入誤差的產(chǎn)生.于是,對物理量進(jìn)行重新標(biāo)度的另一個原則是:避免過于小或過于大的數(shù)量級出現(xiàn)在計算機(jī)模擬中[4].

    為了得到具有簡潔形式的運動方程,本文并不是直接對力學(xué)中常用的基本物理量進(jìn)行重新標(biāo)度,再導(dǎo)出其它物理量.在分子動力學(xué)模擬中,選取一些反映相互作用和運動特征的物理量,首先對它們進(jìn)行恰當(dāng)?shù)闹匦聵?biāo)度.進(jìn)而把這些重新標(biāo)度的物理量作為新的基本物理量[5,6],由新的基本物理量導(dǎo)出其它物理量、標(biāo)度因子和運動方程.

    現(xiàn)在對物理量重新標(biāo)度,并以星號*表示重新標(biāo)度后物理量的數(shù)值.從后文可以看到,一方面物理量的數(shù)值的數(shù)量級不會遠(yuǎn)遠(yuǎn)偏離1,另一方面運動方程經(jīng)過了無量綱處理.

    根據(jù)式(2),容易得到-ε是分子間相互作用能的最小值,這個最小值在r=21/6σ處.參數(shù)ε決定了兩個粒子處于平衡位置的相互作用能,參數(shù)σ決定了兩個粒子處于平衡位置的距離.規(guī)定σ為長度單位,ε為能量單位,得到變換關(guān)系:

    l=l*σ

    (10)

    E=E*ε

    (11)

    考慮粒子i與粒子j間相互作用的時間.在其質(zhì)心系中,將兩粒子所構(gòu)成的系統(tǒng)的總能量E表示為E=aε,則a為一個無量綱的參量.由該系統(tǒng)能量守恒[7]得到

    (12)

    (13)

    上式為系統(tǒng)運動時長t關(guān)于兩粒子間距離rij的微分表達(dá)式.容易得到兩粒子間由最小距離r1處首次運動到最大距離r2處經(jīng)歷的時間T[7]為

    (14)

    圖1 的數(shù)值模擬結(jié)果

    粒子間相互作用的時間也可以從以下的角度考慮.當(dāng)一個粒子的振動范圍很小時,可以對系統(tǒng)中其它粒子對該粒子的作用做線性近似.考慮粒子i在單個粒子的勢場中的運動.粒子i的位置用ri表示,考察粒子i的運動周期.由式(8),粒子i在勢場中受到的作用力f(ri)為

    (15)

    則r0=21/6σ處是粒子i振動的平衡位置.將f(ri)在r0處做泰勒展開,保留線性項,得到

    (16)

    其中x=ri-r0.于是粒子i的牛頓運動方程為

    (17)

    (18)

    這里時間單位中包含一個常數(shù)48.它的作用是把相互作用力代入牛頓第二定律的時候,得到的運動方程形式更加簡單.

    (19)

    于是m*=48,即重新標(biāo)度后質(zhì)量m具有m*=48的數(shù)值大小.

    由式(7)、(9)、(10)、(18),得到重新標(biāo)度后的牛頓運動方程,其表示為

    (20)

    (21)

    這里h*是重新標(biāo)度后的基本時間步長.

    (22)

    (23)

    由式(11)、(23),重新標(biāo)度后平均動能表示為

    (24)

    由式(6)、(22)、(24),重新標(biāo)度后標(biāo)度因子β表示為

    (25)

    這里得到的標(biāo)度因子與相關(guān)書籍中的結(jié)果是一致的[1-3].當(dāng)導(dǎo)出了所有相關(guān)的公式以后,為了簡便,忽略所有無量綱物理量的*號.

    3 雙原子分子與量子分子動力學(xué)方法

    原則上,由大量粒子組成的系統(tǒng)都可應(yīng)用分子動力學(xué)方法模擬.這個系統(tǒng)包括的不只是點粒子系統(tǒng),也包括具有內(nèi)部結(jié)構(gòu)的粒子組成的系統(tǒng).如果按照系統(tǒng)所遵循的運動規(guī)律不同,分子動力學(xué)方法可分為經(jīng)典分子動力學(xué)方法和量子分子動力學(xué)方法.經(jīng)典分子動力學(xué)方法中,系統(tǒng)中每個分子各自都服從經(jīng)典運動定律.有時必須考慮系統(tǒng)中的量子效應(yīng),這需要用量子力學(xué)來處理.這種考慮微觀組態(tài)變化影響的分子動力學(xué)模擬稱為量子分子動力學(xué)模擬,也稱為從頭計算的分子動力學(xué)模擬或第一性原理計算[8].量子分子動力學(xué)方法處理的是多電子的薛定諤方程:

    Hψ(r,R)=Eψ(r,R)

    (26)

    其中,r和R分別代表所有電子和原子核的坐標(biāo)[8].基于密度泛函理論的第一性原理的計算方法由Car與Parrinello于1985年提出并用于實際計算.在過去的幾十年中,分子動力學(xué)模擬發(fā)展很快,并在許多領(lǐng)域得到了廣泛的應(yīng)用[8,9].

    根據(jù)分子動力學(xué)方法的基本步驟,系統(tǒng)中相互作用勢函數(shù)的確定尤為關(guān)鍵.在分子動力學(xué)方法中,通常是基于經(jīng)驗數(shù)據(jù)和經(jīng)驗規(guī)律來選取勢函數(shù).除了常用的Lennard-Jones勢,還有Morse勢、Kihara勢、諧振子勢、偶極子勢等.其中,Morse勢極好地描述了雙原子分子的振動,其數(shù)學(xué)形式如下[10]:

    V(r)=D(e-2αx-2e-αx)

    (27)

    在薛定諤方程的數(shù)值求解中,通常以埃為單位來度量所有的長度,以電子伏為單位來度量一切能量[11].在這樣一種單位制下,方程中的物理量和常數(shù)具有較好的數(shù)量級.在所謂的量子分子動力學(xué)方法中,也可以類似地對系統(tǒng)做無量綱化處理.只需要對長度和能量分別進(jìn)行重新標(biāo)度.這種無量綱化處理方式與前文所述是一致的.

    因此,確定好勢函數(shù)和運動方程后,總是可以按照本文的無量綱化處理方式對物理量進(jìn)行重新標(biāo)度,使運動方程更加簡潔.值得注意的是,其中對基本物理量的選取和重新標(biāo)度依賴于具體的物理模型,這具有一定的物理內(nèi)涵.

    4 結(jié)語

    本文首先給出了趨衡過程中速度標(biāo)度因子在國際單位制下的數(shù)學(xué)表達(dá)式.然后,說明了如何把物理量重新標(biāo)度得到無量綱公式,使得在模擬中計算更加方便.本文通過物理量的重新標(biāo)度,給出了時間單位的物理意義,并且給出了教材中常用的標(biāo)度因子相應(yīng)數(shù)學(xué)表達(dá)式的推導(dǎo)過程.文中的方法可以幫助初學(xué)者掌握完整的模擬原理,以便于在面對新的情景時推導(dǎo)相應(yīng)的公式開展計算機(jī)模擬.

    猜你喜歡
    方法系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    WJ-700無人機(jī)系統(tǒng)
    ZC系列無人機(jī)遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    基于PowerPC+FPGA顯示系統(tǒng)
    學(xué)習(xí)方法
    半沸制皂系統(tǒng)(下)
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    欧美成人a在线观看| 日本与韩国留学比较| 可以在线观看的亚洲视频| 91久久精品电影网| av专区在线播放| 国产中年淑女户外野战色| 精品免费久久久久久久清纯| 天堂√8在线中文| 国产精品无大码| 亚洲一区二区三区色噜噜| 日日摸夜夜添夜夜添小说| 精品日产1卡2卡| 亚洲自拍偷在线| 久久精品国产亚洲av涩爱 | 欧美成人一区二区免费高清观看| 亚州av有码| 久久亚洲精品不卡| 蜜臀久久99精品久久宅男| av.在线天堂| av天堂中文字幕网| 久久久久免费精品人妻一区二区| 精品午夜福利视频在线观看一区| 一边摸一边抽搐一进一小说| 舔av片在线| 久久精品久久久久久噜噜老黄 | 国产精品女同一区二区软件| 国产不卡一卡二| 国内精品美女久久久久久| 欧美日韩国产亚洲二区| 又黄又爽又刺激的免费视频.| 国产人妻一区二区三区在| 欧美性猛交黑人性爽| 精华霜和精华液先用哪个| 亚洲经典国产精华液单| 国产探花极品一区二区| 在线国产一区二区在线| 99久久成人亚洲精品观看| 国产精品一二三区在线看| 久久九九热精品免费| 日本免费a在线| 天堂√8在线中文| 丰满的人妻完整版| 午夜福利18| 人妻夜夜爽99麻豆av| 国产精品乱码一区二三区的特点| 看片在线看免费视频| 淫妇啪啪啪对白视频| 99久国产av精品| 国产中年淑女户外野战色| 搡老妇女老女人老熟妇| 97在线视频观看| 久久99热这里只有精品18| 婷婷色综合大香蕉| 波多野结衣高清无吗| 精品熟女少妇av免费看| 国内精品一区二区在线观看| 尾随美女入室| 国产乱人偷精品视频| 欧美极品一区二区三区四区| 国产精品一二三区在线看| 亚洲成人久久爱视频| 毛片女人毛片| 国产精品久久视频播放| 日韩av在线大香蕉| 超碰av人人做人人爽久久| 九色成人免费人妻av| 日日撸夜夜添| 国产一区亚洲一区在线观看| 97碰自拍视频| 免费观看在线日韩| 亚洲欧美日韩东京热| 成熟少妇高潮喷水视频| 国产人妻一区二区三区在| 亚洲婷婷狠狠爱综合网| 日本免费a在线| 国产精品精品国产色婷婷| 长腿黑丝高跟| 国产一区亚洲一区在线观看| 亚州av有码| 美女免费视频网站| 一级黄色大片毛片| 色5月婷婷丁香| 天堂影院成人在线观看| 老司机午夜福利在线观看视频| 男人和女人高潮做爰伦理| 日本免费a在线| 99在线人妻在线中文字幕| 国产私拍福利视频在线观看| 18禁黄网站禁片免费观看直播| av在线观看视频网站免费| 日本五十路高清| 看黄色毛片网站| 男女那种视频在线观看| 国产v大片淫在线免费观看| 女人被狂操c到高潮| 在线免费观看不下载黄p国产| 人妻少妇偷人精品九色| 性色avwww在线观看| 18禁在线无遮挡免费观看视频 | 亚洲欧美日韩高清专用| 久久久久久九九精品二区国产| 免费av不卡在线播放| 直男gayav资源| 最近2019中文字幕mv第一页| 日韩亚洲欧美综合| 好男人在线观看高清免费视频| 99久国产av精品国产电影| 亚洲精品国产av成人精品 | 免费观看精品视频网站| 小蜜桃在线观看免费完整版高清| 老女人水多毛片| 国产大屁股一区二区在线视频| av在线老鸭窝| 又爽又黄a免费视频| 国产综合懂色| 久久久久久久亚洲中文字幕| 成人综合一区亚洲| 联通29元200g的流量卡| 精品国内亚洲2022精品成人| 国产精品野战在线观看| 久久精品国产鲁丝片午夜精品| 最近的中文字幕免费完整| 午夜激情福利司机影院| 亚洲欧美中文字幕日韩二区| 在线免费十八禁| 亚洲欧美清纯卡通| 欧美日韩综合久久久久久| 久久鲁丝午夜福利片| 国产探花极品一区二区| 亚洲最大成人中文| 国产黄片美女视频| 日本-黄色视频高清免费观看| 国内精品一区二区在线观看| 日韩精品中文字幕看吧| 国语自产精品视频在线第100页| 亚洲色图av天堂| 免费av不卡在线播放| 女同久久另类99精品国产91| 一区二区三区高清视频在线| av在线观看视频网站免费| 精品一区二区三区视频在线观看免费| .国产精品久久| 日本成人三级电影网站| 国产毛片a区久久久久| 成人永久免费在线观看视频| 尾随美女入室| 国产伦一二天堂av在线观看| 日韩欧美在线乱码| 中文字幕av在线有码专区| 成人一区二区视频在线观看| 国产欧美日韩精品一区二区| 亚洲欧美成人精品一区二区| av黄色大香蕉| 啦啦啦啦在线视频资源| 最新在线观看一区二区三区| 中文字幕免费在线视频6| 欧美国产日韩亚洲一区| 别揉我奶头 嗯啊视频| 天堂网av新在线| 国产亚洲精品av在线| 亚洲无线观看免费| 91在线观看av| 国产免费男女视频| 黄色一级大片看看| 亚洲第一区二区三区不卡| 久久精品91蜜桃| 国产一级毛片七仙女欲春2| 亚洲真实伦在线观看| 国产私拍福利视频在线观看| 最近视频中文字幕2019在线8| 美女xxoo啪啪120秒动态图| 精品久久久久久久末码| 少妇裸体淫交视频免费看高清| 国产黄a三级三级三级人| 香蕉av资源在线| www.色视频.com| 免费高清视频大片| 成年女人永久免费观看视频| 观看美女的网站| 国模一区二区三区四区视频| 国产精品嫩草影院av在线观看| 色视频www国产| 亚洲天堂国产精品一区在线| 欧美国产日韩亚洲一区| 亚洲不卡免费看| 在线播放国产精品三级| 一级毛片aaaaaa免费看小| 亚洲色图av天堂| 国产精品福利在线免费观看| 亚洲精品在线观看二区| 日本色播在线视频| 搡女人真爽免费视频火全软件 | 成人高潮视频无遮挡免费网站| 欧美绝顶高潮抽搐喷水| 国产免费男女视频| 少妇熟女aⅴ在线视频| 伦理电影大哥的女人| 黄色日韩在线| 国产免费一级a男人的天堂| av中文乱码字幕在线| 欧美日韩精品成人综合77777| av国产免费在线观看| 久久久成人免费电影| 亚洲精品粉嫩美女一区| 久久久久久大精品| 日日摸夜夜添夜夜爱| 亚洲四区av| 国产麻豆成人av免费视频| 久久精品国产99精品国产亚洲性色| 国内久久婷婷六月综合欲色啪| 国产aⅴ精品一区二区三区波| 国产在线男女| 伦精品一区二区三区| 伦理电影大哥的女人| 精品熟女少妇av免费看| 蜜臀久久99精品久久宅男| 亚洲国产精品sss在线观看| 国产成人91sexporn| 精品乱码久久久久久99久播| 长腿黑丝高跟| 狠狠狠狠99中文字幕| 国产亚洲av嫩草精品影院| 51国产日韩欧美| 亚洲久久久久久中文字幕| 亚洲人成网站高清观看| 久久久欧美国产精品| 欧美区成人在线视频| 中文字幕av成人在线电影| 一a级毛片在线观看| 日本精品一区二区三区蜜桃| 久久精品综合一区二区三区| 欧美性感艳星| 国产成人影院久久av| 国产成人aa在线观看| 亚洲人成网站在线播| 日韩一区二区视频免费看| 网址你懂的国产日韩在线| ponron亚洲| 久久精品夜夜夜夜夜久久蜜豆| 国产av不卡久久| 欧美在线一区亚洲| 久99久视频精品免费| 色尼玛亚洲综合影院| 看黄色毛片网站| 看十八女毛片水多多多| 久久久久久九九精品二区国产| avwww免费| 男女边吃奶边做爰视频| 久久天躁狠狠躁夜夜2o2o| 麻豆精品久久久久久蜜桃| 男女那种视频在线观看| 国产精品国产高清国产av| 午夜影院日韩av| 亚洲av美国av| 午夜福利视频1000在线观看| 成人鲁丝片一二三区免费| 成人综合一区亚洲| 亚洲人成网站在线播放欧美日韩| 国产精品无大码| 国产高清视频在线观看网站| 黄色视频,在线免费观看| 毛片一级片免费看久久久久| 99热全是精品| 国产大屁股一区二区在线视频| 看免费成人av毛片| 搡老岳熟女国产| 国产成人福利小说| 欧美成人一区二区免费高清观看| 美女 人体艺术 gogo| 中文在线观看免费www的网站| 中国美白少妇内射xxxbb| 午夜福利在线观看免费完整高清在 | 国产高清有码在线观看视频| 国产中年淑女户外野战色| 久久久久精品国产欧美久久久| 免费看光身美女| 在线观看av片永久免费下载| 麻豆国产97在线/欧美| 色播亚洲综合网| 韩国av在线不卡| 国产在线精品亚洲第一网站| 熟妇人妻久久中文字幕3abv| 成年女人看的毛片在线观看| 欧美日韩精品成人综合77777| 国国产精品蜜臀av免费| 99riav亚洲国产免费| 免费看av在线观看网站| 欧美潮喷喷水| 小说图片视频综合网站| 欧美bdsm另类| 欧美性感艳星| 别揉我奶头~嗯~啊~动态视频| 亚洲久久久久久中文字幕| 人人妻人人澡人人爽人人夜夜 | 级片在线观看| 亚洲欧美日韩东京热| 国产乱人视频| 国产精品国产三级国产av玫瑰| 校园人妻丝袜中文字幕| 一进一出抽搐动态| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av天美| 亚洲国产日韩欧美精品在线观看| 欧美日韩综合久久久久久| 在线观看一区二区三区| 观看美女的网站| 97超级碰碰碰精品色视频在线观看| 亚洲天堂国产精品一区在线| 在线观看66精品国产| 欧美日韩精品成人综合77777| 女生性感内裤真人,穿戴方法视频| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 在线天堂最新版资源| 久久精品国产清高在天天线| 听说在线观看完整版免费高清| 日韩欧美精品免费久久| 一级av片app| 亚洲中文日韩欧美视频| 99久久久亚洲精品蜜臀av| 亚洲熟妇中文字幕五十中出| 色播亚洲综合网| 99久国产av精品| 欧美+日韩+精品| 亚洲真实伦在线观看| 九九热线精品视视频播放| 韩国av在线不卡| 搡老熟女国产l中国老女人| 99久久久亚洲精品蜜臀av| 亚洲综合色惰| 日日干狠狠操夜夜爽| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 久久精品91蜜桃| 国产精品久久久久久精品电影| 99国产精品一区二区蜜桃av| 欧美丝袜亚洲另类| 免费观看精品视频网站| 亚洲电影在线观看av| 久久久成人免费电影| 国产黄色小视频在线观看| 久久久久性生活片| 国产探花在线观看一区二区| 高清午夜精品一区二区三区 | 小说图片视频综合网站| 久久精品久久久久久噜噜老黄 | 丝袜喷水一区| 婷婷精品国产亚洲av| 久久人人精品亚洲av| 舔av片在线| 国产中年淑女户外野战色| 国产激情偷乱视频一区二区| 午夜福利在线观看免费完整高清在 | 天堂av国产一区二区熟女人妻| 亚洲第一区二区三区不卡| 国产精品一及| 国产成人影院久久av| 高清毛片免费看| 成人特级黄色片久久久久久久| 色吧在线观看| 亚洲四区av| 免费观看精品视频网站| 亚洲精品亚洲一区二区| 国产一区二区三区在线臀色熟女| 久久精品国产亚洲网站| 欧美高清成人免费视频www| 日韩精品青青久久久久久| 99九九线精品视频在线观看视频| 国产美女午夜福利| 99久久精品一区二区三区| 精品人妻偷拍中文字幕| 欧美一区二区国产精品久久精品| 国产黄色视频一区二区在线观看 | 日韩精品有码人妻一区| 国模一区二区三区四区视频| 成人美女网站在线观看视频| 国产女主播在线喷水免费视频网站 | 亚洲av二区三区四区| 国产不卡一卡二| 我的女老师完整版在线观看| 欧美zozozo另类| 久久久国产成人精品二区| 中文资源天堂在线| 一级a爱片免费观看的视频| 亚洲高清免费不卡视频| 国产久久久一区二区三区| 中文亚洲av片在线观看爽| 亚洲熟妇熟女久久| 精品欧美国产一区二区三| 国产激情偷乱视频一区二区| 国产精品福利在线免费观看| 一区二区三区免费毛片| 亚洲人成网站在线播| .国产精品久久| 91麻豆精品激情在线观看国产| 久久久久国内视频| 亚洲在线自拍视频| 日本-黄色视频高清免费观看| 亚洲熟妇中文字幕五十中出| 国产麻豆成人av免费视频| 亚洲欧美成人精品一区二区| 国产真实伦视频高清在线观看| 国产高清有码在线观看视频| 直男gayav资源| 熟女电影av网| 日韩三级伦理在线观看| 成人性生交大片免费视频hd| 精品久久久久久久人妻蜜臀av| 日本与韩国留学比较| av视频在线观看入口| 真实男女啪啪啪动态图| 老熟妇仑乱视频hdxx| 欧美性猛交╳xxx乱大交人| 成人亚洲精品av一区二区| 变态另类成人亚洲欧美熟女| 99在线人妻在线中文字幕| 欧美3d第一页| 一级黄色大片毛片| 精品午夜福利在线看| 桃色一区二区三区在线观看| 久久鲁丝午夜福利片| 高清日韩中文字幕在线| 日韩亚洲欧美综合| 亚洲第一区二区三区不卡| 露出奶头的视频| 蜜桃久久精品国产亚洲av| 亚洲,欧美,日韩| 久久人人精品亚洲av| 最近中文字幕高清免费大全6| 精品久久久久久成人av| 亚洲欧美精品自产自拍| 亚洲av一区综合| 久久久久久久久久久丰满| 国产高清视频在线播放一区| 我要看日韩黄色一级片| 亚洲国产精品久久男人天堂| 国内精品美女久久久久久| 欧美性猛交╳xxx乱大交人| 国产大屁股一区二区在线视频| 日韩av不卡免费在线播放| 狠狠狠狠99中文字幕| 亚洲第一电影网av| 人妻夜夜爽99麻豆av| 欧美一区二区亚洲| 亚洲人成网站在线播| 人妻夜夜爽99麻豆av| 一级毛片久久久久久久久女| 亚洲久久久久久中文字幕| 免费观看在线日韩| 校园人妻丝袜中文字幕| 亚洲欧美日韩高清在线视频| 2021天堂中文幕一二区在线观| 欧美成人a在线观看| 黄色欧美视频在线观看| 村上凉子中文字幕在线| 一个人观看的视频www高清免费观看| 久久久精品大字幕| 少妇丰满av| 国产精品一区www在线观看| 噜噜噜噜噜久久久久久91| 亚洲精品国产av成人精品 | av天堂在线播放| 午夜福利在线观看吧| 亚洲欧美成人精品一区二区| 日本精品一区二区三区蜜桃| 国产单亲对白刺激| 亚洲第一区二区三区不卡| 日韩欧美免费精品| 天堂影院成人在线观看| 亚洲乱码一区二区免费版| 久久久久九九精品影院| 日日摸夜夜添夜夜添av毛片| av免费在线看不卡| 午夜激情福利司机影院| 三级经典国产精品| 淫秽高清视频在线观看| 色尼玛亚洲综合影院| 久久久午夜欧美精品| 精品久久久久久久末码| 国产精品人妻久久久久久| 免费av毛片视频| 亚洲成a人片在线一区二区| 国产视频一区二区在线看| 欧美一区二区亚洲| 国产精品一区二区免费欧美| 国产午夜精品久久久久久一区二区三区 | 男女做爰动态图高潮gif福利片| 国产美女午夜福利| 精品不卡国产一区二区三区| av在线天堂中文字幕| 国产成人a区在线观看| 国产 一区精品| 简卡轻食公司| av天堂中文字幕网| 偷拍熟女少妇极品色| 国产成人a区在线观看| 婷婷色综合大香蕉| 亚洲人与动物交配视频| 国产高清不卡午夜福利| 欧美一区二区亚洲| 人妻丰满熟妇av一区二区三区| 少妇被粗大猛烈的视频| 99热6这里只有精品| 久久精品91蜜桃| 欧美激情在线99| 看黄色毛片网站| 国产精品免费一区二区三区在线| 国产视频内射| 麻豆久久精品国产亚洲av| 亚洲欧美日韩卡通动漫| 日韩一区二区视频免费看| 久久久成人免费电影| 干丝袜人妻中文字幕| 12—13女人毛片做爰片一| 性插视频无遮挡在线免费观看| 国产亚洲精品综合一区在线观看| 丰满人妻一区二区三区视频av| 晚上一个人看的免费电影| 日韩人妻高清精品专区| 春色校园在线视频观看| 天天一区二区日本电影三级| 国产精品三级大全| 国产 一区 欧美 日韩| 亚洲成人精品中文字幕电影| 成人性生交大片免费视频hd| 99视频精品全部免费 在线| 国产高清有码在线观看视频| 欧美高清成人免费视频www| 综合色av麻豆| 免费av观看视频| 亚洲精品粉嫩美女一区| 中文字幕人妻熟人妻熟丝袜美| 深夜a级毛片| 最近2019中文字幕mv第一页| 永久网站在线| 中国美女看黄片| 婷婷亚洲欧美| 美女xxoo啪啪120秒动态图| 麻豆一二三区av精品| 国产高清视频在线观看网站| 最近最新中文字幕大全电影3| 搡女人真爽免费视频火全软件 | 国产av麻豆久久久久久久| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 欧美bdsm另类| 久久精品国产99精品国产亚洲性色| 嫩草影视91久久| 麻豆成人午夜福利视频| 午夜福利在线观看免费完整高清在 | 国产免费一级a男人的天堂| 日韩高清综合在线| .国产精品久久| 91久久精品国产一区二区成人| 大香蕉久久网| 久久午夜福利片| av福利片在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品女同一区二区软件| 国产av麻豆久久久久久久| 欧美精品国产亚洲| 男女那种视频在线观看| 91麻豆精品激情在线观看国产| 亚洲人成网站在线播| 美女免费视频网站| 日本黄色片子视频| 97碰自拍视频| 成人性生交大片免费视频hd| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 女人十人毛片免费观看3o分钟| 黄色一级大片看看| 欧美zozozo另类| 色综合色国产| av视频在线观看入口| 搞女人的毛片| 一夜夜www| 成人永久免费在线观看视频| 黄色视频,在线免费观看| 真人做人爱边吃奶动态| 床上黄色一级片| 成人一区二区视频在线观看| 国产高清视频在线观看网站| 日韩大尺度精品在线看网址| 亚洲av美国av| 精品久久久久久久久久免费视频| 丰满的人妻完整版| 午夜久久久久精精品| 十八禁网站免费在线| 午夜免费男女啪啪视频观看 | 久久久国产成人精品二区| 自拍偷自拍亚洲精品老妇| 韩国av在线不卡| 国产三级中文精品| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 国产片特级美女逼逼视频| 两个人的视频大全免费| 国产人妻一区二区三区在| 成人av一区二区三区在线看| 亚洲成人久久爱视频| 日韩人妻高清精品专区| 精品无人区乱码1区二区| 国产真实伦视频高清在线观看| 99久国产av精品| 久久精品久久久久久噜噜老黄 | 男女啪啪激烈高潮av片| 国产美女午夜福利| 男人舔奶头视频| 在线看三级毛片| 色在线成人网| 欧美成人精品欧美一级黄| 国产 一区精品| 晚上一个人看的免费电影| 国产高清激情床上av|