Morphology of Polarization States in Strained Ferroelectric Films

Ana\"is Sen\'e; Anna Razumnaya; Igor Lukyanchuk; L\'eo Boron; Svitlana Kondovych; Yuri Tikhonov

REVIEW 3 major objections 5 minor 47 references

Morphology of Polarization States in Strained Ferroelectric Films

T0 review · 3 major / 5 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper claims that in strained PbTiO3 films the equilibrium polarization texture is captured by a lowest-harmonic ansatz, and that compressive strain stabilizes a rotated vortex phase while small tensile strain stabilizes a rotated helix

desk verdict A useful analytical shortcut for strained-PbTiO3 phase diagrams, but the predicted h' helix region is contradicted by the paper's own phase-field simulation, so treat the quantitative boundaries cautiously. read the letter →

arxiv 2509.06508 v1 pith:Y45BBUIJ submitted 2025-09-08 cond-mat.mtrl-sci cond-mat.mes-hall

Léo Boron , Anaïs Sené , Yuri Tikhonov , Anna Razumnaya , Igor Lukyanchuk , Svitlana Kondovych This is my paper · ORCID

classification cond-mat.mtrl-scicond-mat.mes-hall PACS 77.80.-e

keywords ferroelectricthinfilmstopologicalpolarizationtexturesvortexphasesoft-domainapproachGinzburg-Landau-DevonshireepitaxialstrainPbTiO3diagrams

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper claims that the vortex-like polarization patterns seen in strained ferroelectric films can be captured by a single three-amplitude Fourier ansatz. Working with the Ginzburg-Landau-Devonshire free energy and strain-renormalized coefficients for PbTiO3, the authors impose a divergence-free condition, volume-average, and minimize. The resulting phase diagrams show that under compressive misfit strain a diagonally oriented vortex phase replaces the uniform out-of-plane c-phase, while small tensile strain stabilizes a rotated helix; ultrathin films or larger tensile strain return uniform in-plane phases. The value is practical: a compact analytic model gives temperature-strain-thickness phase boundaries that full simulations and experiments can test without heavy computation.

What carries the argument

The soft-domain ansatz (Eq. 14) is the central object: a lowest-harmonic sinusoidal polarization field P = (Pa sin(pi z/2a_f) sin(pi x/d) + P1, P2, P3 cos(pi z/2a_f) cos(pi x/d)). The electrostatic constraint div P = 0 collapses into the relation Pa = gamma P3 with gamma = d/2a_f, which eliminates surface and bulk bound charges. Volume averaging with standard trigonometric mean values turns the full GLD functional into a polynomial in three amplitudes whose coefficients are the strain-renormalized effective-functional values, with a 45-degree rotation in-plane to generate diagonal vortex states. This reduction is what makes analytic phase diagrams possible.

What would settle it

Take a PbTiO3 film at thickness 12 nm, T = 300 K, and compressive misfit strain u_m = -0.01; compute the equilibrium with phase-field or atomistic simulation and expand the out-of-plane polarization profile across the thickness in Fourier modes. If the fundamental cosine/sine mode is not dominant, or if the vortex array does not align along an in-plane diagonal, the ansatz's phase diagram is falsified.

Watch

Extended reading notes

Core claim

Within the soft-domain framework the polarization texture is written as the lowest Fourier harmonics: an in-plane vortex component proportional to sin(pi z/2a_f) sin(pi x/d), a uniform in-plane part, and an out-of-plane component proportional to cos(pi z/2a_f) cos(pi x/d). Enforcing vanishing divergence fixes the vortex amplitude Pa = gamma P3, with gamma = d/2a_f, so the whole state is controlled by three amplitudes plus one geometric ratio. After volume-averaging the strain-renormalized free energy over the domain cell, minimization yields phase diagrams in which the rotated variants v' and h' appear where the uniform c- and r-phases would have been electrostatically forbidden. The same an

Load-bearing premise

The ansatz forces the equilibrium texture to be a single sinusoidal harmonic in x and z and to be exactly divergence-free; if higher harmonics, multiple wave vectors, or temperature-dependent profile changes alter the energy ordering, the predicted phase boundaries could be wrong.

