{"id":"7f5e565c-95fd-473c-b20d-6d9ebf033c25","arxiv_id":"2509.06508","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A single-mode variational ansatz yields phase diagrams of vortex, helix, wave, and uniform polarization states in strained PbTiO3 films as functions of temperature, strain, and thickness.","lead":"This paper develops an analytical approximation for the complex swirled polarization patterns that form in strained ferroelectric films, and maps out where each pattern is stable. It offers a fast alternative to expensive computer simulations for guiding experiments on ferroelectric nanostructures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-mode ansatz truncation: phase diagrams not validated against multi-mode equilibrium states; h' region already contradicted by paper's own phase-field result.","rationale":"The reader identified the single-harmonic ansatz as the weakest assumption, and I agree. The paper's phase diagrams are variational upper bounds within a restricted ansatz; if the true minimum lies outside this class, the phase boundaries and even phase ordering could change. The paper provides only three sparse phase-field validation points, and at u_m≈0 the simulation already produces a labyrinthine wh-state rather than the predicted h' helix, which is a concrete internal indicator that the truncation matters. However, this does not by itself prove the central claim is false; it motivates a targeted computational check. Running a phase-field simulation at a point where h' is predicted would settle whether the single-mode state is actually stable. Since the paper is otherwise a reasonable exploratory variational study and the concern is about quantitative predictive power rather than internal inconsistency, the CONDITIONAL verdict remains appropriate. My read does not change the reader's verdict.","tokens_in":14545,"tokens_out":8429,"duration_ms":90093,"concrete_test":"At a point solidly inside the predicted h' region (e.g., u_m=+0.002, 2a_f=12 nm, T=300 K), run a phase-field simulation (as in Sec. IV B) starting from both a random initial condition and the analytic h' ansatz (14). Relax to equilibrium and compare the resulting texture and free energy. If the h' ansatz is not a local minimum, or if a multi-mode/mixed state has lower free energy, the h' phase is not the true equilibrium and the phase diagram needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The phase diagrams in Fig. 3 are computed by minimizing the GLD free energy within the single-wave-vector ansatz (14), which fixes a 1D sinusoidal texture in (x,z) and enforces div P=0. This ansatz excludes higher harmonics and superpositions of wave vectors, yet the paper itself notes in Sec. V.A that bubble/skyrmion states require a star of wave vectors. More directly, the phase-field simulation in Sec. IV.B at u_m≈0 (Fig. 4b) does not relax to the predicted h' helix but to a labyrinthine wh-type texture, indicating that the excluded modes are competitive exactly where the phase diagram predicts h'. Thus the central claim that v' replaces the uniform c-phase and h' appears near small tensile strain is established only within a restricted variational class; the true GLD equilibrium could have different phase boundaries or phase ordering.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14758,"tokens_out":7682,"duration_ms":85129,"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":[{"comment":"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.","section":"Sec. IV.B (Fig. 4b) vs Sec. IV.A (Fig. 3a,c)"},{"comment":"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.","section":"Sec. III.A–C, Eq. (14)"},{"comment":"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.","section":"Sec. IV.B (Fig. 4c) and Sec. III.B/IV.A (aa phase)"}],"minor_comments":[{"comment":"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.","section":"Sec. II.C"},{"comment":"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.","section":"Fig. 4b caption"},{"comment":"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.","section":"Sec. IV.A"},{"comment":"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.","section":"Sec. V.A"},{"comment":"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é').","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper has clear merit as a variational tool, but the central phase-diagram claim for h' is directly contradicted by the authors' own full-GLD phase-field simulation at u_m ≈ 0. This is not a matter of style; it is a load-bearing discrepancy that must be addressed by either providing evidence that the single-mode ansatz is sufficient (e.g., multi-mode convergence or targeted phase-field tests in each phase region) or by explicitly reframing the phase diagrams as idealized variational upper bounds. I therefore recommend major revision rather than rejection. The paper's reliance on the authors' earlier soft-domain papers is not circular, since the GLD coefficients are taken from the literature and not fitted to the target phases."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nThe 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.\n\nWho 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.","headline":"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.","tokens_in":665,"tokens_out":897,"would_cite":true,"duration_ms":36624,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["77.80.-e"],"model":"deepseek-v4-flash","headline":"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","keywords":["ferroelectric thin films","topological polarization textures","vortex phase","soft-domain approach","Ginzburg-Landau-Devonshire","epitaxial strain","PbTiO3","phase diagrams"],"falsifier":"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.","tokens_in":14489,"feed_emoji":"🌀","tokens_out":5867,"duration_ms":64054,"temperature":0.7,"pith_summary":"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.","feed_headline":"Strain tilts ferroelectrics into vortex and helix phases","feed_subtitle":"A compact analytic model maps how temperature, strain, and thickness choose between uniform and swirling polarization in PbTiO3 films.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the strain-renormalized GLD coefficients for PbTiO3 and the uniform-phase baseline that nonuniform states must compete with.","marker":"[20]"},{"why":"Establishes the soft-domain assumption that the spatial profile of the inhomogeneous structure stays fixed while amplitude grows on cooling.","marker":"[27]"},{"why":"Introduces the soft-domain expansion in ferroelectric superlattices that this work extends.","marker":"[23]"},{"why":"Applies soft-domain modeling to ferroelectric domains in thin films and superlattices, grounding the ansatz.","marker":"[25]"},{"why":"Provides the gradient stiffness coefficients for PbTiO3 used in the gradient-energy term.","marker":"[28]"},{"why":"Early prediction of flux-closure vortices, the class of textures the ansatz is designed to describe.","marker":"[16]"}],"fun_headline_variants":["Strain twists ferroelectric films into vortex and helix phases","How strain sculpts topological vortices in ferroelectric films","Strain spins polarization into vortex states in ferroelectrics","Compact theory reveals vortex patterns in strained PbTiO3 films","Vortex and helix phases emerge in strained ferroelectric films"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Strain twists ferroelectric films into vortex and helix phases","How strain sculpts topological vortices in ferroelectric films","Strain spins polarization into vortex states in ferroelectrics","Compact theory reveals vortex patterns in strained PbTiO3 films","Vortex and helix phases emerge in strained ferroelectric films"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000855,"raw_usage":{"total_tokens":3495,"prompt_tokens":630,"completion_tokens":2865,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":374,"completion_tokens_details":{"reasoning_tokens":2795}},"tokens_in":374,"tokens_out":2865,"duration_ms":27203,"temperature":1.0,"reasoning_tokens":2795,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T23:28:50.897054+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the strain-renormalized GLD coefficients for PbTiO3 and the uniform-phase baseline that nonuniform states must compete with."},{"cited_title":"De Guerville, I","cited_arxiv_id":null,"evidence_quote":"Establishes the soft-domain assumption that the spatial profile of the inhomogeneous structure stays fixed while amplitude grows on cooling."},{"cited_title":"Pertsev, A","cited_arxiv_id":null,"evidence_quote":"Introduces the soft-domain expansion in ferroelectric superlattices that this work extends."},{"cited_title":"Stephanovich, I","cited_arxiv_id":null,"evidence_quote":"Applies soft-domain modeling to ferroelectric domains in thin films and superlattices, grounding the ansatz."},{"cited_title":"Baudry, I","cited_arxiv_id":null,"evidence_quote":"Provides the gradient stiffness coefficients for PbTiO3 used in the gradient-energy term."},{"cited_title":"Nahas, S","cited_arxiv_id":null,"evidence_quote":"Early prediction of flux-closure vortices, the class of textures the ansatz is designed to describe."}],"review_version":1}