---
title: Aubry Transition in Incommensurate Systems
url: https://www.emergentmind.com/topics/aubry-transition
type: topic
---

# Aubry Transition in Incommensurate Systems

Searching arXiv for recent and foundational papers on the Aubry transition and closely related Aubry–André work.
The Aubry transition denotes a sharp change between an unpinned, sliding state and a pinned state in systems with competing incommensurate length scales. In its original formulation, it arises in the one-dimensional Frenkel–Kontorova model, where increasing substrate corrugation drives an incommensurate chain from a superlubric phase with zero static friction to a pinned phase with finite static friction, accompanied by breaking of analyticity of the hull function and the opening of a phason gap [1510.07585]. Subsequent work has shown that the concept extends well beyond the classical zero-temperature one-dimensional setting: in two-dimensional incommensurate monolayers it becomes a first-order structural and frictional transition that survives at finite temperature and terminates at a critical point [1705.06111], while in quasiperiodic tight-binding models the same name is also used for the Aubry–André localization transition between extended and localized eigenstates [2508.08255].

## 1. Classical formulation in the Frenkel–Kontorova picture

The canonical setting is the Frenkel–Kontorova model, a chain of particles connected by springs and placed in a sinusoidal substrate potential. A standard form is
\[
H = \sum_i \left[\frac{k}{2}(x_{i+1}-x_i-a_0)^2 + U_0 \cos\!\left(\frac{2\pi x_i}{a_s}\right) \right],
\]
with \(a_0\) the natural lattice spacing of the chain, \(a_s\) the substrate period, and \(U_0\) the corrugation amplitude [1802.09075]. The essential control parameter is the mismatch ratio
\[
\rho = \frac{a_0}{a_s},
\]
which distinguishes commensurate and incommensurate configurations [1802.09075].

For an incommensurate chain at \(T=0\), Aubry showed that increasing the substrate corrugation drives a transition between an unpinned and a pinned phase. Below a critical corrugation, the chain is in a translationally invariant sliding phase with no Peierls–Nabarro barriers and zero static friction; above it, the lattice potential overcomes the chain stiffness, the atoms reorganize toward lattice minima and avoid lattice maxima, and finite static friction appears [1510.07585]. In the infinite one-dimensional Frenkel–Kontorova chain, this Aubry transition is continuous and is conventionally described as a transition by breaking of analyticity of the hull function [1802.09075].

A standard diagnostic is the disorder parameter introduced in later formulations of the Aubry problem. In one dimension it measures whether particles can access the vicinity of substrate maxima. In the sliding phase, maxima remain accessible; in the pinned phase, a neighborhood of maxima becomes strictly unoccupied [1705.06111]. Closely related signatures are the onset of a phonon gap at \(q=0\), the appearance of finite static friction, and the loss of analyticity of the hull function [1802.09075]. In the infinite incommensurate limit, the hull function becomes a nowhere-analytic fractal staircase, or Cantorus, above the transition [1510.07585].

The same framework underlies the language of superlubricity. Below the transition, force cancellation across the incommensurate interface yields vanishing static friction. Above it, the interface becomes structurally pinned even though it remains incommensurate [2501.11061]. This association between analyticity breaking and frictional locking is central to later experimental and theoretical generalizations.

## 2. Finite chains, quantum effects, and experimental observation in trapped ions

A direct experimental realization of Aubry physics in a finite system was achieved with chains of laser-cooled \(^{174}\mathrm{Yb}^+\) ions in a linear Paul trap subject to a periodic optical lattice [1510.07585]. In this setting, the system is described by a generalized Frenkel–Kontorova–Tomlinson model with interatomic stiffness \(g\), external stiffness \(K=m\omega_0^2\), and optical lattice period \(a=185\,\mathrm{nm}\) [1510.07585]. The relevant mismatch condition for a finite chain is not irrationality in the strict infinite-chain sense, but cancellation of the net lattice force on the unperturbed chain,
\[
\sum_j \sin(2\pi x_{j,0}+\phi)=0 \quad \text{for any } \phi,
\]
which defines maximal mismatch and the superlubric case [1510.07585].

