Emergence of rigid Polycrystals from atomistic Systems with general Interactions

Leonard Kreutz; Timo Ziereis

arxiv: 2604.19239 · v1 · submitted 2026-04-21 · ❄️ cond-mat.mes-hall · math-ph· math.MP

Emergence of rigid Polycrystals from atomistic Systems with general Interactions

Leonard Kreutz , Timo Ziereis This is my paper

Pith reviewed 2026-05-10 02:32 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall math-phmath.MP

keywords polycrystalsdiscrete-to-continuum limitGamma-convergencegrain boundariesatomistic modelsorientation fieldsrigid interactions

0 comments

The pith

Atomistic systems with rigid interactions converge to polycrystal models with energy only on grain boundaries.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that particle systems favoring local isometry to a lattice develop polycrystalline structures in the large-scale limit. It uses Gamma-convergence to prove that the atomistic energy approaches a continuum functional where energy is supported solely on interfaces separating regions of constant orientation. These interfaces represent grain boundaries whose energy depends on the two orientations and the interface normal. The rigidity of the interactions prevents low-energy transition layers between solids, causing the interface energy to equal twice the solid-vacuum energy.

Core claim

What carries the argument

The Γ-limit of the atomistic energy, a functional on piecewise constant orientation fields whose density on grain boundaries equals twice the solid-vacuum transition energy.

Load-bearing premise

The rigid interactions are such that any interpolating boundary layer between two different solid orientations costs at least as much energy as two separate solid-vacuum transitions.

What would settle it

Numerical minimization of the atomistic energy for two grains of different orientations showing a thin interpolating layer with total energy less than twice a single grain-vacuum interface would disprove the energy decomposition.

Figures

Figures reproduced from arXiv: 2604.19239 by Leonard Kreutz, Timo Ziereis.

**Figure 1.** Figure 1: Schematic Illustration of a competitor X for the cell problem on Q ν ρ in the definition of φ. Also, this illustrates a configuration X ∈ Adm(z+,z−) ε,λ (Q ν ρ). Note that the scaling factor between the lattice spacing and the boundary thickness does not match the definition. The limiting functional E : P C(R d ; Z) → [0, +∞) is defined as E(u) := ˆ Ju φ(u +, u−, νu) dHd−1 (2.16) . In view of (2.11), funct… view at source ↗

**Figure 2.** Figure 2: This is a schematic picture of a transition from two constant values considered in the construction of the upper bound for d = 2. The dark and light grey regions are regions where the configuration coincides with the lattice and vacuum, respectively. Proof of Theorem 2.6(ii). We split the proof into several steps. Step 1: Energy representation via φvac. Let u ∈ PC(R d ; Z). We then have E(u) = X∞ i=1 ˆ ∂∗G… view at source ↗

**Figure 3.** Figure 3: Illustration of the different regions in the definition of Yε. The dark grey region is D ε ∪ ∂ − λεQ ν , the grey region is (R ν 1,δ \ Q ν rkε−1 ) \ (∂ + λεQ ν ∪ ∂ − λεQ ν ), the light grey region is S ε kε and the white region is (R+ δ,ε \ Q ν rkε )∪∂ + λεQ ν . The dashed lines enclose R ν 1,δ. Energy estimate on Aε 1 : We claim that there exists a C > 0 such that Eε(Yε, Aε 1 (5.18) ) ≤ Cδ . By (5.14) we … view at source ↗

**Figure 4.** Figure 4: The definition of XT Note that for all j ∈ Zν S,T and all x ∈ XT ∩ Qν S (x ± j,S) we have XT ∩ Brint (x) = X ± j,S ∩ Brint (x), whenever x /∈ A ± j := (Qν S (x ± j,S) \ Qν S−2rint (x ± j,S)) ∩ {y | |⟨y − x ± j,S, ν⟩| ≤ 2rint}. However, we have #(XT ∩ A ± j ) ≤ Ld ((A ± j )1) ≤ CSd−2 by Lemma 4.1(v) and therefore, using (8.4) and (8.5) respectively, we obtain E1(XT , Qν S (x + j,S)) ≤ E1(X + j,S, Qν S (x + … view at source ↗

