Zeta Phases in Materials Systems

• Zeta phases are a system-specific label applied to chemically and crystallographically distinct states across oxygen, nitrogen, carbides, alloys, and other materials.
• They are defined by unique bonding, symmetry, and stability criteria identified via first-principles calculations, X-ray diffraction, and spectroscopy.
• Their study offers insights into phase transformation, metastability, and nucleation, aiding the development of advanced ceramics and alloy systems.

Searching arXiv for recent and relevant papers on zeta phases across materials systems. In contemporary condensed-matter, alloy, and ceramic literature, the designation $\zeta$ is applied to several distinct phases whose symmetry, bonding, dimensionality, and stability are strongly system dependent. Reported examples include high-pressure solid oxygen, a bi-layered precipitate in Al–Ag, vacancy-ordered layered carbides and nitrides, a metallic few-layer tellurium allotrope, molecular nitrogen polytypes near megabar pressures, and a non-icosahedral boron allotrope; in colloid science, by contrast, phase behavior has also been analyzed for spheres with a constant zeta-potential, which is an electrostatic boundary condition rather than a crystal phase [(Elatresh et al., 2023); (Zhang et al., 2017); (Nykiel et al., 25 Aug 2025); (Wang et al., 2018); (Goncharov et al., 17 Apr 2026); (Chuvashova et al., 2017); (Smallenburg et al., 2010)]. This suggests that “zeta phase” is a conventional, system-specific phase label rather than a universal structural class.

1. Taxonomic scope and nomenclature

The literature assigns the $\zeta$ label to phases that are chemically and crystallographically unrelated. In the available reports, the label spans molecular solids under extreme compression, coherent precipitates in alloys, vacancy-ordered layered ceramics, and low-dimensional elemental allotropes. The commonality is terminological rather than structural [(Elatresh et al., 2023); (Zhang et al., 2017); (Nykiel et al., 25 Aug 2025); (Wang et al., 2018); (Goncharov et al., 17 Apr 2026); (Chuvashova et al., 2017); (Smallenburg et al., 2010)].

System	$\zeta$ designation	Defining feature
Solid oxygen	candidate and predicted room-temperature $\zeta$ phase	molecular O$_2$ structures including Pm, Pnma, P2$_1$/m, P6$_3$/mmc
Al–Ag alloy	$\zeta$ precipitate	alternating Ag-rich and Al-rich bilayers on $\{111\}_{\rm Al}$
Layered carbides/nitrides	$\zeta$-like phases	vacancy-ordered $\zeta$0 layered hexagonal/trigonal structures
Tellurium few-layers	$\zeta$1 allotrope	metallic P4/mmm layered phase with square-net sublayers
Molecular nitrogen	$\zeta$2-N$\zeta$3 and $\zeta$4-N$\zeta$5	monoclinic C2/c molecular phases near 1 Mbar
Boron	$\zeta$6-B	orthorhombic Cmce $\zeta$7-Ga-type non-icosahedral allotrope
Colloidal spheres	constant zeta-potential	electrostatic charging condition governing fluid/bcc/fcc phase diagrams

A recurrent misconception is that all $\zeta$8 phases form a homologous family. The reported data do not support that interpretation. In one case $\zeta$9 denotes a metallic molecular oxygen phase; in another it denotes a coherent AgAl precipitate; in another it denotes vacancy-ordered carbides and nitrides with 16 stacking motifs; and in boron it denotes an $\zeta$0-Ga-type three-dimensional covalent framework. The terminology therefore functions locally within each phase diagram.

2. High-pressure elemental $\zeta$1 phases

In solid oxygen, first-principles structure searches and molecular dynamics identify several crystalline candidates for the $\zeta$2 phase at 0 K, with space groups Pnma, P2$\zeta$3/m, Pm, and P6$\zeta$4/mmc. The calculations used USPEX coupled to DFT–PBE relaxations and HSE06 enthalpy corrections. Relative to the $\zeta$5 phase, taken as C2/m O$\zeta$6, the relevant thermodynamic quantity is

