🔬 mesoscale physics

Surface Phase-Field-Crystal-Helfrich model for out-of-plane deformations in thin crystalline sheets with lattice mismatch

This paper presents an extended Surface Phase-Field-Crystal-Helfrich model that incorporates spatially varying lattice spacing to simulate lattice mismatch in thin crystalline sheets, validating the approach against classical theories to demonstrate how localized compressive stresses drive out-of-plane deformations.

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

Published 2026-07-15

📖 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

Imagine a thin, flexible sheet of crystal, like a microscopic piece of origami made of atoms. Usually, when you push on such a sheet, it might crumple or bend. But what if the atoms inside the sheet don't quite agree on how far apart they should be? Some parts want to be squeezed tight, while others want to stretch out. This "lattice mismatch" creates a tug-of-war inside the material that can force the sheet to buckle and pop up into the third dimension, creating 3D shapes from a 2D surface.

In this study, researchers Emma Radice, Ingo Nitschke, Marco Salvalaglio, and Axel Voigt built a new digital playground to watch this drama unfold. They created a "multiscale description," which is a fancy way of saying they built a computer model that can see both the tiny, individual atoms and the big, smooth curves of the sheet at the same time. They combined two powerful ideas: the Phase-Field-Crystal (PFC) model, which tracks the atomic dance, and the Helfrich model, which describes how the sheet bends like a piece of paper.

The Secret Ingredient: A Shifting Grid
The big innovation here is that they allowed the "grid" of the crystal to change size in different spots. Think of a checkerboard where some squares are naturally smaller and some are naturally larger. In the real world, this happens in "heterostructures" (layers of different materials stuck together). In their model, they mimicked this by letting the equilibrium distance between atoms vary across the sheet. This creates "eigenstrain," which is like a built-in memory of stress that the material can't get rid of just by relaxing.

Testing the Model: The Bouncy Sheet and the Invisible Bubble
Before showing off their new toy, the team had to make sure it didn't break the laws of physics. They ran two specific tests:

The Squeeze Test: They took a perfect crystal sheet and squished it from one side (uniaxial compression). They compared their simulation results to the classic "Föppl–von Kármán" equations (the old-school math for bending sheets). The results matched perfectly, even when they added complex corrections for how the surface curves.
The Invisible Bubble Test: They simulated an "Eshelby inclusion," which is like dropping a tiny, invisible bubble of different material into the sheet. In classical physics, this bubble creates a uniform stress field inside it. Their simulation reproduced this exact pattern, proving their model could handle internal stress correctly without the sheet bending yet (when they set the bending rigidity, $\beta$ , to a high value of 20).

The Main Discovery: When Stress Makes the Sheet Pop
Once they were sure the model worked, they turned down the "stiffness" of the sheet (lowering $\beta$ to 0.5) to see what happens when the material is more flexible.

They found that when the sheet is flexible, it doesn't just sit there stressed; it relieves that stress by buckling up or down.

The Bulge: When they introduced a mismatched region (the "inclusion"), the sheet would bulge out.
The Defect Dance: They also introduced "dislocations," which are like missing or extra rows of atoms in the crystal. They found that these dislocations act like little anchors. If the stress around a dislocation is compressive (squeezing), the sheet forms a bulge right there.
The Interaction: When they combined a mismatched region with dislocations, the story got interesting. The dislocations provided an extra way for the material to release its stress. Instead of the whole sheet buckling into a huge, smooth wave, the presence of defects made the buckling pattern a bit more irregular and reduced the overall height of the bumps.

What They Didn't Find (And What They Didn't Say)
It is important to note what this paper doesn't do. They did not solve the problem for every possible shape or material. They specifically looked at circular inclusions and rectangular domains. They also noted that while their model works well for small-angle grain boundaries (where the crystal is rotated just a tiny bit, like 5 degrees), reproducing height modulation for large-angle boundaries (where the crystal is twisted wildly) is still a challenge that needs more work. They didn't claim to have solved the behavior of all polycrystalline metal foils, though they suggested their method could be applied to them in the future.

The Bottom Line
The researchers successfully demonstrated that by combining atomic detail with macroscopic bending physics, they can simulate how thin crystalline sheets deform. Their simulations show that localized stress from lattice mismatches and defects drives the sheet to buckle out of the plane. The presence of defects acts as a release valve, changing how the sheet bends and reducing the amplitude of the buckling compared to a perfect crystal.

In the end, this isn't just about math; it's about understanding how the tiny imperfections in a material's atomic structure dictate the big, visible shapes it takes. Whether it's a flexible electronic device or a biological membrane, knowing how these sheets pop and buckle helps us design better, more resilient materials. As the authors suggest, this approach opens the door to studying more complex, dynamic systems where the "mismatch" changes over time, but for now, they have firmly established the rules of the game for these flexible, mismatched crystal sheets.

Technical Summary: Surface Phase-Field-Crystal-Helfrich Model for Out-of-Plane Deformations