read the original abstract

We investigate the formation of polycrystalline structures in a class of particle systems. The atomistic energy is modeled as a sum of particle energies that favor atoms being locally isometric to a reference lattice. The discrete frame invariant energy allows for particle configurations in which no underlying lattice is assumed a priori. We prove a discrete-to-continuum limit for configurations with finite surface-energy scaling by means of $\Gamma$-convergence. The resulting continuum theory is described by piecewise constant fields encoding the local orientation of the configuration. The limiting energy is concentrated on grain boundaries, corresponding to the interfaces between regions where the microscopic configuration has constant orientation. The associated energy density depends on the orientations of the two grains as well as on the normal to the interface. Due to our assumptions on the rigid interactions, solid-solid phase transitions with interpolating boundary layers are not energetically favorable; the energy density therefore decomposes into twice the energy density for solid-vacuum transitions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a Gamma-convergence from frame-invariant atomistic rigid energies to a continuum polycrystal model whose grain-boundary energy reduces to twice the solid-vacuum interface energy under the stated interaction assumptions.

read the letter

The core contribution is a discrete-to-continuum Gamma-convergence result for particle systems whose energy favors local isometry to some reference lattice without fixing that lattice in advance. The limit objects are piecewise-constant orientation fields, and the energy concentrates on the interfaces between regions of different orientation. The authors show that, for their class of interactions, any solid-solid transition that uses an interpolating layer costs strictly more than two separate solid-vacuum surfaces, so the limiting density factors into twice the vacuum-interface density. That explicit decomposition is new for this setting and removes the need to track orientation-dependent solid-solid costs separately. The proof strategy follows the usual compactness-plus-liminf-plus-limsup route for finite-surface-energy sequences, which is clean on paper. The frame-invariance and lack of a pre-assumed lattice are genuine extensions beyond many earlier lattice-based Gamma-convergence results in the same program. The main limitation is that everything rests on the interaction assumptions that make interpolating solid-solid layers expensive. If those assumptions are narrow, the result applies only to a restricted family of potentials and the claimed decomposition may not survive for more general rigid interactions. The abstract does not give concrete examples of admissible potentials or show how much room there is before the exclusion fails, so the practical scope remains unclear. The handling of finite-surface-energy configurations also needs to be checked in detail to confirm that no other low-energy transition mechanisms sneak through. This is a technical variational-analysis paper aimed at readers who already work on discrete-to-continuum limits for polycrystals or grain boundaries. Someone looking for a rigorous justification of a simple continuum model from particles will find the theorem useful; a reader outside that niche will not. It is worth sending to referees because the statement is new, the method is standard but applied cleanly, and the only real question is whether the interaction hypotheses are stated sharply enough and whether the proof closes all the gaps.

Referee Report

2 major / 1 minor

Summary. The manuscript proves a discrete-to-continuum Γ-convergence result for atomistic energies that favor local isometry to a reference lattice without assuming an a priori lattice structure. For sequences with finite surface-energy scaling, the limit consists of piecewise-constant orientation fields, with the continuum energy concentrated on grain boundaries whose density depends on the two orientations and the interface normal. Under the stated assumptions on the rigid interactions, solid-solid transitions with interpolating boundary layers are energetically unfavorable, so the limiting energy density equals twice the solid-vacuum interface energy.

Significance. If the assumptions on the interaction potentials hold, the result supplies a rigorous variational justification for continuum polycrystal models arising from general discrete systems, extending lattice-based analyses to frame-invariant settings. The explicit energy decomposition and the identification of the limit objects as piecewise-constant fields constitute a clear mathematical contribution that could be tested against atomistic simulations.

major comments (2)

