{"id":"90993472-1acb-4e57-95e9-4ba6900a1152","arxiv_id":"2608.11067","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Spontaneous striped order in nematic dimers is explained by flexoelectric-elastic competition, and high-frequency field-induced bands match the inertial conduction instability.","lead":"In two flexible-core liquid crystal dimers, the authors mapped pattern formation that occurs spontaneously and under electric fields, and proposed a simple energy-balance explanation for the spontaneous stripes. The work tests a predicted high-frequency instability and is relevant to researchers studying twist-bend nematics and bent-core soft materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sec. III.B's zero-crossing is arithmetically inconsistent: stated constants give dq≈4, not dq≈10.","rationale":"I re-derived Eqs. (5)-(7) from the director ansatz (1). The elastic energy integrates to Eq. (6); the flexoelectric half-wave dipole p=4de3φ0/π and the dipole-dipole interaction give Eq. (5). Setting E_fluct=0 with the stated constants and ξ=d yields a quartic in x=dq: 4.70×10^-14 x^4 = 3.93×10^-13 x^2 + 5.81×10^-12. The positive root is x≈4.0. At x=10, the left side is 4.70×10^-10 vs right 4.5×10^-11, so flexo dominates by ~10×. Thus the paper's statement 'It crosses 0 when dq≈10' is not reproducible. This matters because the entire explanatory claim for the striped ground state—and the specific predicted wavelength λ≈0.63d—hangs on this number. If the zero-crossing is actually at dq≈4, the model predicts λ≈1.6d, which conflicts with the observed 'a few µm' stripes in cells of d=5–10 µm (observed λ≈0.63d). The reading that the model is only approximate does not rescue the specific quantitative agreement claim. The screening assumption is also a valid concern, but this arithmetic check is more fundamental and directly testable. The authors should either correct the numerical evaluation, use measured constants for 5O9/6O9 (not borrowed from other dimers), or remove the quantitative λ claim. Because the issue is fixable, I recommend CONDITIONAL rather than outright rejection, but the condition is non-negotiable.","tokens_in":17956,"tokens_out":16927,"duration_ms":136758,"concrete_test":"Compute E_flex from Eq. (5) and E_elast from Eq. (6) at dq=10 using k33=1 pN, k22=1.5 pN, e3=10 pC/m, ε0=8.854×10^-12 F/m, ξ=d. If |E_flex| is not equal to E_elast (it is ≈10× larger), the claimed zero-crossing at dq≈10 is refuted; solve the balance equation to find the actual dq (≈4.0).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claim of the SGS model is that E_flex+E_elast crosses zero at dq≈10, giving λ≈0.63d. Substituting the paper's own numbers (k33=1 pN, k22=1.5 pN, e3=10 pC/m, ξ=d) into Eqs. (5) and (6) and solving gives dq≈4.0, not 10. At dq=10, |E_flex|/E_elast≈10.4, so the cancellation has long since occurred; the zero-crossing happens near dq≈4 (λ≈1.6d), which does not match the observed 0.63d. To force a crossing at dq=10 one would need e3≈3 pC/m (with the stated k's) or k's≈10 pN (with e3=10 pC/m), neither of which is reported. This is a purely algebraic inconsistency, independent of the screening assumption flagged by the reader. Since the paper's headline agreement between model and observation rests on this number, the central claim is not supported as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports an experimental study of two dielectrically negative twist-bend nematic dimers, P5OBD9 and P6OBD9. It describes a spontaneously formed, quasiperiodic striped ground state in the uniaxial nematic phase and attributes this state to a balance between a negative flexoelectric interaction and the positive elastic energy of azimuthal director fluctuations (Sec. III.B). The paper further documents a broad phenomenology of field-induced patterns: surface-localized electroconvection and Bobylev–Pikin flexoelectric domains in DC fields, transient surface and volume instabilities at sub-hertz frequencies, an oblique-roll → bimodal-grid → normal-roll → chevron sequence at low frequencies, and wide normal rolls above 100 kHz, which are interpreted as the inertial conduction instability with an extracted effective elastic constant K≈1 pN.","tokens_in":18175,"tokens_out":12972,"duration_ms":113245,"significance":"The paper's main value lies in its rich experimental phenomenology and the schematic U–f phase diagram for patterned states in flexible-core twist-bend nematogens. The use of 90° twist cells to distinguish surface and bulk effects is a clear strength, as is the identification of polarity-sensitive surface electroconvection and the nucleation-and-growth character of flexoelectric domains. The comparison with the Pikin–Chigrinov inertial mode is a plausible consistency check, although it is not a parameter-free test. However, the central quantitative claim of the striped-ground-state model is undermined by an arithmetic inconsistency that must be resolved before the model can be accepted.","major_comments":[{"comment":"The claimed zero-crossing at dq≈10 is not supported by the paper's own parameters. Substituting ξ=d, k33=1 pN, k22=1.5 pN and e3=10 pC/m into Eqs. (5) and (6) gives E_flex=−(4 d e3² (dq)³)/(π⁶ ε₀) φ₀² and E_elast=(π d/8)[k33(dq)+k22π²/(dq)]φ₀². Equating these two terms yields dq≈4.0, not 10; the corresponding wavelength is λ≈1.6d, not 0.63d. At dq=10 the flexoelectric term is already about 10 times larger than the elastic term, so no cancellation occurs there. To force a zero-crossing at dq=10 with the stated elastic constants one would need e3≈3 pC/m, while retaining e3=10 pC/m would require elastic constants of order 10 pN; neither is reported. This removes the model's quantitative agreement with the observed stripe wavelength.","section":"III.B, Eqs. (5)–(7)"},{"comment":"The model discards surface-charge contributions solely because neighboring half-waves are oppositely charged and assumed to be screened by free ions. This assumption is not tested or quantified. If screening is incomplete, the depolarizing field opposing the flexoelectric polarization will reduce or even eliminate the negative interaction term in Eq. (5), changing the zero-crossing condition and hence the predicted wavelength. The authors should either estimate the screening length relative to d and the stripe period, or show that the model's conclusions are robust to partial screening.","section":"III.B, after Eq. (4)"},{"comment":"The extraction of K≈1 pN assumes that σ⊥ represents both σ and σ||−σ⊥, because homeotropic cells were unavailable. This equality is not measured, so the near-linearity of U_NR with f√(ε⊥/σ⊥) is at least partly built into the assumed form of Eq. (9). The identification of the high-frequency wide domains with the inertial conduction instability is plausible and interesting, but it should be presented as a consistency check under an explicit auxiliary assumption, not as a determination of K.","section":"III.F, Eq. (9) and Fig. 20"}],"minor_comments":[{"comment":"The flexoelectric coefficient e3 has units C/m; the sentence 'e3 is about 10 pC/m 2' should read '10 pC/m'.","section":"III.B, text near Eq. (5)"},{"comment":"The phrase 'twit-bend nematogen' contains a typo and should be 'twist-bend nematogen'.","section":"III.A, penultimate paragraph"},{"comment":"The caption describes 'exponential frequency variations' of U_NR and 1/√σ; please clarify whether an exponential fit is actually intended or whether 'monotonic decrease' is the correct description.","section":"Fig. 18 caption"},{"comment":"Threshold-voltage measurements in Figs. 15 and 18 are presented without error bars; an estimate of the uncertainty in U_NR and U_G would strengthen the quantitative claims.","section":"III.F, Fig. 20"},{"comment":"The quoted statement about σ(σ||−σ⊥) depending strongly on frequency would benefit from a page or equation reference to the source, since it is used to justify the decreasing U_NR(f).","section":"III.F, quotation from Ref. [40]"}],"recommendation":"major_revision","confidential_remarks":"The arithmetic error in Sec. III.B is the decisive issue: the model as written contradicts its own stated parameters and does not reproduce the observed stripe wavelength. If the authors can correct the model or show that a different but justified choice of parameters restores the dq≈10 zero-crossing, the substantial experimental material could support publication. If the numbers in Eqs. (5)–(7) are maintained as they stand, the central claim of the paper fails. I recommend major revision rather than rejection because the error is localized and fixable in principle, and because the experimental phenomenology alone is of considerable value."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThe paper is worth a look for the experimental survey alone. Krishnamurthy et al. map out an impressively complete set of patterned states in two dielectrically negative flexible-core dimers: the spontaneous quasiperiodic ground state, surface and bulk electroconvection, Bobylev-Pikin flexo-domains that grow by nucleation and front propagation, the sequence of oblique-roll/grid/normal-roll/chevron, and the high-frequency wide domains that fit the inertial conduction instability reasonably well. The observation that the flexo-domains grow as 2D tactoids is novel, and the demonstration that the very high-frequency threshold decreases with frequency while scaling with f√(ε⊥/σ⊥) is a nice confirmation of the Pikin-Chigrinov inertia mode.\n\nThe problem is the model for the striped ground state (SGS). The paper claims that with the stated constants (k33=1 pN, k22=1.5 pN, e3=10 pC/m, ξ=d), the fluctuation energy E_flex+E_elast crosses zero at dq≈10, giving λ≈0.63d. When I plug those numbers into their own Eqs. (5) and (6), the crossing happens at dq≈4.0, giving λ≈1.6d. At dq=10 the flexo term is about 10 times the elastic term, so the cancellation has long passed. To get a crossing at dq=10 you would need e3≈3 pC/m or elastic constants near 10 pN, neither of which is reported. This is a load-bearing inconsistency: the paper's headline agreement with the observed stripe spacing rests on it. The qualitative idea—flexoelectric-elastic competition driving azimuthal fluctuations—is plausible, but the concrete prediction is not supported as written.\n\nThe high-frequency comparison also has a known soft spot: σ⊥ is used as a proxy for both σ and (σ||−σ⊥), which is a crude simplification, and the extracted K≈1 pN has no error bars. That part is better treated as a consistency check than a quantitative test.\n\nThe experiments themselves seem careful and the systematic observations are valuable. The audience is the twist-bend nematic and bent-core electroconvection community; they will find the phase diagram and many new textures useful. The SGS model needs a corrected calculation (and perhaps a statement about the dielectric constant of the medium, which is ignored). This is a case where a good experimental paper has a weak theoretical centerpiece that can be fixed.\n\nWould I send it to review? Yes, because the experimental content deserves a serious referee and the model article is correctable. But the referee should be asked to verify the arithmetic.\n\nBest.","headline":"Rich experimental phenomenology, but the central SGS model's quantitative claim is contradicted by its own arithmetic.","tokens_in":18721,"tokens_out":13451,"would_cite":false,"duration_ms":105260,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper traces the spontaneous striped ground state in the nematic phase of flexible-core phenyl benzoate dimers to a competition between flexoelectric and elastic energies, predicting a stripe wavelength near $0.63$ times the cell…","keywords":["nematic liquid crystals","twist-bend nematogens","flexoelectricity","striped ground state","electroconvection","inertial conduction instability","phenyl benzoate dimers","director fluctuations"],"falsifier":"Measure the field-free stripe spacing in planar cells of several thicknesses $d$ at a fixed reduced temperature: the model's energy balance crosses zero at $dq\\approx 10$, so the spacing should scale as $\\lambda\\approx 0.63d$ and the pattern should lose its long lifetime in much thinner cells. A second test targets the screening assumption: reduce the free-ion content and the stripes should weaken or vanish if the model's neglect of surface-charge energy is correct, and persist if it is not. For the high-frequency branch, plot $U_{\\rm NR}$ against $f\\sqrt{\\varepsilon_\\perp/\\sigma_\\perp}$ in a material whose dielectric relaxation occurs at a different frequency: the same near-linear relation should hold, with a slope still consistent with $K\\approx 1$ pN.","tokens_in":17796,"feed_emoji":"🔬","tokens_out":17176,"duration_ms":141769,"temperature":0.7,"pith_summary":"This paper tries to establish why the nematic phase of two flexible-core twist-bend dimers (phenyl benzoates with a nonamethylene spacer) is spontaneously striped, even with no field applied. The stripes are azimuthal director undulations, and the mechanism the paper proposes is a competition between flexoelectric energy—electric polarization created by bend distortion—and ordinary elastic energy. With measured elastic and flexoelectric constants, the net free energy of a stripe crosses zero at a wavevector set by $dq\\approx 10$, meaning the stripe wavelength should be near $0.63$ times the cell thickness, close to the observed spacing. The paper also claims that the very high frequency periodic wide domains in the same materials are the inertial conduction instability, with a threshold that falls as frequency rises. A sympathetic reading is that the work unifies the field-free pattern and a long list of field-induced instabilities under one phase diagram rooted in the bend-flexibility and flexoelectricity of bent molecules.","feed_headline":"Nematic stripe spacing set by flexoelectric–elastic balance","feed_subtitle":"The model predicts stripes spaced near 0.63 times cell thickness and groups every field-driven pattern into one phase diagram.","key_machinery":"The load-bearing object is the one-mode azimuthal fluctuation ansatz $\\varphi(x,z)=\\varphi_0\\cos(\\pi z/d)\\sin(qx)$: it puts the maximum director deviation at the midplane, satisfies strong anchoring at the plates, and picks out a wave vector $q$ along the rubbing direction. The argument is carried by the free-energy balance $E_{\\rm fluct}=E_{\\rm flex}+E_{\\rm elast}$. The bend component of the fluctuation produces a flexoelectric polarization (electric polarization induced by a bend distortion) $\\mathbf{P}_{\\rm flex}=\\mathbf{j}\\, e_3 q \\varphi_0\\cos(\\pi z/d)\\cos(qx)$; integrating over half a period and the cell thickness treats each half-wave as a line dipole, giving a negative neighbor interaction $E_{\\rm flex}\\approx -(1/4\\pi\\varepsilon_0)(4de_3\\xi/\\pi)^2/(\\pi/q)^3\\,\\varphi_0^2$. The opposing cost is the bend-plus-twist elastic energy $E_{\\rm elast}=(\\pi\\xi/8)(k_{33}dq + k_{22}\\pi^2/(dq))\\varphi_0^2$. Setting $E_{\\rm fluct}\\approx 0$ with the measured constants makes $dq\\approx 10$, or $\\lambda\\approx 0.63d$. For the high-frequency instability, the named mechanism is the inertial conduction instability, whose threshold is $U_{\\rm NR}=2\\pi f\\sqrt{K\\varepsilon\\varepsilon_0/[\\sigma(\\sigma_\\parallel-\\sigma_\\perp)]}$; the data's near-linear plot of $U_{\\rm NR}$ against $f\\sqrt{\\varepsilon_\\perp/\\sigma_\\perp}$ and the implied $K\\approx 1$ pN are what attach the observations to that mechanism.","core_discovery":"The paper's central claim is that the quasiperiodic striped ground state in these dielectrically negative, bend-flexible nematics is a bulk phenomenon: a long-lived azimuthal modulation of the director with maximum amplitude at the cell midplane. The modulation is described by $\\varphi(x,z)=\\varphi_0\\cos(\\pi z/d)\\sin(qx)$, and its survival is attributed to a near cancellation between the positive elastic cost of bend and twist and a negative interaction energy between neighboring half-waves carrying opposite flexoelectric polarization. Inserting measured values—$k_{33}\\approx 1$ pN, $k_{22}\\approx 1.5$ pN, $e_3\\approx 10$ pC/m—makes the net fluctuation energy pass through zero at $dq\\approx 10$, i.e. $\\lambda\\approx 0.63d$, which the authors say matches the observed stripe spacing. On the electrical side, the paper's second claim is that the wide domains appearing above roughly 100 kHz are the inertial conduction instability: the threshold voltage is nearly linear in $f\\sqrt{\\varepsilon_\\perp/\\sigma_\\perp}$, and the fitted slope gives an effective elastic constant $K\\approx 1$ pN. The paper reads the full sequence of static, low-frequency, and high-frequency patterned states as a single voltage–frequency phase diagram sitting on top of the striped ground state.","pith_inferences":["This model suggests a screening experiment that is not in the paper: systematically doping or purifying the ionic content should modulate the stripe contrast and lifetime, because the negative flexoelectric term depends on free ions neutralizing surface charges.","If the mechanism is generic, then any nematic with sufficiently small $k_{33}$ and large $e_3$—not only twist-bend dimers—should show quasiperiodic azimuthal fluctuations; one could screen candidate compounds by elastic and flexoelectric measurements rather than waiting for stripes to appear.","The same balance, taken to shorter wavelengths, could set the scale of the pretransitional fluctuations as $T\\to T_{\\rm TB}$; the paper does not pursue this, but its zero-energy condition is a natural starting point for such a calculation.","For the high-frequency branch, the prediction could be used to identify the inertial conduction instability in calamitic nematics with suitable dielectric relaxation, which would show that a falling threshold with frequency is not exclusive to bent-core materials."],"forward_implications":["If the ground-state model holds, the stripe wavevector in planar cells is controlled by cell thickness through $dq\\approx 10$, so changing $d$ should move the observed spacing and the pattern should disappear in sufficiently thin cells.","The striped ground state underlies all electrically induced states, giving a single voltage–frequency phase diagram in which static and very low frequency fields select surface electroconvection or volume flexoelectric domains, intermediate frequencies select roll–grid–chevron sequences, and very high frequencies select wide bands.","The high-frequency threshold should continue to fall with frequency wherever $\\varepsilon_\\perp$ relaxes and $\\sigma_\\perp$ rises, and the effective elastic constant extracted from the slope should stay near 1 pN.","The low bend elastic constant and enhanced flexoelectric coefficient, both linked to the molecular bend, become the measurable predictors of whether a nematic will show spontaneous stripes."],"supporting_citations":[{"why":"Supplies the measured bend and twist elastic constants ($k_{33}\\approx 1$ pN, $k_{22}\\approx 1.5$ pN) used as inputs to the energy balance.","marker":"[27]"},{"why":"Shows large effective flexoelectric coefficients in twist-bend dimers, supporting the large $e_3$ premise.","marker":"[28]"},{"why":"Adds flexoelectrooptic evidence that the bend flexo coefficient is enhanced in these dimeric nematics.","marker":"[29]"},{"why":"Reports a very large bend flexo coefficient $e_3$ in a related twist-bend nematic, behind the value-scale used in the model.","marker":"[30]"},{"why":"Documents the nonstandard prewavy wide-domain states in rigid bent-core nematics that the high-frequency comparison extends.","marker":"[34]"},{"why":"Provides the inertia-mode electroconvection treatment whose threshold relation is tested against the measured $U_{\\rm NR}$.","marker":"[35]"},{"why":"Describes the earlier surface-based interpretation of the bent-core striped state that the paper's bulk model replaces.","marker":"[37]"},{"why":"Gives the standard thermal-fluctuation treatment of the director used to set up the stripe fluctuation ansatz.","marker":"[38]"},{"why":"Original inertial conduction instability prediction, used as the theoretical identification of the very high frequency domains.","marker":"[46]"}],"fun_headline_variants":["Nematic stripe spacing from flexoelectric-elastic balance","Flexoelectric-elastic balance sets nematic stripe pattern","Bulk flexoelectric effect creates nematic stripes","Stripe width in bend-flexible nematics: a balance","Flexoelectric-elastic tug shapes nematic instabilities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire energy balance assumes that free ions in the liquid screen the surface charges of neighboring striped half-waves, so those charges cost no energy; if that screening is incomplete, the negative flexoelectric term shrinks and the predicted $0.63d$ spacing no longer follows.","fun_headline_variants_meta":{"raw":{"variants":["Nematic stripe spacing from flexoelectric-elastic balance","Flexoelectric-elastic balance sets nematic stripe pattern","Bulk flexoelectric effect creates nematic stripes","Stripe width in bend-flexible nematics: a balance","Flexoelectric-elastic tug shapes nematic instabilities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000226,"raw_usage":{"total_tokens":1548,"prompt_tokens":1108,"completion_tokens":440,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":724,"completion_tokens_details":{"reasoning_tokens":363}},"tokens_in":724,"tokens_out":440,"duration_ms":4839,"temperature":1.0,"reasoning_tokens":363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:54:33.976326+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the field-free stripe spacing in planar cells of several thicknesses $d$ at a fixed reduced temperature: the model's energy balance crosses zero at $dq\\approx 10$, so the spacing should scale as $\\lambda\\approx 0.63d$ and the pattern should lose its long lifetime in much thinner cells. A second test targets the screening assumption: reduce the free-ion content and the stripes should weaken or vanish if the model's neglect of surface-charge energy is correct, and persist if it is not. For the high-frequency branch, plot $U_{\\rm NR}$ against $f\\sqrt{\\varepsilon_\\perp/\\sigma_\\perp}$ in a material whose dielectric relaxation occurs at a different frequency: the same near-linear relation should hold, with a slope still consistent with $K\\approx 1$ pN.","supporting_citations":[{"cited_title":"Babakhanova, Z","cited_arxiv_id":null,"evidence_quote":"Supplies the measured bend and twist elastic constants ($k_{33}\\approx 1$ pN, $k_{22}\\approx 1.5$ pN) used as inputs to the energy balance."},{"cited_title":"Balachandran, V","cited_arxiv_id":null,"evidence_quote":"Shows large effective flexoelectric coefficients in twist-bend dimers, supporting the large $e_3$ premise."},{"cited_title":"Varanytsia and L","cited_arxiv_id":null,"evidence_quote":"Adds flexoelectrooptic evidence that the bend flexo coefficient is enhanced in these dimeric nematics."},{"cited_title":"Škarabot, N","cited_arxiv_id":null,"evidence_quote":"Reports a very large bend flexo coefficient $e_3$ in a related twist-bend nematic, behind the value-scale used in the model."},{"cited_title":"Wiant, J","cited_arxiv_id":null,"evidence_quote":"Documents the nonstandard prewavy wide-domain states in rigid bent-core nematics that the high-frequency comparison extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the inertia-mode electroconvection treatment whose threshold relation is tested against the measured $U_{\\rm NR}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the earlier surface-based interpretation of the bent-core striped state that the paper's bulk model replaces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the standard thermal-fluctuation treatment of the director used to set up the stripe fluctuation ansatz."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original inertial conduction instability prediction, used as the theoretical identification of the very high frequency domains."}],"review_version":1}