搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

摩擦不对称耦合布朗马达的定向输运

张旭 曹佳慧 艾保全 高天附 郑志刚

引用本文:
Citation:

摩擦不对称耦合布朗马达的定向输运

张旭, 曹佳慧, 艾保全, 高天附, 郑志刚

Investigation on the directional transportation of coupled Brownian motors with asymmetric friction

Zhang Xu, Cao Jia-Hui, Ai Bao-Quan, Gao Tian-Fu, Zheng Zhi-Gang
PDF
HTML
导出引用
  • 本文研究了耦合布朗马达在溶液摩擦对称性破缺条件下的定向输运问题. 详细讨论了溶液阻尼系数比, 周期外力以及耦合作用对布朗粒子定向输运的影响. 研究发现: 溶液摩擦的不对称性不仅能促进摩擦棘轮的定向输运, 且合适的溶液阻尼系数比还可以使耦合棘轮的定向输运达到最强. 此外, 一定摩擦阻尼条件下外力振幅可以使耦合粒子的定向输运呈现多峰结构. 研究还发现合适的自由长度和耦合强度也能增强摩擦棘轮的输运特性. 研究结果不仅能够启发实验上选取合适的摩擦阻尼来优化布朗马达的输运, 也可为信息容量大、反应速度快的纳米器件的制备提供一定的理论参考.
    Biological molecular motors exist in cells widely. They can make use of intracellular free energy to complete all kinds of internal biological transports by transforming chemical energy into mechanical energy. The kind of directional movement of biological molecular motors plays a very important role in intracellular material transportation. In order to study the transport mechanism of molecular motors further, a large number of ratchet models are proposed, such as rocking ratchets and flashing ratchets. By investigating various kinds of ratchets we can not only understand the directional movement mechanism of Brownian particles, but can find suitable conditions in which the performance of Brownian motors’ directional transportation could be enhanced. Meanwhile, the investigation of ratchets could also be applied in manufacturing nanometer devices.At present, the directional transportation of Brownian ratchet has attracted extensive interests of researchers. In general, most friction factors of Brownian ratchet models are considered unit. In fact, the concentration of solutions and cell fluid impurity affect the actual frictional damping conditions, so the real frictional coefficient of Brownian motors is often changed. In addition, lots of experimental studies have shown that the movement of Brownian motors is collectively directed motion and the kind of directional movement is induced by intermolecular coupling interaction. As a result, it is more valuable to investigate the transporting performance of coupled Brownian particles that existed in different frictional damping conditions. In order to enhance the transporting performance of Brownian ratchet in different frictional damping conditions, we discuss how the frictional damping factor influences the directional movement of coupled Brownian particles deeply when Brownian particles drag loads. In this paper, we established the overdamped frictional ratchets, and then we investigated how frictional damping coefficient ratio, coupling strength and external force amplitude affect the transportation of coupled Brownian ratchets. On the basis of the investigation, some interesting results are found. The directional transport of frictional ratchets can be promoted by adjusting the frictional damping factor. Besides, the transportation can obtain the maximum under the appropriate friction factor case. In addition, under certain frictional damping condition, the directional transportation of the friction ratchets present multi-peak structure as the external force amplitude increases. Meanwhile, the appropriate free length and coupling strength can also enhance the transportation characteristic of frictional ratchets. All conclusions obtained in this paper can not be applied in selecting suitable frictional damping conditions experimentally to improve the directional transportation of coupled Brownian ratchets, but they can also be used in developing and manufacturing nanometer devices.
      通信作者: 高天附, tianfugao@synu.edu.cn ; 郑志刚, zgzheng@hqu.edu.cn
    • 基金项目: 国家级-国家自然科学基金(11875135)
      Corresponding author: Gao Tian-Fu, tianfugao@synu.edu.cn ; Zheng Zhi-Gang, zgzheng@hqu.edu.cn
    [1]

    Xie P 2010 Int. J. Biol. Sci. 6 665

    [2]

    舒咬根, 欧阳钟灿 2007 物理 36 735Google Scholar

    Shu Y G, Ouyang Z C 2007 Physics 36 735Google Scholar

    [3]

    Matsuda A, Li J, Brumm P, Adachi T, Inoue Y, Kim T 2019 Biolphys. J. 116 2161Google Scholar

    [4]

    Brown A I, Sivak D A 2019 Europhys. Lett. 126 40004Google Scholar

    [5]

    Grossert C, Leder M, Denisov S, Hänggi P, Weitz M 2016 Nat. Commun. 7 10440Google Scholar

    [6]

    Lee A K, Lee J R, Lee K H 2012 Physica B 407 4298Google Scholar

    [7]

    Wickenbrock A, Cubero D, Wahab N A A, Phoonthong P, Renzoni F 2011 Phys. Rev. E 84 021127Google Scholar

    [8]

    Lade S J 2010 J. Phys. A: Math. Theor. 43 249801Google Scholar

    [9]

    Ethier S N, Lee J L 2018 R. Soc. OpenSci. 5 171685Google Scholar

    [10]

    Vorotnikov D 2013 J. Math. Biol. 68 1677

    [11]

    Roca F, Villaluenga J P G, Dinis L 2014 Europhys. Lett. 107 10006Google Scholar

    [12]

    Bao J D 2018 College Physics 37 11

    [13]

    Huang S N, Zhu W J, Huang X Q, Ai B Q, Li F G 2019 Chin. Phys. B 28 040502Google Scholar

    [14]

    艾保全, 王先菊, 刘国涛, 刘良钢, 郑华 2001 中山大学学报: 自然科学版 40 27

    Ai B Q, Wang X J, Liu G T, Liu L G, Zheng H 2001 ActaSci. Natur. Univ. Sunyatseni 40 27

    [15]

    吕明涛, 延明月, 艾保全, 高天附, 郑志刚 2017 物理学报 66 220501Google Scholar

    Lv M T, Yan M Y, Ai B Q, Gao T F, Zheng Z G 2017 Acta. Phys. Sin. 66 220501Google Scholar

    [16]

    Holmes H R, Böhringer K F 2018 Adv. Colloid. Interface. Sci. 255 18Google Scholar

    [17]

    Latorre J C, Kramer P R, Pavliotis G A 2014 J. Comput. Phys. 257 57Google Scholar

    [18]

    Ryabov A, Holubec V, Yaghoubi M H, Varga M, Foulaadvand M E, Chvosta P 2016 J. Stat. Mech.: Theory Exp. 2016 093202Google Scholar

    [19]

    Radtke M, Patzelt A, Knorr F, Lademann J, Netz R R 2017 Eur. J. Pharm. Biopharm. 116 125Google Scholar

    [20]

    Hojo M, Arai N, Ebisuzaki T 2018 Mol. Simul. 44 523Google Scholar

    [21]

    Pattanayak S, Das R, Kumar M, Mishra S 2019 Eur. Phys. J. E 42 62Google Scholar

    [22]

    延明月, 张旭, 刘晨昊, 黄仁忠, 高天附, 郑志刚 2018 物理学报 67 190501Google Scholar

    Yan M Y, Zhang X, Liu C H, Huang R Z, Gao T F, Zheng Z G 2018 Acta Phys. Sin. 67 190501Google Scholar

    [23]

    Lombardo J, Broadwater D, Collins R, Cebe K, Brady R, Harrison S 2019 Hum. Pathol. 86 129Google Scholar

    [24]

    Toyoshima H 1951 J. Japn. Fores. Soc. 33 203

    [25]

    Toyoshima H 1951 J. Japn. Fores. Soc. 33 83

    [26]

    Pascu M N, Popescu I 2016 Stoch. Proc. Appl. 126 628Google Scholar

    [27]

    Jacka S D, Mijatović A, Širaj D 2014 Stoch. Proc. Appl. 124 1055Google Scholar

    [28]

    Atzberger P J 2006 Phys. Lett. A 351 225Google Scholar

    [29]

    Guo W, Du L C, Liu Z Z, Yang H, Mei D C 2017 Chin. Phys. B 26 010502Google Scholar

    [30]

    von Gehlen S, Evstigneev M, Reimann P 2009 Phys. Rev. E 79 031114Google Scholar

    [31]

    Saikia S, Mahato M C 2010 Physica A 389 4052Google Scholar

    [32]

    von Gehlen S, Evstigneev M, Reimann P 2008 Phys. Rev. E 77 031106Google Scholar

  • 图 1  摩擦对称性破缺耦合粒子的模型示意图

    Fig. 1.  Coupling particles model with broken friction symmetry.

    图 2  耦合粒子的平均速度$\left\langle v \right\rangle$随弹簧自由长度l的变化曲线, 其中γ1 = 5.0, k = 20.0, α = 0.3, A = 1.0, λ = 0.1

    Fig. 2.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the free length l, where ${\gamma _1} = 5.0$, $k = 20.0$, $\alpha = 0.3$, $A = 1.0$, $\lambda = 0.1$.

    图 3  耦合粒子的平均速度$\left\langle v \right\rangle$随耦合强度k的变化曲线, 插图为$A = 20.0$时的变化曲线, 其中${\gamma _1} = 5.0$, $l = 0.25$, $\alpha = 0.3$, $D = 0.35$, $\lambda = 0.1$

    Fig. 3.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the strength of coupling k, the insert figure described the average velocity curve when $A = 20.0$, where ${\gamma _1} = 5.0$, $l = 0.25$, $\alpha = 0.3$, $D = 0.35$, $\lambda = 0.1$.

    图 4  耦合粒子的平均速度$\left\langle v \right\rangle$随外力振幅A的变化曲线, 其中γ1 = 5.0, l = 0.25, k = 20.0, D = 0.35, λ = 0.1

    Fig. 4.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the amplitude A of external force, where ${\gamma _1} = 5.0$, $l = 0.25$, $k = 20.0$, $D = 0.35$, $\lambda = 0.1$.

    图 5  耦合粒子的平均速度$\left\langle v \right\rangle$随摩擦阻尼系数比$\alpha$的变化曲线, 插图为$\lambda = 0.6$时的变化曲线, 其中${\gamma _1} = 5.0$, $k = 20.0$, $l = 0.25$, $A = 20.0$, $D = 0.35$

    Fig. 5.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the frictional damping coefficient ratio $\alpha$, the insert figure described the average velocity curve when $\lambda = 0.6$, where ${\gamma _1} = 5.0$, $k = 20.0$, $l = 0.25$, $A = 20.0$, $D = 0.35$.

  • [1]

    Xie P 2010 Int. J. Biol. Sci. 6 665

    [2]

    舒咬根, 欧阳钟灿 2007 物理 36 735Google Scholar

    Shu Y G, Ouyang Z C 2007 Physics 36 735Google Scholar

    [3]

    Matsuda A, Li J, Brumm P, Adachi T, Inoue Y, Kim T 2019 Biolphys. J. 116 2161Google Scholar

    [4]

    Brown A I, Sivak D A 2019 Europhys. Lett. 126 40004Google Scholar

    [5]

    Grossert C, Leder M, Denisov S, Hänggi P, Weitz M 2016 Nat. Commun. 7 10440Google Scholar

    [6]

    Lee A K, Lee J R, Lee K H 2012 Physica B 407 4298Google Scholar

    [7]

    Wickenbrock A, Cubero D, Wahab N A A, Phoonthong P, Renzoni F 2011 Phys. Rev. E 84 021127Google Scholar

    [8]

    Lade S J 2010 J. Phys. A: Math. Theor. 43 249801Google Scholar

    [9]

    Ethier S N, Lee J L 2018 R. Soc. OpenSci. 5 171685Google Scholar

    [10]

    Vorotnikov D 2013 J. Math. Biol. 68 1677

    [11]

    Roca F, Villaluenga J P G, Dinis L 2014 Europhys. Lett. 107 10006Google Scholar

    [12]

    Bao J D 2018 College Physics 37 11

    [13]

    Huang S N, Zhu W J, Huang X Q, Ai B Q, Li F G 2019 Chin. Phys. B 28 040502Google Scholar

    [14]

    艾保全, 王先菊, 刘国涛, 刘良钢, 郑华 2001 中山大学学报: 自然科学版 40 27

    Ai B Q, Wang X J, Liu G T, Liu L G, Zheng H 2001 ActaSci. Natur. Univ. Sunyatseni 40 27

    [15]

    吕明涛, 延明月, 艾保全, 高天附, 郑志刚 2017 物理学报 66 220501Google Scholar

    Lv M T, Yan M Y, Ai B Q, Gao T F, Zheng Z G 2017 Acta. Phys. Sin. 66 220501Google Scholar

    [16]

    Holmes H R, Böhringer K F 2018 Adv. Colloid. Interface. Sci. 255 18Google Scholar

    [17]

    Latorre J C, Kramer P R, Pavliotis G A 2014 J. Comput. Phys. 257 57Google Scholar

    [18]

    Ryabov A, Holubec V, Yaghoubi M H, Varga M, Foulaadvand M E, Chvosta P 2016 J. Stat. Mech.: Theory Exp. 2016 093202Google Scholar

    [19]

    Radtke M, Patzelt A, Knorr F, Lademann J, Netz R R 2017 Eur. J. Pharm. Biopharm. 116 125Google Scholar

    [20]

    Hojo M, Arai N, Ebisuzaki T 2018 Mol. Simul. 44 523Google Scholar

    [21]

    Pattanayak S, Das R, Kumar M, Mishra S 2019 Eur. Phys. J. E 42 62Google Scholar

    [22]

    延明月, 张旭, 刘晨昊, 黄仁忠, 高天附, 郑志刚 2018 物理学报 67 190501Google Scholar

    Yan M Y, Zhang X, Liu C H, Huang R Z, Gao T F, Zheng Z G 2018 Acta Phys. Sin. 67 190501Google Scholar

    [23]

    Lombardo J, Broadwater D, Collins R, Cebe K, Brady R, Harrison S 2019 Hum. Pathol. 86 129Google Scholar

    [24]

    Toyoshima H 1951 J. Japn. Fores. Soc. 33 203

    [25]

    Toyoshima H 1951 J. Japn. Fores. Soc. 33 83

    [26]

    Pascu M N, Popescu I 2016 Stoch. Proc. Appl. 126 628Google Scholar

    [27]

    Jacka S D, Mijatović A, Širaj D 2014 Stoch. Proc. Appl. 124 1055Google Scholar

    [28]

    Atzberger P J 2006 Phys. Lett. A 351 225Google Scholar

    [29]

    Guo W, Du L C, Liu Z Z, Yang H, Mei D C 2017 Chin. Phys. B 26 010502Google Scholar

    [30]

    von Gehlen S, Evstigneev M, Reimann P 2009 Phys. Rev. E 79 031114Google Scholar

    [31]

    Saikia S, Mahato M C 2010 Physica A 389 4052Google Scholar

    [32]

    von Gehlen S, Evstigneev M, Reimann P 2008 Phys. Rev. E 77 031106Google Scholar

  • [1] 邓红梅, 黄磊, 李静, 陆叶, 李传起. 基于石墨烯加载的不对称纳米天线对的表面等离激元单向耦合器. 物理学报, 2017, 66(14): 145201. doi: 10.7498/aps.66.145201
    [2] 任芮彬, 刘德浩, 王传毅, 罗懋康. 时间非对称外力驱动分数阶布朗马达的定向输运. 物理学报, 2015, 64(9): 090505. doi: 10.7498/aps.64.090505
    [3] 秦天奇, 王飞, 杨博, 罗懋康. 带反馈的分数阶耦合布朗马达的定向输运. 物理学报, 2015, 64(12): 120501. doi: 10.7498/aps.64.120501
    [4] 王飞, 谢天婷, 邓翠, 罗懋康. 系统非对称性及记忆性对布朗马达输运行为的影响. 物理学报, 2014, 63(16): 160502. doi: 10.7498/aps.63.160502
    [5] 戴栋, 王其明, 郝艳捧. 大气压氦气介质阻挡放电中的周期一不对称放电实验研究. 物理学报, 2013, 62(13): 135204. doi: 10.7498/aps.62.135204
    [6] 王飞, 邓翠, 屠浙, 马洪. 耦合分数阶布朗马达在非对称势中的输运. 物理学报, 2013, 62(4): 040501. doi: 10.7498/aps.62.040501
    [7] 王莉芳, 高天附, 黄仁忠, 郑玉祥. 外力作用下反馈耦合布朗棘轮的定向输运. 物理学报, 2013, 62(7): 070502. doi: 10.7498/aps.62.070502
    [8] 陶锋, 陈伟中, 许文, 都思丹. 基于非线性超传导的能流不对称传输现象的研究. 物理学报, 2012, 61(13): 134103. doi: 10.7498/aps.61.134103
    [9] 陈艳, 蒋敏强, 戴兰宏. 金属玻璃温度依赖的拉压屈服不对称研究. 物理学报, 2012, 61(3): 036201. doi: 10.7498/aps.61.036201
    [10] 李体俊. 坐标算符本征矢的表示与不对称投影算符的积分. 物理学报, 2008, 57(7): 3969-3972. doi: 10.7498/aps.57.3969
    [11] 李 微, 赵同军, 郭鸿涌, 纪 青, 展 永. 布朗马达的非均匀高斯跃迁模型. 物理学报, 2004, 53(11): 3684-3689. doi: 10.7498/aps.53.3684
    [12] 缪中林, 陈平平, 蔡玮颖, 李志锋, 徐文兰, 袁先漳, 刘平, 史国良, 陈昌明, 朱德彰, 潘浩昌, 胡军, 李明乾, 陆卫. 组合注入质子导致不对称耦合双量子阱截面混合效应研究. 物理学报, 2001, 50(1): 116-119. doi: 10.7498/aps.50.116
    [13] 张彦鹏, 唐天同, 付盘铭. 四能级系统极化拍频的不对称特性. 物理学报, 1999, 48(2): 242-249. doi: 10.7498/aps.48.242
    [14] 展永, 包景东, 卓益忠, 吴锡真. 布朗马达的定向输运模型. 物理学报, 1997, 46(10): 1880-1887. doi: 10.7498/aps.46.1880
    [15] 王德宁, 潘慧珍. 四层不对称波导的一些近似关系式. 物理学报, 1984, 33(11): 1610-1618. doi: 10.7498/aps.33.1610
    [16] 陈启洲, 胡宁. 奇异粒子衰变的上下不对称问题. 物理学报, 1964, 20(4): 374-377. doi: 10.7498/aps.20.374
    [17] 胡宁. Λ和∑粒子衰变的上下不对称性. 物理学报, 1961, 17(7): 315-320. doi: 10.7498/aps.17.315
    [18] 江安才. 直线式不对称三原分子之振动转动光谱及其势能函数. 物理学报, 1944, 5(1): 49-63. doi: 10.7498/aps.5.49
    [19] 胡宁. 一排等距平行直杆后及方格后激流平均速度及温度之分布. 物理学报, 1944, 5(1): 30-48. doi: 10.7498/aps.5.30
    [20] 胡宁. 圆柱体后及轴对称体后激流平均速度及温度之分布. 物理学报, 1944, 5(1): 1-29. doi: 10.7498/aps.5.1