[Abstract; statement of main theorem (likely Theorem 1.1 or §3)] The central claim that the limiting energy density decomposes into twice the solid-vacuum transition energy rests on the assertion (stated in the abstract and presumably formalized in the hypotheses of the main Γ-convergence theorem) that solid-solid phase transitions with interpolating boundary layers are strictly higher in energy than two separate solid-vacuum interfaces. The precise conditions on the interaction potential that enforce this exclusion are not exhibited with a counter-example or a minimal set of inequalities that would allow verification for a concrete potential; without this, it is unclear whether the liminf inequality produces the claimed continuum energy for all admissible sequences.
[Proof of Γ-convergence (likely §4 or §5)] In the proof of the liminf inequality, the argument that any interpolating layer can be replaced by two solid-vacuum interfaces without increasing energy must be checked against the finite-surface-energy scaling; if a sequence exists whose energy is lower than the proposed limit while still satisfying the discrete rigidity constraints, the identification of the limit objects as piecewise-constant orientation fields would fail.

minor comments (1)

[§2] Notation for the discrete energy and the continuum orientation field should be introduced with a single consistent symbol set early in the paper to avoid later redefinitions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the positive evaluation of its significance. We address each major comment below and will revise the manuscript accordingly to improve clarity.

read point-by-point responses

Referee: [Abstract; statement of main theorem (likely Theorem 1.1 or §3)] The central claim that the limiting energy density decomposes into twice the solid-vacuum transition energy rests on the assertion (stated in the abstract and presumably formalized in the hypotheses of the main Γ-convergence theorem) that solid-solid phase transitions with interpolating boundary layers are strictly higher in energy than two separate solid-vacuum interfaces. The precise conditions on the interaction potential that enforce this exclusion are not exhibited with a counter-example or a minimal set of inequalities that would allow verification for a concrete potential; without this, it is unclear whether the liminf inequality produces the claimed continuum energy for all admissible sequences.

Authors: The conditions ensuring that solid-solid interpolating layers are not energetically favorable are formalized in the hypotheses of Theorem 1.1 (specifically, the rigidity and frame-invariance assumptions on the interaction potentials detailed in Section 2). These prevent low-energy transitions between distinct orientations without passing through vacuum. We agree that a concrete example and a minimal set of inequalities would improve verifiability. In the revision we will add a remark after the statement of the main assumptions that provides (i) an explicit example potential satisfying the hypotheses and (ii) a short derivation showing why any interpolating layer costs at least as much as two solid-vacuum interfaces under those conditions. revision: yes
Referee: [Proof of Γ-convergence (likely §4 or §5)] In the proof of the liminf inequality, the argument that any interpolating layer can be replaced by two solid-vacuum interfaces without increasing energy must be checked against the finite-surface-energy scaling; if a sequence exists whose energy is lower than the proposed limit while still satisfying the discrete rigidity constraints, the identification of the limit objects as piecewise-constant orientation fields would fail.

Authors: The liminf proof in Section 4 proceeds by first establishing discrete rigidity at the finite-surface-energy scaling, which forces the configuration to be close to piecewise-constant orientations separated by interfaces. Any attempted interpolating layer between two solid grains is then shown to violate the rigidity estimates unless its energy is at least that of two independent solid-vacuum transitions; this replacement is performed by cutting the domain along the layer and applying the lower bound separately to each piece. We believe the argument already accounts for the scaling, but to make the replacement step fully explicit we will insert a short auxiliary lemma (or expanded paragraph) that directly compares the energy of an interpolating sequence to the sum of two solid-vacuum energies under the finite-surface-energy assumption. revision: partial

Circularity Check

0 steps flagged

Gamma-convergence proof is self-contained with no circular reduction

full rationale

The paper establishes a discrete-to-continuum limit via Gamma-convergence for atomistic energies favoring local isometry to a reference lattice. The limiting continuum energy is concentrated on grain boundaries between piecewise-constant orientation fields, with the explicit decomposition into twice the solid-vacuum interface energy following directly from the stated assumptions on rigid interactions that penalize interpolating boundary layers. This is a standard variational existence and characterization result; the assumptions are explicit inputs to the theorem rather than outputs derived from the conclusion. No load-bearing step reduces by construction to a fitted parameter, self-definition, or self-citation chain. The derivation chain is independent of the target result and relies on general techniques in the calculus of variations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the specific form of the atomistic energy (sum of local terms favoring isometry to a reference lattice) and on the rigidity assumptions that exclude interpolating layers; no free parameters or invented entities are introduced in the abstract.

