M4N-type Ordered Phases Overview

• M4N-type ordered phases are states characterized by multiple preferred length scales and competing lattice symmetries, driven by multi-mode competition.
• They emerge in systems such as multi-mode phase field crystal models, generalized XY models, and magnetic structures, displaying both conventional and topological order.
• Computational analyses reveal complex interfacial behaviors and phase transitions, providing practical insights into nucleation, coarsening, and defect dynamics.

Searching arXiv for papers relevant to “M4N-type Ordered Phases” and the cited IDs. M4N-type ordered phases, in the available literature, do not appear as a single uniformly standardized category. In explicit usage, the term refers to ordered states with multiple characteristic length scales and competing lattice symmetries in a generic multi-mode phase field crystal model (Mkhonta et al., 2013). Closely related papers do not always use the notation explicitly, but they discuss new ordered phases generated by three competing length scales (Mkhonta et al., 2013), competing ferromagnetic and higher-harmonic couplings (Poderoso et al., 2010), antiferromagnetic and third-order antinematic couplings on a triangular lattice (Lach et al., 2020), symmetry relations between sublattices in bipartite quasicrystals (Ji et al., 27 Nov 2025), nodal magnetic order with d-, g-, or i-wave character (Jungwirth et al., 2024), and three-dimensional Abelian topological order with deconfined point- and loop-like excitations (Sagi et al., 2018). This suggests a broad family of ordered phases in which competition among modes, harmonics, or symmetry operations selects structures not reducible to a single conventional ordered state.

1. Terminological scope and defining features

In the multi-mode phase field crystal literature, M4N-type ordered phases are described as ordered states with multiple characteristic length scales and competing lattice symmetries (Mkhonta et al., 2013). The same paper emphasizes that three competing length scales in a multi-mode phase field crystal model are enough to generate all five Bravais lattices, together with honeycomb, kagome, dimers, and hybrid phases. In this usage, the essential mechanism is reciprocal-space competition among several preferred wave numbers $Q_i$, with phase selection controlled by which wave-vector triads can most efficiently lower the free energy.

Several adjacent literatures broaden this picture without adopting identical terminology. In generalized $XY$ models, higher harmonics and frustration generate additional ordered phases beyond the standard ferromagnetic and nematic ones, including pseudo-ferromagnetic phases $F_1$ and $F_2$ on the square lattice (Poderoso et al., 2010) and a canted antiferromagnetic phase on the triangular lattice (Lach et al., 2020). In magnetic systems, the relevant classification variable is not only net magnetization but also the symmetry of the spin density and the momentum-space nodal structure, as in altermagnets and other nodal magnetically ordered phases (Jungwirth et al., 2024). In 3D interacting-electron constructions, the phrase is only approximate: the relevant phases are fully gapped topologically ordered descendants of a composite Weyl semimetal, with deconfined point-like and line-like excitations and non-trivial mutual statistics (Sagi et al., 2018).

A recurrent misconception is that all such phases are conventional symmetry-breaking states. The cited literature does not support that simplification. Some examples are quasi-long-range ordered phases with separate chiral long-range order (Lach et al., 2020), while others are fully gapped Abelian topological orders rather than ordinary symmetry-breaking phases (Sagi et al., 2018).

2. Multi-mode phase field crystal realizations

The most explicit formulation is the generic multi-mode phase field crystal model, written in terms of a rescaled particle-density field $\psi(\mathbf r,t)$ evolving on diffusive time scales (Mkhonta et al., 2013): $$\mathcal{F} = \int d\vec{r}\, \left\{ \frac{\psi}{2}\left(r+\lambda\prod_{i=0}^{N-1}\left[(Q_i^2+\nabla^2)^2+b_i\right]\right)\psi -\frac{\tau}{3}\psi^3+\frac{\psi^4}{4}\right\},$$ with dynamics

$$\frac{\partial \psi}{\partial t} = \nabla^2\left\{ \left(r+\lambda\prod_{i=0}^{N-1}\left[(Q_i^2+\nabla^2)^2+b_i\right]\right)\psi -\tau\psi^2+\psi^3 \right\}.$$

The core physical content is that the free energy favors density modulations at multiple preferred wave numbers. In Fourier space,

$$\psi(\mathbf r)=\psi_0+\sum_{\mathbf q} A_{\mathbf q}e^{i\mathbf q\cdot\mathbf r},$$

and near freezing the cubic term is especially important. Following the Alexander–McTague logic, the favored crystal is determined by density waves whose wave vectors form closed triads,

$\mathbf q_1+\mathbf q_2+\mathbf q_3=0.$

The ordered phase is therefore selected by the joint action of mode excitation, triad formation, and lattice-symmetry matching.