Editorial extensions

If this is right

Compressive misfit strain stabilizes a diagonally rotated vortex lattice v' rather than the uniform c-phase, with a transition temperature depressed relative to the monodomain value.
Small tensile strain opens a rotated helix h' window; because neighboring states are nearly degenerate, this region should appear as labyrinthine or mixed textures.
Ultrathin films (below about 5 nm) favor uniform in-plane aa polarization under the chosen conditions, while nonuniform vortex states dominate thicker films.
The same three-amplitude machinery gives a morphological classification: v, h, w, and wh states correspond to different combinations of uniform in-plane polarization with vortex arrays.
Phase-field simulations with sinusoidal profiles and near-zero divergence support the soft-domain predictions, including a1/a2 stripe relaxation at large tensile strain.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the single-harmonic ansatz is correct, higher Fourier modes should be small in real films; a direct Fourier analysis of atomistic polarization maps would make this testable.
Extending the plane-wave superposition to a threefold star of wave vectors is the paper's own sketch for bubbles; one could push the same averaged-energy machinery to predict bubble-to-vortex transitions quantitatively.
The near-degeneracy of h, w, hw, and w' states at small strain suggests external fields or slight strain anisotropy can be used to switch among chiral morphologies, an implication the authors leave implicit.
The method's reliance on the monodomain strain renormalization means it should be revisited when full elastic inhomogeneity or flexoelectric coupling is relevant; the paper names these as future extensions.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

This paper delivers something genuinely useful: a compact variational route to phase diagrams for strained PbTiO3 films, built from the authors' earlier soft-domain ansatz and the standard Pertsev-Tagantsev coefficients. No free parameters are fitted to the target phases, the algebra is laid out in enough detail to follow, and the classification of vortex, helix, wave, and combined states is clean. The v'-phase replacing the uniform c-phase under compressive strain is plausible and consistent with the phase-field test they show. That alone is worth having as a fast estimator for experimental design.

The weak spot is the h' region near small tensile strain. Their own phase-field simulation at u_m ≈ 0 does not relax to the predicted regular helix; it gives a labyrinthine wh-like mixture. The paper's response—these states are nearly degenerate and the exact pattern is sensitive to history—is reasonable, but it does not rescue the specific phase boundary. It means the h' stability region is an artifact of the single-mode restriction, not a demonstrated equilibrium of the full GLD model. The same truncation underlies the whole diagram, so while the v' and aa regions are probably robust, the h' region and any boundaries involving it are not quantitatively trustworthy. The sparse phase-field validation (three points, parameters not fully disclosed) reinforces that concern. To their credit, the authors openly flag the discrepancy and discuss the need for multi-k superpositions for bubble states, so the paper is honest about the limitation.

Who should read it: anyone working on ferroelectric topological textures who wants a quick analytic complement to phase-field or atomistic simulations. The paper is clearly written and the method is transparent. For peer review, it deserves a serious referee, but the referee should push on the h' discrepancy and on whether the one-mode ansatz can be relaxed or at least benchmarked more systematically. I would not desk-reject it, but I would expect major revision before the phase diagrams are presented as predictive.

Referee Report

3 major / 5 minor

Summary. The authors extend the soft-domain variational approach to construct phase diagrams of polarization textures in epitaxial PbTiO3 films under misfit strain. Using a single-mode divergence-free ansatz (Eq. 14), they derive an effective GLD free energy with strain-renormalized coefficients, minimize it, and find that a rotated vortex phase v' replaces the uniform c-phase under compression, a rotated helix h' appears near small tensile strain, and the uniform aa-phase is stabilized at larger tensile strain and in ultrathin films. They support the picture with three phase-field simulations and compare the sinusoidal profile with experimental data. The paper is clearly written and provides explicit coefficient tables.

