Anharmonic phonon coupling enabled by local inversion symmetry breaking at domain walls in ferroelastics

Andrey Shalit; Martina Basini; Mattia Udina; Niccolo Sellati; Paolo Barone; Peter Hamm; Seyyed Jabbar Mousavi; Steven L. Johnson; Thomas Feurer; Vivek Unikandanunni

arxiv: 2604.28066 · v2 · submitted 2026-04-30 · ❄️ cond-mat.mtrl-sci

Anharmonic phonon coupling enabled by local inversion symmetry breaking at domain walls in ferroelastics

Seyyed Jabbar Mousavi , Vivek Unikandanunni , Niccolo Sellati , Paolo Barone , Martina Basini , Steven L. Johnson , Andrey Shalit , Peter Hamm

show 2 more authors

Mattia Udina Thomas Feurer

This is my paper

Pith reviewed 2026-05-07 06:18 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci

keywords LaAlO3ferroelastic domain wallsinversion symmetry breakinganharmonic phonon couplingtwo-dimensional Raman-terahertz spectroscopyphonon selection rulestwin domains

0 comments

The pith

Local inversion symmetry breaking at domain walls enables anharmonic phonon coupling in ferroelastic LaAlO3.

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

Ferroelastic materials like LaAlO3 form twin domains through spontaneous symmetry breaking, keeping the bulk crystal centrosymmetric while locally breaking inversion symmetry at the domain walls. This local breaking alters phonon selection rules and allows the Eg phonon to gain finite infrared activity. The paper uses two-dimensional Raman-terahertz spectroscopy to detect cross-peaks that signal both mechanical and electrical anharmonicity between this Eg mode and the A1g Raman-active phonon. A reader would care because it shows how interface features in domain-structured materials control lattice vibrations that bulk measurements overlook. The work highlights a new route to detect subtle symmetry effects through weak anharmonic signals.

Core claim

Direct evidence of anharmonic phonon coupling is reported in ferroelastic LaAlO3 using two-dimensional Raman-terahertz spectroscopy. The observed cross-peaks arise from both mechanical and electrical anharmonicity between the A1g Raman-active phonon and the Eg phonon. The Eg phonon acquires finite infrared activity through local inversion symmetry breaking at ferroelastic domain walls.

What carries the argument

Local inversion symmetry breaking at ferroelastic domain walls, which activates the Eg phonon in the infrared and enables its anharmonic coupling to the A1g Raman mode.

If this is right

Ferroelastic domain walls host anharmonic phonon interactions forbidden in the bulk crystal.
Two-dimensional Raman-terahertz spectroscopy can detect subtle symmetry breaking through intrinsically weak anharmonic signals.
Phonon selection rules change locally at domain walls, altering which modes couple.
Both mechanical and electrical anharmonicity contribute to the observed cross-peaks in such systems.

Where Pith is reading between the lines

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

The same spectroscopy approach could map domain wall density or dynamics by tracking cross-peak strength under varying conditions.
Similar local symmetry effects may appear in other ferroelastic or twinned materials and influence their thermal or dielectric response.
Varying temperature to move domain walls would provide a direct test of whether the coupling strength tracks wall area.

Load-bearing premise

The observed cross-peaks arise specifically from anharmonic coupling enabled by local inversion symmetry breaking at domain walls rather than from bulk effects, experimental artifacts, or other unaccounted mechanisms.

What would settle it

No cross-peaks in spectra from a single-domain LaAlO3 sample or from regions without domain walls would falsify the attribution.

Figures

Figures reproduced from arXiv: 2604.28066 by Andrey Shalit, Martina Basini, Mattia Udina, Niccolo Sellati, Paolo Barone, Peter Hamm, Seyyed Jabbar Mousavi, Steven L. Johnson, Thomas Feurer, Vivek Unikandanunni.

**Figure 1.** Figure 1: FIG. 1. (a) Schematic of the experimental geometry for hy view at source ↗

**Figure 2.** Figure 2: FIG. 2. (a) Experimental 2D Raman-THz time-domain data measured using the RTT pulse sequence for the LAO crystal at view at source ↗

**Figure 3.** Figure 3: FIG. 3. (a,b) Diagrammatic representations of the third-order processes involving (a) three-phonon mechanical anharmonicity view at source ↗