$\zeta$7

Within HSE06, $\zeta$8-O$\zeta$9 remains lower in enthalpy for $\zeta$0 GPa, whereas above $\zeta$1 GPa the Pnma and Pm structures become favored; at 300 K, HSE06-corrected Gibbs free energies place Pm lowest from 111 to 140 GPa. The predicted room-temperature $\zeta$2 phase is therefore the 24-atom Pm structure, which retains molecular O$\zeta$3 bonds of about 1.18–1.19 Å and is metallic (Elatresh et al., 2023).

That identification is supported, but not made trivial, by comparison to experiment. PowderCell-simulated X-ray diffraction at 116 GPa gives the best match for Pm to the broad, overlapping peaks near $\zeta$4. Raman/IR mode counting disfavors the ideal P6$\zeta$5/mmc model because it yields only three Raman-active modes, whereas experiment in $\zeta$6 oxygen shows at least seven modes in the 400–900 cm$\zeta$7 interval. Superconductivity estimates from the McMillan formula give $\zeta$8 K for Pnma, $\zeta$9 K for Pm, $_2$0 K for C2/m, $_2$1 K for P6$_2$2/mmc, and $_2$3 K for P2$_2$4/m, to be compared with the measured 0.6 K; none is exact, but Pm remains compatible in order of magnitude (Elatresh et al., 2023).

Molecular nitrogen exhibits a different use of the same label. At 98 GPa and 293 K, $_2$5-N$_2$6 is monoclinic C2/c with $_2$7 N$_2$8 per cell and lattice parameters $_2$9 Å, $_1$0 Å, $_1$1 Å, $_1$2. Laser heating of $_1$3-N$_1$4 to 1800–2500 K at 78–98 GPa yields a new polytype, $_1$5-N$_1$6, also C2/c but with a tripled $_1$7 axis, $_1$8 Å, and $_1$9. DFT/PBE enthalpies place $_3$0 within 2–3 meV/atom of $_3$1 above 35 GPa and slightly more stable for $_3$2 GPa; Raman spectroscopy reveals additional low-frequency and vibron features, and the phase likely corresponds to the previously reported $_3$3-N$_3$4 (Goncharov et al., 17 Apr 2026).

In boron, $_3$5 denotes a non-molecular, non-icosahedral allotrope synthesized at 115 GPa and 2100 K. $_3$6-B adopts the orthorhombic $_3$7-Ga-type structure, space group Cmce, with $_3$8 and lattice parameters $_3$9 Å, $\zeta$0 Å, $\zeta$1 Å at 115(2) GPa. Each B atom is seven-fold coordinated, with in-plane bond lengths 1.66(1), 1.72(1), and 1.75(1) Å and short interlayer bonds 1.59(1) Å, generating a fully three-dimensional covalent network. The transformation from $\zeta$2-B is accompanied by a $\zeta$3 volume collapse, and the second-order Birch–Murnaghan fit at 115 GPa gives $\zeta$4 GPa with $\zeta$5 Å$\zeta$6 (Chuvashova et al., 2017).

These three elemental cases already show that the $\zeta$7 label covers qualitatively different bonding regimes: molecular metallic oxygen, molecular nitrogen polytypes, and a covalent non-icosahedral boron framework.

3. Precipitate and intermediate $\zeta$8 phases in alloys

In the Al–Ag system, Zhang et al. reported a new precipitate phase $\zeta$9 by scanning transmission electron microscopy. It is a modulated, shear-free layered structure composed of alternating bilayers enriched in Ag or Al. The repeat consists of six Ag-rich and six Al-rich $\{111\}_{\rm Al}$0 planes, giving twelve $\{111\}_{\rm Al}$1 layers per repeat with the chemical modulation

$\{111\}_{\rm Al}$2

and periodicity

$\{111\}_{\rm Al}$3

In a hexagonal description the phase is R$\{111\}_{\rm Al}$4m with $\{111\}_{\rm Al}$5 Å experimentally and $\{111\}_{\rm Al}$6 Å, while the simplest trigonal unit cell is described as P3 and contains 12 sites with overall AgAl composition in the pure-layer model (Zhang et al., 2017).

