Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

The crystallographic quaternions and their product law

T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Quaternion rotations can be multiplied directly in the crystal basis using a metric-tensor product law.

desk verdict Correct algebra, but the novelty is a repackaged inverse metric and the EBSD efficiency claim is not supported by the paper's own numbers. read the letter →

arxiv 2607.16899 v1 pith:SFZ6RJGG submitted 2026-07-18 cond-mat.mtrl-sci

Cyril Cayron This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords quaternionsrotationsmetrictensorcrosscrystalbasiselectronbackscatterdiffractionnon-Cartesiancoordinatescompositionlaw
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to show that the standard quaternion product law, which composes 3D rotations, can be rewritten so that it works directly in the fractional-coordinate basis of a crystal, without transforming to a Cartesian frame. The key move is to replace the ordinary scalar product with the metric tensor and the ordinary cross product with a new 'cross tensor' that equals sqrt(det M) times the inverse of the metric tensor. This yields a product formula that depends only on the lattice parameters and lets orientations be encoded and multiplied in the crystal's own coordinates. The author argues this could make electron backscatter diffraction (EBSD) orientation calculations faster and convention-free.

What carries the argument

The cross tensor X = sqrt(det(M)) M⁻¹, a linear map that turns the ordinary coordinate cross product into the geometric cross product in the crystal basis. It is the cross-product counterpart of the metric tensor and, together with M, allows the scalar and vector parts of the quaternion product to be evaluated entirely in non-Cartesian coordinates.

What would settle it

Take a triclinic cell, encode two rotations as crystallographic quaternions, multiply them by the proposed formula, convert the result to a Cartesian rotation matrix, and compare with the product of their Cartesian matrix forms; any mismatch would refute the product law.

Watch

Extended reading notes

Core claim

At the paper's center is a closed-form product law for rotations encoded as quaternions whose axes are written in the fractional coordinates of an arbitrary crystal cell. For two such quaternions q1 = c1 + s1 u~1 and q2 = c2 + s2 u~2, the composition is q3 = q1q2 = c3 + s3 u~3, with c3 = c1c2 − s1s2 (u~1ᵀ M u~2) and s3 u~3 = c2 s1 u~1 + c1 s2 u~2 + s1 s2 X (u~1 × u~2), where M is the metric tensor of the lattice and X = sqrt(det(M)) M⁻¹. Because X is a function of M alone, the whole composition rule is set by the lattice parameters. The paper argues that this removes the need to shuttle between crystal and Cartesian coordinates during rotation composition.

Load-bearing premise

The practical case for using these quaternions in EBSD rests on an unbenchmarked cost assumption: converting Kikuchi-derived data into Euler angles costs roughly 45 extra floating-point operations per pixel, and spherical indexing can supply orientations directly in crystal coordinates.

Editorial extensions

If this is right

  • EBSD orientation maps could be processed end to end in fractional coordinates, eliminating structure-tensor conversions per pixel.
  • Symmetry recognition in orientation maps simplifies because symmetry axes have coordinates 0 or ±1 in the crystal basis.
  • The same construction transplants to other rotation representations that rely on dot and cross products, such as compact axis-angle forms.
  • The computational cost per product is 46 FLOPs with the crystallographic formula, versus 28 for the Cartesian formula alone, but the paper argues the full pipeline favors the crystallographic route.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The FLOP advantage hinges on the unverified assumption that spherical indexing can deliver orientations directly in crystal coordinates; if that fails, the efficiency argument collapses even though the product law itself is unaffected.
  • Since X and M depend only on the phase, they can be precomputed once per material, making the per-product cost effectively that of a few matrix-vector products, a point the paper leaves implicit.
  • The construction suggests a general recipe: any rotation formula written with ordinary dot and cross products can be lifted to a non-Cartesian basis by substitution (dot becomes M-dot, cross becomes X-cross), potentially simplifying other crystallographic calculations.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper proposes a generalization of unit quaternions to non-orthogonal crystallographic bases. It introduces a 'cross tensor' X = [b×c, c×a, a×b], proves X = sqrt(det M) M^{-1} for the metric tensor M, and uses M and X to write the composition law for quaternions directly in the crystal basis (Eqs. 18-20), together with analogous Rodrigues and vector-rotation formulas (Eqs. 24-25). The final section argues that this construction could make EBSD orientation calculations more efficient by avoiding Cartesian structure-tensor round trips.

Significance. If the stated qualifications are added, the paper gives a correct and useful coordinate formula: the central identity X = sqrt(det M) M^{-1} and the product law (18)-(20) reproduce the standard rotation composition in fractional coordinates, and the derivation in the appendices is self-contained and not circular. The extension to Rodrigues vectors is a useful by-product. The practical EBSD claim, however, is not established, and some stated properties of the metric tensor are overstated; these issues do not invalidate the core algebra but require substantial revision of Section 4.