In that experiment, the Aubry transition was identified simultaneously through static and dynamical observables. Hysteresis loops in the ion position as a function of support position reveal both the opening of position gaps and the onset of stick–slip friction. The heights of the hysteresis loops reconstruct the static gaps in the allowed ion positions, while the separation between the slipping events equals twice the static friction force required to pull the chain over the Peierls–Nabarro barrier [1510.07585]. Below the critical lattice depth, the chain follows the global minimum smoothly and no hysteresis is present; above it, bistability, finite static friction, and stick–slip appear [1510.07585].

The same paper showed that the measured critical depths in finite mismatched chains coincide with the calculated values for the Aubry transition in an infinite golden-ratio chain with the corresponding values of \(g/K\) [1510.07585]. Finite size therefore rounds the transition into a symmetry-breaking crossover, but the onset of friction and analyticity breaking remains sharp enough to define an operational critical depth.

Quantum effects in the Aubry transition were later analyzed for trapped ions in an optical lattice using path-integral Monte Carlo [2008.12699]. There the system is a finite chain of ions with long-range Coulomb repulsion in a harmonic trap plus a periodic optical potential,
\[
\begin{split}
\mathcal{H}=&\sum_i
 \Bigg\{
 \frac{p_i^2}{2m} +\frac{1}{2}m\omega_0^2 \left(x_i-\frac{a}{2\pi}\Phi\right)^2 \\
 &\quad+ V \left[1+\cos\left(\frac{2\pi}{a} x_i\right)\right]
 \Bigg\} +\sum_{i<j}\frac{e^2}{\left | x_i-x_j\right|} .
\end{split}
\]
Quantum tunneling modifies the finite-chain Aubry crossover and provides new signatures that are robust against thermal and finite-size effects [2008.12699]. In particular, the Binder cumulant of the central ion’s position distribution and the total energy difference \(\Delta E(K)=E(K)-E(0)\) display sharp changes near the crossover, and the hull function develops gaps as pinning sets in [2008.12699]. This work emphasizes that the universality class of the quantum Aubry transition remains an open question [2008.12699].

## 3. Two-dimensional generalization: finite-temperature structural and frictional transition

A major extension of Aubry physics concerns two-dimensional incommensurate interfaces. Using a model inspired by colloid monolayers in an optical lattice, simulations showed that a sharp Aubry transition persists in two dimensions at finite temperature, but with properties fundamentally different from the one-dimensional case [1705.06111]. The model consists of classical particles interacting via a screened Coulomb potential,
\[
V(r) = \frac{Q}{r}\exp(-r/\lambda_D),
\]
on a triangular substrate potential
\[
W(\mathbf{r}) = -W_0\frac{2}{9}\left[\frac{3}{2}
    + 2\cos\frac{2\pi x}{a_l}\cos\frac{2\pi y}{\sqrt{3}a_l}
    + \cos\frac{4\pi y}{\sqrt{3}a_l}\right]
   \equiv W_0\, w(\mathbf{r}),
\]
with overdamped Langevin dynamics at temperature \(T\) [1705.06111].

In this system, incommensurability arises from both lattice spacing mismatch and rotational misalignment. The weak-coupling equilibrium configuration is characterized by a nonzero Novaco–McTague angle, and the mismatch produces a moiré superlattice made of nearly commensurate domains separated by a honeycomb network of soliton lines [1705.06111]. The two-dimensional Aubry transition occurs when increasing corrugation causes the local moiré domains to rotate from the Novaco–McTague angle toward near alignment with the substrate and to become locally commensurate [1705.06111].

