Nonlinear Flexoelectric Phonon Coupling in Ferroics

• Nonlinear flexoelectric phonon coupling is the trilinear interaction between strain gradients, order parameter fields, and macroscopic polarization in ferroic materials.
• It comprises symmetric flexoelectric-like and antisymmetric Dzyaloshinskii–Moriya-like components that alter phonon dispersion and trigger spatially modulated phases.
• Experimental and computational methods such as DFPT and LGD theory quantify the coupling tensors, predicting critical instabilities in complex oxides and van der Waals ferroics.

Nonlinear flexoelectric phonon coupling describes the interaction between strain gradients or structural order parameter fields and macroscopic polarization in ferroic materials, mediated by phonon modes and governed by nonlinear energy invariants. This coupling has both symmetric flexoelectric-like and antisymmetric Dzyaloshinskii–Moriya-like components and manifests in macroscopic polarization, phonon dispersion, instability criteria, and spectral response, with significant implications for the phase behavior and transport properties of complex oxides and van der Waals ferroics (Stengel, 2023, Morozovska et al., 23 Jul 2025, Morozovska et al., 8 Mar 2025, Morozovska et al., 2017, Morozovska et al., 2016, Morozovska et al., 2015, Morozovska et al., 2010).

1. Mathematical Formulation and Tensor Structure

Nonlinear flexoelectric phonon coupling is encapsulated by a four-rank tensor $W_{ij\alpha\beta}$ relating local order-parameter modulations (e.g., tilt, polarization, antiferrodistortive modes) and their gradients to the induced macroscopic polarization: $$P_i = W_{ij\alpha\beta}\;i\,q_j\,\phi_\alpha\,\phi_\beta$$ Here, $\phi_\alpha$ is a generalized order parameter, $q_j$ is the phonon wavevector, and $W_{ij\alpha\beta}$ admits a decomposition: $$W_{ij\alpha\beta} = K_{ij\alpha\beta} + \zeta (\delta_{i\alpha} \delta_{j\beta} - \delta_{i\beta} \delta_{j\alpha})$$

• The symmetric part $K_{ij\alpha\beta}$ represents flexoelectric-like coupling that is sensitive to boundary and gauge conditions.
• The antisymmetric component $\zeta$ is Dzyaloshinskii–Moriya-like, boundary-condition independent, and exists as a genuine bulk invariant.

The total coupling produces physical terms analogous to Lifshitz invariants in continuum free energies (Stengel, 2023): $$E_{\text{flexo}} \sim \frac{1}{2} K_{\alpha\beta\gamma\lambda}\left[\partial_\beta P_\alpha \phi_\gamma \phi_\lambda - P_\alpha \partial_\beta (\phi_\gamma \phi_\lambda)\right]$$

$$E_{\text{DM}} \sim \zeta\,P\cdot[\phi(\nabla\cdot\phi) - (\phi\cdot\nabla)\phi]$$

These invariants mediate trilinear coupling between polarization, order parameter fields, and their gradients.

2. First-Principles and Analytical Calculation Methods

The coupling tensors $W$ and $K$ are accessible via density-functional perturbation theory (DFPT), based on the expansion of the dynamical matrix with respect to long-wavelength phonon amplitudes and gradients (Stengel, 2023): $$\Phi(q) \simeq \Phi^{(0)} - i q_j \Phi^{(1,j)} + \mathcal{O}(q^2)$$ Uniform distortions $\phi^0_\beta$ induce a change in the first-order term, and projections onto polar modes yield $W_{ij\alpha\beta}$. All symmetry-allowed components are resolved automatically for cubic (or material-specific) symmetry, as illustrated concretely for SrTiO$_3$ where $K_{11}, K_{12}, K_{44}$ and $\zeta$ are calculated from LDA/ABINIT. Robustness against boundary conditions distinguishes $\zeta$ as a true bulk parameter, while $K$ requires flexoelastic renormalization for meaningful macroscopic response (Stengel, 2023).

In continuum analytical treatments (LGD theory), the nonlinear flexoelectric term appears in the free-energy functionals as $f_{ijkl} u_{ij} \partial_k P$ and its dynamic analog $M \dot U \dot P$ (Morozovska et al., 2017, Morozovska et al., 2016, Morozovska et al., 23 Jul 2025). These directly contribute to the equations of motion, dynamical matrices, and dispersion relations.

3. Impact on Phonon Dispersion and Spectral Properties

Flexoelectric phonon coupling fundamentally modifies soft optic and acoustic phonon dispersions in ferroics (Morozovska et al., 23 Jul 2025, Morozovska et al., 8 Mar 2025, Morozovska et al., 2015). The key outcome is the hybridization and repulsion ("pushing away") between optical and acoustic branches, non-diagonalization of the generalized susceptibility tensor, and broadening of the $k$-spectrum:

• The acoustic branch acquires a strong $f^2 k^4$ term and may soften at finite wavevector $k_{\text{cr}}$ for $f > f_{\text{cr}}$ (critical flexoelectric coefficient), signaling an instability towards a spatially modulated phase (SMP) (Morozovska et al., 2017, Morozovska et al., 2016).
• Flexocoupling-induced "ferrons", collective dipolar fluctuations, emerge as additional branches with spectral density $S_A(k) \propto 1/[\omega_A^2(k)]$, diverging at the soft-mode instability (Morozovska et al., 8 Mar 2025).
• Nonlinear contributions via Landau coefficients and electrostriction further renormalize gaps and mixing.