major comments (3)
  1. [§4 and Table 1] The efficiency comparison is not reliable. The paper counts a 3×3 symmetric matrix-vector product as 6M+3A, but symmetry reduces storage, not arithmetic; the operation count is 9M+6A. With correct counts, the crystallographic product costs at least 58 FLOPs, not the reported 46. In addition, item (b) states 'additional 45 FLOPs (9 M and 6 A)', which is arithmetically inconsistent (9+6=15), and even accepting the paper's counts the comparison is 46 vs 43 in favor of the Cartesian route. The claimed net advantage therefore rests entirely on the unsubstantiated assertion that spherical indexing returns orientations 'directly in the crystal basis'; the cited implementations return orientation matrices in an orthogonal crystal frame. This section needs a corrected benchmark or should be reframed as a qualitative discussion.
  2. [§4(d), §2.1] The statement that the metric tensor is 'independent of any convention' is incorrect. The matrix M is the Gram matrix of the chosen basis vectors (a,b,c); its entries change with the unit-cell setting, axis order, and handedness. What is convention-independent is the underlying bilinear form, not the matrix used in Eqs. (19)-(20). Since the product law is written in a particular crystal basis, it is basis-dependent in the same way as any coordinate formula. This overstatement should be corrected.
  3. [Eq. (16), §2.2] The identity X = sqrt(det M) M^{-1} is valid only for a right-handed basis. The derivation uses a* = (b×c)/V with V = det(a,b,c) and then sets V = sqrt(det M), conflating signed determinant with positive volume. For a left-handed basis, the correct cross tensor is X = -sqrt(det M) M^{-1}; the metric tensor alone does not encode the orientation of the basis. The manuscript should either explicitly assume the standard right-handed crystallographic setting or carry a signed volume throughout. Without this, the claim that the product law depends uniquely on the metric is not general.
minor comments (5)
  1. [§2.1] The text calls M the transformation 'from the reciprocal basis to the direct basis', but the derivation v* = M v shows that M maps direct-basis coordinates to reciprocal-basis coordinates, i.e. M = [B_c -> B_c^*]. Please correct the wording.
  2. [§4] The split of FLOP counts for the Cartesian quaternion product is inconsistent: '4 M and 4 A' plus '12 M and 6 A' sums to 16M+10A, not the stated 16M+12A. The correct split is 4M+3A (angular) and 12M+9A (vector).
  3. [Appendix 1] The sentence 'with c_i = cos(p_jk) and c_i = sin(p_jk)' should read 'and s_i = sin(p_jk)'.
  4. [Appendix 2] The text says 'u_ij = -u_ji' immediately after defining u_ij as an angle between vectors; if the angles are ordinary positive angles, the relation should be u_ij = u_ji. If signed angles are intended, this should be stated.
  5. [§3] The code is described as 'available upon request'; for reproducibility, a permanent repository or supplementary material would be preferable.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the crystallographic quaternion product is a coordinate-form version of the standard quaternion product, and the cross tensor is derived from the metric tensor.

full rationale

The central derivation is self-contained and not circular. Equations (18)-(20) are the standard quaternion composition law re-expressed in a non-orthonormal basis: the scalar product is written with the metric tensor (eq. 12) and the cross product with the cross tensor (eq. 13). The cross tensor itself is derived from the reciprocal-basis relation X = sqrt(det M) M^{-1} (eq. 16), not assumed. The Cartesian composition law is proved in Appendix 1 from spherical laws, and those laws are proved in Appendix 2 from the identity between the direct and reciprocal metric tensors. No fitted parameter is extracted from the target result, and no prediction is statistically forced. The only self-citation (Cayron, 2006) is used to label the composition-of-reflections construction as a groupoid law; the accompanying geometric proof does not depend on that citation, so it is not load-bearing. The Section 4 EBSD-efficiency argument does rest on an unverified external premise about spherical indexing returning orientations directly in the crystal basis, and its FLOP accounting contains an internal inconsistency (45 vs 9M+6A), but these are correctness or applicability concerns, not circularity. The core algebraic derivation therefore has no circular step; the score of 1 reflects only the minor, non-load-bearing self-citation.

Assumptions & free parameters 0 free parameters · 6 assumptions · 1 invented entities

No free parameters fitted to data — lattice parameters are physical inputs. The contribution rests on standard quaternion algebra and reciprocal-lattice identities; the only genuinely ad hoc premise is the EBSD FLOP accounting in Section 4. The 'cross tensor' is a repackaged scaled inverse metric.