axioms (2)

domain assumption The atomistic energy is a sum of particle energies that favor atoms being locally isometric to a reference lattice.
Stated in the abstract as the modeling choice that enables the discrete-to-continuum analysis.
domain assumption The discrete frame-invariant energy allows particle configurations without an underlying lattice assumed a priori.
Explicit modeling hypothesis that distinguishes the setting from classical lattice-based approaches.

pith-pipeline@v0.9.0 · 5463 in / 1492 out tokens · 31982 ms · 2026-05-10T02:32:50.745696+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

[1]

Allinger

N.L. Allinger . Molecular structure: understanding steric and electronic effects from molecular mechanics . John Wiley & Sons (2010)

work page 2010
[2]

Ambrosio, N

L. Ambrosio, N. Fusco, D. Pallara . Functions of bounded variation and free discontinuity problems . Oxford University Press (2000)

work page 2000
[3]

B\'etermin, L

L. B\'etermin, L. De Luca, M. Petrache . Crystallization to the Square Lattice for a Two-Body Potential . Arch.\ Ration.\ Mech.\ Anal.\ 181 (2021), 987--1053

work page 2021
[4]

Blanc, M

X. Blanc, M. Lewin . The crystallization conjecture: a review . EMS Surv.\ Math.\ Sci.\ 2 (2015), 225--306

work page 2015
[5]

A. Braides . -convergence for Beginners . Oxford University Press, Oxford 2002

work page 2002
[6]

Braides, S

A. Braides, S. Conti, A. Garroni . Density of polyhedral partitions . Calc.\ Var.\ Partial Differential Equations 56 (2017), Paper No. 28

work page 2017
[7]

Chambolle, L

A. Chambolle, L. Kreutz . Crystallinity of the Homogenized Energy Density of Periodic Lattice Systems . Multiscale Model. Simul. 21 (2023), 34--79

work page 2023
[8]

Dal Maso

G. Dal Maso . An introduction to -convergence . Birkh \"a user, Boston Basel Berlin 1993

work page 1993
[9]

Dolbilin, J

N. Dolbilin, J. Lagarias, M. Senechal . Multiregular Point Systems . Discrete\ Comput.\ Geom.\ 20 (1998), 477--498

work page 1998
[10]

De Luca, G

L. De Luca, G. Del Nin . A Crystallization Result in Two Dimensions for a Soft Disc Affine Potential . Anisotropic Isoperimetric Problems and Related Topics. INdAM 2022. Springer INdAM Series 62 (2024), 201--212

work page 2022
[11]

De Luca, G

L. De Luca, G. Friesecke . Crystallization in two dimensions and a discrete Gauss–Bonnet Theorem . J.\ Nonlinear Sci.\ 28 (2017), 69--90

work page 2017
[12]

W. E, D. Li . On the crystallization of 2D hexagonal lattices . Comm.\ Math.\ Phys.\ 286 (2009), 1099--1140

work page 2009
[13]

C Evans, R

L. C Evans, R. F. Gariepy . Measure theory and fine properties of functions . CRC Press, Boca Raton London New York Washington, D.C. 1992

work page 1992
[14]

Farmer, S

B. Farmer, S. Esedoḡlu, P. Smereka . Crystallization for a Brenner-like Potential . Commun.\ Math.\ Phys.\ 349 (2017), 1029--1061

work page 2017
[15]

Friedrich, L

M. Friedrich, L. Kreutz . Crystallization in the hexagonal lattice for ionic dimers . Math.\ Models Methods Appl.\ Sci.\ 29 (2019), 1853--1900

work page 2019
[16]

Friedrich, L

M. Friedrich, L. Kreutz . Finite crystallization and Wulff shape emergence for ionic compounds in the square lattice . Nonlinearity 33 (2020), 1240--1296

