{"id":"54c7ba35-be58-4368-bd1b-623bda485e9a","arxiv_id":"2411.14115","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A systematic study of 27 fluorinated liquid crystals shows that specific fluorination patterns, rather than maximal fluorination or large dipole moments, best stabilize the ferroelectric nematic phase.","lead":"This paper synthesizes 27 related liquid crystal molecules and maps how the placement of fluorine atoms determines whether they form a ferroelectric nematic phase, a fluid that is spontaneously polarized. The authors find that specific, not maximum, fluorination patterns give the most stable polar phases, and link this to the pattern of electric charge on the molecular surface.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central design rule rests on an undefined 'oscillation' metric; no quantitative link between 1D ESP shape and polar-phase stability is established.","rationale":"The paper delivers a substantial experimental dataset: 27 homologues, full phase characterization by DSC, POM, X-ray scattering, and polarization measurements, with openly deposited data. That part is credible and valuable. The reader's CONDITIONAL verdict is appropriate. My concern differs slightly from the reader's weakest_assumption: I agree that the 0.0004-isovalue static-DFT radial averaging is a modeling premise, but the sharper problem is that even within that protocol the central descriptor 'oscillatory' is never operationalized. The 1D ESP curves are inspected and described as 'almost sinusoidal' or 'lacking a clear oscillatory structure', but no feature of those curves is measured, thresholded, or regressed against phase behavior. Consequently the central claim — that surface-charge oscillation, not dipole magnitude, is the key design variable — is not yet testable, even in principle. A quantitative index and an out-of-sample check would settle it. This strengthens the condition the reader already attached to the verdict, but does not change the verdict itself. The SI characterization errors (e.g., duplicated 13C data for compound 8, mislabeled compound numbers in Figure 5a) are not load-bearing for the central claim, though they reinforce the need for careful revision before publication.","tokens_in":32553,"tokens_out":4139,"duration_ms":40735,"concrete_test":"Define a quantitative oscillation index, e.g. the summed squared amplitudes of the first three Fourier harmonics of the 1D ESP profile along the long axis normalized by the mean, or a spatial-uniformity measure of positive/negative ESP patches. Compute it for all 27 homologues and test rank correlation with T_NF-I and T_NX-N (treating absence of a polar phase as censored). Also recompute the ESP profiles at isovalues 0.0002 and 0.0008 and with an alternative functional such as ωB97X-D; if the index ranking among compounds 1-9 changes materially, the proposed rule is an artifact of the chosen protocol.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that molecules with 'a more oscillatory distribution of electrons across their surfaces' have a higher propensity to form polar nematic phases, and that this explains the observed X.Y.1 trend. As presented, this is a qualitative visual classification of the 1D radially averaged ESP curves (Figure 4, SI Figure S3), not a measurable quantity. No algorithm, threshold, or scalar index distinguishes 'oscillatory/uniform' from 'non-oscillatory/pronoounced' curves; the comparison between 2.2.Z (polar) and 2.0.Z (non-polar) is made by eye after the phase behavior is known. Because the central design principle is stated as a general rule, it must be possible to compute a number from an ESP surface that predicts propensity to polar order. The current definition also depends on choices in SI 1.6: the 0.0004 isodensity, static B3LYP-GD3BJ geometry, and the radial-averaging/rescaling procedure. There is no test that the qualitative ranking is robust to these choices. Until the oscillation rule is operationalized and shown to correlate with polar-phase stability out-of-sample, the headline claim is not falsifiable. This is load-bearing because if another isovalue, conformer, or density functional changes which compounds look 'oscillatory', the proposed design rule evaporates.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the synthesis and characterization of 27 systematically fluorinated biphenyl benzoate liquid crystals, identifying ferroelectric nematic (NF) and antiferroelectric nematic (NX) phases in several homologues. The key experimental observation is that compounds with the X.Y.1 fluorination pattern exhibit higher polar-phase transition temperatures than their more fluorinated analogues despite having smaller total dipole moments. Using DFT-computed electrostatic potential (ESP) surfaces, the authors argue that molecules with a more oscillatory radial ESP distribution along the long axis have a higher propensity to form polar nematic phases, and they support this with bimolecular potential energy surface calculations. The abstract presents the ESP-oscillation criterion as a general design rule for future ferroelectric nematogens.","tokens_in":32765,"tokens_out":10734,"duration_ms":93772,"significance":"If the central claim is established, the paper would provide a valuable molecular design principle: radial surface-charge pattern, rather than dipole-moment magnitude, determines polar nematic stability. The experimental dataset itself is a strength: a systematic 27-compound series with phase assignments supported by POM, DSC, current-response measurements, and X-ray scattering for representative compounds, and the X.Y.1 trend is an empirical observation that stands independently of the computational interpretation. The DFT calculations are standard and no parameters are fitted to the phase data. However, the headline ESP-oscillation rule is currently qualitative and retrospective, which limits its immediate utility as a predictive design rule. The paper has the potential to be an important contribution to the chemistry of ferroelectric nematics once the computational claim is made testable.","major_comments":[{"comment":"The central design principle is expressed entirely in qualitative terms: curves are described as \"oscillat[ing] almost sinusoidally\", \"more pronounced, lacking a clear oscillatory structure\", and \"spatially uniform\", with no algorithm, threshold, or scalar index defining \"oscillatory\" versus \"non-oscillatory\". Because the abstract states this as a general rule (\"molecules possessing a more oscillatory distribution of electrons across their surfaces possessing a higher propensity to form polar nematic phases\"), the rule should be operationalized so that an independent reader can compute the descriptor from the 1D ESP data and test the correlation with polar-phase propensity. Without such a metric, the central claim is not falsifiable.","section":"Results and Discussion, Figure 4"},{"comment":"The ESP descriptor depends on the chosen electron-density isovalue (0.0004), the static B3LYP-GD3BJ/cc-pVTZ geometry, and the radial-averaging/rescaling procedure. The manuscript reports no test that the qualitative distinction between \"oscillatory\" and \"uniform\" curves is robust to these choices, for example by varying the isovalue over a reasonable range or re-optimizing with a different functional. Since the design rule is presented as a general molecular-level principle, the absence of a robustness check leaves open the possibility that the classification is an artifact of the calculation protocol.","section":"SI section 1.6"},{"comment":"The classification of compounds into \"oscillatory\" versus \"non-oscillatory\" is performed on the same 27 molecules whose phase behavior is already known, and the text explicitly selects examples based on that knowledge (e.g., 2.2.Z vs 2.0.Z). This is a retrospective rationalization rather than a predictive test. The authors should either state this limitation explicitly or provide an out-of-sample validation (e.g., predict the phase behavior of a newly synthesized compound) to support the claim that the criterion is a \"design strategy\".","section":"Results and Discussion, \"Inspection of these plots...\""},{"comment":"The statement that \"for none of 1-27 does the global minima in complexation energy for the parallel packed molecules become lower than antiparallel\" is presented as a general result, but the manuscript shows data only for 2.2.2 in Figure 5. If this claim is used to support the entropy-based argument for parallel packing, the complexation energies, or at least the differences between parallel and antiparallel global minima, should be tabulated for all compounds in the SI.","section":"Results and Discussion, Figure 5"},{"comment":"The claim that X.Y.1 homologues have higher polar-phase transition temperatures \"whilst possessing smaller values of µ\" is based on total dipole magnitude. However, the longitudinal component, which can be computed as μ·cos(θ) from Table S3, is actually slightly larger for the Z=1 compounds than for their Z=2 counterparts in each X.Y series (e.g., 2.2.2 vs 2.2.1: about 11.52 D vs 11.57 D; 1.2.2 vs 1.2.1: about 10.64 D vs 10.70 D). Since the paper frames the design in terms of the longitudinal dipole moment, the authors should report longitudinal components and show that the empirical trend is not explained by them, or else avoid making the \"smaller dipole moment\" claim based on total magnitude.","section":"Figure 3f and Table S3"}],"minor_comments":[{"comment":"Several figure references are incorrect: \"Figure 1b(ii)\", \"Figure 1g\", and \"Figure 1h\" should refer to panels in Figure 3, and the ESP panel referenced as \"Figure 2a\" after the radial-averaging description should be Figure 4a.","section":"Throughout"},{"comment":"The caption labels compound 11 as \"2.1.1\", but Table S1 and the main text identify compound 11 as having the 0.2.1 fluorination pattern.","section":"SI, Figure S9"},{"comment":"The caption states \"Molecular structure of 2∙2∙2 (27)\", but compound 27 is 0.0.0 in Table S1; the structure shown is compound 1 (2.2.2).","section":"Figure 5 caption"},{"comment":"The text says \"the conformation of the transition temperature by differential scanning calorimetry\"; \"conformation\" should read \"confirmation\".","section":"Results and Discussion"},{"comment":"Reference [45] (Hess et al., LINCS) is a molecular dynamics constraint algorithm, not a DFT or dispersion-correction method; it appears to be an incorrect citation for the B3LYP-GD3BJ method, which should instead cite Grimme et al. for GD3BJ and Dunning for cc-pVTZ.","section":"References"},{"comment":"Several width entries are malformed, e.g., \"0.5.09\" and \"0.4.99\", and should be decimal values (0.509, 0.499).","section":"Table S3"},{"comment":"The phrase \"due to it's the potential\" is ungrammatical and should read \"due to its potential\".","section":"Introduction"},{"comment":"The word \"complimented\" in \"We have complimented our synthetic efforts\" should be \"complemented\".","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The experimental work is systematic and the phase assignments appear sound; the concern is the unsupported generality of the ESP-oscillation rule as currently presented. A revision that operationalizes the metric, provides a robustness check, and either validates the rule out of sample or explicitly frames it as a retrospective hypothesis would make the paper suitable for publication. The citation and labeling issues are minor and easily fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you work on ferroelectric nematics. The experimental core is solid: 27 homologues, systematic fluorination at every position on a biphenyl benzoate core, full phase characterization by POM, DSC, current response, and X-ray for representative compounds, plus a complete transition temperature table. The dataset alone is a contribution. The observation that Z=1 fluorination patterns give more stable polar phases than maximal fluorination, despite smaller dipole moments, is new and convincingly documented. That's the part I'd trust.\n\nThe soft spot is exactly what the stress-test flagged. The 'oscillatory ESP' rule is not actually defined. The paper shows 1D radially averaged ESP curves and says some are more oscillatory/uniform than others, but there is no metric, no threshold, no out-of-sample test. The comparison is made by eye after the phase behavior is known. The ESP isovalue (0.0004), the static DFT geometry, and the radial-averaging procedure are reasonable choices, but there is no evidence the qualitative ranking survives changing them. As stated, the headline claim is not falsifiable. That is a real weakness, and it is load-bearing only for the explanatory story. The Z=1 trend itself stands independently of the ESP interpretation.\n\nA few smaller issues: compound 8's 13C data in the SI is duplicated from compound 5, figure S9 mislabels compound 11 as 2.1.1 instead of 0.2.1, and Table S3 has some obvious typos. These are presentation errors, not fatal. The bimolecular PES refinement uses five minima per packing mode, which is a free choice, but it is not fitted to the phase data, so it is not circular.\n\nWho is this for? Synthetic chemists and soft matter experimentalists who want concrete fluorination patterns to try in new NF candidates. The paper deserves a serious referee. The recommendation would be: publish the dataset and the experimental trends, but send it back for a quantitative operationalization of the ESP-oscillation criterion and a robustness check across isovalues or functionals, or at the least an explicit statement that it is a qualitative heuristic. I would not block publication over the ESP rule, but I would not let the 'ground rules' phrasing stand as is.","headline":"A genuinely useful 27-compound fluorination screen with a plausible Z=1 design trend, but the central 'oscillatory ESP' rule is still a qualitative after-the-fact descriptor rather than a testable criterion.","tokens_in":33313,"tokens_out":1491,"would_cite":true,"duration_ms":16559,"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":"Specific fluorination patterns, not maximal fluorination or dipole size, decide which liquid crystals form ferroelectric nematic phases; the controlling feature is how evenly the surface charge oscillates along the molecule.","keywords":["ferroelectric nematic","liquid crystal","fluorination","electrostatic potential","molecular design","polar nematic phase","antiferroelectric nematic","density functional theory"],"falsifier":"Prepare a 2.2.1 analogue in which the outer-ring fluorine is moved from the position adjacent to the ester carbonyl to the opposite side of the benzoate ring; the spatial-uniformity model predicts a drop in T_NF-I, while a pure dipole-moment model predicts little change, so the comparison would separate the two mechanisms. Alternatively, compute the 1D ESP for known ferroelectric nematogens outside this scaffold and check whether the oscillation-uniformity ranking tracks which compounds actually form polar phases.","tokens_in":32326,"feed_emoji":"⚡","tokens_out":8460,"duration_ms":79417,"temperature":0.7,"pith_summary":"Ferroelectric nematics are liquids that flow yet carry a macroscopic electric polarisation comparable to solid ferroelectrics, and this paper asks which molecular features actually favour that polar order. To answer, the authors synthesised all twenty-seven fluorination patterns of a biphenyl benzoate core and measured which ones form polar nematic or antiferroelectric nematic phases. The central finding is that maximal fluorination is not optimal: homologues bearing one fluorine on the outer benzoate ring (X.Y.1 patterns) show the highest polar-phase transition temperatures even though their dipole moments are smaller than those of more fluorinated analogues. Electronic-structure calculations trace this to the molecule's surface charge: when the radially averaged electrostatic potential along the long axis oscillates evenly between positive and negative regions, parallel neighbours pack more favourably and polar order is stabilised. A sympathetic reader would take this as a concrete design rule for future ferroelectric nematics: engineer the surface-charge pattern, not just the dipole moment.","feed_headline":"Specific fluorination, not the biggest dipole, wins in polar nematics","feed_subtitle":"Across 27 homologues, a uniform oscillatory surface-charge pattern—not the largest dipole—stabilises polar order.","key_machinery":"The central object is the molecule's longitudinal charge-density wave: the radially averaged electrostatic potential (ESP) plotted along the long molecular axis at an electron density isovalue of 0.0004, computed from DFT (B3LYP-GD3BJ/cc-pVTZ) geometries. This 1D reduction of the 3D ESP is the diagnostic that separates the polar-phase-forming homologues (uniform, near-sinusoidal oscillations) from the apolar ones (large, non-oscillatory charge variation). The argument is carried further by rigid bimolecular potential-energy scans, which map the complexation energy of two molecules translated over x/y/z in parallel and antiparallel orientations, and by Interaction Region Indicator (IRI) isosurfaces that visualise the non-covalent contacts; together they show offset π-π stacking dominates and that parallel packing offers multiple comparable-energy minima while antiparallel packing has repulsive regions.","core_discovery":"Against the expectation that more fluorination means a bigger dipole and therefore a more stable ferroelectric nematic phase, the paper reports that in a family of twenty-seven biphenyl benzoates the polar phases (NF and NX) appear only in the most fluorinated subset (1-9), and within that subset the most stable polar phases are not the most fluorinated molecules. Removing one fluorine from the outer benzoate ring to make the 2.2.1 homologue raises the NF-isotropic transition from a monotropic 133.5°C for 2.2.2 to an enantiotropic 145.4°C, while lowering the calculated dipole moment; the same X.Y.1 advantage repeats across the 2.2, 2.1 and 1.2 series. The authors' explanation is that the relevant molecular property is the 1D radially averaged electrostatic potential at the 0.0004 electron-density isosurface, which should oscillate uniformly along the long axis so that parallel side-by-side molecules see complementary positive and negative regions. DFT-based bimolecular potential-energy scans and Interaction Region Indicator analysis support the picture: antiparallel packing carries repulsive regions that the parallel arrangement avoids, the dominant interaction is offset π-π stacking rather than point-dipole attraction, and increased fluorination brings the parallel and antiparallel energy minima closer together. The paper's central claim is therefore that the propensity to form polar nematic phases is governed by the spatial uniformity of the surface-charge oscillation, not by the magnitude of the molecular dipole moment.","pith_inferences":["If the surface-charge-oscillation rule transfers to other molecular cores, then screening candidate ferroelectric nematogens by the uniformity of their 1D ESP wave, without synthesising them first, could accelerate discovery of room-temperature NF materials; this is an inference, not a paper claim.","The paper assumes free rotation about the long axis when radially averaging the ESP, so the rule is most confidently applied to rod-like mesogens with low rotational barriers; for bent or bulky cores, 3D ESP anisotropy may need to be included.","The entropic argument that several comparable parallel packing minima offset an enthalpic penalty suggests a testable design strategy: adding substituents that create multiple nearly-degenerate parallel contacts could further stabilise polar order, a direction the paper gestures toward but does not demonstrate.","Quantifying the ESP wave (amplitude, wavelength, phase) and correlating it with T_NF-I across a broader dataset could turn the qualitative oscillation criterion into a quantitative descriptor for machine-learning-guided materials design."],"forward_implications":["Synthetic strategy: for this biphenyl benzoate scaffold, the X.Y.1 fluorination pattern is the target, not maximal fluorination, because it gives the highest polar-phase transition temperatures.","Prediction metric: the 1D radially averaged ESP plot can serve as a computational pre-screen for whether a candidate rod-like molecule is likely to form a polar nematic phase.","Phase control: single-fluorine edits switch materials between ferroelectric nematic (NF), antiferroelectric nematic (NX), and ordinary nematic (N), so the phase sequence is finely tunable.","Interaction picture: since offset π-π stacking rather than point-dipole attraction dominates the lateral interactions, design should optimise surface-charge complementarity rather than only adding polar groups.","Dipole de-emphasis: among chemically similar compounds, longitudinal dipole magnitude is not a reliable guide to polar-phase stability."],"supporting_citations":[{"why":"Supplies the model in which oscillatory longitudinal electrostatic potential makes parallel side-by-side packing attractive; the paper applies it to all 27 homologues.","marker":"[35]"},{"why":"Introduces the archetypal ferroelectric nematic RM734 whose rational design the present study extends.","marker":"[1]"},{"why":"Reports the archetypal DIO ferroelectric nematic and provides the standard for polarisation measurement and phase identification.","marker":"[2]"},{"why":"Gaussian 16, the software used for the DFT geometries, ESP surfaces, and bimolecular potential-energy scans.","marker":"[42]"},{"why":"Provides the Interaction Region Indicator method used to identify offset π-π stacking as the dominant non-covalent interaction.","marker":"[46]"},{"why":"Osipov's analysis of dipole-dipole interactions is cited to show that purely electrostatic surface-charge models cannot fully account for NF formation, which the design rule must respect.","marker":"[25]"}],"fun_headline_variants":["Fluorine pattern, not dipole size, stabilizes polar nematics","Surface charge rhythm, not dipole, sets polar nematic stability","Oscillating surface charge predicts robust ferroelectric nematics","Uniform charge waves, not max fluorination, win in polar order"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The design rule rests on the assumption that a single molecule's computed surface charge, averaged around its long axis as if the molecule spins freely and ignores its neighbours, temperature, and shape changes, captures what actually stabilises polar order in the dense liquid.","fun_headline_variants_meta":{"raw":{"variants":["Fluorine pattern, not dipole size, stabilizes polar nematics","Surface charge rhythm, not dipole, sets polar nematic stability","Oscillating surface charge predicts robust ferroelectric nematics","Uniform charge waves, not max fluorination, win in polar order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000335,"raw_usage":{"total_tokens":1910,"prompt_tokens":1053,"completion_tokens":857,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":669,"completion_tokens_details":{"reasoning_tokens":782}},"tokens_in":669,"tokens_out":857,"duration_ms":8305,"temperature":1.0,"reasoning_tokens":782,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:30:58.253179+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Prepare a 2.2.1 analogue in which the outer-ring fluorine is moved from the position adjacent to the ester carbonyl to the opposite side of the benzoate ring; the spatial-uniformity model predicts a drop in T_NF-I, while a pure dipole-moment model predicts little change, so the comparison would separate the two mechanisms. Alternatively, compute the 1D ESP for known ferroelectric nematogens outside this scaffold and check whether the oscillation-uniformity ranking tracks which compounds actually form polar phases.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the model in which oscillatory longitudinal electrostatic potential makes parallel side-by-side packing attractive; the paper applies it to all 27 homologues."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the archetypal ferroelectric nematic RM734 whose rational design the present study extends."},{"cited_title":"Nishikawa, K","cited_arxiv_id":null,"evidence_quote":"Reports the archetypal DIO ferroelectric nematic and provides the standard for polarisation measurement and phase identification."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gaussian 16, the software used for the DFT geometries, ESP surfaces, and bimolecular potential-energy scans."},{"cited_title":"Lu and Q","cited_arxiv_id":null,"evidence_quote":"Provides the Interaction Region Indicator method used to identify offset π-π stacking as the dominant non-covalent interaction."}],"review_version":1}