REVIEW 3 major objections 5 minor 48 references
Patterned states in the nematic phase of flexible-core phenyl benzoate dimers
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read 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…
desk verdict Rich experimental phenomenology, but the central SGS model's quantitative claim is contradicted by its own arithmetic. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [III.B, Eqs. (5)–(7)] 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.
- [III.B, after Eq. (4)] 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.
- [III.F, Eq. (9) and Fig. 20] 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.
minor comments (5)
- [III.B, text near Eq. (5)] The flexoelectric coefficient e3 has units C/m; the sentence 'e3 is about 10 pC/m 2' should read '10 pC/m'.
- [III.A, penultimate paragraph] The phrase 'twit-bend nematogen' contains a typo and should be 'twist-bend nematogen'.
- [Fig. 18 caption] 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.
- [III.F, Fig. 20] 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.
- [III.F, quotation from Ref. [40]] 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).
Circularity Check
No circular derivation: the SGS model's material constants are external measurements and the INR test applies the external Pikin–Chigrinov formula; only non-load-bearing self-citations appear, so the score is a low 2. A separate, non-circular arithmetic inconsistency undermines the claimed dq≈10 crossing.
full rationale
No circular step is exhibited. For the striped ground state, k33≈1 pN, k22≈1.5 pN and e3≈10 pC/m are adopted from independent published measurements on other twist-bend dimers (Refs. 27-30), and Eqs. (5)-(7) are used to solve E_fluct=0 for dq; the observed stripe wavelength is not inserted as an adjustable parameter, so the comparison λ≈0.63d versus the observed value is an externally constrained prediction rather than a fitted parameter renamed as a result. For the very high frequency mode, Eq. (9) is the external Pikin-Chigrinov inertial-conduction formula, and the slope of UNR versus f√(ε⊥/σ⊥) is used to extract an effective K≈1 pN as an output; the target pattern is not defined by that fit. Self-citations (e.g., Refs. 17-19, 31, 33, 37, 41) are contextual and do not carry the derivation, so they do not make the argument circular. I flag two issues that affect correctness rather than circularity: first, substituting the paper's stated constants into Eqs. (5)-(7) gives a zero crossing near dq≈4, not the claimed dq≈10, so the model-observation agreement in Sec. III.B is arithmetically unsupported as written; second, the neglect of surface-charge depolarization (Sec. III.B, just before Eq. 5) rests on an unvalidated screening assumption. Neither issue is a self-consistent reduction of the prediction to its inputs.
Assumptions & free parameters
free parameters (1)
- K (effective elastic constant) =
≈1 pN
assumptions (5)
- domain assumption Flexoelectric and elastic constants from CB7CB-like dimers (e3≈10 pC/m², k33≈1 pN, k22≈1.5 pN) apply to 5O9 and 6O9.
- ad hoc to paper Surface charge contributions to the flexoelectric energy can be ignored because neighboring half-waves have opposite surface charge densities screened by free ions.
- domain assumption The director distortion in the striped ground state is described by φ(x,z)=φ0 cos(πz/d) sin(qx) with no y-dependence.
- domain assumption σ⊥ can be used in place of both σ and (σ||-σ⊥), and ε=ε⊥, in Eq. (9).
- standard math Standard nematic elasticity and flexoelectricity formalisms apply.
Cite this review
Pith. "Pith review of Patterned states in the nematic phase of flexible-core phenyl benzoate dimers." pith.science (2026). https://pith.science/paper/UBWRJTVA
@misc{pith2026260811067,
author = {Pith},
title = {Pith review of: Patterned states in the nematic phase of flexible-core phenyl benzoate dimers},
year = {2026},
howpublished = {\url{https://pith.science/paper/UBWRJTVA}},
note = {Machine review of arXiv:2608.11067}
}
read the original abstract
This study deals with both spontaneously-formed and electrically-induced structures observed in the nematic phase of two dielectrically negative twist-bend nematogens. In planar cells, the nematic layers are inhomogeneous, existing in a quasiperiodic ground state that involves essentially azimuthal director deviations. This phenomenon is understood using a simple model based on the relative flexoelectric and elastic contributions to free-energy. In an external electric field, a variety of instabilities are obtained. In an increasing static field, for example, the initial periodic surface electroconvective instability is followed by the volume periodic flexoelectric instability. Uncommonly, the growth of the latter takes place via nucleation and front propagation. Due to the low bend elastic deformation cost, flexoelectric bands progressively distort as they narrow (mediated by edge dislocations) under an increasing field to eventually form fanlike objects. In the megahertz region, in each half cycle following 0 V, the electroconvective and flexoelectric instabilities appear transiently with the latter setting in at a higher voltage. Above a few hertz, flexoelectric instability ceases; the sequence of patterned states then is oblique roll to bimodal-grid to normal-roll to chevron. Above 100 kHz, periodic wide bands oriented normal to the rubbing axis are obtained at a threshold voltage that decreases with increasing frequency. Our measurement of threshold voltage corresponding to different frequencies (or electrical conductivity and permittivity values) in this regime agrees well with the earlier theoretical predictions relating to the inertial conduction instability.
Reference graph
Works this paper leans on
-
[1]
R. B. Meyer, in Molecular Fluids, edited by R. Balian and G. Weil (Gordon and Breach, New York, 1976), pp. 271–373
work page 1976
- [2]
-
[3]
M. Cestari, S. Diez -Berart, D. A. Dunmur , A. Ferrarini, M. R. de la Fuente, D. J. B. Jackson, D. O. Lopez, G. R. Luckhurst, M. A. Perez-Jubindo, R. M. Richardson, J. Salud, B. A. Timimi, and H. Zimmermann, Phys. Rev. E 84, 031704 (2011)
work page 2011
-
[4]
V. P. Panov, M. Nagaraj, J. K. Vij, Yu. P. Panarin, A. Kohlmeier, M. G. Tamba, R. A. Lewis, and G. H. Mehl, Phys. Rev. Lett. 105, 167801 (2010)
work page 2010
-
[5]
P. K. Challa, V. Borshch, O. Parri, C. T. Imrie, S. N. Sprunt, J. T. Gleeson, O. D. Lavrentovich, and A. Jákli, Phys. Rev. E 2014, 89, 060501R (2014)
work page 2014
-
[6]
V. P. Panov, S. P. Sreenilayam, Y. P. Panarin, J. K. Vij, C. J. Welch, and G. H. Mehl, Nano Lett. 12, 7515-7519 (2017)
work page 2017
-
[7]
D. Chen, J. H. Porada, J. B. Hooper, A. Klittnick, Y. Shen, M. R.Tuchband, E. Korblova, D. Bedrov, D. M.Walba, M. A. Glaser, J. E. Maclennan , and N. A. Clark, Proc. Natl. Acad. Sci. USA 110, 15931-15936 (2013)
work page 2013
-
[8]
V. Borshch, Y.-K. Kim, J. Xiang, M. Gao, A. Jákli, V.P. Panov, J.K. Vij, C.T. Imrie, M.G. Tamba, G.H. Mehl, and O.D. Lavrentovich, Nature Comm. 4, 2635 (2013). 30
work page 2013
Show all 48 references
-
[9]
D. Chen, M. Nakata, R. Shao, M. R. Tuchband, M. Shuai, U. Baumeister, W. Weissflog, D. M. Walba, M. A. Glaser, J. E. Maclennan , and N. A. Clark, Phys. Rev. E 89, 022506 (2014)
2014
-
[10]
R. R. R. de Almeida, C. Zhang, O. Parri, S. N. Sprunt, and A. Jákli, Liq. Cryst. 41, 1661- 1667 (2014)
2014
-
[11]
V. P. Panov, R. Balachandran, M. Nagaraj, J. K. Vij, M. G. Tamba, A. Kohlmeier , and G. H. Mehl, Appl. Phys. Lett. 99, 261903 (2011)
2011
-
[12]
V. P. Panov, R. Balachandran, J. K. Vij, M. G. Tamba, A. Kohlmeier , and G. H. Mehl, Appl. Phys. Lett. 101, 234106 (2012)
2012
-
[13]
Meyer, G
C. Meyer, G. R. Luckhurst, and I. Dozov, Phys. Rev. Lett. 89, 067801 (2013)
2013
-
[14]
Xiang, S
J. Xiang, S. V. Shiyanovskii, C. Imrie, and O. D. Lavrentovich, Phys. Rev. Lett . 112, 217801 (2014)
2014
-
[15]
Meyer, Liq
C. Meyer, Liq. Cryst. 43, 2144-2162 (2016)
2016
-
[16]
S. P. Sreenilayam, V. P. Panov, J. K. Vij, and G. Shanker, Liq. Cryst. 44, 244-257 (2017)
2017
-
[17]
K. S. Krishnamurthy, M. B. Kanakala, C. V. Yelamaggada, and M. Kleman, Soft Matter 14, 5393- (2018)
2018
-
[18]
Kleman and K
M. Kleman and K. S. Krishnamurthy, Phys. Rev. E 98, 032705 (2018)
2018
-
[19]
K. S. Krishnamurthy, D. S. Shankar Rao, M. B. Kanakala, C. V. Yelamaggada, and M. Kleman, Soft Matter 16, 7479-7491 (2020)
2020
-
[20]
M. A. Osipov, Liq. Cryst. 51, 948-960 (2024). https://doi.org/10.1080/02678292.2024.2321464
2024
-
[21]
Jákli, O
A. Jákli, O. D. Lavrentovich, and J. V. Selinger, Rev. Mod. Phys. 90, 045004 (2018). https://doi.org/10.1103/RevModPhys.90.045004
2018 doi
-
[22]
V. P. Panov, J-K. Song, G. H. Mehl, and J. K. Vij, Crystals 11, 621 (2021). https://doi.org/10.3390/cryst11060621
2021 doi
-
[23]
M. R. Tuchband, M. Shuai, K. A. Graber, D. Chen, C. Zhu, L. Radzihovsky, A. Klittnick, L. Foley, A. Scarbrough, J. H. Porada, M. Moran, J. Yelk, J. B. Hooper, X. Wei, D. Bedrov, C. Wang, E. Korblova, D. M. Walba, A. Hexemer, J. E. Maclennan, M. A. Glaser, and N. A. Clark, Crys...
2024 doi
-
[24]
Adlem, M
K. Adlem, M. Čopič, G. R. Luckhurst, A. Mertelj, O. Parri, R. M. Richardson, B. D. Snow, B. A. Timimi, R. P. Tuffin and D. Wilkes, Phys. Rev. E 88, 022503 (2013)
2013
-
[25]
C.-J. Yun, M. R. Vengatesan , J. K. Vij and J. -K. Song, Appl. Phys. Lett. 106, 173102 (2015)
2015
-
[26]
B. R. Hernández, N. Sebastián, M. R. de la Fuente, D. O. López, S. Diez-Berart, J. Salud, M. B. Ros, D. A. Dunmur, G. R. Luckhurst and B. A. Timimi, Phys. Rev. E 92, 062505 (2015)
2015
-
[27]
Babakhanova, Z
G. Babakhanova, Z. Parsouzi, S. Paladugu, H. Wang, Yu. A. Nastishin, S. V. Shiyanovskii, S. Sprunt and O. D. Lavrentovich, Phys. Rev. E 96, 062704 (2017)
2017
-
[28]
Balachandran, V
R. Balachandran, V. P. Panov, Yu. P. Panarin, M. G. Tamba, G. H. Mehl, J. K. Song , and J. K. Vij, Mol. Cryst. Liq. Cryst. 611, 65-70 (2015)
2015
-
[29]
Varanytsia and L
A. Varanytsia and L. –C. Chien, Sci. Reports 7, 41333 (2017). DOI: 10.1038/srep41333
2017 doi
-
[30]
Škarabot, N
M. Škarabot, N. J. Mottram, S. Kaur, C. T. Imrie, E. Forsyth, J. M. D. Storey, R. Mazur, W. Piecek, and L. Komitov, ACS Omega 7, 9785−9795 2022 (2022)
2022
-
[31]
K. S. Krishnamurthy, N. B. Palakurthy, and C. V. Yelamaggad, J. Phys. Chem. B 121, 5447–5454 (2017)
2017
-
[32]
Oswald, G
P. Oswald, G. Poy, and K. S. Krishnamurthy, Phys. Rev. E 106, 024705 (2022)
2022
-
[33]
K. S. Krishnamurthy, S. Krishna Prasad , D. S. Shankar Rao , C. V. Yelamaggad, and I. Dozov, Phys. Rev. E 111, 015402 (2025)
2025
-
[34]
Wiant, J
D. Wiant, J. T. Gleeson, N. Eber, K. Fodor-Csorba, A. Jakli, and T. Toth-Katona, Phys. Rev. E 72, 041712 (2005)
2005
-
[35]
L. M. Blinov and V. G. Chigrinov, Electrooptic Effects in Liquid Crystal Materials (Springer, New York, 1996)
1996
-
[36]
R. J. Mandle, Soft Matter 12, 7883 (2016)
2016
-
[37]
K. S. Krishnamurthy, M. B. Kanakala, C. V. Yelamaggad, and N. V. Madhusudana, J. Phys. Chem. B 123, 1423-1231(2019)
2019
-
[38]
P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Clarendon Press, Oxford, 1998)
1998
-
[39]
R. N. Thurston, J. Cheng, R. B. Meyer, and G. D. Boyd, J. Appl. Phys. 56, 263-272 (1984). 32
1984
-
[40]
nonstandard
to hold good qualitatively even if e1 – e3 and εa0, provided |μ|= |εoεak/ (e1–e3)2| <1, so that the domain density L=1/wF is linear in U. In the two dimers studied here, we do find L to increase linearly with voltage (Fig. 11). B C D E 20 µm A P 8.1 V a F H G I B C D E 8.1 V...
-
[41]
S. A. Pikin, Structural Transformations in Liquid Crystals (Gordon and Breach, New York, 1991)
1991
-
[42]
Tadapatri, K
P. Tadapatri, K. S. Krishnamurthy, and W. Weissflog, Soft Matter 8, 1202–1214 (2012)
2012
-
[43]
M. May, W. Schöpf, I. Rehberg, A. Krekhov and Á. Buka, Phys. Rev. E 78, 046215 (2008). DOI: https://doi.org/10.1103/PhysRevE.78.046215
2008 doi
-
[44]
N. Éber, L. O. Palomares, P. Salamon, A. Krekhov, and Á. Buka, Phys. Rev. E 86, 021702 (2012). DOI: https://doi.org/10.1103/PhysRevE.86.021702
2012 doi
-
[45]
Salamon, N
P. Salamon, N. Éber, A. Krekhov, and Á. Buka, Phys. Rev. E 87, 032505 (2013). DOI: https://doi.org/10.1103/PhysRevE.87.032505
2013 doi
-
[46]
A. N. Trufanov, L. M. Blinov, and M. Barnik, Sov. Phys. JETP 51, 314-319 (1980)
1980
-
[47]
S. A. Pikin and V. G. Chigrinov, Sov. Phys. JETP 51, 123-126 (1980)
1980
-
[48]
Paparini and E
S. Paparini and E. G. Virga, Phys. Rev. E 103, 022707 (2021). DOI: https://doi.org/10.1103/PhysRevE.103.022707
2021 doi
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.