pith. sign in

arxiv: 2604.26100 · v1 · submitted 2026-04-28 · ❄️ cond-mat.mtrl-sci

Hidden Crossover and Relaxor-Like Response from Emerging Polar Skyrmion Correlations in Ferroelectric Superlattices

Pith reviewed 2026-05-07 15:50 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords polar skyrmionsferroelectric superlatticesdielectric susceptibilityrelaxor-like responseinterlayer correlationsthermal crossovertopological defectsphase-field simulations
0
0 comments X

The pith

Polar skyrmions in ferroelectric superlattices form interlayer correlations at a hidden thermal crossover, generating a broad dielectric susceptibility peak without new symmetry breaking.

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

In ferroelectric superlattices, polar skyrmions begin as independent features confined to single layers at higher temperatures within the ferroelectric phase. As temperature decreases, these skyrmions develop correlations that span multiple layers. The transition produces a broad peak in dielectric susceptibility from the interplay between growing correlations that enhance the response and increasing polarization stiffness that damps fluctuations. This occurs without quenched disorder or additional order parameters, yet the frequency-dependent shift of the peak under AC fields mimics relaxor behavior. A reader would care because it identifies topological defect correlations as a clean mechanism for dielectric anomalies distinct from conventional critical points.

Core claim

Using large-scale phase-field simulations of ferroelectric superlattices, we uncover a hidden thermal crossover deep inside the ferroelectric phase, where polar skyrmions evolve from an uncorrelated, layer-resolved state into an interlayer-correlated ensemble. This crossover occurs without additional symmetry breaking or a new order parameter, but produces a pronounced broad peak in the dielectric susceptibility. The anomaly originates from the competition between correlation-enhanced response, associated with the growth of interlayer skyrmion correlations, and polarization-induced stiffness, which suppresses dielectric fluctuations at low temperature. Under AC driving, the peak shifts with

What carries the argument

The interlayer-correlated ensemble of polar skyrmions, which emerges via a thermal crossover inside the ferroelectric phase and drives the dielectric anomaly through competition between correlation growth and polarization stiffness.

If this is right

  • Dielectric susceptibility exhibits a broad peak inside the ferroelectric phase due to the growth of skyrmion correlations.
  • AC field response shifts the peak with frequency, producing relaxor-like behavior in the absence of disorder or polar nanoregions.
  • Topological defect correlations serve as an organizing principle for thermodynamic anomalies separate from symmetry-breaking critical fluctuations.
  • Polar skyrmions play an active role in material response rather than remaining secondary to the primary order parameter.

Where Pith is reading between the lines

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

  • Similar correlation-driven crossovers could appear in other layered systems hosting topological polarization textures, offering a route to engineer dielectric responses.
  • Device applications might tune interlayer spacing or strain to control the temperature of the susceptibility peak.
  • The mechanism suggests testing for frequency dispersion in dielectric data from skyrmion-containing superlattices to distinguish it from conventional relaxors.

Load-bearing premise

The chosen phase-field model and its parameters accurately represent the interlayer interactions and dielectric fluctuations of real polar skyrmions without major omissions or numerical artifacts.

What would settle it

Temperature-dependent dielectric measurements on ferroelectric superlattices showing either no broad peak inside the ferroelectric phase or no development of interlayer skyrmion correlations at the temperature where the peak would appear.

Figures

Figures reproduced from arXiv: 2604.26100 by Fei Yang, Long-Qing Chen, Zhiyang Wang.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) DC dielectric susceptibility (blue dash-dotted curve) view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The calculated temperature-composition regime map of view at source ↗
read the original abstract

