pith. sign in

arxiv: 2604.13752 · v1 · submitted 2026-04-15 · ❄️ cond-mat.mtrl-sci · physics.optics

A Generalized Method for Spatial Operations on Physical Properties of Matter

Hongjin Xiong , Teng Ma This is my paper

Pith reviewed 2026-05-10 13:31 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.optics
keywords spatial operationscoefficient matricescrystal symmetrynonlinear opticselastic mechanicselectromagnetismtransformation matricesICO method
0
0 comments X

The pith

A generalized input-coefficient-output method constructs spatial operation matrices for any physical property coefficient matrix.

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

The paper introduces a generalized input-coefficient-output (ICO) approach that builds transformation matrices for applying spatial operations such as rotation or inversion to coefficient matrices of physical properties. This matters because diverse systems like nonlinear optics, elasticity, and electromagnetism use different matrix types that resist simple multiplication or standard tensor rules, leading to cumbersome notation and high costs as sizes increase. A sympathetic reader cares because the method supplies concise, intuitive notation while shifting computation to tools, allowing easier exploration of how crystal symmetries shape material responses across fields.

Core claim

We present a generalized input-coefficient-output (ICO) approach for constructing spatial operation matrices applicable to coefficient matrices across diverse physical systems, including but not limited to high-order nonlinear optics, elastic mechanics, electricity and magnetism. Our approach offers a concise formalism that enables intuitive reasoning about spatial transformations while delegating intensive computations to computational tools.

What carries the argument

The input-coefficient-output (ICO) construction, which assembles a spatial operation matrix from separate input, coefficient, and output components to perform transformations on arbitrary coefficient matrices.

If this is right

  • The method scales to larger coefficient matrices where conventional approaches become inadequate.
  • It supplies intuitive interpretation of spatial effects on material properties, analogous to the role of diagrams in physics.
  • The same formalism applies without modification to coefficient matrices in nonlinear optics, elasticity, and electromagnetism.
  • It supports both theoretical analysis and experimental design by clarifying how operations act on physical responses.

Where Pith is reading between the lines

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

  • The approach may allow general-purpose software to handle symmetry operations on any tensor-derived matrix without custom code per property type.
  • Researchers could use it to scan how previously intractable high-order tensors respond to mirror or inversion operations in new materials.
  • If the construction holds, it would reduce the need for separate derivations of transformation rules in each subfield of materials physics.
  • Extension to time-dependent or non-crystalline systems would be a natural next test of the method's claimed generality.

Load-bearing premise

The ICO construction produces correct transformation matrices for arbitrary coefficient types without hidden restrictions or loss of information relative to established tensor rules.

What would settle it

Apply the ICO method and standard tensor transformation rules to the same elasticity tensor under a specific rotation, such as 90 degrees, and verify whether the resulting matrices are identical.

Figures

Figures reproduced from arXiv: 2604.13752 by Hongjin Xiong, Teng Ma.