Its formation pathway is explicitly intermediate. Samples oil-quenched and aged at 200 $\{111\}_{\rm Al}$7C for 7 days develop large $\{111\}_{\rm Al}$8–GP zones of about 20–25 nm. In situ heating at 150–200 $\{111\}_{\rm Al}$9C inside the TEM shows that after as little as 3 min at 200 $\zeta$0C, $\zeta$1 ordering appears within the GP zone. On longer times, typically more than 10–20 min or at higher temperature, newly nucleated $\zeta$2 (Al$\zeta$3Ag, HCP) plates form at the edge of $\zeta$4 regions and grow at the expense of $\zeta$5, completing the sequence $\zeta$6 (Zhang et al., 2017).

The thermodynamic interpretation is that $\zeta$7 is a local energy minimum. DFT formation energies per Ag atom give $\zeta$8 meV/Ag for isolated substitutional Ag, $\zeta$9 meV/Ag for an infinite Ag monolayer on $\zeta$00, and a minimum separation energy of about $\zeta$01 meV/Ag when two Ag layers are separated by two Al layers along $\zeta$02. Periodic repetition of that Ag$\zeta$03Al$\zeta$04 motif yields the $\zeta$05 structure, with $\zeta$06 meV/Ag. The broader energy sequence is solid solution $\zeta$07 small Ag clusters $\zeta$08 planar clusters $\zeta$09 $\zeta$10 meV) $\zeta$11 $\zeta$12 meV) $\zeta$13 $\zeta$14 meV), so $\zeta$15 is more stable than $\zeta$16 but less stable than $\zeta$17 (Zhang et al., 2017).

The metastability has a direct thermodynamic basis. The strain-energy estimate

$\zeta$18

gives about 3 meV/atom using $\zeta$19 GPa, $\zeta$20, and $\zeta$21 Å$\zeta$22. Configurational entropy favors the more disordered $\zeta$23 state, with a Bragg–Williams upper bound

$\zeta$24

per atom for Al–40 at% Ag, corresponding to $\zeta$25 meV/atom at 200 $\zeta$26C. Because this is comparable in magnitude to the enthalpic gain $\zeta$27 meV, $\zeta$28 remains metastable rather than terminally stable (Zhang et al., 2017).

4. Layered and low-dimensional $\zeta$29 phases

A different and more systematic use of the $\zeta$30 designation appears in hexagonal, layered carbides and nitrides studied as ultra-high temperature ceramics. In that setting, $\zeta$31- and $\zeta$32-type structures have general formula $\zeta$33 with $\zeta$34, where $\zeta$35 is a transition metal and $\zeta$36 is C or N. The defining structural feature of the vacancy-ordered $\zeta$37 phase is that every second M–M layer contains a plane of X vacancies. The reported prototype space groups are P6$\zeta$38mc, R$\zeta$39m, P$\zeta$40m1, and P3m1, and stacking is indexed by the Jagodzinski–Wyckoff descriptors $\zeta$41 and $\zeta$42; for example, $\zeta$43-Ta$\zeta$44C$\zeta$45 is described by hcch-hcch-hcch (Nykiel et al., 25 Aug 2025).

High-throughput DFT finds substantial chemical breadth for this motif. Across 11 transition metals and both C and N, 67 previously unreported hexagonal, layered materials lie within 0.1 eV/atom of the convex hull, and 9 lie exactly on the hull. Selected examples on or near the hull are Ta$\zeta$46C$\zeta$47, Nb$\zeta$48C$\zeta$49, and Zr$\zeta$50C$\zeta$51 with $\zeta$52 eV/atom, Nb$\zeta$53C$\zeta$54 at $\zeta$55 eV/atom, Mo$\zeta$56C at $\zeta$57 eV/atom, and V$\zeta$58C$\zeta$59 at $\zeta$60 eV/atom. The authors classify phases as stable if $\zeta$61, near-stable if $\zeta$62 eV/atom, and synthesizable if $\zeta$63 eV/atom (Nykiel et al., 25 Aug 2025).

