Turing patterns on non-fluctuating surfaces under mechanical stresses

This paper numerically investigates Turing patterns on non-fluctuating lattices using Finsler geometry to model mechanical stress, demonstrating that these static systems exhibit stress-induced pattern responses analogous to those observed on fluctuating membranes.

The Gist

Imagine you are looking at a zebra's stripes or the intricate swirls on a seashell. For a long time, scientists have believed these patterns are created by a chemical dance: two substances, an "activator" that encourages growth and an "inhibitor" that stops it, spreading out and reacting with each other. This is known as a Turing Pattern.

Usually, when scientists simulate this on a computer, they imagine the surface (like a fish's skin) is wiggly and flexible, like a rubber sheet. The chemicals move around, and the surface itself ripples and changes shape.

The Big Twist in This Paper
This paper asks a different question: What if the surface is completely rigid and still?

Think of a seashell or a zebra's skin not as a wiggly rubber sheet, but as a fixed grid of tiles, like a chessboard or a mosaic. The "tiles" (representing pigment cells) are stuck in place; they cannot move around. The only thing that changes is the color of the tile (the chemical concentration) and the direction of stress applied to the grid.

The researchers wanted to see if these "frozen" patterns could still react to being stretched or squeezed, just like the wiggly rubber sheets do.

The Secret Ingredient: The "Internal Compass"
To make this rigid grid behave like a real material, the scientists introduced a hidden variable they call an "Internal Degree of Freedom" (IDOF).

• The Analogy: Imagine every single tile on your chessboard has a tiny, invisible compass needle stuck to it.
• How it works: Even though the tile itself can't move, this compass needle can spin. When you stretch the whole board (like pulling a rubber band), these needles try to align themselves with the stretch.
• The Result: The direction these needles point changes how the chemicals (the activator and inhibitor) interact. If the needles point one way, the chemicals spread easily in that direction; if they point another way, they spread differently. This creates the "anisotropic" (direction-dependent) patterns we see in nature.

The Experiment: Stretching the Grid
The team ran computer simulations on three types of grids:

1. 2D Square Grid: Like a checkerboard.
2. 2D Triangular Grid: Like a honeycomb.
3. 3D Cube Grid: Like a block of dice.

They applied a "stretch" to these grids (making them longer in one direction and thinner in another) and watched what happened to the patterns.

What They Found

1. Rigid vs. Wiggly: Surprisingly, the patterns on the rigid, fixed grids behaved almost exactly like the patterns on the wiggly, flexible membranes studied in previous research.
2. The Stress Response: When they stretched the grid, the patterns reoriented themselves.
   • In one type of model, the stripes lined up parallel to the stretch (like lines drawn on a rubber band being pulled).
   • In another model, they lined up perpendicular to the stretch (like the rungs of a ladder being pulled apart).
3. The "Stress Relaxation" Discovery: This is the most fascinating part. The researchers calculated something called "entropy" (a measure of disorder or freedom). They found that at a specific point of stretching, the system reached a state of maximum entropy.
   • The Metaphor: Imagine you are holding a spring. You pull it tight, and it fights back. But at a certain point, the spring "relaxes" its internal tension. The paper suggests that even on a rigid grid where nothing moves, the internal "compass needles" can rearrange themselves to relieve stress, just like a flexible membrane would.

The Bottom Line
This paper proves that you don't need a wiggly, moving surface to create complex biological patterns. Even if the cells are stuck in a rigid grid (like pigment cells in a shell), the internal "direction" of the material is enough to make the patterns react to mechanical forces.

It's like saying you don't need a wiggly dance floor to create a dance; if the dancers (the chemicals) have a strong sense of direction (the compass needles), they can still create a beautiful, responsive pattern even if they are standing on a solid, unmoving floor.

What the Paper Does NOT Claim

• It does not claim this explains how to cure diseases.
• It does not claim this can be used to build new materials in a factory (yet).
• It strictly focuses on the mathematical and numerical proof that rigid grids can mimic the behavior of flexible biological membranes regarding pattern formation and stress response.

Technical Summary

1. Problem Statement

The paper addresses the formation of Turing Patterns (TPs) in biological systems (e.g., zebrafish stripes, sea shell patterns) and non-living materials. While traditional models assume TPs arise from reaction-diffusion (RD) processes on fluctuating membranes where cell movement and surface deformation play a role, recent biological evidence (e.g., in zebrafish pigment cells) suggests that cell movement may not be strictly necessary for pattern formation.

The core challenge is to model TPs on non-fluctuating (fixed) lattices where vertex positions are static, yet the system must still exhibit:

• Anisotropic diffusion: Direction-dependent diffusion coefficients ($D_x \neq D_y$).
• Mechanical response: The ability to react to external mechanical stresses (stretching/compression) and exhibit stress relaxation phenomena.
• Thermodynamic consistency: The ability to calculate entropy and free energy, which is traditionally difficult on fixed lattices because the partition function lacks integration over positional degrees of freedom.

2. Methodology

The authors extend a hybrid simulation model previously developed for fluctuating membranes to fixed lattices using Finsler Geometry (FG) modeling.

A. Model Framework

• Lattices: Simulations are performed on 2D square, 2D triangular, and 3D cubic lattices. Vertex positions ($\vec{r}$) are fixed (non-fluctuating), but the system includes an Internal Degree of Freedom (IDOF), denoted as $\vec{\tau}$, which represents the local direction of mechanical stress or polymer orientation.
• Reaction-Diffusion (RD) Equations: The system uses the FitzHugh-Nagumo (FN) equations for an activator ($u$) and an inhibitor ($v$). The Laplacian operator $\Delta$ is modified to depend on the IDOF $\vec{\tau}$, introducing anisotropy:
  $$\frac{\partial u}{\partial t} = D_u \Delta(\vec{\tau}) u + f(u,v)$$
  $$\frac{\partial v}{\partial t} = D_v \Delta(\vec{\tau}) v + g(u,v)$$
• Hamiltonian and Energy: The total energy includes:
  1. Gaussian Bond Potential (GBP): $H_1 = \sum \Gamma^G_{ij}(\vec{\tau}) \ell_{ij}^2$. On fixed lattices, bond lengths $\ell_{ij}$ are constant. However, the intensive part $\Gamma^G_{ij}$ depends on $\vec{\tau}$, allowing the system to store elastic energy and respond to deformation.
  2. IDOF Correlation: $H_\tau = \sum (1 - (\vec{\tau}_i \cdot \vec{\tau}_j)^2)$, promoting alignment of stress directions.
  3. RD Terms: $H_u$ and $H_v$ representing the reaction-diffusion energy.

B. Simulation Technique

• Hybrid Monte Carlo (MC):
  • RD Evolution: The variables $u$ and $v$ are updated using discrete time-stepping of the RD equations.
  • IDOF Update: The internal variables $\vec{\tau}$ are updated using Metropolis Monte Carlo simulations based on the Hamiltonian $H(\vec{r}, \vec{\tau})$.
• Stress and Entropy Calculation: A key methodological innovation is the derivation of a stress formula and entropy for fixed lattices. The authors demonstrate that by taking the limit of vanishing fluctuations ($\epsilon_R \to 0$) from a fluctuating model, one can define a well-posed stress tensor and entropy on a fixed lattice without positional integration.

C. Deformation Protocol

The lattices are subjected to uniaxial deformation (stretching along $x$, compressing along $y$) while maintaining constant area (2D) or volume (3D). The deformation ratio is denoted by $R$.

3. Key Contributions

1. Fixed Lattice Validity: The paper proves that Turing patterns can be successfully modeled on non-fluctuating lattices. This supports the hypothesis that pigment cell movement is not a prerequisite for pattern formation in organisms like zebrafish; the patterns can emerge from the interaction of chemical fields and internal stress directions on a static substrate.
2. Mechanical Anisotropy via IDOF: The study demonstrates that the Finsler geometry approach successfully induces anisotropic diffusion coefficients ($D_x \neq D_y$) solely through the orientation of the IDOF $\vec{\tau}$, even when the underlying lattice geometry is static.
3. Thermodynamic Consistency on Fixed Lattices: The authors successfully define and calculate entropy and stress on fixed lattices. They show that the stress formula derived for fluctuating membranes remains valid in the limit of zero fluctuation, allowing for the study of stress relaxation.
4. Dimensional Universality: The model is validated across 2D (square/triangular) and 3D (cubic) geometries, showing that the mechanical response of TPs is independent of dimensionality.

4. Key Results

• Pattern Orientation:

  • The direction of the Turing patterns aligns with the external mechanical stress.
  • Model 1 (Sine-type Finsler metric): Patterns align parallel to the tensile stress direction.
  • Model 2 (Cosine-type Finsler metric): Patterns align perpendicular to the tensile stress direction.
  • This behavior matches observations in fluctuating membrane models, confirming the robustness of the FG approach.

• Stress Relaxation and Entropy:

  • The system exhibits a peak in entropy density ($s$) at a specific coupling coefficient $\lambda_c \approx 0.88$.
  • When the lattice is deformed ($R=1.2$), the entropy decreases, and the stress increases.
  • Crucially, the system exhibits stress relaxation: when the external force is released, the system returns to a state of higher entropy, mimicking the behavior of soft biological materials. The peak position of the entropy remains invariant with respect to the deformation ratio $R$.

• Energy Localization:

  • Under deformation, energy localizes along specific axes depending on the model type. For example, in Model 1, energy concentrates along the direction of extension, making the perpendicular direction the "easy axis" for deformation.

• Comparison with Fluctuating Models:

  • The responses of the fixed lattice models (directional dependence of diffusion, stress-entropy relations) are nearly identical to those of fluctuating membrane models reported in previous literature. This validates the fixed lattice approximation as a computationally efficient and physically accurate alternative.

5. Significance

• Biological Insight: The findings provide strong theoretical support for the idea that biological patterns (like zebrafish stripes) can be generated by local chemical interactions and internal stress orientations without requiring the physical migration of pigment cells. This simplifies the biological requirements for pattern formation.
• Computational Efficiency: By removing the need to simulate vertex fluctuations (which requires complex integration over positions), the fixed lattice model significantly reduces computational cost while retaining the essential physics of anisotropic diffusion and mechanical response.
• Material Science Application: The model offers a framework for designing synthetic materials with tunable Turing patterns that respond to mechanical stress, applicable to soft robotics, self-organizing materials, and nanotechnology.
• Theoretical Advancement: The successful definition of entropy and stress on fixed lattices bridges a gap in statistical mechanics, showing that thermodynamic properties can be rigorously derived even when positional degrees of freedom are frozen, provided an internal degree of freedom (like stress direction) is active.