Inversion-Symmetry Order Parameter

• Inversion-symmetry order parameter is a quantitative descriptor that measures the extent of inversion breaking, defined to vanish in centrosymmetric phases.
• It is computed via methods tailored to specific systems, such as mass dimerization in quasicrystals, macroscopic polarization in ferroelectrics, and symmetry indicators in crystalline band theory.
• This parameter underpins the realization of topologically protected states, vibrational anomalies, and multi-component phase transitions across various condensed matter platforms.

The inversion-symmetry order parameter is a fundamental quantitative descriptor for the extent and nature of inversion-symmetry breaking in condensed matter systems. Its formulation varies with physical context—spanning lattice models, vibrational networks, crystalline and topological phases, ferroelectrics, quasicrystals, and higher-order topological insulators—yet it is always constructed to be odd under the inversion operation and to vanish in genuinely centrosymmetric phases. Below, its mathematical definitions, physical manifestations, topological implications, and connections to major theoretical and experimental paradigms are detailed across canonical systems and methodologies.

1. Mathematical Formulations Across Physical Settings

Quasicrystalline Lattices: Mass Dimerization Parameter

In two-dimensional quasicrystal lattices, inversion symmetry is quantified via a mass-dimerization parameter $\delta$ associated with two inequivalent sublattices, A and B. Assigning site masses as $m_A = m_0(1+\delta)$, $m_B = m_0(1-\delta)$, the discrete Hamiltonian possesses an explicit inversion-breaking term proportional to $\delta$: $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$ Under inversion, A $\leftrightarrow$ B and hence $\delta \mapsto -\delta$, establishing $\delta$ as an odd, bona fide inversion-symmetry order parameter. The value $\delta=0$ signals perfect inversion symmetry; $\delta \ne 0$ quantifies its breaking (Beli et al., 2023).

Lattice Polarization: Macroscopic Displacement Order Parameter

For classical statistical models of polar displacive transitions, such as ferroelectric-like lattices, inversion breaking is measured by the macroscopic polarization $m_A = m_0(1+\delta)$0, evaluated as the magnitude of the site-averaged displacement vector: $m_A = m_0(1+\delta)$1 where $m_A = m_0(1+\delta)$2 are local discrete displacements (e.g., $m_A = m_0(1+\delta)$3 on a square or cubic lattice). $m_A = m_0(1+\delta)$4 enforces global inversion symmetry, while $m_A = m_0(1+\delta)$5 signals spontaneous symmetry breaking and macroscopic polarization (Wolpert et al., 2017).

Local Symmetry in Glasses and Defective Crystals

The local inversion-symmetry breaking (ISB) order parameter $m_A = m_0(1+\delta)$6, introduced to identify the origin of the boson peak and nonaffine mechanical response in disordered solids, is defined via unbalanced affine-force fields on each atom: $m_A = m_0(1+\delta)$7 where $m_A = m_0(1+\delta)$8 collects noncentrosymmetry-induced affine forces under strain, and the normalization ensures $m_A = m_0(1+\delta)$9, with $m_B = m_0(1-\delta)$0 for perfect centrosymmetry, $m_B = m_0(1-\delta)$1 for a maximally asymmetric network (Milkus et al., 2016).

Inversion Indicators in Band Theory

For centrosymmetric crystalline insulators (and superconductors), the inversion order parameter is encoded in symmetry indicators computed from inversion (parity) eigenvalues $m_B = m_0(1-\delta)$2 at time-reversal-invariant momenta (TRIMs): $m_B = m_0(1-\delta)$3 with $m_B = m_0(1-\delta)$4 the counts of occupied Bloch states with parity $m_B = m_0(1-\delta)$5 at $m_B = m_0(1-\delta)$6. This $m_B = m_0(1-\delta)$7 index fully classifies inversion-protected higher-order topological phases, trivial insulators, and Weyl semimetal constraints. In superconductors, the analogous symmetry indicator is constructed from the difference in BdG-band representations between physical and atomic-limit configurations, expressed as a $m_B = m_0(1-\delta)$8-valued invariant in $m_B = m_0(1-\delta)$9 dimensions (Kim et al., 2018, Skurativska et al., 2019).

Multi-Order-Parameter Structural Transitions

In systems such as Cd$\delta$0Re$\delta$1O$\delta$2, inversion breaking involves multiple coupled order parameters, e.g., $\delta$3 (primary odd), $\delta$4 (primary even), and $\delta$5 (secondary odd), each transforming as distinct irreducible representations under the parent symmetry. The full structural distortion is parametrized as: $\delta$6 accompanied by a Landau free energy with trilinear invariants, allowing the secondary order to be induced when both primaries condense (Norman, 2019).

2. Physical Significance and Inversion-Breaking Mechanisms

Across these systems, the inversion-symmetry order parameter quantifies a wide variety of physical regimes:

• In quasicrystals, finite $\delta$7 opens a full spectral gap and partitions the system into domains of distinct "valley-Hall" topological insulating character. Domain walls between regions of opposite $\delta$8 host topologically protected interface states, even in the absence of periodicity (Beli et al., 2023).
• In polar materials and models, $\delta$9 directly signals global inversion symmetry breaking, concomitant with the emergence of ferroelectricity (Wolpert et al., 2017).
• In amorphous and defective systems, $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$0 is directly linked to the presence of the boson peak in the vibrational density of states and to nonaffine softening of the elastic shear modulus, unifying the framework for disorder-induced anomalies (Milkus et al., 2016).
• In crystalline band insulators, the quantized inversion indicator $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$1 and its higher-order generalizations (e.g., $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$2 for $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$3 superconductors) partition phases into trivial, weak, and higher-order topological classes (Kim et al., 2018, Skurativska et al., 2019).
• In multi-component phase transitions, as in Cd$$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$4Re$$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$5O$$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$6, an interplay of coupled order parameters yields unconventional criticality and "Landau-violating" transitions, with measurable tensorial signatures in nonlinear optical experiments (Norman, 2019).

