{"id":66988,"date":"2025-08-08T06:29:13","date_gmt":"2025-08-08T06:29:13","guid":{"rendered":"https:\/\/www.newsbeep.com\/us\/66988\/"},"modified":"2025-08-08T06:29:13","modified_gmt":"2025-08-08T06:29:13","slug":"self-sustained-frictional-cooling-in-active-matter","status":"publish","type":"post","link":"https:\/\/www.newsbeep.com\/us\/66988\/","title":{"rendered":"Self-sustained frictional cooling in active matter"},"content":{"rendered":"<p>Active granular particles governed by dry friction<\/p>\n<p>Active systems governed by dry friction are experimentally explored by utilizing 3D-printed particles placed on a vibrating plate. These plastic objects are designed as cylindrical particles with seven legs attached that are tilted in the same direction (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>a and Methods). The asymmetry of the legs generates active motion when these particles are placed on a vibrating plate, activated by an electromagnetic shaker. Indeed, particles jumping on the plate move in the direction where legs are tilted for a short time, while the direction of motion is generally randomized after a long time due to plate and particle imperfections. The quasi-two-dimensional dynamics of each particle are governed by activity (self-propelled speed) and inertia since particles are macroscopic objects.<\/p>\n<p>As shown in ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Antonov, A. P., Caprini, L., Ldov, A., Scholz, C. &amp; L&#xF6;wen, H. Inertial active matter with Coulomb friction. Phys. Rev. Lett. 133, 198301 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR48\" id=\"ref-link-section-d79981332e629\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a>, the single particle dynamics of these granular objects is subject to dry friction in a range of sufficiently low shaker amplitudes A (see Methods for details). This friction is generated by the contact between plastic legs and plate and impedes the particle motion. The amplitude\u2019s increase enhances the active speed compared to dry friction forces: this allows particles to switch from a Brownian (arrested-like) regime, where dry friction dominates and keeps the particle almost arrested, to a dynamical regime where particles move with a typical speed. In the latter regime, a single particle typically accelerates for a few seconds before changing its direction of motion<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Antonov, A. P., Caprini, L., Ldov, A., Scholz, C. &amp; L&#xF6;wen, H. Inertial active matter with Coulomb friction. Phys. Rev. Lett. 133, 198301 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR48\" id=\"ref-link-section-d79981332e636\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a>.<\/p>\n<p>The principle of self-sustained frictional cooling<\/p>\n<p>In order to get systematic insight into the self-sustained frictional cooling, we consider two active particles governed by dry friction which are characterized by a unique collisional mechanism that does not have an equivalent in wet systems governed by Stokes friction (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>a). Indeed, in the Stokes friction case, a moving active particle is able to push an arrested object, so that both start moving together. By contrast, in the dry friction case, interactions reduce the particle velocity and both particles are at rest after the collision as a result of dry friction. This mechanism is confirmed by proof-of-concept numerical simulations of an activated and a resting particle colliding elastically, i.e. without the loose of kinetic energy due to the collision (see Methods for details). However, after a collision, particles with dry friction lose kinetic energy much faster than particles governed by Stokes friction (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>b). Correspondingly, the center of mass of the system is nearly arrested in the former case while it moves almost linearly with time in the latter case (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>c). Due to fluctuations in the active force, after a transient period, particles start to regain the kinetic energy and move away from each other. However, in a high-density system, particles typically undergo frequent collisions and lack long free runs to restore their kinetic energy. As a consequence, our proof-of-concept analysis focuses on the minimum kinetic energy reached during a collision, which reflects the typical conditions encountered in the collective regime and provides information on the presented cooling mechanism. Here, we identify a range of activity where dry friction generates configurations slower than Stokes friction (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>d). In this range, particles exhibit self-sustained frictional cooling, while larger activities generate heated particles. The unique collisional mechanism responsible for the cooling effect is purely induced by dry friction and activity, and does not depend on the specific collision rule adopted. This is verified in an additional numerical study reported in the Supplementary Information (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">SI<\/a>), where elastic collisions are replaced by an exclusive volume Weeks-Chandler-Andersen (WCA) potential (see Methods for details) or by partially inelastic collisions. Therefore, in the subsequent numerical study, we consider particles interacting via the WCA potential.<\/p>\n<p>Fig. 2: Collisional mechanism for self-sustained frictional cooling.<a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62626-9\/figures\/2\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig2\" src=\"https:\/\/www.newsbeep.com\/us\/wp-content\/uploads\/2025\/08\/41467_2025_62626_Fig2_HTML.png\" alt=\"figure 2\" loading=\"lazy\" width=\"685\" height=\"630\"\/><\/a><\/p>\n<p>a Sketch of a typical binary collision between a moving, activated particle (red) and an arrested particle (blue). A different collision scenario is predicted for Stokes and dry friction. In the latter case, both particles are arrested, while in the former case they continue moving. b, c Averaged mean kinetic energy, \u3008Ekin(t)\u3009 (b), and center of mass of two particles (c) as a function of time t. The activated (left) particle in both cases has a typical activation speed v0. The kinetic energy is normalized by the initial activation energy \\({E}_{{{{\\rm{kin}}}}}(0)=m{v}_{0}^{2}\/2\\). In both cases, after an initial drop in kinetic energy to a minimum due to the collision, the system starts to regain kinetic energy from fluctuations in the active force. d Averaged minimal kinetic energy during the cooling process as a function of activity. For low activity, dry friction cools down the system more than Stokes friction, allowing the kinetic temperature to almost approach zero. The crossover from self-sustained frictional cooling to heating is marked by a vertical dashed line dividing the activity region where cooling by dry friction becomes less effective compared to the reference system with Stokes friction. The simulation protocol to generate the elastic collisions in the case of dry friction and Stokes friction is described in the Methods.<\/p>\n<p>Cooled, mixed and heated phases<\/p>\n<p>Experimentally, collective phenomena are explored by placing N active granular particles on the plate, at packing fraction \u03a6\u00a0=\u00a0Nd2\/D2\u00a0=\u00a00.5, where D and d are the plate\u2019s and the particles\u2019 diameters, respectively. We let the system evolve at large shaker amplitude to reach a configuration where particles are randomly placed on the plate. Successively, we sharply decrease the shaker\u2019s amplitude to the desired value tuning the active speed compared to the dry friction. For low amplitude conditions (corresponding to low activity), particles rarely move until to form a cluster where they are almost arrested, as outlined by the steady-state temporal evolution (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>a and Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM3\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>) and the instantaneous kinetic energy of each particle. We refer to this almost arrested dynamical state as a frictional arrested cooled phase. This dynamical feature is confirmed by plotting the distribution of the particle speed p(v) that is characterized by a narrow peak close to zero (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>d) and by a short-tail for non-zero speeds. This tail corresponds to particles that slightly move within the dense cluster. However, these particles are almost immediately stopped (cooled down) by the caging effect due to neighboring particles, as shown in Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM4\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>, where the particles are colored according to their speed.<\/p>\n<p>Fig. 3: Cooled, mixed and heated phases.<a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62626-9\/figures\/3\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig3\" src=\"https:\/\/www.newsbeep.com\/us\/wp-content\/uploads\/2025\/08\/41467_2025_62626_Fig3_HTML.png\" alt=\"figure 3\" loading=\"lazy\" width=\"685\" height=\"449\"\/><\/a><\/p>\n<p>a\u2013c Snapshots at two different times to outline the time evolution of configurations for cooled (A\u00a0=\u00a018.66\u00a0\u00b1\u00a00.08\u2009\u03bcm), mixed (A\u00a0=\u00a018.88\u00a0\u00b1\u00a00.09\u2009\u03bcm) and heated (A\u00a0=\u00a021.56\u00a0\u00b1\u00a00.09\u2009\u03bcm) phases, respectively. The snapshots in the upper row are the initial configurations, with a corresponding to that in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>g, while the lower row reports the snapshot after 30 seconds of time evolution. In the lower row, the trajectories of two highlighted tracers (one inside and one outside the cluster) are displayed for each phase. Particles with low mobility are represented by orange trajectories while those with high mobility are shown in dark green. The evolution reveals that the cluster remains stable in the cooled and mixed phases but becomes transient in the heated phase due to the significantly higher particle mobility. d\u2013f Velocity probability distributions p(v), shown for the corresponding phases in (a)\u2013(c). The colored solid lines represent experimental data, while the black dashed lines correspond to results obtained by simulating the dynamics (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>) in a confining circular arena. Parameters of the simulations are given in Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Tab1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>. Both the cooled and mixed phases exhibit a velocity peak near zero, with the mobile particles outside the cluster in the mixed state contributing to the distribution tail at higher velocities. In contrast, the heated phase shows a complete shift of the velocity peak away from zero.<\/p>\n<p>By increasing the shaker\u2019s amplitude and, thus, the activity compared to dry friction, particles outside the cluster start moving as evidenced by plotting the kinetic energy per particle and the typical particle trajectories (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>b). In this case, a large cluster formed by cooled particles coexists with fast particles with large kinetic energy (see the steady-state Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM3\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>). While the former particles are part of a cooled cluster, the latter particles define a heated phase. Consequently, in this intermediate regime, the cooled and heated phases coexist. This is confirmed by plotting the speed distribution p(v) (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>e) which shows the coexistence between a narrow peak at zero, which is generated by frictional arrested (cooled) particles in the dense cluster, and a long tail reaching large velocity values due to heated particles in the dilute region. The self-sustained frictional cooling mechanism operates in both regimes but is ultimately weaker in the dilute region, where collisions are infrequent. As a result, it cannot effectively stop the particles, which are instead heated by activity, as shown in the Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM4\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>.<\/p>\n<p>Finally, by increasing the shaker amplitude A (further reducing the dry friction), almost all the particles move fast and form unstable aggregates which continuously break and reform (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>c and Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM3\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>). Given the large kinetic energy value per particle, this regime can be identified as a heated phase, characterized by dynamical clustering. In this phase, the speed distribution p(v) shows a peak at large velocity, being determined by hot particles only (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>f and Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM4\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). In this heated regime, the self-sustained frictional cooling is still present but it is too weak to arrest the particles which continuously move being heated by activity. However, the particle acceleration observed at the single-particle level<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Antonov, A. P., Caprini, L., Ldov, A., Scholz, C. &amp; L&#xF6;wen, H. Inertial active matter with Coulomb friction. Phys. Rev. Lett. 133, 198301 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR48\" id=\"ref-link-section-d79981332e957\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a> is not easily discernible in the collective dynamics due to frequent interparticle collisions.<\/p>\n<p>The cooled, mixed, and heated phases depicted in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a> demonstrate configurations already in the steady state. The corresponding Supplementary Movies\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM3\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a> and <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM4\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a> begin after a long transient, lasting several minutes, to ensure that the system has reached this regime. However, experiments are initialized from a loosely packed configuration and spontaneously evolve toward the aforementioned phases, as shown in Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM5\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>, which captures the system from the very moment the shaker is turned on.<\/p>\n<p>Model for self-sustained frictional cooling<\/p>\n<p>Active Brownian particles with Stokes friction typically used in wet active matter cannot reproduce the cooled and heated phases experimentally observed. Indeed, these models show dynamical clustering at large self-propelled speed (large P\u00e9clet number) compared to our system, and, specifically, cannot reproduce the cooled cluster of almost frictional arrested objects experimentally observed. By contrast, the clustering of our experimental system does not originate from the blocking effect caused by high motility and volume exclusion but intuitively arise from the competition between dry friction and the caging effect due to neighboring particles.<\/p>\n<p>To examine the suggested frictional mechanism and reproduce cooled and heated phases, we perform a numerical study based on inertial active particles with mass m, subject to dry friction. Particles evolve with two-dimensional inertial dynamics for the velocity \\({{{{\\bf{v}}}}}_{i}={\\dot{{{{\\bf{x}}}}}}_{i}\\)<\/p>\n<p>$$m{\\dot{{{{\\bf{v}}}}}}_{i}=-{{{\\boldsymbol{\\sigma }}}}({{{{\\bf{v}}}}}_{i})+\\sqrt{2K}{{{{\\boldsymbol{\\xi }}}}}_{i}(t)+f{{{{\\bf{n}}}}}_{i}+{{{{\\bf{F}}}}}_{i}+{{{{\\bf{F}}}}}_{i}^{w},$$<\/p>\n<p>\n                    (1)\n                <\/p>\n<p>where \u03bei are Gaussian white noises with unit variance and zero average and the constant K determines the noise strength. This noise is due to imperfections in the particle shape and the surface of the plate and is generated by the small vertical motion of the vibrobot due to the oscillating plate<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 33\" title=\"Scholz, C., Jahanshahi, S., Ldov, A. &amp; L&#xF6;wen, H. Inertial delay of self-propelled particles. Nat. Commun. 9, 5156 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR33\" id=\"ref-link-section-d79981332e1236\" rel=\"nofollow noopener\" target=\"_blank\">33<\/a>. The term fni models the active force whose evolution follows the active Ornstein-Uhlenbeck dynamics<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Szamel, G. Self-propelled particle in an external potential: Existence of an effective temperature. Phys. Rev. E 90, 012111 (2014).\" href=\"#ref-CR52\" id=\"ref-link-section-d79981332e1248\">52<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Maggi, C. et al. Generalized energy equipartition in harmonic oscillators driven by active baths. Phys. Rev. Lett. 113, 238303 (2014).\" href=\"#ref-CR53\" id=\"ref-link-section-d79981332e1248_1\">53<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Keta, Y.-E., Jack, R. L. &amp; Berthier, L. Disordered collective motion in dense assemblies of persistent particles. Phys. Rev. Lett. 129, 048002 (2022).\" href=\"#ref-CR54\" id=\"ref-link-section-d79981332e1248_2\">54<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 55\" title=\"Gupta, D., Klapp, S. H. L. &amp; Sivak, D. A. Efficient control protocols for an active Ornstein-Uhlenbeck particle. Phys. Rev. E 108, 024117 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR55\" id=\"ref-link-section-d79981332e1251\" rel=\"nofollow noopener\" target=\"_blank\">55<\/a>. In Eq. (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>), the constant f determines the amplitude of the activity and sets the typical speed of an active granular particle while the stochastic term ni evolves as an Ornstein-Uhlenbeck process<\/p>\n<p>$$\\dot{{{{\\bf{n}}}}}(t)=-\\frac{{{{\\bf{n}}}}(t)}{\\tau }+\\sqrt{\\frac{2}{\\tau }}{{{\\boldsymbol{\\eta }}}}(t),$$<\/p>\n<p>\n                    (2)\n                <\/p>\n<p>and determines the direction of the single-particle motion. In Eq. (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>), \u03b7(t) is a Gaussian white noise with zero average and unit variance and \u03c4 represents the persistence time of the active particle. In our experimental setup, dissipation during collisions is negligible<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 38\" title=\"Caprini, L., Breoni, D., Ldov, A., Scholz, C. &amp; L&#xF6;wen, H. Dynamical clustering and wetting phenomena in inertial active matter. Commun. Phys. 7, 343 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR38\" id=\"ref-link-section-d79981332e1397\" rel=\"nofollow noopener\" target=\"_blank\">38<\/a>. Therefore, interactions are modeled by a conservative force Fi derived from a Weeks-Chandler-Andersen potential (see Methods), which accounts for volume exclusion. In addition, a repulsive force \\({{{{\\bf{F}}}}}_{i}^{w}\\), derived from a harmonic potential truncated in its minimum, represents the confinement imposed by the arena (see Methods for further details). Finally, in Eq. (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>), dry friction is included through the term \u00a0\u2212\u00a0\u03c3(vi) which points in the opposite direction compared to the particle velocity and reads<\/p>\n<p>$${{{\\boldsymbol{\\sigma }}}}({{{\\bf{v}}}})={\\Delta }_{C}\\hat{{{{\\bf{v}}}}},$$<\/p>\n<p>\n                    (3)\n                <\/p>\n<p>where \\(\\hat{{{{\\bf{v}}}}}\\) is the normalized velocity vector which is equal to zero if v\u00a0=\u00a00. This expression models the dynamic dry (Coulomb) friction which decelerates an object already in motion and uniquely depends on the velocity direction via the constant friction coefficient \u0394C. The model (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>) is the paradigm to study dry friction in particle dynamics<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"de Gennes, P.-G. Brownian motion with dry friction. J. Stat. Phys. 119, 953&#x2013;962 (2005).\" href=\"#ref-CR56\" id=\"ref-link-section-d79981332e1562\">56<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Touchette, H., der Straeten, E. V. &amp; Just, W. Brownian motion with dry friction: Fokker-Planck approach. J. Phys. A: Math. Theor. 43, 445002 (2010).\" href=\"#ref-CR57\" id=\"ref-link-section-d79981332e1562_1\">57<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Lequy, T. &amp; Menzel, A. M. Stochastic motion under nonlinear friction representing shear thinning. Phys. Rev. E 108, 064606 (2023).\" href=\"#ref-CR58\" id=\"ref-link-section-d79981332e1562_2\">58<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 59\" title=\"Plati, A., Puglisi, A. &amp; Sarracino, A. Thermodynamic uncertainty relations in the presence of non-linear friction and memory. J. Phys. A: Math. Theor. 57, 155001 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR59\" id=\"ref-link-section-d79981332e1565\" rel=\"nofollow noopener\" target=\"_blank\">59<\/a>, with direct applications in Brownian motors<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 60\" title=\"Baule, A. &amp; Sollich, P. Singular features in noise-induced transport with dry friction. Europhys. Lett. 97, 20001 (2012).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR60\" id=\"ref-link-section-d79981332e1569\" rel=\"nofollow noopener\" target=\"_blank\">60<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 61\" title=\"Manacorda, A., Puglisi, A. &amp; Sarracino, A. Coulomb friction driving Brownian motors. Commun. Theor. Phys. 62, 505 (2014).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR61\" id=\"ref-link-section-d79981332e1572\" rel=\"nofollow noopener\" target=\"_blank\">61<\/a> and passive granular particles<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 62\" title=\"Lema&#xEE;tre, A., Mondal, C., Procaccia, I. &amp; Roy, S. Stress correlations in frictional granular media. Phys. Rev. B 103, 054110 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#ref-CR62\" id=\"ref-link-section-d79981332e1576\" rel=\"nofollow noopener\" target=\"_blank\">62<\/a>. We remark that the suggested model does not include static friction, since numerical checks confirm that its inclusion does not qualitatively alter the observed phenomena (see\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">SI<\/a>).<\/p>\n<p>By simulating the dynamics (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>) in the experimental conditions \u2013 e.g., same number of particles and arena size \u2013 we obtain a qualitative agreement with experimental results. By increasing the particle activity f compared to the dry friction coefficient \u0394C, we observe cooled, mixed, and heated phases as revealed by the Supplementary Movie\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM6\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>. This qualitative match is confirmed by monitoring the speed distribution p(v), which is sharply peaked around zero in the cooled phase (compare solid and dashed lines in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>d) and exhibits a broad shape with a peak at large speed in the heated phase (compare solid and dashed lines in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>f). While in these cases we obtain an excellent agreement, p(v) in the mixed phase shows a shorter tail compared to experiments (compare solid and dashed lines in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>e). This occurs because, in simulations, particles in the dilute region of the mixed phase are typically slower as compared to experiments.<\/p>\n<p>Kinetic phase diagram<\/p>\n<p>Cooled and heated phases observed in experiments are numerically reproduced by simulating the dynamics (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>) in a square box of size L with periodic boundary conditions (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>). This numerical study shows that different temperature phases &#8211; specifically, the mixed phase and the cooled cluster &#8211; are generated by dry friction and caging effects from neighboring particles. In those phases, the particles are mostly stopped by dry friction and move due to rare fluctuations of the active force and translational noise. When particles are in close contact (during interactions), the effect of those rare fluctuations is almost suppressed by the caging imposed by the nearly immobile neighboring particles. This scenario dominates inside the cooled clusters and in the dense region of the mixed phase. However, in the latter case, the particle density outside the cluster is low, thereby still allowing particle to exhibit notable mobility as they rarely interact with each other. This physical mechanism generates the self-sustained frictional cooling experimentally observed.<\/p>\n<p>Fig. 4: Kinetic phase diagram.<a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62626-9\/figures\/4\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig4\" src=\"https:\/\/www.newsbeep.com\/us\/wp-content\/uploads\/2025\/08\/41467_2025_62626_Fig4_HTML.png\" alt=\"figure 4\" loading=\"lazy\" width=\"685\" height=\"267\"\/><\/a><\/p>\n<p>a Phase diagram in the plane of reduced activity f0 and packing fraction \u03a6 with a color gradient denoting the mode particle speed (points). Background colors are used to distinguish between different phases: cooled (blue), mixed (pink), and heated (red) phases. The cooled phase occurs when most of the particles are frictional arrested and remain within the cluster, whereas in a heated phase, particles are highly mobile and are not significantly slowed down by the frictional forces. A state where cooled and heated phases coexist is referred to as a mixed phase. Details on the transition line between these states are discussed in the Methods. b\u2013d Snapshots of cooled, mixed, and heated phases. The color gradient denotes the particle speed (red for high and blue for low speeds). The stars above each snapshot indicate the corresponding parameters f0,\u00a0\u03a6 in the phase diagram a. e\u2013g Probability distribution of the velocity p(v) for the cooled (e), mixed (f) and heated (g) phases. For the mixed state, we separately highlight the distributions for particles inside (blue) and outside (red) of the cooled cluster. Particles within the cluster exhibit characteristics of the cooled phase, while those outside behave like particles in the heated phase. In all cases p(v) exhibits a single peak vm (dashed vertical line, mode speed), with the value of this peak depicted in the phase diagram a.<\/p>\n<p>Within this numerical study, we overcome experimental limitations, such as finite-size effects and boundaries of the plate. As explained in the methods, the dynamics are mainly governed by the reduced activity f0\u00a0=\u00a0f\/\u0394C, which quantifies the active force effect compared to the dry friction: The larger f0, the smaller the dry friction or equivalently, the lower the activity. Simulations start from homogeneous configurations and then analyzed when the system reaches the steady state (see Methods for details). This numerical study is performed across various reduced activities and packing fractions \u03a6\u00a0=\u00a0N\u03c0\u03c32\/(4L2). In this way, we extend the experimental study by systematically exploring a broad range of densities that in equilibrium systems show gas-like configurations and high-density homogeneous liquids.<\/p>\n<p>Temperature phases are systematically explored in a phase diagram in the plane of reduced activity f0 and packing fraction \u03a6 with colors representing the typical particle speed, i.e. the mode of the velocity distribution (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>a). Specifically, for low reduced activity f0 (large dry friction), simulations confirm the cooled phase observed in experiments, where particles are almost arrested in cluster structures (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>b). These particles are characterized by low temperature, i.e. low values of kinetic energy, as revealed by plotting the instantaneous kinetic energy per particle. These dynamical properties reflect onto the speed distribution p(v) which shows a narrow peak at vanishing speed (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>e). This shape of the distribution qualitatively agrees with the experimental one (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>d). A detailed discussion on the choice of mode speed as a parameter to distinguish between the cooled and heated phases is provided in the\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">SI<\/a>.<\/p>\n<p>By increasing the reduced activity f0, the arrested-like cluster is surrounded by a heated phase, consisting of fast particles (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>c) with large kinetic energy as in experiments. In this regime, the speed distribution p(v) has a peak at zero and a long tail for large velocities (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>f). The former is generated by slow particles in the cluster, while the latter is due to fast particles in the heated phase. This interpretation results from the direct calculation of p(v) inside and outside the cluster, which is in qualitative agreement with the experimental results (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>e).<\/p>\n<p>Finally, a further increase of f0 (lower values of dry friction) completely suppresses the almost arrested cluster and generates a heated phase with fast particles (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>d). Again, this qualitative picture is confirmed by measuring p(v) which displays a broad shape with a peak at large speed as in experiments (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>f).<\/p>\n<p>Our phase diagram reveals that the cooled phase is promoted by high packing fractions \u03a6. Indeed, in a dense system, the caging effect is enhanced, thereby strengthening the self-sustained frictional cooling. For low \u03a6, the cooled phase is directly followed by the heated phase, with a smooth crossover reflecting onto the smooth change of the particle skewness. Specifically, in this regime, particles interact rarely and therefore the transition from cooled to heated phases takes place when the activity exceeds the dry friction coefficient for a single particle, i.e. for f0\u00a0\u2248\u00a01. By increasing \u03a6, the transition line shifts to larger activities f0 \u2273 1 because interactions are more frequent and collisions &#8211; governed by the self-sustained frictional cooling mechanism &#8211; on average slow down the particles, thereby favoring the cooled phase. For large \u03a6, specifically above the threshold value \u03a6 &gt; 0.4, the transition is anticipated by the mixed phase, where heated and cooled particles coexist and demix. This scenario resembles a first-order phase transition, occurring out-of-equilibrium, entirely due to the competition between dry friction, caging and activity.<\/p>\n<p>Structure analysis: homogeneous, clustered and solid structures<\/p>\n<p>We combine dynamic and static information to discuss the structural properties of cooled and heated phases. Both are characterized by homogeneous, clustered, and solid-like structures depending on reduced activity f0 and packing fraction \u03a6 (see the phase diagram, Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>a). To distinguish these configurations, we study the distribution Prob(\u03d5) of the local packing fraction \u03d5 (the particle area divided by the area of its Voronoi cell), and, specifically, its skewness (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>f\u2013i), i.e. the degree of asymmetry compared to the average packing fraction (see Methods for details). The homogeneous phase is characterized by an almost symmetric Prob(\u03d5) (vanishing skewness), while the presence of clusters (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>c) induces long tails in Prob(\u03d5) (negative skewness) and shifts the main peak compared to the average packing fraction.<\/p>\n<p>Fig. 5: Structure analysis.<a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62626-9\/figures\/5\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig5\" src=\"https:\/\/www.newsbeep.com\/us\/wp-content\/uploads\/2025\/08\/41467_2025_62626_Fig5_HTML.png\" alt=\"figure 5\" loading=\"lazy\" width=\"685\" height=\"244\"\/><\/a><\/p>\n<p>a Phase diagram illustrating the four possible structural states. The thick solid lines are identical to those shown in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>. Phases are distinguished by comparing the distribution of the local packing fraction \u03d5 and the one of the orientational order parameter \u03c86 (see Methods for details). b, e Snapshots for the different structural phases. Particles are colored according to their speed, i.e. red and blue for high and low speed, respectively. The stars above each snapshot indicate the corresponding parameters f0,\u00a0\u03a6 in the phase diagram a. f\u2013i Probability distributions of the local packing fraction Prob(\u03d5) (violet) and the hexagonal local order parameter Prob(\u03c86) (green). The vertical dashed lines mark the average packing fraction value \u03a6. The snapshots (b)\u2013(e) correspond to the distribution (f)\u2013(i), respectively.<\/p>\n<p>Both in the cooled and heated phases, clustering is intuitively favored by large global packing fraction values \u03a6. However, counterintuitively, clustering occurs non-monotonically with the reduced activity f0 revealing a sharp change when the system switches from cooled to heated phases. Indeed, activity favors cluster formation (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>c) in the cooled phase but promotes homogeneous configurations (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>b) in the heated phase. This non-monotonic behavior is absent in overdamped wet active matter where clustering is always favored by the increase of the active speed and requires speeds at least an order of magnitude larger than our activity. These differences are caused by the origin of the mechanism leading to cluster formation: Indeed, in our system, clustering does not occur because particles block each other but because of self-sustained frictional cooling generated by dry friction and caging. Specifically, in the heated phase, the faster the particle, the larger the probability that dry friction is overcome. When this happens particles likely leave the clusters and the system shows a homogeneous phase. By contrast, in the cooled phase, only particles with high enough speed have the capability of moving until they collide with other particles and eventually remain stuck in a cluster.<\/p>\n<p>In agreement with our expectations, the mixed state, i.e. the coexistence of the cooled cluster and the heated homogeneous state (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>d), is characterized by a bimodal packing fraction distribution Prob(\u03d5) (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>h). This phase coexistence differs from motility-induced phase separation typical of overdamped active matter because i) occurs at small activity values and ii) the two coexisting phases have a different temperature: large temperature for heated configurations and low temperature for cooled configurations. Starting from the mixed phase (pink region in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>a), the phase coexistence is suppressed both for large and low activity values. Indeed, when f0 is increased, a larger number of particles have the capability of moving until the cooled phase is completely suppressed and dynamical clustering is recovered. By contrast, when f0 is low, almost all the particles are almost arrested and the cooled phase dominates over the heated one (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>e). We call this regime, a cooled solid. Indeed, in this phase, the bimodality of Prob(\u03d5) is suppressed (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>i).<\/p>\n<p>To distinguish between cooled clustering and cooled solid, we monitor the distribution of the hexatic order parameter Prob(\u03c86) (see the methods for the definition of \u03c86). This observable is close to the unit for solid-like configurations and returns smaller values otherwise. The solid phase can be identified when Prob(\u03c86) shows a peak at \u03c86\u00a0\u2248\u00a01 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>i), while we consider cooled clustering (or homogeneous) those configurations such that the peak of Prob(\u03c86) occurs at smaller \u03c86 values (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>f, g). Finally, as expected, Prob(\u03c86) is bimodal in the mixed phase (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62626-9#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>h).<\/p>\n","protected":false},"excerpt":{"rendered":"Active granular particles governed by dry friction Active systems governed by dry friction are experimentally explored by utilizing&hellip;\n","protected":false},"author":2,"featured_media":66989,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[49],"tags":[1159,1160,199,79,48978,17232],"class_list":["post-66988","post","type-post","status-publish","format-standard","has-post-thumbnail","category-physics","tag-humanities-and-social-sciences","tag-multidisciplinary","tag-physics","tag-science","tag-soft-materials","tag-statistical-physics"],"_links":{"self":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts\/66988","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/comments?post=66988"}],"version-history":[{"count":0,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts\/66988\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/media\/66989"}],"wp:attachment":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/media?parent=66988"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/categories?post=66988"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/tags?post=66988"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}