pith. sign in

arxiv: 1906.08111 · v1 · pith:2Z2CVN3Vnew · submitted 2019-06-19 · ❄️ cond-mat.soft

Distinguishing noisy crystalline structures using bond orientational order parameters

Jan Haeberle , Matthias Sperl , Philip Born This is my paper

Pith reviewed 2026-05-25 20:00 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords bond orientational order parameterscrystal structure identificationneighborhood definitionnoise robustnesssoft matterFCCBCC
0
0 comments X

The pith

A simple change to the neighborhood definition makes bond orientational order parameters continuous and accurate for noisy crystal structures.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Bond orientational order parameters, introduced by Steinhardt and colleagues, characterize local particle arrangements in soft matter. Common ways to define a particle's neighbors produce ambiguous results when thermal noise or disorder is added, making it hard to tell FCC-based from BCC-based structures apart. The authors test these standard definitions and show they lose the ability to separate the structures under noise. They introduce one straightforward adjustment to how neighbors are selected. With this adjustment the order parameters stay continuous and recover the ability to identify the correct crystal type even in the presence of noise.

Core claim

The paper establishes that a straightforward modification to the neighborhood definition used when computing bond orientational order parameters produces values that remain continuous and allow reliable distinction between common crystal structures such as FCC- and BCC-based arrangements even after noise is introduced.

What carries the argument

The modified neighborhood definition applied to the calculation of bond orientational order parameters, which removes the ambiguity that arises with standard neighbor selections under noise.

If this is right

  • The order parameters become continuous functions of particle positions even when noise is present.
  • Standard crystal structures can be identified correctly without the failures observed with earlier neighbor choices.
  • The same parameters apply across a range of noise amplitudes without case-by-case retuning.
  • Local structure analysis in soft-matter simulations gains reliability when particles deviate from perfect lattice sites.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same neighborhood adjustment might be checked on additional lattices such as HCP or simple cubic to see whether the improvement generalizes.
  • Automated structure-recognition pipelines that currently discard noisy frames could retain more data if this definition is adopted.
  • The approach could be combined with existing order-parameter libraries to reduce misclassification rates in colloidal or granular experiments.

Load-bearing premise

The modified neighborhood definition works without introducing new artifacts or requiring separate adjustments for different noise levels and crystal types.

What would settle it

A direct test showing that the modified parameters still fail to separate FCC from BCC structures once moderate noise is added would falsify the central claim.

Figures

Figures reproduced from arXiv: 1906.08111 by Jan Haeberle, Matthias Sperl, Philip Born.

