🔬 mesoscale physics

Wrinkling dynamics of thin crystalline sheets and its dependency on topological defects

Using a Surface Phase-Field-Crystal-Helfrich model, this study reveals that while crystalline defects do not alter the scaling laws of wrinkling coarsening, they significantly influence the onset, wavelength, and pattern irregularity of wrinkles by promoting a higher density of wrinklons and modifying long-term morphological evolution.

Original authors: Emma Radice, Ingo Nitschke, Marco Salvalaglio, Axel Voigt

Published 2026-10-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Emma Radice, Ingo Nitschke, Marco Salvalaglio, Axel Voigt

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Thin sheets of material, whether a piece of paper or a microscopic film, have a peculiar habit when squeezed. If they are free to move, they simply buckle and fold. But if they are stuck to a soft, sticky surface, they cannot fold over themselves. Instead, they buckle upward, forming a landscape of parallel ridges and valleys known as wrinkles. For decades, scientists have understood that these patterns are not random; they follow the laws of elasticity, growing in size and smoothing out over time until they reach a stable, regular state. This behavior is not just a curiosity of physics; it is the foundation for technologies ranging from flexible electronics that can bend without breaking to surfaces designed to harvest light more efficiently. However, this classical understanding assumes the material is perfectly uniform, like a flawless sheet of glass. Real materials, especially those at the microscopic scale, are rarely perfect. They contain tiny imperfections in their atomic structure, such as missing atoms or misaligned rows of crystals, which act as topological defects. The question remains: how do these microscopic flaws change the way a sheet wrinkles when it is compressed?

A team of researchers at the Technical University of Dresden set out to answer this by simulating the behavior of thin crystalline sheets under pressure. They did not use physical materials in a lab, but rather a sophisticated computer model that treats the sheet as a collection of atoms arranged in a hexagonal pattern, similar to the structure of graphite. This model allowed them to watch, in slow motion, how the sheet deforms when squeezed from one side. Crucially, their simulation could track two different types of imperfections simultaneously. The first type exists within the crystal lattice itself, such as dislocations where the atomic rows are misaligned, and grain boundaries where different crystal orientations meet. The second type appears in the wrinkles themselves, known as wrinklons, which are defects in the pattern of the ridges, much like a kink in a line of marching soldiers. By introducing varying numbers of these atomic defects into their virtual sheets, the researchers could observe how the presence of a flawed crystal structure influenced the formation and evolution of the surface wrinkles.

The results revealed a nuanced story about the relationship between the microscopic and the macroscopic. When the researchers squeezed a perfect, defect-free crystal, it developed a neat, regular pattern of wrinkles that grew larger and more uniform over time, exactly as classical physics predicts. The size of these wrinkles increased at a specific, predictable rate. When they introduced atomic defects, such as dislocations, into the sheet, the story changed in specific ways, but not in the way one might expect. The initial size of the wrinkles was slightly different, and the patterns that emerged were far less orderly. Instead of neat, parallel lines, the surface became a chaotic landscape of irregular ridges. The atomic defects acted as anchors, pinning the wrinkles in place and preventing them from smoothing out into a perfect pattern. This led to a much higher density of wrinklons, the kinks and breaks in the wrinkle lines, which persisted for much longer than they did in the perfect sheets.

Despite this visible chaos, the researchers found that the fundamental rule governing how the wrinkles grew over time remained unchanged. Whether the sheet was perfect or riddled with hundreds of atomic defects, the average size of the wrinkles still increased at the same mathematical rate. The defects did not break the law of how the wrinkles coarsen; they simply altered the scenery in which that law played out. The presence of defects delayed the moment when the wrinkles began to grow significantly, and it made the process of releasing the built-up stress much less efficient. In sheets with many defects, the stress remained trapped in the material for longer, and the final pattern was a messy, irregular texture rather than a clean, functional surface. The study suggests that while the underlying physics of how wrinkles grow is robust, the quality and regularity of the final pattern are entirely dependent on the microscopic health of the material.

This distinction is vital for anyone looking to engineer these materials for real-world use. If the goal is to create a flexible sensor or a light-harvesting surface, the regularity of the pattern is often just as important as the fact that it wrinkles at all. The research indicates that controlling the atomic perfection of a film is a critical step in ensuring the resulting surface behaves as intended. A sheet that is slightly flawed might still wrinkle, but it will do so in a way that is unpredictable and potentially less useful. The study confirms that while the broad strokes of wrinkling are universal, the fine details are dictated by the hidden architecture of the crystal. By understanding how these tiny atomic flaws interact with the large-scale bending of the material, scientists can better predict how thin films will behave, moving beyond idealized theories to account for the messy reality of the materials we actually build with.

Technical Summary: Wrinkling Patterns and Topological Defect Dynamics in Crystalline Films

Problem Statement
Thin crystalline sheets under uniaxial compression typically undergo a stress-driven instability, transitioning from a flat state to a wrinkled state. While classical elasticity theory predicts regular sinusoidal wrinkles with a specific equilibrium wavelength and coarsening dynamics, the emergence of these patterns is sensitive to initial perturbations. In real materials, the presence of topological defects within the crystalline lattice (such as dislocations and grain boundaries) introduces complex interactions between in-plane elastic relaxation and out-of-plane deformations. The specific impact of these crystalline defects on the onset of wrinkling, the evolution of wrinkle patterns (including the formation of "wrinklons"), and the long-term coarsening regime remains poorly understood within a unified framework that resolves atomic-scale structures.

