低惯性与高惯性活性粒子混合体系中的相分离现象.pdf
《低惯性与高惯性活性粒子混合体系中的相分离现象.pdf》由会员分享,可在线阅读,更多相关《低惯性与高惯性活性粒子混合体系中的相分离现象.pdf(11页珍藏版)》请在咨信网上搜索。
1、封面文章低惯性与高惯性活性粒子混合体系中的相分离现象*王晶焦阳田文得陈康(苏州大学物理科学与技术学院,软凝聚态物理及交叉研究中心,苏州215006)(2023年 5月 16 日收到;2023年 7月 28 日收到修改稿)自然界中许多活性物质处于高黏度的环境中,研究者们大多采用过阻尼模型来进行描述.近年来,惯性对活性物质运动行为的影响受到了越来越多的关注.本文通过郎之万动力学模拟,研究低惯性活性粒子和高惯性活性粒子组成的混合体系中的相分离现象.研究结果表明,与单一组分的低惯性粒子体系或高惯性粒子体系相比,中等密度的混合体系有利于相分离现象的发生,而高密度的混合体系不利于相分离的发生.此外,中等密
2、度下,在大量低惯性粒子与少量高惯性粒子组成的混合体系中,由于受到高惯性粒子的撞击,低惯性粒子的速度增大,这有利于低惯性粒子与低惯性粒子形成聚集,从而发生相分离.但是随着低惯性粒子的减少和高惯性粒子的增加,高惯性粒子与高惯性粒子之间的弹性碰撞会抑制相分离现象的产生.本研究有助于深入认识活性物质的集体运动行为,为材料的设计与控制提供参考.关键词:活性物质,相分离,惯性PACS:05.40.Jc,05.40.a,81.16.Dn,82.20.WtDOI:10.7498/aps.72.202307921引言活性物质13是一类可以从周围环境中获取能量来实现自主运动的物质.它们可以展现出独特的非平衡现象4
3、10,例如巨大的密度涨落、动态斑图、相分离、涡旋等,因而在科学界备受关注.多数活性物质体系11,如细菌、细胞、微泳体(micro-swimmer)、自驱动胶体粒子等,都处于高黏度的环境中,受到的流体阻力远大于粒子的惯性作用,研究者们大多采用过阻尼的模型来进行描述,如过阻尼的活性布朗粒子(activeBrownianparticles,ABP)模型1214.近年来,越来越多的研究者们开始关注惯性对活性物质运动的影响15.这类物质处于低黏度的环境中,如振动的机器人16、Runner、飞行的昆虫以及微型飞行器.它们可以用欠阻尼的郎之万动力学模型来描述,也称为活性惯性粒子(activeinertial
4、particles,AIP)模型17.已有的研究结果表明,在一定的密度下,AIP系统中也可以发生相分离1720.但是,与 ABP 系统相比,AIP 系统中惯性粒子之间具有强烈的碰撞反弹效应,不利于粒子的聚集1820,从而抑制相分离的发生.Dai 等17研究发现,惯性的存在导致体系中的团簇比较稀疏,并且在活性惯性团簇态和均匀无规则流体态之间存在一个较大的过渡区域,在该区域,粒子的局部密度很高,但没有形成六角密排.此外,在上述活性惯性粒子的体系中,尽管粒子的形状是各向同性的,但在相分离的稳态和均匀无规流体态之间依然出现了有趣的振荡行为,这一现象说明在活性团簇的内部存在动量转移的现象,*国家自然科学
5、基金(批准号:21774091,21674078)资助的课题.通信作者.E-mail:通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.19(2023)190501190501-1va这是惯性所诱导的集体运动.另外,Mandal 等21研究了 ABP 系统和 AIP 系统中粒子的热力学温度.在 ABP 系统中,粒子的热力学温度没有明显的差别,但是在 AIP 系统,活性团簇区域与均匀无规流体区域的粒子具有明显的热力学温差.这是由于在高密度的团簇区域,惯性粒子之间存在强烈的弹性碰撞,导致粒子将永
6、远不会拥有沿其取向的速度;而在低密度的流体区域,粒子之间的碰撞较少,粒子可以持续加速,直到它们几乎都达到自身的驱动速度 .因此,就粒子的动力学温度而言,团簇区是“冷”的,流体区是“热”的.在现实情况下,活性物质的生存环境比较复杂2226,许多体系中物质的种类并不是单一的,而是由多种成分混合而成.少数关于活性物质混合体系的研究表明,体系的成分及其比例会对相分离产生影响.在活性与非活性粒子的混合体系27中,研究者们发现,这种混合体系不利于活性诱导的相分离现象的产生,体系中形成的团簇在边缘是活性的,内部是非活性的.与单一体系相比,团簇的分裂与融合比较频繁.在低活性与高活性粒子混合体系的模拟28研究中
7、,当两种粒子的活性差较大时,在致密的团簇内部可以观察到高速区域和低速区域.与高速粒子和低速粒子的单一体系相比,高速粒子进入致密相受到抑制,而低速粒子进入致密相得到促进,随着低速粒子活性的增加,这些影响逐渐减弱.在活性粒子和非活性棒状粒子混合体系的模拟29研究中,研究者们发现该体系中也可以发生相分离,但是由于棒状粒子的形状各向异性,导致体系的相分离能力较弱,形成的团簇比较稀疏.受上述研究的启发,构建了一个二元混合体系,该体系由低惯性和高惯性的两种活性粒子组成.感兴趣的是,与单一的低惯性活性粒子体系或高惯性活性粒子体系相比,混合体系中的相分离现象是更容易发生还是会被抑制?这一研究将为我们在复杂环境
8、下调控活性物质的集体行为提供基础30,31.I=m/1/Dr本文通过郎之万动力学模拟,首先研究了单个活性惯性粒子的运动行为(补充材料 1(online).结果表明,单粒子的运动行为由粒子运动的惯性特征时间()和活性力的弛豫时间()共同决定.然后研究了低惯性和高惯性活性粒子混合体系中的相分离现象.研究发现,与单一组分的低惯性粒子体系或高惯性粒子体系相比,中等密度的混合体系有利于相分离现象的发生,而高密度的混合体系不利于相分离的发生.此外,在混合体系中,两种粒子的比例对体系的相分离程度起着至关重要的影响.在大量低惯性粒子与少量高惯性粒子组成的混合体系中,由于受到高惯性粒子的撞击,低惯性粒子的速度增
9、大,这有利于相分离现象的产生.但是随着高惯性粒子比例的增大,高惯性粒子与高惯性粒子之间的弹性碰撞占据主导地位,这种激烈的碰撞会抑制相分离现象.2体系与模型N1N2N=N1+N2通过郎之万动力学模拟研究了二元混合体系中的相分离现象.二元混合体系由 个低惯性粒子和 个高惯性粒子组成,体系中的总粒子数为N,.两种粒子的直径均为 (图 1(b).所有粒子之间的相互作用遵循纯排斥的 Weeks-Chandler-Anderson(WCA)势,其表达式为UWCA(r)=4(/r)12(/r)6+,rrc,0,r rc,(1)rrc=21/6其中 为相互作用强度,为两粒子之间的距离,为相互作用势的截断距离.
10、粒子的运动遵循郎之万动力学方程,表达式为m ri=iUi ri+Fa ui(t)+2kBTi(t),(2)i=2Dri(t).(3)riUiikBTFava=Fa/i ui=(cosi,sini)ib uiDr1/Dri(t)i(t)其中 和 分别是第 个粒子的位置矢量和总势能.为平动摩擦系数,为标度温度.所有粒子自驱动力(活性力)的大小均相同,为 .粒子自驱动速度的大小也相同,为 .但自驱动速度的方向是独立的,第 个粒子的自驱动方向为,其中 为粒子取向与 x 轴之间的夹角.自驱动方向 变化的快慢由旋转扩散系数 所决定,代表粒子定向驱动的持久时间.和 为单位方差的高斯白噪声.方程(2)描述了粒
11、子在平动自由度上的运动,方程(3)描述了粒子自驱动方向的耦合旋转动力学.m1=m12/=10Fa=10kBT=1Dr=0.001m1m2=100 m2=100m1本文使用约化单位,以 ,和 作为质量、长度和能量的基本单位.因此,时间单位为 ,速度单位为 ,力的单位为 .在模拟中,设置 ,.低惯性粒子的质量为 ,高惯性粒子的质量为().粒子的质量决定粒子运物理学报ActaPhys.Sin.Vol.72,No.19(2023)190501190501-2I=m/lI=IvaLxLy=200 200=N2/(4LxLy)动的惯性特征时间 和惯性特征长度.模拟盒子采用周期性的边界条件,其大小为 .体系
12、中粒子的面密度为.3结果与讨论3.1 相图与机制N1:N2研究不同组分比例()下,低惯性和高惯性活性粒子混合体系中的相分离现象.图 1(a)中的相边界是依据粒子局部密度和格点局部密度=0.1N1:N2=10:0N1:N2=0:10=0.3 0.4判断的(具体内容见 3.2 节).在低密度()的情况下,单一组分的低惯性粒子体系()、高惯性粒子体系()以及两种粒子组成的混合体系均呈现均匀无规流体态(图 1).在中等密度()下,单一组分的体系呈现均匀无规流体态;然而,在大量低惯性粒子和少量高惯性粒子组成的混合体系中却出现了明显的相分离现象,这说明该组分比例有利于粒子的聚集.在高密度()的情况下,单一
13、组分的体系和混合体系中都可以产生宏观相分离现象.1/Dr为了探究相分离现象发生的物理机制,通过最简单的粒子对碰撞模型研究了活性惯性粒子间的碰撞(补充材料 24(online).当两个低惯性粒子发生碰撞(图 2(a)后,粒子会在排斥力的作用下分开一定的距离,但不会立即远离,粒子相互远离所用的时间由活性力的弛豫时间()决定.如果在此时间内,有其他粒子靠近、并包围粒子对,则这些粒子就会在该位置形成聚集.当两个高惯性粒子发生碰撞(图 2(c)后,由于强烈的弹性碰撞,粒子对迅速分离,较难形成聚集.当低惯性粒子与高惯性粒子发生碰撞(图 2(b)后,低惯性粒子由10091827364554637281901
14、0120.60.50.40.30.20.1(a)0.60.50.40.30.20.110091827364554637281901012(b)1=12=100图1(a)模拟相图,黑色空心正方形代表均匀无规流体相,黑色实心圆形代表均匀无规流体与团簇的共存相,黑色虚线表示两相的相边界,背景颜色代表混合体系中的粒子组分比(颜色与模拟相图相对应);(b)快照相图,粒子的颜色代表其质量(红色粒子为低惯性粒子,蓝色粒子为高惯性粒子)Fig.1.(a)Simulated phase diagram,the black hollowsquarerepresentsthehomogeneousrandomflu
15、idphase,theblacksolidcirclerepresentsthecoexistingphaseoftheho-mogeneousrandomfluidandthecluster,theblackdashedlinerepresentsthephaseboundaryofthetwophases,whichis divided according to the local density distribution,thebackgroundcolorrepresentstheparticlecompositionratiointhemixingsystem,andthecolor
16、iscorrespondingtothesimulatedphasediagram;(b)phasediagramofsnapshoots,thecolorofparticlesrepresentstheirmass,theredparticlesare low inertia particles,and the blue particles are highinertiaparticles.(a)(b)(c)m1m2图2粒子碰撞示意图.红色球表示质量为 的低惯性粒子,蓝色球表示质量为 的高惯性粒子,深色球表示碰撞前的粒子,浅色球表示碰撞后的粒子,黑色箭头表示碰撞前粒子的运动轨迹,浅红色箭头表
17、示碰撞后低惯性粒子的运动轨迹,浅蓝色箭头表示碰撞后高惯性粒子的运动轨迹m1m2Fig.2.Schematic diagram of collisions between particles,theredspheresrepresentlowinertiaparticleswith ,thebluespheresrepresenthighinertialparticleswith ,thedarkspheresrepresenttheparticlesbeforethecollision,thelightspheresrepresenttheparticlesafterthecollision,
18、theblackarrowsrepresentthetrajectoryoftheparticlesbeforethecollision,thelightredarrowsrepresentthetrajectoryofthelowinertiaparticlesafterthecollision,thelightbluearrows represent the trajectory of the highly inertialparticlesafterthecollision.物理学报ActaPhys.Sin.Vol.72,No.19(2023)190501190501-3于受到强烈的撞击
19、,迅速改变运动方向,而高惯性粒子的运动方向几乎保持不变.N1:N2=0.3=0.6N1:N2在低惯性粒子与高惯性粒子组成的混合体系中,这 3 种碰撞都会发生.并且,随着粒子组分比例()的变化,这 3 种碰撞发生的概率也会随之变化.中等密度()下,在大量低惯性粒子和少量高惯性粒子组成的混合体系中,高惯性粒子与低惯性粒子之间的碰撞会使得低惯性粒子的速度增大,从而导致低惯性粒子与低惯性粒子的碰撞概率增大,有利于粒子的聚集,促进相分离现象的产生.而在少量低惯性粒子和大量高惯性粒子组成的混合体系中,高惯性粒子间强烈的弹性碰撞占主导地位,使得高惯性粒子的速度较低,难以诱导低惯性粒子产生聚集.在高密度()的
20、混合体系中,低惯性粒子的减少与高惯性粒子的增加(减小)依然不利于相分离现象的发生(图 5(b).3.2 结构P(n)P(i,j)P(n)Ladd 10Snnnn体系中是否发生相分离是根据粒子局部密度分布 和格点局部密度分布 判断的.首先计算粒子的局部密度分布 ,计算方法如下.根据二维平面上粒子的位置得到沃罗诺伊图(Voronoidiagram),图中包含多个沃罗诺伊单元,每个单元中仅包含一个粒子,并且该单元内任意一点到该粒子的距离比到其他粒子的距离都要近.另外,在划分沃罗诺伊单元之前,需要将模拟区域内边界附近的粒子(粒子距离模拟平面边界的垂直距离 )复制到模拟区域的外侧,使得模拟区域内边界附近
21、粒子的沃罗诺伊单元满足周期性边界条件.表示第 个粒子所在沃罗诺伊单元的面积.表示第 个粒子的粒子局部密度,表达式为n=2/(4Sn).(4)=0.3P(n)N1:N2=9:1N1:N2=8:2P(n)P(n)在中等密度()下,单一组分的低惯性粒子体系和高惯性粒子体系的粒子局部密度分布 均呈现单峰状(图 4(a),表明体系处于均匀无规流体态.然而,在 和 的混合体系中,粒子局部密度分布 呈现双峰状,表明体系处于均匀无规流体与团簇的共存态.这说明,与单一组分的体系相比,中等密度的混合体系有利于相分离现象的发生.并且,第 1 个峰(代表低密度的气态区域)的峰值较高,=0.6P(n)第 2 个峰(代表
22、高密度的团簇区域)的峰值较低.这表明在发生相分离的混合体系中,少量粒子形成聚集,而大量粒子仍然处于气态区.在高密度()的情况下,与低惯性或高惯性的体系相比,混合体系中 的第 2 个峰值较低(图 4(b),这说明高密度的混合体系会抑制相分离现象的发生.=0.3P(n)=0.6P(n)P(n)N1:N2P(i,j)1 1(i,j)在中等密度()的体系中,粒子局部密度分布 可以很好地反映体系的相分离程度.但是,在高密度()的体系中,由于气态区域的粒子数量稀少,的第一个峰较低,导致 较难全面反映粒子比例()的变化对相分离的影响.所以,进一步计算了格点局部密度分布 .该分布是在沃罗诺伊图的基础上,将二维
23、模拟平面划分为 的格点单元.格点 的密度为(i,j)=noverlapnwn,(5)wn(i,j)P(i,j)P(n)=0.3N1:N2=8:2P(n)P(i,j)=0.6P(i,j)P(i,j)其中,为格点 与沃罗诺伊单元相交区域的面积.因此,代表空间权重的密度分布,而 代表粒子权重的密度分布.在中等密度()、的混合体系中,与粒子局部密度分布 类似,格点局部密度分布 也呈现双峰状(图 5(a),但其第 2 个峰的峰值较低,表明体系中形成的团簇较小.在高密度()的体系中,格点局部密度分布 呈现双峰状(图 5(b),双峰分别代表局部密度低的气态区域和局部密度高的团簇区域(图 3(e).与单一组分
24、的体系相比,混合体系中 的双峰由两侧向中间位置移动,这说明混合体系不利于相分离现象的发生.P(bn)bnP(bn)=0.3bn=6bn=5bn=7由沃罗诺伊图可知,在高密度的团簇内部,粒子沃罗诺伊单元的面积很小(图 3(a),且形状为规则的正六边形,表明团簇内部粒子有序排列.然而,在低密度的气态区域,粒子沃罗诺伊单元的面积很大(图 3(b),且形状不规则,表明气态区域的粒子无规排列.根据这一特点,计算了近邻粒子数分布.为第 n 个粒子所在沃罗诺伊单元的顶点个数,也称为近邻粒子数.由近邻粒子数分布 可知,在中等密度()下,近邻粒子数为 的粒子占比略高于 和 的粒子(图 6(a).这表明,在中等密
25、度下,体系中形成的团簇较小,缺陷较多.并且,与单一组分的体系相物理学报ActaPhys.Sin.Vol.72,No.19(2023)190501190501-4=0.6比,在大量低惯性粒子与少量高惯性粒子组成的混合体系中团簇内部粒子的排列更加密集.这说明,中等密度下,少量高惯性粒子的存在可以诱导低惯性粒子产生聚集,进而发生相分离.在高密度()下,体系中都可以发生相分离,形成致密的大团簇.然而,与单一组分的体系相比,混合体系不利于相分离现象的发生(图 6(b).这是因为,在高密度下,混合体系的团簇内部,高惯性粒子与低惯性粒子之间的碰撞较为频繁,使得团簇容易分裂,从而抑制相分离现象的发生.0408
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 惯性 活性 粒子 混合 体系 中的 分离 现象
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。