Zirconolite-2M: Structure and Controversies

• Zirconolite-2M is a crystalline polytype (CaZrTi2O7) used as a standard in nuclear waste immobilization studies, despite ongoing crystallographic controversies.
• Recent synchrotron diffraction and DFT analyses reveal that the ambient phase may be reclassified as triclinic zirconolite-2TR rather than the classic monoclinic model.
• Defect studies and diffusion simulations highlight its varied transport properties and the significant impact of radiation damage on both crystalline and amorphous states.

Searching arXiv for zirconolite-2M and closely related zirconolite studies to ground the article in the cited literature. Zirconolite-2M is the conventional name for the ambient crystalline zirconolite polytype commonly associated with the composition $\mathrm{CaZrTi_2O_7}$ in nuclear-waste, defect-physics, and high-pressure studies. In the literature summarized here, however, the status of that designation is not uniform. Several studies investigate “zirconolite” or monoclinic zirconolite without explicitly assigning the 2M polytype (Yang et al., 2015, Mulroue et al., 2012, Diver et al., 2020), whereas a recent high-pressure diffraction and DFT study argues that the ambient structure previously assigned as monoclinic zirconolite-2M is better described as a triclinic phase, termed zirconolite-2TR, with zirconolite-2M then corresponding to a pressure-induced $C2/c$ structure above $14.7$ GPa (Errandonea et al., 5 Aug 2025). Accordingly, Zirconolite-2M occupies a dual role in current research: it remains the canonical reference phase in much waste-form literature, but its exact ambient-pressure crystallographic assignment has become an active subject of revision.

1. Definition, nomenclature, and scope

Zirconolite is treated throughout the cited literature as a technologically important titanate phase with nominal composition $\mathrm{CaZrTi_2O_7}$, especially in the context of SYNROC and actinide immobilization (Yang et al., 2015, Mulroue et al., 2012, Diver et al., 2020). The label “zirconolite-2M” is not universal in these papers. One molecular-dynamics study of diffusion explicitly states that it does not use the label “zirconolite-2M,” although the crystalline zirconolite reference is reasonably interpreted as the standard ordered zirconolite used in atomistic modeling (Yang et al., 2015). A first-principles defect study likewise investigates monoclinic zirconolite derived from Rossell’s X-ray structure but does not explicitly state “zirconolite-2M” (Mulroue et al., 2012). A radiation-damage study on crystalline and amorphous zirconolite also does not provide enough crystallographic detail to verify that its crystal is specifically zirconolite-2M (Diver et al., 2020).

This terminological asymmetry is significant. In older and phase-generic work, “zirconolite” often functions as a practical shorthand for the common monoclinic reference structure. By contrast, the 2025 high-pressure study explicitly compares zirconolite-2M, -4M, -3O, and -3T, and reinterprets ambient-pressure zirconolite-2M as a triclinic structure, $P\bar{1}$, which the authors rename zirconolite-2TR (Errandonea et al., 5 Aug 2025). This suggests that “zirconolite-2M” can mean either the traditional ambient reference in applied literature or, in stricter crystallographic usage, the monoclinic $C2/c$ member of a pressure-dependent structural sequence.

The ideal formula $\mathrm{CaZrTi_2O_7}$ is retained across these studies even when chemistry is modified. Examples include a model with 10% U substitution on the Zr site in diffusion simulations (Yang et al., 2015), Nd-doped precursor glasses designed for zirconolite crystallization (Loiseau et al., 2011), and highly non-ideal meteoritic zirconolite compositions enriched in Nb, Fe, Si, Th, U, and Pb (Sakuma et al., 16 Mar 2026). Thus, zirconolite-2M is best understood as the polytype-specific endpoint of a broader zirconolite structural-chemical family rather than an isolated fixed composition.

2. Crystal-structure assignment and recent reinterpretation

The most explicit crystallographic treatment of zirconolite-2M is the synchrotron XRD and DFT study of high-pressure structural stability (Errandonea et al., 5 Aug 2025). That work begins from the previously assigned monoclinic $C2/c$ model for zirconolite-2M, with unit-cell parameters $a = 12.4413(2)$ Å, $b = 7.2700(1)$ Å, $C2/c$0 Å, $C2/c$1, and $C2/c$2, and notes that one quarter of the Ti atoms occupy a Wyckoff $C2/c$3 position with 50% statistical occupation. Rietveld refinement of this model at $C2/c$4 GPa gives $C2/c$5, $C2/c$6, and $C2/c$7, with weak reflections that cannot be indexed by $C2/c$8 (Errandonea et al., 5 Aug 2025).

