Intermediate scattering function of colloids in a periodic laser field

Abstract

We investigate the dynamics of individual colloidal particles in a one-dimensional periodic potential using the intermediate scattering function (ISF) as a key observable. We elaborate a theoretical framework and derive formally exact analytical expressions for the ISF. We introduce and analyze a generalized ISF with two wave numbers to capture correlations in a periodic potential beyond the standard ISF. Relying on Bloch's theorem for periodic systems and, by solving the Smoluchowski equation for an overdamped Brownian particle in a cosine potential, we evaluate the ISF by numerically solving for the eigenfunctions and eigenvalues of the expression. We apply time-dependent perturbation theory to expand the ISF and extract low-order moments, including the mean-square displacement, the time-dependent diffusivity, and the non-Gaussian parameter. Our analytical results are validated through Brownian-dynamics simulations and experiments on 2D colloidal systems exposed to a light-induced periodic potential generated by two interacting laser beams.

---

1 Introduction

Colloidal particles suspended in fluids serve as a powerful model system in soft matter and biological physics. Their microscopic size makes them easy to manipulate, allowing researchers to explore fundamental physical principles in a controlled environment.^1–7 A major breakthrough in this field was the invention of optical tweezers by Ashkin et al., which enabled the precise trapping and manipulation of individual particles using focused laser beams.^8–10 Building on this advancement, scientists have since developed laser interference patterns to create structured external potentials, offering a unique way to study how particles behave under spatially periodic forces. These optical potentials arise from the interaction between light and particles with a refractive index different from that of the surrounding medium, creating highly tunable energy landscapes.^11,12

In order to understand how colloidal particles move within these periodic potentials, researchers have analyzed the probability distribution, the first-passage time, and low-order moments, such as the mean-square displacement, the long-time diffusivity, and the non-Gaussian parameter.^4,6,13,14 Although these observables provide valuable information, most theoretical studies have focused only on diffusion coefficients in the short- and long-time limits.^15–17 More recent work has expanded on these studies by analytically exploring tilted washboard potentials,^18,19 comparing theoretical predictions with experimental results,²⁰ and investigating memory effects in such systems.^21,22 These low-order moments are useful for identifying deviations from free Brownian motion, revealing effects such as trapping and non-Gaussianity.^23,24 However, they offer only a limited perspective of the system dynamics, as they do not capture the full spatial-temporal evolution of particle motion. To achieve a more complete description, it is necessary to adopt a framework that incorporates both spatial and temporal correlations.

One such framework emerges naturally in Markovian systems, where the future state depends only on the present and not on the past. In these systems, all relevant dynamical information is contained in the propagator, which describes the probability of a particle transitioning between states over a given time interval. Although the propagator provides a full statistical description of particle motion, it is often challenging to access experimentally. A more practical and experimentally accessible alternative is the intermediate scattering function (ISF), which encodes both spatial and temporal correlations in particle motion. Unlike traditional low-order moments, the ISF offers a more comprehensive characterization of the system and can be directly measured using techniques such as dynamic light scattering,^25,26 differential dynamic microscopy,²⁷ and single-particle tracking.²⁸

Furthermore, upon a small-wave-vector expansion of the ISF the low-order moments are recovered. Analytical derivations of the ISF have been achieved for various systems, including anisotropic active Brownian particles,^23,29 Brownian circle swimmers in gravitational fields,^30–32 anisotropically diffusing colloidal dimers,³³ and run-and-tumble agents.^34,35 Notably, experimental validations have been conducted for colloidal dimers and active colloids.

In this work, we extend the traditional ISF by introducing a generalized version that incorporates two wave numbers, allowing us to investigate the correlations in periodic systems in reciprocal space. The central question we address is twofold. First, we develop a theoretical framework to gain analytic insight into the dynamics of colloidal particles in general periodic potentials. Second, we compare our predictions with experimental results, exploring a spatio-temporal regime that has not yet been fully investigated.

To develop the theoretical framework, we solve the Smoluchowski equation for a single, overdamped Brownian particle in a cosine potential. Reformulating the problem as a Hermitian Schrödinger equation allows expressing the solutions in terms of eigenvalues and eigenfunctions.³⁶ Taking advantage of the systems periodicity, we apply Bloch's theorem, which provides a systematic way to expand the ISF using time-dependent perturbation theory and extract key dynamical moments. These methods are based on previous work^31,32,37 and have been successfully applied to bistable periodic potentials.^38,39 To validate our theoretical predictions, we compared our theoretical framework with Brownian dynamics simulations.

Experimentally, we track individual colloidal particles confined between two walls in a two-dimensional configuration, subjected to a periodic potential created by the interference of two laser beams. By systematically varying the laser power, we control the potential amplitude and record particle trajectories under different conditions. The generalized ISF, along with its associated moments, is then extracted from these experimental data and directly compared with our theoretical predictions, enabling a precise evaluation of the Brownian motion description in periodic optical fields. By integrating theoretical modeling, experimental measurements, and computational simulations, this work aims to provide a comprehensive understanding of colloidal motion in structured environments. Our approach bridges the gap between fundamental stochastic dynamics and experimentally accessible observables, offering new insights into the interplay between Brownian motion and periodic potentials.

This work is organized as follows. In Section 2 the experimental materials and methods are described. In Section 3, we introduce the analytical model, and derive the theoretical framework and observables. Readers more interested in the results rather than the theoretical derivation may skip this section and jump directly to Section 4, where we present our findings, compare experimental results with theoretical predictions, and discuss their implications. Finally, we summarize our findings, concluding with an outlook on future research directions in this field in Section 5.

2 Material and methods

This section describes the employed numerical methods and the used experimental setup to study colloidal dynamics in periodic potentials. We describe the processes of preparing two-dimensional colloidal samples and explain our custom-built optical setup that generates periodic laser fields. Also, the experimental protocol, including data acquisition parameters and particle tracking methodology, is described. Our experimental design allows us to systematically investigate the behavior of colloidal particles under varying potential strengths while maintaining high spatial and temporal resolution, which is essential for a rigorous comparison with the theoretical framework presented in Section 3.

2.1 Numerical methods

To complement our analytical and experimental approach, we conducted Brownian-dynamics simulations of single particles in a one-dimensional cosine potential. We implemented the Euler–Maruyama method⁴⁰ for numerically solving the differential equation of motion. For efficient sampling of observables, we used the extended order-n algorithm of Frenkel and Smit, which enables equidistant sampling on a logarithmic time scale.^41–43

To numerically solve the equation of motion, as it will be explained in detail in the theory Section 3, we employed a spectral method. Specifically, we expanded the equation in its Fourier basis to obtain its operator form. Numerically, we truncated the expansions to form a finite-dimensional matrix eigenvalue problem. We used the scipy.linalg.eig function from the SciPy library⁴⁴ to compute the eigenvalues and eigenvectors of the resulting matrix. The convergence of our results was verified by systematically increasing the truncation order until the eigenvalues stabilized to the desired precision.

2.2 Sample preparation

The experiments were performed with two-dimensional (2D) colloidal systems to ensure controlled particle motion mainly in the horizontal plane. We prepared a dilute suspension of monodisperse polystyrene sulfate latex particles (radius 1.5 μm, 4% polydispersity, Thermo Fisher Scientific, batch number 1660463) in ultra-pure water (resistivity 18.2 MΩ cm, Purelabs Flex, Elga). To build the sample cells chamber, we first attached a rectangular coverslip (Thickness No. 1, 24 × 50 mm, VWR 631-0146) to a microscope slide using UV-curing adhesive (NOA61, Norland Products Inc.). We then carefully deposited 2.5 μL of the particle suspension at the center of this base coverslip. A smaller square coverslip (Thickness No. 1.5, 22 × 22 mm, VWR 631-0125) was gently placed on top to create a thin chamber. We sealed the edges with UV-curing adhesive to prevent evaporation during long observation periods. To ensure enough space between the coverslips, particles with a radius of 2 μm were added to the colloidal suspension to act as spacers. These spacers ensure a sufficient gap between coverslips that counteracts the pull of capillary forces, preventing the coverslips from coming close to each other and squeezing the main particles. The resulting configuration reached an area fraction of 0.01. The completed cell was positioned on the microscope stage for observation.

