Spin Model for Flat-Foldability Diagnostics

• The paper introduces a framework where the flat-foldability of origami is determined by finding a frustration-free, zero-energy ground state in a spin-glass Hamiltonian.
• It systematically translates layer-ordering into discrete Ising spins and integrates two-spin, three-spin, and four-spin interactions to enforce crease constraints.
• The approach bridges combinatorial origami analysis with statistical mechanics, offering actionable diagnostics and insights for tuning origami metamaterials.

A spin model for flat-foldability diagnostics is an analytical framework that maps the global flat-foldability of origami crease patterns onto the zero-energy ground-state problem of a spin-glass Hamiltonian on a (hyper)graph. Such models encode the layer-ordering and blocking relationships between overlapping origami facets as Ising spins or related discrete variables. Forbidden layer-orderings, which would violate the physical requirement that facets do not intrude upon one another in the folded state, are assigned positive energy contributions. The existence of a frustration-free ground state (minimum energy zero) of the spin Hamiltonian is necessary and sufficient for global flat-foldability. This approach connects combinatorial origami foldability analysis to paradigms in statistical mechanics and combinatorial optimization (Nakajima, 2024, Assis, 2017).

1. Formalization of Spin Variables and Layer-Ordering

In the fully generalized spin model construction, one begins with a locally flat-foldable origami crease pattern partitioned into $N$ polygonal facets. From the pre-folded diagram (i.e., the geometric planar projection constructed by flattening all local folds and ignoring mountain/valley assignments), overlapping pairs of facets $(j,k)$ ($j < k$) are identified. To each pair, an Ising variable $\sigma_{j,k} \in \{+1, -1\}$ is assigned to encode the relative stacking:

• $\sigma_{j,k} = +1$ means $k$ is above $j$; $\sigma_{j,k} = -1$ means $k$ is below $j$.
• Consistency under relabeling is imposed: $(j,k)$0.

To systematize indices, signs, and account for variable redundancy, the reduced set

$(j,k)$1

is defined such that $(j,k)$2. The full configuration space thus has cardinality $(j,k)$3 for $(j,k)$4 facets (Nakajima, 2024).

2. Construction of the Layer-Ordering Hamiltonian

The key structure is a Hamiltonian $(j,k)$5 defined as a sum of positive semi-definite energy penalties $(j,k)$6 for forbidden local layer orderings. This Hamiltonian contains three classes of terms:

• Two-spin (pairwise) terms $(j,k)$7: Penalize configurations where a facet $(j,k)$8 intrudes between $(j,k)$9 and $j < k$0 in series along adjacent creases.

$j < k$1

• Four-spin (quartic) terms $j < k$2: Penalize forbidden intrusions when creases coincide geometrically.

$j < k$3

• Three-spin cycle terms $j < k$4: Exclude cyclic ordering among fully overlapping facet triples.

$j < k$5

The full Hamiltonian is then a sum over all necessary interactions as dictated by the geometry of the pre-folded diagram:

$j < k$6

(Nakajima, 2024).

3. Diagnostic Algorithm for Flat-Foldability

The spin model enables a concrete algorithmic procedure for diagnosing global flat-foldability:

1. Input: An arbitrary straight-crease pattern with $j < k$7 internal vertices.
2. Local flat-foldability check: Apply Kawasaki’s theorem at each vertex (even degree, alternating angle sum zero). If any vertex fails, local flat-foldability is rejected.
3. Pre-folded diagram construction: Simulate all local folds; identify overlapping facet pairs $j < k$8.
4. Spin Introduction: Assign independent Ising spins $j < k$9 for each facet pair.
5. Hamiltonian Assembly:
   • Add $\sigma_{j,k} \in \{+1, -1\}$0 for every series pair of creases.
   • Add $\sigma_{j,k} \in \{+1, -1\}$1 for each set of coincident creases.
   • Add $\sigma_{j,k} \in \{+1, -1\}$2 for overlapping triples.
6. Interaction Graph Reduction: Prune all dangling spins (leaves) and split into connected components (in $\sigma_{j,k} \in \{+1, -1\}$3 time).
7. Ground-State Search: For each hypergraph component, search for a ground-state spin configuration using, e.g., exhaustive search ($\sigma_{j,k} \in \{+1, -1\}$4 assignments), branch-and-bound, or simulated annealing.
8. Diagnosis: If $\sigma_{j,k} \in \{+1, -1\}$5 for all components, report "flat-foldable;" otherwise, report "not flat-foldable" (Nakajima, 2024).

