pith. sign in

arxiv: 2512.14633 · v1 · pith:4EWE57KQnew · submitted 2025-12-16 · ❄️ cond-mat.mtrl-sci

On the origin of the unusual strain morphologies and polar Moir\'e patterns in twisted ferroelectrics

Sergey Prosandeev , Charles Paillard , Laurent Bellaiche This is my paper

Pith reviewed 2026-05-21 17:23 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords twisted ferroelectricsBaTiO3 bilayersshear strainMoiré patternselectric dipolesvortices and antivorticesacoustic forcesfirst-principles calculations
0
0 comments X

The pith

In twisted BaTiO3 bilayers, acoustic forces generate shear strain standing waves whose gradient couples to electric dipoles to form a Moiré pattern of vortices and antivortices.

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

This paper performs density functional theory calculations on twisted ferroelectric bilayers to trace the source of their complex shear strain shapes and polar Moiré structures. It shows that forces tied to acoustic motions create standing waves in the shear strain. These waves produce a strain gradient that interacts with the electric dipoles, resulting in interpenetrated arrays of vortices and antivortices. The analysis also indicates that forces linked to optical phonons play a lesser role in shaping these patterns.

Core claim

Density functional theory calculations, together with a decomposition of forces into acoustic and optical contributions, show that forces acting mostly on acoustic-related motions produce standing waves of the shear strain. Such waves naturally generate a self-organization of the shear strains and a peculiar gradient of these strains. A Moiré dipole pattern consisting of interpenetrated arrays of vortices and antivortices made of the electric dipoles then mostly arises due to the coupling of this gradient of the shear strain with the electric dipoles, while forces acting on optical-phonon motions play a smaller role.

What carries the argument

Decomposition of forces into acoustic-related and optical-phonon-related contributions that isolates the dominant mechanism producing standing shear-strain waves and the resulting dipole pattern.

If this is right

  • Acoustic forces are the primary driver of the self-organized shear strain morphologies in twisted ferroelectric bilayers.
  • The gradient of shear strain couples directly with electric dipoles to create the observed Moiré pattern of vortices and antivortices.
  • Optical-phonon forces contribute to the polar vortices and antivortices but only to a limited degree.
  • Comparable strain and dipole patterns are expected in other twisted ferroelectric materials when similar acoustic mechanisms operate.

Where Pith is reading between the lines

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

  • Controlling the twist angle may offer a way to engineer specific dipole vortex patterns for potential use in ferroelectric devices.
  • External electric fields could be tested to see whether they modulate the acoustic waves and thereby alter the Moiré dipole arrangement.
  • The strain-gradient mechanism may connect to topological features seen in other two-dimensional layered systems.

Load-bearing premise

Splitting the forces into those linked to acoustic motions versus those linked to optical motions correctly identifies what mainly causes the shear strain waves and the dipole patterns.

What would settle it

A simulation that applies the same twist but uses a different force decomposition or boundary conditions and fails to produce the observed standing shear-strain waves or the Moiré dipole pattern would challenge the proposed origin.

Figures

Figures reproduced from arXiv: 2512.14633 by Charles Paillard, Laurent Bellaiche, Sergey Prosandeev.