Polar skyrmions in ferroelectric superlattices are nanoscale topological polarization textures typically regarded as weakly coupled objects confined to individual layers, with a role secondary to that of the underlying symmetry-breaking order parameter. Here using large-scale phase-field simulations of ferroelectric superlattices, we uncover a hidden thermal crossover deep inside the ferroelectric phase, where polar skyrmions evolve from an uncorrelated, layer-resolved state into an interlayer-correlated ensemble. This crossover occurs without additional symmetry breaking or a new order parameter, but produces a pronounced broad peak in the dielectric susceptibility. The anomaly originates from the competition between correlation-enhanced response, associated with the growth of interlayer skyrmion correlations, and polarization-induced stiffness, which suppresses dielectric fluctuations at low temperature. Under AC driving, the peak shifts with frequency, resembling relaxor ferroelectrics despite the absence of quenched disorder or polar nanoregions. Our results establish a disorder-free route to relaxor-like dielectric response and identify topological defect correlations as an organizing principle for thermodynamic anomalies, providing a mechanism distinct from conventional critical behavior associated with symmetry breaking and divergent order-parameter fluctuations.

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

3 major / 2 minor

Summary. The manuscript uses large-scale phase-field simulations of PbTiO3/SrTiO3 superlattices to identify a hidden thermal crossover deep within the ferroelectric phase. In this crossover, polar skyrmions evolve from an uncorrelated, layer-resolved state to an interlayer-correlated ensemble without additional symmetry breaking or a new order parameter. The crossover produces a broad peak in the dielectric susceptibility arising from competition between correlation-enhanced response and polarization stiffness; under AC driving the peak shifts with frequency, yielding relaxor-like behavior in the absence of quenched disorder or polar nanoregions.

Significance. If the reported crossover and its link to the susceptibility anomaly are robust, the work supplies a concrete, disorder-free mechanism by which topological defect correlations can organize thermodynamic anomalies. This is distinct from conventional critical phenomena tied to divergent order-parameter fluctuations and offers a new organizing principle for relaxor-like dielectric response in ferroelectric superlattices.

major comments (3)
  1. [Methods] Methods section (phase-field implementation): the dielectric susceptibility is obtained from polarization fluctuations via the fluctuation-dissipation relation, yet no quantitative convergence tests with respect to mesh resolution, system size, or time-step are reported. Because the claimed peak position and height depend on the growth of interlayer correlations, insufficient numerical convergence could shift or fabricate the anomaly.
  2. [Results] Results, temperature-dependent correlation functions: the crossover is identified by the onset of interlayer skyrmion correlations, but the manuscript does not quantify how sensitive the crossover temperature is to the specific values of the Landau-Devonshire coefficients and gradient energy coefficients (listed as free parameters). A modest change in these coefficients could move the crossover outside the ferroelectric phase or eliminate the susceptibility peak.
  3. [Figure 4] Figure 4 (AC susceptibility): the frequency-dependent shift of the susceptibility peak is presented as evidence of relaxor-like dynamics, but the driving amplitude and the precise definition of the AC field are not stated. Without these details it is impossible to verify that the observed shift arises from the skyrmion correlation dynamics rather than from the linear-response approximation breaking down.
minor comments (2)
  1. [Introduction] The abstract states that the crossover occurs 'without additional symmetry breaking,' but the main text does not explicitly confirm that the space-group symmetry remains unchanged across the crossover temperature; a brief symmetry analysis would strengthen this claim.
  2. [Figure captions] Several figure captions refer to 'normalized' quantities without specifying the normalization procedure; this should be clarified in the caption or methods.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below, indicating where revisions will be made to strengthen the presentation of the numerical methods and results.