2.3 Experimental setup

Our experimental setup consists of a laser field generator and an imaging system, as illustrated in Fig. 1. The laser field generator creates the periodic potential experienced by the colloidal particles. A 532 nm laser (Cobolt 05-01 Samba, 1.5 W) beam passes through a beam expander and is directed by mirrors (M1 and M2) to a Köster prism (KP). This prism splits the incoming beam into two parallel beams of equal intensity, which are then focused by lens L1 and directed by a dichroic mirror (D1) into the sample plane. At the sample plane, the beams interfere to generate the laser field inducing a periodic potential, with a periodicity of L = (4.15 ± 0.07) μm. By adjusting the position of the Köster prism, we can fine-tune this periodicity. Additionally, by controlling the laser power, we can precisely modulate the amplitude of the potential experienced by the particles. The relationship between the amplitude of the potential as a function of laser power (calibration) is presented in Appendix A.

Fig. 1 Illustration of the experimental setup. The laser field generator utilizes a 532 nm laser beam to create a periodic light field through a Köster prism (KP). The two parallel beams from KP are focused on the sample plane using lens L1 and their interference creates the fringe pattern. The imaging system, consisting of a microscope with a 20× objective and a CMOS camera, captures the sample's particle trajectories.

The imaging system is an inverted bright-field microscope (Nikon Ti-E) with a 20× objective (Nikon Plan Flour, 0.5 NA). We illuminate the sample with a blue LED (Thorlabs M455L4) and capture images using a CMOS camera (Mako U-130B) at a resolution of 1280 × 1024 pixels, with a pixel pitch of 0.24 μm per px. To prevent laser light from reaching the camera sensor, we employed a second dichroic mirror (D2) that redirects the laser beams to a beam dump (BD), with any residual laser light being filtered out by a notch filter (NF). This optical arrangement allows us to simultaneously apply the periodic potential and observe particle dynamics with high precision.

2.4 Experimental realization

We conducted all experiments in a temperature-controlled laboratory environment maintained at 20 °C. Before measurements, samples were allowed to equilibrate for 48 hours to ensure thermal stability. For each experiment, we recorded 180 000 images over 2 hours, using a frame rate of 25 Hz and an exposure time of 1 ms. These parameters represent a balance between temporal resolution and the need to collect sufficient data for robust statistical analysis. We collected between 400 and 1500 individual particle trajectories for each laser power from several measurements. At the highest laser power, the pressure radiation pushed approximately 20% of particles out of the field of view, requiring us to sample multiple regions within each sample cell to obtain statistically meaningful datasets. To extract particle positions with high precision, we employed a custom MATLAB-based tracking algorithms adapted from the work of D. Blair and E. Dufresne⁴⁵ together with the Michalet algorithm,⁴⁶ achieving a localization uncertainty of approximately 3 nm. Trajectories are openly available in Zenodo at https://doi.org/10.5281/zenodo.14931759.

This exceptional spatial resolution is crucial for accurately capturing the subtle dynamics of particles within the periodic potential, particularly when comparing experimental results with theoretical predictions. The observables are computed for each measurement and averaged over all the measurements. Results in Section 4, represent the mean values and the error bars their standard deviation.

3 Theory and observables

In this section, we introduce the model and the theoretical framework. Further, we derive the generalized ISF and low-order moments. From the expansion of the ISF, we compute exact expressions of the mean-square displacement, diffusivity, and non-Gaussian parameter.

3.1 Model

We model the interaction of the colloid and the interfering laser beam by a one-dimensional periodic potential U(x) = U(x + L) with period L. We assume a simple cosine potential

U(x) = U₁ cos(Q₁x), (1)

where Q₁ = 2π/L is the wave number associated with the period L. We introduce the dimensionless amplitude u = U₁/k_BT where k_BT is the thermal energy scale. The colloid is assumed to undergo Brownian motion in the presence of the spatially periodic force −∂_xU(x). Then the equation of motion for its position x(t) is provided by the Langevin equation

ẋ(t) = DQ₁u sin(Q₁x(t)) + η(t). (2)

Here, D denotes the short-time diffusion coefficient such that D/k_BT is the mobility as derived by the Einstein relation. The stochastic term η(t) represents a centered Gaussian white noise with a delta-correlated variance 〈η(t)η(t′)〉 = 2Dδ(t − t′).

We identify three quantities which set the characteristic units of the system: The period L of the potential is the fundamental unit of length. The time a free particle needs to diffuse the distance of a period is L²/D and will be used as a time unit. Energies are compared to the thermal energy k_BT. Therefore, the problem displays the dimensionless amplitude u as a single control parameter.

3.2 Smoluchowski equation

To make analytic progress, we rely on the Smoluchowski equation which is equivalent to the Langevin description in eqn (2) in the case of equilibrium dynamics. Here, the fundamental quantity is the propagator [Doublestruck P] := [Doublestruck P](x,t|x₀), defined as the conditional probability of finding the particle at position x at time t, given its initial position was x₀ at time zero. The initial condition is therefore provided by [Doublestruck P](x,t = 0|x₀) = δ(x − x₀). By standard methods³⁶ one derives the Smoluchowski equation governing the time evolution of the propagator

∂_t[Doublestruck P] = −∂_x[DQ₁u sin(Q₁x)[Doublestruck P]] + D∂_x²[Doublestruck P] =: Ω[Doublestruck P]. (3)

A stationary solution is provided by

with a normalization factor Z₁. Since in the infinite system no stationary solution exists, we choose Z₁ as the normalizing factor for a unit cell

Explicitly

where I_ν(·) denotes the modified Bessel function of the first kind to order ν.

We replace our infinite system with a finite system by dividing it into N ∈ [Doublestruck N] unit cells of length L and apply periodic boundary conditions. The limit N → ∞ will be performed at the end of the calculations.

In order to find non-trivial solutions of eqn (3) we first perform a separation ansatz

[Doublestruck P] = E(t)ψ(x), (7)

and find the solution for the time-dependent part immediately as E(t) = e^−λt. For the position-dependent part, we have to solve an eigenvalue equation. As the Smoluchowski operator Ω is non-Hermitian, we distinguish between right and left eigenfunctionsand Ω^† is the adjoint operator with respect to the scalar product

Left and right eigenfunctions to different eigenvalues are orthonormal, satisfying

Therefore, only the product of left and right eigenstates to identical eigenvalue is normalized. Furthermore, the eigenfunctions are complete, fulfilling the condition

By comparing eqn (3) and (8) for the stationary state, we identify that the eigenfunction ψ^R₀(x) to eigenvalue zero has to be proportional to the equilibrium distribution ψ^R₀(x) ∝ p^eq(x). We choose

ψ^R₀(x) = p^eq(x), ψ^L₀(x)* = 1, (12)

which fulfills the normalization conditions eqn (5) and (10).

The formal solution of the Smoluchowski, eqn (3), can be written as [Doublestruck P] = e^Ωtδ(x − x₀). We can insert the completeness relation, eqn (11), and apply the eigenvalue equation, eqn (8), and obtain

3.3 Schrödinger form

It is favorable to transform the Smoluchowski equation, eqn (3), into a Schrödinger-like equation using the ‘gauge transform’³⁶

[Doublestruck P] = e^{−U(x)/2k_BT}[Doublestruck P]₀. (14)

A straightforward calculation reveals that this yields in general to

By specializing to the potential given in eqn (1), we find

We note that the operator [script L]₀ is Hermitian with respect to the standard scalar product in eqn (9). Equivalently to the procedure of finding solutions of the Smoluchowski operator, to find the non-trivial solution of eqn (16) we perform again a separation ansatz

