HomeCompanion pagesIrrational Lattice

Irrational Lattice

A new primitive: deterministic, algebraic “colored noise” for lattices. Coordinates in ℚ(√D) deform ℤᵈ and snap back via a nearest-lattice map — provably aperiodic, spectrally tunable, integer-arithmetic clean.

Open the interactive lab Watch the video

Algebraic Colored Lattice Fields: A New Primitive for Structured Aperiodic Geometry

The Problem with Lattices

Lattices are too clean. The same properties that make them analytically tractable — periodicity, translational symmetry, uniform spacing — are precisely what makes them poor substrates for anything that is supposed to look physical, noisy, or generic. A pure lattice carries no information. Every point is equivalent to every other point under translation, and that uniformity is both a mathematical gift and a physical lie.

The obvious fix is to add noise. But random perturbations introduce a different problem: they are structureless in a different way. White noise has no correlations, no spectral color, no algebraic backbone. Smooth random fields are better, but they still require a probability measure, a choice of distribution, and they lose the exact arithmetic that makes lattices computationally tractable.

What we actually want is something in between: a perturbation that is deterministic, algebraically grounded, provably aperiodic, and spectrally tunable. Not random. Not periodic. Colored noise for fields.

This is the problem that motivates the construction described here.


The Core Idea: Algebraic Deformation and Nearest-Lattice Snapping

Start with a base lattice, say $L = \mathbb{Z}^d$. The goal is to define a family of point sets $X_\varepsilon$ that are:

  1. Infinitesimally close to $L$ as $\varepsilon \to 0$,
  2. Provably aperiodic — no nonzero translation $t$ satisfies $X_\varepsilon + t = X_\varepsilon$,
  3. Algebraically structured — coordinates live in a number field, not in a probability space,
  4. Spectrally tunable — the "color" of the noise can be shaped by design.

The mechanism has two steps.

Step one: algebraic deformation. For each lattice point $x \in \mathbb{Z}^d$, define a displacement vector

$$\eta(x) = a(x) + b(x)\sqrt{D}$$

where $a(x), b(x) \in \mathbb{Q}^d$ are rational vectors computed from $x$, and $D$ is a fixed square-free integer. The coordinates of $\eta(x)$ live in the quadratic number field $K = \mathbb{Q}(\sqrt{D})$. The deformed position is

$$T_\varepsilon(x) = x + \varepsilon , \eta(x).$$

Step two: nearest-lattice snapping. Project back to the original lattice using a nearest-point map under the $\ell^\infty$ norm:

$$S_\varepsilon(x) = \operatorname{argmin}{y \in \mathbb{Z}^d} |T\varepsilon(x) - y|_\infty.$$

The effective point set is then

$$X_\varepsilon = { S_\varepsilon(x) : x \in \mathbb{Z}^d }.$$

This is the construction. It is deterministic, algebraic, and produces a point set that is geometrically close to the original lattice while being provably free of translational symmetry.


Why the Irrational Component Is the Key

The critical constraint is that the irrational part of $\eta(x)$ must be bounded away from zero:

$$|b(x)| \geq c > 0 \quad \text{for all } x \in \mathbb{Z}^d.$$

This is not a technical nicety — it is the mechanism that kills periodicity.

Suppose there were a nonzero translation $t$ such that $X_\varepsilon + t = X_\varepsilon$. Tracing this back through the snapping map, it would require that for each $x \in \mathbb{Z}^d$ there exists $y \in \mathbb{Z}^d$ such that

$$t = (y - x) + \varepsilon(\eta(y) - \eta(x)).$$

The left side is a fixed vector. The right side has an irrational component $\varepsilon(b(y) - b(x))\sqrt{D}$. For this to hold for all $x$ simultaneously, the irrational parts would have to cancel — which requires $b(y) = b(x)$ for all relevant pairs, collapsing the irrational component to zero. But the lower bound $|b(x)| \geq c$ forbids this.

The minimum irrational component is a rigidity condition. It prevents the configuration from collapsing into any rationally periodic pattern, regardless of how the snapping map rearranges points.


Relationship to Known Structures

This construction does not live in a vacuum. It touches several well-studied mathematical objects, but it is not identical to any of them.

Cut-and-Project Sets and Quasicrystals

The closest relatives are quasicrystals constructed via the cut-and-project method. In that framework, one starts with a periodic lattice in a higher-dimensional space, applies an irrational linear projection to a lower-dimensional subspace, and retains only the points that fall within a bounded "window." The result is an aperiodic but highly ordered point set — Penrose tilings, Ammann-Beenker tilings, and their higher-dimensional generalizations all arise this way.

The algebraic structure is similar: quasicrystals built from $\mathbb{Q}(\sqrt{5})$ (the golden ratio field) have the same flavor of "irrational coordinates with algebraic backbone" as the construction here. The key difference is geometric. Cut-and-project sets are not infinitesimal perturbations of a lattice — they are genuinely different point sets, with different local configurations and different density. The construction here stays infinitesimally close to $\mathbb{Z}^d$ and snaps back to it, rather than projecting away from it.