For $N=3$, the model realizes all five 2D Bravais lattices: triangular, square, rectangular, rhombic, and oblique (Mkhonta et al., 2013). With the specific choice

$XY$0

the simulations produce Tri0, Tri1, Tri2 triangular phases, Honeycomb, Kagome, Rectangular, Dimer, and Intermediate hybrid states. The same paper reports other choices of $XY$1 that generate square, rhombic, oblique, and square-dimer structures. The stability logic is explicit: if one mode dominates strongly, simpler one-mode structures like triangular order appear; if two modes are balanced, two-mode structures like honeycomb, kagome, or rhombic order appear; if all three modes participate comparably, more complex or lower-symmetry states appear, including oblique and hybrid phases.

This construction is also dynamical. The paper examines non-equilibrium transitions such as

$XY$2

and reports nucleation, growth and coarsening, formation and annihilation of topological defects, and competition between different ordered domains. The same framework is used to compute elastic properties, especially for the dimer phase, whose smaller shear modulus relative to the triangular phase is attributed to dimer rotation.

3. Generalized $XY$3 models and frustration-induced hybrid phases

A second major realization comes from generalized $XY$4 models with competing harmonics. On the square lattice, the Hamiltonian

$XY$5

interpolates between ordinary ferromagnetic $XY$6 order and $XY$7-fold pseudonematic order (Poderoso et al., 2010). For $XY$8, the low-temperature structure remains the familiar paramagnet–nematic–ferromagnet sequence, with KT, Ising, or 3-state Potts criticality depending on $XY$9. For $F_1$0, the phase diagram changes qualitatively and two new pseudo-ferromagnetic phases emerge, $F_1$1 and $F_1$2. $F_1$3 has one preferred spin orientation, while $F_1$4 has four preferred spin orientations with different weights in the $F_1$5 example. Both have quasi-long-range order and broken reflection symmetry.

On the triangular lattice, geometric frustration produces a distinct extension of this pattern. The generalized model with AFM and third-order antinematic couplings is

$F_1$6

with

$F_1$7

Monte Carlo simulations find three quasi-long-range ordered phases: AFM, CAFM, and AN3 (Lach et al., 2020). The CAFM phase is the central new result. It appears at low temperature, is wedged between the AFM and AN3 phases, and contains simultaneous AFM and AN3 QLRO. It is absent in the analogous AN2 model.

The order parameters are

$F_1$8

with $F_1$9 detecting magnetic order and $F_2$0 detecting third-order nematic order. The standard and generalized staggered chiralities are denoted $F_2$1 and $F_2$2. A particularly important result is that all three QLRO phases also have true long-range order of both standard and generalized chiralities, and that both chiralities vanish simultaneously at their own second-order transitions slightly above the magnetic and nematic order-disorder transition temperatures. The high-temperature transitions to the paramagnet are BKT-type, the AFM–CAFM phase transition belongs to the weak Ising universality class, and the AN3–CAFM phase transition belongs to the weak three-state Potts universality class. The chiral and magnetic or nematic sectors are therefore decoupled.

4. Magnetic classifications: Néel, altermagnetic, and nodal ordered phases

Magnetic ordered phases provide a third setting in which the notion of M4N-type order is used or inferred. In half-filled Hubbard models on 2D bipartite quasicrystals, sign-problem-free projector quantum Monte Carlo yields Néel ordered states in all studied cases, but the magnetic class is not fixed a priori (Ji et al., 27 Nov 2025). The paper proposes a general criterion: if the two sublattices are related by inversion, the Néel state is AFM; if they are related by another point-group operation, it is AM; if no point-group symmetry relates them, the Néel state is generally FM. Explicit examples are AM in the two $F_2$3-symmetric Thue-Morse quasicrystals, AFM in the $F_2$4-symmetric Thue-Morse quasicrystal, and FM in the $F_2$5-symmetric Penrose and $F_2$6-symmetric Ammann-Beenker quasicrystals.

A broader symmetry-based classification appears in the review of altermagnets and nodal magnetically ordered phases (Jungwirth et al., 2024). Altermagnets are defined as collinear magnets with zero net magnetization, spin-split electronic bands away from symmetry-protected nodes, and alternating-sign even-parity spin splitting in momentum space. They feature d-, g-, or i-wave magnetic ordering, with spin-degenerate nodes and a characteristic alternation of spin polarization. The spin symmetry classes are written as

$F_2$7

Within this framework, 4 spin Laue groups correspond to d-wave, 4 to g-wave, and 2 to i-wave altermagnets, with 2, 4, and 6 nodal surfaces, respectively, crossing the $F_2$8-point.

