Framework for Quasiperiodic Interfaces: Proximal Coincidence Point Set and Computation

This paper presents a unified theoretical and computational framework based on proximal coincidence point set (PCPS) theory that bridges mathematical quasiperiodicity with classical crystallography to rigorously model and resolve complex quasiperiodic interface structures, including non-crystallographic symmetries in various grain boundaries and phase boundaries.

The Gist

Imagine you are trying to fit two different patterns of floor tiles together at a seam.

If the tiles are simple squares and you rotate one by a perfect 90 degrees, they line up perfectly. This is like a periodic crystal, where everything repeats in a neat, predictable grid. Scientists have known how to model these "perfect" seams for a long time.

But what happens if you rotate the tiles by a weird angle, like 30 degrees, or if the tiles are slightly different sizes? The patterns won't line up perfectly. Instead of a neat grid, you get a messy, jumbled seam where the tiles never quite repeat the same way twice. In the world of atoms, this is called a quasiperiodic interface. It's a "perfectly imperfect" boundary that exists in real materials but has been incredibly hard to model because it doesn't follow simple repeating rules.

This paper introduces a new, unified way to understand and calculate these messy atomic seams. Here is the breakdown of their solution using simple analogies:

1. The Problem: The "Mismatched Puzzle"

Think of two large crowds of people (the atoms) standing in two different rooms (the bulk materials). They are walking toward a doorway (the interface).

• Old Models: Scientists used to assume the crowds were perfectly synchronized. If they weren't, they just faked it by forcing the patterns to repeat, which is like pretending a jigsaw puzzle fits when it clearly doesn't. This led to errors and missed important details.
• The Reality: In real life, the atoms don't care about our neat grids. They just want to get as close to their neighbors as possible without crashing into them. This creates a complex, non-repeating pattern that looks like a beautiful, chaotic mosaic.

2. The Solution: The "Proximal Coincidence" (PCPS)

The authors created a new theory called the Proximal Coincidence Point Set (PCPS).

• The Analogy: Imagine you have two flashlights shining on a wall. One is from the left room, one from the right. You don't need the light beams to hit the exact same spot to call them a match. You just need them to be close enough (within a small "tolerance" distance).
• The Magic: The PCPS theory says, "If an atom from the left side and an atom from the right side are close enough to shake hands, they form a special pair." By mapping all these "handshake" pairs, a hidden, beautiful pattern emerges. It turns out that even though the atoms are messy, their "handshake points" follow a strict mathematical rule based on cutting and projecting.

The "Cut-and-Project" Metaphor:
Imagine a 3D loaf of bread (the atomic structure). If you slice it straight down, you get a simple 2D pattern (a periodic crystal). But if you slice the loaf at a weird, irrational angle, the pattern on the slice looks random and never repeats.

• The authors realized that these messy 2D interfaces are actually just slices of a higher-dimensional, perfect 6D crystal.
• Their PCPS model is the mathematical "knife" that slices this 6D crystal at the perfect angle to reveal the 2D interface pattern we see in real life.

3. The Computer Simulation: The "Landau-Brazovskii" Engine

Knowing the theory is one thing; calculating it is another. To simulate this, they used a tool called the Landau-Brazovskii (LB) model.

• The Analogy: Think of the interface as a calm lake. The atoms are like ripples on the water. The LB model is a super-accurate weather forecast for these ripples. It calculates the energy of every possible ripple pattern and finds the one that is the most stable (the calmest water).
• The Innovation: Usually, computers struggle with quasiperiodic patterns because they try to force them into a repeating box (like trying to wrap a sphere in a square box). The authors used a Projection Method to solve this. Instead of forcing the pattern into a box, they projected the calculation onto a "donut" shape (a mathematical torus). This allows the computer to handle the infinite, non-repeating nature of the interface without losing accuracy.

4. What They Found (The Results)

Using this new framework, they simulated three types of atomic seams and found some cool things:

