pith. sign in

arxiv: 2510.15703 · v2 · submitted 2025-10-17 · ❄️ cond-mat.mtrl-sci

Lattice excitations with finite polarization and magnetization

Mike Pols , Carl P. Romao , Dominik M. Juraschek This is my paper

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

classification ❄️ cond-mat.mtrl-sci
keywords multiferronsdynamical multiferroicityLiNbO3ferroelectric polarizationmagnetizationlattice excitationsmultipolonsmultiferroics
0
0 comments X p. Extension

The pith

Multiferrons are lattice excitations that tilt polarization and generate magnetization through dynamical multiferroicity.

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

The paper introduces multiferrons as quasiparticles that act as elementary excitations carrying both electric polarization and magnetization. These excitations arise when lattice vibrations produce a tilt and elliptical precession of the polarization while simultaneously creating a net magnetization. First-principles calculations for LiNbO3 show that the multiferron polarization points perpendicular to the material's equilibrium ferroelectric polarization and the magnetization points parallel to it. Multiferrons also possess net quadrupole and octupole moments that the authors call multipolons. A reader would care because this mechanism links electric and magnetic responses at the level of lattice vibrations without external fields.

Core claim

Multiferrons are elementary excitations of the ferroelectric order that possess both electric and magnetic character. They produce a tilt and elliptical precession of the polarization and create a magnetization through dynamical multiferroicity. First-principles calculations in LiNbO3 show the electric polarization of multiferrons to be perpendicular to the equilibrium ferroelectric polarization while the magnetization is parallel to it. Multiferrons further carry net electric and magnetic quadrupole and octupole moments termed multipolons that could couple to internal multipolar degrees of freedom or external probes such as neutrons.

What carries the argument

Multiferrons, the elementary excitations with both electric and magnetic character produced by dynamical multiferroicity in the lattice.

Load-bearing premise

Dynamical multiferroicity produces a net magnetization from lattice excitations in LiNbO3 without additional fitting parameters or external fields.

What would settle it

A first-principles calculation or measurement showing zero net magnetization parallel to the ferroelectric polarization after coherent excitation of the relevant lattice modes in LiNbO3 would falsify the claim.

Figures

Figures reproduced from arXiv: 2510.15703 by Carl P. Romao, Dominik M. Juraschek, Mike Pols.