assumptions (6)
  • standard math Unit quaternions represent rotations, and Hamilton's product law composes them (with the norm-multiplicativity property).
    Section 1.1: rotations encoded as q = c + sũ and composed by (2)-(4). Standard algebra, used as the starting point of the generalization.
  • domain assumption Reciprocal-basis identity: b×c = V a*, c×a = V b*, a×b = V c*, with V = det(S) = √det(M).
    Section 2.2 (eq. 14): the entire cross tensor rests on this standard crystallographic dual-basis relation; if the metric's determinant were not the cell volume, X = V M⁻¹ would fail.
  • domain assumption M = SᵀS and det(M) = V² for any choice of Cartesian structure tensor S; M is symmetric positive definite.
    Sections 2.1-2.2 (eqs. 11, 15): required to express the cross tensor and its determinant solely through lattice parameters.
  • standard math Spherical laws of cosines and sines, proved in Appendix 2 from metric identities (M_u* = M_u⁻¹).
    Appendix 1 uses these laws to prove the quaternion product law; Appendix 2's proof is standard spherical trigonometry expressed through metric tensors.
  • ad hoc to paper EBSD pipeline accounting: converting Kikuchi-derived orientation data to Euler angles costs ~45 FLOPs per pixel; spherical indexing yields orientations directly in the crystal basis.
    Section 4(b): estimated, unbenchmarked, and internally inconsistent in the stated M/A split (45 FLOPs but only 9 M + 6 A listed). The whole 'more effective for EBSD' conclusion depends on it.
  • domain assumption The metric tensor is 'independent of any convention' while the structure tensor is not.
    Sections 2.1 and 4(d): overstated — metric components depend on the chosen unit-cell setting; only the physical tensor is convention-invariant. This premise underlies the claimed software advantage.
invented entities (1)
  • Cross tensor X = [b×c, c×a, a×b] = √det(M) M⁻¹ independent evidence
    purpose: To compute the cross product of two vectors directly in a non-orthonormal crystal basis (eq. 13), the cross-product analogue of the metric tensor.
    Not a new physical entity: it is a renamed, scaled inverse metric. Independent handle: eq. (13) is checkable against standard crystallographic cross-product computations for any lattice, and it reduces to the textbook reciprocal-basis identity (eq. 14). The paper's claim of absence in the literature is unsubstantiated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The crystallographic quaternions and their product law." pith.science (2026). https://pith.science/paper/SFZ6RJGG

@misc{pith2026260716899,
  author       = {Pith},
  title        = {Pith review of: The crystallographic quaternions and their product law},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SFZ6RJGG}},
  note         = {Machine review of arXiv:2607.16899}
}
read the original abstract

Unit quaternions are widely used in science to encode rotations because the quaternion product is more efficient than matrix product and more stable than Rodrigues product to calculate the composition of two rotations. However, quaternions in their usual form refers to a Cartesian basis; they cannot be used in crystallography as they are. The usual way to solve this issue to apply back-and-forth coordinate changes from the crystal basis to a Cartesian basis attached to the crystal with the help of the structure tensor. Here, we show that actually quaternions can be used directly in the crystal basis by generalizing the quaternion product law. In that aim, we introduced a matrix that we called cross tensor. It allows the calculation of the cross product in the crystal basis, a bit like the metric tensor allows it for the scalar product. We also show the cross tensor is proportional to the inverse of the metric tensor. The formula of crystallographic quaternion product is then given; it depends uniquely on the metric tensor. The application of the crystallographic quaternions to Electron Back Scatter Diffraction is discussed.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. The crossmetric tensor and the geometrical meaning of the imaginary numbers

    cond-mat.mtrl-sci 2026-08 conditional novelty 3.0 of 10

    A unit quaternion is represented as a pair of oriented planes; for Cartesian bases i, j, k become pairs of cube faces, which is said to explain Hamilton's multiplication rules.

Reference graph

Works this paper leans on

3 extracted references · cited by 1 Pith paper

  1. [1]

    L., Wright, S

    Adams, B. L., Wright, S. I. & Kunze, K. (1993). Metall Trans A 24, 819–831. Britton, T. B., Jiang, J., Guo, Y., Vilalta-Clemente, A., Wallis, D., Hansen, L. N., Winkelmann, A. & Wilkinson, A. J. (2016). Materials Characterization 117, 113–126. Cayron, C. (2006). Acta Crystallogr A Found Crystallogr 62, 21–40. Day, A. P. (2008). Journal of Microscopy 230, ...

    1993

  2. [3]

    The composition law of quaternions and the spherical laws A rotation is the composition of two mirror planes (tail, head) such that these planes intersect along the rotation axis, and the angle between them is the semiangle of the rotation angle. The composition of two rotations can thus be geometrically constructed by decomposing them such that the head ...

    2006

  3. [590]

    & Field, D

    Morawiec, A. & Field, D. P. (1996). Philosophical Magazine A 73, 1113–1130. Nolze, G. (2015). Cryst. Res. Technol. 50, 188–201. Schwarzer, R. A. (1997). Micron 28, 249–265. Stewart, G. W. (1973). Introduction to matrix computations New York: Academic Press. Wilkinson, A. J. & Hirsch, P. B. (1997). Micron 28, 279–308. Winkelmann, A., Cios, G., Tokarski, T....

    1996

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.