Significance. If the phase diagrams were reliable, this would be a valuable compact predictive tool for strained ferroelectric films, complementing heavy numerical simulations. Strengths: the calculations use literature GLD coefficients with no fitted parameters; the ansatz enforces div P = 0 exactly; explicit coefficient tables make the results reproducible; the phase-field simulations provide a benchmark; and the experimental profile comparison in Fig. 5b is a strong positive point. The main caveat is that the phase diagrams are variational results within a restricted single-mode class, and the paper's own phase-field test at u_m ≈ 0 does not reproduce the predicted h' phase, so the central phase-stability claims need further support or qualification.

major comments (3)

[Sec. IV.B (Fig. 4b) vs Sec. IV.A (Fig. 3a,c)] The phase-field simulation at u_m ≈ 0, which lies inside the predicted h' region of Fig. 3, relaxes to a labyrinthine wh-like state rather than the regular h' helix. The paper's attribution to a 'transient regime with almost degenerate energies' (Sec. V.B) does not clarify whether h' is an equilibrium phase of the full GLD model. This is a direct test of the variational ansatz, and it fails in the very region where one of the two central claims — the emergence of h' — is made. Please either (i) demonstrate that seeding h' in the phase-field model relaxes to a stable h' state with the same parameters, or (ii) revise the phase diagrams and central claims to describe h' as metastable or as an upper-bound variational result, adding a corresponding caveat in the abstract and conclusions.
[Sec. III.A–C, Eq. (14)] The entire phase diagram is obtained by minimizing within the one-mode ansatz (14), which fixes a sinusoidal z-profile and a single wave vector. Higher harmonics and superpositions of wave vectors are excluded, although the paper itself notes in Sec. V.A that bubble states require a star of wave vectors. Since phase boundaries are determined by comparing energies of different ansatz states, omitted modes can change not only the boundaries but also the ordering of phases. The profile-shape assumption is stated in Sec. III.A ("in the spirit of Landau theory") but is not tested away from T_c; the phase diagrams extend to room temperature and below, where nonlinearities are strong. Please add a convergence check (e.g., include a third harmonic in the ansatz) or perform phase-field relaxations at several representative points in each phase region of Fig. 3 to validate the phase boundaries.
[Sec. IV.B (Fig. 4c) and Sec. III.B/IV.A (aa phase)] At tensile misfit strain u_m = +0.01, the phase-field result is a parquet-like a1/a2 domain arrangement, not the uniform aa-phase shown in the phase diagrams of Fig. 3. The paper acknowledges an effective equivalence on long length scales, but if the phase diagram labels an actual texture, this is another discrepancy; if the aa-phase is a coarse-grained description, this should be stated explicitly wherever the phase diagram is presented (Fig. 3, Sec. IV.A). Without this clarification, readers may overinterpret the tensile side of the phase diagram.

minor comments (5)

[Sec. II.C] The gradient coefficients appear rounded inconsistently: with G1111 = 2.77, G1122 = 0, G1212 = 1.38, one obtains K1 = K2 = 2.76 and K_a = 0.01, whereas the text gives K1 = 2.77 and states K_a = 0. Please reconcile the numerical values or qualify the statement that the anisotropy term 'vanishes' for PbTiO3.
[Fig. 4b caption] The caption labels the texture 'Wave-helix wh-phase', but the main text states that the simulation produces a labyrinthine mixture closer to a distorted wh-phase. Since this figure is the key comparison, the caption should reflect the actual result (e.g., 'labyrinthine wh-like pattern') to avoid confusing readers.
[Sec. IV.A] The minimization with respect to the structural parameter γ is not explicitly described. Please state whether all phase diagrams include minimization over γ and over the amplitudes, and whether multiple local minima were checked.
[Sec. V.A] The mapping of experimental observations to the unprimed v, h, w, and wh phases should be qualified: those experiments are performed in superlattices with different boundary conditions, and the computed phases in Fig. 3 are primarily the rotated v' and h' phases. The match may therefore be more qualitative than the current wording suggests.
[Throughout] Several minor typographical issues: 'F erroelectric' in the Fig. 1 caption, inconsistent use of the film-thickness notation '2a_f' in figure captions, and a few accent/rendering issues in the author list (e.g., 'Léo Boron', 'Anaïs Sené').

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: phase diagrams are variational energy minimizations within an explicitly stated ansatz, benchmarked by phase-field simulation, not fitted to the predicted phases.