**Figure 4.** Figure 4: FIG. 4. (a) Evolution across a DW (at X=0 view at source ↗

**Figure 6.** Figure 6: FIG. 6. (a) Sketch of the contributions to the 2D response view at source ↗

**Figure 7.** Figure 7: FIG. 7. Raman pulse fluence dependence of the (3.7 THz, 1 view at source ↗

**Figure 9.** Figure 9: FIG. 9. (a) Diagrammatic representation of the ISRS-like view at source ↗

**Figure 10.** Figure 10: FIG. 10. (a) Supercell used for simulating two symmetric view at source ↗

**Figure 11.** Figure 11: FIG. 11. Measured THz radiation from the 50 view at source ↗

read the original abstract

In ferroelastic materials, spontaneous symmetry breaking leads to the formation of twin domains. Although the bulk crystal typically remains centrosymmetric, inversion symmetry can be locally broken at the domain walls, potentially changing phonon selection rules and enabling local anharmonic phonon coupling. Here we report direct evidence of such anharmonic coupling in ferroelastic LaAlO$_3$ using two-dimensional Raman-terahertz spectroscopy. We attribute the cross-peaks observed in the two-dimensional spectra to both mechanical and electrical anharmonicity between the $A_{1g}$ Raman-active phonon and the $E_g$ phonon, which acquires finite infrared activity through local inversion symmetry breaking at ferroelastic domain walls. These findings provide new insight into the complex lattice dynamics of ferroelastic materials and highlight the potential of two-dimensional Raman-terahertz spectroscopy to uncover subtle symmetry breaking through the detection of intrinsically weak anharmonic signals.

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.

Desk Editor's Note private letter to a colleague

The paper detects cross-peaks in 2D Raman-THz spectra of LaAlO3 and attributes them to anharmonic phonon coupling at domain walls, but lacks a quantitative check on whether the small wall volume can explain the signal strength.

read the letter

The paper reports cross-peaks in 2D Raman-THz spectra of LaAlO3 that they tie to anharmonic phonon coupling between the A1g Raman mode and the Eg mode, with the latter gaining IR activity only through local inversion symmetry breaking at ferroelastic domain walls. That is the central observation they present as new. The work applies two-dimensional spectroscopy to isolate weak anharmonic signals that linear measurements would miss, and the assignment to both mechanical and electrical anharmonicity is a reasonable way to frame the coupling. The abstract states the claim cleanly and the spectra are shown as direct evidence. The method itself is the real contribution here: it offers a route to probe subtle local symmetry changes in centrosymmetric bulk crystals without spatial resolution at the wall scale. For groups already using nonlinear THz or Raman techniques on complex oxides, this is a logical next step. The data quality appears standard for the technique, and the choice of LaAlO3 is appropriate given its well-characterized twin domains. The main limitation is the missing link to intensity. Domain walls are atomically thin, so their volume fraction is typically well below 1 percent for micron-scale domains. The paper does not report a measurement of domain density on the actual sample, a comparison to a detwinned reference crystal, or a back-of-the-envelope estimate of the local enhancement factor needed to produce the observed cross-peak strength from that small volume. Without those pieces, bulk anharmonicity, surface contributions, or experimental artifacts remain plausible alternatives. The interpretation therefore rests on the assumption that only the walls enable the Eg IR activity, but the evidence for that exclusivity is not yet quantitative. This is for condensed-matter researchers working on ferroelastics or on ultrafast spectroscopy of symmetry-breaking systems. A reader in that niche would find the experimental approach useful and would want to see the full methods and any domain-structure characterization. I would send it to peer review. The observation is new enough that referees can request the missing controls and decide whether the domain-wall attribution holds.

Referee Report

1 major / 1 minor

Summary. The manuscript reports direct evidence from two-dimensional Raman-terahertz spectroscopy on ferroelastic LaAlO3 for anharmonic phonon coupling. Cross-peaks are attributed to mechanical and electrical anharmonicity between the A1g Raman-active phonon and the Eg phonon, with the latter acquiring finite IR activity due to local inversion symmetry breaking at ferroelastic domain walls.

Significance. If substantiated, the result would be significant for showing how atomically thin domain walls can activate otherwise forbidden anharmonic couplings in centrosymmetric bulk crystals. It also illustrates the sensitivity of 2D Raman-THz spectroscopy for detecting weak signals from local symmetry breaking, with potential extension to other ferroelastic and multiferroic systems.

major comments (1)

[Discussion] Discussion section (cross-peak origin): The central attribution of observed cross-peaks to domain-wall-enabled anharmonicity is not supported by any estimate of the domain-wall volume fraction in the measured crystal or the expected cross-peak intensity scaling with that fraction. In LaAlO3, domain walls are atomically thin and the volume fraction for typical micron-scale domains is ≪1%; without a quantitative consistency check or a detwinned/single-domain reference measurement, bulk anharmonicity, surface effects, or experimental artifacts cannot be excluded. This is load-bearing for the claim in the abstract.

minor comments (1)

[Abstract] The abstract and introduction would benefit from a brief statement of the sample domain structure (e.g., average domain size or wall density) to contextualize the expected signal strength.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for the positive evaluation of the work's significance and for the detailed major comment. We respond point by point below.

read point-by-point responses

Referee: [Discussion] Discussion section (cross-peak origin): The central attribution of observed cross-peaks to domain-wall-enabled anharmonicity is not supported by any estimate of the domain-wall volume fraction in the measured crystal or the expected cross-peak intensity scaling with that fraction. In LaAlO3, domain walls are atomically thin and the volume fraction for typical micron-scale domains is ≪1%; without a quantitative consistency check or a detwinned/single-domain reference measurement, bulk anharmonicity, surface effects, or experimental artifacts cannot be excluded. This is load-bearing for the claim in the abstract.

Authors: We agree that the current manuscript lacks a quantitative estimate of domain-wall volume fraction and intensity scaling, which is a valid concern. In the revised version we will add a dedicated paragraph in the Discussion section providing this analysis. Optical micrographs of our LaAlO3 crystals show typical domain widths of 5–20 μm; assuming domain-wall thicknesses of ~0.5 nm yields a volume fraction of order 10^{-3}–10^{-4}. We will estimate the local enhancement of the mechanical and electrical anharmonic coefficients required to reproduce the observed cross-peak amplitudes and show that an enhancement of 10^2–10^3 is sufficient and physically plausible given the local loss of inversion symmetry. Bulk anharmonicity is ruled out because the Eg mode is strictly IR-inactive in the centrosymmetric bulk, precluding direct THz coupling. Surface contributions are negligible given the penetration depths of the THz and optical beams. We do not possess a detwinned single-domain reference sample, but the mode-specific character of the cross-peaks (A1g–Eg) is difficult to explain by artifacts or surfaces. We will revise the abstract and discussion to reflect these quantitative considerations and to note the absence of a single-domain control as a limitation. revision: partial

standing simulated objections not resolved

We do not have access to a detwinned or single-domain LaAlO3 reference sample for a control measurement.

Circularity Check

0 steps flagged

No circularity; experimental attribution is interpretive, not tautological

full rationale

The paper is a purely experimental report using 2D Raman-THz spectroscopy to observe cross-peaks in LaAlO3 and attribute them to anharmonic coupling between A1g and Eg phonons enabled by local inversion-symmetry breaking at domain walls. No equations, parameter fitting, or derivation chain appear in the provided text. The central claim is a physical interpretation of spectral data grounded in known phonon selection rules and material symmetry properties; it does not reduce to self-definition, fitted inputs renamed as predictions, or load-bearing self-citations. The argument is therefore self-contained against external benchmarks and receives the default non-finding.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on experimental observation and standard interpretation of phonon modes and selection rules in LaAlO3; no free parameters, new entities, or ad-hoc axioms are introduced in the abstract.

axioms (1)

domain assumption Standard phonon mode assignments (A1g Raman-active, Eg) and selection-rule expectations for centrosymmetric LaAlO3.
The attribution of cross-peaks relies on these established mode properties from prior literature on the material.

pith-pipeline@v0.9.0 · 5487 in / 1208 out tokens · 109241 ms · 2026-05-07T06:18:19.795562+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages

[1]

E. K. Salje, Annual Review of Materials Research42, 265 (2012)

work page 2012
[2]

Nataf, M

G. Nataf, M. Guennou, J. Gregg, D. Meier, J. Hlinka, E. Salje, and J. Kreisel, Nature Reviews Physics2, 634 (2020)

work page 2020
[3]

Arlt, Journal of Materials Science25, 2655 (1990)

G. Arlt, Journal of Materials Science25, 2655 (1990)

work page 1990
[4]

Bueble, K

S. Bueble, K. Knorr, E. Brecht, and W. W. Schmahl, Surface science400, 345 (1998)

work page 1998
[5]

E.-J. Guo, R. Desautels, D. Keavney, M. A. Roldan, B. J. Kirby, D. Lee, Z. Liao, T. Charlton, A. Herklotz, T. Zac Ward, et al., Science advances5, eaav5050 (2019)

work page 2019
[6]

Kustov, I

S. Kustov, I. Liubimova, and E. Salje, Applied Physics Letters112(2018)

work page 2018
[7]

Hayward, S

S. Hayward, S. Redfern, and E. Salje, Journal of Physics: Condensed Matter14, 10131 (2002)

work page 2002
[8]

S. A. Hayward, F. D. Morrison, S. A. T. Redfern, E. K. H. Salje, J. F. Scott, K. S. Knight, S. Tarantino, A. M. Glazer, V. Shuvaeva, P. Daniel, M. Zhang, and M. A. Carpenter, Phys. Rev. B72, 054110 (2005)

work page 2005
[9]

M. V. Abrashev, A. P. Litvinchuk, M. N. Iliev, R. L. Meng, V. N. Popov, V. G. Ivanov, R. A. Chakalov, and C. Thomsen, Phys. Rev. B59, 4146 (1999)

work page 1999
[10]

Basini, V

M. Basini, V. Unikandanunni, F. Gabriele, M. Cross, A. M. Derrico, A. X. Gray, M. C. Hoffmann, F. Forte, 6 M. Cuoco, and S. Bonetti, Physical Review Letters136, 156902 (2026)

work page 2026
[11]

Y. Liu, A. Frenkel, G. Garrett, J. Whitaker, S. Fahy, C. Uher, and R. Merlin, Physical review letters75, 334 (1995)

work page 1995
[12]

Yokota, S

H. Yokota, S. Matsumoto, E. Salje, and Y. Uesu, Physical Review B98, 104105 (2018)

work page 2018
[13]

Yokota and Y

H. Yokota and Y. Uesu, Journal of Applied Physics129 (2021)

work page 2021
[14]

Yokota, C

H. Yokota, C. Haines, S. Matsumoto, N. Hasegawa, M. A. Carpenter, Y. Heo, A. Marin, E. Salje, and Y. Uesu, Physical Review B102, 104117 (2020)

work page 2020
[15]

E. K. Salje, M. Alexe, S. Kustov, M. C. Weber, J. Schiemer, G. F. Nataf, and J. Kreisel, Scientific re- ports6, 27193 (2016)

work page 2016
[16]

J. Suda, O. Kamishima, J. Kawamura, T. Hattori, and T. Sato, in Journal of Physics: Conference Series, Vol. 150 (IOP Publishing, 2009) p. 052249

work page 2009
[17]

Hortensius, D

J. Hortensius, D. Afanasiev, A. Sasani, E. Bousquet, and A. Caviglia, npj Quantum Materials5, 95 (2020)

work page 2020
[18]

M. J. Neugebauer, D. M. Juraschek, M. Savoini, P. En- geler, L. Boie, E. Abreu, N. A. Spaldin, and S. L. John- son, Phys. Rev. Res.3, 013126 (2021)

work page 2021
[19]

Ciardi, A

G. Ciardi, A. Berger, P. Hamm, and A. Shalit, The Jour- nal of Physical Chemistry Letters10, 4463 (2019)

work page 2019
[20]

Shalit, S

A. Shalit, S. J. Mousavi, and P. Hamm, The Journal of Physical Chemistry B125, 581 (2021)

work page 2021
[21]

S. J. Mousavi, A. Berger, P. Hamm, and A. Shalit, The Journal of Chemical Physics156, 174501 (2022)

work page 2022
[22]

A. A. Maradudin and R. F. Wallis, Phys. Rev. B2, 4294 (1970)

work page 1970
[23]

R. F. Wallis and A. A. Maradudin, Phys. Rev. B3, 2063 (1971)

work page 2063
[24]

L. B. Humphreys, Phys. Rev. B6, 3886 (1972)

work page 1972
[25]

F¨ orst, C

M. F¨ orst, C. Manzoni, S. Kaiser, Y. Tomioka, Y. Tokura, R. Merlin, and A. Cavalleri, Nature Physics7, 854 (2011)

work page 2011
[26]

Subedi, A

A. Subedi, A. Cavalleri, and A. Georges, Phys. Rev. B 89, 220301 (2014)

work page 2014
[27]

von Hoegen, R

A. von Hoegen, R. Mankowsky, M. Fechner, M. F¨ orst, and A. Cavalleri, Nature555, 79 (2018)

work page 2018
[28]

Khalsa, N

G. Khalsa, N. A. Benedek, and J. Moses, Phys. Rev. X 11, 021067 (2021)

work page 2021
[29]

T. G. Blank, K. A. Grishunin, K. Zvezdin, N. Hai, J. Wu, S.-H. Su, J.-C. Huang, A. Zvezdin, and A. V. Kimel, Physical review letters131, 026902 (2023)

work page 2023
[30]

Petretto, S

G. Petretto, S. Dwaraknath, H. PC Miranda, D. Win- ston, M. Giantomassi, M. J. Van Setten, X. Gonze, K. A. Persson, G. Hautier, and G.-M. Rignanese, Scientific data 5, 1 (2018)

work page 2018
[31]

Miranda, Phonon website: Visualize phonon vibra- tional modes,https://henriquemiranda.github.io/ phononwebsite/(2025), accessed: 2025-11-19

H. Miranda, Phonon website: Visualize phonon vibra- tional modes,https://henriquemiranda.github.io/ phononwebsite/(2025), accessed: 2025-11-19

work page 2025
[32]

A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder, and K. a. Persson, APL Materials1, 011002 (2013)

work page 2013
[33]

Sidler and P

D. Sidler and P. Hamm, The Journal of Chemical Physics 150, 044202 (2019)

work page 2019
[34]

Sidler and P

D. Sidler and P. Hamm, The Journal of Chemical Physics 153, 044502 (2020)

work page 2020
[35]

J. F. Scott, Phys. Rev.183, 823 (1969)

work page 1969
[36]

Lloyd-Hughes, S

J. Lloyd-Hughes, S. Jones, E. Castro-Camus, K. Doig, and J. MacManus-Driscoll, Optics letters39, 1121 (2014)

work page 2014
[37]

S. J. Mousavi, M. F. Biggs, J. A. Johnson, P. Hamm, and A. Shalit, The Journal of Chemical Physics161, 134502 (2024)

work page 2024
[38]

T. P. Martin and L. Genzel, Phys. Status Solidi B61, 493 (1974)

work page 1974
[39]

D. M. Juraschek and S. F. Maehrlein, Phys. Rev. B97, 174302 (2018)

work page 2018
[40]

Giorgianni, M

F. Giorgianni, M. Romani, P. Puphal, M. Isobe, L. Spitz, M. C. Guidi, C. Vicario, and M. Udina, Phys. Rev. B 111, 205138 (2025)

work page 2025
[41]

Sertcan, S

B. Sertcan, S. J. Mousavi, M. Iannuzzi, and P. Hamm, The Journal of Chemical Physics158, 014203 (2023)

work page 2023
[42]

Fiore, N

J. Fiore, N. Sellati, M. Udina, and L. Benfatto, arXiv preprint arXiv:2509.25060 (2025)

work page arXiv 2025
[43]

Proietto, A

F. Proietto, A. Milloch, P. Franceschini, M. Azarm, N. Sellati, R. Mishra, P. C. Moen, S. Smit, M. Bluschke, M. Greven, et al., arXiv preprint arXiv:2603.29713 (2026)

work page arXiv 2026
[44]

A. M. Glazer, Acta Crystallographica Section A31, 756 (1975)

work page 1975
[45]

Gopalan and D

V. Gopalan and D. B. Litvin, Nature Materials10, 376 (2011)

work page 2011
[46]

Bellaiche and J

L. Bellaiche and J. ´I˜ niguez, Phys. Rev. B88, 014104 (2013)

work page 2013
[47]

Barone, D

P. Barone, D. Di Sante, and S. Picozzi, Phys. Rev. B89, 144104 (2014)

work page 2014
[48]

Salje, Ferroelectrics104, 111 (1990)

E. Salje, Ferroelectrics104, 111 (1990)

work page 1990
[49]

Shukla, A

S. Shukla, A. Shalit, and P. Hamm, The Journal of Chem- ical Physics162, 034502 (2025)

work page 2025
[50]

Duchi, S

M. Duchi, S. Shukla, A. Shalit, and P. Hamm, The Jour- nal of Chemical Physics155, 174201 (2021)

work page 2021
[51]

K. M. Farrell, J. S. Ostrander, A. C. Jones, B. R. Yakami, S. S. Dicke, C. T. Middleton, P. Hamm, and M. T. Zanni, Opt. Express28, 33584 (2020)

work page 2020
[52]

Y. Feng, I. Vinogradov, and N.-H. Ge, Opt. Express27, 20323 (2019)

work page 2019
[53]

Savolainen, S

J. Savolainen, S. Ahmed, and P. Hamm, Proceedings of the National Academy of Sciences110, 20402 (2013)

work page 2013
[54]

P. E. Bl¨ ochl, Phys. Rev. B50, 17953 (1994)

work page 1994
[55]

Kresse and D

G. Kresse and D. Joubert, Phys. Rev. B59, 1758 (1999)

work page 1999
[56]

Kresse and J

G. Kresse and J. Furthm¨ uller, Computational Materials Science6, 15 (1996)

work page 1996
[57]

Kresse and J

G. Kresse and J. Furthm¨ uller, Phys. Rev. B54, 11169 (1996)

work page 1996
[58]

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

work page 1996
[59]

Baroni, S

S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Gian- nozzi, Rev. Mod. Phys.73, 515 (2001)

work page 2001
[60]

Giannozzi, O

P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavaz- zoni, D. Ceresoli, M. Cococcioni, et al., Journal of physics: Condensed matter29, 465901 (2017)

work page 2017
[61]

𝑓#−𝑓$ 0246𝑓!(THz) (b)𝑄%&,(! 𝑄&,*

P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, et al., Journal of physics: Condensed mat- ter21, 395502 (2009). 7 END MATTER Instrument Response Function—The absolute value of the instrument response function in the frequency do- main for the 2D-RTT pulse sequence is shown in ...

work page 2009
[62]

domain. The order parameter of the tilted structure can be constructed from the static rotational momenta of AlO6 octahedra as [45–47]: Φ(R) = (−1)i+j+k X l=1,6 ˆrl × ˆr′ l (15) wherei, j, kare integer numbers and the lattice vector R=ia X +ja Y +ka Z describes the B-site position in the supercell. Unit vector ˆrl (ˆr′ l) represents oxygen po- sitions mea...

work page
[63]

We adopted PBE exchange-correlation functional and used PAW pseudopotentials taken from the Quantum ESPRESSO pseudopotential data base

as implemented in Quantum ESPRESSO [60, 61]. We adopted PBE exchange-correlation functional and used PAW pseudopotentials taken from the Quantum ESPRESSO pseudopotential data base. Given the com- putational load of simulating the IR vibrational spectra of the 320-atom supercell, we performed lattice-dynamics calculations for the rhombohedral bulk phase of...

work page

[1] [1]

E. K. Salje, Annual Review of Materials Research42, 265 (2012)

work page 2012

[2] [2]

Nataf, M

G. Nataf, M. Guennou, J. Gregg, D. Meier, J. Hlinka, E. Salje, and J. Kreisel, Nature Reviews Physics2, 634 (2020)

work page 2020

[3] [3]

Arlt, Journal of Materials Science25, 2655 (1990)

G. Arlt, Journal of Materials Science25, 2655 (1990)

work page 1990

[4] [4]

Bueble, K

S. Bueble, K. Knorr, E. Brecht, and W. W. Schmahl, Surface science400, 345 (1998)

work page 1998

[5] [5]

E.-J. Guo, R. Desautels, D. Keavney, M. A. Roldan, B. J. Kirby, D. Lee, Z. Liao, T. Charlton, A. Herklotz, T. Zac Ward, et al., Science advances5, eaav5050 (2019)

work page 2019

[6] [6]

Kustov, I

S. Kustov, I. Liubimova, and E. Salje, Applied Physics Letters112(2018)

work page 2018

[7] [7]

Hayward, S

S. Hayward, S. Redfern, and E. Salje, Journal of Physics: Condensed Matter14, 10131 (2002)

work page 2002

[8] [8]

S. A. Hayward, F. D. Morrison, S. A. T. Redfern, E. K. H. Salje, J. F. Scott, K. S. Knight, S. Tarantino, A. M. Glazer, V. Shuvaeva, P. Daniel, M. Zhang, and M. A. Carpenter, Phys. Rev. B72, 054110 (2005)

work page 2005

[9] [9]

M. V. Abrashev, A. P. Litvinchuk, M. N. Iliev, R. L. Meng, V. N. Popov, V. G. Ivanov, R. A. Chakalov, and C. Thomsen, Phys. Rev. B59, 4146 (1999)

work page 1999

[10] [10]

Basini, V

M. Basini, V. Unikandanunni, F. Gabriele, M. Cross, A. M. Derrico, A. X. Gray, M. C. Hoffmann, F. Forte, 6 M. Cuoco, and S. Bonetti, Physical Review Letters136, 156902 (2026)

work page 2026

[11] [11]

Y. Liu, A. Frenkel, G. Garrett, J. Whitaker, S. Fahy, C. Uher, and R. Merlin, Physical review letters75, 334 (1995)

work page 1995

[12] [12]

Yokota, S

H. Yokota, S. Matsumoto, E. Salje, and Y. Uesu, Physical Review B98, 104105 (2018)

work page 2018

[13] [13]

Yokota and Y

H. Yokota and Y. Uesu, Journal of Applied Physics129 (2021)

work page 2021

[14] [14]

Yokota, C

H. Yokota, C. Haines, S. Matsumoto, N. Hasegawa, M. A. Carpenter, Y. Heo, A. Marin, E. Salje, and Y. Uesu, Physical Review B102, 104117 (2020)

work page 2020

[15] [15]

E. K. Salje, M. Alexe, S. Kustov, M. C. Weber, J. Schiemer, G. F. Nataf, and J. Kreisel, Scientific re- ports6, 27193 (2016)

work page 2016

[16] [16]

J. Suda, O. Kamishima, J. Kawamura, T. Hattori, and T. Sato, in Journal of Physics: Conference Series, Vol. 150 (IOP Publishing, 2009) p. 052249

work page 2009

[17] [17]

Hortensius, D

J. Hortensius, D. Afanasiev, A. Sasani, E. Bousquet, and A. Caviglia, npj Quantum Materials5, 95 (2020)

work page 2020

[18] [18]

M. J. Neugebauer, D. M. Juraschek, M. Savoini, P. En- geler, L. Boie, E. Abreu, N. A. Spaldin, and S. L. John- son, Phys. Rev. Res.3, 013126 (2021)

work page 2021

[19] [19]

Ciardi, A

G. Ciardi, A. Berger, P. Hamm, and A. Shalit, The Jour- nal of Physical Chemistry Letters10, 4463 (2019)

work page 2019

[20] [20]

Shalit, S

A. Shalit, S. J. Mousavi, and P. Hamm, The Journal of Physical Chemistry B125, 581 (2021)

work page 2021

[21] [21]

S. J. Mousavi, A. Berger, P. Hamm, and A. Shalit, The Journal of Chemical Physics156, 174501 (2022)

work page 2022

[22] [22]

A. A. Maradudin and R. F. Wallis, Phys. Rev. B2, 4294 (1970)

work page 1970

[23] [23]

R. F. Wallis and A. A. Maradudin, Phys. Rev. B3, 2063 (1971)

work page 2063

[24] [24]

L. B. Humphreys, Phys. Rev. B6, 3886 (1972)

work page 1972

[25] [25]

F¨ orst, C

M. F¨ orst, C. Manzoni, S. Kaiser, Y. Tomioka, Y. Tokura, R. Merlin, and A. Cavalleri, Nature Physics7, 854 (2011)

work page 2011

[26] [26]

Subedi, A

A. Subedi, A. Cavalleri, and A. Georges, Phys. Rev. B 89, 220301 (2014)

work page 2014

[27] [27]

von Hoegen, R

A. von Hoegen, R. Mankowsky, M. Fechner, M. F¨ orst, and A. Cavalleri, Nature555, 79 (2018)

work page 2018

[28] [28]

Khalsa, N

G. Khalsa, N. A. Benedek, and J. Moses, Phys. Rev. X 11, 021067 (2021)

work page 2021

[29] [29]

T. G. Blank, K. A. Grishunin, K. Zvezdin, N. Hai, J. Wu, S.-H. Su, J.-C. Huang, A. Zvezdin, and A. V. Kimel, Physical review letters131, 026902 (2023)

work page 2023

[30] [30]

Petretto, S

G. Petretto, S. Dwaraknath, H. PC Miranda, D. Win- ston, M. Giantomassi, M. J. Van Setten, X. Gonze, K. A. Persson, G. Hautier, and G.-M. Rignanese, Scientific data 5, 1 (2018)

work page 2018

[31] [31]

Miranda, Phonon website: Visualize phonon vibra- tional modes,https://henriquemiranda.github.io/ phononwebsite/(2025), accessed: 2025-11-19

H. Miranda, Phonon website: Visualize phonon vibra- tional modes,https://henriquemiranda.github.io/ phononwebsite/(2025), accessed: 2025-11-19

work page 2025

[32] [32]

A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder, and K. a. Persson, APL Materials1, 011002 (2013)

work page 2013

[33] [33]

Sidler and P

D. Sidler and P. Hamm, The Journal of Chemical Physics 150, 044202 (2019)

work page 2019

[34] [34]

Sidler and P

D. Sidler and P. Hamm, The Journal of Chemical Physics 153, 044502 (2020)

work page 2020

[35] [35]

J. F. Scott, Phys. Rev.183, 823 (1969)

work page 1969

[36] [36]

Lloyd-Hughes, S

J. Lloyd-Hughes, S. Jones, E. Castro-Camus, K. Doig, and J. MacManus-Driscoll, Optics letters39, 1121 (2014)

work page 2014

[37] [37]

S. J. Mousavi, M. F. Biggs, J. A. Johnson, P. Hamm, and A. Shalit, The Journal of Chemical Physics161, 134502 (2024)

work page 2024

[38] [38]

T. P. Martin and L. Genzel, Phys. Status Solidi B61, 493 (1974)

work page 1974

[39] [39]

D. M. Juraschek and S. F. Maehrlein, Phys. Rev. B97, 174302 (2018)

work page 2018

[40] [40]

Giorgianni, M

F. Giorgianni, M. Romani, P. Puphal, M. Isobe, L. Spitz, M. C. Guidi, C. Vicario, and M. Udina, Phys. Rev. B 111, 205138 (2025)

work page 2025

[41] [41]

Sertcan, S

B. Sertcan, S. J. Mousavi, M. Iannuzzi, and P. Hamm, The Journal of Chemical Physics158, 014203 (2023)

work page 2023

[42] [42]

Fiore, N

J. Fiore, N. Sellati, M. Udina, and L. Benfatto, arXiv preprint arXiv:2509.25060 (2025)

work page arXiv 2025

[43] [43]

Proietto, A

F. Proietto, A. Milloch, P. Franceschini, M. Azarm, N. Sellati, R. Mishra, P. C. Moen, S. Smit, M. Bluschke, M. Greven, et al., arXiv preprint arXiv:2603.29713 (2026)

work page arXiv 2026

[44] [44]

A. M. Glazer, Acta Crystallographica Section A31, 756 (1975)

work page 1975

[45] [45]

Gopalan and D

V. Gopalan and D. B. Litvin, Nature Materials10, 376 (2011)

work page 2011

[46] [46]

Bellaiche and J

L. Bellaiche and J. ´I˜ niguez, Phys. Rev. B88, 014104 (2013)

work page 2013

[47] [47]

Barone, D

P. Barone, D. Di Sante, and S. Picozzi, Phys. Rev. B89, 144104 (2014)

work page 2014

[48] [48]

Salje, Ferroelectrics104, 111 (1990)

E. Salje, Ferroelectrics104, 111 (1990)

work page 1990

[49] [49]

Shukla, A

S. Shukla, A. Shalit, and P. Hamm, The Journal of Chem- ical Physics162, 034502 (2025)

work page 2025

[50] [50]

Duchi, S

M. Duchi, S. Shukla, A. Shalit, and P. Hamm, The Jour- nal of Chemical Physics155, 174201 (2021)

work page 2021

[51] [51]

K. M. Farrell, J. S. Ostrander, A. C. Jones, B. R. Yakami, S. S. Dicke, C. T. Middleton, P. Hamm, and M. T. Zanni, Opt. Express28, 33584 (2020)

work page 2020

[52] [52]

Y. Feng, I. Vinogradov, and N.-H. Ge, Opt. Express27, 20323 (2019)

work page 2019

[53] [53]

Savolainen, S

J. Savolainen, S. Ahmed, and P. Hamm, Proceedings of the National Academy of Sciences110, 20402 (2013)

work page 2013

[54] [54]

P. E. Bl¨ ochl, Phys. Rev. B50, 17953 (1994)

work page 1994

[55] [55]

Kresse and D

G. Kresse and D. Joubert, Phys. Rev. B59, 1758 (1999)

work page 1999

[56] [56]

Kresse and J

G. Kresse and J. Furthm¨ uller, Computational Materials Science6, 15 (1996)

work page 1996

[57] [57]

Kresse and J

G. Kresse and J. Furthm¨ uller, Phys. Rev. B54, 11169 (1996)

work page 1996

[58] [58]

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

work page 1996

[59] [59]

Baroni, S

S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Gian- nozzi, Rev. Mod. Phys.73, 515 (2001)

work page 2001

[60] [60]

Giannozzi, O

P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavaz- zoni, D. Ceresoli, M. Cococcioni, et al., Journal of physics: Condensed matter29, 465901 (2017)

work page 2017

[61] [61]

𝑓#−𝑓$ 0246𝑓!(THz) (b)𝑄%&,(! 𝑄&,*

P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, et al., Journal of physics: Condensed mat- ter21, 395502 (2009). 7 END MATTER Instrument Response Function—The absolute value of the instrument response function in the frequency do- main for the 2D-RTT pulse sequence is shown in ...

work page 2009

[62] [62]

domain. The order parameter of the tilted structure can be constructed from the static rotational momenta of AlO6 octahedra as [45–47]: Φ(R) = (−1)i+j+k X l=1,6 ˆrl × ˆr′ l (15) wherei, j, kare integer numbers and the lattice vector R=ia X +ja Y +ka Z describes the B-site position in the supercell. Unit vector ˆrl (ˆr′ l) represents oxygen po- sitions mea...

work page

[63] [63]

We adopted PBE exchange-correlation functional and used PAW pseudopotentials taken from the Quantum ESPRESSO pseudopotential data base

as implemented in Quantum ESPRESSO [60, 61]. We adopted PBE exchange-correlation functional and used PAW pseudopotentials taken from the Quantum ESPRESSO pseudopotential data base. Given the com- putational load of simulating the IR vibrational spectra of the 320-atom supercell, we performed lattice-dynamics calculations for the rhombohedral bulk phase of...

work page