[Doublestruck P]₀ = E(t)Ψ(x), (17)

with E(t) = e^−λt and for the position-dependent part, we have to solve an eigenvalue equation

[script L]₀Ψ_λ(x) = −λΨ_λ(x), (18)

where λ represents the eigenvalue and Ψ_λ(x) the eigenfunction. We note that the transformation of the eigenfunctions of the Schrödinger operator to the Smoluchowski operator eqn (8) is provided by

As [script L]₀ is an Hermitian operator, the eigenvalues are real and two eigenfunctions with different eigenvalue are orthonormal, equivalently to eqn (10). The eigenfunctions are complete, fulfilling the completeness relation, similarly to eqn (11). The eigenfunction to eigenvalue zero can be easily found using eqn (12) and (19)

3.4 Bloch representation

For periodic systems we make use of the Bloch representation of the eigenfunctions in terms of the discrete wave vector q ∈ (2π/NL)[Doublestruck Z], with −π/L < q ≤ π/L and the discrete band index n. The wave vector can be restricted to the first Brillouin zone (BZ) because of the periodicity of the systems. The eigenfunctions are of the form

Ψ_nq(x) = e^iqxu_nq(x), (21)

where the Bloch function u_nq(x) = u_nq(x + L) is periodic, and the wave functions obtain a second index q, representing the wave number. The associated eigenvalue will be denoted by λ_nq. The orthonormality relation for the Bloch functions can be expressed assee Apendix C for the detailed derivation. Also, the completeness relation still holdswhere x and x₀ are within the same unit cell. Finally, we obtain the probability density using eqn (13) and transforming the eigenfunction of the Smoluchowski operator to the eigenfunctions of the Schrödinger operator with eqn (19). Lastly, inserting the Bloch form of the eigenfunctions we find the probability densitywhere in the last step the thermodynamic limit N → ∞ was performed and therefore, the sum over the wave numbers can be replaced by an integral over the Brillouin zone (BZ).

3.5 Dirac notation

To make further progress, we introduce the compact Dirac notation where we rely on the isomorphism between the periodic square-integrable functions u_nq(x) ∈ L²[0,L] and the abstract kets |u_nq〉 in a separable Hilbert space [script letter H]. We introduce the generalized position basis |x〉, such that u_nq(x) = 〈x|u_nq〉. In particular, the set of associated eigenstates is orthonormaland completefor a given wave number q. From the completeness relation in real spacewe infer 〈x|x₀〉 = δ(x₀ − x). From eqn (25) the (over-) completeness relation for the basis states |x〉 follows

For the actual computation of the eigenfunctions we use the Fourier modes as orthonormal basis {|ν〉: ν ∈ [Doublestruck Z]} in [script letter H] with real-space representation . It is favorable to express the Bloch functions u_nq(x) in terms of their Fourier decomposition, and we express our eigenmodes in Dirac notation

where the corresponding Fourier coefficients are provided by the integral

This is possible because all Bloch functions are lattice periodic.

3.6 Intermediate scattering function: definition and properties

We aim to analyze the characteristic function of the random displacements Δx(t) := x(t) − x(0), which corresponds to the self-ISF and provides full spatio-temporal resolution of the particle dynamics. Making use of the periodicity of the system, we introduce the Bravais lattice, defined as Λ := {nL: n ∈ [Doublestruck Z]}. The reciprocal lattice is similarly defined as Λ* := {Q_μ = (2πμ/L): μ ∈ [Doublestruck Z]}. Any wave vector k can be uniquely decomposed as k = q + Q_μ, where q ∈ BZ lies within the first Brillouin zone (BZ) and Q_μ ∈ Λ* is a reciprocal lattice vector. We define the generalized ISF

F_μν(q,t) := 〈e^{−i(q+Q_μ)x(t)}e^{i(q+Q_ν)x(0)}〉, (31)

with the mode indices μ,ν ∈ [Doublestruck Z]. In the following we refer to this quantity as the ISF. The brackets 〈⋯〉 indicate an ensemble average. The conventional ISF to wave number q + Q_μ then corresponds to F(q + Q_μ,t) := F_μμ(q,t) and probes diagonal correlations in reciprocal space. The off-diagonal elements F_μν(q,t), μ ≠ ν, encode Umklapp-processes where wave vectors differ by a reciprocal lattice vector Q_ν − Q_μ. These Umklapp-processes enter by scattering from the periodic modulation, while they are absent in the homogeneous case. Thus, F_μν(q,t) corresponds to a matrix-valued correlation function: for a set of complex numbers b_μ(q) ∈ [Doublestruck C] the weighted sumforms an autocorrelation function, with . In particular, autocorrelation functions display non-negative spectra.⁴⁷ Notably, the diagonal elements, F_μμ(q,t), are autocorrelation functions. For the case of purely relaxational dynamics, for example, Smoluchowski dynamics, autocorrelation functions are completely monotone, i.e., all time derivatives exist and satisfy the inequality (−1)^m∂^m_t〈b(q,t)*b(q,0)〉 ≥ 0 for m ∈ [Doublestruck N]₀ and t > 0,⁴⁸ ensuring monotonically decaying, non-oscillatory behavior. In contrast, the individual off-diagonal elements, μ ≠ ν, can exhibit non-monotonic behavior, including local minima or maxima.

Further, we note that the wave number q is identical in both exponentials and that the diagonal elements, F_μμ(q,t), correspond to the conventional ISF evaluated at wave vector q + Q_μ. In translationally invariant systems only the diagonal elements are non-vanishing, since shifting the trajectory of a particle by an arbitrary displacement leads to an equally likely trajectory. For our case, the discrete symmetry is reflected in the sense that a common shift x(t) ↦ x(t) + R for all times by a lattice vector R ∈ Λ leaves the ISF invariant. The brackets 〈⋯〉 indicate an ensemble average.

Using the conditional probability density, the ISF of eqn (31) can be expressed as

Here, we used that without restriction, the initial position of the particle can be chosen to be in a definite cell and is sampled from the equilibrium distribution for this single cell.

Reversely, we can compute the probability density by the backward Fourier transform

see Appendix D for the derivation. The previous relation also reveals that the conventional ISF is not sufficient to reconstruct the full probability density unless the system is fully translationally invariant.

The explicit form of eqn (33) in terms of the Bloch functions is obtained by inserting eqn (24) and simplifies to express the stationary solution in terms of the eigenfunctions of the Schrödinger operator, using eqn (25). After some algebra and rearranging the terms, we find the final expression

see Appendix D for more details.

To determine the functions u_nq(x) we introduce the operator [script L]_q for which the eigenvalue equation

[script L]_qu_nq(x) = −λ_nqu_nq(x), (36)

holds. Straightforward substitution leads towhich can be abbreviated as [script L]_q = [script L]₀ + δ[script L]_q. Here, the operator [script L]₀ (first three terms) encodes the interaction of the diffusive particle with the potential and the q-dependent operator δ[script L]_q (last two terms) contains a drift or advection-like term, linear in the wave number, and a diffusion-like term proportional to q².

The matrix representation of the operator, is given by

and

〈μ|δ[script L]_q|ν〉 = (−4πqDν/L − Dq²)δ_μ,ν. (39)

The matrix [script L]₀ is a Hermitian matrix and pentadiagonal in the Fourier basis, i.e., it has non-zero elements only on the main diagonal and the two diagonals above and below it. The q-dependent matrix is diagonal in the Fourier basis. The eigenvalue problem [script L]_q|u_nq〉 = −λ_nq|u_nq〉 is computed numerically by diagonalizing the (truncated) Hermitian matrix

