REVIEW 2 major objections 3 minor 41 references
Structural Chirality from Short-Range Order in Heteroanionic Materials
T0 review · 2 major / 3 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The ordering of two anion species over an achiral parent lattice can by itself force structural chirality; for ReO3-type oxyfluorides such as NbO2F, every configuration satisfying two short-range ordering rules is chiral, with no atomic…
desk verdict A genuinely new configurational mechanism for chirality, with a solid combinatorial proof and a credible but surrogate-limited NbO2F demonstration. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing device is a combinatorial encoding: each cation with cis-MX4Y2 coordination has exactly one fluoride-free octahedral axis, so the entire 3×3×3 configuration can be coded by cation labels (x, y, or z). Period-three OOF ordering is then equivalent to demanding that each <001> chain contain exactly one cation labelled with the chain's direction, which turns the problem into a depth-first search over 10,752 valid labellings that reduce to twelve chiral orbits. The finite-temperature demonstration is carried by two order parameters—inter-chain coherence, which measures whether chain patterns hold a fixed relative phase, and the pseudoscalar chirality $\chi = |E_+|^2 - |E_-|^2$, built from the wavevector $\mathbf{k} = (1/3,1/3,1/3)$ Fourier amplitudes of the three anion sublattices—with the first-order transition located by the equal-weight coexistence point and Binder cumulant of $\chi$.
What would settle it
A single-crystal or electron-diffraction experiment on NbO2F below about 494 K could settle it: the predicted chiral phase sharpens the {hk1/3}* superlattice features from diffuse sheets into spots, so a sample that keeps diffuse sheets at low temperature would rule out the predicted equilibrium chiral ground state.
Extended reading notes
Core claim
The central claim is a structural theorem and its chemical demonstration. In an achiral cubic parent with a 2:1 anion ratio, the rules 'each cation sees two minority anions cis' and 'each <001> chain follows a period-three X–X–Y pattern' jointly remove every mirror, inversion, and rotoinversion from the stabiliser of every allowed configuration; the twelve symmetry-inequivalent configurations that satisfy both rules on the 3×3×3 cell all belong to Sohncke space groups (one P3_121, one P3_1, three C2, seven P1). Neither rule is chiral alone, and the result survives any displacive subgroup of the parent. In NbO2F the most coherent of these chiral configurations, P3_121 with Nb–F–Nb helices winding about <111>, is the ground state by DFT and by cluster-expansion searches, and the equilibrium chiral phase persists to a first-order transition at about 494 K. The ordering that sets in below the transition is the relative phase of the period-three chains; the chirality is an improper order parameter that appears only when that interlacing locks in.
Load-bearing premise
The chemical conclusion rests on the assumption that both short-range ordering rules are genuinely preferred in the real crystal and that the DFT-trained cluster-expansion model, accurate to about 1.4 meV/atom, correctly ranks configurations separated by as little as 2.2 meV/atom.
Editorial extensions
If this is right
- Any ReO3-type MX2Y system in which both ordering rules are active has the same twelve chiral orbits; TaO2F, TiOF2, and oxynitride perovskites such as SrTaO2N are named in-scope candidates.
- The chiral phase of NbO2F is the equilibrium state up to about 494 K; above that temperature the material keeps cis coordination and OOF chains but loses the fixed phase relationship between chains, so the chirality vanishes.
- Since chirality is an improper order parameter triggered by inter-chain coherence, it cannot be made to order on its own; any sample that is not coherently chain-interlaced will be achiral even if both local rules hold.
- The two phases are distinguishable in diffraction: {hk1/3}* superlattice features appear as diffuse sheets in the achiral phase and sharpen toward spots once the chains lock, providing a direct experimental fingerprint.
Reading between the lines
- Editorial: the design rule this suggests is broader than the paper's examples—any heteroanionic compound with two species, 2:1 stoichiometry, and an octahedral framework that enforces both cis preference and period-three chain preference should be forced chiral, so the mechanism could be screened purely from these two local energetic biases.
- Editorial: because the transition is first order, growing the chiral phase likely requires annealing below the transition temperature to let chains re-phase; quenching from high temperature would trap the achiral state, and two enantiomeric helices may form domains, which would affect how optical activity is measured in practice.
- Editorial: the combinatorial core—two jointly imposed rules eliminating all improper operations—might extend to other parent sublattices and stoichiometries; a search for pairs of local rules whose simultaneous imposition kills every mirror in the stabiliser could reveal further configurational chiral mechanisms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript argues that structural chirality in inorganic crystals can arise purely from configurational anion ordering, without displacive symmetry breaking. For ReO3-type MX2Y materials with cis-MX4Y2 coordination and period-three [X–X–Y] chain ordering, the authors use a depth-first search to enumerate all configurations on the 3×3×3 supercell, obtaining 10,752 labellings that reduce to twelve symmetry orbits, and prove that each orbit belongs to a Sohncke space group. For NbO2F, direct DFT energies place the maximally symmetric P3_121 chiral orbit lowest among the twelve, and a cluster expansion trained on 96 DFT configurations is used in Monte Carlo searches and Wang–Landau simulations, yielding a predicted first-order chiral-to-achiral transition at Tc = 494 ± 1 K. The paper concludes that configurational ordering is a chemically realizable route to structural chirality and identifies heteroanionic materials as candidates for chiral functional materials.
Significance. The combinatorial theorem is a genuinely new and clean result: it shows that two local anion-ordering rules force chirality on the smallest commensurate cell, and the one-to-one labelling argument in Appendix A is convincing and machine-checkable in spirit. This mechanism is distinct from displacive chirality and is of clear interest to the heteroanionic-materials community. The paper also provides explicit falsifiable predictions (ground-state space group, first-order transition, diffraction signatures) and deposits code and data, which strengthens its reproducibility. However, the significance of the NbO2F demonstration depends on the accuracy of the cluster expansion, whose cross-validated RMSE (1.4 meV/atom) is comparable to the energy spread among the twelve chiral orbits (2.2 meV/atom) and to the latent heat of the transition (≈3.5 meV/atom). If the CE mis-ranks low-energy achiral configurations, the predicted P3_121 ground state and Tc could change, and the 'chemically realizable route' claim would be weakened, even though the theorem itself would remain valid.
major comments (2)
- [Section IV C, Methods B] The global ground-state claim for NbO2F rests on a cluster expansion with cross-validated RMSE of 1.4 meV/atom, while the twelve chiral orbits span only 2.2 meV/atom (Fig. 6). The search over 6×6×6 supercells for lower-energy achiral configurations uses this surrogate model, and the stopping rule in Methods B ('stopped predicting configurations lower in energy than these') is self-referential: a CE that systematically mis-ranks the low-energy basin will not predict lower-energy configurations even if they exist. Since the combinatorial theorem applies only to the 3×3×3 supercell, the paper's assertion that P3_121 is the ground state ('on three counts') depends entirely on the CE search to exclude achiral cis+OOF configurations on larger cells. I recommend adding direct DFT energies for the lowest-energy achiral configurations identified by the CE on larger supercells, or for a targeted set of plausible achiral cis+OOF configurations, to close this gap. Without such a validation, the central claim that configurational ordering is a chemically realizable route to chirality in NbO2F is not fully supported.
- [Section V, Methods D] The reported Tc = 494 ± 1 K is the statistical uncertainty from the finite-size extrapolation only; it does not include any contribution from the cluster-expansion model error. The transition is driven by inter-chain phase ordering, with energy differences of order a few meV/atom (the latent heat is about 3.5 meV/atom), which is the same order as the CE's 1.4 meV/atom RMSE. The ARDR posterior over ECIs provides a natural way to propagate this uncertainty into the density of states and Tc; at minimum, a sensitivity analysis (e.g., perturbing ECIs within their posterior and recomputing Tc) is needed to know whether the prediction '494 K' is meaningful. The current statement '494±1 K' overstates the precision of a model whose underlying energy errors are an order of magnitude larger than the energy scale of the transition.
minor comments (3)
- [Figure 4c] The Binder cumulant minima for L=15 and L=18 are omitted from the figure to keep the smaller sizes visible; please include them in a supplementary figure or in an inset, so the reader can verify the deepening trend that is claimed in the text.
- [Section IX C] The list of Monte Carlo trial moves includes 'cell reflections' that map a configuration to its mirror image. Since the paper emphasizes that the chiral phase is incompatible with mirror symmetry, please clarify whether such moves are used only as unphysical equilibration aids and how they are reconciled with detailed balance when the two enantiomers are distinct thermodynamic states.
- [Abstract and Section II] The term 'OOF' is used freely after the abstract but is only defined in Section II; please define it at first use in the abstract or introductory paragraph.
Circularity Check
The enumeration theorem is self-contained, but the Monte-Carlo 'confirmation' of the P3_121 ground state restates the cluster-expansion training stop rule.
-
fitted input called prediction
[Section IV C (Ground-state search) and Section IX B (Cluster expansion)]
"This continued until the expansion both reproduced the DFT ordering of the low-lying configurations (which include the twelve chiral orbits) and stopped predicting configurations lower in energy than these. ... P3 121 is therefore the ground state on three counts: the enumeration admits only the twelve orbits, the density-functional energies place it lowest, and the Monte Carlo search finds nothing beneath it."
The cluster expansion was iteratively trained until it stopped predicting configurations lower in energy than the twelve chiral orbits. The later 'Monte Carlo search finds nothing beneath it' therefore restates that stopping condition rather than providing an independent test: the surrogate was selected to satisfy exactly the property being verified. This does not touch the Section III enumeration theorem or the direct DFT ranking of the twelve orbits, but it means the claimed exclusion of lower-energy achiral configurations on 6x6x6 cells is a fitted-input confirmation, not an independent prediction.
full rationale
The central combinatorial claim is not circular: Appendix A gives a one-to-one labelling-to-configuration mapping, and the stabiliser computation for each of the twelve orbits is a direct symmetry calculation with no fitted parameters. The DFT energies of the twelve chirality orbits are first-principles data, and the statement that P3_121 is the lowest of those twelve is independent of the cluster expansion. The finite-temperature transition at T_c = 494 K is also a genuine model prediction: no experimental chirality temperature is used as input, and the quoted uncertainty is a posterior over finite-size extrapolation. The one circular element is the ground-state search: the cluster expansion was trained until it both reproduced the DFT ordering of the chiral orbits and predicted no lower-energy configuration, and the later Monte-Carlo quenches are then cited as a third count that 'finds nothing beneath it.' That is a fitted input being reported as a prediction. It weakens the global-ground-state and Tc chemical-demonstration claims, though not the theorem; hence a partial circularity score of 6.
Assumptions & free parameters
free parameters (2)
- Cluster expansion effective cluster interactions (ECIs) =
Not enumerated numerically in the paper
- Cluster expansion hyperparameters (pair/triplet/quadruplet cutoffs and ARDR relevance threshold) =
9 Å, 5 Å, 5 Å; threshold not stated
assumptions (5)
- standard math A crystal structure is chiral if and only if its space group is Sohncke, containing no improper symmetry operations.
- domain assumption The ReO3-type framework is treated as a fixed achiral scaffold with anion sites occupied by two species; enumeration ignores atomic displacements.
- domain assumption The local ordering rules cis-MX4Y2 and period-three OOF chains are simultaneously active and sufficient to describe NbO2F anion order.
- domain assumption PBEsol DFT and the fitted cluster expansion accurately resolve energy differences of order a few meV/atom among anion configurations.
- domain assumption First-order finite-size scaling Tc(L) = Tc(infinity) + a/L^3 applies, and the replica-exchange Wang-Landau walkers are ergodic over anion configuration space.
Cite this review
Pith. "Pith review of Structural Chirality from Short-Range Order in Heteroanionic Materials." pith.science (2026). https://pith.science/paper/24NPM76R
@misc{pith2026260804841,
author = {Pith},
title = {Pith review of: Structural Chirality from Short-Range Order in Heteroanionic Materials},
year = {2026},
howpublished = {\url{https://pith.science/paper/24NPM76R}},
note = {Machine review of arXiv:2608.04841}
}
abstract
Established routes to structural chirality in inorganic crystals depend on symmetry-lowering atomic displacements of an achiral parent structure. We show that chirality can instead be forced by the ordering of two anion species over the sites of an achiral parent, independent of atomic displacements. ReO$_3$-type oxyfluorides with a 2:1 anion stoichiometry, such as NbO$_2$F, exhibit two ordering patterns: cis octahedral coordination and period-three anion-chain ordering, both rooted in the off-centring of $d^0$ cations. By direct enumeration, we prove that every configuration combining both orderings on the smallest commensurate cell is chiral, belonging to a Sohncke space group. For the specific case of NbO$_2$F, density-functional theory calculations predict that the ground state is the maximally symmetric chiral configuration. Wang-Landau Monte Carlo simulations of a DFT-trained cluster-expansion model predict an equilibrium chiral phase up to a first-order transition at 494 K. These results establish configurational ordering as a chemically realisable route to structural chirality, and mark heteroanionic materials with analogous ordering chemistry as candidates for a new class of chiral functional materials.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
A. M. Glazer and K. Stadnicka, On the origin of optical activity in crystal structures, J. Appl. Crystallogr.19, 108 (1986)
work page 1986
-
[2]
E. Bousquet, M. Fava, Z. Romestan, F. G´ omez-Ortiz, E. E. McCabe, and A. H. Romero, Structural chirality and related properties in periodic inorganic solids: Re- view and perspectives, J. Phys.: Condens. Matter37, 163004 (2025)
work page 2025
-
[3]
F. G´ omez-Ortiz, A. Zabalo, A. M. Glazer, E. E. McCabe, A. H. Romero, and E. Bousquet, Structural chirality and natural optical activity across theαtoβphase transi- tion in SiO 2 and AlPO 4 from first principles, J. Appl. Crystallogr.59, 225 (2026)
work page 2026
-
[4]
F. G´ omez-Ortiz, A. H. Romero, and E. Bousquet, Path- ways to crystal chirality: An algorithm to identify displacive chiral phase transitions, Phys. Rev. B112, 014107 (2025)
work page 2025
-
[5]
N. Charles, R. J. Saballos, and J. M. Rondinelli, Struc- tural diversity from anion order in heteroanionic materi- als, Chem. Mater.30, 3528 (2018)
work page 2018
-
[6]
Neither of the two ordering rules produces chirality on its own
and sevenP1. Neither of the two ordering rules produces chirality on its own. OOF ordering alone leaves mirror-symmetric ar- rangements available (a mirror perpendicular to a chain through its fluoride maps such a configuration onto it- self), and cis coordination alone likewise leaves arrange- ments with mirrors that exchange the two fluorides at a catio...
- [7]
-
[8]
T. Yamada and H. Koizumi, Czochralski growth of Ag4P2O7 crystals, J. Cryst. Growth64, 558 (1983)
work page 1983
Show all 41 references
-
[9]
M. V. Talanov, V. B. Shirokov, and V. M. Talanov, Anion order in perovskites: A group-theoretical analysis, Acta Crystallogr. A72, 222 (2016)
2016
-
[10]
F. J. Brink, R. L. Withers, and L. Nor´ en, An elec- tron diffraction and crystal chemical investigation of oxy- gen/fluorine ordering in niobium oxyfluoride, NbO 2F, J. Solid State Chem.166, 73 (2002)
2002
-
[11]
R. L. Withers, F. Brink, Y. Liu, and L. Nor´ en, Cluster chemistry in the solid state: Structured diffuse scatter- ing, oxide/fluoride ordering and polar behaviour in tran- sition metal oxyfluorides, Polyhedron26, 290 (2007)
2007
-
[12]
Dabachi, M
J. Dabachi, M. Body, C. Galven, F. Boucher, and C. Leg- ein, Preparation-dependent composition and O/F order- ing in NbO2F and TaO2F, Inorg. Chem.56, 5219 (2017)
2017
-
[13]
Legein, B
C. Legein, B. J. Morgan, A. G. Squires, M. Body, W. Li, M. Burbano, M. Salanne, T. Charpentier, O. J. Borkiewicz, and D. Dambournet, Correlated anion dis- order in heteroanionic cubic TiOF 2, J. Am. Chem. Soc. 146, 21889 (2024)
2024
-
[14]
L. K. Frevel and H. W. Rinn, The crystal structure of NbO2F and TaO 2F, Acta Crystallogr.9, 626 (1956)
1956
-
[15]
Kunz and I
M. Kunz and I. D. Brown, Out-of-center distortions around octahedrally coordinatedd 0 transition metals, J. Solid State Chem.115, 395 (1995)
1995
-
[16]
P. S. Halasyamani and K. R. Poeppelmeier, Noncen- trosymmetric oxides, Chem. Mater.10, 2753 (1998)
1998
-
[17]
Pauling, The principles determining the structure of complex ionic crystals, J
L. Pauling, The principles determining the structure of complex ionic crystals, J. Am. Chem. Soc.51, 1010 (1929)
1929
-
[18]
C. R. Morelock, B. K. Greve, M. Cetinkol, K. W. Chap- man, P. J. Chupas, and A. P. Wilkinson, Role of anion site disorder in the near zero thermal expansion of tan- talum oxyfluoride, Chem. Mater.25, 1900 (2013)
2013
-
[19]
M. Yang, J. Or´ o-Sol´ e, J. A. Rodgers, A. B. Jorge, A. Fuertes, and J. P. Attfield, Anion order in perovskite oxynitrides, Nat. Chem.3, 47 (2011)
2011
-
[20]
J. P. Attfield, Principles and applications of anion order in solid oxynitrides, Cryst. Growth Des.13, 4623 (2013)
2013
-
[21]
J. K. Harada, N. Charles, K. R. Poeppelmeier, and J. M. Rondinelli, Heteroanionic materials by design: Progress toward targeted properties, Adv. Mater.31, 1805295 (2019)
2019
-
[22]
D. A. Keen and A. L. Goodwin, The crystallography of correlated disorder, Nature521, 303 (2015)
2015
-
[23]
N. A. Benedek and C. J. Fennie, Hybrid improper fer- roelectricity: A mechanism for controllable polarization- magnetization coupling, Phys. Rev. Lett.106, 107204 (2011)
2011
-
[24]
Kresse and J
G. Kresse and J. Furthm¨ uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)
1996
-
[25]
Kresse and D
G. Kresse and D. Joubert, From ultrasoft pseudopoten- tials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)
1999
-
[26]
J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, Phys. Rev. Lett.100, 136406 (2008)
2008
-
[27]
˚Angqvist, W
M. ˚Angqvist, W. A. Mu˜ noz, J. M. Rahm, E. Fransson, C. Durniak, P. Rozyczko, T. H. Rod, and P. Erhart, ICET – a Python library for constructing and sampling alloy cluster expansions, Adv. Theory Simul.2, 1900015 (2019)
2019
-
[28]
Fransson, F
E. Fransson, F. Eriksson, and P. Erhart, Efficient con- struction of linear models in materials modeling and applications to force constant expansions, npj Comput. Mater.6, 135 (2020)
2020
-
[29]
Pedregosa, G
F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cour- napeau, M. Brucher, M. Perrot, and ´E. Duchesnay, Scikit-learn: Machine learning in Python, J. Mach. Learn. Res.12, ...
2011
-
[30]
Kleiven, J
D. Kleiven, J. Akola, A. A. Peterson, T. Vegge, and J. H. Chang, Training sets based on uncertainty estimates in the cluster-expansion method, J. Phys. Energy3, 034012 (2021)
2021
-
[31]
Wang and D
F. Wang and D. P. Landau, Efficient, multiple-range ran- dom walk algorithm to calculate the density of states, Phys. Rev. Lett.86, 2050 (2001)
2001
-
[32]
Wang and D
F. Wang and D. P. Landau, Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram, Phys. Rev. E64, 056101 (2001)
2001
-
[33]
Vogel, Y
T. Vogel, Y. W. Li, T. W¨ ust, and D. P. Landau, Generic, hierarchical framework for massively parallel Wang–Landau sampling, Phys. Rev. Lett.110, 210603 (2013)
2013
-
[34]
B. J. Morgan, mchammer-pt: Replica-exchange orches- trators for icet/mchammer ensembles,https://github. com/bjmorgan/mchammer-pt(2026), software
2026
-
[35]
R. E. Belardinelli and V. D. Pereyra, Fast algorithm to calculate density of states, Phys. Rev. E75, 046701 (2007)
2007
-
[36]
R. E. Belardinelli and V. D. Pereyra, Wang–Landau al- gorithm: A theoretical analysis of the saturation of the error, J. Chem. Phys.127, 184105 (2007)
2007
-
[37]
Binder, Finite size scaling analysis of ising model block distribution functions, Z
K. Binder, Finite size scaling analysis of ising model block distribution functions, Z. Phys. B: Condens. Matter43, 119 (1981)
1981
-
[38]
Binder and D
K. Binder and D. P. Landau, Finite-size scaling at first- order phase transitions, Phys. Rev. B30, 1477 (1984)
1984
-
[39]
B. J. Morgan, chainorder: Order parameters for chain- ordered mixed-anion structures,https://github.com/ bjmorgan/chainorder(2026), software
2026
-
[40]
B. J. Morgan, hofmann: Ball-and-stick crystal- structure figures with Matplotlib,https://github.com/ bjmorgan/hofmann(2026), software
2026
-
[41]
Beckett, J
G. Beckett, J. Beech-Brandt, K. Leach, Z. Payne, A. Simpson, L. Smith, A. Turner, and A. Whiting, ARCHER2 service description (2024)
2024
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.