Figure 1
Figure 1. Figure 1: FIG. 1. (Color online) Forces centered on the Ti ions in the bottom atomic layer of the BTO [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (Color online) Map in a (x,y) plane of the shear strain morphology from Eq. (7) (panel a), [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
read the original abstract

Density functional theory calculations are conducted to understand and reveal the origin of the complex shear strain morphology and of the polar Moir\'e topological pattern recently observed in twisted BaTiO$_3$ bilayers. Our first-principles calculations, along with an original analysis of them allowing the decomposition of forces into the acoustic and optical contributions, point out to the occurrence of forces mostly acting on the {\it acoustic-related} motions to produce the standing waves of the shear strain. Such acoustic waves naturally generate a striking self-organization of the shear strains, and hence create a peculiar gradient of these shear strains. A Moir\'e dipole pattern, consisting of the interpenetrated arrays of vortices and antivortices made of the electric dipoles, then mostly arises due to the coupling of this gradient of the shear strain with the electric dipoles. Furthermore, other forces, namely acting on the motions associated with the {\it optical phonons}, could also play a role in the formation of these polar vortices and antivortices, but at a smaller extent.

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 uses density functional theory to examine the origin of unusual shear-strain morphologies and polar Moiré patterns observed in twisted BaTiO3 bilayers. It introduces an original force decomposition separating acoustic-related and optical-phonon-related contributions, concluding that forces acting primarily on acoustic motions generate standing shear-strain waves; the resulting strain gradient then couples to electric dipoles to produce the interpenetrated vortex-antivortex Moiré dipole pattern, while optical contributions play a secondary role.

Significance. If the force decomposition can be shown to cleanly isolate the dominant mechanism, the work would supply a first-principles mechanistic account of how strain gradients drive polar topology in twisted ferroelectrics. The combination of DFT with a mode-based force analysis is a positive feature that could inform design of moiré ferroelectric heterostructures.

major comments (2)
  1. [force decomposition analysis] The central attribution of the shear-strain waves and subsequent dipole pattern to acoustic-related forces rests on an unspecified decomposition procedure. The manuscript does not state whether the projection uses eigenvectors of the untwisted bilayer dynamical matrix, a supercell phonon-mode decomposition, or a real-space filter, nor does it report tests of stability under small changes in reference structure or cutoff. Without this information the claim that acoustic forces 'mostly' dominate cannot be quantitatively assessed.
  2. [results on dipole pattern formation] The manuscript states that the gradient of the shear strain couples to the electric dipoles to form the Moiré vortex-antivortex pattern, yet no explicit calculation or correlation function is shown that quantifies the relative strength of this coupling versus other possible contributions (e.g., direct piezoelectric or flexoelectric terms). A concrete decomposition of the dipole field into strain-gradient-driven and other components is needed to support the 'mostly arises due to' statement.
minor comments (2)
  1. [methods] The abstract and main text should explicitly list the exchange-correlation functional, plane-wave cutoff, k-point sampling, and convergence criteria employed in the DFT calculations.
  2. [force decomposition analysis] Notation for the acoustic and optical force components should be defined once and used consistently; currently the distinction is introduced only descriptively.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable feedback on our manuscript. We address each major comment below and have revised the manuscript to provide the requested details and supporting analyses.

read point-by-point responses

  1. Referee: [force decomposition analysis] The central attribution of the shear-strain waves and subsequent dipole pattern to acoustic-related forces rests on an unspecified decomposition procedure. The manuscript does not state whether the projection uses eigenvectors of the untwisted bilayer dynamical matrix, a supercell phonon-mode decomposition, or a real-space filter, nor does it report tests of stability under small changes in reference structure or cutoff. Without this information the claim that acoustic forces 'mostly' dominate cannot be quantitatively assessed.

    Authors: We agree that the original manuscript lacked sufficient detail on the force decomposition. The procedure projects the DFT-computed forces onto the acoustic and optical eigenvectors of the dynamical matrix calculated for the untwisted BaTiO3 bilayer reference structure, performed within the supercell. We have added a dedicated subsection in the Methods section that specifies the projection formula, the mode classification criteria, and explicit tests of robustness under small shifts in the reference atomic positions and force cutoff. These additions confirm the dominance of the acoustic contribution and are supported by new supplementary figures. revision: yes

  2. Referee: [results on dipole pattern formation] The manuscript states that the gradient of the shear strain couples to the electric dipoles to form the Moiré vortex-antivortex pattern, yet no explicit calculation or correlation function is shown that quantifies the relative strength of this coupling versus other possible contributions (e.g., direct piezoelectric or flexoelectric terms). A concrete decomposition of the dipole field into strain-gradient-driven and other components is needed to support the 'mostly arises due to' statement.

    Authors: We acknowledge that the original text did not include a quantitative decomposition of the dipole contributions. We have now performed an explicit decomposition by isolating the dipole field generated from the computed shear-strain gradient via the flexoelectric response and comparing it against the total DFT dipole field as well as the direct piezoelectric term. Correlation analysis shows that the strain-gradient coupling is the leading mechanism. This decomposition, together with supporting plots, has been added to the Results section and a new figure in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central claim follows from DFT analysis without reduction to inputs by construction.

full rationale

The derivation relies on density functional theory calculations of twisted BaTiO3 bilayers, followed by decomposition of forces into acoustic and optical contributions to identify shear-strain waves and their coupling to dipoles. No step equates a reported prediction or pattern to a parameter fitted from the same data, nor does any load-bearing premise reduce to a self-citation chain or ansatz smuggled via prior work by the same authors. The force decomposition is presented as an original analysis of the computed results rather than a definitional tautology, and the Moiré dipole pattern is attributed directly to the gradient of the computed shear strains. This is a standard first-principles workflow with independent content.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on standard DFT assumptions such as the Born-Oppenheimer approximation and chosen exchange-correlation functional; no ad-hoc free parameters or new entities are introduced in the abstract.

axioms (1)
  • standard math Born-Oppenheimer approximation separating electronic and nuclear motion
    Invoked implicitly by all DFT calculations of atomic forces and phonons.

pith-pipeline@v0.9.0 · 5727 in / 1275 out tokens · 56957 ms · 2026-05-21T17:23:02.836634+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

28 extracted references · 28 canonical work pages

  1. [1]

    S. Carr, D. Massatt, S. Fang, P. Cazeaux, M. Luskin, E. Kaxiras, Twistronics: Manipulating the electronic properties of two-dimensional layered structures through their twist angle. Phys. Rev. B.95, 075420 (2017)

    work page 2017

  2. [2]

    Oudich, X

    M. Oudich, X. Kong, T. Zhang, C. Qiu, and Y. Jing, Engineered moir´ e photonic and phononic superlattices. Nature Materials23, 1169 (2024)

    work page 2024

  3. [3]

    F. He, Y. Zhou, Z. Ye, S.-H. Cho, J. Jeong, X. Meng, and Y. Wang, Moir´ e Patterns in 2D Materials: A Review, ACS Nano15, 5944 (2021)

    work page 2021

  4. [4]

    Y. Cao, D. Rodan-Legrain, O. Rubies-Bigorda, J. Min Park, K. Watanabe, T. Taniguchi, P. Jarillo-Herrero, Tunable correlated states and spin-polarized phases in twisted bilayer-bilayer graphene, Nature583, 215 (2020)

    work page 2020

  5. [5]

    Y. Jia, J. Yu, J. Liu, J. Herzog-Arbeitman, Z. Qi, H. Pi, N. Regnault, H. Weng, B. A. Bernevig, and Q. Wu, Moir´ e fractional Chern insulators. I. First-principles calculations and continuum models of twisted bilayer MoTe 2, Phys. Rev. B109, 205121 (2024)

    work page 2024

  6. [6]

    Herzog-Arbeitman, Y

    J. Herzog-Arbeitman, Y. Wang, J. Liu, P. M. Tam, Z. Qi, Y. Jia, D. K. Efetov, O. Vafek, N. Regnault, H. Weng, Q. Wu, B. A. Bernevig, and J. Yu, Moir´ e fractional Chern insulators. II. First-principles calculations and continuum models of rhombohedral graphene superlattices, Phys. Rev. B109, 205122 (2024). 12

    work page 2024

  7. [7]

    E. Y. Andrei , D. K. Efetov , P. Jarillo- Herrero, A. H. MacDonald, K. F. Mak , T. Senthil, E. Tutuc, A. Yazdani, and A. F. Young, The marvels of moir´ e materials, Nat. Rev. Mat.6, 201 (2021)

    work page 2021

  8. [8]

    X. Yan, Z. Zheng, V. K. Sangwan, J. H. Qian, X. Wang, S. E. Liu, K. Watanabe, T. Taniguchi, S.-Y. Xu, P. Jarillo-Herrero, Q. Ma, and M. C. Hersam, Moir´ e synaptic transistor with room- temperature neuromorphic functionality. Nature624, 551 (2023)

    work page 2023

  9. [9]

    Fan, Moir´ e turns neuromorphic

    W. Fan, Moir´ e turns neuromorphic. Nature Materials23, 1164 (2024)

    work page 2024

  10. [10]

    S´ anchez-Santolino, V

    G. S´ anchez-Santolino, V. Rouco, S. Puebla, H. Aramberri, V. Zamora, M. Cabero, F. A. Cuellar, C. Munuera, F. Mompean, M. Garcia-Hernandez, A. Castellanos-Gomez, J. ´I˜ niguez, C. Leon, and J. Santamaria, A 2D ferroelectric vortex pattern in twisted BaTiO3 freestanding layers, Nature (London)626, 529 (2024)

    work page 2024

  11. [11]

    S. M. Kogan, Flexoelectric polarization of barium strontium titanate in the paraelectric state, Sov. Phys. Solid State5, 2069 (1964)

    work page 2069

  12. [12]

    Zubko, G

    P. Zubko, G. Catalan, A. K. Tagantsev, Flexoelectric Effect in Solids, Ann. Rev. Mater. Res. 43, 387 (2013)

    work page 2013

  13. [13]

    Prosandeev, C

    S. Prosandeev, C. Paillard, and L. Bellaiche, Understanding and controlling dipolar Moir´ e pattern in ferroelectric perovskite oxide nanolayers. Phys. Rev. B111, L180103 (2025)

    work page 2025

  14. [14]

    S. Lee, D. J. P. de Sousa, B. Jalan, T. Low, Moir´ e polar vortex, flat bands, and Lieb lattice in twisted bilayer BaTiO 3, Sci. Adv.10, 293 (2024)

    work page 2024

  15. [15]

    J. M. Kosterlitz and D. J. Thouless, Ordering, metastability and phase transitions in two- dimensional systems, J. Phys. C6, 1181 (1973)

    work page 1973

  16. [16]

    V. L. Berezinskii, Destruction of Long-range Order in One-dimensional and Two-dimensional Systems Possessing a Continuous Symmetry Group. II. Quantum Systems, ZhETF61, 1144 (1972) [Sov. Phys. JETP34, 610 (1972)]

    work page 1972

  17. [17]

    W. Kohn, L. J. Sham, Self-consistent equations including exchange and correlation effects. Phys. Rev.140, 1133 (1965)

    work page 1965

  18. [18]

    Kresse and J

    G. Kresse and J. Hafner. Ab initio molecular-dynamics simulation of the liquid-metal- amorphous-semiconductor transition in germanium. Phys. Rev. B49, 14251 (1994)

    work page 1994

  19. [19]

    Furthmuller

    G Kresse and J. Furthmuller. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science6, 15 (1996)

    work page 1996

  20. [20]

    P. E. Blochl. Projector augmented-wave method. Phys. Rev. B50, 17953 (1994). 13

    work page 1994

  21. [21]

    From ultrasoft pseudopotentials to the projector augmented-wave method

    G Kresse and D Joubert. From ultrasoft pseudopotentials to the projector augmented-wave method. Physical Review B59, 1758 (1999)

    work page 1999

  22. [22]

    J. P. Perdew, K. Burke and M. Ernzerhof. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett.77, 3865 (1996)

    work page 1996

  23. [23]

    R. P. Feynman, Forces in Molecules. Phys. Rev.56, 340 (1939)

    work page 1939

  24. [24]

    Zhong, D

    W. Zhong, D. Vanderbilt, and K. M. Rabe, Phase Transitions in BaTiO3 from First Principles. Phys. Rev. Lett.73, 1861 (1994)

    work page 1994

  25. [25]

    Zhong, D

    W. Zhong, D. Vanderbilt, and K. M. Rabe, First-principles theory of ferroelectric phase tran- sitions for perovskites: The case of BaTiO 3. Phys. Rev. B52, 6301 (1995)

    work page 1995

  26. [26]

    Bellaiche, A

    L. Bellaiche, A. Garc ˜Aa, and D. Vanderbilt, Finite-Temperature Properties of Pb(Zr1−xTix)O3 Alloys from First Principles. Phys. Rev. Lett.84, 5427 (2000)

    work page 2000

  27. [27]

    R. D. King-Smith and D. Vanderbilt, First-principles investigation of ferroelectricity in per- ovskite compounds. Phys. Rev. B49, 5828 (1994)

    work page 1994

  28. [28]

    Bastogne, F

    L. Bastogne, F. G´ omez-Ortiz, S. Anand, and P. Ghosez, Dynamical manipulation of polar topologies from acoustic phonon excitations, Nano Lett.24, 13783 (2024). 14

    work page 2024