full rationale

The paper's derivation chain is: GLD functional with literature coefficients → soft-domain ansatz (14) → volume-averaged effective free energy (Tables II and III) → minimization → phase diagrams (Fig. 3). The ansatz is an explicit variational restriction, not a hidden result, and the phase-diagram predictions follow from minimizing the energy over that trial family using external literature coefficients (Pertsev–Tagantsev, Wang et al.). No free parameter is tuned to reproduce the v′ or h′ phases, and the paper does not claim the ansatz is exact; it acknowledges in Sec. V.A that bubble/skyrmion states require a star of wave vectors and in Sec. V.B that the h′ region is actually realized as labyrinthine mixed states in phase-field simulation (Fig. 4b). The self-citations to Refs. [23,25,27] and [33] are historical lineage for the soft-domain method rather than load-bearing evidence: the present paper re-derives the effective coefficients and compares with an independent numerical implementation. The single-harmonic approximation is a validity limitation, not a circular step, because the conclusion does not reduce to an input by construction. Thus no circularity is present, and the appropriate score is 0.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The model contains no fitted free parameters: the ansatz amplitudes and the geometric parameter gamma are variational and determined by energy minimization, and all material coefficients are taken from prior literature. The main postulates are the single-mode ansatz shape and the use of strain-renormalized uniform coefficients for nonuniform textures. No new physical entities are introduced.

assumptions (6)

domain assumption GLD free energy with Pertsev-Tagantsev strain-renormalized coefficients correctly describes uniform and nonuniform polarization in strained PbTiO3 films.
Sec. II D and III C: the P-T coefficients from Ref [20] are used for all textures, including nonuniform ones, even though the P-T approach was derived for uniform states.
ad hoc to paper The polarization texture is exactly the single-mode ansatz Eq. (14) with divergence-free constraint.
Sec. III B: this reduces the functional to three amplitudes; higher harmonics and other wave vector combinations are excluded.
domain assumption The soft-mode spatial profile remains unchanged as temperature is lowered below Tc.
Sec. III A: 'we assume that close to the transition the spatial profile of the inhomogeneous structure remains essentially unchanged,' citing Ref [27].
standard math Volume averaging of powers of sine and cosine (Eq. 16) over the domain cell is a valid reduction of the thermodynamic energy.
Sec. III C: averages used to obtain effective energy.
standard math The 45-degree rotation of coordinates and coefficients (Table III) preserves the relevant physics.
Sec. III C: transformation (17) and subsequent energy minimization.
domain assumption Enforcing div P=0 fully removes the depolarizing energy, so no explicit electrostatic energy term is needed.
Sec. II E and III B: the divergence-free condition is used to eliminate bound charges; any residual depolarizing field from higher-order effects is neglected.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Morphology of Polarization States in Strained Ferroelectric Films." pith.science (2026). https://pith.science/paper/Y45BBUIJ

@misc{pith2026250906508,
  author       = {Pith},
  title        = {Pith review of: Morphology of Polarization States in Strained Ferroelectric Films},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Y45BBUIJ}},
  note         = {Machine review of arXiv:2509.06508}
}

abstract

Ferroelectric thin films under epitaxial strain exhibit a variety of vortex-like topological polarization textures. To analyze them, we build on the Ginzburg-Landau-Devonshire framework and extend the previously introduced soft-domain approach. This formulation provides a compact variational theoretical description of polarization morphologies in strained PbTiO$_3$ films. It yields phase diagrams as a function of temperature, strain, and thickness, and clarifies the morphological structure of emergent topological states. The method is computationally efficient and offers practical guidance for experimental studies of ferroelectric nanostructures.

Figures

Figures reproduced from arXiv: 2509.06508 by the authors.

