{"id":"ce1c6054-5e63-4ef9-9941-970d39fecc90","arxiv_id":"2601.09160","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Closed-form extensions of Stoney's curvature equation for films with flexoelectric and piezoelectric coupling, in open/closed circuit and uniform/graded cases.","lead":"This paper derives new equations that predict how much a film-substrate plate bends when the film couples mechanical strain to electric fields through piezoelectricity or flexoelectricity. The formulas extend the classic 1909 Stoney equation, used to measure film stress from curvature, to electromechanically active films.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The open/closed circuit contrast rests on switching between two flexoelectric energy forms that differ by a boundary term; without a boundary-consistent single energy, the central curvature formulas are not established.","rationale":"The reader's weakest assumption correctly identifies the most load-bearing issue: the paper uses two different flexoelectric enthalpy densities for the two electrical boundary conditions, although these forms differ only by a boundary term. In a well-posed continuum theory, direct and converse flexoelectricity are two manifestations of one coupling; the choice of E vs ∇E in the energy is a gauge convention fixed by boundary terms, not by whether the circuit is open or closed. Since the paper omits those surface terms and does not show gauge equivalence, the central formulas—especially the claimed difference between open- and closed-circuit curvature and strain—are not established. This is more central than the secondary substrate-electrically-inert assumption because it directly affects the headline generalization and the open/closed comparison. The recommended verdict remains CONDITIONAL: the concern does not prove the results are wrong, but it shows the derivation is incomplete in a way that could change the formulas. A corrected derivation with a single energy and consistent boundary conditions might produce similar results, but that must be checked before quantitative metrology is based on Eqs. (51)–(54) and (94)–(97).","tokens_in":20955,"tokens_out":10507,"duration_ms":102655,"concrete_test":"Repeat the §4.1 closed-circuit derivation using one energy density for both circuits—e.g., use Eq. (5) for the closed-circuit optimization and include the proper boundary surface term at z=hs/2 and z=hs/2+hf, or equivalently solve the full Euler–Lagrange equations with fixed φ boundaries for the closed case and D·n=0 for the open case. Compare the resulting ε0 and κ to Eqs. (51)–(54). A minimal check: set e31=0, V≠0; Eq. (5) gives a nonzero flexoelectric torque proportional to μ12 V κ, while Eq. (40) gives zero. If the minimized curvature changes, the converse-side formulas are mis-modeled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim depends on treating Eqs. (5) and (6) as two independent material models, one for the open-circuit ('direct') case and one for the closed-circuit ('converse') case. They are not independent: ψ_conv − ψ_dir = μ_lijk(ε_ij E_l),k, a pure divergence. Integrated over the film, the difference is a surface term. The paper minimizes the volume enthalpy (49) without including this flexoelectric surface term, and never shows that the same bulk energy plus consistent boundary conditions reproduces both circuit cases. In a fixed-voltage problem, the boundary term is not generally zero; it is the flexoelectric surface work that must be added to (or subtracted from) the bulk enthalpy. A concrete failure appears in Case IV: for e31=0, the closed-circuit converse enthalpy (40) vanishes because E_z=V/hf is uniform, so the paper concludes the flexoelectric effect disappears in the closed circuit (Eq. 57). But using the single standard direct-form energy −μ ε' E, a bent film has a nonzero bulk term −2μ12 κ E_z even when E_z is uniform. The predicted open/closed contrast is therefore likely an artifact of the gauge choice rather than physics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives closed-form extensions of Stoney's equation for a film/substrate bilayer in which the film is piezoelectric and flexoelectric. Assuming Kirchhoff kinematics, transverse isotropy, and axial symmetry, it solves the one-dimensional electrostatic field for open- and closed-circuit conditions, integrates the enthalpy over the thickness, and minimizes with respect to stretching strain and curvature. Explicit formulas are given for uniform properties (Eqs. 51–54) and for a linear mismatch-strain gradient (Eqs. 94–97), followed by parametric plots and polarization expressions. The stated goal is to enable curvature-based extraction of film electromechanical properties.","tokens_in":21282,"tokens_out":20929,"duration_ms":231642,"significance":"The linear-elastic, axisymmetric framework is natural, and the paper is explicit about its variational principle, uses no fitted parameters, and provides closed-form expressions. If the closed-circuit branch is correct, the work is a practical extension of Stoney metrology to piezo/flexo films. However, the central open/closed contrast is tied to an unproven choice of enthalpy density; the most striking result—that flexoelectricity disappears in a closed circuit for e31=0, Eq. (57)—is not supported by a boundary-consistent calculation. The direct/open-circuit part and the elastic/piezoelectric reductions are plausible, but the closed-circuit formulas need revision.","major_comments":[{"comment":"The two flexoelectric enthalpy densities are related by ψ_conv − ψ_dir = (µ_ijlk ε_ij E_l),k, a pure divergence. The paper adopts ψ_conv for the closed circuit and ψ_dir for the open circuit, but minimizes the volume integral (49) without including the corresponding surface term. This is not invariant after eliminating the electrostatic field and minimizing over κ, because κ changes ε at the electrode surfaces. A concrete failure is Case IV: for e31=0 the closed-circuit solution (32) gives E_z=V/h_f, so ψ_conv=0 and Eq. (57) follows, whereas the paper's own direct form (5) gives −2µ12 ε_rr,z E_z = 2µ12 κ V/h_f ≠ 0 in the same state. Thus the disappearance of the flexoelectric effect in the closed circuit, and more generally the open/closed contrast in §4.1, is not established unless the missing surface work is included or shown to vanish. Please re-derive with a single energy and consist","section":"§2 and §4, Eqs. (5)–(6), (49), (57)"},{"comment":"The statement that the stress and electric displacement expressions hold irrespective of the choice of enthalpy density is true for the bulk constitutive laws, but it does not justify replacing the volume enthalpy in Eq. (49). Because the two forms differ by a boundary term, the reduced plate energy is changed by a κ-dependent surface contribution. This affects the closed-circuit formulas: Eq. (40) is proportional to e31, so it vanishes when e31→0 and leads to Eq. (57); using the direct form (5) instead gives a leading closed-circuit flexoelectric term proportional to µ12 V κ/h_f. The open/closed contrast is therefore a consequence of the chosen energy gauge rather than of the physics, unless the authors can prove that the boundary term cancels identically for the stated boundary conditions.","section":"§4, Eqs. (12)–(13), (40), (51)–(52)"}],"minor_comments":[{"comment":"The dielectric enthalpy density for the direct case is missing a factor: −(1/2)k33(E^Dir_z)^2 should contain the coefficient 2e31/k33, not e31. The integrated result in Eq. (46) appears consistent with the corrected density, but the displayed formula is misleading.","section":"Eq. (37)"},{"comment":"Eq. (61) uses k instead of κ for the curvature in the strain expression; Eq. (73) writes the substrate elastic enthalpy with M^f instead of M^s. These typos should be fixed.","section":"Eqs. (61), (73)"},{"comment":"The statement that the electric displacement in the substrate is zero because piezoelectricity and flexoelectricity are absent is not correct as written: for a dielectric substrate, D = kE. Either the substrate is intended to be conductive or the electric field is assumed to be screened; this assumption should be stated explicitly.","section":"§3, after Eq. (27)"},{"comment":"The open-circuit electric field (33) is derived for e31 ≠ 0 and contains factors of 1/e31. The e31 → 0 limit used in Case IV is well-defined for the integrated quantities, but this should be stated explicitly to avoid confusion about division by zero.","section":"§4, Eq. (33) and Case IV"},{"comment":"The final formulas are very long and no intermediate algebra or symbolic verification is provided. Given the typos in nearby equations, the authors should supply a supplementary derivation or a machine-checked symbolic verification of these expressions.","section":"§4.1 and §5.1, Eqs. (51)–(54), (94)–(97)"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern lands: the closed-circuit branch rests on the choice of Eq. (6), and Eq. (57) is the clearest symptom. I do not see misconduct or a novelty problem; the issue is a fixable modeling inconsistency. If the author can justify ψ_conv as the correct fixed-voltage enthalpy including surface terms, or re-derive the closed-circuit formulas from a boundary-consistent single-energy principle, the paper could become acceptable. Otherwise the open/closed contrast and all closed-circuit formulas should be revised."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper extends Freund's Stoney-type analysis to films with piezoelectric and flexoelectric coupling, and that extension is real. The closed-form expressions for curvature and stretching strain in the uniform and graded-mismatch cases (Eqs. 51-54 and 94-97) are not in the cited literature, and the derivation from a variational principle is transparent. The parameter studies are thorough. If the central results hold, this is a useful tool for thin-film metrology.\n\nThe central results do not all hold in their present form. The author uses two flexoelectric enthalpy densities, direct (Eq. 5) and converse (Eq. 6), and assigns one to the open-circuit case and the other to the closed-circuit case. These two expressions differ by a pure divergence, so they are not independent material models; they differ by a surface term. The paper minimizes the volume enthalpy without including that surface term. That is exactly the step that needs justification, and it is the load-bearing step.\n\nThe red flag is Case IV. With e31=0, the closed-circuit electric field is uniform, so the converse enthalpy (Eq. 40) vanishes and the paper concludes that flexoelectricity has no effect on curvature in the closed circuit (Eq. 57). But the same physics, written with the direct form, has a nonzero bulk term proportional to mu12*kappa*Ez even when Ez is uniform. The open/closed contrast is therefore likely an artifact of the gauge choice, not a real effect. This is not a cosmetic issue; it compromises the closed-circuit formulas (51), (52), (94), (95), and the discussion that depends on them.\n\nThere are also several typos (Eqs. 37, 61, 73) that make verification harder. None of these are fatal by themselves, but they add to the sense that the algebra has not been checked symbolically.\n\nThe open-circuit results and the overall framework are worth taking seriously. The fix is to redo the closed-circuit derivation with a single flexoelectric energy and the proper boundary/surface contribution—or to explicitly justify why the two volume forms are each valid under the stated electrical boundary condition. Without that, the quantitative formulas should not be used for metrology.\n\nThis is exactly the kind of paper that deserves peer review: important enough, novel enough, and the flaw is identifiable and fixable. I would send it to a competent referee with instructions to check the boundary-term consistency and to verify the algebra in a computer algebra system.","headline":"A genuinely new Stoney extension for flexoelectric films, but the closed-circuit formulas rest on a questionable swap between two energy forms that differ by a boundary term.","tokens_in":21643,"tokens_out":4621,"would_cite":false,"duration_ms":50383,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives explicit closed-form generalizations of Stoney's equation that predict curvature and stretching strain of a piezoelectric-flexoelectric thin film on an elastic substrate, opening curvature metrology to electromechanical f","keywords":["flexoelectricity","piezoelectricity","Stoney's equation","thin films","substrate curvature","electromechanical coupling","strain gradient","film-substrate bilayer"],"falsifier":"A finite-element solution of the full coupled equilibrium equations using one flexocoupling tensor and explicit surface terms—rather than the paper's two separate enthalpy densities—computed for the same geometry and constants as Figures 3 to 5 would settle the matter: if its predicted stretching strain and curvature do not match Eqs. (51)–(54) in both open and closed circuits, the two-enthalpy treatment is not equivalent.","tokens_in":20835,"feed_emoji":"📐","tokens_out":11022,"duration_ms":480215,"temperature":0.7,"pith_summary":"This paper generalizes Stoney's classical relation between film stress and substrate curvature to thin films that are both piezoelectric and flexoelectric. It derives closed-form formulas for the midplane stretching strain and curvature of a film-substrate bilayer, for uniform and linearly graded film properties, in both open-circuit (direct flexoelectric) and closed-circuit (applied-voltage, converse) configurations. If correct, these formulas let experimenters extract electromechanical film constants from curvature measurements in the same way residual stress is measured today, and they predict how applied voltage bends the bilayer. A notable conclusion is that the electric polarization in the film varies linearly through the thickness, and that Stoney's original estimates increasingly overestimate stretching strain and curvature as the film becomes thicker and stiffer relative to the substrate.","feed_headline":"Stoney's equation gains piezo- and flexoelectric terms","feed_subtitle":"New formulas link substrate curvature to film electromechanical constants in open and closed circuits.","key_machinery":"The carrying mechanism is a total enthalpy functional for the film–substrate bilayer: elastic, dielectric, piezoelectric, and flexoelectric energy densities integrated over volume, with flexoelectricity represented in two forms—a direct form coupling strain gradient to electric field (voltage from deformation) and a converse form coupling strain to electric-field gradient (deformation from voltage). The argument proceeds by requiring this functional to be stationary with respect to the midplane stretching strain ε0 and the curvature κ, which yields two coupled algebraic equations whose solution is the closed-form generalization. In the non-uniform case, the same stationarity machinery is app","core_discovery":"The paper's central claim is that electromechanical coupling in a thin film can be incorporated into explicit, energy-minimized expressions for the bilayer's midplane stretching strain ε0 and curvature κ. For uniform film properties it reports closed-form formulas (Eqs. 51–54) that connect ε0 and κ to the film's piezoelectric constant e31, flexocoupling constant µ12, permittivity k33, elastic moduli, film and substrate thicknesses, mismatch strain, and applied voltage, with separate branches for closed-circuit (converse) and open-circuit (direct) conditions; a linearly varying mismatch strain leads to the analogous formulas Eqs. (94)–(97). The paper shows these expressions reduce to the pure","pith_inferences":["The difference between open- and closed-circuit curvature of the same sample could be used to identify e31 and µ12 separately without knowing the residual stress, because voltage and boundary conditions enter the formulas through different terms—a protocol the paper motivates but does not spell out.","Because the direct and converse flexoelectric enthalpy densities differ only by an integration-by-parts boundary term, the predicted open/closed contrast should be cross-checked against a formulation with one flexoelectric tensor plus consistent surface terms before taking the quantitative difference at face value.","The same energy-minimization framework could accommodate finite deformations, anisotropic or inhomogeneous substrates, and non-axisymmetric buckling, with the present formulas serving as the linear benchmark.","The linear-through-thickness polarization prediction could be tested directly on ultrathin films via surface-potential or second-harmonic measurements, which would also test the constitutive assumption that flexoelectricity couples to strain gradient rather than to strain itself."],"forward_implications":["Curvature metrology can be extended from residual stress to electromechanical properties: measuring κ and ε0 in open- and closed-circuit configurations gives access to e31 and µ12 as well as the mismatch strain.","An applied voltage bends the bilayer, with the bending direction and magnitude set by the film/substrate thickness and stiffness ratios—so the formulas double as a design rule for flexoelectric actuators on compliant substrates.","In the closed-circuit/converse case with no piezoelectricity, flexoelectricity alone leaves the curvature and stretching strain unchanged; the open-circuit/direct case does change them, so comparing the two configurations separates flexoelectric from piezoelectric effects.","For films with a linear mismatch gradient, the derived expressions let a single curvature/strain measurement be decomposed into average and gradient contributions to the film stress.","The polarization is linear through the film thickness and tied to both ε0 and κ, giving an independent surface measurement that can corroborate the curvature-derived constants."],"fun_headline_variants":["Stoney's equation gets a flexo-piezo upgrade","Circuit matters: Stoney's equation for flexoelectric films","Beyond Stoney: electromechanical curvature formulas","Stoney's law now bends for piezo- and flexoelectric films","New Stoney equations for open and closed circuits in thin films"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that direct and converse flexoelectricity can be modeled by two distinct enthalpy densities that differ by an integration-by-parts boundary term; if a single flexoelectric tensor together with the correct boundary conditions changes the open- and closed-circuit curvature predictions, the paper's central contrast between the two configurations would not survive.","fun_headline_variants_meta":{"raw":{"variants":["Stoney's equation gets a flexo-piezo upgrade","Circuit matters: Stoney's equation for flexoelectric films","Beyond Stoney: electromechanical curvature formulas","Stoney's law now bends for piezo- and flexoelectric films","New Stoney equations for open and closed circuits in thin films"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000674,"raw_usage":{"total_tokens":2860,"prompt_tokens":656,"completion_tokens":2204,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":400,"completion_tokens_details":{"reasoning_tokens":2120}},"tokens_in":400,"tokens_out":2204,"duration_ms":19707,"temperature":1.0,"reasoning_tokens":2120,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T06:21:45.656158+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A finite-element solution of the full coupled equilibrium equations using one flexocoupling tensor and explicit surface terms—rather than the paper's two separate enthalpy densities—computed for the same geometry and constants as Figures 3 to 5 would settle the matter: if its predicted stretching strain and curvature do not match Eqs. (51)–(54) in both open and closed circuits, the two-enthalpy treatment is not equivalent.","supporting_citations":[],"review_version":1}