Abstract
Collective particle transport across periodic energy landscapes is ubiquitously present in many condensed matter systems spanning from vortices in high-temperature superconductors, frictional atomic sliding, driven skyrmions to biological and active matter. Here we report the emergence of fast solitons propagating against a rotating optical landscape. These experimentally observed solitons are stable cluster waves that originate from a coordinated particle exchange process which occurs when the number of trapped microparticles exceeds the number of potential wells. The size and speed of individual solitons rapidly increase with the particle diameter as predicted by theory and confirmed by numerical simulations. We show that when several solitons coexist, an effective repulsive interaction can stabilize their propagation along the periodic potential. Our experiments demonstrate a generic mechanism for cluster-mediated transport with potential applications to condensed matter systems on different length scales.
Similar content being viewed by others
Introduction
The driven motion of particles through periodic structures is intensively studied with the aim to understand and control non-equilibrium processes in biology, chemistry and physics1,2,3. Particleâparticle and particleâsubstrate interactions yield a wide variety of phenomena as directional locking4,5,6, stepwise increases of ratchet currents7, dynamic mode locking8, or pinning-depinning transition with a complex alternation between static and dynamic modes9,10.
In a recent theoretical work11 it was predicted that a collection of hard spheres driven across a periodic potential could produce solitons at high densities, i.e., stable propagating cluster waves. The predicted effect is appealing, since it allows to achieve a net transportation in dense systems even when single particles cannot surmount the energetic barriers of the potential, thus invoking for an experimental realization.
Here we experimentally show that in a rotating periodic potential with a number of particles larger than the number of potential wells, solitons emerge due to collective effects. In general, solitons are solitary waves propagating without distortion12. They have been found in a variety of systems from fermionic superfluids13 and Bose-Einstein condensate14,15, to macroscopic cracks16, mechanical metamaterials17,18,19, and ocean waves20. In our experiments, the solitons are composed of clusters that continuously break and reform by releasing and accepting particles via a periodic exchange particle process. Solitons are commonly known from a continuous description of the nonlinear dynamics of waves, as, e.g., in the Sine-Gordon, Korteweg-De Vries, or nonlinear Schrödinger equations and others21. In contrast, in our experiments the solitons emerge from many-body effects and are observed at the single particle level. Moreover, the colloidal solitons form despite of negligible inertia in the particle motion, i.e. in the limit of fully overdamped Brownian dynamics. They reach a speed much higher than that of the driving potential and even move against the driving direction.
Particle transport against an external bias is a rather unusual effect. It occurs under special circumstances, for example, in the presence of strong confinement22,23,24, multi-particle interactions25 or when the imposed potential is modulated in time26,27. For the soliton transport reported here, particle motion is still in the driving direction, but the spontaneously forming collective excitations in form of localized clusters propagate fast against the external bias.
By extending a theoretical model based on a recent prediction11, we explain the observed speed, growth, and direction of motion of the solitons. We further strengthen our results by complementing them with Brownian dynamics simulations carried out using experimental parameters.
The underlying mechanism of soliton formation is generic and can be extended to other driven periodic systems under overcrowding conditions, such as particles flowing through porous media28,29, vortices in high-Tc superconducting heterostructures30,31,32, frictional atomic sliding33,34, driven skyrmions35,36 ultracold atoms37, photonic38,39,40, and active matter41,42,43,44,45 systems.
Results
Soliton observation
Our setup, schematically shown in Fig. 1a, is made of a closed fluidic cell filled with a water dispersion of spherical polystyrene particles having a diameter of Ïâ=â4âμm46. We confine N of these particles in a rotating one-dimensional periodic potential which is created by an infrared laser beam rapidly steered across Mâ=â27 equispaced positions along a circle of radius R. More details can be found in the Methods section. Once the laser is scanned trough the array of optical traps, it generates an effective periodic potential along the azimuthal direction with wavelength λâ=â2ÏR/M, which rotates at a constant angular velocity Ï, Fig. 1a.
The generated optical potential strongly confines the particles in the radial direction. Indeed, particle displacements along this direction have a Gaussian distribution with a mean at the ring radius Râ=â20âμm and a standard deviation smaller than 0.01âRâ=â0.2âμm, i.e., they are negligible compared to the particle motion along the ring. Thus, we consider particle motion to take place in a traveling-wave potential
where xâ=âRÏ with Ï the azimuthal angle along the ring. A single trapped colloidal particle is dragged clockwise by this translating wave. For low Ï, it moves at a constant tangential speed ÏR. By increasing Ï, the mean speed decreases because the particle looses its synchronized phase with the moving potential due to the viscous drag47.
The synchronization becomes increasingly restored with larger N due to an effective potential barrier enhancement by hydrodynamic interactions48. For N = M the particle motion is fully synchronized with a mean particle velocity ÏR, as demonstrated in the Supplementary Movie 1.
A different type of collective excitation occurs for Nâ>âM. In this situation of overcrowding, after a short transitory period, the particles generate localized clusters that stably propagate along the ring without dispersion, Fig. 1b (Supplementary Movie 2). This cluster formation is not obvious in light of negligible attractive interactions between the colloidal particles in our experiments.
Let us consider the case of one extra particle, Nâ=âMâ+â1, which generates a double-occupied trap. The excess particle displaces the colloids in neighboring traps from the preferred position close to the potential minimum, forming an extended defect, which appears as an almost compact cluster composed of particles that are nearly in contact. The defect propagates as a dispersion-free solitary wave, or soliton, along the optical ring.
Surprisingly, it moves counter-clockwise along the ring, i.e., backwards against the external driving, and the speed of propagation is much faster than that of the rotating optical traps. The backward propagation occurs independent of the sense of rotation of the optical landscape, which means it is clockwise if the traps are rotated counter-clockwise and vice versa. It is not caused by the viscous drag only. In the case of negligible potential barriers, the viscous drag cannot lead to backward motion, but at most to a vanishing mean displacement of the particles along the ring.
Soliton stability and propagation
To understand the backward movement of a cluster, we must consider a wavelike cooperative movement, where the particles forming the cluster change. Specifically, this change is occurring due to particles detaching from the cluster at its back end and attaching at its front end, as shown in Fig. 1c. Back and front end refer to the direction of the cluster motion. However, why do particle stay together without attractive interactions and why does their cooperative movement appear as a stable soliton?
To answer these questions, let us view the motion in a frame corotating with the optical traps, for which we show particle trajectories in Fig. 1d. If the particle number N is equal to the number M of traps, each particle in this frame stays at the bottom of a potential well. For Nâ=âMâ+â1, coherent particle movements between the wells are present due to the soliton propagation, see the parts of the particle trajectories marked in red.
In the corotating frame, the particles are driven by a mean flow of the surrounding fluid that tries to move them clockwise, corresponding to the positive x-direction. The effective force of the flow-driving follows from a coordinate transformation to the comoving frame: \({x}^{{\prime} }=x+\omega Rt\), yielding \({{{{{{{\rm{d}}}}}}}}{x}^{{\prime} }/{{{{{{{\rm{d}}}}}}}}t={{{{{{{\rm{d}}}}}}}}x/{{{{{{{\rm{d}}}}}}}}t+\omega R\). This amounts to a constant drag force Fâ=âÏR/μ acting on the particles in the corotating frame, which locally tilts the periodic potential. For the parameters used in our experiments, this tilting is always undercritical and the effective barrier between traps much larger than the thermal energy kBT. This means that a single particle could hardly move from one optical trap to the next due to thermal excitation. In the many-particle system with filling factor one (Nâ=âM), where the system is perfectly ordered, the drag force F does not generate sustainable motion, see Fig. 1d.
In the overcrowded system, propagating clusters of particles can occur, which are not formed by attractive interactions but are stabilized by the external forces \({F}^{{{{{{{{\rm{ext}}}}}}}}}({x}^{{\prime} })=\omega R/\mu -\partial U({x}^{{\prime} })/\partial {x}^{{\prime} }\), where \(U({x}^{{\prime} })=({U}_{0}/2)\cos (2\pi {x}^{{\prime} }/\lambda )\) is the potential in the corotating frame. The condition on the external forces to hold n particles together in an n-cluster is11
Here, \({F}_{j}^{{{{{{{{\rm{ext}}}}}}}}}={F}^{{{{{{{{\rm{ext}}}}}}}}}({x}_{j}^{{\prime} })\) is the external force on the jth particle in the cluster. If inequalities (2) hold, a fragmentation of the cluster is impossible. This is because for a fragmentation to occur into a left subcluster of the first i particles and a right subcluster of the (nâââi) remaining particles, the average force \(\mathop{\sum }\nolimits_{j=1}^{i}{F}_{j}^{{{{{{{{\rm{ext}}}}}}}}}/i\) on the left subcluster would have to be smaller than the average force \(\mathop{\sum }\nolimits_{j=i+1}^{n}{F}_{j}^{{{{{{{{\rm{ext}}}}}}}}}/(n-i)\) on the right subcluster.
For the particles to stay together in the cluster during its motion, the inequalities (2) must be obeyed for all points in some interval. If this interval would span a full period λ, then the cluster would move without changing its size. However, the soliton dynamics are more complex, as the size of a cluster changes within λ. This leads to different types of solitons.
Soliton types
As the potential barrier U0 is much larger than kBT, formation and propagation of solitons can be understood by considering the limit of vanishing noise. Theoretically, the solitons consist of periodically repeating movements of clusters with different size in this limit11,49, for which the equations of motions in the corotating frame are
Introducing scaled dimensionless coordinates and time, \({x}_{i}^{{\prime} }\to {y}_{i}={x}_{i}^{{\prime} }/\lambda\), tâââλ2t/(ÏμU0), and a dimensionless driving force fâ=âλÏR/(ÏμU0)â=âFλ/(ÏU0), these equations take the form
Solving Eq. (4) subject to the force conditions (2) and an initial condition with one double-occupied potential well, we find that after a transient time, periodic motions of two soliton types appear as limit cycles: an A type soliton given by two subintervals of the movements of an n- and (nâ+â1)-cluster during one period [n-(nâ+â1)-soliton], and a B type soliton given by four subintervals of cluster movements [n-(nâ+â1)-(nâ+â2)-(nâ+â1)-soliton].
The decrease and increase of a cluster size in the sequences is by detachments and attachments of a single particle to the cluster. Exact expressions for the position of and times between attachment and detachment events are derived in Supplementary Discussion 1. Knowing these quantities, we can in particular calculate mean sizes and mean angular velocities of the solitons.
Type B solitons occur only in narrow regimes of particles sizes and driving forces. Typical solitons are of type A. Within one period of the motion of an A type soliton of core size n, a single particle detaches at the back end of an (nâ+â1)-cluster and attaches at the front end of an n-cluster. Back and front end refer to the direction of cluster motion in the corotating frame, i.e., to the direction of positive x for fâ>â0.
Figure 2a shows type A solitons observed in the experiments. For Ïâ=â0.67λ, sequences of 3- and 2-cluster movements occur, corresponding to a 3-2-soliton. For Ïâ=â0.76λ, we see the occurrence of a 4-3-soliton. The clusters in these solitons are not formed by particles that are exactly in contact, because the noise in the experiments is weak but not vanishing. Due to thermal fluctuations, particles cannot stay in contact during their motion. In our analysis of both experiments and simulations at finite noise, we identify a cluster as a sequence of neighboring particles, where the empty gap between two neighboring particles is smaller than a cutoff-distance 5âÃâ10â3λ. The number n of particles in the corresponding sequence gives the cluster size. For defining the soliton position, we select the pair of particles in the cluster whose positions are in the same optical trap, as shown in Fig. 1a. The soliton coordinate is taken to be the center of mass position [angle coordinate Ï(t)] of the two particles belonging to that pair.
Figure 2b shows the time intervals of alternating n- and (nâ+â1)-cluster motion. The duration of two consecutive time intervals is close to the average time period λ/vsolâ=âλ/ÏsolR, reflecting the underlying periodicity of the soliton motion. In the presence of noise, the motion is no longer perfectly periodic in time, but the power spectrum of the cluster position should still reflect the underlying periodicity by exhibiting a peak at a frequency vsol/λ. Indeed, we find this peak in the spectrum in both experiments [top image in Fig. 2d] and numerical simulations [bottom image in Fig. 2d]. To filter out the trivial peak at frequency vsol/L, we have used the relative position of the clusters to the next potential minimum in the calculation of the spectra.
In Fig. 2c we show the fraction of time a soliton stays in the n-cluster state in dependence of the particle diameter. These results have been calculated analytically in the zero-noise limit, see Supplementary Discussion 1. They demonstrate in particular that the clusters of the solitons become larger with increasing Ï, in agreement with the experimental observations.
Soliton size and speed
In Fig. 3a we show that the mean cluster size ãnã increases rapidly with Ï/λ, as the perturbation induced by a double-occupied trap spreads over longer distances. In the experiments (empty symbols), the wavelength λ was varied according to λâ=â2ÏR/M with Mâ=â21-27. We can estimate ãnã by the following reasoning: An n-cluster covers a space nÏ. The remaining (Nââân) particles not belonging to the n-cluster are close to potential minima and thus distributed over the length (Nââân)λ. The two lengths should fill the ring of length Mλâ=â(Nâââ1)λ, yielding (Nââân)λâ+ânÏâââ(Nâââ1)λ, i.e., nÏâââ(nâââ1)λ. As this reasoning can be applied to any cluster in a soliton mode, we obtain
where α is a prefactor of order unity. In Supplementary Note 4 we show that this approximate expression indeed captures the result of an exact calculation of ãnã in the zero-noise limit. The experimental results for ãnã in Fig. 3a can be well fitted to Eq. (5) with αâ=â0.83.
Figure 3b shows the normalized mean angular soliton speed Ïsol/Ï as a function of the scaled particle diameter Ï/λ, where we varied the potential wavelength λ in the experiments, as in Fig. 3b. The observed speed raises nonlinearly with Ï/λ, reaching a maximum value of vsolâ=â34.7âμmâsâ1 for λâ=â4.65âμm, which is about five times higher than the speed of the optical traps moving in the opposite direction (ÏRâ=â7.2âμmâsâ1). This means that the soliton velocity relative to the optical traps is almost six times higher.
The soliton mean velocity vsol can be estimated by a scaling argument in the corotating frame. After the detachment and attachment of a particle, the mean distance moved by the cluster is one wavelength λ. The time needed for an attachment and detachment is of the order of the mean distance (λâââÏ) between two particles divided by the velocity ÏR. Accordingly, vsol should be proportional to λÏR/(λâââÏ), i.e. vsolâ=âβÏR/(λâââÏ)â=âβÏRâãnã, where β is a constant. This scaling argument is corroborated by an exact calculation of soliton velocities in the zero-noise limit outlined in Supplementary Discussion 1. In the laboratory frame, this gives
Figure 3b shows that the mean soliton velocity in the experiment can be well described by Eq. (6) with \({\beta }_{\exp }=1.02\). The simulated data can be well fitted with βthâ=â1.31. This value is very close to the one obtained when fitting Eq. (6) to the the exact analytical calculation in the zero-noise limit, see Supplementary Discussion 1. The difference between \({\beta }_{\exp }\) and βth can be explained by the fact that an ideal traveling sinusoidal wave is modeled in the simulations, while in the experiments there are always some small deviations from the ideal behavior, e.g., deviations from the exact sinusoidal form of the potential and from a purely one-dimensional particle motion. Due to these imperfections, we have rescaled the simulated soliton velocities in the corotating frame by the factor \({\beta }_{\exp }/{\beta }_{{{{{{{{\rm{th}}}}}}}}}\). The correspondingly rescaled simulated data are shown in Fig. 3a.
Many interacting colloidal solitons
We find that the number of solitons present in our system is equal to the overcrowding, defined as the difference NâââM. Figure 4a demonstrates the increase of the number of solitons with the overcrowding for two particle diameters. Our observations suggest that solitons tend to repel each other and propagate at a well-defined mean distance. To quantify this effect, we analyzed distributions of soliton distances.
For two solitons, the distribution Ï(ÎÏ) of distances ÎÏ(t)â=ââ£Ï2(t)âââÏ1(t)⣠between the positions Ï1(t) and Ï2(t) of two solitons is shown in Fig. 4 for experiment (b, c) and simulation (d, e). In both cases, Ï(ÎÏ) exhibits several peaks with a Gaussian-like envelope centered around a mean \(\overline{{{\Delta }}\varphi }\simeq \pi\). Such fine structure appears due to the fact that the solitons have preferential positions close to the potential minima, while other positions are less likely. Accordingly, the distance distribution shows peaks separated by 2Ï/M. The localized envelope points to an effective repulsive soliton-soliton interaction, which tends to keep the soliton positions at a maximal distance Ï. Note that the difference between the experimental (Fig. 4b, c) and simulation data (Fig. 4d, e) is only in the widths of the angle distributions, but not in the mean value, which demonstrates the presence of an effective solitonâsoliton interaction.
For steady states characterized by three propagating solitons, their position can be ordered at each time t such that Ï1(t)â<âÏ2(t)â<âÏ3(t). Thus we show in Fig. 4c, e the distributions Ï(ÎÏ12) and Ï(ÎÏ13) of the distances ÎÏ12â=âÏ2(t)âââÏ1(t) and ÎÏ13â=âÏ3(t)âââÏ1(t). Here the envelopes of Ï(ÎÏ12) and Ï(ÎÏ13) are centered close to 2Ï/3 and 4Ï/3, respectively.
Discussion
We have experimentally observed solitary cluster waves in a highly crowded system of driven Brownian particles along a periodic potential. The solitons are robust, their number can be controlled by the overcrowding, and they propagate against the rotation direction at a speed up to six times higher than the angular velocity of the rotating wells. Viewed in a frame comoving with the potential wells, the experiment realizes particle motion through a periodic potential driven by a constant drag force. The solitons then move in the force direction and provide a formidable way to transport matter at the microscale in confined space.
Cluster-mediated transport realized in our work may appear similar to the dynamics of kinks and anti-kinks predicted by the Frenkel-Kontorova (FK) model50, as, for example, seen in nanotribology33,34 and diffusion of crowdions51,52 and voidions53,54. In recent experimental realizations on the colloidal length scale, kinks were produced by sliding electrostatically repulsive microspheres across a static, two-dimensional potential55,56. The mechanism behind the kink formation in the FK model appears, however, to be fundamentally different from the mechanism of soliton formation reported here, which is given by the conditions in Eq. (2) between the external forces. Due to these conditions, particles can keep together during cluster propagation. This mechanism leads to cluster waves that, once excited, are robust and persist during all experimental time.
Cluster-mediated transport is widespread in many artificial and biological systems including, for example, microfluidic channels or vein networks. The rich physics unveiled in an overcrowded single-file system driven across a periodic potential opens different perspectives for further research.
Methods
Experimental system
Our colloidal suspension consists of spherical polystyrene particles having a diameter of Ïâ=â4âμm (CML, Molecular Probes, 4% w/v in water) that are dispersed in deionized water (MilliQ). This suspension is confined within a fluidic cell, Fig. 1a, which is composed of two coverslips separated byâââ100âμm and placed on the stage of a custom-built optical microscope. The particles are left sediment close to the bottom of the cell due to density mismatch, and float there showing a small diffusion coefficient Dâ=â0.1295âμm2âsâ1, as measured from the mean squared displacement in absence of the optical potential. The system dynamics are visualized by using a Nikon 40âÃâmicroscope objective (plan Apo) illuminated by a light emitting diode, while video recording is performed at 60âHz with a complementary metal oxide camera (Ximea MQ003MG-CM).
We create a rotating one-dimensional periodic potential by passing an infrared laser beam (Manlight ML5-CW-P/TKS-OTS) through a pair of acousto-optic deflectors (AA Optoelectronics DTSXY-400-1064). The beam is characterized by a wavelength 1064âμm and power Pâ=â3âW. The AODs have an input frequency in the range from 60 to 90âMHz. It is produced via a two-channel radio frequency wave generator (DDSPA2X-D431b-34) which is addressed by a digital output card (National Instruments cDAQ NI-9403) with a refresh frequency of 150âkHz. More technical details on the experimental system can be found in a related work47, which investigates the dynamics of a number of particles in the low filling regime, Nâ<âM,46.
Numerical simulations
In the experiments, the particles perform a Brownian motion in the time-dependent optical potential generated by the rotating laser beam. As shown in ref. 47, these dynamics can be described by that of hard-spheres with diameter Ï, where the center of mass position ri of each particle i moves according to the Langevin equation
Here Uopt(r,ât) is the time-dependent optical potential, μâ=âD/kBT is the particle mobility, and ζi(t) are Gaussian white noise processes with ãζi(t)ãâ=â0 and \(\langle {\zeta }_{i\alpha }(t){\zeta }_{j\beta }({t}^{{\prime} })\rangle=2D{\delta }_{ij}{\delta }_{\alpha \beta }\delta (t-{t}^{{\prime} })\). Because the particle confinement in the radial direction is very strong and the radius of curvature is much greater than the mean distance between the colloidal particles, the motion can be restricted to one dimension in a traveling-wave potential \(U(x,t)=({U}_{0}/2)\cos (2\pi x/\lambda+M\omega t)\) with xâ=âRÏ, where Ï is the azimuthal angle47,57,58. Note that Ï is the angular velocity of the rotation of the optical traps, which implies that the traveling wave has the frequency MÏ for M traps. The impact of hydrodynamic interactions is effectively tackled by the recently discovered potential barrier enhancement48 and the condition â£riââârjâ£ââ¥âÏ implied by the hard-sphere interaction is treated by the method proposed in ref. 59 and additionally by a recently developed cluster algorithm60. The two simulations methods gave the same results.
Data availability
The data that support the findings of this study are available from the corresponding author upon request. Parts of the originally submitted manuscript have been incorporated in a doctoral thesis46 available at http://hdl.handle.net/10803/688857. Source data are provided with this paper.
Code availability
All computer codes are available from the corresponding authors upon request.
References
Reimann, P. Brownian motors: noisy transport far from equilibrium. Phys. Rep. 361, 57 (2002).
Hänggi, P. & Marchesoni, F. Artificial Brownian motors: controlling transport on the nanoscale. Rev. Mod. Phys. 81, 387 (2009).
Reichhardt, C. & Reichhardt, C. J. O. Depinning and nonequilibrium dynamic phases of particle assemblies driven over random and ordered substrates: a review. Rep. Prog. Phys. 80, 026501 (2016).
Reichhardt, C. & Nori, F. Phase locking, devilâs staircases, Farey trees, and Arnold tongues in driven vortex lattices with periodic pinning. Phys. Rev. Lett. 82, 414 (1999).
Korda, P. T., Taylor, M. B. & Grier, D. G. Kinetically lockedin colloidal transport in an array of optical tweezers. Phys. Rev. Lett. 89, 128301 (2002).
Bohlein, T. & Bechinger, C. Experimental observation of directional locking and dynamical ordering of colloidal monolayers driven across quasiperiodic substrates. Phys. Rev. Lett. 109, 058301 (2012).
Tierno, P. & Fischer, T. M. Excluded volume causes integer and fractional plateaus in colloidal ratchet currents. Phys. Rev. Lett. 112, 048302 (2014).
Juniper, M. P., Straube, A. V., Besseling, R., Aarts, D. G. A. L. & Dullens, R. P. Microscopic dynamics of synchronization in driven colloids. Nat. Commun. 6, 7187 (2015).
Reichhardt, C. & Reichhardt, C. J. O. Random organization and plastic depinning. Phys. Rev. Lett. 103, 168301 (2009).
Tierno, P. Depinning and collective dynamics of magnetically driven colloidal monolayers. Phys. Rev. Lett. 109, 198304 (2012).
Antonov, A. P., Ryabov, A. & Maass, P. Solitons in overdamped Brownian dynamics. Phys. Rev. Lett. 129, 080601 (2022).
Drazin, P. G. & Johnson, R. S. Solitons: An Introduction (Cambridge University Press, Cambridge, 1993).
Yefsah, T. et al. Heavy solitons in a fermionic superfluid. Nature 499, 426 (2013).
Strecker, K. E., Partridge, G. B., Truscott, A. G. & Hulet, R. G. Formation and propagation of matter-wave soliton trains. Nature 417, 150 (2002).
Nguyen, J. H. V., Luo, D. & Hulet, R. G. Formation of matter-wave soliton trains by modulational instability. Science 356, 422 (2017).
Sharon, E., Cohen, G. & Fineberg, J. Propagating solitary waves along a rapidly moving crack front. Nature 410, 68 (2001).
Chen, B. G., Upadhyaya, N. & Vitelli, V. Nonlinear conduction via solitons in a topological mechanical insulator. Proc. Natl Acad. Sci. USA 111, 13004 (2014).
Deng, B., Wang, P., He, Q., Tournat, V. & Bertoldi, K. Metamaterials with amplitude gaps for elastic solitons. Nat. Commun. 9, 3410 (2018).
Zhang, Y., Li, B., Zheng, Q. S., Genin, G. M. & Chen, C. Q. Programmable and robust static topological solitons in mechanical metamaterials. Nat. Commun. 10, 5605 (2019).
Costa, A. et al. Soliton turbulence in shallow water ocean surface waves. Phys. Rev. Lett. 113, 108501 (2014).
Ablowitz, M. J. Nonlinear Dispersive Waves: Asymptotic Analysis and Solitons, Cambridge Texts in Applied Mathematics (Cambridge University Press, Cambridge, 2011).
Eichhorn, R., Reimann, P. & Hänggi, P. Brownian motion exhibiting absolute negative mobility. Phys. Rev. Lett. 88, 190601 (2002).
Matthias, S. & Müller, F. Asymmetric pores in a silicon membrane acting as massively parallel Brownian ratchets. Nature 424, 53 (2003).
Ros, A. et al. Absolute negative particle mobility. Nature 436, 928 (2005).
McDermott, D., Reichhardt, C. J. O. & Reichhardt, C. Collective ratchet effects and reversals for active matter particles on quasi-one-dimensional asymmetric substrates. Soft Matter 12, 8606 (2016).
Lee, S.-H., Ladavac, K., Polin, M. & Grier, D. G. Observation of flux reversal in a symmetric optical thermal ratchet. Phys. Rev. Lett. 94, 110601 (2005).
Schwemmer, C., Fringes, S., Duerig, U., Ryu, Y. K. & Knoll, A. W. Experimental observation of current reversal in a rocking Brownian motor. Phys. Rev. Lett. 121, 104102 (2018).
Gerber, G., Rodts, S., Aimedieu, P., Faure, P. & Coussot, P. Particle-size-exclusion clogging regimes in porous media. Phys. Rev. Lett. 120, 148001 (2018).
Bizmark, N., Schneider, J., Priestley, R. D. & Datta, S. S. Multiscale dynamics of colloidal deposition and erosion in porous media. Sci. Adv. 6, eabc2530 (2020).
Blatter, G., Feigelâman, M. V., Geshkenbein, V. B., Larkin, A. I. & Vinokur, V. M. Vortices in high-temperature superconductors. Rev. Mod. Phys. 66, 1125 (1994).
Villegas, J. E. et al. A superconducting reversible rectifier that controls the motion of magnetic flux quanta. Science 302, 1188 (2003).
Wu, H. et al. The field-free Josephson diode in a van der Waals heterostructure. Nature 604, 653 (2022).
Vanossi, A., Manini, N., Urbakh, M., Zapperi, S. & Tosatti, E. Colloquium: Modeling friction: from nanoscale to mesoscale. Rev. Mod. Phys. 85, 529 (2013).
Hod, O., Meyer, E., Zheng, Q. & Urbakh, M. Structural superlubricity and ultralow friction across the length scales. Nature 563, 485 (2018).
Mühlbauer, S. et al. Skyrmion lattice in a chiral magnet. Science 323, 915 (2009).
Reichhardt, C., Reichhardt, C. J. O. & MiloseviÄ, M. V. Statics and dynamics of skyrmions interacting with disorder and nanostructures. Rev. Mod. Phys. 94, 035005 (2022).
Bloch, I. Quantum coherence and entanglement with ultracold atoms in optical lattices. Nature 453, 1016 (2008).
Wang, P. et al. Localization and delocalization of light in photonic Moirâe lattice. Nature 577, 42 (2020).
Christodoulides, D. N., Lederer, F. & Silberberg, Y. Discretizing light behaviour in linear and nonlinear waveguide lattices. Nature 424, 817 (2003).
Jürgensen, M., Mukherjee, S., Jörg, C. & Rechtsman, M. C. Quantized fractional Thouless pumping of solitons. Nat. Phys. 19, 420â426 (2023).
Bechinger, C. et al. Active particles in complex and crowded environments. Rev. Mod. Phys. 88, 045006 (2016).
Reichhardt, C. J. O. & Reichhardt, C. Ratchet effects in active matter systems. Ann. Rev. Cond. Matt. Phys. 8, 51 (2017).
Bhattacharjee, T. & Datta, S. S. Bacterial hopping and trapping in porous media. Nat. Commun. 10, 2075 (2019).
Meacock, O. J., Doostmohammadi, A., Foster, K. R., Yeomans, J. M. & Durham, W. M. Bacteria solve the problem of crowding by moving slowly. Nat. Phys. 17, 205 (2021).
Dehkharghani, A., Waisbord, N. & Guasto, J. S. Selftransport of swimming bacteria is impaired by porous microstructure. Commun. Phys. 6, 18 (2023).
Cereceda-López, E. Non-Equilibrium Dynamics of Driven and Confined Colloidal Systems (Phd Thesis, University of Barcelona, 2023).
Lips, D. et al. Hydrodynamic interactions hinder transport of flow-driven colloidal particles. Soft Matter 18, 8983 (2022).
Cereceda-López, E. et al. Hydrodynamic interactions can induce jamming in flow-driven systems. Phys. Rev. Lett. 127, 214501 (2021).
Antonov, A. P., Vorác, D., Ryabov, A. & Maass, P. Collective excitations in jammed states: ultrafast defect propagation and finite-size scaling. New J. Phys. 24, 093020 (2022).
Braun, O. M. & Kivshar, Y. S. Nonlinear dynamics of the Frenkel-Kontorova model. Phys. Rep. 306, 1 (1998).
Paneth, H. R. The mechanism of self-diffusion in alkali metals. Phys. Rev. 80, 708 (1950).
Derlet, P. M., Nguyen-Manh, D. & Dudarev, S. L. Multiscale modeling of crowdion and vacancy defects in bodycentered- cubic transition metals. Phys. Rev. B 76, 054107 (2007).
Matsukawa, Y. & Zinkle, S. J. One-dimensional fast migration of vacancy clusters in metals. Science 318, 959 (2007).
van der Meer, B., van Damme, R., Dijkstra, M., Smallenburg, F. & Filion, L. Revealing a vacancy analog of the crowdion interstitial in simple cubic crystals. Phys. Rev. Lett. 121, 258001 (2018).
Bohlein, T., Mikhael, J. & Bechinger, C. Observation of kinks and antikinks in colloidal monolayers driven across ordered surfaces. Nat. Mater. 11, 126 (2012).
Brazda, T., Julya, C. & Bechinger, C. Experimental observation of shapiro-steps in colloidal monolayers driven across time-dependent substrate potentials. Soft Matter 13, 4024 (2017).
Juniper, M. P. N., Straube, A. V., Aarts, D. G. A. L. & Dullens, R. P. A. Colloidal particles driven across periodic optical-potential-energy landscapes. Phys. Rev. E 93, 012608 (2016).
Juniper, M. P. N., Besseling, R., Aarts, D. G. A. L. & Dullens, R. P. A. Acousto-optically generated potential energy landscapes: potential mapping using colloids under flow. Opt. Express 20, 28707 (2012).
Scala, A. Event-driven Langevin simulations of hard spheres. Phys. Rev. E 86, 026709 (2012).
Antonov, A. P., Schweers, S., Ryabov, A. & Maass, P. Brownian dynamics simulations of hard rods in external fields and with contact interactions. Phys. Rev. E 106, 054606 (2022).
Acknowledgements
This project has received funding from the European Research Council (ERC) under the European Unionâs Horizon 2020 research and innovation programme (grant agreement no. 811234). A.A., A.R., and P.M. gratefully acknowledge financial support by the Czech Science Foundation (Project No. 20-24748J) and the Deutsche Forschungsgemeinschaft (Project No. 432123484). P.T. acknowledges support from the Ministerio de Ciencia e Innovació (Project No. PID2022-137713NB-C22), the Agència de Gestió dâAjuts Universitaris i de Recerca (Project No. 2021 SGR 00450) and the Generalitat de Catalunya under Program âICREA Acadèmiaâ.
Author information
Authors and Affiliations
Contributions
E.C. performed the experiments and analyzed the experimental data. A.A. performed numerical simulations, analyzed experimental data, and derived analytical formulas presented in the supplementary information. A.R., P.M., and P.T. supervised the work. All authors discussed the results and commented on the manuscript at all stages.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cereceda-López, E., Antonov, A.P., Ryabov, A. et al. Overcrowding induces fast colloidal solitons in a slowly rotating potential landscape. Nat Commun 14, 6448 (2023). https://doi.org/10.1038/s41467-023-41989-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-41989-x