Lowering the symmetry through a $C2/c$9 group-subgroup relation yields a substantially better fit. The triclinic model gives $14.7$0, $14.7$1, and $14.7$2, explains weak low-angle reflections, and removes the problematic partial Ti occupancy because all Ti atoms fully occupy their Wyckoff positions (Errandonea et al., 5 Aug 2025). The authors therefore propose that the ambient-pressure phase should be called zirconolite-2TR rather than zirconolite-2M. At $14.7$3 GPa, the triclinic structure is assigned to space group $14.7$4 with $14.7$5 Å, $14.7$6 Å, $14.7$7 Å, $14.7$8, $14.7$9, $\mathrm{CaZrTi_2O_7}$0, and $\mathrm{CaZrTi_2O_7}$1 (Errandonea et al., 5 Aug 2025).

The same study describes the structure in terms of planes of corner-sharing CaO$\mathrm{CaZrTi_2O_7}$2 and ZrO$\mathrm{CaZrTi_2O_7}$3 polyhedra running perpendicular to [001], separated by hexagonal tungsten bronze-type Ti layers containing TiO$\mathrm{CaZrTi_2O_7}$4 octahedra and TiO$\mathrm{CaZrTi_2O_7}$5 trigonal bipyramids (Errandonea et al., 5 Aug 2025). This preserves the familiar zirconolite framework motif while replacing the older disordered monoclinic average structure with a lower-symmetry ordered variant.

DFT strengthens that reinterpretation. The previously assigned ambient-pressure monoclinic $\mathrm{CaZrTi_2O_7}$6 structure shows imaginary phonon frequencies, whereas the proposed triclinic $\mathrm{CaZrTi_2O_7}$7 structure has all positive phonon branches and lower enthalpy below 28 GPa (Errandonea et al., 5 Aug 2025). A plausible implication is that much of the historical use of “zirconolite-2M” in ambient-pressure discussions may correspond, structurally, to what this study reclassifies as zirconolite-2TR.

3. Relationship to monoclinic zirconolite in defect and waste-form studies

Despite the reinterpretation just noted, several earlier studies remain directly relevant to Zirconolite-2M because they analyze monoclinic or conventional crystalline zirconolite as the parent ordered state from which defect production, amorphization, and transport are discussed. The first-principles intrinsic-defect study states that zirconolite has a monoclinic crystal structure taken from Rossell’s X-ray diffraction work and describes rows of 6-fold and 4-fold coordinated Ti–O polyhedra separated by alternating Ca and Zr layers (Mulroue et al., 2012). Although the paper does not explicitly name zirconolite-2M, that structural description is fully consistent with the common monoclinic zirconolite form generally associated with 2M.

The defect topology is crystallographically consequential. The idealized structure contains distinct Ti environments and oxygen sites of differing connectivity, including 2-fold and 3-fold coordinated oxygen (Mulroue et al., 2012). These low-coordinate oxygen sites later become central to defect physics because holes associated with cation vacancies localize on them (Mulroue et al., 2012). The same study also notes that experimental studies often report a 5-fold coordinated Ti absent in the ideal stoichiometric structure, and that such 5-fold Ti appears naturally in the calculations when Ti/Zr antisite-like substitutions are introduced (Mulroue et al., 2012). This suggests that polyhedral multiplicity observed experimentally may partly reflect off-stoichiometry rather than the ideal ordered crystal alone.

The diffusion study treats the crystalline reference more schematically. It does not provide lattice parameters, space group, coordination numbers, or polytype distinctions such as 2M, 3O, or 3T (Yang et al., 2015). Instead, it frames crystalline zirconolite as an ordered network with equivalent interatomic distances and local diffusion pathways, in contrast to the amorphous state, which has a broad distribution of interatomic distances and angles (Yang et al., 2015). For Zirconolite-2M specifically, the key point is methodological: much of the atomistic waste-form literature uses a standard ordered zirconolite crystal as the reference, even when the polytype is not named explicitly.

The precursor-glass study is even more indirect. It addresses Nd-doped glasses designed for preparation of zirconolite $\mathrm{CaZrTi_2O_7}$8-based glass-ceramics and shows that the parent glass contains Zr-centered short-range order with similarities to zirconolite, including $\mathrm{CaZrTi_2O_7}$9-fold oxygen coordination at $P\bar{1}$0 Å and likely second neighbors around $P\bar{1}$1 Å (Loiseau et al., 2011). However, it does not establish the crystallized phase as zirconolite-2M. Its relevance lies in precursor chemistry and nucleation rather than polytype identification.

