Efficient algorithm for simulating particles in real quasiperiodic environments

TL;DR

This paper presents an efficient algorithm based on the Generalized Dual Method for simulating particle dynamics in quasiperiodic environments, significantly reducing computational resources needed for large-scale simulations.

Contribution

The authors introduce a novel algorithm that enables efficient simulation of particles in quasiperiodic lattices without periodic approximations, applicable to any cut-and-project generated lattice.

Findings

The algorithm maintains constant computation time and memory with increasing distance R.

Free path length distribution depends on the rank of the quasiperiodic system, not symmetry.

Distribution combines exponential decay and power-law behavior, approaching exponential for high rank.

Abstract

We introduce an algorithm based on Generalized Dual Method (GDM) to efficiently study the dynamics of a particle in quasiperiodic environments without the need to use periodic approximations or to save the information of the vertices that make up the quasiperiodic lattice. We show that the computation time and the memory required to find the tile in which a particle is located as a function of the distance $R$ to the center of symmetry remains constant in our algorithm, while using the GDM directly both quantities go like $R^2$.This allows us to perform realistic simulations with low consumption of computational resources. The algorithm can be used to study any quasiperiodic lattice that can be produced by the cut-and-project method. Using this algorithm, we have calculated the free path length distribution in quasiperiodic Lorentz gases reproducing previous results and for systems with high symmetries at the Boltzmann-Grad limit. We have found for the Boltzmann-Grad limit, that the distribution of free paths depends on the rank $r$ of the quasiperiodic system and not on its symmetry. The distribution as a function of the free path length $l$ appears to be a combination of exponential decay and a power-law behavior. The latter seems to become important only for probabilities less than $(2^{r-2} r (r+1))^{-1}$, showing an exponential decaying free-path length distribution for $r \rightarrow \infty$, similar to what is observed in disordered systems.