Figure 1
Figure 1. Figure 1: FIG. 1. Ferrons in LiNbO [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Magnetization induced by circular and elliptical exc [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Electric polarization and magnetization dynamics ind [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Time-averaged quadrupole tensors of radial polar [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

Ferrons are a type of quasiparticle corresponding to elementary excitations of the ferroelectric order. Analogously to how magnons modulate and transport magnetization, ferrons modulate and transport electric polarization. Here, we introduce multiferrons as elementary excitations with both electric and magnetic character. Multiferrons lead to a tilt and elliptical precession of the polarization and at the same time create a magnetization through the mechanism of dynamical multiferroicity. Using first-principles calculations for LiNbO$_3$, we show that the electric polarization of multiferrons is perpendicular to the equilibrium ferroelectric polarization, whereas the magnetization is parallel to it. Our calculations further demonstrate that multiferrons carry net electric and magnetic quadrupole and octupole moments, which we term multipolons. These multipolons could couple to internal multipolar degrees of freedom, for example in altermagnets, or to external probes such as neutrons, leading to potentially experimentally observable phenomena following coherent or thermal excitation of multiferrons.

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

2 major / 2 minor

Summary. The manuscript introduces multiferrons as lattice excitations possessing both ferroelectric and magnetic character. Using first-principles calculations on LiNbO3, it reports that multiferrons produce an electric polarization perpendicular to the equilibrium ferroelectric polarization and a magnetization parallel to it, arising from tilt and elliptical precession via dynamical multiferroicity. The excitations are further shown to carry net electric and magnetic quadrupole and octupole moments, termed multipolons, with potential coupling to internal or external probes.

Significance. If the computational results are robust, the work would provide a concrete microscopic mechanism linking lattice dynamics to induced magnetization in a non-magnetic ferroelectric, extending the concept of dynamical multiferroicity to quasiparticle excitations and suggesting new multipolar observables.

major comments (2)
  1. [Methods] Methods section (computational details for LiNbO3): the procedure used to extract a time-averaged net magnetization from the multiferron dynamics is not specified. It is unclear whether the reported parallel magnetization is obtained by explicit time integration over the precession trajectory, by an effective static model, or by a single-snapshot evaluation; this detail is load-bearing for the central claim that dynamical multiferroicity produces a finite M without external bias.
  2. [Results] Results section on polarization and magnetization directions: the statement that the multiferron polarization is perpendicular and magnetization parallel to the equilibrium ferroelectric axis relies on the output of the first-principles run, yet no convergence tests, k-point sampling, or error estimates on the computed moments are provided. Without these, it is difficult to assess whether the reported directions are numerically stable or sensitive to the chosen supercell or exchange-correlation functional.
minor comments (2)
  1. [Introduction] The introduction of the terms 'multiferron' and 'multipolon' would benefit from a brief comparison to existing literature on hybrid magnon-phonon or electromagnon modes to clarify novelty.
  2. [Figures] Figure captions describing the precession and multipole moments should explicitly state the coordinate system relative to the equilibrium polarization axis for clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major comment below and have revised the manuscript to improve clarity and provide additional supporting information.

read point-by-point responses

  1. Referee: [Methods] Methods section (computational details for LiNbO3): the procedure used to extract a time-averaged net magnetization from the multiferron dynamics is not specified. It is unclear whether the reported parallel magnetization is obtained by explicit time integration over the precession trajectory, by an effective static model, or by a single-snapshot evaluation; this detail is load-bearing for the central claim that dynamical multiferroicity produces a finite M without external bias.

    Authors: We agree that the extraction procedure requires explicit clarification. The reported magnetization is obtained by explicit time integration of the instantaneous magnetization generated via dynamical multiferroicity (M(t) proportional to the cross product of the polarization vector and its time derivative) over multiple periods of the elliptical precession, followed by a time average. We have added a dedicated paragraph in the revised Methods section describing the integration scheme, the numerical time step used, and the averaging formula. revision: yes

  2. Referee: [Results] Results section on polarization and magnetization directions: the statement that the multiferron polarization is perpendicular and magnetization parallel to the equilibrium ferroelectric axis relies on the output of the first-principles run, yet no convergence tests, k-point sampling, or error estimates on the computed moments are provided. Without these, it is difficult to assess whether the reported directions are numerically stable or sensitive to the chosen supercell or exchange-correlation functional.

    Authors: We acknowledge that convergence information was not included in the original submission. The calculations were performed with a converged 6×6×6 k-point mesh on a 2×2×2 supercell using the PBE functional; additional tests with denser meshes and larger supercells confirm that the perpendicular and parallel directions remain stable within numerical noise. We have added a new subsection to the Supplementary Information that reports these convergence tests together with estimated uncertainties on the induced moments obtained from the standard deviation over the dynamical trajectory. revision: yes

Circularity Check

0 steps flagged

No circularity: first-principles outputs are independent of inputs

full rationale

The paper's central results—the perpendicular electric polarization and parallel magnetization of multiferrons in LiNbO3—are presented as direct outputs of first-principles calculations rather than definitions, fits, or self-citations. The abstract and derivation invoke dynamical multiferroicity as the underlying mechanism but do not reduce the reported directions or multipole moments to a fitted parameter or prior self-citation by construction. No equation is shown to equal its own input, and the computational methodology supplies independent, falsifiable content against external benchmarks. This is the expected non-finding for a DFT-based study.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 2 invented entities

The work rests on standard density-functional theory for lattice dynamics plus the dynamical-multiferroicity mechanism; no new free parameters are explicitly introduced in the abstract, but the computational setup implicitly contains typical DFT choices such as exchange-correlation functional and k-point sampling.

axioms (1)
  • domain assumption Dynamical multiferroicity generates net magnetization from ionic displacements in a ferroelectric lattice
    Invoked to link the lattice excitations to finite magnetization without additional external fields.
invented entities (2)
  • multiferron no independent evidence
    purpose: Quasiparticle carrying both polarization and magnetization
    Newly defined excitation combining ferroelectric and magnetic character
  • multipolon no independent evidence
    purpose: Higher-order electric and magnetic moments carried by multiferrons
    Term introduced for quadrupole and octupole moments

pith-pipeline@v0.9.0 · 5703 in / 1382 out tokens · 35545 ms · 2026-05-18T06:20:28.630776+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

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

  1. [1]

    16 THz [ 15, 16] (SM Note 1). We therefore express the polarization in terms of the phonon amplitude Q as Pph = ZQ/V c, (2) arXiv:2510.15703v1 [cond-mat.mtrl-sci] 17 Oct 2025 2 P x P y Ferron In-plane ferron P z P x P y (a) [111] P 0 P ph (b) (c) [111] P 0 P ph Multiferron [111] P 0 P ph V P P P P M P P M FIG. 1. Ferrons in LiNbO 3. (a) Excitation of A1 m...

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  2. [2]

    ( 1) yields an equivalent double-well potential in terms of Q, which provides a natural starting point for analyzing lat- tice excitations in LiNbO 3

    into Eq. ( 1) yields an equivalent double-well potential in terms of Q, which provides a natural starting point for analyzing lat- tice excitations in LiNbO 3. In the following, we will de- rive the net electric polarization and magnetization gen- erated by multiferrons. Close to the minima at ± P0, the potential energy of a mode can be described using a ...

    work page

  3. [3]

    1a), whereas excitations of the E modes lead to an increase and tilting of the total polarization (Fig

    Excitations of the A1 modes decrease the magnitude of the total polarization (Fig. 1a), whereas excitations of the E modes lead to an increase and tilting of the total polarization (Fig. 1b), which could be referred to as Higgs- and Goldstone-like ferrons, respectively [ 5]. B. Magnetization When the degenerate E modes are excited circularly, the resultin...

    work page

  4. [4]

    Because the out-of-plane magnetization and radial magnetization both scale with the magnitude of the ionic charge current, the contri- bution of each atom to the radial magnetization can be expressed as |Mrad,n |= | Mn| |P0| |Pph,n | (7) where |Pph,n | and |Mn| denote the atomic contribu- tions to the phonon-induced polarization and magnetiza- 3 M rad P 0...

    work page

  5. [5]

    The temporal profile and polarization of the driving field govern the am- plitude and phase of the coherent phonon response

    describes the external driving force from a laser pulse, where Zα is the mode effective charge that deter- mines the coupling to the electric field E. The temporal profile and polarization of the driving field govern the am- plitude and phase of the coherent phonon response. The functional form used to model the laser pulse is given in SM Note 4. In Fig. 3, w...

    work page

  6. [6]

    A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, Magnon spintronics, Nat. Phys 11, 453 (2015)

    work page 2015

  7. [7]

    S. M. Rezende, A. Azevedo, and R. L. Rodr ´ ıguez-Su´ arez, Introduction to antiferromagnetic magnons, J. Appl. Phys. 126, 151101 (2019)

    work page 2019

  8. [8]

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

    work page 2022

  9. [9]

    G. E. W. Bauer, P. Tang, R. Iguchi, and K.-i. Uchida, Magnonics vs. Ferronics, J. Magn. Magn. Mater. 541, 168468 (2022)

    work page 2022

  10. [10]

    Bauer, P

    G. Bauer, P. Tang, R. Iguchi, J. Xiao, K. Shen, Z. Zhong, T. Yu, S. Rezende, J. Heremans, and K. Uchida, Polar- ization transport in ferroelectrics, Phys. Rev. Appl. 20, 050501 (2023)

    work page 2023

  11. [11]

    X.-H. Zhou, C. Cai, P. Tang, R. L. Rodr ´ ıguez-Su´ arez, S. M. Rezende, G. E. W. Bauer, and T. Yu, Surface Fer- ron Excitations in Ferroelectrics and Their Directional Routing, Chinese Phys. Lett. 40, 087103 (2023)

    work page 2023

  12. [12]

    Tang and G

    P. Tang and G. E. W. Bauer, Electric analog of magnons in order-disorder ferroelectrics, Phys. Rev. B 109, L060301 (2024)

    work page 2024

  13. [13]

    K. Shen, P. Tang, X. Chen, Y. Gao, Y. Fan, Z. Guo, Y. Wei, H. Jiang, X. Zhang, M. Wang, P. He, W. Shi, J. Han, Y. Wu, J. Shen, Q. Liu, G. E. W. Bauer, and M. Liu, Observation of ferron transport in ferroelectrics (2025), 2505.24419 [physics]

    work page arXiv 2025

  14. [14]

    J. Choe, T. Handa, C.-Y. Huang, A. K. Liston, J. Cox, J. Stensberg, Y. Hong, D. G. Chica, D. Xu, F. Tay, 6 S. Husremovic, V. d. S. L. Avelar, E. A. Arsenault, Z. Zhang, J. McIver, D. N. Basov, M. Delor, X. Roy, and X.-Y. Zhu, Observation of Coherent Ferrons (2025), 2505.22559 [cond-mat]

    work page arXiv 2025

  15. [15]

    A. N. Morozovska, E. A. Eliseev, O. V. Bereznikov, M. Y. Yelisieiev, G.-D. Zhao, Y. Zhu, V. Gopalan, L.-Q. Chen, J.-M. Hu, and Y. M. Vysochanskii, Flexocoupling- induced phonons and ferrons in van der Waals ferro- electrics, Phys. Rev. B 112, 014110 (2025)

    work page 2025

  16. [16]

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

    work page 2025

  17. [17]

    G. E. W. Bauer, P. Tang, R. Iguchi, J. Xiao, K. Shen, Z. Zhong, T. Yu, S. M. Rezende, J. P. Heremans, and K. Uchida, A Perspective on Ferrons (2023), 2302.12985 [cond-mat]

    work page arXiv 2023

  18. [18]

    S. C. Abrahams, H. J. Levinstein, and J. M. Reddy, Fer- roelectric lithium niobate. 5. Polycrystal X-ray diffrac- tion study between 24 ◦ and 1200 ◦ C, J. Phys. Chem. Solids 27, 1019 (1966)

    work page 1966

  19. [19]

    Dhar and A

    A. Dhar and A. Mansingh, Optical properties of reduced lithium niobate single crystals, J. Appl. Phys. 68, 5804 (1990)

    work page 1990

  20. [20]

    Mankowsky, A

    R. Mankowsky, A. von Hoegen, M. F¨ orst, and A. Caval- leri, Ultrafast Reversal of the Ferroelectric Polarization, Phys. Rev. Lett. 118, 197601 (2017)

    work page 2017

  21. [21]

    Henstridge, M

    M. Henstridge, M. F¨ orst, E. Rowe, M. Fechner, and A. Cavalleri, Nonlocal nonlinear phononics, Nat. Phys. 18, 457 (2022)

    work page 2022

  22. [22]

    D. M. Juraschek, M. Fechner, A. V. Balatsky, and N. A. Spaldin, Dynamical multiferroicity, Phys. Rev. Mater. 1, 014401 (2017)

    work page 2017

  23. [23]

    D. M. Juraschek and N. A. Spaldin, Orbital magnetic mo- ments of phonons, Phys. Rev. Mater. 3, 064405 (2019)

    work page 2019

  24. [24]

    Cheng, T

    B. Cheng, T. Schumann, Y. Wang, X. Zhang, D. Bar- balas, S. Stemmer, and N. P. Armitage, A Large Effective Phonon Magnetic Moment in a Dirac Semimetal, Nano Lett. 20, 5991 (2020)

    work page 2020

  25. [25]

    R. M. Geilhufe, V. Juriˇ ci´ c, S. Bonetti, J.-X. Zhu, and A. V. Balatsky, Dynamically Induced Magnetism in KTaO3, Phys. Rev. Res. 3, L022011 (2021)

    work page 2021

  26. [26]

    C. Xiao, Y. Ren, and B. Xiong, Adiabatically induced orbital magnetization, Phys. Rev. B 103, 115432 (2021)

    work page 2021

  27. [27]

    Y. Ren, C. Xiao, D. Saparov, and Q. Niu, Phonon Mag- netic Moment from Electronic Topological Magnetiza- tion, Phys. Rev. Lett. 127, 186403 (2021)

    work page 2021

  28. [28]

    Baydin, F

    A. Baydin, F. G. G. Hernandez, M. Rodriguez-Vega, A. K. Okazaki, F. Tay, G. T. Noe, I. Katayama, J. Takeda, H. Nojiri, P. H. O. Rappl, E. Abramof, G. A. Fiete, and J. Kono, Magnetic Control of Soft Chiral Phonons in PbTe, Phys. Rev. Lett. 128, 075901 (2022)

    work page 2022

  29. [29]

    F. G. G. Hernandez, A. Baydin, S. Chaudhary, F. Tay, I. Katayama, J. Takeda, H. Nojiri, A. K. Okazaki, P. H. O. Rappl, E. Abramof, M. Rodriguez-Vega, G. A. Fiete, and J. Kono, Observation of interplay between phonon chirality and electronic band topology, Science Advances 9, eadj4074 (2023)

    work page 2023

  30. [30]

    R. M. Geilhufe and W. Hergert, Electron Magnetic Mo- ment of Transient Chiral Phonons in KTaO 3, Phys. Rev. B 107, L020406 (2023)

    work page 2023

  31. [31]

    Zhang, Y

    X.-W. Zhang, Y. Ren, C. Wang, T. Cao, and D. Xiao, Gate-Tunable Phonon Magnetic Moment in Bilayer Graphene, Phys. Rev. Lett. 130, 226302 (2023)

    work page 2023

  32. [32]

    Shabala and R

    N. Shabala and R. M. Geilhufe, Phonon Inverse Faraday Effect from Electron-Phonon Coupling, Phys. Rev. Lett. 133, 266702 (2024)

    work page 2024

  33. [33]

    Klebl, A

    L. Klebl, A. Schobert, M. Eckstein, G. Sangiovanni, A. V. Balatsky, and T. O. Wehling, Ultrafast Pseudomagnetic Fields from Electron-Nuclear Quantum Geometry, Phys. Rev. Lett. 134, 016705 (2025)

    work page 2025

  34. [34]

    Mustafa, C

    H. Mustafa, C. Nnokwe, G. Ye, M. Fang, S. Chaudhary, J.-A. Yan, K. Wu, C. J. Cunningham, C. M. Hemesath, A. J. Stollenwerk, P. M. Shand, E.-H. Yang, G. A. Fiete, R. He, and W. Jin, Origin of Large Effective Phonon Magnetic Moments in Monolayer MoS 2, ACS Nano 19, 11241 (2025)

    work page 2025

  35. [35]

    Chen, X.-W

    W. Chen, X.-W. Zhang, Y. Su, T. Cao, D. Xiao, and S.- Z. Lin, Gauge theory of giant phonon magnetic moment in doped Dirac semimetals, Phys. Rev. B 111, 035126 (2025)

    work page 2025

  36. [36]

    Paiva, M

    C. Paiva, M. Fechner, and D. M. Juraschek, Dynamically Induced Multiferroic Polarization, Phys. Rev. Lett. 135, 066702 (2025)

    work page 2025

  37. [37]

    F. Wang, X. Liu, H. Sun, H. Wang, S. Murakami, L. Zhang, H. Zhang, and D. Xing, Ab Initio The- ory of Large Phonon Magnetic Moments Induced by Electron-Phonon Coupling in Magnetic Materials (2025), 2503.10160 [cond-mat]

    work page arXiv 2025

  38. [38]

    J. Luo, T. Lin, J. Zhang, X. Chen, E. R. Blackert, R. Xu, B. I. Yakobson, and H. Zhu, Large effective magnetic fields from chiral phonons in rare-earth halides, Science 382, 698 (2023)

    work page 2023

  39. [39]

    Basini, M

    M. Basini, M. Pancaldi, B. Wehinger, M. Udina, V. Unikandanunni, T. Tadano, M. C. Hoffmann, A. V. Balatsky, and S. Bonetti, Terahertz electric-field-driven dynamical multiferroicity in SrTiO 3, Nature 628, 534 (2024)

    work page 2024

  40. [40]

    C. S. Davies, F. G. N. Fennema, A. Tsukamoto, I. Raz- dolski, A. V. Kimel, and A. Kirilyuk, Phononic switching of magnetization by the ultrafast Barnett effect, Nature 628, 540 (2024)

    work page 2024

  41. [41]

    M. F. Biggs, Sin-hang, Ho, A. Alejandro, M. Lutz, C. D. Moss, and J. A. Johnson, Ultrafast Faraday Rotation Probe of Chiral Phonon-Polaritons in LiNbO 3 (2025), 2507.22232 [cond-mat]

    work page arXiv 2025

  42. [42]

    Castro, C

    M. Castro, C. Saji, G. Saez, P. Vergara, S. Allende, and A. S. Nunez, Elementary theory of Magnetoferrons: Bringing magnons and ferrons together in multiferroic systems (2024), 2412.15796 [cond-mat]

    work page arXiv 2024

  43. [43]

    M. A. Castro, C. Saji, G. Saez, P. Vergara, S. Allende, and A. S. Nunez, Phenomenological theory of electro- magnons in multiferroic systems, Phys. Rev. B 111, 214401 (2025)

    work page 2025

  44. [44]

    Shitade, H

    A. Shitade, H. Watanabe, and Y. Yanase, Theory of or- bital magnetic quadrupole moment and magnetoelectric susceptibility, Phys. Rev. B 98, 020407 (2018)

    work page 2018

  45. [45]

    Bhowal and N

    S. Bhowal and N. A. Spaldin, Ferroically Ordered Mag- netic Octupoles in d-Wave Altermagnets, Phys. Rev. X 14, 011019 (2024)

    work page 2024

  46. [46]

    K. Hart, R. Sutcliffe, G. Refael, and A. Paramekanti, Phonon-Driven Multipolar Dynamics in a Spin-Orbit Coupled Mott Insulator, Phys. Rev. Lett. 134, 246701 (2025)

    work page 2025

  47. [47]

    Sutcliffe, K

    R. Sutcliffe, K. Hart, S. Chaudhary, and A. Paramekanti, Pseudo-chiral phonon splitting from octupolar magnetic order (2025), 2506.18978 [cond-mat]

    work page arXiv 2025

  48. [48]

    Urru, J.-R

    A. Urru, J.-R. Soh, N. Qureshi, A. Stunault, B. Roessli, H. M. Rønnow, and N. A. Spaldin, Neutron scattering 7 from local magnetoelectric multipoles: A combined the- oretical, computational, and experimental perspective, Phys. Rev. Res. 5, 033147 (2023)

    work page 2023

  49. [49]

    H. Ueda, A. Nag, C. P. Romao, M. Garc ´ ıa-Fern´ andez, K.- J. Zhou, and U. Staub, Chiral Phonons in Polar LiNbO 3 (2025), 2504.03330 [cond-mat]

    work page arXiv 2025

  50. [50]

    P. Tang, R. Iguchi, K.-i. Uchida, and G. E. W. Bauer, Thermoelectric Polarization Transport in Ferroelectric Ballistic Point Contacts, Phys. Rev. Lett. 128, 047601 (2022). Supplemental Material: Multiferrons: lattice excitations with finite polarization and magnetization Mike Pols, 1, ∗ Carl P. Romao, 2, † and Dominik M. Juraschek 1, ‡ 1Department of Appli...

    work page 2022

  51. [51]

    Computational details S3 Density functional theory (DFT) S3 Phonon calculations S3 Ferroelectric polarization S5

    work page

  52. [52]

    Phonon potential energy surfaces (PES) S7 Transforming degenerate PES to cylindrical coordinates S7 Phonon anharmonicity S9

    work page

  53. [53]

    Single-phonon per unit cell properties S12 Single-phonon per unit cell calculations S12 (Multi)ferron properties S12

    work page

  54. [54]

    Laser pulses S14 Functional form S14 Laser energy S14

    work page

  55. [55]

    Polarization and magnetization dynamics S16 High-frequency modes S16 Tunability of multiferrons S17

    work page

  56. [56]

    Multipole moments S18 Quadrupole moments S18 Octupole moments S20 References S27 S2

    work page

  57. [57]

    We use the standard PA W pseudopotentials with the valence electron configurati ons Li (2s 1), Nb (4p 65s14d4), and O (2s 22p4) [ 5]

    COMPUT A TIONAL DET AILS Density functional theory (DFT) We calculate the properties of LiNbO 3 using density functional theory (DFT) as im- plemented in the Vienna Ab-initio Simulation Package ( VASP) [ 1–4]. We use the standard PA W pseudopotentials with the valence electron configurati ons Li (2s 1), Nb (4p 65s14d4), and O (2s 22p4) [ 5]. Exchange-corre...

    work page

  58. [58]

    , 4) are the anharmonic coupling coefficients

    PHONON POTENTIAL ENERGY SURF ACES (PES) T ransforming degenerate PES to cylindrical coordinates We model the potential energy surface (PES) of two degenerat e E modes in LiNbO 3 as V (Qa, Q b) = ω 2 a 2 Q2 a + ω 2 b 2 Q2 b + a1Q3 a + a2Q2 aQb + a3QaQ2 b + a4Q3 b (S2) where Qa and Qb are the mode amplitudes of the degenerate modes, ω a = ω b ≡ ω is the pho...

    work page

  59. [59]

    As such, we use a value of ⟨Q2⟩ = ℏ ω to compute the polarization associated with a single phonon per unit cell with Eq

    SINGLE-PHONON PER UNIT CELL PROPERTIES Single-phonon per unit cell calculations Treating phonons as a quantum mechanical harmonic oscillat or, we can relate the mean- squared displacement of a mode ⟨Q2⟩ to the angular frequency ω of that mode using ladder operators, a and a†, as ⟨Q2⟩ = ⟨n|Q2|n⟩ = ℏ 2ω ⟨n| ( a + a†) 2 |n⟩ = ℏ 2ω ⟨n|a2 + ( a†) 2 + 2a†a + 1|...

    work page

  60. [60]

    54 THz mode, which has the largest in-plane polarization, it i s tilted by ϕ = 0. 18◦. S12 TABLE S4. Calculated phonon frequencies ν0, mode effective charges |Z|, and net single-phonon per unit cell polarization |P| of the A1 and E modes in LiNbO 3. Single-phonon per unit cell magnetization |M| of circularly polarized E modes in LiNbO 3. Symmetry Mode νDFT...

    work page

  61. [61]

    The electric field is typ ically represented as a vector with components that vary in time

    LASER PULSES F unctional form To model a laser pulse and its interaction with a material, we describe the time-dependent electric field generated by the pulse. The electric field is typ ically represented as a vector with components that vary in time. Here, we describe a laser pu lse propagating along the z-axis by E(t) = ˜E0 √ 2 exp      − (t − t0)2 (...

    work page

  62. [62]

    Ex (t) = Ey (t) = ˜E0√ 2 cos (ω 0t)

    For a linearly polarized laser pulse ( ϕ = 0), the x- and y-components of the electric S14 field are identical, i.e. Ex (t) = Ey (t) = ˜E0√ 2 cos (ω 0t). The resulting instantaneous intensity is |E(t)|2 = ˜E2 0 cos2(ω 0t). Averaging over a complete cycle of the laser pulse yields ⟨|E(t)|2⟩ = ˜E2 0 ⟨cos2(ω 0t)⟩ = ˜E2 0 2

    work page

  63. [63]

    The resulting instantaneous intensity is constant in time: |E (t) |2 = ˜E2 0 2 [ cos2 (ω 0) + sin2 (ω 0t) ] = ˜E2 0 2

    Circular excitations ( ϕ = ± π/ 2) have out-of-phase x- and y-components for the elec- tric field, for instance Ex (t) = ˜E0√ 2 cos (ω 0t) and Ey (t) = ˜E0√ 2 sin (ω 0t). The resulting instantaneous intensity is constant in time: |E (t) |2 = ˜E2 0 2 [ cos2 (ω 0) + sin2 (ω 0t) ] = ˜E2 0 2

    work page

  64. [64]

    Over a full cy- cle, it averages to ⟨|E(t)|2⟩ = ˜E2 0 2

    Elliptical excitations (0 < |ϕ| < π/ 2) have x- and y-components for the electric field of Ex (t) = ˜E0√ 2 cos (ω 0t) and Ey (t) = ˜E0√ 2 cos (ω 0t + ϕ), yielding an instanta- neous intensity that can be rewritten using trigonometric i dentities to: |E (t) |2 = ˜E2 0 2 [cos2 (ω 0t) + cos2 (ω 0t + ϕ)] = ˜E2 0 4 [2 + cos(2ω 0t) − cos(2ω 0t + 2ϕ)]. Over a ful...

    work page

  65. [65]

    32 THz, we also present the dynamics of modes 17 and 18 in Fig

    POLARIZA TION AND MAGNETIZA TION DYNAMICS High-frequency modes To demonstrate that the polarization and magnetization dyn amics shown in the main text are not unique to modes 4 and 5 at 4 . 32 THz, we also present the dynamics of modes 17 and 18 in Fig. S5. These modes have a frequency of 10 . 54 THz, as evidenced by the more rapid fluctuations in transien...

    work page

  66. [66]

    We used Eq

    MUL TIPOLE MOMENTS Quadrupole moments For completeness, we here show all quadrupole tensor compon ents for the radial polar- ization Pph and radial magnetization Mrad. We used Eq. (10) and Eq. (11) to determine the quadrupole contributions shown in Fig. S7 and Fig. S8. -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 P ij ( C Å / cm 2)(a) xx ×102 (b) xy ×102 (c) xz ×102 = ...

    work page

  67. [67]

    Kresse and J

    G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993)

    work page 1993

  68. [68]

    Kresse and J

    G. Kresse and J. Hafner, Phys. Rev. B 49, 14251 (1994)

    work page 1994

  69. [69]

    Kresse and J

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

    work page 1996

  70. [70]

    Kresse and J

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

    work page 1996

  71. [71]

    Kresse and D

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

    work page 1999

  72. [72]

    J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scus eria, L. A. Constantin, X. Zhou, and K. Burke, Phys. Rev. Lett. 100, 136406 (2008)

    work page 2008

  73. [73]

    H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976)

    work page 1976

  74. [74]

    S. C. Abrahams, H. J. Levinstein, and J. M. Reddy, J. Phys. Chem. Solids 27, 1019 (1966)

    work page 1966

  75. [75]

    A. Togo, J. Phys. Soc. Jpn. 92, 012001 (2023)

    work page 2023

  76. [76]

    A. Togo, L. Chaput, T. Tadano, and I. Tanaka, J. Phys.: Condens. Matter 35, 353001 (2023)

    work page 2023

  77. [77]

    Ridah, M

    A. Ridah, M. D. Fontana, and P. Bourson, Phys. Rev. B 56, 5967 (1997)

    work page 1997

  78. [78]

    M. R. Chowdhury, G. E. Peckham, and D. H. Saunderson, J. Phys. C: Solid State Phys. 11, 1671 (1978)

    work page 1978

  79. [79]

    Kokanyan, D

    N. Kokanyan, D. Chapron, and M. D. Fontana, Appl. Phys. A 117, 1147 (2014)

    work page 2014

  80. [80]

    R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 1651 (1993)

    work page 1993

Showing first 80 references.