Elastic and thermal indicators further motivate the ceramic interpretation. Representative values include $\zeta$64 GPa, $\zeta$65 GPa, $\zeta$66 GPa, $\zeta$67 GPa, $\zeta$68 GPa, $\zeta$69 GPa, and $\zeta$70 K for $\zeta$71-Ta$\zeta$72C$\zeta$73. The final melting temperature estimates average a Lindemann model with parsimonious neural networks, and five newly stable candidates are predicted above 2500 K: Nb$\zeta$74C$\zeta$75 at $\zeta$76 K, several Mo$\zeta$77C stackings at $\zeta$78–2750 K, and Zr$\zeta$79C$\zeta$80 at $\zeta$81 K. Additional reported values are $\zeta$82 K for V$\zeta$83C$\zeta$84, $\zeta$85 K for Nb$\zeta$86C$\zeta$87, and $\zeta$88 K for Ta$\zeta$89C$\zeta$90 (Nykiel et al., 25 Aug 2025).

Few-layer tellurium supplies yet another low-dimensional $\zeta$91 allotrope. Wang et al. identified a $\zeta$92 phase consisting, at the monolayer level, of three Te sublayers stacked along $\zeta$93, each sublayer forming a perfect square net. The symmetry is tetragonal P4/mmm, with $\zeta$94–3.15 Å and vacuum spacing chosen above 15 Å in slab calculations. The monolayer is metallic, and the phase remains more stable than the monolayer $\zeta$95 phase by about 29 meV/Te; the bilayer is about 35 meV/Te more stable than the corresponding $\zeta$96 bilayer. The energy difference with $\zeta$97 decreases with thickness and vanishes at four layers, i.e. 12 sublayers, so the crossover thickness is approximately 4 layers (Wang et al., 2018).

The electronic and interlayer characteristics of Te $\zeta$98 are likewise distinctive. The phase is metallic in PBE+SOC, with bands crossing $\zeta$99 along $\zeta$00–X, X–M, and M–$\zeta$01, and SOC does not open a gap. Strong interlayer coupling gives layer-dependent quantum-well states; the average interlayer stabilization relative to the monolayer is about $\zeta$02 meV/Te for the bilayer, $\zeta$03 meV/Te for the trilayer, and $\zeta$04 meV/Te for the four-layer system. Charge doping further stabilizes $\zeta$05, and for $\zeta$06 e/Te the $\zeta$07–$\zeta$08 crossover thickness increases beyond 4 layers (Wang et al., 2018).

5. Thermodynamic and structural criteria of identification

Across these systems, identification of a $\zeta$09 phase is primarily a problem of relative stability under constrained thermodynamic variables. In high-pressure oxygen, the central quantities are enthalpy and Gibbs free energy relative to $\zeta$10-O$\zeta$11,

$\zeta$12

$\zeta$13

with $\zeta$14 obtained from the MD-derived vibrational density of states. In layered carbides and nitrides, the corresponding stability measure is the formation enthalpy

$\zeta$15

followed by convex-hull analysis against all known competing compounds for each $\zeta$16–$\zeta$17 pair. In nitrogen, enthalpies are again referenced directly to $\zeta$18-N$\zeta$19 through

$\zeta$20

These are formally different constructions, but all are used to decide whether a $\zeta$21 phase is equilibrium, near-stable, or metastable under specified conditions (Elatresh et al., 2023, Nykiel et al., 25 Aug 2025, Goncharov et al., 17 Apr 2026).

Finite-temperature and kinetic effects are often decisive. The oxygen study explicitly combines PBE molecular dynamics, HSE06 thermodynamic perturbation, and vibrational entropy at 300 K, changing the preferred $\zeta$22 assignment from a near-degeneracy at 0 K to a specific Pm equilibrium phase above about 111 GPa. In Al–Ag, by contrast, the enthalpic preference for ordering into $\zeta$23 is offset by configurational entropy, while formation also depends on vacancy flux and the pre-existence of sufficiently large $\zeta$24 GP zones. A plausible implication is that the same $\zeta$25 label can refer either to an equilibrium phase boundary or to a kinetic waypoint embedded in a transformation sequence (Elatresh et al., 2023, Zhang et al., 2017).

