In the rocksalt IV–VI chalcogenides AQ (A = Ge, Sn, Pb; Q = S, Se, Te), Waghmare, Spaldin, Kandpal, and Seshadri (Phys. Rev. B 67, 125111 (2003)) showed that the longitudinal–transverse optical (LO–TO) phonon splitting at the zone center, computed in the undistorted high-symmetry phase, tracks the observed tendency of the cation to off-center better than other single quantities (Born effective charges, dielectric constants, electron localization lengths), because it reflects both the polar structural instability and the screening associated with incipient metallicity. The authors themselves stated that this indicator was established only within the nine studied compounds, at theoretical lattice constants, without spin–orbit coupling in all electronic-structure stages, and at zero temperature. The open problem is to determine whether the LO–TO splitting is a generally valid, quantitative predictor of structural (off-centering) versus metallic instability in other chemistries, other crystal structures, and at finite temperature — including the volume and relativistic corrections needed for quantitative prediction of transition temperatures. Later work has deepened the finite-temperature physics of the same family and extended lone-pair concepts to other chemistries using other descriptors, but no systematic validation or refutation of the LO–TO splitting indicator beyond the original nine compounds has been reported in the literature.

Background

The compounds AQ with a divalent group-14 cation A = Ge, Sn, or Pb and a chalcogen Q = S, Se, or Te (the IV–VI chalcogenides) crystallize near the cubic rocksalt structure. Some of them (GeTe, and SnTe at low temperature) spontaneously distort to a rhombohedral, polar (ferroelectric) phase in which the cation shifts off the center of its octahedral coordination site; others (GeS, GeSe, SnS, SnSe) adopt layered low-symmetry structures; PbS, PbSe, and PbTe remain cubic. The cation carries a nominally inert valence s² electron pair (the 'lone pair'), and its tendency to localize into a directed lobe — the stereochemical activity of the lone pair — is the chemical origin of this off-centricity. The traditional rationalizations (on-site cation s–p mixing after Orgel, and the VSEPR rules of valence-shell electron-pair repulsion) proved incomplete: the anion p states are essential to forming the lone-pair lobe.

The relevant first-principles quantities come from density-functional theory (DFT) linear-response (density-functional perturbation theory), which computes phonon frequencies and dielectric responses without supercells. The Born effective charge Z∗Z^{*} is the change of macroscopic polarization per unit displacement of an atom; in these covalent, highly polarizable compounds it is anomalously large (values far above the formal ionic charge +2), a signature of dynamic covalency rather than of a definite structural instability. The high-frequency (optical) dielectric constant ε∞\varepsilon_{\infty} measures purely electronic screening. Because long-range electrostatic dipole–dipole forces act only on longitudinal displacements, the longitudinal-optical (LO) and transverse-optical (TO) phonons at the zone center split by an amount that scales as ΔωLO−TO∝Z∗2/(ε∞M)\Delta\omega_{\mathrm{LO-TO}} \propto Z^{*2}/(\varepsilon_{\infty}M) (M the reduced mass) — a single number that combines the strength of the polar coupling with the electronic screening that suppresses it. A ferroelectric (polar) structural instability appears as a TO phonon that softens to zero frequency; the competing instability is electronic delocalization (incipient metallicity), which grows as the gap narrows and the electronic localization measure l2l^{2} increases along the series S → Se → Te.

Within this framework, Waghmare, Spaldin, Kandpal, and Seshadri (Phys. Rev. B 67, 125111 (2003)) surveyed the nine rocksalt-structure-related compounds with plane-wave pseudopotential linear response and with linear muffin-tin orbital (LMTO, an all-electron band-structure method) orbital-resolved analysis. They found that the two instabilities compete: the polar one (soft TO mode, lone-pair localization) and the metallic one (delocalization, large ε∞\varepsilon_{\infty}). Individually, Z∗Z^{*} is anomalously large even in compounds that stay cubic, ε∞\varepsilon_{\infty} and l2l^{2} grow with metallicity rather than with distortion tendency, and the screened (longitudinal) effective charges capture only part of the trend. The computed LO–TO splitting, however, decreases monotonically from GeS to PbTe exactly in the order of decreasing observed off-centering tendency, and they proposed ΔωLO−TO\Delta\omega_{\mathrm{LO-TO}} computed in the high-symmetry phase as a single ab initio indicator of proximity to the competing structural and metallic instabilities.