The time evolution of the ISF is encoded in the eigenvalues and eigenfunctions of the operator [script L]_q. The eigenvalues λ_nq form continuous bands as the number of cells goes to infinity, N → ∞, see Fig. 2. All eigenvalues are non-negative, and the only zero eigenvalue is in the center of the Brillouin zone at the lowest band. Only the lowest bands are significantly affected by the potential U(x). For λ ≳ π²Du/L² the bands correspond to a particle freely diffusing without underlying potential modulation. Albeit the bands come very close at the edges of the Brillouin zone, we checked numerically that they do not touch. For symmetric potentials the avoided crossing theorem does not apply and in principle bands can cross. For the simple cosine potential, one can actually show that all eigenvalues for q = 0 except for λ₀₀ = 0 are twofold degenerate, see Appendix C.

Fig. 2 Eigenvalues λ_nq of the operator [script L]_q for U₁ = 1.0k_BT are shown, for the five lowest bands in the first Brillouin zone.

Finally, the ISF in eqn (35) can be conveniently expressed in a spectral representation using the Fourier basis

see Appendix D for the derivation. The previous relation reveals that is completely monotone.

Of particular interest is the conventional ISF with wave vector in the first Brillouin zone, F(q,t) := F₀₀(q,t), which simplifies upon exploiting the completeness of the Fourier basis

Since all eigenvalues are strictly larger than zero, except for λ₀₀ = 0 in the lowest band and in the center of the Brillouin zone, F(q,t) decays to zero for large time, t → 0 for q ≠ 0.

3.6.1 Symmetry relations. The symmetries of the generalized ISF are derived following the same arguments as in ref. 49. For equilibrium dynamics, time inversion symmetry holds and the ISF is even in time. Time-translation invariance then reveals the relations

F_μν(q,t) = F_μν(q, −t) = F_νμ(q,t)*. (43)

For symmetric potentials U(−x) = U(x), space-inversion symmetry implies

F_μν(q,t) = F_μν(q,t)* = F_−μ,−ν(−q,t). (44)

In particular, the ISF is a real quantity and symmetric upon interchanging its mode indices. At the edge of the Brillouin zone we have the additional relation

F_μν(π/L,t) = F_μ+1,ν+1(−π/L,t) = F_{−(μ+1),−(ν+1)}(π/L,t). (45)

3.6.2 Short- and long-time limits. For the simple cosine potential, eqn (1), the short-time limit of the ISF can be calculated explicitly

From eqn (35) and , one infers that for q = 0 the ISF displays a non-vanishing long-time limit

F_μν(0,t → ∞) = 〈e^−iQ_μx(t)〉〈e^iQ_νx(0)〉. (47)

The factorization of the limit can be interpreted as the system being ergodic. For the simple cosine potential, eqn (1), the limit can be calculated explicitly

An equivalent formal expression for the long-time limit follows from eqn (41)

Since and employing the Jacobi–Anger expansion,⁵⁰ the Fourier coefficients 〈μ|u₀₀〉 can be calculated explicitly

With Neumann's addition theorem for modified Bessel functions⁵⁰ the sums in eqn (49) can be performed and we recover eqn (48). Note that eqn (46) and (48) represent static quantities that only depend on the potential amplitude.

3.7 Mean-squared displacement and non-Gaussian parameter

The goal of this subsection is to elaborate the low-order moments of the fluctuating variable Δx(t). They can be derived from the conventional ISF, F(q,t) = F₀₀(q,t), which is the characteristic function of the random displacements. As moments of odd order vanish by symmetry, we compute only the even moments. The lowest nontrivial is the time-dependent mean-squared displacement 〈[Δx(t)]²〉. We also define the time-dependent diffusion coefficient via the derivative of the mean-squared displacement (MSD)

The long-time dynamics is diffusive, in particular, the long-time limit D_∞ := D(t → ∞) > 0 is finite and defines the long-time diffusion constant. An analytic expression for arbitrary periodic potentials is known.^{15,16,18,22,51} In particular, for the simple cosine potential it evaluates to

The exponential suppression of the diffusion constant reflects Kramers's rule for hopping over a potential barrier.³⁶ A convenient measure to discuss deviations from simple diffusion is the non-Gaussian parameter^47,52

The derivation for extracting the moments from the ISF is similar to ref. 31 and 32 with only the essential steps presented here. The key idea is to apply perturbation theory for small wave numbers q and compare it to the Taylor series of the ISF, which yields the moments

where we already exploited that all odd-order moments vanish. The ISF, eqn (42), will be expanded in powers of q by [script L]_q = [script L]₀ − Dq² + δ[script L]¹_q with δ[script L]¹_q = 2iqD∂_x. We rely on the Dyson representation

Replacing the time-evolution operator in the integral on the right-hand side iteratively generates the Born series, see ref. 32. The main simplification steps are to make use of the fact that e^{[script L]₀t}|u_n0〉 = |u_n0〉 and also 〈u_n0|e^{[script L]₀t} = 〈u_n0| and to insert complete basis sets, eqn (26), for q = 0. Occurring integrals can be formally evaluated and the terms are simplified to obtain the final result, which is similar to the result in ref. 32, but slightly changed for our operator and eigenvectors. We find the formal expression

with
and,

Here, all sums over n, m, p include all bands and therefore formally the expression causes divisions by zero if a band index corresponds to the lowest band or two band indices correspond to the same band. In both cases the corresponding numerators also vanish. The appearance of the zero divisors can be avoided in the first place by handling these case separately before performing the integrals in the simplification steps. Here we follow a different route to keep the expressions simple by analytically continuing the expression for the case of zero numerators/denominators.

A further complication arises in the case of a simple cosine potential, since all eigenvalues λ_n0 are, additionally, twofold degenerate, except for the ground state causing additional zero divisors. However, as in degenerate perturbation theory, one can choose basis states such that the matrix elements of δ[script L]¹_q coupling different states to the same eigenvalues vanish. Since [script L]¹_q anticommutes with space inversion, only states of different parity couple, however, because u_n0 are either even or odd, no zero divisors occur.

The low-order cumulants of the random variable Δx(t) are generated upon expanding the logarithm of the ISF in powers of the wave number q

To order O(q²) we find the mean-square displacement as first nonvanishing cumulant

where no contributions from n = 0 as intermediate state arises since the transition matrix element vanisheswhere the last equality follows by integration by parts and observing that u₀₀(x) is real. The time-dependent diffusion coefficient, eqn (51), can be computed using the time-derivative of eqn (59).

For the fourth cumulant we need to collect terms of order O(q⁴) and we obtain

〈[Δx(t)]⁴〉 − 3〈[Δx(t)]²〉² = F̃₄(t) − 3[F̃₂(t)]². (62)

For completeness, let us argue explicitly that all odd powers in q in the expansion of F(q,t), eqn (56), vanish. A term linear in q corresponding to a mean drift would involve the matrix element 〈u₀₀|δ[script L]¹_q|u₀₀〉 which is shown to vanish in eqn (60). This vanishing of the mean drift is, of course, a general property in equilibrium. For any symmetric potential all odd moments vanish. The expansion of F(q,t) in eqn (56) generates a chain 〈u₀₀|δ[script L]¹_q|u_n0〉 〈u_n0|…|u_p0〉 〈u_p0|δ[script L]¹_q|u₀₀〉 in eqn (57) of products of matrix elements. By parity and the property of the operator δ[script L]¹_q a matrix element is non-vanishing only if the states are of different parity. Thus, for the chain to yield a non-vanishing contribution, the first intermediate state has to be odd, the second even, and so on. However, the last state is the ground state again which is even. Therefore, only even powers of q are generated.

4 Comparison of the theoretical framework with experimental and simulation results

In this section, we analyze the observables and compare the theoretical predictions with the experimental results. First, we examine the diagonal and off-diagonal elements of the ISF for various wave numbers and potential amplitudes. Next, we discuss the mean-square displacement (MSD), time-dependent diffusivity, and the non-Gaussian parameter for different potential amplitudes. To facilitate the comparison with the experimental results, we restore units; in particular, we use the amplitude of the modulation U₁, eqn (1), as a control parameter. The presented experimental results correspond to the mean of multiple measurements, with error bars indicating the standard deviation. Each measurement consisted of a sufficiently large number of recorded trajectories, from which the observable of interest was computed as an average. To enable comparison with theoretical predictions, we established the relationship between the potential amplitude U₁ and laser power (LP) (see Appendix A for details). Additionally, length and time scales were calibrated by determining the short-time diffusion coefficient D and the characteristic period L from the experimental data. We determined from the initial slope of the mean-square displacement an average short-time diffusion coefficient of (0.050 ± 0.002) μm² s⁻¹.