This framework is provably complete, with correctness reduced to finding (or not) a zero-energy ground state of the constructed Hamiltonian.

4. Frustration, Cycles, and Flat-Foldability Obstructions

A central feature of the diagnostic method is its correspondence between frustrated cycles in the (hyper)graph of spin interactions and obstructions to flat-foldability. All $\sigma_{j,k} \in \{+1, -1\}$6, $\sigma_{j,k} \in \{+1, -1\}$7, and $\sigma_{j,k} \in \{+1, -1\}$8 couplings are $\sigma_{j,k} \in \{+1, -1\}$9, making the following criteria exact:

• For any even-length cycle of 2-spin edges $\sigma_{j,k} = +1$0, non-satisfiability (frustration) is detected if $\sigma_{j,k} = +1$1.
• For composite cycles including 4-spin hyperedges, the product of all $\sigma_{j,k} = +1$2’s and $\sigma_{j,k} = +1$3’s along the chain is checked; a value $\sigma_{j,k} = +1$4 signals frustration.
• Presence of a single frustrated loop (or hyper-loop) implies $\sigma_{j,k} = +1$5 and hence, global flat-foldability fails.

Explicit illustrative examples for six-spin cycles and coupled cycles with four-spin terms are provided in (Nakajima, 2024), solidifying the correspondence between spin-glass frustration and origami folding infeasibility.

5. Connections to Exactly Solvable Vertex Models and Tiling Cases

For certain classes of quadrilateral origami tilings (e.g., Miura-ori, Barreto’s Mars), spin models admit exact solvable constructions mapped to vertex models from statistical mechanics:

• Crease-Spin Models: Edges receive $\sigma_{j,k} = +1$6 for mountain/valley; faces may receive three-state colors modeling layer ordering.
• Vertex Model Mapping: Locally flat-foldable configurations correspond to the eight allowed Maekawa/Kawasaki-compliant spin patterns (odd 8-vertex model); the partition function is exactly solvable via free-fermion and Pfaffian/dimer methods.
• Layer Order and Coloring: Layer order can be mapped to Baxter's three-coloring model, where adjacent face colors encode stacking; critical points can be analyzed in terms of fugacities for 'defect' configurations.
• Relaxing Flat-Foldability: Allowing all-crease (even 8-vertex) configurations leads to the 16-vertex model; phase transitions and tunability differ from strictly flat-foldable cases.

This spin-model treatment for quadrilateral tilings enables exact calculation of defect densities, analytic free energies, and explicit relations to mechanical properties like the elastic modulus (Assis, 2017).

6. Computational Complexity and Practical Considerations

• The decision problem for global flat-foldability, when formulated as finding a zero-energy ground state of the constructed Hamiltonian, is NP-hard in the worst case (by Bern–Hayes reduction).
• Preprocessing steps (local theorem checks, hypergraph sparsification, componentization) are polynomial, $\sigma_{j,k} = +1$7.
• Failure modes are entirely captured by the existence of frustrated cycles/hypercycles.
• For practically relevant origami tilings (e.g., Miura-ori), exactly solvable model frameworks facilitate deeper analysis of tunability and phase transitions, including lattice-gas interpretations and explicit connections to mechanical moduli (Nakajima, 2024, Assis, 2017).

7. Comparative Analysis and Applicability to Origami Metamaterials

Spin model diagnostics provide both a universal test for global flat-foldability and a fertile computational/statistical mechanics framework for exploring origami design spaces:

• Miura-ori and trapezoid tilings exhibit sharp phase transitions in defect density at critical fugacity, corresponding to sudden loss of crease order; Barreto's Mars allows continuous tunability without phase transition, implying greater mechanical stability but less "switchability."
• The defect density for Miura-ori as a function of lattice-gas fugacity $\sigma_{j,k} = +1$8 is given explicitly by

$\sigma_{j,k} = +1$9

where $k$0 is the complete elliptic integral of the first kind (Assis, 2017).

• Experimentally, the effective in-plane elastic modulus $k$1 is linear in $k$2, linking spin model diagnostics and material tunability.

This suggests that spin model frameworks not only diagnose flat-foldability but offer a quantitative, theoretically robust pathway to analyze and engineer the tunability of origami metamaterials, subject to the intrinsic combinatorial and physical constraints of flat-folding (Nakajima, 2024, Assis, 2017).