Three limitations of that study define the open question. First, harmonic phonon frequencies — especially the soft TO modes — are highly sensitive to the cell volume, so calculations at the theoretical (local-density-approximation, a standard DFT exchange-correlation choice) lattice constants cannot be expected to give quantitatively accurate phonon frequencies or ferroelectric transition temperatures; the study was a zero-temperature one. Second, the LMTO stages omitted explicit spin–orbit coupling (SOC, the interaction between the electron spin and the orbital motion that is large for heavy elements), which for the heavier cations and anions can shift band positions and width and hence the delicate balance among the symmetry-adapted cation s state (S+S^{+}), anion p state (P−P^{-}), and cation p state (P+P^{+}) that controls the lone-pair formation and the approach to metallicity. Third, the indicator was calibrated on the same nine compounds whose experimental structures were already known, so its predictive power on genuinely independent materials — other chemistries, other crystal structures, or at finite temperature — was untested. These three items jointly ask whether the indicator can be turned from a same-family correlation into a validated quantitative predictor for materials design.

Problem Statement

Determine whether the LO–TO phonon splitting ΔωLO−TO\Delta\omega_{\mathrm{LO-TO}} computed from first principles in the undistorted high-symmetry phase is a generally valid and quantitative indicator that distinguishes compounds prone to a polar structural instability (cation off-centricity, lone-pair stereochemical activity) from compounds dominated by incipient metallicity — beyond the nine IV–VI rocksalt chalcogenides (A = Ge, Sn, Pb; Q = S, Se, Te) for which it was originally calibrated. Concretely: (1) test the indicator on independent compound sets spanning other chemistries and other centrosymmetric parent crystal structures, stating its predictive accuracy, its failure modes (for example where octahedral tilting — collective rotation of the corner-sharing octahedra — a Peierls-like distortion (a unit-cell-doubling, gap-opening lattice distortion driven by electron–phonon coupling), so-called metavalent bonding (a proposed electron-sharing regime intermediate between covalent and metallic bonding, invoked for phase-change chalcogenides), or other competing instabilities dominate, or in low-dimensional structures where the longitudinal–transverse splitting is not defined by the same three-dimensional electrostatics), and whether it performs better than, or is redundant with, computing a full phonon spectrum; (2) extend the indicator to finite temperature — connecting it to anharmonic phonon theory so that transition temperatures, not only zero-temperature tendencies, can be predicted quantitatively, with explicit treatment of the demonstrated sensitivity of harmonic frequencies to the exchange-correlation functional and to cell volume (experimental versus theoretical lattice constants); and (3) establish quantitatively how spin–orbit coupling modifies the balance among the cation s (S+S^{+}), anion p (P−P^{-}), and cation p (P+P^{+}) states that controls both the lone-pair formation and the approach to metallicity in heavy-element compounds, and hence how it shifts the indicator. The overarching objective is a validated, calibrated, single-quantity first-principles predictor usable for materials design — or a demonstrated, evidence-backed statement of the class of materials for which it fails.