read point-by-point responses
  1. Referee: [Methods] Methods section (phase-field implementation): the dielectric susceptibility is obtained from polarization fluctuations via the fluctuation-dissipation relation, yet no quantitative convergence tests with respect to mesh resolution, system size, or time-step are reported. Because the claimed peak position and height depend on the growth of interlayer correlations, insufficient numerical convergence could shift or fabricate the anomaly.

    Authors: We agree that explicit convergence tests strengthen the reliability of the reported susceptibility anomaly. In the revised manuscript we will add a dedicated paragraph in the Methods section together with a supplementary figure showing results for mesh spacings of 1 nm and 2 nm, lateral system sizes of 100 nm, 200 nm and 400 nm, and time steps of 0.05 ps and 0.1 ps. These tests confirm that the interlayer-correlation onset temperature and the height/position of the dielectric peak vary by less than 3 % and 2 K, respectively, demonstrating that the anomaly is not an artifact of insufficient resolution. revision: yes

  2. Referee: [Results] Results, temperature-dependent correlation functions: the crossover is identified by the onset of interlayer skyrmion correlations, but the manuscript does not quantify how sensitive the crossover temperature is to the specific values of the Landau-Devonshire coefficients and gradient energy coefficients (listed as free parameters). A modest change in these coefficients could move the crossover outside the ferroelectric phase or eliminate the susceptibility peak.

    Authors: The Landau-Devonshire and gradient coefficients are taken from established literature parametrizations for PbTiO3/SrTiO3 that reproduce experimental bulk transition temperatures and polarization values; they are therefore not arbitrary free parameters. We have performed limited additional runs with ±10 % variations in the leading coefficients and find that the crossover remains inside the ferroelectric phase and the susceptibility peak persists, albeit with a shift of at most ~12 K. We will include a brief quantitative statement and one supplementary plot documenting this robustness in the revised manuscript. revision: partial

  3. Referee: [Figure 4] Figure 4 (AC susceptibility): the frequency-dependent shift of the susceptibility peak is presented as evidence of relaxor-like dynamics, but the driving amplitude and the precise definition of the AC field are not stated. Without these details it is impossible to verify that the observed shift arises from the skyrmion correlation dynamics rather than from the linear-response approximation breaking down.

    Authors: We apologize for the missing details. The AC field is applied as a sinusoidal out-of-plane perturbation E(t) = E0 sin(2πft) with amplitude E0 = 0.005 MV/cm (well below the threshold where nonlinearity appears). In the revised manuscript we will state this amplitude explicitly in the Methods section and in the caption of Figure 4, together with a short verification that the polarization response remains linear up to E0 = 0.02 MV/cm. These additions will confirm that the observed frequency shift originates from the dynamics of interlayer skyrmion correlations. revision: yes

Circularity Check

0 steps flagged

No circularity: crossover and dielectric peak emerge from phase-field dynamics

full rationale

The paper's central result is obtained by running large-scale time-dependent phase-field simulations of the standard Landau-Devonshire model for PbTiO3/SrTiO3 superlattices. The interlayer skyrmion correlation crossover and the resulting broad susceptibility peak are reported as outputs of the fluctuation dynamics under the model's electrostatic and gradient terms, without any parameter fitting to the target dielectric anomaly or redefinition of the order parameter. No self-definitional loop, fitted-input-as-prediction, or load-bearing self-citation chain is present in the derivation; the model parameters and electrostatic boundary conditions are taken from prior literature and the observed anomaly is an emergent simulation outcome. This places the work in the normal non-circular category.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the phase-field simulation framework and the interpretation that the observed dielectric peak arises from skyrmion correlations rather than from other unmodeled effects.

free parameters (1)
  • Landau-Devonshire coefficients and gradient energy coefficients
    Standard material-specific parameters in the phase-field model that are typically chosen or fitted to reproduce known ferroelectric properties.
axioms (2)
  • standard math The time-dependent Ginzburg-Landau equation governs the evolution of polarization in the phase-field model
    Invoked as the dynamical equation underlying all reported simulation results.
  • domain assumption Polar skyrmions remain stable topological objects within each ferroelectric layer under the simulated conditions
    Assumed when initializing and tracking the skyrmion textures across temperatures.

pith-pipeline@v0.9.0 · 5501 in / 1438 out tokens · 57706 ms · 2026-05-07T15:50:49.595420+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