work page 2020
[17]

Friedrich, L

M. Friedrich, L. Kreutz . A Proof of Finite Crystallization via Stratification . J Stat Phys 190 (2023), 199

work page 2023
[18]

Friedrich, L

M. Friedrich, L. Kreutz, B. Schmidt . Emergence of rigid polycrystals from atomistic systems with Heitmann-Radin sticky disk energy . Arch.\ Ration.\ Mech.\ Anal.\ 240 (2021), 627--698

work page 2021
[19]

Friedrich, L

M. Friedrich, L. Kreutz, U. Stefanelli . Crystallization in the Winterbottom shape and sharp fluctuation laws . Preprint at arxiv:2509.05642 https://arxiv.org/pdf/2509.05642

work page arXiv
[20]

C. S. Gardner, C. Radin . The infinite-volume ground state of the Lennard-Jones potential . J.\ Stat.\ Phys.\ 20 (1979), 719--724

work page 1979
[21]

Heitmann, C

R. Heitmann, C. Radin . The ground state for sticky disks . J.\ Stat.\ Phys.\ 22 (1980), 281--287

work page 1980
[22]

J. Lee . Introduction to Smooth Manifolds . Springer (2012)

work page 2012
[23]

Mainini, P

E. Mainini, P. Piovano, U. Stefanelli . Finite crystallization in the square lattice . Nonlinearity 27 (2014), 717--737

work page 2014
[24]

Mainini, U

E. Mainini, U. Stefanelli . Crystallization in carbon nanostructures . Comm.\ Math.\ Phys.\ 328 (2014), 545--571

work page 2014
[25]

C. Radin . The ground state for soft disks . J.\ Stat.\ Phys.\ 26 (1981), 365--373

work page 1981
[26]

T. Schmidt . Strict interior approximation of sets of finite perimeter and functions of bounded variation . Proc. Am. Math. Soc. 143 (2015), 2069--2084

work page 2015
[27]

F. Theil . A proof of crystallization in two dimensions . Comm.\ Math.\ Phys.\ 262 (2006), 209--236

work page 2006
[28]

A. Vince . Periodicity, quasiperiodicity and Bieberbach's theorem on crystallographic groups . Am.\ Math.\ Mon.\ 104 (1997), 27--35

work page 1997

[1] [1]

Allinger

N.L. Allinger . Molecular structure: understanding steric and electronic effects from molecular mechanics . John Wiley & Sons (2010)

work page 2010

[2] [2]

Ambrosio, N

L. Ambrosio, N. Fusco, D. Pallara . Functions of bounded variation and free discontinuity problems . Oxford University Press (2000)

work page 2000

[3] [3]

B\'etermin, L

L. B\'etermin, L. De Luca, M. Petrache . Crystallization to the Square Lattice for a Two-Body Potential . Arch.\ Ration.\ Mech.\ Anal.\ 181 (2021), 987--1053

work page 2021

[4] [4]

Blanc, M

X. Blanc, M. Lewin . The crystallization conjecture: a review . EMS Surv.\ Math.\ Sci.\ 2 (2015), 225--306

work page 2015

[5] [5]

A. Braides . -convergence for Beginners . Oxford University Press, Oxford 2002

work page 2002

[6] [6]

Braides, S

A. Braides, S. Conti, A. Garroni . Density of polyhedral partitions . Calc.\ Var.\ Partial Differential Equations 56 (2017), Paper No. 28

work page 2017

[7] [7]

Chambolle, L

A. Chambolle, L. Kreutz . Crystallinity of the Homogenized Energy Density of Periodic Lattice Systems . Multiscale Model. Simul. 21 (2023), 34--79

work page 2023

[8] [8]

Dal Maso

G. Dal Maso . An introduction to -convergence . Birkh \"a user, Boston Basel Berlin 1993

work page 1993

[9] [9]

Dolbilin, J

N. Dolbilin, J. Lagarias, M. Senechal . Multiregular Point Systems . Discrete\ Comput.\ Geom.\ 20 (1998), 477--498