Several observables establish the transition. A two-dimensional disorder parameter,
\[
\Psi = \frac{N_s}{N_p},
\]
counts the fraction of particles in a geometrically defined region around substrate maxima where the potential exceeds the saddle-point value. In the unpinned phase, particles occupy these high-energy regions with finite probability; in the pinned phase, those regions become effectively forbidden and \(\Psi\) drops sharply [1705.06111]. A local orientation angle,
\[
\theta_{\rm loc} = \left\langle\frac{1}{M}\sum_{\langle i,j\rangle}
   \bmod\!\left(\theta_{ij},\frac{\pi}{3}\right)\right\rangle,
\]
shows a simultaneous drop from the Novaco–McTague angle to nearly zero [1705.06111].

Unlike the \(1D\) Aubry transition, the two-dimensional transition is first order. Simulations show discontinuous jumps in \(\Psi\), \(\theta_{\rm loc}\), substrate energy, and interparticle energy, together with hysteresis and a finite coexistence region [1705.06111]. The free-energy balance is written as
\[
\Delta U + \Delta W - T\,\Delta S = \Delta F = 0,
\]
and the strongly oblique first-order line in the \(T\)–\(W_0\) plane implies that the unpinned phase has larger entropy [1705.06111]. This first-order Aubry line terminates at a finite-temperature critical point,
\[
W_{0c} \simeq 0.44,\quad T_c \simeq 0.11,
\]
marked by a susceptibility peak in the substrate-energy fluctuations [1705.06111].

This scenario was then realized experimentally in a two-dimensional colloidal monolayer on an optical lattice at room temperature [1802.09075]. The colloidal monolayer, with mismatch ratio \(\delta=a_l/a_c=0.84\), exhibits a sharp drop in mobility at a critical corrugation \(U_c \approx 34\,k_B T_0\), together with the onset of finite static friction [1802.09075]. Structural observables such as the fraction of tilted particles and a disorder parameter defined from occupation of the repulsive regions of the substrate unit cell show the same transition from tilted, superlubric domains to aligned, pinned domains [1802.09075]. The coexistence of pinned and unpinned regions, together with hysteresis in simulations, identifies the transition as first order in \(2D\) [1802.09075].

## 4. Friction, superlubricity, and load-induced locking in realistic interfaces

From its origin, the Aubry transition has been inseparable from friction. In the Frenkel–Kontorova picture, the sliding phase has vanishing static friction because no Peierls–Nabarro barriers exist, whereas the pinned phase exhibits a finite depinning threshold [1510.07585]. In the ion-chain experiment, the onset of friction and the opening of position gaps were shown to be the dynamic and static aspects, respectively, of the same Aubry transition [1510.07585]. In the two-dimensional colloidal monolayer, the Aubry line likewise separates a superlubric regime with \(F_s \approx 0\) from a pinned regime where static friction rises rapidly with corrugation, and the frictional upswing disappears as temperature approaches the finite-temperature critical point [1705.06111].

Recent work has pushed this connection into fully three-dimensional crystal interfaces. Simulations of Au(111) twist grain boundaries showed that the load-free moiré is smooth and superlubric at incommensurate twists, but increasing load induces a first-order structural transformation in which the highest-energy AA moiré nodes are removed and replaced by commensurate mini-domains [2501.11061]. This is described as an Aubry-type transition because the highest-energy local configurations disappear, in direct analogy with the elimination of particles at substrate maxima in the classical disorder-parameter picture [2501.11061].

The corresponding Landau description uses the mini-domain size \(Q\) as an Aubry order parameter and writes a free energy per moiré cell,
\[
\begin{aligned}
F(Q) &= F_1\,f(Q/a_0) + F_2\,[1-f(Q/a_0)] \\
F_1 &= 3 E_\mathrm{AA} \\
F_2 &= 3 E_\mathrm{vac}\, a\, Q^{-1} + \frac{\tau}{2}\, a^2 Q + \frac{G \theta^2}{2}\, a Q^{2}.
\end{aligned}
\]
Under load, this produces a first-order twist–load transition line terminating at a critical point [2501.11061]. Frictionally, the transformation causes a superlubric-locked transition with a huge friction jump, stick–slip, and irreversible plastic flow [2501.11061]. That study makes explicit that Aubry pinning in realistic crystal interfaces can become so strong that sliding relocates from the grain boundary to a nearby pristine crystal plane [2501.11061].