4.1 Diagonal ISF, μ = ν

The diagonal ISF is computed from our experimental data relying on the definition given in eqn (31) for μ = ν. We have checked that its imaginary part is in deed negligible reflecting the mirror symmetry of the potential. We compare the experimental results to the numerical ones relying on the spectral representation in eqn (41). Last, all results are corroborated by Brownian-dynamics simulations. The ISF is displayed in Fig. 3 for a range of different wave numbers Q_μ ∈ Λ* and q ∈ BZ for three distinct amplitudes U₁.

Fig. 3 Diagonal ISF for various mode indices μ = ν and wave numbers q and three different amplitudes U₁, (a)–(c). Full colored lines correspond to the theory and colored squares to the experimental results. Black circles represent the simulation results. The dotted lines represent the harmonic approximation.

We first focus on the behavior for q ≠ 0, where all ISF eventually relax to zero. For moderate potential amplitudes U₁ ≈ k_BT, the potential is not high enough to significantly inhibit hopping between different potential valleys, yielding a single-step relaxation. For larger amplitudes, a two-step process occurs. The particle initially freely diffuses with short-time diffusion coefficient D until the potential forces become dominant. For U₁ ≳ k_BT the motion occurs essentially at the bottom of the potential, which can be approximated by a harmonic well

The particle then locally equilibrates on the time scale of the harmonic relaxation time, τ = (L²/4π²D)(k_BT/U₁), and the ISF saturates at a plateau value, see Fig. 3. The ISF within the harmonic approximation can be calculated explicitly, see Appendix E. For large barriers, Fig. 3(c), the harmonic approximation quantitatively describes the relaxation to the plateau value for wave numbers resolving smaller length scale than a period L. The relaxation from the plateau occurs on a much larger time scale provided by Kramers’ theory τ_K ∝ exp(2U₁/k_BT). Once the particle overcomes the barrier and reaches additional minima, the ISF eventually decays to zero. For small wave numbers and long times, the hydrodynamic regime is reached F₀₀(q,t) = exp(−D_∞q²t). In this regime, the wave numbers only resolve the motion over many periods at time scales much larger than Kramers’ escape time. Our analytical predictions as well as the simulation results within the Smoluchowski picture of a simple cosine potential show excellent agreement with the experimental results.

For the wave number in the center of the BZ, q = 0, the ISF does not decay to zero in the long-time limit t → ∞, rather approaches a finite value, as computed in eqn (48). In Appendix A we used this feature to calibrate the laser power of the experiment to the theoretical k_BT value.

4.2 Off-diagonal ISF, μ ≠ ν

For the off-diagonal ISF, μ ≠ ν, we present results for different modes corresponding to |μ − ν| = 1 and |μ − ν| = 2. The results are shown for three different wave numbers q and the amplitude U₁ = 0.99k_BT, see Fig. 4. The general symmetries of the ISF for symmetric potentials are summarized in eqn (43) and (44).

Fig. 4 Off-diagonal ISF for ν ≠ μ and U₁ = 0.99k_BT for different mode indices (μ, ν) and wave numbers q, (a)–(c). Full colored lines correspond to the theory and squares to the experimental results. Black circles represent the simulation results. In panel (a) at the edge of the Brillouin zone, q = π/L, for (μ, ν) = (1, −1) and (0, −2) coincide, as well as the ones for (1, 0) and (−1, −2). In panel (c) for wave numbers q → 0 the curves of (μ, ν) and (−ν, −μ) approach each other, where μ, ν ∈ [Doublestruck Z].

For wave numbers at the edge of the BZ q = π/L additional symmetries of the ISF, eqn (45), hold. For example, the curves for (μ, ν) = (1, −1) and (0, −2), or (1,0) and (−1, −2) are identical.

For q = 0 both the initial value and the long-time limit are non-zero. As soon as q ≠ 0 the curves decay to zero for long times. If |μ − ν| is odd, the initial value is negative and for even differences the initial value of the ISF is positive, eqn (46). And also according to eqn (48) for odd (or even) values, the long-time limit is negative (or positive, respectively). In contrast to the diagonal elements of the ISF, we find no longer strictly monotone behavior but minima and maxima, see eqn (32).

Results from the experimental data of the off-diagonal ISF closely follow the shape of the analytical predictions and simulation results at long times. However, a clear deviation is observed at short and intermediate times, as shown in Fig. 4. We attribute these discrepancies to experimental factors not implemented in the theoretical model. We consider the primary factor to be the spatial inhomogeneities in the periodicity and amplitude of the light-induced potential across the field of view, which slightly deviates from the ideal cosine form assumed in the theoretical model. As a result, individual particles effectively experience slightly different periodicities and amplitudes across the field of view. A detailed characterization of the amplitude and periodicity in the field of view is shown in Appendix B. Additionally, confinement effects due to the two-dimensional nature of the system and inertial effects not considered in the theoretical framework are factors that might also contribute to these discrepancies. Despite these differences, experimental results for the indices μ and ν with equal |μ − ν| collapse into a common intercept (short-time limit), capturing the expected phenomenology from the theoretical predictions given by eqn (46), though with a noticeable shift respect to theory.

Theoretically, each Fourier component of the ISF is defined in terms of a well-defined lattice wavenumber 2π/L. However, in experiments, achieving this level of precision is challenging because of variations in how individual particles interact with the optical field. To allow a more accurate comparison between theory and experiment, a normalization procedure is applied, an approach that is further examined in Fig. 5.

Fig. 5 Off-diagonal ISF for different amplitudes U₁ for wave number q = π three different mode indices (μ, ν), (a)–(c). Full colored lines correspond to the theory and squares to the experimental results. Black circles represent the simulation results and dotted lines the harmonic approximation.

Fig. 5 illustrates the normalized ISF, F_μν(q,t)/F_μν(q,0), for different potential amplitudes U₁/k_BT, three different μ, ν combinations and, for all of them, qL = π. The choice of qL = π is particularly insightful, as it balances, for the measurements time window, sensitivity to both free diffusion and potential-induced localization, providing a clear distinction between different transport regimes. As it can be seen in Fig. 5, after normalization, the agreement between experimental results, the theoretical framework and computer simulations is remarkable, confirming that the model effectively captures full description of the system.

The most striking feature in the off-diagonal elements is the emergence of a maximum at intermediate times, τ ≪ t ≪ τ_K as seen in Fig. 5(b) and (c). For larger times, the curves decay to zero for q ≠ 0. From the harmonic approximation we anticipate the development of a plateau, whose value can be determined from eqn (77). If (Q_μ + q)(Q_ν + q) < 0, the plateau corresponds to a maximum, which is nicely approached for large potential amplitudes, see Fig. 5(b) and (c). If (Q_μ + q)(Q_ν + q) > 0, the curves look similar to the diagonal elements of the ISF, where a simple plateau emerges, see Fig. 5(a). The slowing down of the relaxation towards the plateau or maximum as the potential amplitude grows, is captured as well by τ ∝ 1/U₁.

Finally, it is important to emphasize that, although experimental factors cause deviations at short and intermediate times, normalization effectively accounts for these variations, leading to excellent agreement between theory, simulations and experiments.

4.3 Low-order moments of Δx(t): MSD, diffusivity and non-Gaussian parameter