Methodology
The authors employ a Surface Phase-Field-Crystal-Helfrich (SPFCH) model to investigate the wrinkling dynamics of thin crystalline sheets. This approach is chosen over classical elasticity theories (e.g., Föppl–von Kármán equations) because the PFC model operates on atomic length and diffusive time scales, naturally incorporating crystalline order and elasticity, while allowing for the explicit initialization of topological defects.

Model Formulation: The SPFCH model couples the evolution of the order parameter $\psi$ (representing atomic density and crystalline order) and the surface geometry $S$ (representing out-of-plane deformations). It utilizes a $(L^2, H^{-1})$ surface gradient flow to minimize a free energy functional $F = F_L + F_P + F_b$ , where $F_L$ is the Swift-Hohenberg elastic energy, $F_P$ is the thermotropic part, and $F_b$ is the Helfrich bending energy.
Dynamics: The system is modeled with finite mobility coefficients ( $M_\psi$ and $M_X$ ). The mobility $M_X$ implicitly accounts for viscous dissipation in the surrounding medium (substrate), allowing the model to capture overdamped dynamics without explicitly simulating a viscous substrate.
Simulation Setup: Numerical simulations are performed on a rectangular domain with periodic boundary conditions using a Fourier pseudo-spectral method on GPUs. The study compares four distinct initial conditions:
1. A perfect hexagonal crystal.
2. Crystals with varying densities of isolated dislocations ( $N_\psi = 8, 16, 60, 100$ ), which are manually inserted by initializing the density field and imposing a singularity in the phase to create specific Burgers vectors.
3. A polycrystalline sheet with small- and large-angle grain boundaries, generated by placing multiple crystal seeds with random orientations.
Analysis: The authors track the evolution of mean wavelength ( $\bar{\lambda}$ ), root-mean-square amplitude ( $W_{RMS}$ ), Gaussian curvature ( $K_{RMS}$ ), von Mises stress ( $\bar{\sigma}_{VM}$ ), the number of crystalline defects ( $N_\psi$ ), and the number of topological defects in the wrinkle pattern known as wrinklons ( $N_h$ ). Wrinklons are identified by demodulating the height field and counting phase windings.

Key Contributions and Results
The study reveals a complex coupling between atomic-scale crystalline defects and macroscopic wrinkling patterns:

Onset Wavelength: The numerical simulations validate the theoretical prediction for the onset wavelength in defect-free crystals (Eq. 3.6). The presence of crystalline defects causes minor quantitative deviations in the initial wrinkling wavelength but does not fundamentally alter the theoretical scaling.
Coarsening Scaling: A primary finding is that crystalline defects do not affect the scaling law of the coarsening regime. Across all systems (perfect crystals, dislocated crystals, and polycrystals), the wavelength and roughness evolve according to a $t^{1/4}$ power law, consistent with previous studies on viscous substrates.
Pattern Irregularity and Wrinklons: While the scaling law remains universal, the morphology of the wrinkles is significantly altered by defects.
- Higher densities of crystalline defects (dislocations and grain boundaries) lead to the nucleation of a higher density of wrinklons (topological defects in the wrinkle pattern).
- These systems exhibit more irregular, persistent wrinkle patterns that deviate from the ideal alignment perpendicular to the compression axis.
- In polycrystalline systems, the onset of the coarsening regime is delayed, and the growth of wrinkle amplitude is hindered compared to defect-free systems.
Stress Relaxation Mechanisms: The interplay between in-plane and out-of-plane relaxation is critical.
- In defect-free crystals, stress relaxation is dominated by the formation of perfectly aligned wrinkles, leading to a significant drop in von Mises stress at late times.
- In defective crystals, stress relaxation is less efficient. The presence of crystalline defects and the resulting high density of wrinklons constrain out-of-plane deformation. Consequently, residual stress remains higher in defective systems at long timescales.
- Crystalline defects themselves migrate and annihilate (e.g., dislocation dipoles reacting), contributing to initial plastic relaxation, but the remaining defects ultimately determine the residual stress state.

Significance and Claims
The paper claims that while classical elasticity theory correctly predicts the universal scaling laws for wrinkling coarsening, it fails to capture the morphological details and long-term evolution dictated by the material's crystalline microstructure.

The authors assert that crystalline defects primarily control the morphology and relaxation efficiency of wrinkled thin sheets rather than the underlying coarsening laws. Specifically, defects promote irregular patterns and hinder the system's ability to reach a low-stress, perfectly ordered state. The authors suggest that for applications relying on the functional properties of wrinkled surfaces (such as flexible electronics, strain sensors, or light harvesting), the specific arrangement of crystalline defects and the resulting density of wrinklons are crucial factors that determine mechanical and functional performance. The work provides a foundation for understanding how atomic-scale processes couple with macroscopic deformation, offering a basis for potentially tailoring the morphological evolution of crystalline thin films through microstructural control.

Technical Summary: Wrinkling Patterns and Topological Defect Dynamics in Crystalline Films

More like this