Quantitatively, in van der Waals CuInP$_2$S$_6$, $f_{55}^{\text{cr}} \sim 5$ V at $T = 293$ K, and the modulation period at instability is $\lambda_{\text{mod}} \sim 20$ nm (Morozovska et al., 8 Mar 2025). The frequency of acoustic flexophonons and ferrons can vanish at nonzero $k$ and threshold electric field, enabling experimental extraction of flexoelectric tensor components (Morozovska et al., 23 Jul 2025).

4. Instabilities and Spatially Modulated Phases

Nonlinear flexoelectric phonon coupling is the mechanism underlying the appearance of spatially modulated (incommensurate) polar phases in ferroics (Morozovska et al., 2017, Morozovska et al., 2016). The critical coupling $f_{\text{cr}}(T)$ for SMP onset is determined by temperature, order parameter stiffness, strain, and higher-order gradient coefficients: $$f_{\text{cr}}^2(T) = a_s v + c g + 2 \sqrt{c g a_s v}$$ For $f \geq f_{\text{cr}}$, the homogeneous state bifurcates to a sinusoidal order parameter modulation $P(x) = A(T)\sin(k_{\text{cr}} x)$, with $A(T) \propto (T_{\text{IC}} - T)^{1/2}$. The instability produces a square-root dispersion near $k_{\text{cr}}$, a hallmark observed in neutron and Raman scattering experiments (Morozovska et al., 2017).

The generalized susceptibility in the SMP regime exhibits a kink at $T_{\text{IC}}$ rather than a Curie-like divergence. Inclusion of squired elastic strain gradient terms yields a temperature-dependent upper bound for flexoelectric coupling strength, superseding previous static bulk limits (Morozovska et al., 2016).

5. Generalized Susceptibility and Correlation Functions

Joint static and dynamic flexoelectric couplings induce off-diagonal components in the generalized susceptibility $\chi_{ij}(k, \omega)$, proportional to convolutions of spontaneous polarization and flexocoupling tensors (Morozovska et al., 2015). The static susceptibility spectrum is broadened and correlation radii $R_{ij}$ reduced: $$R_{ij}^2 = \frac{g_{ij} C_{ij} - f_{ij}^2}{C_{ij} \alpha_{ij}}$$ Ferroelectric nonlinearity and electrostriction can further tune the broadening or narrowing of $\chi(k)$; non-degenerate TO/TA branches in the ferroelectric state also reflect mode splitting and phonon repulsion.

The fluctuation–dissipation theorem yields order-parameter correlation functions with full tensorial structure, encapsulating all static, dynamic, and anisotropic contributions (Morozovska et al., 2015).

6. Experimental Signatures and Applications

Flexoelectric phonon coupling drives observable phenomena in ferroelectrics, multiferroics, and van der Waals ferrielectrics:

• Softening/dip of the acoustic phonon branch at finite $k$ in Brillouin or neutron scattering (Morozovska et al., 2017, Morozovska et al., 23 Jul 2025).
• Giant enhancement of low-energy Raman/neutron intensity at softening points (Morozovska et al., 23 Jul 2025).
• Finite-$k$ superlattice peaks in diffraction and SHG microscopy, indicative of SMP phases.
• Enhanced pyroelectric and electrocaloric responses sourced by flexophonon/flexoferron contributions; $\Pi \sim 10^{-4}$–$10^{-3}$ C/m$^2$K at low $T$ for CuInP$_2$S$_6$ (Morozovska et al., 8 Mar 2025).
• Nonlinear and field-tunable electromechanical response in mixed ionic–electronic conductors, with measurable bias and frequency harmonics (Morozovska et al., 2010).

For Pb(Zr$_{0.4}$Ti$_{0.6}$)O$_3$ and PbTiO$_3$, fitting the phonon dispersion and susceptibility spectra requires inclusion of both static $f_{ij}$ and dynamic $M_{ij}$ flexoelectric tensors for quantitative agreement with experiment (Morozovska et al., 2015).

7. Boundary Condition Sensitivity and Bulk Invariants

The symmetric (flexoelectric-like) coupling tensor $K$ is generally sensitive to the choice of electrostatic boundary conditions and gauge conventions, complicating the definition of "bulk" response (Stengel, 2023). In contrast, the antisymmetric (DM-like) scalar $\zeta$ remains invariant under changes in boundary conditions and represents a genuine transverse bulk property. Only after flexoelastic renormalization does $K$ contribute unambiguously to the physically relaxed macroscopic polarization.

This distinction has implications for first-principles extraction, transport modeling, and interpretation of experiments where boundary and sample geometry effects are significant.

---

Nonlinear flexoelectric phonon coupling unifies trilinear gradient interactions in ferroic crystals, establishes the criteria for modulated phase formation, controls phonon spectrum topology, and governs temperature, field, and symmetry dependencies of macroscopic polarization and dielectric response. It is a vital ingredient in the emerging understanding of nanostructured ferroics, flexo-engineering, and field-tunable dielectric and caloric functionality (Stengel, 2023, Morozovska et al., 23 Jul 2025, Morozovska et al., 8 Mar 2025, Morozovska et al., 2017, Morozovska et al., 2016, Morozovska et al., 2015, Morozovska et al., 2010).