Perturbed Lattices and Delone Sets

A perturbed lattice in the standard sense is a point set of the form ${x + \delta_x : x \in \mathbb{Z}^d}$ where the displacements $\delta_x$ are small and often random. If the perturbation is Lipschitz with small constant, the result is a Delone set: uniformly discrete and relatively dense, bi-Lipschitz equivalent to the original lattice.

The construction here is a perturbed lattice in this sense, but with two non-standard features: the perturbation is algebraically structured (not random), and the snapping step introduces a nonlinear quantization that is not present in the standard theory. The bi-Lipschitz equivalence still holds for small $\varepsilon$, but the mechanism is different.

Algebraic Number Fields in Geometry

The use of quadratic number fields to generate aperiodic structure is well-established in the theory of substitution tilings, Pisano sequences, and algebraic quasicrystals. What is unusual here is the role of the number field: rather than generating the point set directly, it generates the deformation field that is then snapped back to the lattice. The algebraic structure is one level removed from the geometry, which gives more flexibility in tuning the spectral properties of the resulting field.


The Field Perspective: Deterministic Colored Noise

The most important reframing is this: the construction is not primarily about the point set $X_\varepsilon$. It is about the field $\eta : \mathbb{Z}^d \to K^d$.

Think of $\eta$ as a noise field on the lattice. It is:

A concrete example: define

$$\eta(x) = \sum_{k=1}^{K} \left( r_k \sin(2\pi \alpha_k \cdot x) + s_k \sqrt{D} \cos(2\pi \beta_k \cdot x) \right)$$

where $\alpha_k, \beta_k$ are algebraically independent irrational vectors and $r_k, s_k$ are rational coefficients. By choosing the $\alpha_k$ and $\beta_k$, you control which spatial frequencies are present in the field. By choosing the $r_k$ and $s_k$, you control the amplitude at each frequency. The result is a deterministic field with a designed power spectrum — colored noise, but algebraic.

The perturbed lattice $X_\varepsilon$ is then just one geometric manifestation of this field. The field is the primary object; the geometry is derived from it.


Computational Structure: Why Quadratic Rationals Are Efficient

The choice of $K = \mathbb{Q}(\sqrt{D})$ is not just mathematically convenient — it is computationally optimal for large-scale parallel generation.

Each value of $\eta(x)$ is a pair of integers $(a(x), b(x))$ representing $a(x) + b(x)\sqrt{D}$. If $\eta(x)$ is built from linear forms in the lattice coordinates — say $a(x) = A \cdot x$ and $b(x) = B \cdot x$ for integer matrices $A, B$ — then computing $\eta(x)$ for every point in a large grid requires only integer dot products. There are no transcendental function calls, no floating-point instability, no per-site branching.

This is SIMD-friendly and GPU-friendly in the strongest sense: the computation at each lattice site is identical, involves only integer arithmetic, and has no data dependencies between sites. A large field of $\eta$ values can be generated in a single parallel pass over the lattice.

The $\sqrt{D}$ constant is global — it appears once, not per site. If you need actual floating-point geometry (for the snapping step or for visualization), you convert at the last moment. The algebraic structure lives in integer space throughout.

This is a significant practical advantage over smooth random fields, which require pseudorandom number generation (inherently sequential or at least stateful) and floating-point arithmetic throughout.


What the Visualization Will Look Like

If you plot the total perturbation magnitude $|\eta(x)|$ over a large region of $\mathbb{Z}^2$, the result will be deeply strange — in a structured way.

The most likely visual character is something close to the large-scale moiré patterns visible in Penrose tilings. This is not a coincidence; the mechanism is the same.

Why Penrose Moiré Appears

Penrose tilings arise from projecting a 5-dimensional cubic lattice along directions defined by the golden ratio $\phi = (1+\sqrt{5})/2$. When you zoom out on a Penrose tiling, you see ghostly bands and interference patterns — quasi-stripes that almost repeat but never quite do. These are the shadow of the 5-dimensional periodicity projected onto 2 dimensions. The "moiré" is the interference of incommensurate algebraic gratings.

The perturbation field $\eta(x)$ built from $\mathbb{Q}(\sqrt{D})$ has the same structure: it is the interference of incommensurate algebraic directions on a square lattice. The $\sqrt{D}$ component introduces a grating that is incommensurate with the integer lattice, and the interaction of these gratings produces large-scale interference patterns.

The differences from Penrose moiré are:

The visual result should feel like "Penrose moiré constrained to a square grid, with adjustable spectral weight." Long-range quasi-stripes, algebraic interference, no true periodicity, no randomness.

Fourier-Bohr Spectrum

