1
|
Toivonen E, Kaipainen J, Molkkari M, Keski-Rahkonen J, Klages R, Räsänen E. Anomalous diffusion in the square soft Lorentz gas. Phys Rev E 2025; 111:014216. [PMID: 39972858 DOI: 10.1103/physreve.111.014216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2024] [Accepted: 12/24/2024] [Indexed: 02/21/2025]
Abstract
We demonstrate and analyze anomalous diffusion properties of point-like particles in a two-dimensional system with circular scatterers arranged in a square lattice and governed by smooth potentials, referred to as the square soft Lorentz gas. Our numerical simulations reveal a rich interplay of normal and anomalous diffusion depending on the system parameters. To describe diffusion in normal regimes, we develop a unit cell hopping model that, in the single-hop limit, recovers the Machta-Zwanzig approximation and converges toward the numerical diffusion coefficient as the number of hops increases. Anomalous diffusion is characterized by quasiballistic orbits forming Kolmogorov-Arnold-Moser islands in phase space, alongside a complex tongue structure in parameter space defined by the interscatterer distance and potential softness. The distributions of the particle displacement vector show notable similarities to both analytical and numerical results for a hard-wall square Lorentz gas, exhibiting Gaussian behavior in normal diffusion and long tails due to quasiballistic orbits in anomalous regimes. Our work thus provides a catalog of key dynamical system properties that characterize the intricate changes in diffusion when transitioning from hard billiards to soft potentials.
Collapse
Affiliation(s)
- Esko Toivonen
- Tampere University, Computational Physics Laboratory, P.O. Box 600, FI-33014 Tampere, Finland
| | - Joni Kaipainen
- Tampere University, Computational Physics Laboratory, P.O. Box 600, FI-33014 Tampere, Finland
- VTT Technical Research Centre of Finland, Integrated Computational Materials Engineering Group, Ltd, FI-02044 Espoo, Finland
| | - Matti Molkkari
- Tampere University, Computational Physics Laboratory, P.O. Box 600, FI-33014 Tampere, Finland
| | - Joonas Keski-Rahkonen
- Tampere University, Computational Physics Laboratory, P.O. Box 600, FI-33014 Tampere, Finland
- Harvard University, Department of Physics, Cambridge, Massachusetts 02138, USA
| | - Rainer Klages
- Queen Mary University of London, Centre for Complex Systems, School of Mathematical Sciences, Mile End Road, London E1 4NS, United Kingdom
- London Mathematical Laboratory, 8 Margravine Gardens, London W6 8RH, United Kingdom
| | - Esa Räsänen
- Tampere University, Computational Physics Laboratory, P.O. Box 600, FI-33014 Tampere, Finland
| |
Collapse
|
2
|
Bystrik YS, Denisov S. Asymptotic densities of planar Lévy walks: A nonisotropic case. Phys Rev E 2021; 104:064131. [PMID: 35030855 DOI: 10.1103/physreve.104.064131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2021] [Accepted: 12/06/2021] [Indexed: 06/14/2023]
Abstract
Lévy walks are a particular type of continuous-time random walks which results in a super-diffusive spreading of an initially localized packet. The original one-dimensional model has a simple schematization that is based on starting a new unidirectional motion event either in the positive or in the negative direction. We consider two-dimensional generalization of Lévy walks in the form of the so-called XY model. It describes a particle moving with a constant velocity along one of the four basic directions and randomly switching between them when starting a new motion event. We address the ballistic regime and derive solutions for the asymptotic density profiles. The solutions have a form of first-order integrals which can be evaluated numerically. For specific values of parameters we derive an exact expression. The analytic results are in agreement with the results of finite-time numerical samplings.
Collapse
Affiliation(s)
- Yu S Bystrik
- Institute of Applied Physics, National Academy of Sciences of Ukraine, Petropavlivska Street 58, 40000 Sumy, Ukraine
| | - S Denisov
- Department of Computer Science, Oslo Metropolitan University, N-0130 Oslo, Norway and NordSTAR-Nordic Center for Sustainable and Trustworthy AI Research, Pilestredet 52, Oslo N-0166, Norway
| |
Collapse
|
3
|
Propagation of waves in high Brillouin zones: Chaotic branched flow and stable superwires. Proc Natl Acad Sci U S A 2021; 118:2110285118. [PMID: 34580232 DOI: 10.1073/pnas.2110285118] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/16/2021] [Indexed: 11/18/2022] Open
Abstract
We report unexpected classical and quantum dynamics of a wave propagating in a periodic potential in high Brillouin zones. Branched flow appears at wavelengths shorter than the typical length scale of the ordered periodic structure and for energies above the potential barrier. The strongest branches remain stable indefinitely and may create linear dynamical channels, wherein waves are not confined directly by potential walls as electrons in ordinary wires but rather, indirectly and more subtly by dynamical stability. We term these superwires since they are associated with a superlattice.
Collapse
|
4
|
Hu SY, Chu JJ, Shelley MJ, Zhang J. Lévy Walks and Path Chaos in the Dispersal of Elongated Structures Moving across Cellular Vortical Flows. PHYSICAL REVIEW LETTERS 2021; 127:074503. [PMID: 34459633 DOI: 10.1103/physrevlett.127.074503] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Revised: 03/19/2021] [Accepted: 07/20/2021] [Indexed: 06/13/2023]
Abstract
In cellular vortical flows, namely arrays of counterrotating vortices, short but flexible filaments can show simple random walks through their stretch-coil interactions with flow stagnation points. Here, we study the dynamics of semirigid filaments long enough to broadly sample the vortical field. Using simulation, we find a surprising variety of long-time transport behavior-random walks, ballistic transport, and trapping-depending upon the filament's relative length and effective flexibility. Moreover, we find that filaments execute Lévy walks whose diffusion exponents generally decrease with increasing filament length, until transitioning to Brownian walks. Lyapunov exponents likewise increase with length. Even completely rigid filaments, whose dynamics is finite dimensional, show a surprising variety of transport states and chaos. Fast filament dispersal is related to an underlying geometry of "conveyor belts." Evidence for these various transport states is found in experiments using arrays of counterrotating rollers, immersed in a fluid and transporting a flexible ribbon.
Collapse
Affiliation(s)
- Shi-Yuan Hu
- Applied Mathematics Lab, Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- Department of Physics, New York University, New York, New York 10003, USA
| | - Jun-Jun Chu
- School of Physics Science and Engineering, Tongji University, Shanghai 200092, China
| | - Michael J Shelley
- Applied Mathematics Lab, Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- Center for Computational Biology, Flatiron Institute, New York, New York 10010, USA
| | - Jun Zhang
- Applied Mathematics Lab, Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- Department of Physics, New York University, New York, New York 10003, USA
- NYU-ECNU Institute of Physics at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
5
|
Barkai E, Burov S. Packets of Diffusing Particles Exhibit Universal Exponential Tails. PHYSICAL REVIEW LETTERS 2020; 124:060603. [PMID: 32109131 DOI: 10.1103/physrevlett.124.060603] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Revised: 10/23/2019] [Accepted: 11/20/2019] [Indexed: 06/10/2023]
Abstract
Brownian motion is a Gaussian process described by the central limit theorem. However, exponential decays of the positional probability density function P(X,t) of packets of spreading random walkers, were observed in numerous situations that include glasses, live cells, and bacteria suspensions. We show that such exponential behavior is generally valid in a large class of problems of transport in random media. By extending the large deviations approach for a continuous time random walk, we uncover a general universal behavior for the decay of the density. It is found that fluctuations in the number of steps of the random walker, performed at finite time, lead to exponential decay (with logarithmic corrections) of P(X,t). This universal behavior also holds for short times, a fact that makes experimental observations readily achievable.
Collapse
Affiliation(s)
- Eli Barkai
- Physics Department, Bar-Ilan University, Ramat Gan 5290002, Israel
| | - Stanislav Burov
- Physics Department, Bar-Ilan University, Ramat Gan 5290002, Israel
| |
Collapse
|
6
|
Zarfaty L, Peletskyi A, Barkai E, Denisov S. Infinite horizon billiards: Transport at the border between Gauss and Lévy universality classes. Phys Rev E 2019; 100:042140. [PMID: 31770907 DOI: 10.1103/physreve.100.042140] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2019] [Indexed: 06/10/2023]
Abstract
We consider transport in two billiard models, the infinite horizon Lorentz gas and the stadium channel, presenting analytical results for the spreading packet of particles. We first obtain the cumulative distribution function of traveling times between collisions, which exhibits nonanalytical behavior. Using a renewal assumption and the Lévy walk model, we obtain the particles' probability density. For the Lorentz gas, it shows a distinguished difference when compared with the known Gaussian propagator, as the latter is valid only for extremely long times. In particular, we show plumes of particles spreading along the infinite corridors, creating power-law tails of the density. We demonstrate the slow convergence rate via summation of independent and identically distributed random variables on the border between Lévy and Gauss laws. The renewal assumption works well for the Lorentz gas with intermediate-size scattering centers, but fails for the stadium channel due to strong temporal correlations. Our analytical results are supported with numerical samplings.
Collapse
Affiliation(s)
- Lior Zarfaty
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Alexander Peletskyi
- Institute of Physics, University of Augsburg, Universitätsstrasse 1, 86135 Augsburg, Germany
| | - Eli Barkai
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Sergey Denisov
- Department of Computer Science, Oslo Metropolitan University, 0130 Oslo, Norway
| |
Collapse
|
7
|
Klages R, Gallegos SSG, Solanpää J, Sarvilahti M, Räsänen E. Normal and Anomalous Diffusion in Soft Lorentz Gases. PHYSICAL REVIEW LETTERS 2019; 122:064102. [PMID: 30822076 DOI: 10.1103/physrevlett.122.064102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Revised: 11/02/2018] [Indexed: 06/09/2023]
Abstract
Motivated by electronic transport in graphenelike structures, we study the diffusion of a classical point particle in Fermi potentials situated on a triangular lattice. We call this system a soft Lorentz gas, as the hard disks in the conventional periodic Lorentz gas are replaced by soft repulsive scatterers. A thorough computational analysis yields both normal and anomalous (super)diffusion with an extreme sensitivity on model parameters. This is due to an intricate interplay between trapped and ballistic periodic orbits, whose existence is characterized by tonguelike structures in parameter space. These results hold even for small softness, showing that diffusion in the paradigmatic hard Lorentz gas is not robust for realistic potentials, where we find an entirely different type of diffusion.
Collapse
Affiliation(s)
- Rainer Klages
- Queen Mary University of London, School of Mathematical Sciences, Mile End Road, London E1 4NS, United Kingdom
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstraße 36, 10623 Berlin, Germany
- Institute for Theoretical Physics, University of Cologne, Zülpicher Straße 77, 50937 Cologne, Germany
| | - Sol Selene Gil Gallegos
- Queen Mary University of London, School of Mathematical Sciences, Mile End Road, London E1 4NS, United Kingdom
| | - Janne Solanpää
- Computational Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Mika Sarvilahti
- Computational Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Esa Räsänen
- Computational Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| |
Collapse
|