REVIEW 3 major objections 4 minor 45 references
Symmetry breaking in Prussian Blue Analogues via growth--guided local ordering of hexacyanometallate vacancies
T0 review · 3 major / 4 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read Growth-guided vacancy ordering reduces the local symmetry of cubic Prussian Blue analogue crystals to tetragonal.
desk verdict The diffuse-scattering evidence for local tetragonal vacancy order in a cubic-average Prussian blue is solid and new; the growth-guided mechanism is plausible but the axis-to-growth mapping is under-documented. 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 central object is the anisotropic pair-correlation model of vacancy ordering. Quantitatively, it is the three-dimensional difference pair distribution function (3D-DeltaPDF) refined from single-crystal diffuse scattering, which yields pair probabilities of Co(CN)_6 occupancy at each lattice vector; the key result is the contrast between correlation length of about 18 unit cells along [001] and about 8 along [100] and [010]. The mechanism is carried by a Monte Carlo Hamiltonian with a nearest-neighbour repulsion J1 enforcing electroneutrality and a split second-neighbour interaction J2,x = J2,y = -0.0135 versus J2,z = -0.11, whose anisotropy produces the observed tetragonal diffuse scattering.
What would settle it
Measure diffuse scattering from a Mn[Co]-PBA crystal whose growth direction is known independently, for example by in-situ growth observation or by facet indexing: if the long-correlation axis does not coincide with the growth direction, or if a crystal grown along [111] shows no trigonal diffuse pattern, the mechanism fails. Alternatively, anneal a freshly grown crystal and check whether the diffuse-scattering anisotropy disappears; if the correlations become cubic, the ordering was not locked in during growth.
Extended reading notes
Core claim
The central discovery is that in Mn[Co]-Prussian Blue Analogues the Co(CN)_6 vacancies are non-randomly arranged: they avoid each other at first neighbours and correlate positively along [001], with a correlation length of about 18 unit cells along the growth direction versus only about 8 unit cells perpendicular to it. This anisotropy gives the local structure tetragonal 4/mmm Laue symmetry while Bragg diffraction shows the average structure is cubic m-3m, and it manifests optically as birefringence. The authors identify the long-correlation axis with the growth direction and show in Monte Carlo simulations that a single anisotropic second-neighbour interaction (J2,z about six times J2,xy) reproduces the observed diffuse scattering. Growing crystals with {111} faces produces additional domains with trigonal (bar-3m) local symmetry. The general claim is that any disordered crystal grown irreversibly locks in growth-direction-specific defect correlations, providing a route to symmetry engineering by controlling growth direction.
Load-bearing premise
The long-correlation axis of the vacancies is identified with the crystal growth direction from the cube morphology, and the authors assume that once deposited, vacancies cannot rearrange; if that axis identification or that locking-in is wrong, the 'growth-guided' conclusion does not follow even though the tetragonal local order itself is real.
Editorial extensions
If this is right
- If growth direction is controlled, the local symmetry, and thus the optical, transport, and mechanical properties, of disordered crystals can be engineered without changing their composition.
- Cube-shaped Mn[Co]-PBA crystals are six pyramidal tetragonal twins sharing a single cubic average structure, with defect-rich twin boundaries at the cube edges that explain the previously observed preferential etching.
- Pore connectivity in these Prussian Blue analogues becomes anisotropic, with vacancy channels connecting better within planes perpendicular to the growth direction than along it.
- Growing along a <111> direction should produce trigonal (bar-3m) local domains, which the optical data already indicate and which future diffuse-scattering tomography could confirm.
- The mechanism should generalize to other out-of-equilibrium disordered crystals, such as electrochemically grown or flux-grown materials, provided post-growth rearrangement is suppressed.
Reading between the lines
- The same mechanism may underlie unexplained birefringence or anisotropy in other nominally cubic framework materials where disorder is high but diffuse scattering has not been measured.
- A direct test of the growth-guided claim would be to grow the same material along two different directions and compare the symmetry of the diffuse scattering; if the local symmetry tracks the growth direction, the mechanism is confirmed.
- If the locking-in hypothesis is correct, the correlation length along the growth direction could serve as a readout of growth kinetics, distinguishing layer-by-layer from island growth modes.
- The principle may extend well beyond insoluble salts to molecular and colloidal crystals whenever rearrangement is slow compared with the growth timescale, potentially making growth-direction control a general symmetry-engineering tool.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a combined single-crystal diffuse scattering, 3D-ΔPDF, and Mueller polarimetry study of Mn[Co]-Prussian Blue Analogues. It shows that although the average structure refines as cubic Fm\bar{3}m, the diffuse scattering from an optically selected pyramidal domain is tetragonal, with vacancy correlations decaying much more slowly along one cubic axis (~18 unit cells) than along the perpendicular axes (~8 unit cells). A Monte Carlo model with anisotropic second-neighbor interactions reproduces the pattern. The authors interpret this as 'Growth-Guided Local Ordering': the slow-correlation axis is the growth direction, the vacancy order is locked in during layer-by-layer growth, and the local Laue symmetry is reduced to 4/mmm while the average structure remains cubic. Growth along <111> is claimed to yield trigonal domains. The central experimental finding is the anisotropic local vacancy order; the causal growth mechanism is an interpretation that goes beyond the final-state data.
Significance. If the growth-direction attribution can be established, this is a significant result: it would provide a mechanism for symmetry control in disordered crystals and a concrete experimental demonstration in an important materials family. The paper has real strengths: the anisotropic diffuse scattering is observed directly and is not an artifact of the Monte Carlo model; the 3D-ΔPDF refinement separates substitutional and displacive correlations; the average-structure analysis is careful; Mueller polarimetry provides independent optical evidence; and the data and scripts are made available. The Monte Carlo model is transparent and explicitly labeled as approximate. The main weakness is that the causal 'growth-guided' and 'locked-in' claims are not uniquely supported by the static final-state measurements, because the relation between the diffuse-scattering axis and the external crystal morphology is not experimentally documented.
major comments (3)
- [Real-Space Model / Methods] The assignment of the long-correlation direction to the crystal growth direction is not experimentally documented. The Methods describe cutting a pyramidal domain under a polarizing microscope and gluing it to a loop, followed by indexing with XDS/Dials, but no step records how the optical or morphological orientation of the cut domain was transferred to the diffractometer coordinate frame. Because the average structure is cubic, indexing does not fix the physical orientation of the unique diffuse-scattering axis relative to the crystal habit: any cubic crystal can be reindexed so that the slow axis is [001]. Therefore the sentence 'The [001] direction, characterized by longer vacancy correlations, corresponds to the crystal growth direction' is an assumption, not a measured relation. This assumption is load-bearing for the 'Growth-Guided' mechanism claimed in the title and abstract. A co-registration experiment (for example, mounting a crystal with a known face on the diffractometer, or comparing several crystals whose growth axes are known from morphology) is needed; without it, the causal wording should be softened to a hypothesis.
- [Discussion and outlook] The claim that 'the insoluble nature of the material ensures that once the building blocks attach during growth, they do not rearrange, thus preserving the growth-induced anisotropy' is asserted without any kinetic test. No annealing, aging, dissolution/regrowth, or growth-rate variation experiment is reported. The static diffuse scattering data are equally consistent with anisotropic short-range order established during growth and with an equilibrium or post-growth reorganization that also produces a tetragonal local axis. To support the locking-in step of the mechanism, the paper needs an experimental control (for example, measuring diffuse scattering after thermal treatment, or comparing crystals grown under conditions that change the growth rate or post-growth time). Without such a test, the abstract's claim that the reduced symmetry 'persists in the final structure' through a growth process is not uniquely established.
- [Abstract / Symmetry Control Through Growth Direction] The statement that growth along [111] 'produces domains with trigonal symmetry' overstates what is shown. The Symmetry Control section acknowledges that the small size and overlapping domains 'currently prevent us from performing detailed diffuse scattering analysis'; the evidence for the {111} domains is Mueller polarimetry only, and the trigonal symmetry is inferred from a subgroup analysis ('likely produces domains'). The abstract should be phrased in terms of consistency with trigonal local symmetry, or the corresponding diffuse scattering measurement should be provided.
minor comments (4)
- [3D-ΔPDF Analysis / Figure 2c] The quoted correlation lengths of approximately 18 and 8 unit cells are presented without uncertainties or a discussion of how resolution and background subtraction affect the decay profiles; since the ratio is used as a quantitative measure of anisotropy, please add error estimates or explicitly state that these are approximate values.
- [Throughout] There are several typographical errors: 'margnianlly' should be 'marginally' (Average Structure section), 'untwined' should be 'untwinned' (Methods), 'withing' should be 'within' (Real-Space Model), 'the the scaling' should be 'the scaling' (Methods), 'Data avalibility' should be 'Data availability' (heading), and 'Views and options' should probably be 'Views and opinions' (Acknowledgements).
- [Figure 2 caption] The parenthetical expression '(h+1, k+1, 0, with h,k = odd)' is awkward and should be rewritten for clarity.
- [Methods / Monte Carlo modelling] The mean-field derivation of the Monte Carlo interaction parameters and the displacive correlation analysis are attributed to reference 23, a 'Manuscript in preparation (2025)'; because these items are needed to reproduce the model, please provide a preprint, a deposition, or an expanded description in the Methods.
Circularity Check
No significant circularity: the tetragonal local ordering is read directly from measured single-crystal diffuse scattering and 3D-ΔPDF refinement, not generated by the Monte Carlo model; the only self-referential element is a non-load-bearing deferral of mean-field parameter details to an unpublished companion manuscript.
full rationale
The paper's central result, that Mn[Co]-PBA retains cubic m-3m average symmetry while the vacancy distribution has tetragonal 4/mmm local symmetry, is established from direct experimental observations: single-crystal diffuse scattering shows tetragonal intensity distributions (Figure 2a), and the 3D-ΔPDF refinement yields anisotropic vacancy-vacancy correlations with longer range along the unique axis (Table S6 and Figure 2b-c). These data are not produced by the Monte Carlo model; rather, the model is fitted to them, as the text states that "to reproduce the observed diffuse scattering pattern, we found that the strength of J2,z must be set approximately six times of J2,xy", and the Methods explicitly disclaim any pretense of an optimal model ("the aim is not to represent an optimal model"). The causal attribution to growth rests on the empirical identification of the long-correlation axis with the growth direction and on a qualitative layer-by-layer argument, which is an interpretation rather than a derivation from the model. The only self-citation of note is reference 23, an in-preparation manuscript by overlapping authors, for the mean-field derivation of the Monte Carlo interaction parameters and for the displacive-correlation analysis; this citation is not load-bearing because neither the tetragonal symmetry conclusion nor the observed correlation anisotropy depends on those parameters. No equation in the paper reduces a predicted quantity to a fitted input, and no uniqueness theorem or ansatz is imported from prior same-author work. The main substantive weakness is evidentiary, not circular: the mapping of the diffuse-scattering long axis to the morphological growth direction is asserted rather than documented, and the locking-in mechanism is not tested by annealing or time-resolved experiments.
Assumptions & free parameters
free parameters (4)
- J1 =
0.605
- J2,x = J2,y =
-0.0135
- J2,z =
-0.11
- Monte Carlo temperature T =
1.3
assumptions (5)
- domain assumption Bragg diffraction measures only the cubic average structure; diffuse scattering is needed to see local ordering.
- domain assumption The diffuse signal can be cleanly separated into substitutional and displacive correlations, and the substitutional part is unbiased by the punch-and-fill removal of Bragg peaks.
- domain assumption The measured tetragonal c-axis corresponds to the crystal growth direction identified from morphology.
- domain assumption Post-growth rearrangement is negligible because the material is insoluble.
- domain assumption Mean-field approximation for the interaction parameters, with details in an unpublished manuscript.
Cite this review
Pith. "Pith review of Symmetry breaking in Prussian Blue Analogues via growth--guided local ordering of hexacyanometallate vacancies." pith.science (2026). https://pith.science/paper/FJOEAGPK
@misc{pith2026250205936,
author = {Pith},
title = {Pith review of: Symmetry breaking in Prussian Blue Analogues via growth--guided local ordering of hexacyanometallate vacancies},
year = {2026},
howpublished = {\url{https://pith.science/paper/FJOEAGPK}},
note = {Machine review of arXiv:2502.05936}
}
abstract
We report Growth--Guided Local Ordering, a novel mechanism of symmetry reduction in disordered crystals. This mechanism operates through the directional ordering of point defects during crystal growth, where defect correlations develop preferentially along the growth direction, resulting in reduced symmetry that persists in the final structure through the spatial distribution of defects. We demonstrate this phenomenon in Mn[Co]-Prussian Blue Analogues, disordered cyanide crystals containing numerous Co(CN)$_6$ vacancies. Single crystal diffuse scattering reveals pronounced anisotropy in vacancy distribution: strong correlations along [001] growth direction contrast with weak correlations perpendicular to it. This local ordering reduces the Laue symmetry to tetragonal $4/mmm$, evident in properties such as birefringence, while the average structure retains cubic $m\bar 3m$ symmetry. When growth proceeds along [111] direction, the same mechanism produces domains with trigonal symmetry. Because this mechanism relies on fundamental aspects of crystal growth rather than specific material properties, it offers a general strategy for symmetry control in disordered crystals. Crucially, it transforms the complex task of altering crystal symmetry into the more manageable challenge of controlling growth direction, achievable through various established techniques such as the use of surfactants during crystallization.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
Nye, J. F. Physical properties of crystals: their representation by tensors and matrices (Oxford Univer- sity Press, 1985)
work page 1985
-
[2]
Nonlinear optics on ferroic materials (John Wiley & Sons, 2023)
Fiebig, M. Nonlinear optics on ferroic materials (John Wiley & Sons, 2023)
work page 2023
-
[3]
Bloom, B. P., Paltiel, Y ., Naaman, R. & Waldeck, D. H. Chemical Reviews 124, 1950–1991 (2024)
work page 2024
-
[4]
V ojta, M.Advances in Physics 58, 699–820 (2009)
work page 2009
- [5]
-
[6]
& Jona-Lasinio, G.Physical Review 122, 345 (1961)
Nambu, Y . & Jona-Lasinio, G.Physical Review 122, 345 (1961)
work page 1961
-
[7]
Higgs, P. W. Physical Review Letters 13, 508 (1964)
work page 1964
- [8]
Show all 45 references
-
[9]
Annual Review of Materials Science 28, 1–27 (1998)
Goodenough, J. Annual Review of Materials Science 28, 1–27 (1998). 8/19
1998
-
[10]
Benedek, N. A. & Hayward, M. A. Annual Review of Materials Research 52, 331–355 (2022)
2022
-
[11]
Bostr"om, H. L. CrystEngComm 22, 961–968 (2020)
2020
-
[12]
L., Senn, M
Bostr"om, H. L., Senn, M. S. & Goodwin, A. L. Nature Communications 9, 2380 (2018)
2018
-
[13]
W., Southon, P
Chapman, K. W., Southon, P. D., Weeks, C. L. & Kepert, C. J. Chemical Communications 3322–3324 (2005)
2005
-
[14]
K., Thallapally, P
Motkuri, R. K., Thallapally, P. K., McGrail, B. P. & Gho- rishi, S. B. CrystEngComm 12, 4003–4006 (2010)
2010
-
[15]
Y ., Rodríguez-García, B
Goberna-Ferrón, S., Hernández, W. Y ., Rodríguez-García, B. & Galán-Mascarós, J. R. ACS Catal. 4, 1637–1641 (2014)
2014
-
[16]
Guo, S. et al. Industrial & Engineering Chemistry Research 59, 13831–13840 (2020)
2020
-
[17]
& Goodwin, A
Cattermull, J., Pasta, M. & Goodwin, A. L. Materials Horizons 8, 3178–3186 (2021)
2021
-
[18]
J., Pasta, M
Cattermull, J., Roth, N., Cassidy, S. J., Pasta, M. & Good- win, A. L. Journal of the American Chemical Society 145, 24249–24259 (2023)
2023
-
[19]
& Hashimoto, K
Sato, O., Iyoda, T., Fujishima, A. & Hashimoto, K. Science 272, 704–705 (1996)
1996
-
[20]
S., Mathonière, C
Aguilà, D., Prado, Y ., Koumousi, E. S., Mathonière, C. & Clérac, R. Chemical Society Reviews 45, 203–224 (2016)
2016
-
[21]
& Ohkoshi, S.-i
Stefánczyk, O. & Ohkoshi, S.-i. Chem. Eur. J. 25, 15963– 15977 (2019)
2019
-
[22]
Simonov, A. et al. Nature 578, 256–260 (2020)
2020
-
[23]
& Arkadiy, S
Kholina, Y ., Chernyshov, D., Warren, M. & Arkadiy, S. Manuscript in preparation (2025)
2025
-
[24]
& Zeng, Y
Hu, M., Jiang, J.-S. & Zeng, Y . Chemical Communications 46, 1133–1135 (2010)
2010
-
[25]
Welberry, T. R. Diffuse X-ray Scattering and Models of Disorder (Oxford University Press, 2022)
2022
-
[26]
Welberry, T., Miller, G. H. & Pickard, D. Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 367, 175–192 (1979)
1979
-
[27]
W., Quintel, A
Hulliger, J., Roth, S. W., Quintel, A. & Bebie, H. Journal of Solid State Chemistry 152, 49–56 (2000)
2000
-
[28]
& Fl"orsheimer, M
Rechsteiner, P., Hulliger, J. & Fl"orsheimer, M. Chemistry of Materials 12, 3296–3300 (2000)
2000
-
[29]
Hulliger, J. et al. The crystallization of polar, channel- type inclusion compounds: Property-directed supramolec- ular synthesis (1997)
1997
-
[30]
Roth, S. W. et al. Advanced Materials 10, 1543–1546 (1998)
1998
-
[31]
& Chernyshov, D
Dyadkin, V ., Pattison, P., Dmitriev, V . & Chernyshov, D. Journal of Synchrotron Radiation 23, 825–829 (2016)
2016
-
[32]
MiTeGen, LTD, info@mitegen.com. P.O. Box 3867, Ithaca, NY 14852
-
[33]
Acta Crystallographica Section D: Biological Crystallography 66, 125–132 (2010)
Kabsch, W. Acta Crystallographica Section D: Biological Crystallography 66, 125–132 (2010)
2010
-
[34]
Winter, G. et al. Acta Crystallographica Section D: Structural Biology 74, 85–97 (2018)
2018
-
[35]
Meerkat: program for reciprocal space reconstruction
Simonov, A. Meerkat: program for reciprocal space reconstruction. https://github.com/aglie/meerkat (2019)
2019
-
[36]
Sheldrick, G. M. Acta Crystallographica Section A: Foundations of Crystallography 64, 112–122 (2008)
2008
-
[37]
V ., Bourhis, L
Dolomanov, O. V ., Bourhis, L. J., Gildea, R. J., Howard, J. A. & Puschmann, H. Journal of Applied Crystallography 42, 339–341 (2009)
2009
-
[38]
Gorfman, S. et al. Journal of Applied Crystallography 54, 914–923 (2021)
2021
-
[39]
& Zhang, N
Gorfman, S., Spirito, D., Zhang, G., Detlefs, C. & Zhang, N. Acta Crystallographica Section A: Foundations and Advances 78, 158–171 (2022)
2022
-
[40]
& Gorfman, S
Biran, I., Bosak, A., Ye, Z.-G., Levin, I. & Gorfman, S. Applied Physics Letters 124 (2024)
2024
-
[41]
& Steurer, W.Journal of Applied Crystallography 47, 1146–1152 (2014)
Simonov, A., Weber, T. & Steurer, W.Journal of Applied Crystallography 47, 1146–1152 (2014)
2014
-
[42]
& Simonov, A.Zeit
Weber, T. & Simonov, A.Zeit. Krist. (2012)
2012
-
[43]
& Dušek, M
Petˇríˇcek, V ., Palatinus, L., Plášil, J. & Dušek, M. Zeitschrift für Kristallographie-Crystalline Materials 238, 271–282 (2023)
2023
-
[44]
Chipman, R., Lam, W. S. T. & Young, G.Polarized Light and Optical Systems. CRC Press (2018)
2018
-
[45]
thermal diffuse scattering-like
Lu, S.-Y . & Chipman, R. A. JOSA A 13, 1106–1113 (1996). Data avalibility The raw data along with the scripts used to produce the figures in this paper are shared on the follow- ing doi:10.3929/ethz-b-000719385, https://www.research- collection.ethz.ch/handle/20.500.11850/7193...
1996 doi
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.