A distinct but conceptually related realization appears in an active mechanical setting. A snake-like robot undulating through a channel with a bichromatic array of hemispherical obstacles experiences a drag landscape approximating a one-dimensional Aubry–André potential. When the landscape is strictly periodic, the robot traverses the channel ballistically; when the landscape is sufficiently aperiodic, it becomes trapped and fails to exit [2507.03715]. The average travel distance, the mean-squared displacement exponent, and the distribution of travel distances reproduce the change from extended to localized transport familiar from the Aubry–André transition [2507.03715]. Although this is not the classical pinned–sliding Aubry transition of the Frenkel–Kontorova type, it illustrates the broader principle that incommensurate structure plus sufficient amplitude can suppress transport [2507.03715].

## 5. Aubry–André localization transition and its generalizations

In the spectral-theory and condensed-matter literature, “Aubry transition” often refers to the localization transition of the Aubry–André or Aubry–André–Harper model. The canonical Hermitian model,
\[
H_{\text{AAH}\psi)_n = t\left(\psi_{n+1}+\psi_{n-1}\right) + 2\lambda\cos\big(2\pi(\Phi n+\theta)\big)\,\psi_n,
\]
exhibits a sharp metal–insulator transition controlled by Aubry–André self-duality [2508.08255]. For irrational \(\Phi\), all eigenstates are extended for \(|\lambda/t|<1\), all are localized for \(|\lambda/t|>1\), and all are critical at \(|\lambda/t|=1\) [2508.08255]. Equivalent conventions in other papers place the critical point at \(\lambda_c=2J\) when the Hamiltonian is written with hopping \(J\) and quasiperiodic potential amplitude \(\lambda\) [1802.08859].

This localization transition has several generalizations in the supplied literature. Periodic modulation of the Aubry–André potential produces a Floquet localization–delocalization transition controlled by drive frequency and amplitude. In the noninteracting driven model,
\[
H(t)=H_0 + A\cos(\omega t)\sum_i\cos(2\pi\beta i+\phi)\ket{i}\bra{i},
\]
delocalization occurs at sufficiently low frequency and sufficiently large amplitude, with a low-frequency threshold
\[
A_c = \lambda - 2J
\]
and a characteristic frequency scale \(\hbar\omega_c \approx 2\lambda\) [1802.08859].

Non-Hermitian generalizations alter both spectral and dynamical critical behavior. For the non-Hermitian Aubry–André–Harper model with asymmetric hopping,
\[
i \frac{d \psi_n}{dt} = J_R \psi_{n+1} + J_L \psi_{n-1} + 2V \cos(2\pi\alpha n)\,\psi_n,
\]
the critical potential is \(V_c=J_L\), and the localization–delocalization transition becomes discontinuous not only in the diffusion exponent but also in the ballistic velocity [2102.09214]. In the \(\mathcal{PT}\)-symmetric non-Hermitian AAH model with complex onsite potential,
\[
E \psi_n = J(\psi_{n+1}+\psi_{n-1}) + V_0 e^{-2\pi i \alpha n} \psi_n,
\]
the metal–insulator transition at \(V_0=J\) coincides with \(\mathcal{PT}\)-symmetry breaking and a topological change in the complex spectrum; the localization length in the insulating phase remains energy-independent,
\[
\xi(E)=\frac{1}{\log(V_0/J)},
\]
exactly as in the Hermitian Aubry–André model [1908.03371].