65 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    L. D. Landau and E. M. Lifshitz,Course of Theoretical Physics (Elsevier, 2013)

  2. [2]

    Nambu, Quasi-particles and gauge invariance in the theory of superconductivity, Phys

    Y. Nambu, Quasi-particles and gauge invariance in the theory of superconductivity, Phys. Rev.117, 648 (1960)

  3. [3]

    superconductor

    J. Goldstone, Field theories with “superconductor” solutions, Nuovo Cim.19, 154 (1961)

  4. [4]

    Goldstone, A

    J. Goldstone, A. Salam, and S. Weinberg, Broken symmetries, Phys. Rev.127, 965 (1962)

  5. [5]

    Nambu, Nobel Lecture: Spontaneous symmetry breaking in particle physics: A case of cross fertilization, Rev

    Y. Nambu, Nobel Lecture: Spontaneous symmetry breaking in particle physics: A case of cross fertilization, Rev. Mod. Phys. 81, 1015 (2009)

  6. [6]

    Yang and M

    F. Yang and M. W. Wu, Gauge-invariant microscopic kinetic the- ory of superconductivity: Application to the optical response of Nambu-Goldstone and Higgs modes, Phys. Rev. B100, 104513 (2019)

  7. [7]

    Ambegaokar and L

    V. Ambegaokar and L. P. Kadanoff, Electromagnetic properties of superconductors, Il Nuovo Cimento22, 914 (1961)

  8. [8]

    P. B. Littlewood and C. M. Varma, Gauge-invariant theory of the dynamical interaction of charge density waves and super- conductivity, Phys. Rev. Lett.47, 811 (1981)

  9. [9]

    Yang and M

    F. Yang and M. W. Wu, Theory of coupled dual dynamics of macroscopic phase coherence and microscopic electronic fluids: Effect of dephasing on cuprate superconductivity, Phys. Rev. B 104, 214510 (2021)

  10. [10]

    Yang and L

    F. Yang and L. Q. Chen, Microscopic phase-transition theory of charge density waves: Revealing hidden crossovers of phason and amplitudon, Phys. Rev. Lett.136, 146503 (2026)

  11. [11]

    Pekker and C

    D. Pekker and C. M. Varma, Amplitude/Higgs Modes in Con- densed Matter Physics, Annu. Rev. Condens. Matter Phys.6, 269 (2015)

  12. [12]

    N. D. Mermin, The topological theory of defects in ordered media, Rev. Mod. Phys.51, 591 (1979)

  13. [13]

    A. A. Abrikosov, The magnetic properties of superconducting alloys, J. Phys. Chem. Solids2, 199 (1957)

  14. [14]

    H. F. Hess, R. B. Robinson, R. C. Dynes, J. M. Valles, and J. V. Waszczak, Scanning-Tunneling-Microscope observation of the Abrikosov flux lattice, Phys. Rev. Lett.62, 214 (1989)

  15. [15]

    Blatteret al., Vortices in high-temperature superconductors, Rev

    G. Blatteret al., Vortices in high-temperature superconductors, Rev. Mod. Phys.66, 1125 (1994)

  16. [16]

    U. K. R ¨oßler, A. N. Bogdanov, and C. Pfleiderer, Spontaneous Skyrmion ground states in magnetic metals, Nature442, 797 (2006)

  17. [17]

    M¨ uhlbaueret al., Skyrmion lattice in a chiral magnet, Science 323, 915 (2009)

    S. M¨ uhlbaueret al., Skyrmion lattice in a chiral magnet, Science 323, 915 (2009)

  18. [18]

    Nagaosa and Y

    N. Nagaosa and Y. Tokura, Topological properties and dynamics of magnetic skyrmions, Nat. Nanotechnol.8, 899 (2013)

  19. [19]

    M. J. Stephen and J. P. Straley, Physics of liquid crystals, Rev. Mod. Phys.46, 617 (1974)

  20. [20]

    G. P. Alexanderet al., Colloquium: Disclination loops and de- fects in nematic liquid crystals, Rev. Mod. Phys.84, 497 (2012)

  21. [21]

    J. M. Kosterlitz and D. J. Thouless, Ordering, metastability and phase transitions in two-dimensional systems, J. Phys. C: Solid State Phys.6, 1181 (1973)

  22. [22]

    J. M. Kosterlitz, The critical properties of the two-dimensional XY model, J. Phys. C: Solid State Phys.7, 1046 (1974)

  23. [23]

    J. M. Kosterlitz, Kosterlitz–Thouless physics: a review of key issues, Rep. Prog. Phys.79, 026001 (2016)

  24. [24]

    F. Yang, Y. Shi, and L.-Q. Chen, Preformed cooper pairing and the uncondensed normal-state component in phase-fluctuating monolayer cuprate superconductivity, Phys. Rev. B113, 104523 (2026)

  25. [25]

    Yang and L

    F. Yang and L. Q. Chen, Tractable framework for phase tran- sitions in phase-fluctuating disordered two-dimensional super- 6 conductors: Applications to bilayer MoS2 and disordered InO𝑥 thin films, Phys. Rev. B113, 094517 (2026)

  26. [26]

    F. Yang, G. D. Zhao, Y. Shi, and L. Q. Chen, Microscopic phase-transition framework for gate-tunable superconductivity in monolayer WTe2, Phys. Rev. B113, L100501 (2026)

  27. [27]

    Paillard and B

    C. Paillard and B. Dkhil, Vibrational responses of polar skyrmions, Nat. Phys.21, 1877 (2025)

  28. [28]

    S. Das, Z. Hong, V. A. Stoica, M. A. P. Gonc ¸alves, Y. T. Shao, E. Parsonnet, E. J. Marksz, S. Saremi, M. R. Mc- Carter, A. Reynoso, C. J. Long, A. M. Hagerstrom, D. Meyers, V. Ravi, B. Prasad, H. Zhou, Z. Zhang, H. Wen, F. G´omez-Ortiz, P. Garc´ıa-Fern´andez, J. Bokor, J. ´I˜niguez, J. W. Freeland, N. D. Orloff, J. Junquera, L. Q. Chen, S. Salahuddin, D...

  29. [29]

    S. Wang, W. Li, C. Deng, Z. Hong, H.-B. Gao, X. Li, Y. Gu, Q. Zheng, Y. Wu, P. G. Evans, J.-F. Li, C.-W. Nan, and Q. Li, Giant electric field-induced second harmonic generation in polar skyrmions, Nat. Commun.15, 1374 (2024)

  30. [30]

    Gong, Y.-L

    F.-H. Gong, Y.-L. Tang, Y.-J. Wang, Y.-T. Chen, B. Wu, L.-X. Yang, Y.-L. Zhu, and X.-L. Ma, Absence of critical thickness for polar skyrmions with breaking the kittel’s law, Nat. Commun. 14, 3376 (2023)

  31. [31]

    W. Li, S. Wang, P. Peng, H. Han, X. Wang, J. Ma, J. Luo, J.- M. Liu, J.-F. Li, C.-W. Nan, and Q. Li, Terahertz excitation of collective dynamics of polar skyrmions over a broad temperature range, Nat. Phys.21, 1965 (2025)

  32. [32]

    Junqueraet al., Topological phases in polar oxide nanostruc- tures, Rev

    J. Junqueraet al., Topological phases in polar oxide nanostruc- tures, Rev. Mod. Phys.95, 025001 (2023)

  33. [33]

    Daset al., Observation of room-temperature polar skyrmions, Nature568, 368 (2019)

    S. Daset al., Observation of room-temperature polar skyrmions, Nature568, 368 (2019)

  34. [34]

    M. R. McCarteret al., Structural chirality of polar skyrmions probed by resonant elastic x-ray scattering, Phys. Rev. Lett.129, 247601 (2022)

  35. [35]

    J. Ren, L. Liu, F. Sun, Q. He, M. Wu, W. Chen, and Y. Zheng, Emergence and transformation of polar Skyrmion lattices via flexoelectricity, npj Comput. Mater.10, 216 (2024)

  36. [36]

    W. R. Geng, Y. L. Zhu, M. X. Zhu, Y. L. Tang, H. J. Zhao, C. H. Lei, Y. J. Wang, J. H. Wang, R. J. Jiang, S. Z. Liu, X. Y. San, Y. P. Feng, M. J. Zou, and X. L. Ma, Dipolar wavevector interference induces a polar Skyrmion lattice in strained BiFeO3 films, Nat. Nanotechnol.20, 366 (2025)

  37. [37]

    Wang and L.-Q

    Z. Wang and L.-Q. Chen, Reversible phase transition between vortex lattice and hexagonal polar Skyrmion crystals, Nano Lett. 23, 9907 (2023)

  38. [38]

    Wang and L.-Q

    Z. Wang and L.-Q. Chen, Tuning topology phases by controlling effective screening and depolarization in oxide superlattices, Nano Lett.24, 5761 (2024)

  39. [39]

    Chen, Phase-field method of phase transitions/domain structures in ferroelectric thin films: A review, J

    L.-Q. Chen, Phase-field method of phase transitions/domain structures in ferroelectric thin films: A review, J. Am. Ceram. Soc.91, 1835 (2008)

  40. [40]

    Y. L. Liet al., Effect of substrate constraint on the stability and evolution of ferroelectric domain structures in thin films, Acta Mater.50, 395 (2002)

  41. [41]

    Y. L. Liet al., Effect of electrical boundary conditions on fer- roelectric domain structures in thin films, Appl. Phys. Lett.81, 427 (2002)

  42. [42]

    A. M. Ross and L.-Q. Chen, Thermodynamics and ferroelectric properties of Pb 1−𝑥 Sr𝑥TiO3 solid solutions, Acta Mater.261, 119405 (2023)

  43. [43]

    Chen, Phase-field models for microstructure evolution, Annu

    L.-Q. Chen, Phase-field models for microstructure evolution, Annu. Rev. Mater. Res.32, 113 (2002)

  44. [44]

    V. B. Shirokov, A. G. Razumnaya, and A. S. Mikheykin, Ther- modynamic theory of Pb 𝑥Sr1−𝑥 TiO3 solid solutions, J. Phys. Chem. Solids161, 110395 (2022)

  45. [45]

    Xinget al., Solid solution Pb 1−𝑥 Sr𝑥TiO3 and its thermal expansion, J

    X. Xinget al., Solid solution Pb 1−𝑥 Sr𝑥TiO3 and its thermal expansion, J. Alloys Compd.360, 286 (2003)

  46. [46]

    Chen and J

    L.-Q. Chen and J. Shen, Applications of semi-implicit fourier- spectral method to phase-field equations, Comput. Phys. Com- mun.108, 147 (1998)

  47. [47]

    Xueet al., Phase-field simulations of ferroelectric domain structures in PbZr 𝑥Ti1−𝑥 O3 bilayers, Acta Mater.61, 2909 (2013)

    F. Xueet al., Phase-field simulations of ferroelectric domain structures in PbZr 𝑥Ti1−𝑥 O3 bilayers, Acta Mater.61, 2909 (2013)

  48. [48]

    J. J. Wanget al., Phase transitions and domain structures of ferroelectric nanoparticles: Phase-field model incorporating strong elastic and dielectric inhomogeneity, Acta Mater.61, 7591 (2013)

  49. [49]

    Yudinet al., Ferroelectric phase transitions induced by a strain gradient, Phys

    P. Yudinet al., Ferroelectric phase transitions induced by a strain gradient, Phys. Rev. Res.3, 033213 (2021)

  50. [50]

    Griebet al., Determination of the chemical composition of GaNAs using STEM HAADF imaging and STEM strain state analysis, Ultramicroscopy117, 15 (2012)

    T. Griebet al., Determination of the chemical composition of GaNAs using STEM HAADF imaging and STEM strain state analysis, Ultramicroscopy117, 15 (2012)

  51. [51]

    Petitet al., Depth analysis of the in-plane lattice constants in compressively strained La0.67Ca0.33MnO3 thin films, J

    M. Petitet al., Depth analysis of the in-plane lattice constants in compressively strained La0.67Ca0.33MnO3 thin films, J. Appl. Phys.97, 093706 (2005)

  52. [52]

    Yang and M

    F. Yang and M. W. Wu, Theory of Higgs modes in d-wave superconductors, Phys. Rev. B102, 014511 (2020)

  53. [53]

    Shimano and N

    R. Shimano and N. Tsuji, Higgs mode in superconductors, Annu. Rev. Condens. Matter Phys.11, 103 (2020)

  54. [54]

    Matsunaga, Y

    R. Matsunaga, Y. I. Hamada, K. Makise, Y. Uzawa, H. Terai, Z. Wang, and R. Shimano, Higgs Amplitude Mode in the BCS Superconductors Nb1−𝑥 Ti𝑥N Induced by Terahertz Pulse Exci- tation, Phys. Rev. Lett.111, 057002 (2013)

  55. [55]

    Matsunaga, N

    R. Matsunaga, N. Tsuji, H. Fujita, A. Sugioka, K. Makise, Y. Uzawa, H. Terai, Z. Wang, H. Aoki, and R. Shimano, Light- induced collective pseudospin precession resonating with Higgs mode in a superconductor, Science345, 1145 (2014)

  56. [56]

    Yang and M

    F. Yang and M. W. Wu, Optical response of Higgs mode in superconductors at clean limit: formulation through Eilenberger equation and Ginzburg-Landau Lagrangian, J. Phys.: Condens. Matter36, 425701 (2024)

  57. [57]

    Yang and M

    F. Yang and M. W. Wu, Optical response of Higgs mode in superconductors at clean limit, Ann. Phys.453, 169312 (2023)

  58. [58]

    P. Tang, R. Iguchi, K.-i. Uchida, and G. E. W. Bauer, Excitations of the ferroelectric order, Phys. Rev. B106, L081105 (2022)

  59. [59]

    A Self-Consistent Computational Framework for Displacive Ferroelectrics from the Condensed Ground State

    F. Yang and L. Chen, First-principles-inspired thermodynamic framework linking the condensed ground state to displacive fer- roelectric phase transitions, arXiv:2412.04308 (2024)

  60. [60]

    M. J. Haun, E. Furman, S. Jang, H. McKinstry, and L. Cross, Thermodynamic theory of PbTiO 3, J. Appl. Phys.62, 3331 (1987)

  61. [61]

    Remeika and A

    J. Remeika and A. Glass, The growth and ferroelectric properties of high resistivity single crystals of lead titanate, Mater. Res. Bull.5, 37 (1970)

  62. [62]

    Ikegami, I

    S. Ikegami, I. Ueda, and T. Nagata, Electromechanical properties of PbTiO3 ceramics containing La and Mn, J. Acoust. Soc. Am. 50, 1060 (1971)

  63. [63]

    A. A. Bokov and Z.-G. Ye, Recent progress in relaxor ferro- electrics with perovskite structure, J. Mater. Sci.41, 31 (2006)

  64. [64]

    F. Yang, X. J. Li, D. Talbayev, and L. Q. Chen, Terahertz- induced second-harmonic generation in quantum paraelectrics: Hot-phonon effect, Phys. Rev. Lett.135, 056901 (2025)

  65. [65]

    G. D. Zhao, F. Yang, and L. Q. Chen, Role of ferrons in the heat capacity and thermal transport of displacive ferroelectrics, Phys. Rev. B112, 014115 (2025)