4. Intrinsic defects and local relaxation modes

The DFT study of intrinsic defects provides the most detailed microscopic account of the defect landscape relevant to monoclinic zirconolite and, by inference, to the conventional Zirconolite-2M reference (Mulroue et al., 2012). Calculations were performed with VASP using PBE, PAW pseudopotentials, a $P\bar{1}$2 supercell containing 176 ions, a 600 eV plane-wave cutoff, and $P\bar{1}$3-point sampling corresponding to $P\bar{1}$4 minimum reciprocal-space sampling (Mulroue et al., 2012). Stable interstitial sites were located using ab initio random structure search with 100 random insertions per interstitial species (Mulroue et al., 2012).

A central result is the pronounced site and charge-state dependence of both vacancies and interstitials. Neutral oxygen interstitials are structurally flexible: in the 88-ion cell the lowest-energy configuration is a $P\bar{1}$5 dumbbell with O–O separation $P\bar{1}$6, whereas in the larger 176-ion cell the extra oxygen increases two 6-fold Ti to 7-fold and converts the 4-fold Ti environments mostly to 5- or 6-fold (Mulroue et al., 2012). The singly charged oxygen interstitial binds to a 4-fold Ti polyhedron to create a defective 5-fold coordinated Ti environment (Mulroue et al., 2012). Ca and Ti interstitials preferentially occupy the $P\bar{1}$7 channels running through zirconolite, while the lowest-energy Zr interstitial configuration is a substitutional-Zr plus Ti-interstitial complex in which added Zr displaces a lattice 6-fold Ti and the displaced Ti enters the $P\bar{1}$8 channel (Mulroue et al., 2012). This indicates an intrinsic tendency toward Ti–Zr cation disordering under defect accommodation.

Vacancy relaxations are especially rich. Neutral Zr vacancies form an $P\bar{1}$9-like molecule with O–O bond length $C2/c$0, close to free $C2/c$1, while neutral and singly charged vacancies in one Ti environment also form $C2/c$2 molecules with O–O distances $C2/c$3 and $C2/c$4 (Mulroue et al., 2012). Charged Zr vacancies do not form $C2/c$5, but instead convert neighboring Ti polyhedra and delocalize holes over low-coordinate oxygen (Mulroue et al., 2012). Ca vacancies leave holes localized on 2-fold oxygen sites and do not generate gap states (Mulroue et al., 2012). Collectively, these results show that oxygen rebonding, Ti polyhedral reconfiguration, and oxygen-site hole trapping are basic structural response modes of the crystal.

The Frenkel pair energies quantify defect susceptibility. Reported values are $C2/c$6, $C2/c$7, $C2/c$8, and $C2/c$9 eV for different oxygen Frenkel pairs, $\mathrm{CaZrTi_2O_7}$0 eV for Ca Frenkel, $\mathrm{CaZrTi_2O_7}$1 eV for Zr Frenkel, and $\mathrm{CaZrTi_2O_7}$2, $\mathrm{CaZrTi_2O_7}$3, and $\mathrm{CaZrTi_2O_7}$4 eV for the three Ti Frenkel environments (Mulroue et al., 2012). The oxygen defects therefore show the strongest dependence on local chemical environment, and some oxygen Frenkel pairs are very inexpensive (Mulroue et al., 2012). This supports a picture in which oxygen-sublattice disorder is comparatively easy to produce and may contribute materially to radiation-induced amorphization of the crystalline parent phase.

5. Transport, amorphization, and radiation-damage evolution

A major theme in zirconolite research is that the long-term behavior relevant to waste immobilization cannot be inferred from the pristine crystal alone. In classical MD simulations of solid-state diffusion, the authors compared crystalline zirconolite, high-density amorphous zirconolite with the same density as the crystal, and low-density amorphous zirconolite with 5% lower density (Yang et al., 2015). Diffusion coefficients were extracted from the Einstein relation,

$\mathrm{CaZrTi_2O_7}$5

and Arrhenius fits used

$\mathrm{CaZrTi_2O_7}$6

The central result is that amorphization alone causes a major increase in diffusion, even at constant density (Yang et al., 2015).