Further extensions include disorder and hybrid universality. In the non-Hermitian disordered Aubry–André model with nonreciprocal hopping and random onsite disorder, both the quasiperiodic control parameter \(\delta\) and the disorder strength \(\Delta\) are relevant directions at the non-Hermitian AA critical point, leading to a new universality class and a hybrid scaling law in the overlap between Aubry–André and Anderson critical regions [2405.15220]. The non-Hermitian interpolating Aubry–André–Fibonacci chain exhibits a cascade of delocalization transitions as the potential is deformed from the Aubry–André limit toward the Fibonacci limit, with self-similar critical modes and only a few plateaux surviving under nonreciprocal hopping [2104.06035].

The many-particle ground-state version of the Aubry–André transition also acquires a richer arithmetic structure. For the interacting Aubry–André model, the critical behavior depends on a Diophantine equation relating the incommensurate frequency and the filling fraction,
\[
Q_k M_k - P_k N_k = \pm N_F,
\]
which generalizes the dependence of the single-particle critical properties on the continued-fraction expansion of the incommensurate frequency [2001.04954]. Numerical evidence suggests that nearest-neighbor interactions may be irrelevant at at least some of these critical points, so the Diophantine classification can survive in the interacting case [2001.04954].

## 6. Open directions, reinterpretations, and broader significance

Several recent papers recast the Aubry transition by changing the microscopic setting while preserving the central competition between incommensurate structure and coupling. One such direction concerns the shape of the substrate potential itself. In a modified Frenkel–Kontorova model with periodic potential
\[
V_\alpha(x) = \frac{\cosh\alpha\,\cos(2\pi x) - 1}{\cosh\alpha - \cos(2\pi x)},
\]
the Aubry transition can occur with very small structural distortions if the local substrate force becomes very strong in narrow regions while the potential remains weak elsewhere [2405.01878]. In that model, the critical coupling \(K_c(\alpha)\) tends to zero as \(\alpha\to 0\), the phason gap opens in the pinned phase, yet the high-energy part of the phonon spectrum remains close to that of an undistorted chain [2405.01878]. This suggests that pinned incommensurate phases with small observable distortions need not contradict an Aubry mechanism [2405.01878].

Another active direction concerns dissipation and nonequilibrium open-system transport. In the Aubry–André–Harper chain coupled to a non-Markovian bath, the role of bath memory depends qualitatively on the side of the localization transition [2605.22966]. In the extended phase, memory reshapes the dynamical generator and produces transport patterns that cannot be reduced to a simple rescaling of time, whereas in the localized phase the bath activates motion between localized states and memory mainly renormalizes the dynamical timescales [2605.22966]. This suggests that localization acts as a filter of non-Markovian effects [2605.22966].

Across these formulations, several points recur. First, the Aubry transition is not restricted to a single model, but identifies a class of incommensurate pinning or localization phenomena controlled by the competition of length scales and couplings. Second, dimensionality matters: the one-dimensional pinned–sliding transition is continuous at \(T=0\), whereas realistic two-dimensional interfaces can exhibit first-order, finite-temperature transitions with coexistence and critical endpoints [1705.06111]. Third, the transition is inseparable from transport: it manifests as the onset of static friction in tribology and as localization in quasiperiodic spectral problems. Fourth, structural degrees of freedom absent in the original Frenkel–Kontorova chain—rotation of moiré domains, twist elasticity, superconducting pairing, non-Hermiticity, or bath memory—do not abolish Aubry physics but can profoundly alter its order, universality, and observables.

In that sense, the Aubry transition remains both a specific phenomenon and a unifying concept. In one branch of the literature it describes the passage from superlubric sliding to structural pinning at incommensurate interfaces [1510.07585]; in another it denotes the self-duality-driven localization transition of quasiperiodic quantum lattices [2508.08255]. The supplied work shows that both usages continue to expand, into two-dimensional colloids at ambient temperature [1802.09075], three-dimensional crystal grain boundaries under load [2501.11061], active mechanical transport [2507.03715], and non-Hermitian, interacting, driven, and dissipative quasiperiodic systems [2405.15220].

Source: https://www.emergentmind.com/topics/aubry-transition