Problem Statement
Thin, flexible crystalline sheets exhibit complex mechanical behaviors driven by the coupling of in-plane elasticity and out-of-plane deformations. A comprehensive description of these systems must account for the crystalline structure, its defects (such as dislocations), and lattice mismatches found in heterostructures. While previous approaches have utilized atomistic descriptions, continuum modeling, or multiscale methods, there is a need for a framework that retains atomistic detail while capturing macroscopic bending effects and spatially varying equilibrium lattice parameters. Specifically, understanding how localized lattice eigenstrain (mimicking lattice mismatch) and defects interact to drive out-of-plane deformation remains a critical challenge.

Methodology
The authors extend the surface Phase-Field-Crystal (sPFC) model by coupling it with a macroscopic Helfrich energy description to account for bending effects. The core of the methodology involves:

Model Formulation: The system is defined by a free energy functional $F = F_L + F_P + F_b$ $F = F_{L} + F_{P} + F_{b}$ , where $F_L$ $F_{L}$ and $F_P$ $F_{P}$ represent the sPFC energy (characterizing the lattice structure and phase space) and $F_b$ $F_{b}$ represents the Helfrich bending energy.
- A key extension allows for a spatially varying equilibrium lattice parameter $q(x)$ , enabling the representation of localized lattice eigenstrain to mimic lattice mismatch in heterostructures.
- The model employs a height observer parameterization $X(t, x, y) = [x, y, h(t, x, y)]^T$ to describe surface deformations.
Mathematical Framework: The evolution of the system is governed by $(L^2, H^{-1})$ -surface gradient flows. To ensure energy dissipation and mass conservation ( $\frac{d}{dt}\int_S \psi dS = 0$ ) on evolving surfaces, the authors utilize the Truesdell time derivative and Truesdell gauge of surface independence. This approach addresses limitations in previous derivations regarding local expansion or contraction of the surface.
Numerical Implementation: The coupled equations for the density field $\psi$ and the height field $h$ are solved using a Fourier pseudo-spectral method with a semi-implicit time-stepping scheme. Nonlinear terms are evaluated in physical space and transformed via Fast Fourier Transforms (FFT).

Key Contributions

Extended Model: The development of a multiscale sPFC-Helfrich model that incorporates spatially varying equilibrium lattice spacing, allowing for the simulation of lattice mismatch and eigenstrain within a crystalline sheet.
Theoretical Refinement: The rigorous derivation of surface gradient flows using the Truesdell gauge, which introduces additional terms (e.g., $\psi \Delta_S h$ ) not present in prior simplified models.
Validation: The model is validated against two classical analytical benchmarks:
- Föppl–von Kármán (FvK) Equations: Under uniaxial compression, the model's predictions for buckling amplitudes match analytical solutions, including higher-order corrections.
- Eshelby's Inclusion Problem: For spherical inclusions with eigenstrain, the model reproduces the uniform stress field inside the inclusion and the decay outside, consistent with classical elasticity theory when bending rigidity is high.

Results
The validated framework was applied to study the interplay between localized lattice misfit, defects, and out-of-plane deformation:

Bending Rigidity Effects: Reducing the bending rigidity ( $\beta$ ) introduces a size dependence in the out-of-plane deformation. Unlike the classical Eshelby problem (or high $\beta$ cases), lower rigidity allows the sheet to bend, providing an additional relaxation mechanism that partially releases stress within the inclusion.
Defect-Driven Deformation: Inclusions rotated relative to the matrix (creating small-angle grain boundaries with isolated dislocations) induce out-of-plane bulges. These bulges form at dislocations where the hydrostatic stress is compressive. The direction of the bulge depends on initial perturbations, but the interaction between neighboring dislocations influences the final height profile.
Combined Effects: In simulations combining eigenstrain (lattice mismatch) and rotation-induced dislocations, the stress release from the inclusion superposes with the bulging at dislocations. The results indicate that defects act as an additional source of strain release, reducing the overall buckling amplitude compared to defect-free cases with the same eigenstrain. Consequently, the inclusion experiences smaller compressive stresses on average, while the interface with the surrounding matrix experiences larger tensile stresses.

Significance and Claims
The paper claims that the developed sPFC-Helfrich model provides a robust framework for predicting and controlling the morphology and mechanical response of 2D crystalline sheets. By successfully coupling atomistic detail with macroscopic bending, the model elucidates the competition between defect formation (plastic deformation) and out-of-plane deformation as relaxation mechanisms for internal stresses. The authors note that while the current work focuses on specific validation cases and simple heterostructures, the approach is applicable to larger systems and more complex structures, such as polycrystalline metal foils. The paper modestly concludes that further assessment is required to reproduce height modulation in the presence of large-angle grain boundaries and under dynamic conditions, which will be the subject of future work.

Technical Summary: Surface Phase-Field-Crystal-Helfrich Model for Out-of-Plane Deformations

More like this