In the crystal, only oxygen diffuses on accessible MD timescales; in amorphous zirconolite, all species diffuse: O, Ca, Zr, Ti, and U (Yang et al., 2015). Oxygen has fitted activation energy $\mathrm{CaZrTi_2O_7}$7 eV in the crystal and $\mathrm{CaZrTi_2O_7}$8 eV in amorphous zirconolite at crystal density, while cation diffusion in amorphous zirconolite has activation energies $\mathrm{CaZrTi_2O_7}$9 eV for Ca, $C2/c$0 eV for Zr, $C2/c$1 eV for Ti, and $C2/c$2 eV for U (Yang et al., 2015). Reducing amorphous density by 5% does not appreciably change these activation energies within error, but increases the Arrhenius prefactor for cations by factors of about 2–3, whereas oxygen is only weakly affected (Yang et al., 2015). The interpretation is that disorder creates a distribution of local migration pathways, with faster or softer pathways dominating net transport even when average density is unchanged (Yang et al., 2015).

A later large-scale cascade study extends that reasoning from static amorphization to continued irradiation of the amorphous state (Diver et al., 2020). Using DL_POLY_4, a 13,365,000-atom cell, a Buckingham potential mixed with ZBL at short range, and repeated overlapping 70 keV cascades in the same region, the authors compared crystalline and melt-quenched amorphous zirconolite (Diver et al., 2020). The amorphous model was produced by melting at 6000 K for 100 ps, cooling to 300 K at 10 K ps$C2/c$3, and relaxing for 100 ps at 300 K in NPT, yielding a 4% increase in volume relative to the crystal (Diver et al., 2020).

The headline result is that amorphous zirconolite is “softer” than crystalline zirconolite, with overall about 4 times more atoms displaced during a cascade than in the crystal (Diver et al., 2020). Oxygen is the most displaced species in both phases; Ca becomes the second most displaced species in the amorphous state (Diver et al., 2020). The crystal shows an early-time recovery peak around $C2/c$4 ps and is not fully amorphized after 5 cascades in the overlap region, whereas the amorphous state shows no early-time recovery peak and continues to evolve structurally under further irradiation (Diver et al., 2020). Local density mapping in 10 Å cubes reveals regions with up to 20% lower density in irradiated amorphous zirconolite (Diver et al., 2020). The authors suggest that such density inhomogeneities may become important for future developments in nuclear waste storage, and a plausible implication is that the irradiated amorphous state constitutes its own evolving structural regime rather than a terminal, inert disorder state.

6. Chemical flexibility, precursor states, and natural occurrences

Zirconolite’s technological and geological significance derives in large part from its chemical flexibility. In the Nd-doped glass study, the precursor composition belongs to the $C2/c$5 system, with base glass composition $C2/c$6, $C2/c$7, $C2/c$8, $C2/c$9, $a = 12.4413(2)$0, and $a = 12.4413(2)$1, plus up to 2.2 mol.% (10 wt.%) $a = 12.4413(2)$2 (Loiseau et al., 2011). These glasses had already been shown to lead to zirconolite crystallization after nucleation and crystal-growth treatments (Loiseau et al., 2011). Zr occupies a well-defined local site in the glass, with $a = 12.4413(2)$3, $a = 12.4413(2)$4 Å, and likely Ca/Ti and Zr second neighbors around $a = 12.4413(2)$5 Å (Loiseau et al., 2011). The authors explicitly state that this short-range order “could predispose to zirconolite nucleation” (Loiseau et al., 2011).

Ti also shows local environments compatible with zirconolite. ESR detects two types of Ti$a = 12.4413(2)$6 sites; the authors interpret the spectra in terms of a compressed octahedral symmetry highly distorted to yield $a = 12.4413(2)$7, corresponding either to $a = 12.4413(2)$8 or $a = 12.4413(2)$9 symmetry and thus 5- to 6-fold coordination (Loiseau et al., 2011). They note that these Ti coordinations are close to those encountered in zirconolite (Loiseau et al., 2011). Nd behaves differently: EXAFS gives $b = 7.2700(1)$0, $b = 7.2700(1)$1 Å, and $b = 7.2700(1)$2, with no clearly identified second neighbors (Loiseau et al., 2011). The authors conclude that Nd$b = 7.2700(1)$3 cannot act as a nucleating agent for zirconolite crystallization and may hinder nucleation with increasing diffusionnal difficulties (Loiseau et al., 2011). This underscores that zirconolite-forming short-range order in precursor materials resides mainly in the Ti–Zr framework, not necessarily in the rare-earth or actinide site chemistry.