The dynamics previously discussed in terms of the ISF, can also be found in the standard observables. Therefore, the discussion will be kept rather brief. The MSD, Fig. 6(a), is plotted for various amplitudes U₁ of the external potential. For short times, t ≪ (L²/4π²D)(k_BT/U₁), with U₁ ≳ k_BT or for times t ≪ L²/4π²D for very low amplitudes, U₁ ≪ k_BT, we observe the expected linear increase of the MSD, characteristic of free diffusion. For higher potential barriers, however, a plateau emerges around t ≈ (L²/4π²D)(k_BT/U₁), corresponding to the time scale at which the particle becomes temporarily trapped. At longer times, ∝ exp(2U₁/k_BT), the particle eventually overcomes the barrier and resumes diffusion, with a reduced diffusion constant for higher amplitudes compared to lower ones. This behavior is further highlighted in the time-dependent diffusion coefficient, eqn (51). We find that the long-time diffusion coefficient decreases as U₁ increases, see Fig. 6(b). There we also report good agreement with the known values of the long-time diffusion coefficient. As expected, the harmonic approximation captures the behavior of the relaxation towards the plateau increasingly better for higher amplitudes.

Fig. 6 (a) MSD, (b) time-dependent diffusion coefficient, and (c) non-Gaussian parameter for different potential amplitudes U₁. Full colored lines correspond to the theory and squares to the experimental results. Black circles represent the simulation results and dotted lines the harmonic approximation. The dashed lines in (b) display the theoretical long-term limit.

Furthermore, we analyze the non-Gaussianity of the particle displacements using the parameter defined in eqn (53). As expected, for higher barriers, the particle dynamics become increasingly non-Gaussian, see Fig. 6(c). We observe that both very small and very large amplitudes pose challenges in experiments. For small amplitudes, it is difficult to distinguish the dynamics from those of a free particle, as the external potential has little effect. Conversely, for very large amplitudes, the low diffusivity makes it challenging to sample a sufficient number of particles that successfully hop over a barrier within the experimental observation time. However, the experimental results show excellent agreement with both simulations and analytical predictions.

5 Conclusion

In this work, we investigated the dynamics of dilute colloidal suspensions in the presence of a periodic potential. Our approach combined theoretical analysis, Brownian-dynamics simulations, and experimental measurements. By evaluating the spatio-temporal dynamics, we demonstrated that the behavior of a colloid in a periodic potential can be accurately described by analytical solutions of the Smoluchowski equation. This was achieved through the analysis of a generalized ISF, which captures how particle positions are correlated in a periodic system. Low-order moments were derived, with a focus on the MSD, time-dependent diffusion coefficient, and non-Gaussian parameter.

Based on the Smoluchowski equation reformulated in a Hermitian Schrödinger form, we found formal expressions in terms of a spectral-theoretical approach. The eigenfunctions were expressed in Bloch form, to make use of the periodic nature of the system. We found an analytic expression for the generalized ISF, F_μν(q,t). By using the time-dependent perturbation theory and Taylor expansion of the ISF we computed lower-order moments. In our system, without memory effects, the ISF effectively captures the full dynamics of colloidal particles in periodic potentials, revealing both short- and long-time diffusive behavior and trapping at intermediate times.

We performed experiments on 2D dilute colloidal suspensions subjected to a periodic potential generated by two interfering laser beams. Using particle tracking, we obtained particle trajectories and averaged them to extract relevant observables. The laser power was calibrated to its corresponding amplitude value using two theoretical predictions (eqn (48) and (52)). From our experimental data, we extracted the observables of interest and compared them to our analytical solutions and Brownian-dynamics simulations. We compared the results for various strengths of the amplitude of the potential and found excellent agreement between the theoretical description and experiments. The most sensitive observables were the off-diagonal elements of the generalized ISF, where slight differences in the experimental setup were amplified in the curves. It was crucial to ensure equivalent experimental conditions, minimizing variations in periodicity, potential amplitude, and confinement effects. To obtain good agreement with the analytical predictions, a normalization was necessary.

We have carried out a comprehensive test of the underlying fundamental dynamics of colloidal dynamics in a structured environment, combining theory, simulations, and experiments. By analyzing the generalized ISF, we identified new observables with distinct features, including a non-vanishing long-time limit. We provide a detailed theoretical framework and rationalize our findings through a harmonic approximation. Developing an explicit formula for the whole-time dependence of the ISF allowing for a comprehensive description of colloidal motion across all time scales. Furthermore, we introduced a new approach for calibrating the experiment using these observables, offering a reciprocal space alternative to conventional calibration methods. Comparing the results to a harmonic approximation, we confirm that the Brownian particle first diffuses freely, before it is temporarily trapped in the minima of the periodic cosine potential. For these times the dynamics are well approximated by a harmonic potential for large enough amplitudes, and only at longer times it hops over the potential barriers and once again exhibits diffusive behavior.

The analytical and experimental framework presented can be extended to more complex systems. Although this work focused on dilute suspensions, exploring more dense systems would allow us to study particle interactions and many-body effects. A possible other extension is the study of periodic lattices, where higher-dimensional effects and collective behavior become important. Investigating tilted washboard potentials could provide further insight into driven transport. Our framework is also applicable to a wide range of periodic systems beyond simple cosine potentials. Experimentally, the new observables could also be measured using differential dynamic microscopy.

Author contributions

MAES and TF conceived and administrated the project. YMS, ABZB, and MAES conducted the experimental work. RR and TF performed the analytic modeling. RR implemented the numerical analysis. YMS, MAES, and RR performed computer simulations, analyzed the experimental data, and contributed to data visualization. RR drafted the original manuscript, and all authors contributed to the manuscript.

Data availability

Data for this article, including all particle trajectories analyzed in the manuscript, are available at Zenodo https://doi.org/10.5281/zenodo.14931759.

Conflicts of interest

There are no conflicts to declare.

Appendices

A. Calibration

The precise determination of the periodic potential amplitude (U₁/k_BT) as a function of laser power is a fundamental step in experiments involving colloidal particles confined in optical potentials. In this work, we perform both the traditional diffusion-based calibration and a novel approach utilizing the generalized intermediate scattering function (ISF). The conventional approach to calibrating periodic potentials relies on measuring the normalized long-time diffusion coefficient (D_∞/D) as a function of laser power. As shown in the central panel of Fig. 6, the diffusivity of particles decreases as the potential amplitude increases, since energy barriers progressively hinder diffusive motion. For a simple cosine potential, U(x) = U₁ cos(2πx/L), the analytical relationship between the normalized diffusion coefficient and the potential amplitude is given by eqn (52).

Our alternative calibration method utilizes the generalized ISF, specifically its asymptotic behavior, which explicitly depends on the potential amplitude. For a simple cosine potential and q = 0, the long-time limit of the generalized ISF is analytically described by eqn (48), where the dependence on the potential amplitude, U₁, is evident. Specifically, for μ = ν = 1, eqn (48) (F₁₁(0,t)) evaluates the characteristic wave number imposed on the system by the periodic potential, i.e., Q_μ = Q_ν = 2π/L. In the top panel of Fig. 7 we plot F₁₁(0,t), which exhibits an opposite trend to diffusivity, with higher plateaus as laser power increases. It is important to note that this approach utilizes equilibrium correlation properties instead of transport characteristics, providing complementary insights into the system.

Fig. 7 Panel (a) shows the ISF for the wave numbers (μ, ν, qL) = (1, 1, 0) and for different amplitudes U₁. The markers represent experimental results, while the solid lines correspond to analytical predictions, with the amplitudes calibrated to match the experimental data. Panel (b) compares two different calibration methods. One based on the long-time diffusivity D_∞ and the other on the long-time limit of the ISF, F₁₁(0, t → ∞). The laser powers (LP) used in the experiments, for which the calibration was performed, were 130 mW, 250 mW 400 mW, 500 mW, 800 mW, 1000 mW, 1300 mW, and 1550 mW.

The calibration is performed for both methods by extracting plateau values as a function of laser power and solving eqn (48) and (52) for the diffusivity and generalized ISF methods, respectively. The calibration results from both methods are presented in Fig. 7(b), showing the relationship between laser power and dimensionless potential amplitude U₁/k_BT. The remarkable agreement between these independent approaches validates both the theoretical framework and the experimental implementation, as seen in the central panel of Fig. 6 and 7(a). This strong consistency confirms the theoretical predictions and demonstrates the reliability of this method for calibration purposes.

