搜索

x

留言板

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

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

低惯性与高惯性活性粒子混合体系中的相分离现象

王晶 焦阳 田文得 陈康

引用本文:
Citation:

低惯性与高惯性活性粒子混合体系中的相分离现象

王晶, 焦阳, 田文得, 陈康

Phase separation phenomenon in mixed system composed of low- and high-inertia active particles

Wang Jing, Jiao Yang, Tian Wen-De, Chen Kang
PDF
HTML
导出引用
  • 自然界中许多活性物质处于高黏度的环境中, 研究者们大多采用过阻尼模型来进行描述. 近年来, 惯性对活性物质运动行为的影响受到了越来越多的关注. 本文通过郎之万动力学模拟, 研究低惯性活性粒子和高惯性活性粒子组成的混合体系中的相分离现象. 研究结果表明, 与单一组分的低惯性粒子体系或高惯性粒子体系相比, 中等密度的混合体系有利于相分离现象的发生, 而高密度的混合体系不利于相分离的发生. 此外, 中等密度下, 在大量低惯性粒子与少量高惯性粒子组成的混合体系中, 由于受到高惯性粒子的撞击, 低惯性粒子的速度增大, 这有利于低惯性粒子与低惯性粒子形成聚集, 从而发生相分离. 但是随着低惯性粒子的减少和高惯性粒子的增加, 高惯性粒子与高惯性粒子之间的弹性碰撞会抑制相分离现象的产生. 本研究有助于深入认识活性物质的集体运动行为, 为材料的设计与控制提供参考.
    Active matter refers to a class of substance capable of autonomously moving by harnessing energy from its surrounding environment. The substance exhibits unique non-equilibrium phenomenon, and hence has attracted great attention in the scientific community. Many active matters, such as bacteria, cells, micro-swimmers, and self-propelled colloidal particles, operate in viscous environments and their motions are described usually by using overdamped models. Examples include overdamped active Brownian particle (ABP) model for self-propelled colloidal particles in solution and run-and-tumble (RTP) model for swimming bacteria. In recent years, increasing studies focus on the influence of inertia on the behavior of active matter. Vibrating robots, runners, flying insects, and micro-fliers are typical of active systems under the underdamped condition. The motions of these active matters can be modelled by underdamped Langevin equation, known as the active inertial particle (AIP) model. Previous studies have demonstrated that like the scenarios in ABP systems, motility-induced phase separation (MIPS) phenomena also happen in AIP systems under certain density conditions. However, due to the strong collision-and-rebound effect, aggregation of AIP particles and hence the MIPS are impeded. In complex living/application environments, mixture of different active agents is often seen. Some studies on mixed systems of active matter show that the composition is an important quantity, which influences the phase separation phenomena. In this paper, we study the phase separation phenomena in a mixed system composed of low- and high-inertia active particles by underdamped Langevin dynamics simulations. We find that compared with single-component system, the mixed system is unexpectedly favorable for the occurrence of phase separation at a moderate overall concentration and a certain range of component fraction, while unfavorable for phase separation at a high overall concentration. The underlying mechanism is that the presence of a small number of the high-inertia particles could accelerate the motion of the low-inertia particles, thus facilitating their aggregation and promoting the phase separation. However, when the fraction of the high-inertia particles is large, frequent elastic collisions would disturb the aggregation of the low-inertia particles and suppress the occurrence of phase separation. Our results provide a new insight into the collective behavior of active materials and also a reference for their design and applications.
      通信作者: 田文得, tianwende@suda.edu.cn ; 陈康, kangchen@suda.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 21774091, 21674078)资助的课题.
      Corresponding author: Tian Wen-De, tianwende@suda.edu.cn ; Chen Kang, kangchen@suda.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 21774091, 21674078).
    [1]

    Dauchot O, Löwen H 2019 J. Chem. Phys. 151 114901Google Scholar

    [2]

    Ramaswamy S 2010 Annu. Rev. Condens. Matter Phys. 1 323Google Scholar

    [3]

    Ouellette N T 2019 Matter 1 297Google Scholar

    [4]

    Wensink H H, Löwen H 2008 Phys. Rev. E 78 20Google Scholar

    [5]

    Tailleur J, Cates M E 2008 Phys. Rev. Lett. 100 218103Google Scholar

    [6]

    Peruani F, Ginelli F, Bär M, Chaté H 2011 J. Phys. Conf. Ser. 297 12014Google Scholar

    [7]

    Toner J, Tu Y 1995 Phys. Rev. Lett. 75 4326Google Scholar

    [8]

    Stenhammar J, Marenduzzo D, Allen R J, Cates M E 2014 Soft Matter 10 1489Google Scholar

    [9]

    Liu P, Zhu H W, Zeng Y, et al. 2020 Proc. Natl. Acad. Sci. USA 117 11901Google Scholar

    [10]

    Bricard A, Caussin J B, Das D, et al. 2015 Nat. Commun. 6 7470Google Scholar

    [11]

    Volpe G, Buttinoni I, Vogt D, Kümmerer H J, Bechinger C 2011 Soft Matter 7 8810Google Scholar

    [12]

    Romanczuk P, Bär M, Ebeling W, Lindner B, Schimansky-Geier L 2012 Eur. Phys. J. Spec. Top. 202 1Google Scholar

    [13]

    Speck T 2020 Soft Matter 16 2652Google Scholar

    [14]

    Liao G J, Hall C K, Klapp S H L 2020 Soft Matter 16 2208Google Scholar

    [15]

    Scholz C, Jahanshahi S, Ldov A, Löwen H 2018 Nat. Commun. 9 5156Google Scholar

    [16]

    Li W J, Li L F, Shi Q F, Yang M C, Zheng N 2022 Soft Matter 18 5459Google Scholar

    [17]

    Dai C, Bruss I R, Glotzer S C 2020 Soft Matter 16 2847Google Scholar

    [18]

    Su J, Jiang H, Hou Z 2021 New J. Phys. 23 013005Google Scholar

    [19]

    Löwen H 2020 J. Chem. Phys. 152 040901Google Scholar

    [20]

    Deblais A, Barois T, Guerin T, et al. 2018 Phys. Rev. Lett. 120 188002Google Scholar

    [21]

    Mandal S, Liebchen B, Löwen H 2019 Phys. Rev. Lett. 123 228001Google Scholar

    [22]

    Pan J X, Wei H, Qi M J, Wang H F, Zhang J J, Tian W D, Chen K 2020 Soft Matter 16 5545Google Scholar

    [23]

    Shi S J, Li H S, Feng G Q, Tian W D, Chen K 2020 Phys. Chem. Chem. Phys. 22 14052Google Scholar

    [24]

    Qian B S, Tian W D, Chen K 2021 Phys. Chem. Chem. Phys. 23 20388Google Scholar

    [25]

    Li W J, Li L F, Shi Q F, Yang M C, Zheng N 2022 Powder Technol. 407 117671Google Scholar

    [26]

    Yang Q, Zhu H W, Liu P, Liu R, Shi Q F, Chen K, Zheng N, Ye F F, Yang M C 2021 Phys. Rev. Lett. 126 198001Google Scholar

    [27]

    Stenhammar J, Wittkowski R, Marenduzzo D, Cates M E 2015 Phys. Rev. Lett. 114 018301Google Scholar

    [28]

    Kolb T, Kolb T, Klotsa D 2020 Soft Matter 16 1967Google Scholar

    [29]

    Wang Y, Shen Z L, Xia Y Q, Feng G Q, Tian W D 2020 Chin. Phys. B 29 050103Google Scholar

    [30]

    Popkin G 2016 Nature 529 16Google Scholar

    [31]

    Needleman D, Dogic Z 2017 Nat. Rev. Mater. 2 17048Google Scholar

    [32]

    Shan W J, Zhang F, Tian W D, Chen K 2019 Soft Matter 15 4761Google Scholar

    [33]

    Fily Y, Henkes S, Marchetti C M 2014 Soft Matter 10 2132Google Scholar

  • 图 1  (a)模拟相图, 黑色空心正方形代表均匀无规流体相, 黑色实心圆形代表均匀无规流体与团簇的共存相, 黑色虚线表示两相的相边界, 背景颜色代表混合体系中的粒子组分比(颜色与模拟相图相对应); (b)快照相图, 粒子的颜色代表其质量(红色粒子为低惯性粒子, 蓝色粒子为高惯性粒子)

    Fig. 1.  (a) Simulated phase diagram, the black hollow square represents the homogeneous random fluid phase, the black solid circle represents the coexisting phase of the homogeneous random fluid and the cluster, the black dashed line represents the phase boundary of the two phases, which is divided according to the local density distribution, the background color represents the particle composition ratio in the mixing system, and the color is corresponding to the simulated phase diagram; (b) phase diagram of snapshoots, the color of particles represents their mass, the red particles are low inertia particles, and the blue particles are high inertia particles.

    图 2  粒子碰撞示意图. 红色球表示质量为$ {m}_{1} $的低惯性粒子, 蓝色球表示质量为$ {m}_{2} $的高惯性粒子, 深色球表示碰撞前的粒子, 浅色球表示碰撞后的粒子, 黑色箭头表示碰撞前粒子的运动轨迹, 浅红色箭头表示碰撞后低惯性粒子的运动轨迹, 浅蓝色箭头表示碰撞后高惯性粒子的运动轨迹

    Fig. 2.  Schematic diagram of collisions between particles, the red spheres represent low inertia particles with $ {m}_{1} $, the blue spheres represent high inertial particles with ${m_2}$, the dark spheres represent the particles before the collision, the light spheres represent the particles after the collision, the black arrows represent the trajectory of the particles before the collision, the light red arrows represent the trajectory of the low inertia particles after the collision, the light blue arrows represent the trajectory of the highly inertial particles after the collision.

    图 5  中等密度和高密度体系中的格点局部密度分布 (a) ϕ = 0.3; (b) ϕ = 0.6

    Fig. 5.  Local density distribution of lattices in the system with medium density and high density: (a) ϕ = 0.3; (b) ϕ = 0.6.

    图 4  中等密度和高密度体系中的粒子局部密度分布$ P({\phi _n}) $, Peak 1, Peak 2分别表示混合体系($ {N}_{1}:{N}_{2}=8:2 $)中粒子局部密度分布$ P({\phi _n}) $的第1个峰和第2个峰 (a) ϕ = 0.3; (b) ϕ = 0.6

    Fig. 4.  Local density distribution $ P({\phi }_{n}) $ of particles in the system with medium density (a) and high density (b). Peak 1 and Peak 2 represent the first and second peaks of local density distribution $ P({\phi }_{n}) $ of particles in the mixed system of $ {N}_{1}:{N}_{2}=8:2 $, respectively: (a) ϕ = 0.3; (b) ϕ = 0.6.

    图 3  (a), (b)分别为图(c)中团簇内外所选区域的沃罗诺伊图 (红色粒子为低惯性粒子, 蓝色粒子为高惯性粒子); (c)模拟快照(${N_1}:{N_2} = 8:2$, $ \phi =0.6 $); (d)图(c)所对应的粒子沃罗诺伊单元的面积${S_n}$; (e) 图(c)所对应的格点局部密度$ \phi (i, j) $, 图中颜色代表格点$ (i, j) $的密度

    Fig. 3.  (a), (b) Voronoi diagrams of cluster region and gas region in panel (c), respectively (Red particles are the low inertia particle, and the blue particles are the high inertia particle); (c) simulation snapshot (${N_1}:{N_2} = 8:2$, $ \phi =0.6 $); (d) area ${S_n}$ of the Voronoi unit of the particle corresponding to panel (c); (e) lattice local density $ \phi (i, j) $ corresponding to panel (c); the color in the figure represents the density of the lattice $ (i, j) $.

    图 6  中等密度和高密度体系中的近邻粒子数分布 (a) ϕ = 0.3; (b) ϕ = 0.6

    Fig. 6.  Neighbor particle population distribution in the system with medium density and high density: (a) ϕ = 0.3; (b) ϕ = 0.6.

    图 7  中等密度和高密度的体系中, 两种粒子的粒子局部密度分布$P({\phi _n})$及粒子局部密度的平均值$\overline {{\phi _n}} $. 红色表示低惯性粒子; 蓝色表示高惯性粒子; 黄色表示体系中的所有粒子 (a), (c) ϕ = 0.3; (b), (d) ϕ = 0.6

    Fig. 7.  Local density distribution $P({\phi _n})$ and the average $\overline {{\phi _n}} $ of the local density of the particles with different inertias in the system with medium density and high density. Red represents low inertial particles; the blue represents high inertial particles; the yellow represents all particles in the system: (a), (c) ϕ = 0.3; (b), (d) ϕ = 0.6.

    图 8  中等密度和高密度的体系中, 粒子的近邻粒子类型之比, 纵轴表示粒子周围出现同类(或异类)粒子的概率p. 红色实心正方形表示低惯性粒子附近出现低惯性粒子的概率, 红色空心正方形表示低惯性粒子附近出现高惯性粒子的概率; 蓝色实心圆形表示高惯性粒子附近出现低惯性粒子的概率, 蓝色空心圆形表示高惯性粒子附近出现高惯性粒子的概率; 灰色实心三角形和空心三角形分别表示体系中低惯性粒子和高惯性粒子占总粒子数的比值 (a) ϕ = 0.3; (b) ϕ = 0.6

    Fig. 8.  Ratio of the particle’s nearest neighbor particle types in the system with medium density and high density, and the vertical axis shows the probability p of similar (or dissimilar) particles occurring around the particle. Red solid square indicates the probability of low inertial particles near the low inertial particles, and the red hollow square indicates the probability of high inertial particles near the low inertial particles. Blue solid circle indicates the probability of low inertial particles near high inertial particles, and blue hollow circle indicates the probability of high inertial particles near high inertial particles. Gray solid and hollow triangle represent the ratio of low and high inertial particles to the total number of particles in the system, respectively: (a) ϕ = 0.3; (b) ϕ = 0.6.

    图 9  中等密度(a)和高密度(b)的体系中, 团簇的尺寸分布$P(n)$; (c)平均畴尺寸$\Re $

    Fig. 9.  Size distribution $P(n)$ of clusters in the system with medium density (a) and high density (b); (c) average domain size $\Re $.

    图 10  中等密度和高密度的体系中, cluster区域和gas区域中的粒子类型之比(红色代表cluster区域, 蓝色代表gas区域, 灰色代表整个混合体系) (a) ϕ = 0.3; (b) ϕ = 0.6

    Fig. 10.  Ratio of particle types in the cluster and gas regions in the system with medium density and high density (The red represents the cluster region, blue represents the gas region, and gray represents the whole mixed system): (a) ϕ = 0.3; (b) ϕ = 0.6.

    图 11  (a)热力学温度相图, 粒子颜色代表其热力学温度${T_{{\text{eff}}}}$; (b), (c) cluster区域和gas区域中粒子的平均热力学温度${\overline T _{{\text{eff}}}}$, 低惯性粒子的热力学温度参照左侧刻度线, 高惯性粒子的热力学温度参照右侧刻度线

    Fig. 11.  (a) Thermodynamic temperature phase diagram, in which the color of particle represents its thermodynamic temperature ${T_{{\text{eff}}}}$; (b), (c) the mean thermodynamic temperatures ${\overline T _{{\text{eff}}}}$ of the particles in the cluster and gas regions, in panels (b) and (c), the thermodynamic temperature of low inertia particles is based on the left scale line, and the thermodynamic temperature of high inertia particles is based on the right scale line.

  • [1]

    Dauchot O, Löwen H 2019 J. Chem. Phys. 151 114901Google Scholar

    [2]

    Ramaswamy S 2010 Annu. Rev. Condens. Matter Phys. 1 323Google Scholar

    [3]

    Ouellette N T 2019 Matter 1 297Google Scholar

    [4]

    Wensink H H, Löwen H 2008 Phys. Rev. E 78 20Google Scholar

    [5]

    Tailleur J, Cates M E 2008 Phys. Rev. Lett. 100 218103Google Scholar

    [6]

    Peruani F, Ginelli F, Bär M, Chaté H 2011 J. Phys. Conf. Ser. 297 12014Google Scholar

    [7]

    Toner J, Tu Y 1995 Phys. Rev. Lett. 75 4326Google Scholar

    [8]

    Stenhammar J, Marenduzzo D, Allen R J, Cates M E 2014 Soft Matter 10 1489Google Scholar

    [9]

    Liu P, Zhu H W, Zeng Y, et al. 2020 Proc. Natl. Acad. Sci. USA 117 11901Google Scholar

    [10]

    Bricard A, Caussin J B, Das D, et al. 2015 Nat. Commun. 6 7470Google Scholar

    [11]

    Volpe G, Buttinoni I, Vogt D, Kümmerer H J, Bechinger C 2011 Soft Matter 7 8810Google Scholar

    [12]

    Romanczuk P, Bär M, Ebeling W, Lindner B, Schimansky-Geier L 2012 Eur. Phys. J. Spec. Top. 202 1Google Scholar

    [13]

    Speck T 2020 Soft Matter 16 2652Google Scholar

    [14]

    Liao G J, Hall C K, Klapp S H L 2020 Soft Matter 16 2208Google Scholar

    [15]

    Scholz C, Jahanshahi S, Ldov A, Löwen H 2018 Nat. Commun. 9 5156Google Scholar

    [16]

    Li W J, Li L F, Shi Q F, Yang M C, Zheng N 2022 Soft Matter 18 5459Google Scholar

    [17]

    Dai C, Bruss I R, Glotzer S C 2020 Soft Matter 16 2847Google Scholar

    [18]

    Su J, Jiang H, Hou Z 2021 New J. Phys. 23 013005Google Scholar

    [19]

    Löwen H 2020 J. Chem. Phys. 152 040901Google Scholar

    [20]

    Deblais A, Barois T, Guerin T, et al. 2018 Phys. Rev. Lett. 120 188002Google Scholar

    [21]

    Mandal S, Liebchen B, Löwen H 2019 Phys. Rev. Lett. 123 228001Google Scholar

    [22]

    Pan J X, Wei H, Qi M J, Wang H F, Zhang J J, Tian W D, Chen K 2020 Soft Matter 16 5545Google Scholar

    [23]

    Shi S J, Li H S, Feng G Q, Tian W D, Chen K 2020 Phys. Chem. Chem. Phys. 22 14052Google Scholar

    [24]

    Qian B S, Tian W D, Chen K 2021 Phys. Chem. Chem. Phys. 23 20388Google Scholar

    [25]

    Li W J, Li L F, Shi Q F, Yang M C, Zheng N 2022 Powder Technol. 407 117671Google Scholar

    [26]

    Yang Q, Zhu H W, Liu P, Liu R, Shi Q F, Chen K, Zheng N, Ye F F, Yang M C 2021 Phys. Rev. Lett. 126 198001Google Scholar

    [27]

    Stenhammar J, Wittkowski R, Marenduzzo D, Cates M E 2015 Phys. Rev. Lett. 114 018301Google Scholar

    [28]

    Kolb T, Kolb T, Klotsa D 2020 Soft Matter 16 1967Google Scholar

    [29]

    Wang Y, Shen Z L, Xia Y Q, Feng G Q, Tian W D 2020 Chin. Phys. B 29 050103Google Scholar

    [30]

    Popkin G 2016 Nature 529 16Google Scholar

    [31]

    Needleman D, Dogic Z 2017 Nat. Rev. Mater. 2 17048Google Scholar

    [32]

    Shan W J, Zhang F, Tian W D, Chen K 2019 Soft Matter 15 4761Google Scholar

    [33]

    Fily Y, Henkes S, Marchetti C M 2014 Soft Matter 10 2132Google Scholar

  • [1] 金燕, 石子璇, 金奕扬, 田文得, 张天辉, 陈康. 有限多孔介质诱导活性哑铃的聚集行为. 物理学报, 2024, 73(16): 160502. doi: 10.7498/aps.73.20240784
    [2] 贺华丹, 钟琦超, 解文军. 声悬浮条件下双水相液滴的蒸发与相分离. 物理学报, 2024, 73(3): 034304. doi: 10.7498/aps.73.20230963
    [3] 宁鲁慧, 张雪, 杨明成, 郑宁, 刘鹏, 彭毅. 活性浴中惰性粒子形状对有效作用力的影响. 物理学报, 2024, 73(15): 158202. doi: 10.7498/aps.73.20240650
    [4] 石子璇, 金燕, 金奕扬, 田文得, 张天辉, 陈康. 活性三嵌段共聚物的凝胶化转变. 物理学报, 2024, 73(17): 170501. doi: 10.7498/aps.73.20240796
    [5] 高艺雯, 王影, 田文得, 陈康. 空间调制的驱动外场下活性聚合物的动力学行为. 物理学报, 2022, 71(24): 240501. doi: 10.7498/aps.71.20221367
    [6] 刘博阳, 宋文涛, 刘争晖, 孙晓娟, 王开明, 王亚坤, 张春玉, 陈科蓓, 徐耿钊, 徐科, 黎大兵. AlGaN表面相分离的同位微区荧光光谱和高空间分辨表面电势表征. 物理学报, 2020, 69(12): 127302. doi: 10.7498/aps.69.20200099
    [7] 梁燚然, 梁清. 带电纳米颗粒与相分离的带电生物膜之间相互作用的分子模拟. 物理学报, 2019, 68(2): 028701. doi: 10.7498/aps.68.20181891
    [8] 段华, 李剑锋, 张红东. 二维情况下两组分带电囊泡形变耦合相分离的理论模拟研究. 物理学报, 2018, 67(3): 038701. doi: 10.7498/aps.67.20171740
    [9] 纪丹丹, 张劭光. 三区域膜泡相分离模式之间转变的研究. 物理学报, 2018, 67(18): 188701. doi: 10.7498/aps.67.20180828
    [10] 陈雷鸣. 干活性物质的动力学理论. 物理学报, 2016, 65(18): 186401. doi: 10.7498/aps.65.186401
    [11] 任群, 王楠, 张莉, 王建元, 郑亚萍, 姚文静. 调幅分解及形核对相分离作用机理研究. 物理学报, 2012, 61(19): 196401. doi: 10.7498/aps.61.196401
    [12] 张凤国, 周洪强, 张广财, 洪滔. 惯性及弹塑性效应对延性金属材料层裂损伤的影响. 物理学报, 2011, 60(7): 074601. doi: 10.7498/aps.60.074601
    [13] 王强. Bi0.5Ca0.5Mn1-xCoxO3体系中的电荷有序和相分离. 物理学报, 2010, 59(9): 6569-6574. doi: 10.7498/aps.59.6569
    [14] 李美丽, 付兴烨, 孙宏宁, 赵洪安, 李丛, 段永平, 闫元, 孙民华. 高压作用下相分离液体玻璃转变的分子动力学研究. 物理学报, 2009, 58(8): 5604-5609. doi: 10.7498/aps.58.5604
    [15] 李美丽, 张 迪, 孙宏宁, 付兴烨, 姚秀伟, 李 丛, 段永平, 闫 元, 牟洪臣, 孙民华. 二元Lennard-Jones液体的相分离过程及其扩散性质的分子动力学研究. 物理学报, 2008, 57(11): 7157-7163. doi: 10.7498/aps.57.7157
    [16] 刘 锐, 李寅阊, 厚美瑛. 三维颗粒气体相分离现象. 物理学报, 2008, 57(8): 4660-4666. doi: 10.7498/aps.57.4660
    [17] 翟 薇, 王 楠, 魏炳波. 偏晶溶液相分离过程的实时观测研究. 物理学报, 2007, 56(4): 2353-2358. doi: 10.7498/aps.56.2353
    [18] 蒋中英, 郁伟中, 黄彦君, 夏元复, 马淑新. SMMA/SMA共聚物共混物的自由体积的热动态特性与相分离行为的PALS研究. 物理学报, 2006, 55(6): 3136-3140. doi: 10.7498/aps.55.3136
    [19] 张华力, 刘 卫, 李栋才, 吴修胜, 陈初升. La2NiO4+δ体系相分离现象的低频内耗研究. 物理学报, 2004, 53(11): 3834-3838. doi: 10.7498/aps.53.3834
    [20] 冯文强, 诸跃进. 外噪声场对二元混合物相分离的驱动作用. 物理学报, 2004, 53(11): 3690-3694. doi: 10.7498/aps.53.3690
计量
  • 文章访问数:  4341
  • PDF下载量:  235
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-05-16
  • 修回日期:  2023-07-28
  • 上网日期:  2023-08-02
  • 刊出日期:  2023-10-05

/

返回文章
返回