{"id":"bffe79ed-381b-4f60-bb44-3a28724844dd","arxiv_id":"2505.15714","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"In a single material showing both ordinary and ferroelectric nematic phases, splay elasticity stiffens sharply in the polar phase while twist elasticity softens, consistent with electrostatic stiffening by polarization charges.","lead":"This paper measures how hard it is to bend the alignment of a ferroelectric nematic liquid crystal, comparing the nonpolar and polar phases of the same material. It finds that the polar phase resists splay deformation much more strongly while allowing twist deformation more easily, a behavior attributed to electric charges created by the molecular alignment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Twist softening may be an artifact of the Haller-extrapolated chi_m,a, which the paper's own birefringence data show is invalid in the NF phase.","rationale":"The reader's weakest assumption—that chi_m,a is not measured in NF but extrapolated via a Haller fit—is exactly the load-bearing point. The stress-test sharpens it: the paper's own birefringence data contradict the extrapolation, because Delta n deviates upward from the Haller trend in the M phase and jumps at the M-NF transition (Fig. 2 and Section 3.1). Consequently, the extrapolated chi_m,a in NF is likely an underestimate, which directly scales K11 and K22 downward. The twist softening, a key new result, could disappear once a correct chi_m,a is used. The splay stiffening, however, would only become stronger, so the splay part of the claim is robust. The reader's CONDITIONAL verdict is appropriate: the quantitative claims and especially the twist softening need direct verification of chi_m,a or an explicit uncertainty analysis. My assessment does not change the verdict but reinforces the condition. I agree with the reader that this is the weakest assumption; the additional internal inconsistency from the birefringence data makes the concern more concrete. The paper's electrostatic estimate also contains a numerical inconsistency (400 pN from P0 = 6 uC/cm^2, epsilon = 100, lambda_D = 100 nm actually gives about 40 nN; lambda_D = 10 nm gives 400 pN), but this is secondary because the experimental stiffening observation does not rely on the estimate. The central concern remains the experimental determination of chi_m,a in NF.","tokens_in":15100,"tokens_out":8724,"duration_ms":74708,"concrete_test":"Recompute K11 and K22 in the NF phase using chi_m,a rescaled by the measured birefringence: in the N phase calibrate the relation between chi_m,a and Delta n/Delta n0 using Fig. 8b and Fig. 2, then extend this calibration to NF temperatures where Delta n is measured, and recompute K_i from the thresholds in Fig. S2b using K_i = (Bc d/pi)^2 chi_m,a/mu0. If the NF K22 no longer lies below the N-phase trend, or the NF K11 increase changes materially, the Haller extrapolation is the cause. A direct torque-magnetometer measurement of chi_m,a on a uniformly aligned NF sample would settle the question definitively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing concern is the determination of chi_m,a in the NF phase. In Section 3.4, K11 and K22 in NF are computed from magnetic Fréedericksz thresholds in Fig. S2b using chi_m,a extrapolated from the N phase with a Haller fit (chi_m,a = chi_m,a0(1 - T/T*)^beta, with chi_m,a0 = (6 +/- 1) x 10^-6, T* = 359.8 K, beta = 0.40). The paper states that Haller-type behaviour of the orientational order parameter holds across N and NF, but Fig. 2 and the accompanying text show exactly the opposite: birefringence deviates upward in the M phase and jumps at the M-NF transition, so S(T) in NF is higher than the Haller extrapolation. Because chi_m,a is proportional to S at leading order, the extrapolated chi_m,a underestimates the true NF value. Since K_i = (Bc d/pi)^2 chi_m,a/mu0, all NF elastic constants are proportionally underestimated. The reported twist softening is therefore suspect: with a larger, physically justified chi_m,a, K22 in NF moves upward and the softening may vanish. The splay stiffening, by contrast, would only increase, so that part of the claim is more robust. The novelty of the paper is partly the twist softening; if that is an artifact, the central claim is reduced.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a multi-technique study of a liquid crystal mixture exhibiting N, M (antiferroelectric), and NF (ferroelectric nematic) phases, with the goal of comparing the Frank elastic constants of the nonpolar and polar nematic phases. Birefringence, second-harmonic generation, broadband dielectric spectroscopy, and electric/magnetic Fréedericksz transition measurements are combined. The central claims are that the splay elastic constant K11 increases by nearly an order of magnitude in the NF phase, that the twist elastic constant K22 softens significantly in the NF phase, and that the splay stiffening can be attributed to the electrostatic energy of polarization splay, described in the supplementary information by an effective splay constant Keff(k) = K1 + P0^2/(εε0(k^2+κ^2)). The NF-phase elastic constants are obtained from magnetic Fréedericksz thresholds using a diamagnetic anisotropy χm,a extrapolated from the N phase via a Haller-type fit.","tokens_in":15506,"tokens_out":5779,"duration_ms":52872,"significance":"If the results are correct, the paper provides valuable quantitative information on how ferroelectric polar order modifies the mechanical response of a fluid nematic, and it tests a comparatively simple electrostatic model for splay stiffening. The strengths include the use of complementary experimental techniques, the SHG confirmation of polar order, a self-contained derivation in supplementary S3 that uses externally specified parameters rather than fitting the measured elastic constants, and the explicit recognition of the difficulty of measuring χm,a directly in the NF phase. However, the twist-softening claim rests on an extrapolated diamagnetic anisotropy that is contradicted by the paper's own birefringence data, and the absence of error bars on the NF elastic constants makes it difficult to assess the significance of the reported softening. The central splay-stiffening claim is more robust to the extrapolation issue, but the electrostatic estimate contains an apparent arithmetic inconsistency. These issues require substantive revision.","major_comments":[{"comment":"The Haller extrapolation of χm,a into the NF phase is load-bearing. The text states that 'Haller-type behaviour of the orientational order parameter across the whole range of N and NF phases was established by birefringence measurements,' but Fig. 2 shows the opposite: Δn deviates upward in the M phase and jumps at the M–NF transition, so the orientational order parameter in NF lies above the Haller continuation. Since χm,a is proportional to the orientational order parameter at leading order, the extrapolated χm,a is likely an underestimate of the true NF value. Every NF elastic constant Ki = (Bc d/π)^2 χm,a/μ0 is then proportionally underestimated. In particular, the twist softening in Fig. 9b may be an artifact: with a larger, physically motivated χm,a, K22 in the NF phase moves upward and the softening could vanish. The splay stiffening would survive this correction, but the twist claim needs either a direct measurement of χm,a in the NF phase or a sensitivity analysis over the plausible range of χm,a values.","section":"Section 3.4, Fig. 8b and Fig. 2"},{"comment":"No error bars or uncertainty propagation are provided for the elastic constants in the NF phase. The Haller parameters are quoted with uncertainties, e.g., χm,a0 = (6 ± 1) × 10^-6, and the threshold fields in Fig. S2b also carry measurement uncertainty. These propagate directly into K. Without confidence intervals it is not possible to judge whether the reported twist softening is statistically significant, which is essential because the softening is the less robust part of the central claim.","section":"Section 3.4, Fig. 9"},{"comment":"The numerical estimate for the electrostatic contribution is arithmetically inconsistent as written. The paper states that P0 = 6 μC cm^-2, ε = 100, and λD ≈ 100 nm give a correction of 400 pN. Direct evaluation of P0^2 λD^2/(εε0) gives approximately 4 × 10^-8 N = 4 × 10^4 pN, a factor of 100 larger than the quoted value. If a different choice of parameters is intended, that must be stated explicitly; as written, the claim that the electrostatic correction is of 'similar order of magnitude' to the observed stiffening is not supported by the formula and parameters given.","section":"Section 3.4, electrostatic estimate"}],"minor_comments":[{"comment":"The exponent N in the conductive term σDC/(iωε0)^N is not defined in the text; please specify its range and role in the fitting.","section":"Section 2, Eq. (1)"},{"comment":"The text says the splay constant 'exhibits a sharp increase by nearly an order of magnitude,' but the actual K11 values and the temperatures at which they are compared are not quoted. Please state these values so the reader can verify the magnitude of the effect.","section":"Section 3.4, Fig. 9"},{"comment":"The sentence 'A limitation of magnetic field measurements, however, is the need for accurate knowledge of the diamagnetic anisotropy, which is often challenging to determine directly' is important, but the subsequent discussion does not explain how the uncertainty in this extrapolation affects the central claims. A brief quantitative sensitivity statement would be helpful.","section":"Section 3.4, first paragraph"},{"comment":"The phase boundaries (N–M and M–NF) are not marked in Fig. S2b; adding vertical lines would make it easier to identify where the extrapolated χm,a is being used.","section":"Supplementary S2, Fig. S2b"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and important question in ferroelectric nematic research, and the experimental effort is substantial. The main concern is the Haller extrapolation of χm,a into the NF phase, which is contradicted by the paper's own birefringence data; this directly affects the twist-softening claim. The arithmetic error in the electrostatic estimate should also be corrected. I believe the manuscript can be revised to address these points, but the revision needs to include either additional measurements or a carefully argued sensitivity analysis rather than only a change in presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The twist softening is the weak link. It rests on a Haller extrapolation of the diamagnetic anisotropy into the NF phase, and the paper's own birefringence data show that extrapolation is invalid. The splay stiffening, by contrast, survives.\n\nWhat's genuinely new: the same material with N, M, and NF phases, and elastic constants measured in both the nonpolar and ferroelectric nematic phases. The N-phase determination of chi_m,a by combining electric and magnetic Fréedericksz thresholds is clean and carefully done. The electrostatic splay-stiffening calculation in S3 is self-contained, uses external parameter values, and is not fitted to the measured K's—so the central mechanism for splay stiffening is on solid ground.\n\nThe soft spot is load-bearing. Section 3.4 claims that Haller-type behavior of the orientational order parameter holds across N and NF, but Section 3.1 reports birefringence deviating upward in the M phase and jumping at the M–NF transition. Since birefringence tracks S, and chi_m,a is proportional to S, the extrapolated chi_m,a underestimates the true NF value. That means all NF elastic constants are proportionally underestimated. The reported twist softening could vanish once a physically justified chi_m,a is used; the splay stiffening would only increase. There are also no error bars on the K values in Fig. 9, and the electrostatic estimate uses assumed values for P0, epsilon, and lambda_D, giving only order-of-magnitude agreement.\n\nThis deserves serious refereeing—the experimental dataset is valuable and the splay stiffening is robust—but the twist-softening claim needs to be either backed by a direct measurement of chi_m,a in NF or explicitly downgraded to a tentative observation. The authors should also address the internal contradiction about Haller behavior and provide uncertainty analysis for the elastic constants.","headline":"Splay stiffening is real, but the twist softening is likely an artifact of the chi_m,a extrapolation that the paper's own birefringence data contradict.","tokens_in":16039,"tokens_out":3293,"would_cite":false,"duration_ms":28598,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In a material that shows both an ordinary nematic and a ferroelectric nematic phase, magnetic-field thresholds reveal that polar order raises the splay elastic constant by nearly an order of magnitude while softening twist.","keywords":["ferroelectric nematic","elastic constants","splay stiffening","twist softening","Fréedericksz transition","electrostatic screening","polarization bound charge","diamagnetic anisotropy"],"falsifier":"Measure the diamagnetic anisotropy directly in the $N_F$ phase by an independent method (for example torque magnetometry or a geometry that isolates the magnetic torque) and recompute $K_{11}$ and $K_{22}$ from the measured Fréedericksz thresholds; if the corrected values no longer show the near-order-of-magnitude splay increase and the twist decrease, the paper's central mechanical claim is refuted.","tokens_in":14934,"feed_emoji":"🧲","tokens_out":7408,"duration_ms":62864,"temperature":0.7,"pith_summary":"The paper asks how the mechanical elasticity of a liquid crystal changes when it acquires true ferroelectric order, and it answers by measuring the same material in its nonpolar nematic (N) and ferroelectric nematic ($N_F$) phases. Using magnetic Fréedericksz transitions, it finds that entering the $N_F$ phase raises the splay elastic constant $K_{11}$ by nearly an order of magnitude while significantly lowering the twist constant $K_{22}$. The authors trace the splay stiffening to the electrostatic energy of the bound polarization charges that any splay deformation creates, and they connect the twist softening to competing elastic and electrostatic torques that favor twisted director configurations. These results matter because ferroelectric nematics are the first three-dimensional ferroelectric fluids, and knowing how their mechanical response is shaped by electrostatic interactions is central to understanding their textures, defects, and potential uses in fast electro-optic devices.","feed_headline":"Splay rigidity jumps tenfold; twist softens in a polar nematic","feed_subtitle":"Magnetic-field thresholds across the N–NF transition show electrostatic splay stiffening in a room-temperature ferroelectric fluid.","key_machinery":"The load-bearing object is the effective wavevector-dependent splay elastic constant $K_{\\rm eff}(k)=K_1+P_0^2/(\\varepsilon\\varepsilon_0(k^2+\\kappa^2))$, derived in the supplementary information from the screened-Coulomb free energy of bound charges $\\rho=-\\nabla\\cdot\\mathbf{P}$. In the long-wavelength limit it becomes $K_1+P_0^2\\lambda_D^2/(\\varepsilon\\varepsilon_0)$, which for typical values ($P_0=6\\ \\mu\\mathrm{C\\,cm}^{-2}$, $\\varepsilon=100$, $\\lambda_D\\approx100\\ \\mathrm{nm}$) gives a correction of order 400 pN. The measurement machinery is the magnetic Fréedericksz transition: threshold fields $B_c=(\\pi/d)\\sqrt{\\mu_0 K_i/\\chi_{m,a}}$ in splay and twist geometries give $K_{11}$ and $K_{22}$, with the diamagnetic anisotropy $\\chi_{m,a}$ determined in the N phase by combining electric and magnetic thresholds and then extrapolated into the $N_F$ phase by a Haller fit.","core_discovery":"The central claim is that polar order reverses the mechanical hierarchy of the nematic state: in the $N_F$ phase the splay constant grows sharply, by nearly an order of magnitude compared with the adjacent N phase, while the twist constant softens markedly. The evidence comes from the magnetic Fréedericksz transition, where the critical field for splay and twist reorientation of the director is measured optically in planar cells; the same geometry yields the N-phase constants, so the comparison is made on one material across its N, intermediate, and $N_F$ phases. The splay jump is attributed to the electrostatic cost of polarization splay: with $\\mathbf{P}=P_0\\mathbf{n}$, a splay deformation creates bound charge, and screened Coulomb repulsion between those charges adds a wavevector-dependent term to the elastic energy, giving an effective splay constant $K_{\\rm eff}(k)=K_1+P_0^2/(\\varepsilon\\varepsilon_0(k^2+\\kappa^2))$. The twist softening is interpreted through the idea that electrostatic interactions in a polar fluid favor ambidextrous twist deformations, so the Frank twist term is effectively reduced.","pith_inferences":["The extrapolated diamagnetic anisotropy is the main quantitative uncertainty; an independent measurement in the $N_F$ phase could shift both constants, though the qualitative stiffening would survive unless the anisotropy changes by a large factor.","The electrostatic formula suggests a direct experiment the paper does not report: doping the material with an ionic additive should shorten the Debye length and continuously tune $K_{\\rm eff}$, providing a separate check of the mechanism.","Twist softening implies that weak-anchoring or confined $N_F$ samples might spontaneously develop twisted or chiral director fields, connecting this measurement to the helical polar phases mentioned in the introduction.","The wavevector dependence of $K_{\\rm eff}$ means that Fréedericksz thresholds probe only its long-wavelength value; short-wavelength distortions such as defect cores should experience a much larger effective stiffness, which could explain the evolution of the striped textures observed near the transition."],"forward_implications":["Splay deformations in the $N_F$ phase become much more expensive than in the N phase, so polar-aligned cells should resist splay distortions and favor configurations that avoid director divergence.","Twist reorientation becomes easier in the $N_F$ phase, so twist Fréedericksz transitions should occur at lower magnetic fields and twisted textures should appear more readily.","The electrostatic contribution to $K_{\\rm eff}$ depends on the Debye screening length, so ionic content and impurity concentration should measurably alter the apparent splay rigidity.","Because the effective splay constant is wavevector-dependent, the stiffening is strongest at short wavelengths, which bears on the formation of striped textures, conics, and other small-scale director structures in ferroelectric nematics."],"supporting_citations":[{"why":"Supplies the earlier magnetic-Fréedericksz measurement of splay and twist constants in a ferroelectric nematic that this paper extends to a material with an intervening antiferroelectric phase.","marker":"[24]"},{"why":"Provides the flexoelectric and electrostatic free-energy treatment that underlies the effective splay constant formula.","marker":"[42]"},{"why":"Documents conic textures in ferroelectric nematics attributed to splay stiffening, cited as independent evidence for the effect.","marker":"[21]"},{"why":"Predicts helical twist configurations from dipole-dipole interactions, cited to explain the twist softening in the $N_F$ phase.","marker":"[45]"},{"why":"Reports fluid superscreening and polarization following in confined ferroelectric nematics, supporting the screened-electrostatics picture.","marker":"[23]"},{"why":"Provides out-of-plane reorientation measurements in ferroelectric nematics used as a comparison for the magneto-optical threshold behavior.","marker":"[44]"}],"fun_headline_variants":["Splay stiffness jumps tenfold; twist softens in ferroelectric nematic","Polar order flips nematic elasticity: splay stiff, twist soft","Electrostatic splay stiffening, twist softening in polar nematic","Splay up, twist down: ferroelectric nematic mechanics","Nematic elasticity reversed: splay hard, twist weak in polar phase"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on a diamagnetic anisotropy that is extrapolated from the nonpolar nematic phase into the ferroelectric phase rather than measured there, so if polar order changes how strongly the molecules respond to a magnetic field, both the reported splay stiffening and twist softening would be systematically wrong.","fun_headline_variants_meta":{"raw":{"variants":["Splay stiffness jumps tenfold; twist softens in ferroelectric nematic","Polar order flips nematic elasticity: splay stiff, twist soft","Electrostatic splay stiffening, twist softening in polar nematic","Splay up, twist down: ferroelectric nematic mechanics","Nematic elasticity reversed: splay hard, twist weak in polar phase"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000761,"raw_usage":{"total_tokens":3337,"prompt_tokens":861,"completion_tokens":2476,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":477,"completion_tokens_details":{"reasoning_tokens":2378}},"tokens_in":477,"tokens_out":2476,"duration_ms":16529,"temperature":1.0,"reasoning_tokens":2378,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:12:12.230418+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the diamagnetic anisotropy directly in the $N_F$ phase by an independent method (for example torque magnetometry or a geometry that isolates the magnetic torque) and recompute $K_{11}$ and $K_{22}$ from the measured Fréedericksz thresholds; if the corrected values no longer show the near-order-of-magnitude splay increase and the twist decrease, the paper's central mechanical claim is refuted.","supporting_citations":[{"cited_title":"Zavvou, M","cited_arxiv_id":null,"evidence_quote":"Supplies the earlier magnetic-Fréedericksz measurement of splay and twist constants in a ferroelectric nematic that this paper extends to a material with an intervening antiferroelectric phase."},{"cited_title":"Kumari, B","cited_arxiv_id":null,"evidence_quote":"Documents conic textures in ferroelectric nematics attributed to splay stiffening, cited as independent evidence for the effect."},{"cited_title":"Z ˜ρ(k)eik·r dk (2π)3 #","cited_arxiv_id":null,"evidence_quote":"Predicts helical twist configurations from dipole-dipole interactions, cited to explain the twist softening in the $N_F$ phase."},{"cited_title":"Caimi, G","cited_arxiv_id":null,"evidence_quote":"Reports fluid superscreening and polarization following in confined ferroelectric nematics, supporting the screened-electrostatics picture."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides out-of-plane reorientation measurements in ferroelectric nematics used as a comparison for the magneto-optical threshold behavior."}],"review_version":1}