• The "Low-Angle" Seams: When the rotation is small, the atoms form a neat grid of "defects" (like a fence with holes). The size of these holes can be predicted perfectly by their math.
• The "High-Angle" Seams: When the rotation is large, the neat fence breaks down. The pattern becomes a true quasicrystal.
  • The Surprise: In some specific angles (like 30° and 45°), the atoms arrange themselves into 12-sided and 8-sided symmetries.
  • Why it matters: In normal crystals, you can't have 5, 8, or 12-sided symmetry (it's mathematically impossible to tile a floor with pentagons). But in these interfaces, the atoms do it! The authors proved that this happens because the "higher-dimensional" crystal they are slicing has a special symmetry that allows these impossible shapes to appear on the surface.
• The "Fibonacci" Pattern: In some tilt boundaries, the spacing between atoms follows the Fibonacci sequence (1, 1, 2, 3, 5, 8...), a pattern famous in nature (like sunflower seeds) but rarely seen in crystal interfaces.

Summary

This paper is like building a new pair of glasses for materials scientists.

• Before: They looked at messy atomic seams and saw only chaos or had to force them into simple, wrong models.
• Now: They have a unified framework that sees the hidden order in the chaos. They proved that these messy interfaces are actually slices of a perfect, higher-dimensional world.

This doesn't just explain how atoms stick together; it gives engineers a blueprint to design new materials with specific, exotic properties by intentionally creating these "perfectly imperfect" interfaces. It's a bridge between the messy reality of the physical world and the elegant perfection of mathematics.

Technical Summary

1. Problem Statement

Interfaces in crystalline materials (grain boundaries and phase boundaries) are critical to macroscopic properties like fracture resistance and irradiation response. While periodic interfaces are well-understood via the Coincidence Site Lattice (CSL) model, many real-world interfaces exhibit quasiperiodic (non-periodic but ordered) structures.

• Limitations of Existing Models:
  • CSL Model: Fails for general misorientations where lattices do not match periodically.
  • Geometric Approaches (Cut-and-Project, GCSN): While they describe quasiperiodicity, they are often purely geometric, lack a unified mathematical formulation for the transition from bulk to interface, and do not explicitly account for atomic mobility or physical perturbations.
  • Numerical Simulations: Current methods relying on periodic boundary conditions or free surfaces introduce Diophantine approximation errors, destroy long-range quasiperiodicity, and require prohibitive computational resources for large-scale observations.
• Core Challenge: There is a lack of a unified theoretical and computational framework that rigorously explains the origin of quasiperiodic order in interfaces while enabling accurate, large-scale numerical simulation.

2. Methodology

The authors propose a unified framework combining a new theoretical model with an efficient numerical solver.

A. Theoretical Model: Proximal Coincidence Point Set (PCPS)

The authors extend the classical CSL model to a Proximal Coincidence Point Set (PCPS) theory based on the Cut-and-Projection Scheme (CPS).

• Definition: Two atoms from adjoining bulk phases ($\Gamma_1, \Gamma_2$) are considered "proximal" if their midpoint lies near the interface and their mismatch is within a tolerance $\epsilon$.
• Mathematical Formulation:
  • The interface is modeled as a projection from a 6-dimensional lattice ($\mathbb{R}^6$) onto the 2D physical interface plane ($\mathbb{R}^2$).
  • The standard PCPS $\Lambda$ is defined as the projection of lattice points falling within a compact "window" $W$ in the internal space.
  • Physical Realism: The model incorporates atomic mobility by allowing arbitrary displacements of up to $\epsilon/2$ for points in the PCPS. This accounts for the fact that atoms at the interface are not rigidly fixed, bridging the gap between mathematical idealization and physical reality.
• Properties: The resulting point sets possess Finite Local Complexity (FLC), aperiodicity, and repetitivity, satisfying the rigorous mathematical definitions of quasicrystals.

B. Computational Framework: Conserved Landau-Brazovskii (LB) Model