References

  1. U. V. Waghmare; N. A. Spaldin; H. C. Kandpal; Ram Seshadri, 'First-principles indicators of metallicity and cation off-centricity in the IV-VI rocksalt chalcogenides of divalent Ge, Sn, and Pb', Physical Review B 67, 125111 (2003), DOI 10.1103/PhysRevB.67.125111, https://doi.org/10.1103/PhysRevB.67.125111
  2. Emil S. Božin; Christos D. Malliakas; Petros Souvatzis; Thomas Proffen; Nicola A. Spaldin; Mercouri G. Kanatzidis; Simon J. L. Billinge, 'Entropically Stabilized Local Dipole Formation in Lead Chalcogenides', Science 330 (6011), 1660–1663 (2010), DOI 10.1126/science.1192759, https://doi.org/10.1126/science.1192759
  3. Bart J. Kooi; Beatriz Noheda, 'Ferroelectric chalcogenides—materials at the edge', Science 353 (6296), 221–222 (2016), DOI 10.1126/science.aaf9081, https://doi.org/10.1126/science.aaf9081
  4. Leena Aggarwal; Ananya Banik; Shashwat Anand; Umesh V. Waghmare; Kanishka Biswas; Goutam Sheet, 'Local ferroelectricity in thermoelectric SnTe above room temperature driven by competing phonon instabilities and soft resonant bonding', Journal of Materiomics 2 (2), 196–202 (2016), DOI 10.1016/j.jmat.2016.04.001, https://doi.org/10.1016/j.jmat.2016.04.001
  5. M. P. Jiang; M. Trigo; I. Savić; S. Fahy; É. D. Murray; C. Bray; J. Clark; T. Henighan; M. Kozina; M. Chollet; J. M. Glownia; M. C. Hoffmann; D. Zhu; O. Delaire; A. F. May; B. C. Sales; A. M. Lindenberg; P. Zalden; T. Sato; R. Merlin; D. A. Reis, 'The origin of incipient ferroelectricity in lead telluride', Nature Communications 7, 12291 (2016), DOI 10.1038/ncomms12291, https://doi.org/10.1038/ncomms12291
  6. Ronan M. Murphy; Éamonn D. Murray; Stephen Fahy; Ivana Savić, 'Broadband phonon scattering in PbTe-based materials driven near ferroelectric phase transition by strain or alloying', Physical Review B 93, 104304 (2016), DOI 10.1103/PhysRevB.93.104304, https://doi.org/10.1103/PhysRevB.93.104304
  7. Guilherme A. S. Ribeiro; Lorenzo Paulatto; Raffaello Bianco; Ion Errea; Francesco Mauri; Matteo Calandra, 'Strong anharmonicity in the phonon spectra of PbTe and SnTe from first principles', Physical Review B 97, 014306 (2018), DOI 10.1103/PhysRevB.97.014306, https://doi.org/10.1103/PhysRevB.97.014306
  8. Wenhui Wan; Chang Liu; Wende Xiao; Yugui Yao, 'Promising ferroelectricity in 2D group IV tellurides: a first-principles study', Applied Physics Letters 111, 132904 (2017), DOI 10.1063/1.4996171, https://doi.org/10.1063/1.4996171
  9. Y. Okamura; H. Handa; R. Yoshimi; A. Tsukazaki; K. S. Takahashi; M. Kawasaki; Y. Tokura; Y. Takahashi, 'Terahertz lattice and charge dynamics in ferroelectric semiconductor SnxPb1−xTe', npj Quantum Materials 7, 91 (2022), DOI 10.1038/s41535-022-00501-2, https://doi.org/10.1038/s41535-022-00501-2
  10. Ananya Banik; Tanmoy Ghosh; Raagya Arora; Moinak Dutta; Juhi Pandey; Somnath Acharya; Ajay Soni; Umesh V. Waghmare; Kanishka Biswas, 'Engineering ferroelectric instability to achieve ultralow thermal conductivity and high thermoelectric performance in Sn1−xGexTe', Energy & Environmental Science 12 (2), 589–595 (2019), DOI 10.1039/C8EE03162B, https://doi.org/10.1039/C8EE03162B
  11. Geneva Laurita; Douglas H. Fabini; Constantinos C. Stoumpos; Mercouri G. Kanatzidis; Ram Seshadri, 'Chemical tuning of dynamic cation off-centering in the cubic phases of hybrid tin and lead halide perovskites', Chemical Science 8 (8), 5628–5635 (2017), DOI 10.1039/c7sc01429e, https://doi.org/10.1039/c7sc01429e
  12. Kasper Tolborg; Carlo Gatti; Bo B. Iversen, 'Expression and interactions of stereochemically active lone pairs and their relation to structural distortions and thermal conductivity', IUCrJ 7 (3), 480–489 (2020), DOI 10.1107/S2052252520003619, https://doi.org/10.1107/S2052252520003619
  13. Dejan Zagorac; Klaus Doll; J. Christian Schön; Martin Jansen, 'Sterically Active Electron Pairs in Lead Sulfide? An Investigation of the Electronic and Vibrational Properties of PbS in the Transition Region Between the Rock Salt and the α-GeTe-Type Modifications', Chemistry – A European Journal 18 (35), 10929–10936 (2012), DOI 10.1002/chem.201200180, https://doi.org/10.1002/chem.201200180
  14. Shuva Biswas; Sakshi Verma; Subarna Das; Debattam Sarkar; Adrija Ghosh; Diptikanta Swain; Tapas Kumar Maji; Umesh V. Waghmare; Kanishka Biswas, 'Ferroelectric Instability Driven Local Atomic Off-Centering and Low Thermal Conductivity in a Topological Crystalline Insulator Sn3Bi0.7Te4', Chemistry of Materials 38 (15), 8000–8011 (2026), DOI 10.1021/acs.chemmater.6c01598, https://doi.org/10.1021/acs.chemmater.6c01598
  15. Hongyao Xie; Li-Dong Zhao, 'Origin of off-centering effect and the influence on heat transport in thermoelectrics', Materials Futures 3, 013501 (2024), DOI 10.1088/2752-5724/ad1ac0, https://doi.org/10.1088/2752-5724/ad1ac0