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

    過阻尼布朗棘輪的斯托克斯效率研究?

    2017-12-05 02:35:16呂明濤延明月艾保全高天附鄭志剛
    物理學報 2017年22期
    關(guān)鍵詞:斯托克斯棘輪布朗

    呂明濤 延明月 艾保全 高天附 鄭志剛

    1)(沈陽師范大學物理科學與技術(shù)學院,沈陽 110034)2)(華南師范大學物理與電信工程學院,廣州 510006)3)(華僑大學信息科學與工程學院,廈門 361021)

    過阻尼布朗棘輪的斯托克斯效率研究?

    呂明濤1)延明月1)艾保全2)高天附1)?鄭志剛3)

    1)(沈陽師范大學物理科學與技術(shù)學院,沈陽 110034)2)(華南師范大學物理與電信工程學院,廣州 510006)3)(華僑大學信息科學與工程學院,廈門 361021)

    (2017年6月27日收到;2017年8月27日收到修改稿)

    根據(jù)隨機能量理論解析得到阻尼環(huán)境中布朗粒子的概率流和斯托克斯效率,并進一步研究布朗粒子的輸運性能.詳細討論了空間的不對稱性、外偏置力及外勢結(jié)構(gòu)等對棘輪定向輸運的影響.研究發(fā)現(xiàn),合適的外偏置力能使棘輪的定向輸運達到最強.通過調(diào)節(jié)外勢的不對稱性可使棘輪中粒子的運動反向,當選擇合適的空間不對稱性時布朗粒子的反向輸運可獲得最強.此外,一定條件下合適的外勢高度也能增強棘輪輸運,且能使粒子克服黏滯阻力的斯托克斯效率達到最大.所得結(jié)論能夠啟發(fā)實驗上設(shè)計合適的外勢及外偏置來優(yōu)化布朗棘輪的定向輸運性能,并為生物納米器件的研制提供一定的理論參考.

    概率流,黏滯阻力,斯托克斯效率,輸運性能

    1 引 言

    生命體內(nèi)部的細胞活動中分子馬達扮演了重要角色.由于馬達能夠把化學能直接轉(zhuǎn)化為機械能,所以這一獨特的能量轉(zhuǎn)換形式對細胞內(nèi)的物質(zhì)輸送至關(guān)重要.為了研究分子馬達的定向輸運,人們提出了各種模型[1],包括物理上簡單的布朗棘輪模型[2?4].通過研究不同棘輪中布朗粒子的運動情況,不但能了解布朗棘輪產(chǎn)生定向輸運的機制,找到其定向輸運的最優(yōu)條件,而且能夠提高棘輪對輸入到系統(tǒng)中能量的有效利用.

    對于不同棘輪的研究,如何增強棘輪的定向輸運一直廣受學者的關(guān)注[4?7].早期研究中,如Bao等[8]發(fā)現(xiàn)了在黏滯阻力及雙噪聲環(huán)境下布朗粒子的定向輸運可被優(yōu)化,說明黏滯阻力在一定條件下能夠增強布朗粒子的定向輸運.同時,Ai等[9]根據(jù)Fick-Jacobs方程解析得到了二維通道中布朗粒子概率流的表達式,并討論了棘輪熱泵的性能.此外,我們早期也研究了不同時空不對稱度下布朗馬達的定向輸運性能,發(fā)現(xiàn)通過改變外勢的不對稱性粒子的輸運能夠增強[10].此外,文獻[11]最新研究發(fā)現(xiàn),改變外力的變化頻率能夠增強粒子的概率流,且當棘輪輸入合適的外力時棘輪的能量利用率能夠獲得提高.

    隨著棘輪定向輸運問題的深入研究,如何有效利用棘輪系統(tǒng)的能量已引起人們的極大興趣.Sekimoto等[12]首先根據(jù)隨機能量理論提出了不同環(huán)境中棘輪能量的分析方法,并計算了三種不同情況下棘輪的輸入、輸出功.有趣的是,卓益忠等[8]研究發(fā)現(xiàn)通過外偏置力的作用,布朗馬達的能量利用率能夠得到提高.此外大量實驗表明,布朗馬達拖動負載做功時會降低棘輪系統(tǒng)的輸入能[10].同時,實驗結(jié)果還表明,生物體內(nèi)仍存在一類沒有拖動負載卻能產(chǎn)生定向運動的生物分子馬達,并且這類馬達克服阻尼環(huán)境運動時的定向輸運效率非常高[13].因此,研究無負載下馬達的能量利用問題也至關(guān)重要.然而這種情況下,若仍采用傳統(tǒng)熱力學效率的研究方法分析馬達的能量轉(zhuǎn)化問題,該計算方法將失去其物理意義.因此,無負載條件下布朗粒子的定向輸運效率問題目前仍是人們研究的熱點.理論上,根據(jù)Wang的分析,無負載情況下布朗粒子克服黏滯阻力運動時的輸運效率稱為斯托克斯效率[14].目前關(guān)于斯托克斯效率雖有一定研究[11,14],但并不是十分深入.如文獻[15]中研究了時空不對稱勢下棘輪的輸運特性,并對棘輪的熱力學效率和定向輸運效率(斯托克斯效率)分別進行了討論.此外也對耦合布朗粒子的能量轉(zhuǎn)化問題做了初步探討,發(fā)現(xiàn)在一定條件下耦合粒子拖動負載做功時的能量轉(zhuǎn)化效率能夠獲得提高[10].同時文獻[11]對搖擺棘輪中概率流和斯托克斯效率的隨機行為進行了討論.在上述棘輪效率的研究中,人們大都采用數(shù)值方法研究不同棘輪的定向輸運及能量轉(zhuǎn)換問題.然而關(guān)于斯托克斯效率的解析求解至今還不夠深入.因此,本文進一步采用解析方法討論棘輪勢中布朗粒子的定向輸運效率問題.

    根據(jù)Sekimoto[12]隨機能量理論,研究了黏滯阻力作用下布朗粒子的定向輸運及斯托克斯效率問題[16?20].研究發(fā)現(xiàn),在一定不對稱性下隨外偏置力振幅的變化粒子可產(chǎn)生反向運動,且存在合適的外力振幅使概率流達到最強.此外,本文還發(fā)現(xiàn)合適的外勢結(jié)構(gòu)能夠使粒子克服黏滯環(huán)境時的定向輸運效率達到最強.應(yīng)用本文所得結(jié)論實驗上可以啟發(fā)人們設(shè)計合適的外勢及外偏置來增強棘輪的定向輸運性能.

    2 模 型

    本文研究了處在過阻尼環(huán)境下的布朗粒子在空間不對稱周期勢中的運動,同時粒子受周期外力以及熱噪聲的作用,其動力學行為可由如下過阻尼朗之萬方程描述[21,22]:

    其中,x(t)表示粒子t時刻的位置,γ為介質(zhì)阻尼系數(shù),ξ(t)為高斯白噪聲,且滿足〈ξ(t)〉=0,〈ξi(t)ξj(t′)〉=2Dδijδ(t?t′), δ(t)是高斯噪聲的關(guān)聯(lián)函數(shù),D表示噪聲強度,F(t)為t時刻布朗粒子受到的非對稱周期外力:

    上式中A為外力振幅,2π/?為外力周期.此外,U(x)為周期外勢,如圖1所示,其表達式為

    其中a為周期外勢高度,L為外勢的周期且滿足U(x)=U(x+L),?表示外勢的不對稱性,因此通過調(diào)節(jié)?可改變外勢的結(jié)構(gòu).

    圖1 外勢U(x)的結(jié)構(gòu)示意圖Fig.1.The schematic diagram of U(x).

    為討論方便,令V′(x)=?U(x)/?x,根據(jù)文獻[23]理論,與方程(1)對應(yīng)的Fokker-Planck方程和概率密度流守恒方程為

    其中kB是玻爾茲曼常數(shù),T為熱浴溫度,且噪聲強度D=kBT,P(x,t)表示在t時刻在位置x處的概率密度,并滿足P(x,t)=P(x+L,t),j(x,t)為概率流,由上式可得:

    如果F(t)變化很慢即系統(tǒng)為準靜態(tài),則

    根據(jù)上式和已知條件可得系統(tǒng)概率流表達式:

    其中

    由于外力F(t)變化很慢,所以概率流也可表示為[23,25,26]

    根據(jù)Sekimoto隨機能量理論[12],釋放到熱浴中的能量為

    其中ti和tf為積分始末時刻,且積分區(qū)間τ=ti?tf,單位時間內(nèi)勢能的變化為?U,由能量守恒關(guān)系R=?U+E,輸入到系統(tǒng)中的總能量R可表述為

    根據(jù)dU=(?U/?x)dx+(?U/?t)dt, 上式化簡為[12]

    由于外勢的變化與時間t無關(guān),?U(x,t)/?t=0,上式可寫為[23,25]

    在準靜態(tài)下,根據(jù)文獻[24]理論,上式可進一步化簡:

    根據(jù)Mamata理論,斯托克斯效率定義為[11]

    由(16)式可知輸入到系統(tǒng)的能量:

    粒子克服阻尼力做功[15]:

    其中平均速度

    所以斯托克斯效率的解析表達式為

    將(11),(17),(20)及(21)式代入(22)式可得:

    根據(jù)(11)及(23)式,本文主要研究阻尼環(huán)境中布朗棘輪的輸運性能隨系統(tǒng)各參量的變化關(guān)系,且所有物理量已采用無量綱化處理.如無特殊說明,文中參量取L=1.0,γ=1.0,D=0.01,?=1.0,a=1/2π,A=1.0,?=2π.

    3 結(jié)果與討論

    3.1 外力振幅A的影響

    為了解棘輪的定向輸運,研究了概率流隨各參量的變化.首先討論了當外勢不對稱度較小時,即?=1.0時不同噪聲強度下概率流J隨外偏置力振幅A的變化關(guān)系,如圖2(a)所示.當A較小時,概率流為負,且隨A的增加棘輪產(chǎn)生的反向概率流逐漸增大.此時由于外偏置力振幅A較小,所以棘輪不容易產(chǎn)生定向輸運.隨著A的繼續(xù)增加,反向概率流均存在峰值,該現(xiàn)象表明合適的外偏置力振幅能夠促進棘輪的定向輸運.當A繼續(xù)增加時,|J|逐漸減小并趨于零.然而當A大于某一臨界值且增大時,J>0,此時概率流與外偏置力呈類似正比關(guān)系,如圖2(a)插圖所示.說明此時外偏置力較大并占主導(dǎo)作用,因此粒子產(chǎn)生的定向輸運更強.研究發(fā)現(xiàn)隨噪聲強度D的增加,概率流峰值逐漸減小,該現(xiàn)象表明較大的噪聲并不利于粒子的定向輸運.可見除外偏置力,噪聲強度仍能影響粒子的定向輸運.下文將深入討論噪聲強度對概率流的影響.

    圖2 (網(wǎng)刊彩色)?=1.0時(a)概率流J和(b)斯托克斯效率ηs隨外力振幅A的變化Fig.2.(color online)The curves of(a)particle flow J,(b)Stokes efficiency ηsvarying with the amplitude of external force A,where?=1.0.

    由(22)式分析可知,斯托克斯效率ηs反映了布朗粒子克服黏滯阻力產(chǎn)生定向輸運的性能.因此本文研究了ηs隨外偏置力振幅A的變化關(guān)系,如圖2(b)所示.結(jié)果表明隨著A的增加,不同噪聲強度下的斯托克斯效率均逐漸增大且出現(xiàn)峰值.以D=0.02為例,發(fā)現(xiàn)在圖2(a)和圖2(b)中兩組曲線的變化趨勢類似.例如在A=1.0附近,反向概率流最大時相對應(yīng)的斯托克斯效率也接近最大.根據(jù)(22)和(23)式可知,在一定條件下概率流越大,斯托克斯效率也越大,因此當過阻尼棘輪的定向輸運達到最強時,相對應(yīng)地粒子克服黏滯阻力的斯托克斯效率也會增強.此外研究還發(fā)現(xiàn)隨D的增加ηs峰值也會減小,說明過大的熱噪聲不利于棘輪定向輸運效率的提高,下文將繼續(xù)討論熱噪聲對棘輪輸運性能的影響.

    3.2 熱浴環(huán)境中噪聲強度D的影響

    由于噪聲強度對粒子輸運有一定影響,研究了不同不對稱度?下概率流J隨熱噪聲強度D的變化關(guān)系,如圖3(a)所示.研究結(jié)果表明隨D的增加,概率流也會小于零,說明在一定條件下熱噪聲可以使棘輪產(chǎn)生反向輸運,該結(jié)論與文獻[27,28]中所得結(jié)論類似,但產(chǎn)生機制并不完全相同.當D繼續(xù)增加,棘輪存在合適的噪聲強度使概率流達到最強.這一現(xiàn)象說明合適的熱噪聲也能增強棘輪的定向輸運.隨著D的進一步增大,|J|→0.說明當溫度較高時,粒子在做自由擴散運動,所以此時棘輪不容易產(chǎn)生定向輸運.有趣的是隨不對稱度?的增大,反向概率流的峰值逐漸減小.該現(xiàn)象表明在一定條件下,較大的外勢不對稱度會抑制粒子的定向運動.通過(3)式分析,?越大外勢高度越大,因此粒子較難跨越勢壘產(chǎn)生定向輸運.由圖2(a)可知,?=1.0時,概率流可為負.為了深入研究粒子的輸運方向問題,進一步研究了?較小時,如?=0.1,0.01時粒子的輸運行為,如3(a)插圖所示.研究發(fā)現(xiàn),?較小時概率流可為正,說明棘輪的定向輸運對外勢的不對稱度有較強的敏感性,因此當?較大時圖3(a)中J-D會呈現(xiàn)反向概率流.因此該現(xiàn)象與文獻[27]中噪聲誘導(dǎo)棘輪產(chǎn)生流反轉(zhuǎn)的機制并不完全相同.

    如上述討論,噪聲強度能增強棘輪的定向輸運.因此進一步研究了斯托克斯效率ηs隨溫度D的變化關(guān)系,如圖3(b)所示.結(jié)果表明,在一定條件下存在合適的噪聲強度使斯托克斯效率ηs達到最大.對照圖3不難發(fā)現(xiàn),J-D圖像和ηs-D圖像均出現(xiàn)峰值.這是因為合適的熱噪聲強度能夠增強粒子輸運并能提高棘輪的定向輸運效率,所以布朗粒子克服黏滯阻力產(chǎn)生定向輸運的性能也會獲得提高.此外研究還發(fā)現(xiàn)隨著外勢不對稱性的增大,斯托克斯效率的峰值逐漸減小.同樣由(3)式可知這是由于?越大,外勢高度a越大,導(dǎo)致粒子不容易跨越外勢產(chǎn)生定向運動.因此粒子克服黏滯阻力運動時的性能也被減弱.

    圖3 (網(wǎng)刊彩色)A=1.0時(a)概率流J和(b)斯托克斯效率ηs隨熱噪聲強度D的變化Fig.3.(color online)The curves of(a)particle flow J,(b)Stokes efficiency ηsvarying with the noise intensity D,where A=1.0.

    3.3 外勢不對稱參數(shù)?的影響

    為了深入了解外勢不對稱性對棘輪輸運性能的影響,本文研究了不同外力振幅下概率流J隨外勢不對稱度?的變化關(guān)系,如圖4(a)所示.當?較小時概率流為正且逐漸減小.當?達到某一臨界值時,即?=?c(如A=1.6),J=0.隨著?增加J<0且開始反向增大.當?繼續(xù)增大到另一臨界值?=?m時,棘輪的反向輸運最強.最終隨著?增加,|J|→0.產(chǎn)生上述輸運現(xiàn)象主要是由于當外勢不對稱度?<?c時,?越大,外勢高度越大,粒子越難跨越外勢產(chǎn)生正向輸運;當?>?c時,概率流開始反向增加,這是由于隨?的繼續(xù)增加,粒子將開始反向運動;有趣的是當?=?m時反向概率流最大,說明在一定條件下,合適的外勢不對稱性可增強粒子的定向輸運;當?>?m時|J|→0,這是由于較大的外勢不對稱度占主導(dǎo)作用,導(dǎo)致粒子很難形成定向流.此外研究還發(fā)現(xiàn),在一定?范圍內(nèi),J>0時的概率流與外偏置力振幅呈正比例變化,該結(jié)論與圖2(a)中?=1.0時的結(jié)論一致.因此只有在合適的外勢不對稱性和外偏置力共同作用下,棘輪的定向輸運才能增強.

    圖4 (網(wǎng)刊彩色)D=0.01時(a)概率流J和(b)斯托克斯效率ηs隨不對稱參量?的變化Fig.4.(color online)The curves of(a)particle flow J,(b)Stokes efficiency ηsvarying with the barrier asymmetric parameter?,where D=0.01.

    由于棘輪輸運對外勢不對稱度有一定敏感性,討論了斯托克斯效率ηs隨外勢不對稱度?的變化,如圖4(b)所示.隨?的增加不同外偏置力振幅A下的ηs均逐漸減小,如A=1.6,當?增加到某一臨界值?=?c時,ηs→0.隨?繼續(xù)增加,ηs開始逐漸增大,當外勢的不對稱度取某一合適值時,即?=?m,ηs達到峰值.研究發(fā)現(xiàn),?<?c時ηs→0的現(xiàn)象主要是由于在圖4(a)中相對應(yīng)的J隨不對稱度的增加也逐漸減小,根據(jù)類似的討論,由于ηs與概率流呈類似正比關(guān)系,所以ηs逐漸降低并趨于零;斯托克斯效率在?=?m處達到峰值說明合適的外勢不對稱度能夠增強棘輪的定向輸運性能;當?>?m,ηs逐漸減小并趨于零是由于外勢的不對稱度過大,外勢過高,粒子更不容易產(chǎn)生定向運動,所以粒子克服阻尼力的定向輸運效率減小并趨于零.此外研究還發(fā)現(xiàn),ηs的峰值隨外偏置力振幅的增大而減小,這一結(jié)論仍與圖2(b)結(jié)論一致.

    3.4 外勢高度a的影響

    由(3)式分析可知,外勢不對稱性可通過改變外勢的高度影響粒子的定向輸運.因此進一步研究了概率流隨外勢高度a的變化關(guān)系,如圖5(a)所示.研究發(fā)現(xiàn)隨a的增加概率流逐漸減小,當a達到某一臨界值,a=ac時,J=0.當a繼續(xù)增加J<0且反向增大,并在a=am處反向概率流達到峰值.最終隨著a的增加,J→0.結(jié)果表明,以A=3.0為例,當a<ac時外勢高度越小,粒子越容易跨越外勢產(chǎn)生正向輸運;當a>ac時J<0,說明當外勢高度達到某一臨界值后,粒子將反向輸運;有趣的是概率流在a=am處達到峰值,說明外勢高度仍存在一個合適值使粒子的輸運達到最強;當a>am后概率流減小且|J|→0,這是由于外勢過高,布朗粒子很難跨越外勢產(chǎn)生定向輸運.此外,隨著外偏置力振幅A的增加,反向概率流的峰值也逐漸增大,說明一定條件下外偏置力越大,越促進粒子的定向運動,該結(jié)論與圖2(a)中結(jié)論一致.因此,為了增強棘輪的定向輸運,實驗上可設(shè)計合適的外勢結(jié)構(gòu)及外偏置力使布朗粒子的定向輸運達到最強[27].

    由上述分析可知,外勢高度能夠影響粒子輸運.因此進一步研究了不同外偏置力振幅下ηs隨a的變化關(guān)系,如圖5(b)所示.研究發(fā)現(xiàn)ηs隨a的變化仍會出現(xiàn)極值,即在波峰a=am處和波谷a=ac處.結(jié)果表明,仍以A=3.0為例,a<ac時斯托克斯效率逐漸減小,說明一定條件下a越大粒子克服黏滯阻力的定向輸運效率越低;當a>ac時ηs在a=am處出現(xiàn)峰值,說明合適的外勢高度能促使粒子產(chǎn)生較強的定向輸運,從而使粒子克服介質(zhì)阻尼力的輸運效率達到最強;然而當a>am時ηs逐漸減小并趨于零,該現(xiàn)象表明當外勢較高時,粒子想要跨過外勢形成定向輸運的能力減小并逐漸消失.此外,研究還發(fā)現(xiàn)隨外偏置力振幅的增加,粒子定向輸運效率的峰值逐漸增大.該結(jié)論能夠啟發(fā)人們設(shè)計合適的外勢結(jié)構(gòu)來優(yōu)化棘輪的定向輸運性能.

    圖5 (網(wǎng)刊彩色)D=0.1時(a)概率流J和(b)斯托克斯效率隨外勢高度a的變化Fig.5.(color online)The curves of(a)Particle flow J,(b)Stokes efficiency ηsvarying with the barrier height a,where D=0.1.

    4 結(jié) 論

    研究了阻尼環(huán)境中處在時變外偏置力以及噪聲作用下布朗棘輪的定向輸運.深入分析了外勢的空間不對稱度、外偏置力振幅、外勢高度及熱噪聲強度對棘輪定向輸運的影響.研究結(jié)果表明,在合適的外勢結(jié)構(gòu)下,布朗粒子能夠產(chǎn)生反向運動.同時存在合適的空間不對稱度使概率流達到最強,并且當外偏置力振幅為某一合適值時,棘輪定向輸運的能力能夠獲得提高.除此之外,合適的外勢高度也能夠使粒子克服黏滯阻力時的定向輸運效率達到最強,從而有效地提高了棘輪的輸運性能.本文所得結(jié)論不但能夠啟發(fā)人們設(shè)計不同類型的棘輪以增強粒子的定向輸運并提高布朗粒子的輸運效率,還可為醫(yī)學上的藥物輸送及納米器件的控制提供一定的理論指導(dǎo)作用[29?32].

    [1]Xie P 2010Int.J.Biol.Sci.6 665

    [2]Doering C R 1995Nuovo Cimento17 685

    [3]Astumian R D,Bier M 1994Phys.Rev.Lett.72 1766

    [4]Gao T F,Chen J C 2009J.Phys.A:Math.Theor.42 065002

    [5]Ai B Q,He Y F,Zhong W R 2011Phys.Rev.E83 051106

    [6]Zhang H W,Wen S T,Zhang H T,Li Y X,Chen G R 2012Chin.Phys.B21 078701

    [7]Gao T F,Liu F S,Chen J C 2012Chin.Phys.B21 020502

    [8]Zhan Y,Bao J D,Zhuo Y Z 1997Acta Phys.Sin.46 1880(in Chinese)[展永,包景東,卓益忠 1997物理學報461880]

    [9]Ai B Q,He Y F,Li F G,Zhong W R 2013Phys.Rev.E138 154107

    [10]Fan L M,Lü M T,Huang R Z,Gao T F,Zheng Z G 2017Acta Phys.Sin.66 010501(in Chinese)[范黎明,呂明濤,黃仁忠,高天附,鄭志剛2017物理學報66 010501]

    [11]Sahoo M,Jayannavar A M 2017Physica A465 40

    [12]Sekimoto K 1997J.Phys.Soc.Jpn.66 1234

    [13]Parrondo J M R,Cisneros B J D 2002Physics A75 179

    [14]Wang H,Oster G 2002EPL57 134

    [15]Li Y X,Wu X Z,Zhuo Y Z 2000Physica A286 147

    [16]Chueshov I,Kuksin S 2008Physica D237 1352

    [17]Winkler M,Abel M 2015Phys.Rev.E92 063002

    [18]Anders M 2013Phys.Rev.E92 063002

    [19]Sztuk E,Przekop R,Gradoń L 2012Chem.Process Eng.33 279

    [20]Spiechowicz J,Luczka J,Machura L 2016Physics2016 054038

    [21]Zheng Z G,Cross M C,Hu G 2002Phys.Rev.Lett.89 154102

    [22]Li C P,Han Y R,Zhan Y,Hu J J,Zhang L G,Qu J 2013Acta Phys.Sin.62 230051(in Chinese)[李晨璞,韓英榮,展永,胡金江,張禮剛,曲蛟 2013物理學報 62 230051]

    [23]Zeng C H,Wang H 2012Chin.Phys.B21 76

    [24]Ai B Q 2004Ph.D.Dissertation(Guangzhou:Sun Yatsen University)(in Chinese)[艾保全 2004博士學位論文(廣州:中山大學)]

    [25]Kamegawa H,Hondou T,Takagi F 1998Phys.Rev.Lett.80 5251

    [26]Ai B Q,Xie H Z,Liao H Y,Liu L G 2006J.Stat.Mech.50 09016

    [27]Bartussek R,Hanggi P,Lindner B,Schimansky Geier L 1997Physica D109 17

    [28]Kula J,Czernik T,Luczka J 1998Phys.Rev.Lett.80 1377

    [29]Linke H 2002Appl.Phys.A:Mater.Sci.Process.75 167

    [30]Heuvel MGLVD,Dekker C 2007Science317 333

    [31]Braun O M,Kivshar Y S 1998Phys.Rep.306 1

    [32]Landa P S,McClintock P V E 2000Phys.Rep.323 1

    PACS:05.40.–a,05.40.Jc,02.30.Yy,02.50.–rDOI:10.7498/aps.66.220501

    *Project supported by the National Natural Science Foundation of China(Grand Nos.11475022,11347003),the Scienti fic Research Funds of Huaqiao University,China and the Excellent Talents Program of Shenyang Normal University,China(Grand No.91400114005),and the Teaching Research Program of Thermodynamics and Statistical Physics in the Institution of Higher Education,China.

    ?Corresponding author.E-mail:tianfugao@synu.edu.cn

    Stokes efficiency in the overdamped Brownian ratchet?

    Lü Ming-Tao1)Yan Ming-Yue1)Ai Bao-Quan2)Gao Tian-Fu1)?Zheng Zhi-Gang3)

    1)(College of Physical Science and Technology,Shenyang Normal University,Shenyang 110034,China)2)(School of Physics and Telecommunication Engineering,South China Normal University,Guangzhou 510006,China)3)(College of Information Science and Engineering,Huaqiao University,Xiamen 361021,China)

    27 June 2017;revised manuscript

    27 August 2017)

    Molecular motors in life activities of cell are known to operate efficiently.They could convert molecular-scale chemical energy into macroscopic-scale mechanical work with high efficiency.In order to acquire the transport mechanism of the molecular motor,the Brownian ratchet has been proposed to explore the property of directed transport and energy conversion.There are di ff erent kinds of Brownian ratchet models like flashing ratchets,rocking ratchets,and time-asymmetric ratchets and so on.Through investigating the performance of Brownian ratchet moving in periodic potential,the directed transport of ratchet could be explained,and the e ff ective usage of ratchet energy for directed transport could also be improved.Recently,optimizing the transport of Brownian ratchet has aroused the interest of researchers.It is found that the viscous resistance could reinforce the directed transport of the Brownian particle in damping liquid.Meanwhile,a large number of conclusions indicate that the transport of Brownian ratchets would be enhanced if the asymmetry of the potential is changed.Those results show that the in fluences of the external potential and the damping force on the particle flow cannot be neglected.Hence in this paper,the e ff ects of the potential structure and the temperature of heat bath on transport are discussed.

    Furthermore,how to use the ratchet energy e ff ectively has been investigated in recent years.When the Brownian motor operates with load,the input energy is reduced.More importantly,the energy transformation efficiency de fined as the ratio of the useful work done against the load to the input energy is assumed to be a zero value in the absence of load.With the help of stochastic energetic theory proposed by Sekimoto,the Stokes efficiency has been used to explore the performance of the Brownian ratchet.So far,the numerical solution has been used extensively in most theoretical investigations.Nevertheless,in our work,the Stokes efficiency is discussed analytically for explaining the mechanism of directed transport.We consider the transport performance of the Brownian ratchet described by the Fokker Planck equation which is corresponding to the Langevin equation under time-varying external force and thermal noise.Mainly,the e ff ects of potential asymmetry,external force,height of the barrier,and intensity of the thermal noise on transport are discussed in detail.It is found that the transport direction of Brownian ratchet will be reversed under the condition of appropriate potential structures,and the probability current can reach a maximal value by changing the asymmetry of potential.It is worthwhile to point out that the performance of directed transport of the ratchet can be improved when an appropriate amplitude of the external force is applied.Meanwhile,there is an optimal value of the barrier height at which the Stokes efficiency reaches a maximal value and the directed transport of ratchet is enhanced.Through our conclusions,the ratchets of di ff erent structures could be designed for improving the transport property of Brownian motor.And the results have helpful theoretical guidance not only in the aspect of medical delivery but also in the control of nano-devices.

    probability current,viscous resistance,Stokes efficiency,transport property

    10.7498/aps.66.220501

    ?國家自然科學基金(批準號:11475022,11347003)、華僑大學科研啟動費項目和沈陽師范大學優(yōu)秀人才支持計劃 (批準號:91400114005)和高等學校熱力學與統(tǒng)計物理課程教學研究項目資助的課題.

    ?通信作者.E-mail:tianfugao@synu.edu.cn

    猜你喜歡
    斯托克斯棘輪布朗
    基于拉曼散射光動態(tài)校準的分布式光纖溫度傳感系統(tǒng)*
    物理學報(2022年20期)2022-10-27 02:58:36
    Legendary British Climber Joe Brown喬·布朗
    基于抽運-探測法的皮秒反斯托克斯拉曼頻移器的理論研究*
    物理學報(2021年9期)2021-05-14 02:44:50
    基于A-V模型的改進模型及單軸棘輪效應(yīng)預(yù)測
    向前!向前!
    你好,我是布朗熊
    你好,我是布朗熊
    學生天地(2018年30期)2018-10-17 01:33:50
    丹·布朗主要作品
    當代作家(2018年8期)2018-08-31 02:48:44
    基于相干反斯托克斯拉曼散射的二維溫度場掃描測量
    棘輪下錨補償裝置斷線制動研究與探討
    電氣化鐵道(2015年2期)2015-06-29 10:05:39
    青春草国产在线视频 | 99热只有精品国产| 亚洲精品日韩在线中文字幕 | 神马国产精品三级电影在线观看| 变态另类成人亚洲欧美熟女| 午夜精品在线福利| 亚洲18禁久久av| 国产精品国产三级国产av玫瑰| 色综合色国产| 国产精华一区二区三区| 亚洲激情五月婷婷啪啪| 欧美丝袜亚洲另类| 91久久精品电影网| 国产精品蜜桃在线观看 | 欧美日韩乱码在线| 欧美日本视频| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 插逼视频在线观看| 哪里可以看免费的av片| 亚洲久久久久久中文字幕| 日本欧美国产在线视频| 女同久久另类99精品国产91| 欧美日韩国产亚洲二区| 国产综合懂色| 久久6这里有精品| 亚洲精品久久久久久婷婷小说 | 日韩欧美精品免费久久| 日韩强制内射视频| 大型黄色视频在线免费观看| 亚洲av中文av极速乱| 亚洲人成网站在线播放欧美日韩| 在线观看免费视频日本深夜| 免费无遮挡裸体视频| 黄色配什么色好看| 搞女人的毛片| 桃色一区二区三区在线观看| 精品免费久久久久久久清纯| 国产蜜桃级精品一区二区三区| h日本视频在线播放| 99热这里只有精品一区| 成人二区视频| 国产精品国产三级国产av玫瑰| 国产国拍精品亚洲av在线观看| 日本色播在线视频| 国产黄色小视频在线观看| 日韩亚洲欧美综合| 亚洲av免费高清在线观看| 亚洲人成网站在线播| 国产精品不卡视频一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲乱码一区二区免费版| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一级毛片在线| 国产精品电影一区二区三区| 小说图片视频综合网站| 久久久成人免费电影| 午夜福利在线在线| 国产国拍精品亚洲av在线观看| 日韩欧美精品免费久久| 国产日本99.免费观看| 国产又黄又爽又无遮挡在线| 亚洲欧美日韩高清在线视频| 亚洲一区二区三区色噜噜| 少妇丰满av| 色综合站精品国产| 精品一区二区三区视频在线| 亚洲国产精品国产精品| 成年版毛片免费区| 亚洲精品成人久久久久久| 一边亲一边摸免费视频| 久久午夜亚洲精品久久| 国产人妻一区二区三区在| 精品久久久久久成人av| 一区二区三区免费毛片| 最近视频中文字幕2019在线8| 久久综合国产亚洲精品| 男人狂女人下面高潮的视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产精品sss在线观看| 国产精品蜜桃在线观看 | 91久久精品电影网| 国产精品女同一区二区软件| 91精品国产九色| 老熟妇乱子伦视频在线观看| 波多野结衣高清作品| 日韩亚洲欧美综合| 中文资源天堂在线| 亚洲综合色惰| 欧美色视频一区免费| 国产激情偷乱视频一区二区| 在线观看av片永久免费下载| 免费观看人在逋| 国产精品伦人一区二区| 插逼视频在线观看| 激情 狠狠 欧美| 最近2019中文字幕mv第一页| 国产亚洲精品av在线| 国产一区二区三区av在线 | 精品一区二区三区人妻视频| 深爱激情五月婷婷| 非洲黑人性xxxx精品又粗又长| 国产精品麻豆人妻色哟哟久久 | 国产乱人视频| 亚洲欧美精品综合久久99| 看黄色毛片网站| 亚洲av电影不卡..在线观看| 又爽又黄无遮挡网站| 又黄又爽又刺激的免费视频.| 乱人视频在线观看| 九九热线精品视视频播放| 精品一区二区三区视频在线| 久久99热6这里只有精品| 99精品在免费线老司机午夜| 成人亚洲欧美一区二区av| av在线观看视频网站免费| 一边亲一边摸免费视频| 亚洲国产色片| 午夜精品国产一区二区电影 | 国产成人91sexporn| 免费一级毛片在线播放高清视频| 日本黄色片子视频| 免费黄网站久久成人精品| 免费av观看视频| 成年女人看的毛片在线观看| 婷婷色综合大香蕉| 中文字幕熟女人妻在线| 给我免费播放毛片高清在线观看| 成年av动漫网址| 日韩成人伦理影院| 日韩欧美精品v在线| 毛片一级片免费看久久久久| 亚洲人成网站在线观看播放| 午夜a级毛片| 免费观看的影片在线观看| 99精品在免费线老司机午夜| www日本黄色视频网| 在线观看午夜福利视频| 国产精品人妻久久久久久| 国产真实伦视频高清在线观看| 好男人视频免费观看在线| 男女那种视频在线观看| 日本免费一区二区三区高清不卡| 两个人视频免费观看高清| 日韩欧美国产在线观看| 亚洲精品久久久久久婷婷小说 | 少妇裸体淫交视频免费看高清| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| 97在线视频观看| 亚洲精品456在线播放app| 国产精品1区2区在线观看.| 我的女老师完整版在线观看| 国产中年淑女户外野战色| 国产精品99久久久久久久久| 国产一区亚洲一区在线观看| 免费看av在线观看网站| 免费av不卡在线播放| 国产高清视频在线观看网站| 嘟嘟电影网在线观看| 人人妻人人澡人人爽人人夜夜 | 午夜福利高清视频| 免费大片18禁| 内地一区二区视频在线| 亚洲无线观看免费| 日韩av不卡免费在线播放| 白带黄色成豆腐渣| 精品欧美国产一区二区三| 亚洲在线自拍视频| 成人毛片60女人毛片免费| 亚洲真实伦在线观看| 少妇丰满av| 在线观看午夜福利视频| 久久久久久久亚洲中文字幕| 悠悠久久av| 国内揄拍国产精品人妻在线| a级毛片a级免费在线| 亚洲最大成人手机在线| 欧美极品一区二区三区四区| 99视频精品全部免费 在线| 久久精品国产自在天天线| 亚洲无线在线观看| 久久久精品欧美日韩精品| 可以在线观看的亚洲视频| 啦啦啦啦在线视频资源| 少妇的逼水好多| 在线a可以看的网站| 91aial.com中文字幕在线观看| 亚洲第一区二区三区不卡| 久久人人精品亚洲av| 五月玫瑰六月丁香| 久久九九热精品免费| 少妇人妻精品综合一区二区 | 日韩欧美国产在线观看| 亚洲欧美精品自产自拍| 白带黄色成豆腐渣| 国产精品免费一区二区三区在线| 卡戴珊不雅视频在线播放| 亚洲精品456在线播放app| 国产白丝娇喘喷水9色精品| av免费在线看不卡| 中国美女看黄片| 国产色婷婷99| 99久久久亚洲精品蜜臀av| 波多野结衣高清作品| 综合色丁香网| 亚洲精品乱码久久久v下载方式| www.av在线官网国产| 22中文网久久字幕| 99在线人妻在线中文字幕| 国产精品伦人一区二区| 日韩,欧美,国产一区二区三区 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| 舔av片在线| 婷婷色av中文字幕| 中文亚洲av片在线观看爽| 国产黄a三级三级三级人| 99热6这里只有精品| 麻豆一二三区av精品| 国产三级中文精品| 久久婷婷人人爽人人干人人爱| 欧美激情在线99| 丝袜喷水一区| 国产精品一区二区性色av| 成人国产麻豆网| 99国产精品一区二区蜜桃av| 99热这里只有是精品在线观看| 亚洲在久久综合| 婷婷亚洲欧美| 黄片wwwwww| 亚洲婷婷狠狠爱综合网| 天堂网av新在线| .国产精品久久| 国产高清视频在线观看网站| 最近的中文字幕免费完整| 亚洲成av人片在线播放无| 人人妻人人澡欧美一区二区| 免费观看人在逋| 国内久久婷婷六月综合欲色啪| 国产精品99久久久久久久久| 我的老师免费观看完整版| 久久久久久伊人网av| 亚洲av一区综合| 免费看日本二区| 亚洲精华国产精华液的使用体验 | 精品一区二区三区视频在线| 国产欧美日韩精品一区二区| 国产探花极品一区二区| 久久久久久久久久久丰满| 一个人看的www免费观看视频| 亚洲成人精品中文字幕电影| 嫩草影院新地址| 久久精品国产亚洲av天美| 国产黄片美女视频| 长腿黑丝高跟| 国产成人aa在线观看| 国产伦一二天堂av在线观看| 国产免费一级a男人的天堂| av在线天堂中文字幕| 欧美极品一区二区三区四区| 久久精品91蜜桃| 国产国拍精品亚洲av在线观看| 91久久精品国产一区二区成人| 校园春色视频在线观看| 国产极品精品免费视频能看的| 长腿黑丝高跟| 国产老妇伦熟女老妇高清| 成人鲁丝片一二三区免费| 国产亚洲精品av在线| 久久精品久久久久久噜噜老黄 | 国产一区二区三区在线臀色熟女| 蜜桃亚洲精品一区二区三区| 亚洲经典国产精华液单| 黑人高潮一二区| 最近的中文字幕免费完整| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| 狠狠狠狠99中文字幕| 久久精品国产亚洲av天美| 天堂av国产一区二区熟女人妻| 波野结衣二区三区在线| 亚洲av免费在线观看| 欧美xxxx性猛交bbbb| 免费无遮挡裸体视频| 国产黄a三级三级三级人| 国产中年淑女户外野战色| 久久久久久久久久成人| 亚洲乱码一区二区免费版| 午夜福利在线观看吧| 亚洲三级黄色毛片| 又粗又爽又猛毛片免费看| 色综合站精品国产| 国产单亲对白刺激| 白带黄色成豆腐渣| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 日本-黄色视频高清免费观看| 国产一区亚洲一区在线观看| 少妇熟女欧美另类| av在线老鸭窝| 成人亚洲欧美一区二区av| av在线天堂中文字幕| 国产日韩欧美在线精品| 久久久午夜欧美精品| 综合色丁香网| 最近最新中文字幕大全电影3| 小说图片视频综合网站| 国产精品一区二区三区四区久久| 舔av片在线| 美女 人体艺术 gogo| 亚洲在线观看片| 色综合色国产| 免费电影在线观看免费观看| 亚洲成人久久爱视频| 日韩国内少妇激情av| 亚洲国产精品sss在线观看| 嫩草影院入口| 国产精品电影一区二区三区| 97人妻精品一区二区三区麻豆| 成人午夜高清在线视频| 午夜精品国产一区二区电影 | 一边亲一边摸免费视频| 99riav亚洲国产免费| 97超碰精品成人国产| 日韩欧美一区二区三区在线观看| 亚洲国产色片| 免费av不卡在线播放| 最近视频中文字幕2019在线8| 成人综合一区亚洲| 日本与韩国留学比较| 亚洲无线在线观看| 又爽又黄无遮挡网站| 99久久九九国产精品国产免费| 久久久久久伊人网av| 久久午夜福利片| 少妇的逼水好多| 我的女老师完整版在线观看| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 中国美女看黄片| 亚洲国产精品合色在线| 一级毛片电影观看 | 少妇熟女欧美另类| 全区人妻精品视频| 久久人人爽人人爽人人片va| 直男gayav资源| 波多野结衣高清无吗| 免费观看的影片在线观看| 久久鲁丝午夜福利片| 色哟哟·www| 午夜福利成人在线免费观看| 1000部很黄的大片| 欧美性猛交╳xxx乱大交人| 老女人水多毛片| 神马国产精品三级电影在线观看| 国产精品av视频在线免费观看| 人妻少妇偷人精品九色| 国产免费男女视频| 国产成人精品婷婷| 欧美成人精品欧美一级黄| 亚洲欧洲日产国产| 一本一本综合久久| 亚洲欧美成人精品一区二区| 亚洲av男天堂| 国产麻豆成人av免费视频| 美女脱内裤让男人舔精品视频 | 狂野欧美白嫩少妇大欣赏| 国产探花极品一区二区| 少妇的逼好多水| 精品人妻一区二区三区麻豆| 久久人人精品亚洲av| 简卡轻食公司| 国产一级毛片七仙女欲春2| 女人十人毛片免费观看3o分钟| 麻豆一二三区av精品| 国产不卡一卡二| 好男人在线观看高清免费视频| 亚洲,欧美,日韩| 国产成人精品婷婷| 极品教师在线视频| 久久精品人妻少妇| 春色校园在线视频观看| 免费无遮挡裸体视频| 观看美女的网站| 麻豆一二三区av精品| 色吧在线观看| 国产亚洲欧美98| 女同久久另类99精品国产91| 一级二级三级毛片免费看| 国产成人a区在线观看| 久久国产乱子免费精品| 日本一二三区视频观看| 亚洲精华国产精华液的使用体验 | 观看免费一级毛片| 国产真实伦视频高清在线观看| 久99久视频精品免费| 人妻制服诱惑在线中文字幕| 26uuu在线亚洲综合色| 国产美女午夜福利| 国产毛片a区久久久久| 成人漫画全彩无遮挡| 亚洲丝袜综合中文字幕| 99久久中文字幕三级久久日本| 99热这里只有是精品在线观看| 日韩一本色道免费dvd| 国产精品久久久久久av不卡| 国产精品人妻久久久久久| 亚洲精品色激情综合| 亚洲成人久久爱视频| 插逼视频在线观看| 亚洲内射少妇av| 国产精品蜜桃在线观看 | 久久久成人免费电影| 最近的中文字幕免费完整| 国产免费一级a男人的天堂| 一本久久精品| 波野结衣二区三区在线| 中文字幕熟女人妻在线| 看黄色毛片网站| 男女做爰动态图高潮gif福利片| 天天躁日日操中文字幕| 少妇的逼水好多| 国产黄a三级三级三级人| 亚洲av电影不卡..在线观看| 亚洲精品日韩在线中文字幕 | 熟女电影av网| 青青草视频在线视频观看| 亚洲欧美日韩高清在线视频| 搡女人真爽免费视频火全软件| 99riav亚洲国产免费| 久久久久国产网址| 99九九线精品视频在线观看视频| 亚洲欧美成人精品一区二区| 男女边吃奶边做爰视频| 日韩大尺度精品在线看网址| 免费电影在线观看免费观看| 国产av一区在线观看免费| 波野结衣二区三区在线| 又粗又硬又长又爽又黄的视频 | 久久久欧美国产精品| 国产精品一区二区三区四区久久| 淫秽高清视频在线观看| 精品国产三级普通话版| 久久精品久久久久久久性| 久久久欧美国产精品| 26uuu在线亚洲综合色| 国产成人aa在线观看| 天堂中文最新版在线下载 | 国产成人一区二区在线| 男女做爰动态图高潮gif福利片| 99久国产av精品| 国产三级中文精品| 美女cb高潮喷水在线观看| 国产伦在线观看视频一区| 久久久久久大精品| 美女被艹到高潮喷水动态| 男女做爰动态图高潮gif福利片| 欧美成人免费av一区二区三区| 中文字幕av成人在线电影| 亚洲aⅴ乱码一区二区在线播放| 看片在线看免费视频| 热99在线观看视频| 国产精品久久久久久久久免| 深夜a级毛片| 非洲黑人性xxxx精品又粗又长| 国产精品综合久久久久久久免费| 成人无遮挡网站| 久久鲁丝午夜福利片| 我要搜黄色片| 91狼人影院| 国产精品精品国产色婷婷| 国产精品野战在线观看| 又粗又硬又长又爽又黄的视频 | 97超碰精品成人国产| 久久久久久久久久黄片| 亚洲人成网站高清观看| 我的女老师完整版在线观看| 国产在视频线在精品| 99久久中文字幕三级久久日本| 精品久久国产蜜桃| 欧美丝袜亚洲另类| 日韩成人伦理影院| 国产视频首页在线观看| 黄色欧美视频在线观看| 亚洲av成人av| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看 | av女优亚洲男人天堂| 欧美一区二区精品小视频在线| 亚洲最大成人手机在线| 一级黄色大片毛片| 亚洲国产精品成人综合色| 岛国毛片在线播放| 岛国在线免费视频观看| 亚洲,欧美,日韩| 一级二级三级毛片免费看| 久久久a久久爽久久v久久| 91麻豆精品激情在线观看国产| 午夜久久久久精精品| 少妇丰满av| 夜夜夜夜夜久久久久| 日产精品乱码卡一卡2卡三| 男人和女人高潮做爰伦理| 欧美不卡视频在线免费观看| av视频在线观看入口| 特级一级黄色大片| 丝袜喷水一区| 日韩精品有码人妻一区| 美女 人体艺术 gogo| 两个人的视频大全免费| 观看美女的网站| 国产精品不卡视频一区二区| 91久久精品国产一区二区三区| 日韩大尺度精品在线看网址| АⅤ资源中文在线天堂| 国产探花在线观看一区二区| 少妇的逼好多水| 午夜视频国产福利| 国产视频首页在线观看| 成人午夜精彩视频在线观看| 激情 狠狠 欧美| 中文资源天堂在线| 观看美女的网站| 亚洲av二区三区四区| 51国产日韩欧美| 国产不卡一卡二| 啦啦啦韩国在线观看视频| 国产黄片视频在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 色视频www国产| 久久久久久伊人网av| 91精品一卡2卡3卡4卡| 欧美成人一区二区免费高清观看| 国产av在哪里看| 国产又黄又爽又无遮挡在线| 国产精品三级大全| 岛国毛片在线播放| 最近2019中文字幕mv第一页| 亚洲成人中文字幕在线播放| 成人一区二区视频在线观看| 一个人观看的视频www高清免费观看| 成年女人看的毛片在线观看| 国产片特级美女逼逼视频| 国内少妇人妻偷人精品xxx网站| 热99re8久久精品国产| 亚洲欧美日韩无卡精品| 国产v大片淫在线免费观看| 草草在线视频免费看| 女的被弄到高潮叫床怎么办| av国产免费在线观看| 中出人妻视频一区二区| 村上凉子中文字幕在线| 狂野欧美白嫩少妇大欣赏| 久久久久久九九精品二区国产| 午夜福利视频1000在线观看| 日本-黄色视频高清免费观看| 麻豆国产av国片精品| 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| 麻豆久久精品国产亚洲av| 国产真实乱freesex| 亚洲精品自拍成人| 欧美3d第一页| 精品国产三级普通话版| 99热这里只有是精品50| 久久久精品94久久精品| 国产亚洲精品av在线| 国产黄片视频在线免费观看| 看十八女毛片水多多多| 国产淫片久久久久久久久| av免费观看日本| a级一级毛片免费在线观看| 一区二区三区免费毛片| 蜜臀久久99精品久久宅男| 最近最新中文字幕大全电影3| 搡老妇女老女人老熟妇| 日本免费一区二区三区高清不卡| 内地一区二区视频在线| 国产精品精品国产色婷婷| 永久网站在线| 国产精品伦人一区二区| 亚洲欧洲日产国产| 99视频精品全部免费 在线| 国产亚洲欧美98| 亚州av有码| 国产伦在线观看视频一区| 国产精品嫩草影院av在线观看| 欧美+亚洲+日韩+国产| 亚洲av不卡在线观看| 亚洲18禁久久av| 亚洲经典国产精华液单| 桃色一区二区三区在线观看| 天堂√8在线中文| 国产精品一区二区在线观看99 | 日本色播在线视频| 国产伦精品一区二区三区四那| 18禁在线无遮挡免费观看视频| 国产午夜精品久久久久久一区二区三区| 国产探花在线观看一区二区| 国模一区二区三区四区视频| 97人妻精品一区二区三区麻豆| 成人一区二区视频在线观看| 午夜老司机福利剧场| 国产 一区精品| 国产色婷婷99| 好男人视频免费观看在线| 日韩欧美精品v在线| 午夜视频国产福利| 亚洲四区av| 久久99热这里只有精品18| 不卡视频在线观看欧美| 精品无人区乱码1区二区|