B. Light-field spatial characterization

This appendix provides a detailed look at the spatial characteristics of the laser-induced periodic potential, as observed through the motion of particles across the experimental field of view. To determine the spatial periodicity, we follow a similar approach as described in ref. 53. However, for the estimation of the amplitude, we used a more direct estimation. Rather than fitting a sine wave to the derivative of the logarithmic particle density, we extracted the amplitude in each window by identifying its maximum and minimum values. The extracted amplitude of the periodic potential as a function of the position along the field of view is shown in Fig. 8(a). The amplitude displays a clear spatial dependence, which arises from the Gaussian envelope of the beam.^6,54,55 The mean amplitude was measured to U₁ = (1.03 ± 0.08) k_BT, with a coefficient of variation (CV) of 7.76%. The inset histogram shows the distribution of the measured amplitude values. Fig. 8(b) shows the spatial periodicity of the laser field in the same region. In contrast to the amplitude, the periodicity remains remarkably consistent throughout the field of view, as it is unaffected by the Gaussian beam envelope. It maintains an average value of L = (4.15 ± 0.07) μm, with a very low CV of 1.68%. The inset histogram confirms this stability by showing a narrow, normally distributed spread of periodicities. The error bars in both figures represent variations between different experimental runs in the same spatial regions.

Fig. 8 Panel (a): The spatial variation of the extracted potential amplitude U₁ (orange squares) along the field of view, showing a mean value of U₁ = 1.03k_BT (red dashed line) with a standard deviation of 0.08 (orange error bar). Inset: Histogram of the extracted amplitude values fitted with a Gaussian distribution (red curve). Panel (b): The spatial variation of the extracted periodicity L (green squares), with a mean value of L = 4.15 μm (blue dashed line) with a standard deviation of 0.07 (green error bar). Inset: Histogram of the extracted periodicities fitted with a Gaussian distribution (blue curve).

C. Properties of the Bloch functions

In this appendix, we repeat the argument that the Bloch function u_nq(x) to the same wave vector q are orthonormal. Furthermore, we recall some properties of the Bloch functions, in particular, for the case of a symmetric potential.

Orthogonality of Bloch functions. The overlap of two wave function can be simplified using the periodic Bloch functions and splitting the integral in cells

Since

for the discrete wave vectors in the Brillouin zone and because the wave functions are orthonormal 〈Ψ_nq′|Ψ_mq〉 = δ_qq′δ_nm, Bloch functions to the same wave vector are orthogonal

Symmetric band structure. The Hermitian operator [script L]₀ is real. We introduce the eigenfunction Ψ_λ(x) to [script L]₀Ψ_λ(x) = −λΨ_λ(x). Taking the complex conjugate of both sides reveals that Ψ_λ(x)* is again eigenfunction to the same eigenvalue. Bloch theorem states that eigenfunctions can be expressed as Ψ_nq(x) = exp(iqx)u_nq(x) with u_nq(x) periodic. Hence, u_n,−q(x) coincides with u_nq(x)* up to a phase factor. Without restriction, the phase factor can be chosen to be real. Therefore, the corresponding band structure is symmetric with respect to flipping the sign of the wave vector, λ_n,q = λ_n,−q. Furthermore, the eigenfunctions in the center of the Brillouin zone are real u_n0(x) = u_n0(x)*.

Degeneracies for symmetric potentials. We define the parity operator P acting on functions PΨ(x) = Ψ(−x). For a symmetric potential PU(x) = U(−x) = U(x) and the parity operator commutes with [script L]₀. Then, with Ψ_λ(x) eigenfunction, we find λPΨ(x) = P[script L]₀Ψ_λ(x) = [script L]₀PΨ_λ(x). Hence, PΨ_λ(x) = Ψ_λ(−x) is again eigenfunction to the same eigenvalue. For the Bloch representation this implies u_nq(−x) = ±u_n,−q(x). In particular for q = 0 we find u_n0(x) is either even or odd. As a consequence all matrix elements 〈u_m0|[script L]₀|u_n0〉 vanish if the eigenfunctions u_m0(x), u_n0(x) have different parity. In particular, the avoided crossing theorem does not apply, the eigenvalues at the center of the Brillouin zone can be twofold degenerate.

Peculiarity of the cosine potential. In general, band crossings at the center of the BZ are allowed for symmetric potentials since the avoided crossing theorem does not apply, but usually only some band crossings appear while other bands still avoid each other. The simple cosine potential is special in the sense that all eigenvalues to q = 0 are twofold degenerate except for the ground state λ₀₀ = 0.

This property is somewhat hard to see in the representation of the time-evolution operator in the Schrödinger representation in the Fourier basis 〈μ|[script L]₀|ν〉. However, the property can be easily deduced, omitting the gauge transform in the first place, i.e. representing the dynamics in terms of the non-Hermitian matrix

This matrix displays the symmetry 〈−μ|Ω|−ν〉 = 〈μ|Ω|ν〉*. The argument now follows the one of Appendix A of ref. 32. The matrix 〈μ|Ω|ν〉 displays a zero row for μ = 0 and splits into a part with entries for μ > 0, ν ≥ 0 and an identical one for μ < 0, ν ≤ 0.

The only entries preventing the matrix to split into blocks with positive/negative μ, ν are the matrix elements 〈±1|Ω|0〉. However, as 〈l₀| = 〈0| is a left eigenvector to eigenvalue 0, all eigenvectors |r_n0〉 to non-zero eigenvalues have a zero entry in their Fourier representation by orthogonality of eigenvectors 0 = 〈l₀₀|r_n0〉 = 〈0|r_n0〉. Therefore, the blocks with both μ,ν positive do not communicate with the blocks with both indices negative. In particular, one can choose eigenvectors with 〈μ|r_n〉 = 0 for μ ≷ 0 or symmetric and antisymmetric eigenfunctions to the twofold degenerate eigenvalue λ_n0 > 0.

D. ISF and probability density

Starting with the general definition of the ISF, eqn (33), and inserting the expression for the probability density, eqn (24), we obtain

Here in the second equality we used that the Bloch function to wave vector zero at the lowest band is related to the equilibrium density . Furthermore, we observed that the integral vanishes for q ≠ q′ due to the periodicity of the Bloch functions.

We can make further progress by using the Fourier modes as basis functions, eqn (29). By expanding the Bloch functions we obtain

Collecting terms, we find the expression for the ISF

This relation is eqn (41) in the main text.

We also derive how the probability density can be obtained from the intermediate scattering function

This is eqn (34) of the main text.

E. Harmonic approximation

The harmonic approximation of the Langevin equation, eqn (2), is given bywhere the first term on the right-hand side is the restoring force, and the second term is the random force of the Brownian motion. The equation is expressed by the shifted position, x̄ := x − L/2, so that the potential minimum is at the center. The Smoluchowski equation to solve for the propagator [Doublestruck P] := [Doublestruck P](x̄t|x̄₀) is given bywith the harmonic-well relaxation time τ = 1/DuQ₁². The solution is known as the Ornstein–Uhlenbeck process³⁶ and we find the propagatorwith V(t) = Dτ[1 − exp(−2t/τ)], and the stationary solution as the long-time limit of the propagator

We can readily compute the generalized ISF by using the definition of the main text, eqn (33), and by extending the integrals to infinity

Solving the integrals then yields the ISF

We further calculate the ratio

We also readily find the MSD

and time-dependent diffusion coefficient

D(t) = De^−t/τ. (79)

Acknowledgements

We acknowledge the use of AI (ChatGPT, Grammarly, LanguageTool, Perplexity) for its assistance with grammar checking, translations, and text enhancement. This research was funded in part by the Austrian Science Fund (FWF) 10.55776/I5257. For open access purposes, the author has applied a CC BY public copyright license to any author accepted manuscript version arising from this submission. MAES acknowledges funding by Deutsche Forschungsgemeinschaft (DFG) – Project number 459399860 for financial support.

