The electric field gradient tensor as a symmetry-adapted order parameter in Landau theory

Technical Summary: The Electric Field Gradient Tensor as a Symmetry-Adapted Order Parameter in Landau Theory

Problem Statement
For decades, quadrupolar hyperfine spectroscopies (NQR, NMR, TDPAC, Mössbauer) have utilized the electric field gradient (EFG) tensor at a nuclear site as an empirical proxy for order parameters ($\phi$) in structural and electronic phase transitions. Historically, the EFG has been treated as a secondary observable: a "spectator" quantity coupled to a primary order parameter (such as spontaneous strain or polarization) defined elsewhere. Empirical relations, such as $\Delta\nu_Q \propto \langle\phi\rangle$ or $\langle\phi^2\rangle$, were established case-by-case, but the EFG itself lacked a systematic, symmetry-based role within Landau theory. Specifically, there was no general framework to determine a priori whether a specific EFG component constitutes an independent realization of the broken symmetry or merely a secondary response.

Methodology
The authors construct a unified framework combining site-symmetry representation theory, induction over Wyckoff orbits, and Landau free-energy analysis. The methodology proceeds through three logical steps:

1. Mathematical Definition: The EFG is defined as the traceless, symmetric, pure rank-2 ($\ell=2$) tensor of the electrostatic potential's second derivatives at a nuclear site. Crucially, the authors argue that tracelessness is an exact property derived from the decoupling of the trace from the quadrupole interaction, rather than an idealization relying on Laplace's equation in charge-free regions.
2. Representation Theoretical Construction:
   • The EFG tensor at a specific crystallographic site $q$ transforms under the site-symmetry group $G_q$.
   • To apply Landau theory, which requires irreducible representations (irreps) of the parent space group $G_0$, the authors induce the representation from the site group over the entire Wyckoff orbit.
   • This induction yields a decomposition of the EFG tensor into irreps of the parent group.
3. Landau Coupling Analysis: The authors analyze the coupling between the transition order parameter $\phi$ (transforming as irrep $\Gamma_\phi$) and the orbit-adapted EFG channels $v_\alpha$ (transforming as $\Gamma_\alpha$). By minimizing the Landau free energy, they derive the leading-order dependence of $v_\alpha$ on $\phi$.

Key Contributions
The paper establishes that the EFG tensor is not merely a probe but can be a symmetry-adapted realization of the order parameter itself. The central theoretical result is a classification of EFG channels into three distinct symmetry classes:

• Primary Channels (Linear): If the induced EFG representation contains the transition's irrep $\Gamma_\phi$, the corresponding orbit-adapted EFG combination vanishes identically above the critical temperature $T_c$ and grows linearly with $\phi$ below $T_c$ ($v \propto \phi$). These channels inherit the exact critical exponent, sign, and domain structure of the order parameter.
• Secondary Channels (Quadratic): EFG components belonging to other irreps (where $\Gamma_\alpha \neq \Gamma_\phi$ but $\Gamma_\alpha \subset [\Gamma_\phi^2]$) grow quadratically ($v \propto \phi^2$). This recovers the classic empirical $\Delta\nu_Q \propto \phi^2$ relations as a theorem derived from symmetry, rather than an ansatz.
• Forbidden Channels (Zero): EFG combinations that do not appear in the induced representation are symmetry-forbidden and must remain identically zero in equilibrium, regardless of temperature.

The authors emphasize that while the EFG realizes the order parameter (transforming identically to it), it does not drive the thermodynamic instability; the transition is driven by the primary microscopic degree of freedom.

Results and Verification
The framework is validated through two independent "halves": thermodynamic behavior via experimental literature and structural form via first-principles calculations.

1. First-Principles Verification ($\alpha$-Quartz):

   • The authors performed frozen-mode all-electron DFT calculations on $\alpha$-quartz, scanning the amplitude of the soft mode.
   • They decomposed the EFG tensor over the Si Wyckoff orbit into in-phase ($v_+$) and out-of-phase ($v_-$) combinations.
   • Result: The combination matching the distortion irrep ($v_-$) switched on linearly with the distortion amplitude ($R^2 = 0.9999$), while its symmetry-orthogonal partner ($v_+$) was suppressed by two orders of magnitude ($|a_+/a_-| \approx 0.018$). This confirms the orbit-selection prediction of the induced-representation theorem.

2. Experimental Consistency:

   • The framework reinterprets decades of NQR/NMR data. It explains why some transitions show linear critical exponents (primary channels) while others show quadratic behavior (secondary channels or "masqueraded" primary channels where powder averaging hides the linear dependence).
   • AgNa(NO$_2$)$_2$: Oriented-crystal satellite NMR data allows for a channel-resolved test, distinguishing linear and quadratic components.
   • Rb$_2$ZnBr$_4$: The framework explains the coexistence of different exponents ($\beta$ and $2\beta$) in different observables at the normal-to-incommensurate transition.
   • Null Results: The framework provides exclusionary logic; for example, the absence of an EFG anomaly in URu$_2$Si$_2$ at the hidden-order temperature generically excludes quadrupolar hidden order.

3. Proposed Case Study (BaFe$_2$As$_2$):

   • The authors propose a specific test for the nematic transition in BaFe$_2$As$_2$ using the $^{75}$As site.
   • Prediction: The transition irrep ($B_{2g}$) matches the $V_{xy}$ channel. The theory predicts that the in-phase combination of the two As sublattices ($v_+$) will grow linearly with the nematic order parameter, while the out-of-phase combination ($v_-$) must remain exactly zero. This "null prediction" is presented as a falsifiable signature.

Significance and Claims
The paper claims to supply the missing construction that elevates the EFG from a supporting empirical proxy to a principal order-parameter carrier in Landau theory.

• Theoretical Impact: It provides a rigorous symmetry criterion to determine a priori whether a quadrupolar observable is the order parameter or its shadow. It derives the linear/quadratic dichotomy as a consequence of group theory rather than phenomenological assumption.
• Practical Utility: The framework enables mutual calibration between first-principles calculations and Landau phenomenology. The coupling coefficients in the Landau free energy are exactly the derivatives computed by frozen-distortion DFT. Conversely, the symmetry-protected "zero" and parity rules in DFT provide internal consistency checks for hyperfine calculations.
• Modesty of Claims: The authors explicitly state that the EFG realizes the order parameter but does not drive the transition. The framework is restricted to zone-center ($k=0$) transitions of time-even order parameters. It constrains the mean values of the EFG (symmetry zeros) but does not constrain fluctuations (variances), which require dynamic extensions. The quantitative amplitudes remain subject to the approximations of electronic-structure theory, though the symmetry structure (parity and zeros) is exact and functional-independent.

In summary, the paper argues that the EFG tensor, when analyzed through the lens of induced representations over Wyckoff orbits, offers a robust, falsifiable, and symmetry-adapted realization of order parameters, bridging the gap between hyperfine spectroscopy and the fundamental thermodynamics of phase transitions.