**Figure 2.** FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: a shows the vortex phase stabilized under compressive misfit strain um = −0.01. The polarization streamlines form a vortex array structure aligned approximately along the in-plane diagonal, corresponding to the v ′ -phase in the phase diagrams of [PITH_FULL_IMAGE:fi…

**Figure 5.** Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

47 extracted references · 43 canonical work pages

[1]

The realized soft mode corresponds to the maximal critical temperature of the solution of (13), and the as- sociated wave vectork ∗ sets the modulation periodicity. In the spirit of Landau theory, we assume that close to the transition the spatial profile of the inhomogeneous structure,P(r), remains essentially unchanged, while its amplitude grows upon co...

work page
[2]

Junquera, Y

J. Junquera, Y. Nahas, S. Prokhorenko, L. Bellaiche, J. ´I˜ niguez, D. G. Schlom, L.-Q. Chen, S. Salahuddin, D. A. Muller, L. W. Martin,et al., Topological phases in polar oxide nanostructures, Rev. Mod. Phys.95, 025001 (2023)

work page 2023
[3]

S. Das, Z. Hong, M. McCarter, P. Shafer, Y.-T. Shao, D. Muller, L. Martin, and R. Ramesh, A new era in fer- roelectrics, APL Mater.8, 120902 (2020)

work page 2020
[4]

I. A. Lukyanchuk, A. G. Razumnaya, S. Kondovych, Y. A. Tikhonov, B. Khesin, and V. M. Vinokur, Topo- logical foundations of ferroelectricity, Phys. Rep.1110, 1 (2025)

work page 2025
[5]

Luk’yanchuk, A

I. Luk’yanchuk, A. Razumnaya, S. Kondovych, Y. Tikhonov, and V. M. Vinokur, Topological fer- roelectric chirality, Preprint arXiv:2406.19728 (2024)

arXiv 2024
[6]

A. K. Yadav, K. X. Nguyen, Z. Hong, P. Garc ´ ıa- Fern´ andez, P. Aguado-Puente, C. T. Nelson, S. Das, B. Prasad, D. Kwon, S. Cheema,et al., Spatially re- solved steady-state negative capacitance, Nature565, 468 (2019)

work page 2019
[7]

Yadav, C

A. Yadav, C. Nelson, S. Hsu, Z. Hong, J. Clarkson, C. Schlep¨ utz, A. Damodaran, P. Shafer, E. Arenholz, L. Dedon,et al., Observation of polar vortices in oxide superlattices, Nature530, 198 (2016)

work page 2016
[8]

S. Das, Y. Tang, Z. Hong, M. Gon¸ calves, M. Mc- Carter, C. Klewe, K. Nguyen, F. G´ omez-Ortiz, P. Shafer, E. Arenholz,et al., Observation of room-temperature po- lar skyrmions, Nature568, 368 (2019)

work page 2019

Show all 47 references

[9]

Shafer, P

P. Shafer, P. Garc ´ ıa-Fern´ andez, P. Aguado-Puente, A. R. Damodaran, A. K. Yadav, C. T. Nelson, S.-L. Hsu, J. C. Wojde l, J. ´I˜ niguez, L. W. Martin,et al., Emergent chi- rality in the electric polarization texture of titanate su- perlattices, Proc. Natl. Acad. Sci.115, 915 (2018)

2018
[10]

murmurations

directions. This so-calleda 1/a2 domain structure is well-established in strained ferroelectric films [30–32] and represents the optimal configuration to minimize the in- homogeneous part of the elastic energy, which, as already noted, is not fully captured by the soft-domain ...
[11]

Y. Wang, Y. Feng, Y. Zhu, Y. Tang, L. Yang, M. Zou, W. Geng, M. Han, X. Guo, B. Wu,et al., Polar meron lattice in strained oxide ferroelectrics, Nat. Mater.19, 881 (2020)

2020
[12]

Gong, Y.-L

F.-H. Gong, Y.-L. Tang, Y.-L. Zhu, H. Zhang, Y.-J. Wang, Y.-T. Chen, Y.-P. Feng, M.-J. Zou, B. Wu, W.-R. Geng,et al., Atomic mapping of periodic dipole waves in ferroelectric oxide, Sci. Adv.7, eabg5503 (2021)

2021
[13]

Behera, M

P. Behera, M. A. May, F. G´ omez-Ortiz, S. Susarla, S. Das, C. T. Nelson, L. Caretta, S.-L. Hsu, M. R. Mc- Carter, B. H. Savitzky,et al., Electric field control of chirality, Sci. Adv.8, eabj8030 (2022)

2022
[14]

Lichtensteiger, S

C. Lichtensteiger, S. Fernandez-Pena, C. Weymann, P. Zubko, and J.-M. Triscone, Tuning of the depolariza- tion field and nanodomain structure in ferroelectric thin films, Nano Lett.14, 4205 (2014)

2014
[15]

Zhang, L

Q. Zhang, L. Xie, G. Liu, S. Prokhorenko, Y. Nahas, X. Pan, L. Bellaiche, A. Gruverman, and N. Valanoor, Nanoscale bubble domains and topological transitions in ultrathin ferroelectric films, Adv. Mater.29, 1702375 (2017)

2017
[16]

Nahas, S

Y. Nahas, S. Prokhorenko, Q. Zhang, V. Govinden, N. Valanoor, and L. Bellaiche, Topology and control of self-assembled domain patterns in low-dimensional ferro- electrics, Nat. Commun.11, 5779 (2020)

2020
[17]

Nahas, S

Y. Nahas, S. Prokhorenko, J. Fischer, B. Xu, C. Carr´ et´ ero, S. Prosandeev, M. Bibes, S. Fusil, B. Dkhil, V. Garcia,et al., Inverse transition of labyrinthine do- main patterns in ferroelectric thin films, Nature577, 47 (2020)

2020
[18]

I. I. Naumov, L. Bellaiche, and H. Fu, Unusual phase transitions in ferroelectric nanodisks and nanorods, Na- ture432, 737 (2004)

2004
[19]

D. G. Schlom, L.-Q. Chen, C.-B. Eom, K. M. Rabe, S. K. Streiffer, and J.-M. Triscone, Strain tuning of ferroelec- tric thin films, Annu. Rev. Mater. Res.37, 589 (2007)

2007
[20]

D. G. Schlom, L.-Q. Chen, X. Pan, A. Schmehl, and M. A. Zurbuchen, A thin film approach to engineering functionality into oxides, J. Am. Ceram. Soc.91, 2429 (2008)

2008
[21]

Stachiotti and M

M. Stachiotti and M. Sepliarsky, Toroidal ferroelectricity in PbTiO 3 nanoparticles, Phys. Rev. Lett.106, 137601 (2011)

2011
[22]

Pertsev, A

N. Pertsev, A. Zembilgotov, and A. Tagantsev, Effect of mechanical boundary conditions on phase diagrams of epitaxial ferroelectric thin films, Phys. Rev. Lett.80, 1988 (1998)

1988
[23]

Pertsev, A

N. Pertsev, A. Zembilgotov, and A. Tagantsev, Equilib- rium states and phase transitions in epitaxial ferroelectric thin films, Ferroelectrics223, 79 (1999)

1999
[24]

Bratkovsky and A

A. Bratkovsky and A. Levanyuk, Abrupt appearance of the domain pattern and fatigue of thin ferroelectric films, Phys. Rev. Lett.84, 3177 (2000)

2000
[25]

Stephanovich, I

V. Stephanovich, I. Luk’yanchuk, and M. Karkut, Do- main proximity and ferroelectric transition in ferro- 11 paraelectric superlattices, Ferroelectrics291, 169 (2003)

2003
[26]

Stephanovich, I

V. Stephanovich, I. Luk’yanchuk, and M. Karkut, Domain-enhanced interlayer coupling in ferroelec- tric/paraelectric superlattices, Phys. Rev. Lett.94, 047601 (2005)

2005
[27]

De Guerville, I

F. De Guerville, I. Luk’yanchuk, L. Lahoche, and M. El Marssi, Modeling of ferroelectric domains in thin films and superlattices, Mater. Sci. Eng. B120, 16 (2005)

2005
[28]

Baudry, I

L. Baudry, I. A. Luk’yanchuk, and A. Razumnaya, Dynamics of field-induced polarization reversal in thin strained perovskite ferroelectric films with c-oriented po- larization, Phys. Rev. B91, 144110 (2015)

2015
[29]

I. A. Luk’yanchuk, L. Lahoche, and A. Sen´ e, Universal properties of ferroelectric domains, Phys. Rev. Lett.102, 147601 (2009)

2009
[30]

Wang, S.-Q

J. Wang, S.-Q. Shi, L.-Q. Chen, Y. Li, and T.-Y. Zhang, Phase-field simulations of ferroelectric/ferroelastic polar- ization switching, Acta Mater.52, 749 (2004)

2004
[31]

De Gennes and J

P.-G. De Gennes and J. Prost,The physics of liquid crys- tals, 83 (Oxford university press, 1993)

1993
[32]

Pompe, X

W. Pompe, X. Gong, Z. Suo, and J. Speck, Elastic energy release due to domain formation in the strained epitaxy of ferroelectric and ferroelastic films, J. Appl. Phys.74, 6012 (1993)

1993
[33]

V. G. Koukhar, N. A. Pertsev, and R. Waser, Ther- modynamic theory of epitaxial ferroelectric thin films with dense domain structures, Phys. Rev. B64, 214103 (2001)

2001
[34]

S. Li, Y. Zhu, Y. L. Tang, Y. Liu, S. Zhang, Y. Wang, and X. Ma, Thickness-dependent a1/a2 domain evolution in ferroelectric PbTiO3 films, Acta Mater.131, 123 (2017)

2017
[35]

Kondovych, L

S. Kondovych, L. Boron, F. N. Di Rino, M. Sepliarsky, A. G. Razumnaya, A. Sen´ e, and I. A. Lukyanchuk, Surface-tension-induced phase transitions in freestanding ferroelectric thin films, Nano Lett.25, 12987 (2025)

2025
[36]

Logg, K.-A

A. Logg, K.-A. Mardal, G. N. Wells,et al.,Automated Solution of Differential Equations by the Finite Element Method(Springer, Berlin, 2012)

2012
[37]

Geuzaine and J.-F

C. Geuzaine and J.-F. Remacle, Gmsh: A three- dimensional finite element mesh generator with built-in pre- and post-processing facilities, Int. J. Numer. Meth- ods Eng.79, 1039 (2009)

2009
[38]

Janelli and R

A. Janelli and R. Fazio, Adaptive stiff solvers at low ac- curacy and complexity, J. Comput. Appl. Math.191, 246 (2006)

2006
[39]

Balay, S

S. Balay, S. Abhyankar, M. F. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, A. Dener, V. Ei- jkhout, W. D. Gropp, D. Karpeyev, D. Kaushik, M. G. Knepley, D. A. May, L. C. McInnes, R. T. Mills, T. Mun- son, K. Rupp, P. Sanan, B. F. Smith, S. Zampini, H. Zhang, and H. Zha...

2021
[40]

Balay, S

S. Balay, S. Abhyankar, M. F. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, A. Dener, V. Ei- jkhout, W. D. Gropp, D. Karpeyev, D. Kaushik, M. G. Knepley, D. A. May, L. C. McInnes, R. T. Mills, T. Mun- son, K. Rupp, P. Sanan, B. F. Smith, S. Zampini, H. Zhang, and H. Zha...

2021
[41]

Y.-T. Shao, S. Das, Z. Hong, R. Xu, S. Chandrika, F. G´ omez-Ortiz, P. Garc ´ ıa-Fern´ andez, L.-Q. Chen, H. Y. Hwang, J. Junquera,et al., Emergent chirality in a polar meron to skyrmion phase transition, Nat. Commun.14, 1355 (2023)

2023
[42]

Sheng, J

G. Sheng, J. Zhang, Y. Li, S. Choudhury, Q. Jia, Z. Liu, and L. Chen, Domain stability of PbTiO 3 thin films un- der anisotropic misfit strains: Phase-field simulations, J. Appl. Phys.104, 054105 (2008)

2008
[43]

Z. Hong, A. R. Damodaran, F. Xue, S.-L. Hsu, J. Brit- son, A. K. Yadav, C. T. Nelson, J.-J. Wang, J. F. Scott, L. W. Martin,et al., Stability of polar vortex lattice in ferroelectric superlattices, Nano Lett.17, 2246 (2017)

2017
[44]

C. Dai, Z. Hong, S. Das, Y.-L. Tang, L. W. Martin, R. Ramesh, and L.-Q. Chen, Strain effects on stabil- ity of topological ferroelectric polar configurations in (PbTiO3)n/(SrTiO3)n superlattices, Appl. Phys. Lett. 123, 052903 (2023)

2023
[45]

H.-M. Li, H. Zhang, Y.-J. Wang, Y.-L. Tang, Y.-L. Zhu, and X.-L. Ma, Misfit strain-misfit strain phase diagram of (110)-oriented ferroelectric PbTiO3 films: a phase-field study, Microstr.4, 2024004 (2024)

2024
[46]

P. Tong, L. Zhou, K. Du, M. Zhang, Y. Sun, T. Sun, Y. Wu, Y. Liu, H. Guo, Z. Hong,et al., Thermal trig- gering for multi-state switching of polar topologies, Nat. Phys.21, 464 (2025)

2025
[47]

Gregg, Murmurations of electric dipoles, Nat

J. Gregg, Murmurations of electric dipoles, Nat. Phys. 21, 344 (2025)

2025

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.

[1] [1]

The realized soft mode corresponds to the maximal critical temperature of the solution of (13), and the as- sociated wave vectork ∗ sets the modulation periodicity. In the spirit of Landau theory, we assume that close to the transition the spatial profile of the inhomogeneous structure,P(r), remains essentially unchanged, while its amplitude grows upon co...

work page

[2] [2]

Junquera, Y

J. Junquera, Y. Nahas, S. Prokhorenko, L. Bellaiche, J. ´I˜ niguez, D. G. Schlom, L.-Q. Chen, S. Salahuddin, D. A. Muller, L. W. Martin,et al., Topological phases in polar oxide nanostructures, Rev. Mod. Phys.95, 025001 (2023)

work page 2023

[3] [3]

S. Das, Z. Hong, M. McCarter, P. Shafer, Y.-T. Shao, D. Muller, L. Martin, and R. Ramesh, A new era in fer- roelectrics, APL Mater.8, 120902 (2020)

work page 2020

[4] [4]

I. A. Lukyanchuk, A. G. Razumnaya, S. Kondovych, Y. A. Tikhonov, B. Khesin, and V. M. Vinokur, Topo- logical foundations of ferroelectricity, Phys. Rep.1110, 1 (2025)

work page 2025

[5] [5]

Luk’yanchuk, A

I. Luk’yanchuk, A. Razumnaya, S. Kondovych, Y. Tikhonov, and V. M. Vinokur, Topological fer- roelectric chirality, Preprint arXiv:2406.19728 (2024)

arXiv 2024

[6] [6]

A. K. Yadav, K. X. Nguyen, Z. Hong, P. Garc ´ ıa- Fern´ andez, P. Aguado-Puente, C. T. Nelson, S. Das, B. Prasad, D. Kwon, S. Cheema,et al., Spatially re- solved steady-state negative capacitance, Nature565, 468 (2019)

work page 2019

[7] [7]

Yadav, C

A. Yadav, C. Nelson, S. Hsu, Z. Hong, J. Clarkson, C. Schlep¨ utz, A. Damodaran, P. Shafer, E. Arenholz, L. Dedon,et al., Observation of polar vortices in oxide superlattices, Nature530, 198 (2016)

work page 2016

[8] [8]

S. Das, Y. Tang, Z. Hong, M. Gon¸ calves, M. Mc- Carter, C. Klewe, K. Nguyen, F. G´ omez-Ortiz, P. Shafer, E. Arenholz,et al., Observation of room-temperature po- lar skyrmions, Nature568, 368 (2019)

work page 2019