计量
  • 文章访问数:  2721
  • PDF下载量:  79
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-24
  • 修回日期:  2020-03-09
  • 刊出日期:  2020-05-20

摩擦不对称耦合布朗马达的定向输运

  • 1. 沈阳师范大学物理科学与技术学院, 沈阳 110034
  • 2. 华南师范大学物理与电信工程学院, 广州 510006
  • 3. 华侨大学信息科学与工程学院, 厦门 361021
  • 通信作者: 高天附, tianfugao@synu.edu.cn ; 郑志刚, zgzheng@hqu.edu.cn
    基金项目: 国家级-国家自然科学基金(11875135)

摘要: 本文研究了耦合布朗马达在溶液摩擦对称性破缺条件下的定向输运问题. 详细讨论了溶液阻尼系数比, 周期外力以及耦合作用对布朗粒子定向输运的影响. 研究发现: 溶液摩擦的不对称性不仅能促进摩擦棘轮的定向输运, 且合适的溶液阻尼系数比还可以使耦合棘轮的定向输运达到最强. 此外, 一定摩擦阻尼条件下外力振幅可以使耦合粒子的定向输运呈现多峰结构. 研究还发现合适的自由长度和耦合强度也能增强摩擦棘轮的输运特性. 研究结果不仅能够启发实验上选取合适的摩擦阻尼来优化布朗马达的输运, 也可为信息容量大、反应速度快的纳米器件的制备提供一定的理论参考.