3. Computation and Experimental Probes

The explicit computation of the inversion-symmetry order parameter depends on the chosen context:

• In quasicrystal models, $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$7 is prescribed via onsite mass modulation.
• For lattice polarization, $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$8 is evaluated as a Monte Carlo (or experimental) thermal average of local displacements.
• $$H_\mathrm{mass} = \sum_i m_i \omega^2 |w_i|^2, \quad m_i = m_0 [1 + s_i \delta], \quad s_i = \begin{cases} +1 & \text{A-site} \ -1 & \text{B-site} \end{cases}$$9 is computed by analyzing the bond network and atomic positions, extracting nonzero $\leftrightarrow$0 values, and normalizing accordingly.
• Band-structure-based indicators are constructed from inversion eigenvalues at TRIMs obtained from first-principles calculations, often in combination with irreducible representations for multi-symmetry indicators (e.g., combining glide, inversion, and Chern numbers in space groups 13 and 14 (Kim et al., 2018)).
• In structural transitions, order parameters couple directly to measurable quantities such as second harmonic generation (SHG) tensors, lattice distortion amplitudes in scattering experiments, and phonon softening lineshapes.

4. Topological and Phase Classification Implications

The inversion-symmetry order parameter serves as a topologically robust classifier:

• In quasicrystals, its sign determines the topology of the insulating phase, ensuring that spatial variation (domain wall formation) leads to interface-localized (Jackiw–Rebbi) zero modes with decay length $\leftrightarrow$1 (Beli et al., 2023).
• For higher-order topological phases, bulk parity indicators built from inversion eigenvalues at high-symmetry points not only detect insulating versus semimetallic regimes but also the presence and type of edge, hinge, and corner states (Xiong et al., 2020, Skurativska et al., 2019).
• In superconductors, the $\leftrightarrow$2 inversion symmetry indicator partitions first-, second-, and higher-order topological phases, directly predicting boundary Majorana modes, without requiring the explicit BdG spectrum (Skurativska et al., 2019).
• In crystalline insulators with glide, inversion, and time-reversal symmetries, the even-valued $\leftrightarrow$3 identifies higher-order topological insulators with inversion-protected chiral hinge modes, while odd values forbid fully gapped phases (enforcing Weyl nodes) (Kim et al., 2018).

5. Comparative Summary of Key Definitions and Roles

System	Order Parameter	Symmetry Action (P)	Physical Interpretation
Quasicrystals (Beli et al., 2023)	$\leftrightarrow$4	$\leftrightarrow$5	Mass dimerization; controls band inversion; domain-wall modes
Ferroelectric Models (Wolpert et al., 2017)	$\leftrightarrow$6	$\leftrightarrow$7 (under $\leftrightarrow$8)	Average polarization; signals global inversion symmetry breaking
Glasses/Defective Crystals (Milkus et al., 2016)	$\leftrightarrow$9	$\delta \mapsto -\delta$0 (centrosymmetric), $\delta \mapsto -\delta$1 (ISB)	Correlates with boson peak, nonaffine modulus, local symmetry
Band Insulators (Kim et al., 2018)	$\delta \mapsto -\delta$2	Inverts under move across parity inversion	Determines topological class (trivial, HOTI, semimetal)
Superconductors (Skurativska et al., 2019)	$\delta \mapsto -\delta$3, $\delta \mapsto -\delta$4 indicator	Determined by pairing parity $\delta \mapsto -\delta$5	Bulk-boundary Majorana correspondence; higher-order phases
Structural transitions (Norman, 2019)	$\delta \mapsto -\delta$6, $\delta \mapsto -\delta$7, $\delta \mapsto -\delta$8	Assigned by irrep (odd/even under $\delta \mapsto -\delta$9)	Multi-order-parameter criticality; SHG signatures

6. Generalized Lessons and Connections

Several unifying principles emerge:

• The inversion-symmetry order parameter is the coefficient of the only inversion-odd term in the system's Hamiltonian or free energy, and its domain-wall variation is a rigorous guarantee of protected in-gap states or solitonic excitations.
• In periodic crystals, band parity at TRIMs operationalizes the inversion order parameter and bridges symmetry representation theory with bulk, edge, and higher-order topology.
• In aperiodic, disordered, or amorphous systems, local force or displacement-based order parameters provide a robust continuum between crystalline and glassy phases, connecting vibrational anomalies to symmetry breaking directly, even when other order parameters (e.g., bond-orientational) fail.

The inversion-symmetry order parameter thus constitutes a central, multifaceted tool in diagnosing, quantifying, and exploiting inversion-symmetry breaking—whether the goal is to design topologically nontrivial phases, track ferroelectricity, explain vibrational anomalies in disorder, or probe unconventional order in complex crystalline lattices (Beli et al., 2023, Wolpert et al., 2017, Kim et al., 2018, Milkus et al., 2016, Norman, 2019, Skurativska et al., 2019, Xiong et al., 2020).