Natural zirconolite shows still broader chemical departure from the ideal end-member. In the andesitic meteorite Erg Chech 002, zirconolite occurs as needle-, fiber-, and low-aspect-ratio crystals, with the largest grains about 30 $b = 7.2700(1)$4m long and $b = 7.2700(1)$5m wide, mainly in interstitial settings associated with albitic plagioclase, pyroxene, silica, merrillite, baddeleyite, ilmenite, troilite, Fe metal, and unidentified Ti-Zr-bearing phases (Sakuma et al., 16 Mar 2026). The grains are strongly non-ideal, with oxide ranges including CaO 7.97–12.62 wt%, TiO$b = 7.2700(1)$6 19.69–36.64 wt%, ZrO$b = 7.2700(1)$7 27.92–36.00 wt%, Nb$b = 7.2700(1)$8O$b = 7.2700(1)$9 5.02–15.44 wt%, ThO$C2/c$00 0.60–6.86 wt%, and UO$C2/c$01 0.17–2.11 wt% (Sakuma et al., 16 Mar 2026). The paper does not assign a polytype; Raman was featureless, interpreted as recoil-damage-related, so a specific zirconolite-2M identification is not possible (Sakuma et al., 16 Mar 2026).

Nevertheless, the meteoritic case shows that zirconolite can be a robust U–Pb chronometer. Fourteen NanoSIMS spot analyses yield a weighted mean $C2/c$02 ratio of $C2/c$03, corresponding to an age of $C2/c$04 Ma (Sakuma et al., 16 Mar 2026). The authors interpret this age as dating shock metamorphism of the parent asteroid’s crust rather than primary crystallization (Sakuma et al., 16 Mar 2026). For Zirconolite-2M research, this is relevant mainly as evidence that zirconolite-group minerals can preserve chronometric information even when structural polytype resolution is unavailable.

7. Technological significance and open questions

In waste-form science, zirconolite is valued because it can incorporate actinides and remains structurally and chemically consequential over repository-relevant timescales (Yang et al., 2015, Mulroue et al., 2012, Diver et al., 2020). One diffusion study states that under typical waste loading zirconolite is expected to become amorphous after about 1000 years, whereas required operational lifetimes are of order 100,000 to 1,000,000 years (Yang et al., 2015). The practical inference is explicit: predictive models should use diffusion data appropriate to amorphous zirconolite rather than only pristine crystalline zirconolite (Yang et al., 2015). At the same time, cascade simulations show that amorphous zirconolite is not structurally frozen but continues to evolve under irradiation, developing coordination changes and local density inhomogeneities not evident in total pair distribution functions (Diver et al., 2020).

High-pressure work adds a different dimension of structural stability. The triclinic ambient-pressure phase proposed as zirconolite-2TR transforms reversibly at 14.7 GPa to a monoclinic $C2/c$05 phase, with a 1% volume discontinuity and no evidence of amorphization or decomposition up to 30.5 GPa (Errandonea et al., 5 Aug 2025). The low-pressure phase has $C2/c$06 Å$C2/c$07, $C2/c$08 GPa, and $C2/c$09, while the high-pressure monoclinic phase has $C2/c$10 Å$C2/c$11, $C2/c$12 GPa, and $C2/c$13 (Errandonea et al., 5 Aug 2025). Compression is only slightly anisotropic in the low-pressure phase, with principal compressibilities $C2/c$14, $C2/c$15, and $C2/c$16 along $C2/c$17, $C2/c$18, and $C2/c$19, respectively (Errandonea et al., 5 Aug 2025). These data support substantial structural resilience under compression, even as they complicate the historical definition of Zirconolite-2M.

The principal controversy is therefore crystallographic rather than functional. Older defect, diffusion, and waste-form studies rely on a conventional crystalline zirconolite reference that is very likely the common monoclinic 2M form in practical terms, but they generally do not resolve polytype identity explicitly (Yang et al., 2015, Mulroue et al., 2012, Diver et al., 2020). The recent high-pressure and DFT evidence, by contrast, indicates that ambient material long called zirconolite-2M may actually be triclinic zirconolite-2TR (Errandonea et al., 5 Aug 2025). This suggests that future work on defects, transport, radiation damage, and substitution chemistry may need to distinguish more sharply between the historical “2M reference structure” used in atomistic waste-form modeling and the revised crystallographic assignment proposed from modern synchrotron and lattice-dynamical analysis.

A plausible implication is that the most urgent open problem is not whether zirconolite remains a viable waste host, since the cited studies consistently support its relevance, but how precisely its ordered ambient structure should be parameterized when linking crystallography to defect energetics, amorphization pathways, transport coefficients, and long-term performance models.