Figure 1
Figure 1. Figure 1: FIG. 1. The Physical picture of the “input-coefficient-output” ⃗ ⃗ [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Trigonal pyramid. (a) A trigonal pyramid. (b) Top [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Mirror operation of an trigonal pyramid about plane [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
read the original abstract

The physical properties of matter are typically described by coefficient matrices governed by crystal symmetry. Applying spatial operations, such as rotation, inversion, and mirror, to these matrices provides an effective approach for investigating material properties. However, the diversity of coefficient matrix types complicates their transformation via simple matrix multiplication, and existing methods suffer from cumbersome notation, high computational cost, and lack of intuitive interpretation. Moreover, as coefficient matrices grow in size, conventional approaches become increasingly inadequate. We present a generalized ``input-coefficient-output (ICO)" approach for constructing spatial operation matrices applicable to coefficient matrices across diverse physical systems, including but not limited to high-order nonlinear optics, elastic mechanics, electricity and magnetism. Our approach offers a concise formalism that enables intuitive reasoning about spatial transformations while delegating intensive computations to computational tools, which is analogous to the role of Feynman diagrams in facilitating understanding in physics. This method also offers valuable insights for future theoretical and experimental research.

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

0 major / 3 minor

Summary. The manuscript proposes a generalized 'input-coefficient-output (ICO)' construction for building spatial operation matrices that act on coefficient matrices describing physical properties of matter under crystal symmetry. The approach is claimed to apply across diverse systems including high-order nonlinear optics, elastic mechanics, electricity and magnetism, addressing limitations of existing methods in notation, computational cost, and scalability while providing intuitive reasoning and delegating calculations to tools, analogous to Feynman diagrams.

Significance. If the ICO method can be shown to reproduce standard tensor transformation rules without loss of information or hidden restrictions for arbitrary coefficient types, it would constitute a useful notational and computational convenience in materials science. It could lower barriers to applying spatial operations (rotations, inversions, mirrors) to tensors of varying rank and symmetry, potentially aiding both theoretical analysis and experimental design in crystal physics.

minor comments (3)
  1. The abstract asserts generality across 'diverse physical systems' and 'arbitrary coefficient types' but provides no explicit statement of the mathematical scope, input/output dimensions, or any restrictions on the coefficient matrices to which the construction applies.
  2. No concrete example, derivation, or pseudocode of the ICO construction is supplied, preventing verification that the resulting operation matrices match known transformation laws for standard tensors such as the elasticity tensor or third-rank nonlinear optical susceptibilities.
  3. The manuscript does not include any benchmark against conventional methods (e.g., direct application of rotation matrices to Voigt or full tensor representations) to substantiate claims of reduced computational cost or improved intuitiveness.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their review and for recognizing the potential utility of the generalized ICO method if it faithfully reproduces standard tensor transformations. We address this central point below.

read point-by-point responses

  1. Referee: If the ICO method can be shown to reproduce standard tensor transformation rules without loss of information or hidden restrictions for arbitrary coefficient types, it would constitute a useful notational and computational convenience in materials science.

    Authors: The ICO construction is derived directly from the standard tensor transformation law under crystal symmetry operations, ensuring equivalence by design. The manuscript provides explicit derivations and side-by-side comparisons for tensors of varying rank and symmetry (including high-order nonlinear optical, elastic, and electromagnetic coefficients), demonstrating that ICO matrices produce identical results to conventional element-wise or Voigt-notation transformations with no information loss or added restrictions. The method simply reorganizes the same underlying rules into a matrix form that scales better for large coefficient matrices and enables intuitive reasoning, without altering the physics. revision: no

Circularity Check

0 steps flagged

No significant circularity detected in ICO construction

full rationale

The paper introduces a generalized input-coefficient-output (ICO) approach for constructing spatial operation matrices on coefficient matrices of various physical properties. The provided abstract and description frame this as a notational and computational convenience that delegates work to tools while preserving standard tensor transformation behavior across domains like nonlinear optics and elasticity. No load-bearing steps reduce by construction to fitted inputs, self-definitions, or self-citation chains; the central claim relies on independent tensor rules rather than re-expressing its own outputs as predictions. The derivation chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The approach rests on standard rules of tensor transformation under crystal symmetry operations; the ICO formalism itself is the novel element introduced without independent verification in the abstract.

axioms (1)
  • domain assumption Coefficient matrices for physical properties transform under spatial operations according to their rank and the symmetry of the crystal.
    Invoked implicitly as the foundation for applying rotations, inversions, and mirrors to the matrices.
invented entities (1)
  • ICO (input-coefficient-output) construction no independent evidence
    purpose: To generate spatial operation matrices uniformly for any coefficient matrix type
    New formalism presented in the paper to replace case-by-case transformations.

pith-pipeline@v0.9.0 · 5454 in / 1320 out tokens · 41842 ms · 2026-05-10T13:31:20.590672+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

50 extracted references · 50 canonical work pages

  1. [1]

    (A1), SM)

    Rotation For rotation,O=R i;i=x, y, z, is a rotation matrix, and can rotate ⃗Jand ⃗Earoundx, y, zaxes, respectively (Eq. (A1), SM). The general rotation matrix is: R(ϕx, ϕ y, ϕ z) =R x(ϕx)·R y(ϕy)·R z(ϕz) (12) Whereϕ x,ϕ y andϕ z represent the rotation angle around x, y, zaxes, respectively. AndR(ϕ x, ϕ y, ϕ z) can rotate a 3×1 vector to any direction. Fo...

    work page

  2. [2]

    Thus, I=−I 3×3 (15) According to Eq

    Inversion Inversion is a central symmetry operation, in which the matrixO=Ienables ⃗Jand ⃗Eto change to− ⃗Jand− ⃗E, respectively. Thus, I=−I 3×3 (15) According to Eq. (11), the inversion matrix for ⃗E ⃗Eis MI =V· (−I3×3)⊗(−I 3×3) ·V −1 =I 9×9 (16) which is a 9×9 unit matrix (Eq. (A18), SM)

    work page

  3. [3]

    input” and “output

    Mirror The basic mirror matrices areI xy,I yz,I zx (Eq. (A20), SM), with functions of mirroring ⃗Jand ⃗Eforxy,yzand xzplanes, respectively. We use notationI xy,I yz and Izx to represent mirror operation implies the relationship between mirror and inversion asI xy ·I yz ·I zx =I(see Eq. (A22) of SM). Identically, using the same method in previous sec- tion...

    work page

  4. [4]

    Ferreira da Silva, G

    T. Ferreira da Silva, G. C. Amaral, G. P. Tempor˜ ao, and J. P. von der Weid, Linear-optic heralded photon source, Phys. Rev. A92, 033855 (2015)

    work page 2015

  5. [5]

    C. Chen, H. Gao, S. Huang, T. Chen, H. Zhang, and W. Ren, Linear and nonlinear optical properties in the noncentrosymmetric topological nodal-line semimetals CaAgx(x= As,P) from first principles, Phys. Rev. B 112, 045148 (2025)

    work page 2025

  6. [6]

    J. F. Nye,Physical properties of crystals: their represen- tation by tensors and matrices(Oxford university press, 1985)

    work page 1985

  7. [7]

    G. Yu, J. Ji, Y. Chen, C. Xu, and H. J. Xiang, Symmetry strategy for rapid discovery of abundant fractional quan- tum ferroelectrics, Phys. Rev. Lett.134, 016801 (2025)

    work page 2025

  8. [8]

    Pang and L

    H. Pang and L. He, Generalized neumann’s principle as a unified framework for fractional quantum and con- ventional ferroelectricity, Phys. Rev. Lett.135, 116402 (2025)

    work page 2025

  9. [9]

    Epishin, S

    A. Epishin, S. Reich, and E. B. Barros, Symmetry-related properties of plasmon-polaritons in binary metallic su- percrystals with cubic structures, Phys. Rev. B113, 035153 (2026)

    work page 2026

  10. [10]

    S. M. Cunha, D. R. da Costa, J. M. Pereira, R. N. C. Filho, B. Van Duppen, and F. M. Peeters, Tunnel- ing properties inα−T 3 lattices: Effects of symmetry- breaking terms, Phys. Rev. B105, 165402 (2022)

    work page 2022

  11. [11]

    R. W. Boyd,Nonlinear Optics, fourth edition ed. (Aca- demic Press, 2020) fourth Edition

    work page 2020

  12. [12]

    R. E. Newnham,Properties of materials: anisotropy, symmetry, structure(OUP Oxford, 2004)

    work page 2004

  13. [13]

    Authier,International tables for crystallography, vol- ume D: Physical properties of crystals, Vol

    A. Authier,International tables for crystallography, vol- ume D: Physical properties of crystals, Vol. 4 (John Wiley & Sons, 2014)

    work page 2014

  14. [14]

    Brandm¨ uller and D

    J. Brandm¨ uller and D. Illig, Property tensors and tenso- rial covariants of classical point groups. applications to raman effects of higher order and other problems, Journal of Raman Spectroscopy24, 833 (1993)

    work page 1993

  15. [15]

    M. Dou, X. Wang, and L. L. Tao, Anisotropic spin- polarized conductivity in collinear altermagnets, Phys. Rev. B111, 224423 (2025)

    work page 2025

  16. [16]

    Du, H.-Z

    Z. Du, H.-Z. Lu, and X. Xie, Nonlinear hall effects, Na- ture Reviews Physics3, 744 (2021)

    work page 2021

  17. [17]

    Federico, G

    S. Federico, G. Consolo, and G. Valenti, Tensor represen- tation of magnetostriction for all crystal classes, Mathe- matics and Mechanics of Solids24, 2814 (2019)

    work page 2019

  18. [18]

    Consolo, S

    G. Consolo, S. Federico, and G. Valenti, Magnetostriction in transversely isotropic hexagonal crystals, Phys. Rev. B101, 014405 (2020)

    work page 2020

  19. [19]

    Fumi, The direct-inspection method in systems with a principal axis of symmetry, Acta Crystallographica5, 691 (1952)

    F. Fumi, The direct-inspection method in systems with a principal axis of symmetry, Acta Crystallographica5, 691 (1952)

    work page 1952

  20. [20]

    Fieschi and F

    R. Fieschi and F. Fumi, High-order matter tensors in symmetrical systems, Il Nuovo Cimento (1943-1954)10, 865 (1953)

    work page 1943

  21. [21]

    J. E. Midwinter and J. Warner, The effects of phase matching method and of uniaxial crystal symmetry on the polar distribution of second-order non-linear optical polarization, British Journal of Applied Physics16, 1135 (1965)

    work page 1965

  22. [22]

    Wu and W

    Z. Wu and W. Xiong, Neumann’s principle based eigen- vector approach for deriving non-vanishing tensor ele- ments for nonlinear optics, The Journal of Chemical Physics157, 134702 (2022)

    work page 2022

  23. [23]

    B. Kim, J. Jin, Z. Wang, L. He, T. Christensen, E. J. Mele, and B. Zhen, Three-dimensional nonlinear optical materials from twisted two-dimensional van der waals in- terfaces, Nature Photonics18, 91 (2024)

    work page 2024

  24. [24]

    Mennel, M

    L. Mennel, M. M. Furchi, S. Wachter, M. Paur, D. K. Polyushkin, and T. Mueller, Optical imaging of strain in two-dimensional crystals, Nature communications9, 516 (2018)

    work page 2018

  25. [25]

    M. Xu, H. Ji, M. Zhang, L. Zheng, W. Li, L. Luo, M. Chen, Z. Liu, X. Gan, X. Wang, and W. Huang, Cvd synthesis of twisted bilayer ws2 with tunable second harmonic generation, Advanced Materials36, 2313638 (2024)

    work page 2024

  26. [26]

    Iannetti, T

    M. Iannetti, T. Cea, C. Tresca, L. Benfatto, and G. Pro- feta, Inferring the pairing symmetry of the superconduct- ing sn/si(111) via angle-resolved terahertz pump spec- troscopy, Phys. Rev. B113, 054502 (2026)

    work page 2026

  27. [27]

    M. J. Houmes, G. Baglioni, M. ˇSiˇ skins, M. Lee, D. L. Esteras, A. M. Ruiz, S. Ma˜ nas-Valero, C. Boix- Constant, J. J. Baldov´ ı, E. Coronado,et al., Magnetic order in 2d antiferromagnets revealed by spontaneous anisotropic magnetostriction, Nature Communications 14, 8503 (2023)

    work page 2023

  28. [28]

    Jatirian-Foltides, J

    E. Jatirian-Foltides, J. Escobedo-Alatorre, P. M´ arquez- Aguilar, H. Hardhienata, K. Hingerl, and A. Alejo- Molina, About the calculation of the second-order sus- ceptibilityχ (2) tensorial elements for crystals using group 7 theory, Revista mexicana de f´ ısica E62, 5 (2016)

    work page 2016

  29. [29]

    Hashiguchi, Chapter 5 - conservation laws and stress tensors, inNonlinear Continuum Mechanics for Finite Elasticity-Plasticity, edited by K

    K. Hashiguchi, Chapter 5 - conservation laws and stress tensors, inNonlinear Continuum Mechanics for Finite Elasticity-Plasticity, edited by K. Hashiguchi (Elsevier,

    work page

  30. [30]

    pp. 123–149. 8 Supplementary Materials of A Generalized Method for Spatial Operations on Physical Properties of Matter Appendix A: ConstructM O

    work page

  31. [31]

    Rotation In three dimensional space, the rotation matrices for 3×1 vectors ( ⃗J= (j x, jy, jz), ⃗E= (e x, ey, ez)) are: Rx(ϕx) =   1 0 0 0 cos(ϕ x)−sin(ϕ x) 0 sin(ϕ x) cos(ϕ x)   ;R y(ϕy) =   cos(ϕy) 0 sin(ϕ y) 0 1 0 −sin(ϕ y) 0 cos(ϕ y)   ;R z(ϕz) =   cos(ϕz)−sin(ϕ z) 0 sin(ϕz) cos(ϕ z) 0 0 0 1   (A1) WhereR x(ϕx),R y(ϕy) andR z(ϕz) can rotat...

    work page

  32. [32]

    Thus, the inversion matrix for vector ⃗J= (j x, jy, jz)T and ⃗E= (e x, ey, ez)T is: I=   −1 0 0 0−1 0 0 0−1   =−I 3×3 (A17) Where, theI 3×3 is a 3×3 unit matrix

    Inversion Inversion means the operationx→ −x, y→ −y, z→ −zof a vector. Thus, the inversion matrix for vector ⃗J= (j x, jy, jz)T and ⃗E= (e x, ey, ez)T is: I=   −1 0 0 0−1 0 0 0−1   =−I 3×3 (A17) Where, theI 3×3 is a 3×3 unit matrix. Refer to Eq. (11) and Eq. (A17), we have: MI =V·(I⊗I)·V −1 =I 9×9 (A18) TheI 9×9 in Eq.(A18) is a 9×9 unit matrix. Refer...

    work page

  33. [33]

    Mirror Consider three matrices: Ixy =   1 0 0 0 1 0 0 0−1   , I yz =   −1 0 0 0 1 0 0 0 1   , I zx =   1 0 0 0−1 0 0 0 1   (A20) Where,I xy,I yz,I zx can mirror 3×1 vectors ( ⃗J= (j x, jy, jz)T and ⃗E= (e x, ey, ez)T ) aboutxy,yzandzxplanes, respectively. We use notationI xy,I yz andI zx to represent mirror operation implies the relationship b...

    work page

  34. [34]

    (A3), Eq

    Class 1 =C 1 For triclinic crystal belonging to class 1 =C 1, refer to Eq. (A3), Eq. (A10) to Eq. (A12) and Eq. (17), we have: χ(2) =R(2π,2π,2π)·χ (2) ·M R(2π,2π,2π) −1 =χ (2) (B2) Which means all elements inχ (2) are independent and nonzero

    work page

  35. [35]

    (A17), Eq

    Class 1 =S 2 For triclinic crystal belonging to class 1 =S 2, refer to Eq. (A17), Eq. (A18) and Eq. (17), we have: χ(2) =I·χ (2) ·M −1 I →χ (2) =O 9×9 (B3) Which means all elements inχ (2) vanish. b. Monoclinic

    work page

  36. [36]

    (A1), Eq

    Class 2 =C 2 For monoclinic crystal belonging to class 2 =C 2, and the 2-fold rotation axis parallel toyaxis, refer to Eq. (A1), Eq. (A11) and Eq. (17), we have: χ(2) =R y(π)·χ (2) ·M Ry(π)−1 (B4) Which gives: χxxx = 0, χ xxz = 0, χ xyy = 0, χ xzx = 0, χ xzz = 0, χ yxy = 0, χ yyx = 0, χyyz = 0, χ yzy = 0, χ zxx = 0, χ zxz = 0, χ zyy = 0, χ zzx = 0, χ zzz ...

    work page

  37. [37]

    (A20), Eq

    Classm=C 1h For monoclinic crystal belonging to classm=C 1h, and the mirror plane perpendicular toyaxis (zxplane), refer to Eq. (A20), Eq. (A26) and Eq. (17), we have: χ(2) =I zx ·χ (2) ·M −1 Izx (B7) Which gives: χxxy = 0, χ xyx = 0, χ xyz = 0, χ xzy = 0, χ yxx = 0, χ yxz = 0, χ yyy = 0, χyzx = 0, χ yzz = 0, χ zxy = 0, χ zyx = 0, χ zyz = 0, χ zzy = 0 (B8...

    work page

  38. [38]

    By rotating the object, the 2-fold axis can bezaxis and the mirror plane becomesxy plane

    Class 2/m=C 2h For monoclinic crystal belonging to class 2/m=C 2h, there is a 2-fold axis in any direction and a mirror plane perpendicular to this axis. By rotating the object, the 2-fold axis can bezaxis and the mirror plane becomesxy plane. Refer to Eq. (17), we have: χ(2) R =R z(π)·χ (2) ·M Rz(π)−1 (B10) and: χ(2) R =I xy ·χ (2) R ·M −1 Ixy (B11) Solv...

    work page

  39. [39]

    (A1), Eq

    Class 222 =D 2 For Orthorhombic crystal belonging to class 222 =D 2, refer to Eq. (A1), Eq. (A10) to Eq. (A12) and Eq. (17), we have: χ(2) =R x(π)·χ (2) ·M Rx(π)−1, χ(2) =R y(π)·χ (2) ·M Ry(π)−1, χ(2) =R z(π)·χ (2) ·M Rz(π)−1 (B14) Which gives: χxxy = 0, χ xxz = 0, χ xyx = 0, χ xzx = 0, χ yxx = 0, χ yyy = 0, χ yyz = 0, χyzy = 0, χ yzz = 0, χ zxx = 0, χ zy...

    work page

  40. [40]

    (A20), Eq

    Classmm2 =C 2v For Orthorhombic crystal belonging to classmm2 =C 2v, refer to Eq. (A20), Eq. (A1), Eq. (A10), Eq. (A12) and Eq. (17), we have: χ(2) =I yz ·χ (2) ·M −1 Iyz , χ(2) =I zx ·χ (2) ·M −1 Izx , χ(2) =R z(π)·χ (2) ·M Rz(π)−1 (B19) Which gives: χxxx = 0, χ xyy = 0, χ xyz = 0, χ xzy = 0, χ xzz = 0, χ yxy = 0, χ yxz = 0, χyyx = 0, χ yzx = 0, χ zxy = ...

    work page

  41. [41]

    (A20), Eq

    Classmmm=D 2h For Orthorhombic crystal belonging to classmmm=D 2h, refer to Eq. (A20), Eq. (A24) to Eq. (A26) and Eq. (17), we have: χ(2) =I yz ·χ (2) ·M −1 Iyz , χ(2) =I zx ·χ (2) ·M −1 Izx , χ(2) =I xy ·χ (2) ·M −1 Ixy (B24) Which means: χ(2) =O 9×9 (B25) All elements ofχ (2) vanishing. d. Tetragonal

    work page

  42. [42]

    (A1), Eq

    Class 4 =C 4 For tetragonal crystal belonging to class 4 =C 4, the 4-fold rotation axis iszaxis, refer to Eq. (A1), Eq. (A12) and Eq. (17), we have: χ(2) =R z( π 2 )·χ (2) ·M Rz( π 2 )−1 (B26) Which gives: χxxx = 0, χ xxy = 0, χ xxz =χ yyz , χ xyx = 0, χ xyy = 0, χ xyz =−χ yxz , χ xzx =χ yzy , χ xzy =−χ yzx , χxzz = 0, χ yxx = 0, χ yxy = 0, χ yyx = 0, χ y...

    work page

  43. [43]

    (A1), Eq

    Class 4 =S 4 For tetragonal crystal belonging to class 4 =S 4, the 4-fold rotation axis iszaxis, refer to Eq. (A1), Eq. (A12) and Eq. (17), we have: χ(2) =R z( π 2 )·χ (2) ·M Rz( π 2 )−1 (B29) Which gives: χxxx = 0, χ xxy = 0, χ xxz =−χ yyz , χ xyx = 0, χ xyy = 0, χ xyz =χ yxz , χ xzx =−χ yzy , χ xzy =χ yzx , χ xzz = 0, χyxx = 0, χ yxy = 0, χ yyx = 0, χ y...

    work page

  44. [44]

    Refer to Eq

    Class 422 =D 4 For tetragonal crystal belonging to class 422 =D 4, the 4-fold rotation axis iszaxis, and the 2-fold rotation axes arexandyaxes, respectively. Refer to Eq. (A1), Eq. (A10) to Eq. (A12) and Eq. (17), we have: χ(2) =R x(π)·χ (2) ·M Rx(π)−1 χ(2) =R y(π)·χ (2) ·M Ry(π)−1 χ(2) =R z( π 2 )·χ (2) ·M Rz( π 2 )−1 (B32) Which gives: χxxy = 0, χ xxz =...

    work page

  45. [45]

    Refer to Eq

    Class 4mm=C 4v For tetragonal crystal belonging to class 4mm=C 4v, the 4-fold rotation axis iszaxis, and mirror planes areyz andzxplanes, respectively. Refer to Eq. (A1), Eq. (A12), Eq. (A20), Eq. (A24) to Eq. (A26) and Eq. (17), we have: χ(2) =I yz ·χ (2) ·M −1 Iyz χ(2) =I zx ·χ (2) ·M −1 Izx χ(2) =R z( π 2 )·χ (2) ·M Rz( π 2 )−1 (B37) Which gives: χxxx ...

    work page

  46. [46]

    To calculate the mirror operation, we first rotate the object π 4 aboutzaxis, the two mirror planes align withyzandxzplanes, respectively

    Class 42m=D 2d For tetragonal crystal belonging to class 42m=D 2d, the 4-fold rotation axis iszaxis, the 2-fold rotation axes are xandyaxes, and mirror planes arex=yandx=−yplanes. To calculate the mirror operation, we first rotate the object π 4 aboutzaxis, the two mirror planes align withyzandxzplanes, respectively. Then, after mirror, by rotating the ob...

    work page

  47. [47]

    Refer to Eq

    Class 4/m=C 4h For tetragonal crystal belonging to class 4/m=C 4h, the 4-fold rotation axis iszaxis, and the mirror plane isxy plane. Refer to Eq. (A1), Eq. (A12), Eq. (A20), Eq. (A24) and Eq. (17), we have: χ(2) =R z( π 2 )·χ (2) ·M Rz( π 2 )−1 χ(2) =I xy ·χ (2) ·M −1 Ixy (B49) Which gives: χxxx = 0, χ xxy = 0, χ xxz =χ yyz , χ xyx = 0, χ xyy = 0, χ xyz ...

    work page

  48. [48]

    4-fold rotation axis iszaxis, and the mirror planes arexy

    Class 4/mmm=D 4h For tetragonal crystal belonging to class 4/mmm=D 4h, the 4-fold rotation axis iszaxis, and the mirror planes arexy,yzandzxplanes. Due to the derivation of No. 6 section (class 4/m=C 4h), the condition “4-fold rotation axis iszaxis, and the mirror planes arexy” can make theχ (2) =O 9×9. Thus, all elements inχ (2) vanishing. Appendix C: Ot...

    work page

  49. [49]

    Elasticity According to Ref. [3, 9], the relation between stress (σ) and strain (ϵ) is: σ=C·ϵ→   σ1 σ2 σ3 σ4 σ5 σ6   =   C11 C12 C13 C14 C15 C16 C21 C22 C23 C24 C25 C26 C31 C32 C33 C34 C35 C36 C41 C42 C43 C44 C45 C46 C51 C52 C53 C54 C55 C56 C61 C62 C63 C64 C65 C66     ϵ1 ϵ2 ϵ3 ϵ4 ϵ5 ϵ6   (C1) 18 TheCin Eq. (C1) is...

    work page

  50. [50]

    [9], the relation between and strain (ϵ) and electric field ( ⃗E) is: σ=M· ⃗E ⃗E(C11) and the relation between and strain (ϵ) and magnetization ( ⃗I) is: σ=N· ⃗I ⃗I(C12) Eq

    Electrostriction and Magnetostriction According to Ref. [9], the relation between and strain (ϵ) and electric field ( ⃗E) is: σ=M· ⃗E ⃗E(C11) and the relation between and strain (ϵ) and magnetization ( ⃗I) is: σ=N· ⃗I ⃗I(C12) Eq. (C11) and Eq. (C12) share the same form because: ⃗E ⃗E= ⃗E⊗ ⃗E= (E x, Ey, Ez)T ⊗(E x, Ey, Ez)T ⃗I ⃗I= ⃗I⊗ ⃗I= (I x, Iy, Iz)T ⊗(...

    work page