The field $\eta$ has a pure point Fourier-Bohr spectrum — infinitely many "Bragg peaks" at positions determined by the algebraic structure of $K$, with intensities controlled by the coefficients $r_k, s_k$. When you blur or zoom out, these peaks interfere and produce the large-scale modulations visible in the plot.

This is the same mechanism that produces diffraction patterns in physical quasicrystals. The visualization of $|\eta(x)|$ is, in a precise sense, a real-space image of this spectrum.


Properties of the Construction

Collecting the key properties:

Infinitesimal perturbation. For small $\varepsilon$, $X_\varepsilon$ is within Hausdorff distance $O(\varepsilon)$ of $\mathbb{Z}^d$. As $\varepsilon \to 0$, $X_\varepsilon \to \mathbb{Z}^d$.

Bi-Lipschitz equivalence. If $\eta$ is Lipschitz with constant $L$ and $\varepsilon L < 1/2$, then

$$\frac{1}{2}|x - y| \leq |T_\varepsilon(x) - T_\varepsilon(y)| \leq \frac{3}{2}|x - y|$$

for all $x, y \in \mathbb{Z}^d$. The snapped set $X_\varepsilon$ is bi-Lipschitz equivalent to $\mathbb{Z}^d$.

Uniform density. The asymptotic density of $X_\varepsilon$ equals that of $\mathbb{Z}^d$.

Separation. The minimum inter-point distance in $X_\varepsilon$ is bounded below by a constant depending on $\varepsilon$ and the Lipschitz constant of $\eta$.

Provable aperiodicity. No nonzero translation $t$ satisfies $X_\varepsilon + t = X_\varepsilon$, provided the irrational component of $\eta(x)$ is bounded away from zero.

Deterministic. No randomness. The construction is fully explicit and reproducible.

Algebraic. All coordinates live in $\mathbb{Q}(\sqrt{D})$. Exact arithmetic is possible throughout.

Spectrally tunable. The power spectrum of $\eta$ can be shaped by choosing the algebraic frequencies and coefficients.


What This Is Not

It is worth being explicit about what this construction is not, to locate it precisely in the landscape.

It is not a quasicrystal. Quasicrystals are not infinitesimal perturbations of a lattice; they are genuinely different point sets with different local structure. This construction stays close to $\mathbb{Z}^d$ and snaps back to it.

It is not a random perturbed lattice. There is no probability measure, no randomness, no ensemble. The perturbation is fully deterministic.

It is not a substitution tiling. There is no inflation rule, no hierarchical decomposition, no tile alphabet.

It is not a smooth deformation. The snapping step introduces a discontinuous, nonlinear quantization that is not present in smooth deformation theory.

It is a new primitive: an algebraic colored lattice field, with a derived geometry given by the snapped perturbation.


Open Directions

Several natural questions follow from this construction.

Spectral classification. Given a target power spectrum (blue noise, pink noise, hyperuniform), what choice of algebraic frequencies and coefficients achieves it? This is an inverse problem in algebraic harmonic analysis.

Multi-scale stacking. Define

$$\eta(x) = \sum_k \varepsilon_k , \eta_k(x)$$

with different algebraic directions and scales. For $\sum_k \varepsilon_k$ small, this gives a hierarchical aperiodic lattice with multi-scale colored noise. What are the spectral and geometric properties of this stacked construction?

Higher-degree fields. The construction generalizes immediately to number fields of higher degree, $K = \mathbb{Q}(\alpha)$ with $\alpha$ algebraic of degree $n \geq 3$. Higher degree gives more complex quasi-periodic structure and richer interference patterns. What is the relationship between the degree of $K$ and the visual complexity of the resulting field?

Physical coupling. If $\eta(x)$ is interpreted as a physical field — a coherence-cost field, a local energy density, a disorder field — then the snapped geometry $X_\varepsilon$ is the substrate that physical processes see. How do dynamics constrained to $X_\varepsilon$ reflect the underlying algebraic noise field? This is the question of how algebraic structure in the substrate propagates into observable physics.

Formal naming. The object deserves a name. Candidates include: algebraic snap-lattice, infinitesimal algebraic quasicrystal, algebraically quantized aperiodic perturbed lattice. The name should reflect the three essential features: algebraic origin, infinitesimal amplitude, and nearest-lattice quantization.


Summary

The construction described here is a minimal, deterministic, algebraically grounded mechanism for producing a perturbed lattice that is:

The key insight is that the irrational component of a quadratic number field, when used to define a deformation field on a lattice and then snapped back via a nearest-lattice map, produces exactly the kind of structured irregularity that neither random noise nor smooth deformations can provide. The algebraic backbone gives you exact arithmetic and provable aperiodicity. The snapping gives you a point set that stays close to the original lattice. The spectral tuning gives you control over the "color" of the disorder.

This is not a known object in the literature. It is a hybrid of quasicrystal theory, perturbed lattice theory, algebraic number theory, and quantization theory — assembled into a single primitive that is simpler and more computationally tractable than any of its ancestors.