work page 1998

[10] [10]

De Luca, G

L. De Luca, G. Del Nin . A Crystallization Result in Two Dimensions for a Soft Disc Affine Potential . Anisotropic Isoperimetric Problems and Related Topics. INdAM 2022. Springer INdAM Series 62 (2024), 201--212

work page 2022

[11] [11]

De Luca, G

L. De Luca, G. Friesecke . Crystallization in two dimensions and a discrete Gauss–Bonnet Theorem . J.\ Nonlinear Sci.\ 28 (2017), 69--90

work page 2017

[12] [12]

W. E, D. Li . On the crystallization of 2D hexagonal lattices . Comm.\ Math.\ Phys.\ 286 (2009), 1099--1140

work page 2009

[13] [13]

C Evans, R

L. C Evans, R. F. Gariepy . Measure theory and fine properties of functions . CRC Press, Boca Raton London New York Washington, D.C. 1992

work page 1992

[14] [14]

Farmer, S

B. Farmer, S. Esedoḡlu, P. Smereka . Crystallization for a Brenner-like Potential . Commun.\ Math.\ Phys.\ 349 (2017), 1029--1061

work page 2017

[15] [15]

Friedrich, L

M. Friedrich, L. Kreutz . Crystallization in the hexagonal lattice for ionic dimers . Math.\ Models Methods Appl.\ Sci.\ 29 (2019), 1853--1900

work page 2019

[16] [16]

Friedrich, L

M. Friedrich, L. Kreutz . Finite crystallization and Wulff shape emergence for ionic compounds in the square lattice . Nonlinearity 33 (2020), 1240--1296

work page 2020

[17] [17]

Friedrich, L

M. Friedrich, L. Kreutz . A Proof of Finite Crystallization via Stratification . J Stat Phys 190 (2023), 199

work page 2023

[18] [18]

Friedrich, L

M. Friedrich, L. Kreutz, B. Schmidt . Emergence of rigid polycrystals from atomistic systems with Heitmann-Radin sticky disk energy . Arch.\ Ration.\ Mech.\ Anal.\ 240 (2021), 627--698

work page 2021

[19] [19]

Friedrich, L

M. Friedrich, L. Kreutz, U. Stefanelli . Crystallization in the Winterbottom shape and sharp fluctuation laws . Preprint at arxiv:2509.05642 https://arxiv.org/pdf/2509.05642

work page arXiv

[20] [20]

C. S. Gardner, C. Radin . The infinite-volume ground state of the Lennard-Jones potential . J.\ Stat.\ Phys.\ 20 (1979), 719--724

work page 1979

[21] [21]

Heitmann, C

R. Heitmann, C. Radin . The ground state for sticky disks . J.\ Stat.\ Phys.\ 22 (1980), 281--287

work page 1980

[22] [22]

J. Lee . Introduction to Smooth Manifolds . Springer (2012)

work page 2012

[23] [23]

Mainini, P

E. Mainini, P. Piovano, U. Stefanelli . Finite crystallization in the square lattice . Nonlinearity 27 (2014), 717--737

work page 2014

[24] [24]

Mainini, U

E. Mainini, U. Stefanelli . Crystallization in carbon nanostructures . Comm.\ Math.\ Phys.\ 328 (2014), 545--571

work page 2014

[25] [25]

C. Radin . The ground state for soft disks . J.\ Stat.\ Phys.\ 26 (1981), 365--373

work page 1981

[26] [26]

T. Schmidt . Strict interior approximation of sets of finite perimeter and functions of bounded variation . Proc. Am. Math. Soc. 143 (2015), 2069--2084

work page 2015

[27] [27]

F. Theil . A proof of crystallization in two dimensions . Comm.\ Math.\ Phys.\ 262 (2006), 209--236

work page 2006

[28] [28]

A. Vince . Periodicity, quasiperiodicity and Bieberbach's theorem on crystallographic groups . Am.\ Math.\ Mon.\ 104 (1997), 27--35

work page 1997