The review also extends the classification to non-collinear spin densities and predicted p-wave magnets, characterized by

$F_2$9

An objective caveat is that not every candidate realization is settled experimentally: the paper identifies RuO$\psi(\mathbf r,t)$0 as the prototypical metallic altermagnet, though the magnetic ground state remains debated experimentally.

5. Three-dimensional topological order as an alternative usage

In a different meaning, the literature associates M4N-type order with 3D Abelian topological order of the $\psi(\mathbf r,t)$1-$\psi(\mathbf r,t)$2 type (Sagi et al., 2018). The starting point is a stack of alternating thin topological-insulator and normal-insulator layers, with a single Dirac cone of surface electrons on each interface. After applying 2+1D fermionic duality layer by layer, strong inter-layer correlations generate a genuine 3+1D gauge theory of neutral fermions coupled to an emergent $\psi(\mathbf r,t)$3 gauge field: $\psi(\mathbf r,t)$4

The parent state can become a gapless composite Weyl semimetal under a time-reversal-breaking mass term, or a fully gapped phase under dual-fermion pairing,

$\psi(\mathbf r,t)$5

The fully gapped descendants host deconfined point-like and line-like excitations with non-trivial mutual statistics. The point-like excitations are gapped dual fermions, while the line-like excitations are vortex lines of the paired dual-fermion condensate. When a dual fermion loops around a vortex line, it picks up a phase of $\psi(\mathbf r,t)$6. The paper explicitly identifies this as the three-dimensional incarnation of the toric code.

This use of the term is conceptually distinct from the symmetry-breaking examples above. The phases are electronic insulators, but their defining structure is topological rather than conventional magnetic or crystalline order. The paper further discusses symmetry enrichment, including an antiunitary symmetry $\psi(\mathbf r,t)$7 with $\psi(\mathbf r,t)$8, symmetry-protected Majorana surface cones, and a ground-state degeneracy of eight on the 3-torus.

6. Interfaces, compatibility, and computation of ordered-phase structure

The interfacial problem for ordered phases is addressed by a general variational framework based on free-energy minimization with boundary conditions compatible with both bulk phases (Xu et al., 2015). The paper does not explicitly use the term M4N-type, but it treats ordered periodic, quasiperiodic, and modulated phases, including lamellar, cylindrical, spherical, and gyroid morphologies in the Landau–Brazovskii model. The free-energy density is

$\psi(\mathbf r,t)$9

with conserved order parameter $$\mathcal{F} = \int d\vec{r}\, \left\{ \frac{\psi}{2}\left(r+\lambda\prod_{i=0}^{N-1}\left[(Q_i^2+\nabla^2)^2+b_i\right]\right)\psi -\frac{\tau}{3}\psi^3+\frac{\psi^4}{4}\right\},$$0.

The key methodological point is compatibility. In the transverse plane, the admissible function space is constructed so that each bulk phase individually satisfies the boundary conditions. In the normal direction, because the energy involves second derivatives, $$\mathcal{F} = \int d\vec{r}\, \left\{ \frac{\psi}{2}\left(r+\lambda\prod_{i=0}^{N-1}\left[(Q_i^2+\nabla^2)^2+b_i\right]\right)\psi -\frac{\tau}{3}\psi^3+\frac{\psi^4}{4}\right\},$$1 and $$\mathcal{F} = \int d\vec{r}\, \left\{ \frac{\psi}{2}\left(r+\lambda\prod_{i=0}^{N-1}\left[(Q_i^2+\nabla^2)^2+b_i\right]\right)\psi -\frac{\tau}{3}\psi^3+\frac{\psi^4}{4}\right\},$$2 are fixed at the domain boundaries. The interface is then obtained as a local minimizer of the free energy. The paper argues that no artificial pinning constraint is needed for ordered phases, because anisotropy creates local minima that stabilize the interface position.

The resulting energy landscape is wavy rather than flat. For both cylindrical–gyroid and lamellar–gyroid interfaces, the paper reports four local minima within one period. Orientation and translation can produce planar interfaces, non-planar interfaces, or an intermediate ordered state such as a perforated-layer-like structure. This makes interfacial structure part of the ordered-phase problem itself rather than a purely geometric boundary layer.

The combined literature therefore presents M4N-type ordered phases not as one narrowly defined class, but as a family of ordered or topologically ordered states selected by competition among modes, harmonics, chirality sectors, symmetry operations, or emergent gauge fields. A plausible implication is that the term is most useful when it emphasizes the mechanism of phase selection—multi-mode competition, frustration, symmetry-enforced nodal structure, or gauge-theoretic fractionalization—rather than a single universal order parameter.