• <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麻豆av| 久久人妻熟女aⅴ| 人妻夜夜爽99麻豆av| 搡女人真爽免费视频火全软件| 日日撸夜夜添| 免费观看无遮挡的男女| 亚洲美女搞黄在线观看| 国产在线视频一区二区| 国产黄频视频在线观看| 日韩成人av中文字幕在线观看| 欧美xxxx性猛交bbbb| 午夜影院在线不卡| 夫妻性生交免费视频一级片| 我的女老师完整版在线观看| 国产亚洲精品久久久com| 如何舔出高潮| 人妻制服诱惑在线中文字幕| 精品久久久久久久久亚洲| 在线观看国产h片| 国产免费一区二区三区四区乱码| 日日啪夜夜爽| 岛国毛片在线播放| 少妇被粗大猛烈的视频| 综合色丁香网| 晚上一个人看的免费电影| 国产毛片在线视频| 亚洲天堂av无毛| 99九九线精品视频在线观看视频| 妹子高潮喷水视频| 久久鲁丝午夜福利片| 天天躁夜夜躁狠狠久久av| 久久人人爽人人爽人人片va| 亚洲精品日本国产第一区| 久久久国产欧美日韩av| 制服丝袜香蕉在线| 国产成人精品婷婷| 97精品久久久久久久久久精品| 免费人成在线观看视频色| 大话2 男鬼变身卡| 多毛熟女@视频| 亚洲精品日本国产第一区| 国产成人精品一,二区| 大陆偷拍与自拍| 又大又黄又爽视频免费| 午夜老司机福利剧场| 国产深夜福利视频在线观看| 插逼视频在线观看| 欧美精品国产亚洲| 丝袜脚勾引网站| 夜夜爽夜夜爽视频| a级毛片在线看网站| 国产国拍精品亚洲av在线观看| 大香蕉97超碰在线| 午夜影院在线不卡| 美女大奶头黄色视频| 大码成人一级视频| 51国产日韩欧美| 久热这里只有精品99| 人妻夜夜爽99麻豆av| 九色成人免费人妻av| 国产无遮挡羞羞视频在线观看| a级毛片黄视频| 久久国产精品男人的天堂亚洲 | 久久国产精品大桥未久av| 国产国语露脸激情在线看| 高清不卡的av网站| av在线老鸭窝| 亚洲精品乱码久久久久久按摩| 高清av免费在线| 国产片特级美女逼逼视频| 大又大粗又爽又黄少妇毛片口| 亚洲av日韩在线播放| 精品久久蜜臀av无| 免费黄色在线免费观看| 亚洲欧美成人精品一区二区| 这个男人来自地球电影免费观看 | 日韩,欧美,国产一区二区三区| 午夜av观看不卡| 欧美日韩视频高清一区二区三区二| 久久精品久久久久久久性| 亚洲性久久影院| 天堂中文最新版在线下载| 亚洲精品aⅴ在线观看| 91精品一卡2卡3卡4卡| 赤兔流量卡办理| 久久婷婷青草| 乱人伦中国视频| 丁香六月天网| 97超碰精品成人国产| 91精品三级在线观看| 国产成人精品在线电影| 亚洲av日韩在线播放| 欧美日韩在线观看h| 国产国语露脸激情在线看| 久久久精品94久久精品| 99九九在线精品视频| 久久久久国产精品人妻一区二区| 亚洲,欧美,日韩| 亚洲美女黄色视频免费看| 日本爱情动作片www.在线观看| 亚洲精品成人av观看孕妇| 午夜免费观看性视频| 大话2 男鬼变身卡| 多毛熟女@视频| 国产乱来视频区| 麻豆成人av视频| 日韩中字成人| 久久国内精品自在自线图片| 中文字幕久久专区| 这个男人来自地球电影免费观看 | 国产免费现黄频在线看| 少妇的逼好多水| 99热全是精品| .国产精品久久| 18禁裸乳无遮挡动漫免费视频| 夜夜骑夜夜射夜夜干| 精品一区二区三卡| 制服诱惑二区| 18禁在线无遮挡免费观看视频| 涩涩av久久男人的天堂| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 久久精品国产亚洲网站| 免费av不卡在线播放| 国产精品国产三级国产av玫瑰| 精品一区二区三区视频在线| 日韩欧美一区视频在线观看| 黄片无遮挡物在线观看| 插阴视频在线观看视频| 国产不卡av网站在线观看| 街头女战士在线观看网站| 亚洲色图综合在线观看| 一个人免费看片子| 大香蕉久久成人网| av国产久精品久网站免费入址| 丁香六月天网| 婷婷色麻豆天堂久久| 亚洲欧美日韩卡通动漫| 免费少妇av软件| 日韩一本色道免费dvd| 成人18禁高潮啪啪吃奶动态图 | 久久久久国产网址| 桃花免费在线播放| 久久精品夜色国产| 成人亚洲精品一区在线观看| 亚洲伊人久久精品综合| 水蜜桃什么品种好| 黄色欧美视频在线观看| 91aial.com中文字幕在线观看| 91精品国产国语对白视频| 国产亚洲av片在线观看秒播厂| 亚洲第一区二区三区不卡| √禁漫天堂资源中文www| 一本久久精品| a级毛片黄视频| 婷婷色麻豆天堂久久| 精品国产乱码久久久久久小说| 午夜激情久久久久久久| 亚洲欧美日韩卡通动漫| 亚洲欧美清纯卡通| 久久精品夜色国产| 亚洲精品456在线播放app| 内地一区二区视频在线| 国产精品国产av在线观看| 亚洲国产日韩一区二区| 欧美精品一区二区大全| 国产 一区精品| 亚洲精品久久久久久婷婷小说| 各种免费的搞黄视频| 色哟哟·www| 日本爱情动作片www.在线观看| av女优亚洲男人天堂| 国产69精品久久久久777片| 中文字幕制服av| 国产不卡av网站在线观看| 少妇高潮的动态图| 日产精品乱码卡一卡2卡三| 一本色道久久久久久精品综合| 精品少妇久久久久久888优播| av国产精品久久久久影院| 亚洲欧美一区二区三区国产| 国产日韩欧美亚洲二区| 不卡视频在线观看欧美| 午夜日本视频在线| 欧美日韩在线观看h| 一本大道久久a久久精品| 久久精品久久精品一区二区三区| 国产精品免费大片| 欧美性感艳星| 热99国产精品久久久久久7| 乱码一卡2卡4卡精品| 精品一品国产午夜福利视频| 天天操日日干夜夜撸| 视频在线观看一区二区三区| 免费观看在线日韩| 欧美成人午夜免费资源| 免费人成在线观看视频色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| 99国产综合亚洲精品| 亚洲精品aⅴ在线观看| 国产免费一区二区三区四区乱码| 高清在线视频一区二区三区| 一区在线观看完整版| 国产精品人妻久久久影院| 欧美日韩国产mv在线观看视频| 精品99又大又爽又粗少妇毛片| 国产探花极品一区二区| 成人亚洲欧美一区二区av| 免费观看av网站的网址| 美女大奶头黄色视频| 欧美日韩视频精品一区| 又大又黄又爽视频免费| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品色激情综合| 丰满乱子伦码专区| 天堂8中文在线网| 国产成人av激情在线播放 | 高清欧美精品videossex| 国产 精品1| 麻豆精品久久久久久蜜桃| 免费看av在线观看网站| 高清在线视频一区二区三区| 精品一品国产午夜福利视频| 亚洲久久久国产精品| 女人久久www免费人成看片| 国产成人精品在线电影| 成人毛片a级毛片在线播放| 一区二区三区精品91| 亚洲欧美成人综合另类久久久| 大片电影免费在线观看免费| 9色porny在线观看| 热re99久久国产66热| 色婷婷久久久亚洲欧美| 亚洲成人一二三区av| 久久久国产欧美日韩av| 国产成人午夜福利电影在线观看| 亚洲国产日韩一区二区| 亚洲欧美精品自产自拍| 一区二区三区精品91| 亚洲综合色惰| 国产高清国产精品国产三级| 99久久综合免费| 九色亚洲精品在线播放| 狂野欧美白嫩少妇大欣赏| 日本欧美国产在线视频| 亚洲成人手机| 精品卡一卡二卡四卡免费| 欧美xxxx性猛交bbbb| 亚洲不卡免费看| 少妇的逼水好多| 日韩熟女老妇一区二区性免费视频| 中文字幕人妻丝袜制服| 欧美亚洲 丝袜 人妻 在线| 日本vs欧美在线观看视频| 久久ye,这里只有精品| 99九九在线精品视频| 韩国av在线不卡| 老女人水多毛片| 精品少妇黑人巨大在线播放| 日本免费在线观看一区| 2022亚洲国产成人精品| 亚洲美女视频黄频| av又黄又爽大尺度在线免费看| 在线观看免费高清a一片| 如何舔出高潮| 久久狼人影院| 国产成人精品福利久久| 国产精品 国内视频| 国产一区二区在线观看日韩| 91在线精品国自产拍蜜月| 中文字幕精品免费在线观看视频 | 丰满乱子伦码专区| 国产精品国产三级国产专区5o| 男女无遮挡免费网站观看| 午夜福利网站1000一区二区三区| 精品久久久久久电影网| 欧美日韩国产mv在线观看视频| 国产日韩欧美在线精品| 只有这里有精品99| 18禁在线无遮挡免费观看视频| 美女中出高潮动态图| 蜜臀久久99精品久久宅男| 国产精品久久久久久av不卡| 中文字幕久久专区| 国产一区二区三区av在线| 美女内射精品一级片tv| 成年女人在线观看亚洲视频| 日本黄大片高清| 久久久精品94久久精品| 香蕉精品网在线| 免费看光身美女| 一级片'在线观看视频| 精品99又大又爽又粗少妇毛片| 国产白丝娇喘喷水9色精品| www.色视频.com| 久久久久精品性色| 成年人免费黄色播放视频| 我的女老师完整版在线观看| 欧美日韩av久久| 99re6热这里在线精品视频| 七月丁香在线播放| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 在线天堂最新版资源| 又粗又硬又长又爽又黄的视频| 女性被躁到高潮视频| 午夜免费男女啪啪视频观看| 亚洲精品视频女| 青青草视频在线视频观看| 韩国高清视频一区二区三区| 性高湖久久久久久久久免费观看| 亚州av有码| 三级国产精品片| 人人妻人人爽人人添夜夜欢视频| 日韩av在线免费看完整版不卡| 中文字幕av电影在线播放| 丝袜喷水一区| 久久久国产精品麻豆| 日产精品乱码卡一卡2卡三| 亚洲欧洲日产国产| 狂野欧美激情性bbbbbb| 美女视频免费永久观看网站| 欧美激情 高清一区二区三区| 91久久精品电影网| 国精品久久久久久国模美| 丝袜在线中文字幕| 99国产综合亚洲精品| av又黄又爽大尺度在线免费看| 中文字幕精品免费在线观看视频 | 久久毛片免费看一区二区三区| 五月伊人婷婷丁香| 午夜福利影视在线免费观看| 日日摸夜夜添夜夜添av毛片| 一边亲一边摸免费视频| 大又大粗又爽又黄少妇毛片口| 国产精品国产三级国产av玫瑰| 亚洲少妇的诱惑av| 一级毛片 在线播放| 日日啪夜夜爽| 一本色道久久久久久精品综合| av黄色大香蕉| 成人国产av品久久久| 国产一区二区三区av在线| 韩国av在线不卡| 亚洲av中文av极速乱| 亚洲人与动物交配视频| 高清黄色对白视频在线免费看| 免费观看无遮挡的男女| av在线app专区| 只有这里有精品99| 亚洲四区av| 欧美激情极品国产一区二区三区 | 妹子高潮喷水视频| 视频在线观看一区二区三区| 人妻系列 视频| 欧美亚洲 丝袜 人妻 在线| 亚洲精品成人av观看孕妇| 视频在线观看一区二区三区| 丰满饥渴人妻一区二区三| 国产黄色免费在线视频| 久久人人爽人人片av| 最新的欧美精品一区二区| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 国产精品一二三区在线看| 只有这里有精品99| 91精品国产九色| 美女主播在线视频| 精品99又大又爽又粗少妇毛片| 少妇被粗大猛烈的视频| 久久久久久久久久久丰满| 亚洲国产精品专区欧美| 一级爰片在线观看| 国产亚洲午夜精品一区二区久久| 久久久精品免费免费高清| 亚洲内射少妇av| 国产淫语在线视频| 欧美国产精品一级二级三级| 日日摸夜夜添夜夜添av毛片| 国产黄片视频在线免费观看| 亚洲经典国产精华液单| 久久国产亚洲av麻豆专区| 国产不卡av网站在线观看| 精品久久蜜臀av无| 九九爱精品视频在线观看| 2022亚洲国产成人精品| 99久国产av精品国产电影| 精品熟女少妇av免费看| 国产精品一区www在线观看| 22中文网久久字幕| 国产亚洲一区二区精品| 久久精品国产自在天天线| 蜜桃久久精品国产亚洲av| 亚洲国产av影院在线观看| 国产在线一区二区三区精| 一级毛片 在线播放| 一级爰片在线观看| 欧美少妇被猛烈插入视频| kizo精华| 人人妻人人添人人爽欧美一区卜| 国产白丝娇喘喷水9色精品| 久久午夜福利片| 国产男女超爽视频在线观看| 国产精品久久久久成人av| 大陆偷拍与自拍| 亚洲精品国产色婷婷电影| 老司机影院成人| 9色porny在线观看| 久久久久久久亚洲中文字幕| 美女主播在线视频| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 亚洲av男天堂| 亚洲av福利一区| 自线自在国产av| 国产精品国产三级专区第一集| 免费日韩欧美在线观看| 人人妻人人澡人人爽人人夜夜| 18禁在线播放成人免费| 日韩一区二区视频免费看| 黑人高潮一二区| 九九爱精品视频在线观看| 欧美日韩在线观看h| 丝袜脚勾引网站| 免费观看在线日韩| 精品国产国语对白av| 十八禁网站网址无遮挡| 国产精品欧美亚洲77777| 91成人精品电影| 搡女人真爽免费视频火全软件| 亚洲精品中文字幕在线视频| 亚洲精品久久久久久婷婷小说| 在线 av 中文字幕| 国产av码专区亚洲av| 美女cb高潮喷水在线观看| 91久久精品电影网| 妹子高潮喷水视频| 精品一区二区免费观看| 日韩精品免费视频一区二区三区 | 男女边吃奶边做爰视频| 亚洲在久久综合| 亚洲情色 制服丝袜| 午夜免费观看性视频| 国产精品99久久久久久久久| 日韩av在线免费看完整版不卡| 最后的刺客免费高清国语| 欧美激情 高清一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 男女边摸边吃奶| 国产视频首页在线观看| 亚洲av在线观看美女高潮| 在线观看免费日韩欧美大片 | 国产精品成人在线| 少妇 在线观看| 你懂的网址亚洲精品在线观看| 国产亚洲欧美精品永久| 乱人伦中国视频| 一区在线观看完整版| 亚洲精品第二区| 国产熟女欧美一区二区| 秋霞伦理黄片| 最黄视频免费看| 在线观看www视频免费| 久久久亚洲精品成人影院| 中文字幕人妻丝袜制服| 亚洲人成网站在线观看播放| 日韩大片免费观看网站| 三上悠亚av全集在线观看| 黄片无遮挡物在线观看| 久久人人爽人人片av| 久久久亚洲精品成人影院| 久久久久网色| 99视频精品全部免费 在线| 九色成人免费人妻av| 久久国产精品男人的天堂亚洲 | 成人免费观看视频高清| 久久久久精品久久久久真实原创| 少妇 在线观看| 久久久久久久久久久免费av| 街头女战士在线观看网站| 久久久久精品性色| 亚洲av国产av综合av卡| 久久人人爽人人片av| freevideosex欧美| 国产日韩欧美视频二区| av在线观看视频网站免费| 女性被躁到高潮视频| 制服丝袜香蕉在线| 99久国产av精品国产电影| 日本-黄色视频高清免费观看| 一本大道久久a久久精品| 夜夜看夜夜爽夜夜摸| av不卡在线播放| av网站免费在线观看视频| 99热全是精品| 欧美日韩成人在线一区二区| 日本黄色片子视频| 91国产中文字幕| 欧美3d第一页| 成人手机av| a级毛色黄片| 视频在线观看一区二区三区| 美女国产高潮福利片在线看| 天天操日日干夜夜撸| 最近手机中文字幕大全| 国产成人精品在线电影| 日本爱情动作片www.在线观看| 尾随美女入室| 中国国产av一级| 在线观看免费视频网站a站| 亚洲一级一片aⅴ在线观看| 久久久久久人妻| 一区二区三区精品91| 九色亚洲精品在线播放| 欧美xxⅹ黑人| 国模一区二区三区四区视频| 日韩,欧美,国产一区二区三区| 青春草亚洲视频在线观看| 大又大粗又爽又黄少妇毛片口| 女人精品久久久久毛片| 99热这里只有精品一区| 欧美日韩在线观看h| 国产 精品1| 久久免费观看电影| 亚洲国产最新在线播放| 天天影视国产精品| 人人妻人人澡人人看| 亚洲欧美中文字幕日韩二区| www.av在线官网国产| 国国产精品蜜臀av免费| 在线观看免费视频网站a站| av又黄又爽大尺度在线免费看| 国产高清三级在线| 久久人人爽人人爽人人片va| 在线播放无遮挡| 天美传媒精品一区二区| 国产在线视频一区二区| av国产精品久久久久影院| 精品酒店卫生间| 美女福利国产在线| 亚洲在久久综合| 国产高清有码在线观看视频| av有码第一页| 久久韩国三级中文字幕| 国产免费现黄频在线看| 在线观看免费高清a一片| 女人久久www免费人成看片| 好男人视频免费观看在线| 亚洲成色77777| 人妻人人澡人人爽人人| 日本免费在线观看一区| 中文乱码字字幕精品一区二区三区| 国产永久视频网站| 最近的中文字幕免费完整| 成人手机av| 成年av动漫网址| 色吧在线观看| 国产淫语在线视频| 狂野欧美激情性xxxx在线观看| 久久久久久久久久成人| 国产成人免费无遮挡视频| 99精国产麻豆久久婷婷| 亚洲精品日本国产第一区| a级毛片在线看网站| 另类精品久久| 日韩不卡一区二区三区视频在线| 国产亚洲av片在线观看秒播厂| 国产国拍精品亚洲av在线观看| 免费看不卡的av| 久久久久久久精品精品| 欧美最新免费一区二区三区| 国内精品宾馆在线| 日韩av不卡免费在线播放| 亚洲五月色婷婷综合| 日本色播在线视频| 精品人妻一区二区三区麻豆| 青青草视频在线视频观看| 三级国产精品欧美在线观看| 亚洲精品国产av蜜桃| 伦理电影免费视频| 狠狠婷婷综合久久久久久88av| 欧美成人精品欧美一级黄| 免费av中文字幕在线| 免费高清在线观看日韩| 欧美3d第一页| 日本猛色少妇xxxxx猛交久久| videosex国产| 建设人人有责人人尽责人人享有的| 大码成人一级视频| 久久久久久久久久久免费av| 国产视频内射| 一边亲一边摸免费视频| 成人亚洲欧美一区二区av| 蜜桃国产av成人99| 免费av中文字幕在线| 你懂的网址亚洲精品在线观看| 欧美成人精品欧美一级黄| 国产片内射在线| 永久网站在线| 亚洲图色成人| 中国国产av一级| 国产免费一级a男人的天堂| 99久久人妻综合| 人人妻人人添人人爽欧美一区卜| 美女大奶头黄色视频| 乱码一卡2卡4卡精品| 高清不卡的av网站| 国产精品国产三级专区第一集| 日韩视频在线欧美| 99九九在线精品视频| 日本午夜av视频| 国产不卡av网站在线观看| 老司机影院成人|