Figure 1
Figure 1. Figure 1: Radial distribution functions g(r) of a) BCC-based and b) FCC-bsed structures with noise of σ = 0.01d and σ = 0.036d as they are used in the analysis (see Sec. 2 for details of generating the structures). The distance r is scaled by the lattice constant d. The dashed lines indi￾cate the peak positions of the noise-free ideal lattice. The structures with noise still exhibit long-ranged positional correlatio… view at source ↗
Figure 2
Figure 2. Figure 2: The bond orientational order of the ideal noise [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Bond orientational order parameters of a BCC-based structure (orange) and a FCC-based structure (green) [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Bond orientational order parameter analysis of [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the average bond orientational order parameter attained from BCC- and FCC-based structures as [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

The bond orientational order parameters originally introduced by Steinhardt \emph{et. al.} [Phys. Rev. B \textbf{28}, 784 (1983)] are a common tool for local structure characterization in soft matter studies. Recently, Mickel \emph{et. al.} [J. Chem. Phys. \textbf{138}, 044501 (2013)] highlighted problems of the bond orientational order parameters due to the ambiguity of the underlying neighbourhood definition. Here we show the difficulties of distinguish common structures like FCC- and BCC-based structures with the suggested neighbourhood definitions when noise is introduced. We propose a simple improvement to the neighbourhood definition that results in robust and continuous bond orientational order parameters with which we can accurately distinguish crystal structures even when noise is present.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript addresses ambiguities in neighborhood definitions for bond orientational order parameters (BOOP) as identified by Mickel et al., showing that standard choices hinder distinction between FCC- and BCC-based structures under added noise. It proposes a simple modification to the neighborhood rule that is claimed to produce robust, continuous BOOP values enabling accurate structure identification even in noisy conditions.

Significance. If the proposed neighborhood modification performs as described across the tested cases, the work supplies a practical, low-overhead refinement to a standard tool in soft-matter simulations and experiments. The approach directly targets a documented limitation without introducing free parameters or fitted quantities, which strengthens its potential utility for reproducible local-structure analysis.

major comments (2)
  1. [Methods / §2] The central claim of 'accurate distinction' under noise rests on the new neighborhood definition, yet the manuscript provides neither an explicit equation nor algorithmic pseudocode for this rule (cf. abstract and §2). Without this, the reported continuity and robustness cannot be independently verified or reproduced.
  2. [Results] No quantitative metrics (classification accuracy, overlap integrals, or error bars versus noise amplitude) are reported for the distinction between structures; the results appear to rely on qualitative visual comparison of order-parameter distributions. This leaves the generality assumption (effectiveness across noise levels and crystal types without new artifacts) untested in a falsifiable manner.
minor comments (2)
  1. [Figures] Figure captions should explicitly state the noise amplitude and the precise neighborhood cutoff used for each panel to allow direct comparison with the text.
  2. [Results] A short table summarizing the order-parameter values (means and standard deviations) for each structure at representative noise levels would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the work's potential utility. We address each major comment below and will revise the manuscript to improve clarity and rigor.

read point-by-point responses

  1. Referee: [Methods / §2] The central claim of 'accurate distinction' under noise rests on the new neighborhood definition, yet the manuscript provides neither an explicit equation nor algorithmic pseudocode for this rule (cf. abstract and §2). Without this, the reported continuity and robustness cannot be independently verified or reproduced.

    Authors: We agree that the manuscript describes the neighborhood modification in prose but does not supply an explicit equation or pseudocode. This omission limits independent verification. In the revised manuscript we will add the precise mathematical definition of the modified neighborhood rule together with pseudocode in Section 2. revision: yes

  2. Referee: [Results] No quantitative metrics (classification accuracy, overlap integrals, or error bars versus noise amplitude) are reported for the distinction between structures; the results appear to rely on qualitative visual comparison of order-parameter distributions. This leaves the generality assumption (effectiveness across noise levels and crystal types without new artifacts) untested in a falsifiable manner.

    Authors: The present results are conveyed through visual comparison of order-parameter distributions. While these figures illustrate the improvement, we acknowledge that quantitative measures would make the claims more falsifiable. In the revision we will add classification accuracy, overlap integrals between distributions, and error bars plotted against noise amplitude for the tested crystal types. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper proposes an empirical improvement to the neighborhood definition for bond orientational order parameters to address noise-induced ambiguities noted in prior literature (Mickel et al.). No load-bearing step reduces a claimed prediction or result to a fitted input, self-definition, or self-citation chain by construction. The derivation relies on external benchmarks from Steinhardt et al. and Mickel et al. without importing uniqueness theorems or ansatzes from the authors' own prior work. The central claim remains independent and falsifiable against crystal structure data.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the effectiveness of an empirically chosen neighborhood rule whose performance is asserted but not derived from first principles or external benchmarks.

axioms (2)
  • standard math Bond orientational order parameters are computed from the Steinhardt et al. (1983) spherical-harmonic formulation
    The paper takes these definitions as given and modifies only the neighbor list construction.
  • domain assumption Noise is added to otherwise perfect crystalline lattices to simulate realistic conditions
    The abstract states that difficulties appear when noise is introduced.

pith-pipeline@v0.9.0 · 5661 in / 1208 out tokens · 29594 ms · 2026-05-25T20:00:51.120442+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [1]

    Allahyarov, K

    E. Allahyarov, K. Sandomirski, S. U. Egelhaaf, and H. L¨ owen. Crystallization seeds favour crystallization only during initial growth. Nature Communications, 6(1), November 2015

    work page 2015

  2. [2]

    Gasser, E

    U. Gasser, E. R. Weeks, A. Schofield, P. N. Pusey, and D. A. Weitz. Real-Space Imaging of Nucle- ation and Growth in Colloidal Crystallization. Science, 292(5515):258–262, April 2001

    work page 2001

  3. [3]

    J. H. Chu and Lin I. Direct observation of Coulomb crys- tals and liquids in strongly coupled rf dusty plasmas.Phys- ical Review Letters , 72(25):4009–4012, June 1994. Jan Haeberle et al.: Distinguishing noisy crystalline structures using bond orientational order parameters 7

    work page 1994

  4. [4]

    O. Arp, D. Block, A. Piel, and A. Melzer. Dust Coulomb Balls: Three-Dimensional Plasma Crystals. Physical Re- view Letters, 93(16), October 2004

    work page 2004

  5. [5]

    A. P. Nefedov et. al. PKE-Nefedov: plasma crystal experi- ments on the International Space Station. New Journal of Physics, 5(1):33, 2003

    work page 2003

  6. [6]

    Panaitescu and A

    A. Panaitescu and A. Kudrolli. Epitaxial growth of ordered and disordered granular sphere packings. Physical Review E, 90(3), September 2014

    work page 2014

  7. [7]

    T. Aste, M. Saadatfar, A. Sakellariou, and T. J. Senden. Investigating the geometrical structure of disordered sphere packings. Physica A: Statistical Mechanics and its Applications, 339(1-2):16–23, August 2004

    work page 2004

  8. [8]

    P. J. Steinhardt, D. R. Nelson, and M. Ronchetti. Bond- orientational order in liquids and glasses. Physical Review B, 28(2):784–805, July 1983

    work page 1983

  9. [9]

    Lechner and C

    W. Lechner and C. Dellago. Accurate determination of crystal structures based on averaged local bond or- der parameters. The Journal of Chemical Physics , 129(11):114707, September 2008

    work page 2008

  10. [10]

    Mickel, S

    W. Mickel, S. C. Kapfer, G. E. Schr¨ oder-Turk, and K. Mecke. Shortcomings of the bond orientational or- der parameters for the analysis of disordered particulate matter. The Journal of Chemical Physics , 138(4):044501, January 2013

    work page 2013

  11. [11]

    E Schr¨ oder-Turk, R

    G. E Schr¨ oder-Turk, R. Schielein, S. C Kapfer, F. M Schaller, G. W Delaney, T. Senden, M. Saadatfar, T. Aste, and K. Mecke. Minkowski tensors and local structure met- rics: Amorphous and crystalline sphere packings. In AIP Conference Proceedings, volume 1542, pages 349–352. AIP, 2013

    work page 2013

  12. [12]

    P. N. Pusey. The effect of polydispersity on the crystal- lization of hard spherical colloids. Journal de Physique , 48(5):709–712, 1987

    work page 1987

  13. [13]

    P. R. ten Wolde, M. J. Ruiz-Montero, and D. Frenkel. Numerical Evidence for bcc Ordering at the Surface of a Critical fcc Nucleus. Physical Review Letters, 75(14):2714– 2717, October 1995

    work page 1995

  14. [14]

    Rein ten Wolde, M

    P. Rein ten Wolde, M. J. RuizMontero, and D. Frenkel. Numerical calculation of the rate of crystal nucleation in a Lennard-Jones system at moderate undercooling. The Journal of Chemical Physics , 104(24):9932–9947, June 1996

    work page 1996

  15. [15]

    Volkov, M

    I. Volkov, M. Cieplak, J. Koplik, and J. R. Banavar. Molecular dynamics simulations of crystallization of hard spheres. Physical Review E , 66(6), December 2002

    work page 2002

  16. [16]

    Desgranges and J

    C. Desgranges and J. Delhommelle. Crystallization mech- anisms for supercooled liquid Xe at high pressure and temperature: Hybrid Monte Carlo molecular simulations. Physical Review B , 77(5), February 2008

    work page 2008

  17. [17]

    Rietz, C

    F. Rietz, C. Radin, H. L. Swinney, and M. Schr¨ oter. Nu- cleation in Sheared Granular Matter. Physical Review Let- ters, 120(5), February 2018

    work page 2018

  18. [18]

    F. A. Lindemann. ¨Uber die Berechnung molekularer Eigen- frequenzen. Phys. Z. , 11:609–612, 1910

    work page 1910

  19. [19]

    R. W. Cahn. Materials science: Melting from within. Na- ture, 413(6856):582, 2001

    work page 2001

  20. [20]

    C. Rycroft. Voro++: a three-dimensional Voronoi cell li- brary in C++. Technical Report LBNL-1432E, Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States), January 2009

    work page 2009

  21. [21]

    S. Fortune. Voronoi diagrams and Delaunay triangula- tions. In Computing in Euclidean Geometry , volume 4 of Lecture Notes Series on Computing , pages 225–265. World Scientific Publishing Co. Pte. Ltd, 1995

    work page 1995

  22. [22]

    J. P. Troadec, A. Gervois, and L. Oger. Statistics of voronoi cells of slightly perturbed face-centered cubic and hexag- onal close-packed lattices. Europhysics Letters, 42(2):167, 1998

    work page 1998