References

1. W. Poon, P. Pusey and H. Lekkerkerker, Phys. World, 1996, 9, 27.
2. F. Q. Potiguar and R. Dickman, Phys. Rev. E:Stat., Nonlinear, Soft Matter Phys., 2007, 76, 031103.
3. C. Bechinger and E. Frey, Colloids on Patterned Substrates, John Wiley and Sons, Ltd, 2007, ch. 3, pp. 87–158.
4. D. Pérez Guerrero, J. Arauz-Lara, E. Sarmiento and G. Guerrero-Garcia, Front. Phys., 2021, 9, 635269.
5. X.-G. Ma, P.-Y. Lai, B. J. Ackerson and P. Tong, Soft Matter, 2015, 11, 1182–1196.
6. C. Dalle-Ferrier, M. Krüger, R. D. L. Hanes, S. Walta, M. C. Jenkins and S. U. Egelhaaf, Soft Matter, 2011, 7, 2064–2075.
7. R. Castaneda-Priego, Rev. Mex. Fis., 2021, 67, 050101.
8. A. Ashkin, Phys. Rev. Lett., 1970, 24, 156.
9. A. Ashkin, Science, 1980, 210, 1081–1088.
10. A. Ashkin, Proc. Natl. Acad. Sci. U. S. A., 1997, 94, 4853–4860.
11. A. Jonas and P. Zemanek, Electrophoresis, 2009, 29, 4813–4851.
12. D. Gao, W. Ding, M. Nieto-Vesperinas, X. Ding, M. Rahman, T. Zhang, C. Lim and C.-W. Qiu, Light: Sci. Appl., 2017, 6, e17039.
13. M. C. Jenkins and S. U. Egelhaaf, J. Phys.: Condens. Matter, 2008, 20, 404220.
14. R. Capellmann, A. Khisameeva, F. Platten and S. Egelhaaf, J. Chem. Phys., 2018, 148, 114903.
15. S. Lifson and J. L. Jackson, J. Chem. Phys., 1962, 36, 2410–2414.
16. R. Festa and E. G. d'Agliano, Phys. A, 1978, 90, 229–244.
17. R. L. Stratonovich, Topics in the theory of random noise, 1967.
18. P. Reimann, C. Van den Broeck, H. Linke, P. Hänggi, J. M. Rubi and A. Pérez-Madrid, Phys. Rev. E:Stat., Nonlinear, Soft Matter Phys., 2002, 65, 031104.
19. P. Reimann, C. Van den Broeck, H. Linke, P. Hänggi, J. M. Rubi and A. Pérez-Madrid, Phys. Rev. Lett., 2001, 87, 010602.
20. M. Evstigneev, O. Zvyagolskaya, S. Bleil, R. Eichhorn, C. Bechinger and P. Reimann, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2008, 77, 041107.
21. P. Fulde, L. Pietronero, W. R. Schneider and S. Strässler, Phys. Rev. Lett., 1975, 35, 1776–1779.
22. A. V. Straube and F. Höfling, J. Phys. A: Math. Theor., 2024, 57, 295003.
23. C. Kurzthaler, S. Leitmann and T. Franosch, Sci. Rep., 2016, 6, 36702.
24. H. Xu, X. Zheng and X. Shi, Phys. Rev. Res., 2024, 6, 023014.
25. B. J. Berne and R. Pecora, Dynamic light scattering: with applications to chemistry, biology, and physics, Courier Corporation, 2013.
26. J. K. Dhont, An introduction to dynamics of colloids, Elsevier, 1996.
27. R. Cerbino and V. Trappe, Phys. Rev. Lett., 2008, 100, 188102.
28. J. C. Crocker and D. G. Grier, J. Colloid Interface Sci., 1996, 179, 298–310.
29. C. Kurzthaler, C. Devailly, J. Arlt, T. Franosch, W. C. K. Poon, V. A. Martinez and A. T. Brown, Phys. Rev. Lett., 2018, 121, 078001.
30. C. Kurzthaler and T. Franosch, Soft Matter, 2017, 13, 6396–6406.
31. O. Chepizhko and T. Franosch, Phys. Rev. Lett., 2022, 129, 228003.
32. R. Rusch, O. Chepizhko and T. Franosch, Phys. Rev. E, 2024, 110, 054606.
33. D. B. Mayer, E. Sarmiento-Gómez, M. A. Escobedo-Sánchez, J. P. Segovia-Gutiérrez, C. Kurzthaler, S. U. Egelhaaf and T. Franosch, Phys. Rev. E, 2021, 104, 014605.
34. L. Angelani, Europhys. Lett., 2013, 102, 20004.
35. C. Kurzthaler, Y. Zhao, N. Zhou, J. Schwarz-Linek, C. Devailly, J. Arlt, J.-D. Huang, W. C. K. Poon, T. Franosch, J. Tailleur and V. A. Martinez, Phys. Rev. Lett., 2024, 132, 038302.
36. H. Risken, Fokker-Planck equation, Springer, 1996.
37. A. Lapolla, D. Hartich and A. Godec, Phys. Rev. Res., 2020, 2, 043084.
38. A. Asaklil, M. M. Mazroui and Y. Boughaleb, Phys. Condens. Matter, 1999, 10, 91–97.
39. A. Pattanayak, A. Shee, D. Chaudhuri and A. Chaudhuri, New J. Phys., 2024, 26, 083024.
40. G. Maruyama, Rendiconti del Circolo Matematico di Palermo, 1955, vol. 4, pp. 48–90.
41. D. Frenkel and B. Smit, Understanding Molecular Simulation: From Algorithms to Applications, Academic Press, San Diego, 2nd edn, 2002, vol. 2, pp. 2–5.
42. U. Siems, Dissertation, Universität Konstanz, Constance, 2017.
43. U. Siems and P. Nielaba, J. Phys.: Conf. Ser., 2018, 1012, 012011.
44. P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser and J. Bright, et al., Nat. Methods, 2020, 17, 261–272.
45. D. Blair and E. Dufresne, The Matlab Particle Tracking Code Repository, https://site.physics.georgetown.edu/matlab/, visited on 21/01/20.
46. X. Michalet, Phys. Rev. E:Stat., Nonlinear, Soft Matter Phys., 2010, 82, 041914.
47. J.-P. Hansen and I. R. McDonald, Theory of simple liquids: with applications to soft matter, Academic Press, 2013.
48. S. Lang, R. Schilling and T. Franosch, J. Stat. Mech.: Theory Exp., 2013, 2013, P12007.
49. S. Lang, R. Schilling, V. Krakoviack and T. Franosch, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2012, 86, 021502.
50. NIST Digital Library of Mathematical Functions, https://dlmf.nist.gov/, Release 1.2.3 of 2024-12-15, https://dlmf.nist.gov/, ed. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl and M. A. McClain.
51. J. Spiechowicz, I. G. Marchenko, P. Hänggi and J. Łuczka, Entropy, 2023, 25(2), 318.
52. A. Rahman, Phys. Rev., 1964, 136, A405–A411.
53. C. Dieball, Y. Mohebi Satalsari, A. B. Zuccolotto-Bernez, S. U. Egelhaaf, M. A. Escobedo-Sánchez and A. Godec, Soft Matter, 2025, 2170–2181.
54. M. C. Jenkins and S. U. Egelhaaf, J. Phys.: Condens. Matter, 2008, 20, 404220.
55. Q.-H. Wei, C. Bechinger, D. Rudhardt and P. Leiderer, Phys. Rev. Lett., 1998, 81, 2606–2609.

---

Footnotes

† We dedicate this work to the memory of our friend, mentor, and colleague Stefan U. Egelhaaf, with sincere gratitude for his contributions and dedication to the field of soft matter.

‡ R. R. and Y. S. contributed equally to this work.

---