English Abstract

    • 生物分子马达广泛存在于细胞内部, 它们可将生物体内部的化学能转化为机械能, 并利用细胞内的自由能完成各种输运[1,2]. 例如, 细胞内的F-actin马达、肌球蛋白马达、交联蛋白等它们能够产生机械运动[3], 而动力蛋白和驱动蛋白还能拖动囊泡进行步进并完成细胞内的货物运输[4]. 生物分子马达的这种定向运动对于细胞内的物质运输起着非常重要的作用. 为了深入研究分子马达的运动机制, 人们建立了摇摆棘轮[5-8], 闪烁棘轮[9-12]等大量棘轮模型[13]. 通过研究各类棘轮不仅能够深入理解布朗粒子的定向运动行为[14], 找到其定向输运的优化条件[15], 研究结果还可应用于制备纳米器件等领域.

      目前, 关于布朗棘轮定向输运的研究已引起人们广泛的兴趣[16-19]. 如Hojo等[20]通过研究由气泡状中心束缚力驱动的布朗马达, 发现驱动力作用下的马达输运会优于通常单分子马达的输运行为. 此外, Pattanayak等[21]在理论上研究了周期性障碍阵列作用下布朗粒子的输运, 发现障碍物阵列的周期性分布能促进棘轮的定向输运. 同时, Yan等[22]还研究了反馈脉冲棘轮模型, 结果表明粒子间的自由长度和耦合强度都能促进棘轮的定向输运. 然而, 上述各类棘轮的研究中大都考虑的是单位阻尼条件下布朗粒子的运动情况. 由于生物体内细胞液的浓度、杂质等都会影响溶液的实际环境[23-25], 因此分子马达受到的介质阻尼通常是变化的. 此外, 大量实验研究已表明布朗马达的运动都是通过耦合相互作用拖动负载并进行集体定向步进的[26-28]. 由此可见, 研究溶液阻尼对耦合布朗粒子运动的影响更具实际意义. 因此将对不同摩擦阻尼条件下耦合布朗粒子的定向输运展开详细讨论. 通过对不同摩擦阻尼影响的研究, 不仅能深入理解以集体协作形式步进的布朗马达的定向输运特性, 而且对于研究摩擦对称性破缺情况下分子马达的能量转化也具有一定的实验启发.

      迄今为止, 关于溶液阻尼对布朗马达定向输运影响的研究理论上虽有涉及[29], 但并不深入. 如von Gehlen等[30]通过研究无负载作用下摩擦对称性破缺溶液环境中棘轮的输运, 结果发现一定条件下的溶液阻尼会减小棘轮的定向输运, 且通过选择合适的阻尼还能诱导粒子的流反转. 此外, Saikia和Mahato[31]还研究了空间关联摩擦棘轮的输运行为, 结果表明只有在小阻尼条件下布朗粒子才能产生定向运动. 然而, 对于复杂细胞溶液环境中布朗粒子拖动负载时的定向输运问题仍少有研究. 为了在理论上提高不同摩擦阻尼条件下耦合布朗马达拖动负载时的定向输运能力, 深入研究了介质阻尼对反馈耦合布朗粒子定向运动的影响.

      本文通过建立过阻尼摩擦棘轮模型, 研究了阻尼系数比、耦合强度和外力振幅等对耦合布朗棘轮定向输运的影响. 研究发现溶液阻尼并不总是阻碍棘轮的定向输运, 合适的阻尼条件还能使棘轮的定向输运达到最强. 此外, 研究还发现阻尼一定时存在多个外力振幅能使摩擦棘轮的定向输运达到极值, 且合适的外力振幅也可以使棘轮的输运达到最强. 所得结果可用于研究纳米尺度下分子马达的可控制备, 以实现马达做功及能量转换等工作.

    • 主要研究摩擦对称性破缺条件下耦合布朗马达在拖动负载时的运动情况, 其动力学行为可由无量纲化的过阻尼朗之万方程描述[32]:

      $\begin{split} {\gamma _i}{\dot x_i}\left( t \right) =\; & - \partial U\left( {{x_1},{x_2};l} \right)/\partial {x_i} + F\left( t \right) \\ &- \lambda + \sqrt {2{\gamma _i}D} {\xi _i}\left( t \right),\;\;\;(i = 1,2),\end{split}$

      其中t为时间, 两个耦合布朗粒子的位置分别为${x_1}\left( t \right)$, ${x_2}\left( t \right)$. ${\gamma _i}(i = 1, 2)$为第i个粒子的阻尼系数, $\alpha = {\gamma _2}/{\gamma _1}$为两个粒子的溶液阻尼系数之比, 反应的是摩擦阻尼的不对称度. 模型中两个耦合粒子的摩擦系数通常是不同的, 即摩擦不对称. 此外, 虽然方程(1)中两个粒子的摩擦系数${\gamma _i}$即不依赖于时间也不依赖于空间, 但是由于马达两个头部的大小不同, 也会造成粒子受到的摩擦力不等, 因此会导致马达的摩擦对称性破缺. 由于耦合粒子摩擦对称性的破缺会导致一个粒子运动快, 另一个粒子运动慢. 因此, 即使在对称势中耦合棘轮也能够产生定向流. 下文通过调节摩擦阻尼系数比$\alpha$来改变耦合棘轮摩擦的不对称性, 从而研究摩擦棘轮的定向输运行为. 此外, 方程中$F\left( t \right)$为时变外力, $\lambda$为负载, D为热噪声强度, 满足关系$D = {k_{\rm B}}T$, ${k_{\rm B}}$为玻尔兹曼常量, T为溶液环境温度. ${\xi _i}\left( t \right)$为高斯白噪声, 满足如下统计关系:

      $\left\langle {{\xi _i}\left( t \right)} \right\rangle = 0,$

      $ \left\langle {{\xi _i}\left( t \right){\xi _j}\left( {t'} \right)} \right\rangle = {\delta _{ij}}\delta \left( {t - t'} \right),\;\;\;(i,j = 1,2). $

      方程(1)中外势$U\left( {{x_1}, {x_2};l} \right)$的具体形式为

      $\begin{split} & U\left( {{x_1},{x_2};l} \right) = \beta \cdot {V_{\rm{p}}}\left( {{x_i}} \right) + {V_{{\rm{in}}}}\left( {{x_1},{x_2};l} \right),\\ & \qquad \qquad (i = 1,2).\end{split}$

      其中耦合粒子受到的棘轮势${V_{\rm{p}}}\left( {{x_i}} \right)$

      ${V_{\rm{p}}}\left( {{x_i}} \right) = \frac{1}{2}\left[ {1 - \cos \left( {2{\text{π}}{x_i}/L} \right)} \right],$

      L为棘轮势的周期长度. 同时两个粒子还将受到耦合相互作用${V_{{\bf{in}}}}\left( {{x_1}, {x_2};l} \right)$, 其表达式为

      ${V_{\rm in}}\left( {{x_1},{x_2};l} \right) = \frac{k}{2}{\left( {{x_2} - {x_1} - l} \right)^2},$

      式中, k为耦合强度, l为弹簧自由长度, 其结构示意图如图1所示. 图1中曲线代表棘轮势${V_{\rm{p}}}$, 黑色与灰色两个球分别代表耦合粒子${x_1}$${x_2}$, 球的颜色代表粒子受到的不同摩擦力. 两个耦合粒子的相互作用通过弹簧来实现. A, B两组图分别表示棘轮势${V_{\rm{p}}}$处于不同的相位, A组表示耦合粒子处于压缩状态, B组则表示耦合粒子处于伸长状态. 图1中第1行表示棘轮势${V_{\rm{p}}}$存在时, 耦合粒子分别处于压缩和拉伸状态; 第2行表示棘轮势${V_{\rm{p}}}$消失瞬间耦合粒子所处的状态; 第3行表示棘轮势${V_{\rm{p}}}$消失后, 两个耦合粒子恢复到原长状态. 弹簧回到原长的趋势能够促使耦合粒子运动. 由于粒子${x_2}$受到的摩擦力通常小于粒子${x_1}$ ($\alpha < 1.0$), 所以粒子${x_2}$的运动会快于粒子${x_1}$, 进而耦合棘轮能够形成定向输运.

      图  1  摩擦对称性破缺耦合粒子的模型示意图

      Figure 1.  Coupling particles model with broken friction symmetry.

      此外, (4)式中$\beta$为一个重要控制参量, 其具体形式为

      $\beta = \left\{ \begin{aligned} & 1,\qquad G > 0,\\ & 0,\qquad G \leqslant 0. \end{aligned} \right.$

      其中

      $G = \frac{1}{2}\left( { V_{\rm {p}}^\prime \left( {{x_1}} \right) + V_{\rm{p}}^\prime \left( {{x_2}} \right)} \right),$

      表示耦合粒子受到棘轮势的平均作用. 同时, 摩擦棘轮还将受到时变外力$F\left( t \right)$的作用, 其表述为

      $F\left( t \right) = A\sin \left( {\omega t} \right),$

      其中A为外力振幅, $\omega$为圆频率, 外驱动力的周期$\tau = {{2{\text{π}}}}/{\omega }$.

      为了研究负载作用下不同摩擦阻尼对耦合棘轮定向输运的影响, 采用质心平均速度描述耦合棘轮的定向输运行为, 其具体表示为

      $\left\langle v \right\rangle = \frac{1}{2}\left( {\left\langle {{v_1}} \right\rangle + \left\langle {{v_2}} \right\rangle } \right),$

      $ \left\langle {{v_i}} \right\rangle = \mathop {\lim }\limits_{n\tau \to \infty } \frac{1}{{n\tau }}\int_{{t_0}}^{n\tau + {t_0}} {{{\dot x}_i}\left( t \right){\rm{d}}t} ,\;\;i = 1,2. $

      其中$\tau$为周期时间, n为周期数, ${t_0}$为初始时刻, $n\tau$表示耦合粒子的演化时间, $\left\langle \cdot \right\rangle$表示系综平均, $\left\langle {{v_i}} \right\rangle$表示第i个布朗粒子的平均速度.

      采用二阶龙格-库塔算法对方程(1)式进行数值计算, 研究耦合布朗粒子拖动负载时的输运特性随棘轮各参量的变化. 为了得到稳定的系综平均值, 模拟了$1500$个轨道, 每个轨道演化$2000$个周期, 时间步长取$h = 1 \times {10^{ - 3}}$. 文中所有物理量均采用无量纲化参量, 无特殊说明参量取$L = 1.0$, $\omega = {\text{π}}$.

    • 为了研究摩擦棘轮的定向输运特性, 深入讨论了负载作用下摩擦对称性破缺棘轮的输运随各参量变化的行为. 首先, 不同热噪声D条件下, 耦合粒子的质心平均速度随弹簧自由长度l的变化关系如图2所示. 由于外势$U\left( {{x_1}, {x_2};l} \right)$具有平移不变性, 即$U\left( {{x_1}, {x_2};l} \right) \!=\! U\left( {{x_1} \!+\! nL, {x_2} + mL;l - nL \!+\! mL} \right)$, $n, m \in Z$, 且外势还满足反演对称性, 即$U\left( {{x_1}, {x_2};l} \right) = U\left( { - {x_1}, - {x_2}; - l} \right)$. 因此耦合粒子的平均速度随自由长度的变化具有周期性. 这与von Gehlen等[30]研究棘轮输运时得到的结论类似. 为方便起见, 仅研究一个演化周期内耦合粒子的输运随弹簧自由长度的变化情况. 图2所示的研究结果表明, 在一个演化周期内耦合粒子的质心平均速度随自由长度l的变化会出现极值. 如, 当$l = 1.0$时耦合粒子在一个演化周期内弹簧自由长度最长, 此时粒子间的相互作用最弱, 所以耦合布朗粒子容易跨越势垒形成定向运动. 然而, 在一个演化周期内还存在某一个特定的弹簧自由长度使耦合粒子的定向输运达到最小, 类似的分析可知由于此时粒子间的相互作用最强, 使得耦合布朗粒子很难跨越势垒形成定向运动. 由此可见, 通过选取合适的弹簧自由长度可以增强棘轮的输运. 此外, 研究还发现随着热噪声D的增加, 耦合粒子的输运也会增强, 说明一定强度的热噪声也能增强耦合棘轮的定向输运.

      图  2  耦合粒子的平均速度$\left\langle v \right\rangle$随弹簧自由长度l的变化曲线, 其中γ1 = 5.0, k = 20.0, α = 0.3, A = 1.0, λ = 0.1

      Figure 2.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the free length l, where ${\gamma _1} = 5.0$, $k = 20.0$, $\alpha = 0.3$, $A = 1.0$, $\lambda = 0.1$.

    • 通过上文研究已知粒子间的耦合对于棘轮的输运会产生一定影响, 接下来研究了不同外力振幅A作用下耦合强度k对棘轮平均速度的影响, 如图3所示. 研究结果表明, 在小外力振幅条件下, 如$A = 4$时, 随着k的增加, 耦合粒子的几率流单调增加. 然而, 随着外力振幅逐渐变大, 如$A > 4$时, 发现在弱耦合条件下粒子的几率流随k的增加呈非单调的变化关系, 也就是说耦合粒子的定向输运能够产生极值. 有趣的是在较大外力振幅作用下, 如$A \geqslant 15$时, 耦合粒子的几率流会由单峰变成双峰结构, 如图3插图所示. 产生极值的原因主要是由于简谐力与弱耦合作用的相互竞争与协作导致. 由图3还可以发现在振幅$A \geqslant 10$时, 耦合粒子几率流随k变化的曲线都存在峰值, 且随着外力振幅的增加, 曲线峰值对应的优化耦合强度${k_{\rm{opt}}}$也随之增大. 然而, 随着A的增大各曲线对应的优化耦合强度的间隔$\Delta {k_{{\rm{opt}}}}$却越来越小, 因此对于不同A值的曲线能够产生交叉现象. 然而随着耦合强度的继续增大, 即$k \geqslant 30$, 耦合粒子的平均速度逐渐增加并趋于稳定值. 这是因为随着耦合作用的不断增强, 其他外力很难与耦合作用进行竞争, 因此耦合粒子的输运最终会趋于稳定. 由此可见, 在一定外力振幅条件下, 通过选择一个或多个合适的耦合强度能够增强布朗粒子的定向输运. 此外, 还发现随着外力振幅A的增加, 耦合棘轮的几率流整体减小, 这说明较大的外力振幅对耦合棘轮的定向输运还会产生抑制作用. 同时, 图3中外力振幅较大时曲线出现的交叉现象说明外力振幅对耦合棘轮几率流的影响并不是单调的, 所以下文将进一步讨论外力振幅对耦合棘轮定向输运的影响.

      图  3  耦合粒子的平均速度$\left\langle v \right\rangle$随耦合强度k的变化曲线, 插图为$A = 20.0$时的变化曲线, 其中${\gamma _1} = 5.0$, $l = 0.25$, $\alpha = 0.3$, $D = 0.35$, $\lambda = 0.1$

      Figure 3.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the strength of coupling k, the insert figure described the average velocity curve when $A = 20.0$, where ${\gamma _1} = 5.0$, $l = 0.25$, $\alpha = 0.3$, $D = 0.35$, $\lambda = 0.1$.

    • 图3的研究结果已表明外力振幅A对棘轮的输运能够产生抑制作用, 且棘轮的输运随振幅的变化关系并不是单调的. 为了深入研究外力振幅对摩擦棘轮输运的影响, 还详细讨论了不同摩擦阻尼系数比$\alpha$条件下的几率流随A的变化关系. 由图4可知, 随着外力振幅A的增加耦合粒子几率流的变化曲线呈现多峰结构, 且在较大的外力振幅条件下, 如$A > 30$时, 曲线逐渐变得平缓并最终趋于稳定. 产生此现象的原因可由棘轮受到外力作用的分析得到, 即在一定振幅条件下, 如$A < 30$, 此时简谐力会与摩擦阻尼, 耦合作用等其他外力相互协作与竞争, 随着振幅A的增加, 耦合粒子受到的简谐力逐渐增大, 从而棘轮产生的几率流会逐渐变大. 由于几率流的增加会导致耦合粒子受到的黏滞阻尼随之增大. 此外, 由于黏滞阻力还会阻碍耦合粒子的定向运动, 因此当黏滞阻力增大到一定程度时, 耦合棘轮的几率流将会减小, 而减小的几率流又将导致耦合粒子受到的黏滞阻力减小. 当溶液的黏滞阻力减小时耦合粒子的定向输运会再次增强. 因此耦合粒子几率流随外力振幅A的变化曲线能够呈现多个极值. 然而, 随着A的继续增加, 简谐力逐渐变大, 此时耦合粒子受到的主要外力为$F\left( t \right)$, 而其他外力很难再与简谐力竞争, 因此曲线的波动逐渐变小, 进而摩擦棘轮的几率流会变得更加平缓. 上述结果表明, 在一定条件下存在多个合适的外力振幅使摩擦棘轮的质心平均速度达到极值. 此外研究还发现, 随着摩擦阻尼系数比$\alpha$的增加, 摩擦棘轮的输运整体减小, 说明摩擦阻尼对摩擦棘轮的输运也有一定抑制作用. 由此可知, 通过选取合适的外偏置力振幅及溶液阻尼系数比也能增强摩擦棘轮的定向输运. 关于摩擦阻尼对耦合粒子定向输运的影响将在下文进行详细讨论.

      图  4  耦合粒子的平均速度$\left\langle v \right\rangle$随外力振幅A的变化曲线, 其中γ1 = 5.0, l = 0.25, k = 20.0, D = 0.35, λ = 0.1

      Figure 4.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the amplitude A of external force, where ${\gamma _1} = 5.0$, $l = 0.25$, $k = 20.0$, $D = 0.35$, $\lambda = 0.1$.

    • 为了进一步研究溶液阻尼对摩擦棘轮定向输运的影响, 研究了耦合粒子的几率流随溶液阻尼系数比$\alpha = {\gamma _2}/{\gamma _1}$的变化关系. 由图5可得, 随着阻尼系数比的增加, 耦合粒子的定向输运整体减小. 产生这一现象的原因主要是随着$\alpha$的增加, 将导致${\gamma _2}$增加(${\gamma _1}$是常数). 这就意味着$\alpha$越大, 粒子2受到溶液的摩擦阻尼也越大, 因此耦合粒子的几率流会整体变小. 然而, 在一定条件下, 如小负载条件$\lambda < 0.54$, 曲线存在两个不明显的拐点. 随着负载$\lambda$的继续增加, 如$\lambda \geqslant 0.54$, 曲线能够产生两个明显的峰值, 图5中插图仅给出$\lambda = 0.6$时的速度变化曲线. 产生上述峰值现象的原因主要是由于负载较小时, 负载与阻尼的竞争较弱, 随着负载逐渐变大, 如$\lambda \geqslant 0.54$时负载与阻尼的竞争变强, 因而速度曲线的变化会由两个不明显的拐点逐渐变成明显的波峰. 由此可见, 在一定负载条件下通过选取合适的摩擦阻尼系数比也能够增强摩擦棘轮的定向输运.

      图  5  耦合粒子的平均速度$\left\langle v \right\rangle$随摩擦阻尼系数比$\alpha$的变化曲线, 插图为$\lambda = 0.6$时的变化曲线, 其中${\gamma _1} = 5.0$, $k = 20.0$, $l = 0.25$, $A = 20.0$, $D = 0.35$

      Figure 5.  The average velocity $\left\langle v \right\rangle$ of coupled particles varying with the frictional damping coefficient ratio $\alpha$, the insert figure described the average velocity curve when $\lambda = 0.6$, where ${\gamma _1} = 5.0$, $k = 20.0$, $l = 0.25$, $A = 20.0$, $D = 0.35$.

    • 本文主要研究了摩擦对称性破缺条件下耦合布朗粒子的定向输运行为, 同时耦合粒子还会受到外势的反馈作用. 详细讨论了摩擦阻尼系数比、弹簧自由长度、耦合强度以及外力振幅对耦合布朗棘轮定向输运的影响. 研究发现, 一定的自由长度和耦合强度都可以促进棘轮的定向输运. 同时, 有趣地发现在一定条件下摩擦阻尼并不总是阻碍布朗粒子的定向输运, 合适的摩擦阻尼系数比还可增强摩擦棘轮的定向输运特性. 此外, 研究还发现当摩擦阻尼一定时, 耦合粒子的输运随外力振幅的变化呈现多峰结构, 即存在多个外力振幅能使耦合粒子的质心平均速度达到极值, 这也意味着合适的外力振幅可以使摩擦棘轮的定向输运达到最强. 所得结论不仅可应用于研究ATP合成、信号传导等生命活动方面, 还可为微观尺度下纳米器件的研制提供一定的理论指导.

参考文献 (32)

目录

    /

    返回文章
    返回