To numerically realize the PCPS theory, the authors employ a Conserved Landau-Brazovskii (LB) free energy functional.

• Domain Partitioning: The simulation space is divided into two bulk regions ($x < -L, x > L$) and an interfacial transition region ($-L \le x \le L$).
• Spectral Discretization (Projection Method):
  • Instead of truncating the physical domain (which causes errors), the method uses the Fourier-Bohr spectrum derived from the PCPS theory.
  • The interfacial field is expanded using a projection method, mapping the 2D quasiperiodic plane onto a $d$-dimensional torus ($T^d$, where $3 \le d \le 6$).
  • This allows the quasiperiodic order to be resolved with uniform high accuracy across the entire interface without artificial boundary effects.
• Optimization: The energy minimization is solved using an accelerated proximal gradient method, with the bulk phases constrained by Dirichlet boundary conditions and the interface governed by the spectral properties of the PCPS.

3. Key Contributions

1. Unified PCPS Theory: A rigorous mathematical model that generalizes CSL to quasiperiodic interfaces, explicitly incorporating atomic mobility and physical perturbations.
2. Spectral-Computational Bridge: A novel implementation of the projection method within the LB model, enabling the simulation of arbitrary quasiperiodic interfaces with high precision and without Diophantine approximation errors.
3. Mechanistic Explanation of Symmetry: The framework explains why specific non-crystallographic symmetries (e.g., 8-fold, 12-fold) emerge in certain BCC interfaces while others (e.g., 5-fold, 10-fold) do not, based on the arithmetic properties of cyclotomic fields and the dimensionality of the projection.

4. Key Results

The framework was applied to three representative systems:

• BCC [100] Twist Grain Boundaries (GBs):

  • Low-Angle ($\theta \lesssim 15^\circ$): Revealed quasiperiodic dislocation networks with side lengths scaling as $\ell \approx a/\tan\theta$. The patterns are quasiperiodic but subtle, becoming apparent only at large scales.
  • High-Angle ($\theta \gtrsim 15^\circ$): Dislocation networks break down into distinct quasiperiodic patterns.
  • Quasicrystals: At specific angles ($\theta = 30^\circ$ and $45^\circ$), the model predicts and simulates 12-fold and 8-fold rotational symmetries (quasicrystals).
  • Symmetry Constraints: The model mathematically proves that 5-fold and 10-fold symmetries cannot emerge in BCC [100] twist GBs because the internal space projection does not support the required ring of integers ($\mathbb{Z}[\xi_n]$) for $n=5, 10$ given the 4D lattice constraint ($\phi(n)=4$).

• BCC [110] Tilt GBs:

  • Observed generalized Fibonacci sequences in atomic arrangements.
  • The interface patterns follow substitution rules (e.g., $L \to LSL, S \to L$), confirming the quasiperiodic nature of the structural units.

• BCC/FCC Phase Boundaries:

  • Demonstrated that interfaces between different crystal structures (BCC and FCC) are inherently quasiperiodic due to lattice constant mismatches.
  • Low-angle cases showed dislocation networks similar to twist GBs, while high-angle cases exhibited complex quasiperiodic ordering.

5. Significance

• Theoretical Advancement: Provides the first rigorous, unified explanation for the origin of quasiperiodic order in interfaces, moving beyond ad-hoc geometric constructions to a physically motivated model.
• Computational Efficiency: The projection method allows for the simulation of infinite quasiperiodic systems on compact domains, overcoming the limitations of traditional periodic boundary conditions.
• Materials Design: Offers a predictive tool for engineering interfaces with specific symmetries and properties. The insight that realizing $n$-fold quasicrystals requires specific internal space symmetries ($\phi(n)=4$ for 4D projections) guides the design of new quasicrystalline materials.
• General Applicability: The framework is not limited to grain boundaries but can be extended to any incommensurate system, including phase boundaries and heterogeneous interfaces.