Experimental discrimination is correspondingly multimodal. Oxygen relies on agreement among static structure search, X-ray diffraction peak positions, Raman/IR mode multiplicity, and superconducting $\zeta$26. Nitrogen uses single-crystal X-ray diffraction plus Raman signatures, especially additional low-frequency modes and vibron multiplicity, to distinguish $\zeta$27 from $\zeta$28. Boron required single-crystal synchrotron X-ray diffraction under laser-heated diamond-anvil-cell conditions to establish the Cmce $\zeta$29-Ga framework. In the Al–Ag alloy, atomic-resolution STEM and in situ annealing establish both structure and transformation path. These cases show that no single signature is sufficient across all materials; the operational meaning of a $\zeta$30 phase is fixed by a combination of crystallography, spectroscopy, and free-energy analysis (Elatresh et al., 2023, Goncharov et al., 17 Apr 2026, Chuvashova et al., 2017, Zhang et al., 2017).

6. Distinctions, misconceptions, and outstanding problems

An important distinction is between a $\zeta$31 phase and a constant zeta-potential. In suspensions of colloidal spheres, Smallenburg et al. studied phase diagrams under the boundary condition

$\zeta$32

within a Poisson–Boltzmann cell model. The resulting renormalized charge $\zeta$33 and effective screening $\zeta$34 feed a DLVO/Yukawa description whose phase behavior includes fluid, bcc, and fcc domains, as well as re-entrant melting driven by density- and salt-dependent discharge. This is not a crystallographic $\zeta$35 phase; it is a charging protocol whose phase consequences are nonetheless mathematically well defined (Smallenburg et al., 2010).

Within materials science proper, the literature shows several distinct stability classes. In oxygen, Pm is predicted to be the equilibrium room-temperature $\zeta$36 phase above approximately 111 GPa. In Al–Ag, $\zeta$37 is a metastable intermediate between $\zeta$38 GP zones and $\zeta$39. In tellurium, $\zeta$40 is the most stable few-layer form up to a crossover thickness of four layers, with the crossover shifted to larger thickness by doping. In carbides and nitrides, $\zeta$41 identifies a broad vacancy-ordered structural family with both stable and near-stable members. In nitrogen, $\zeta$42 is a polytype closely competing with $\zeta$43 itself. In boron, $\zeta$44 is a high-pressure allotrope reached only under extreme conditions (Elatresh et al., 2023, Zhang et al., 2017, Wang et al., 2018, Nykiel et al., 25 Aug 2025, Goncharov et al., 17 Apr 2026, Chuvashova et al., 2017).

Several unresolved points remain explicit in the reports. For solid oxygen, the Pm assignment still leaves discrepancies such as extra Raman modes and broad XRD peaks, which may arise from mixed-phase coexistence in the 96–124 GPa range and instrumental resolution; higher-pressure, higher-resolution experiments are specifically encouraged. For layered carbides and nitrides, the cited work states that open questions remain about structure, stability, and compositional pervasiveness. For nitrogen, the likely correspondence between $\zeta$45-N$\zeta$46 and the previously reported $\zeta$47-N$\zeta$48 reframes earlier phase identification. For Al–Ag, the recognition of $\zeta$49 as a vacancy-mediated intermediate implies that controlling vacancy flux or in situ ageing conditions may alter precipitation pathways. Taken together, these issues indicate that the term “zeta phase” often marks regions of dense polymorphism, subtle free-energy competition, and nontrivial kinetic trapping rather than settled crystallographic closure (Elatresh et al., 2023, Nykiel et al., 25 Aug 2025, Goncharov et al., 17 Apr 2026, Zhang et al., 2017).