Pith. sign in

REVIEW 3 major objections 6 minor 2 cited by

Phase diagram, confining strings, and a new universality class in nematopolar matter

T0 review · 3 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A minimal free-energy model of coupled nematic and polar order predicts a new universality class and string-confined +1/2 topological defects.

desk verdict The phase diagram and string confinement are worth a look; the 'new universality class' claim is not supported by the evidence in the paper. read the letter →

arxiv 2501.04769 v1 pith:SK6UUYBY submitted 2025-01-08 cond-mat.soft cond-mat.stat-mech

Farzan Vafa , Amin Doostmohammadi This is my paper · ORCID
classification cond-mat.softcond-mat.stat-mech
keywords nematopolarmattercoupledorderparametersphasediagramcriticalexponentstricriticalpointtopologicaldefectstringsconfinementactive
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper studies a minimal free-energy model in which a polar order parameter that prefers isotropy and a nematic order parameter that prefers order interact through a term that favors locking them together. The authors claim that the resulting phase diagram has three phases -- isotropic, nematic, and a nematopolar phase with locked polar and nematic directions -- and that the continuous nematic-to-nematopolar transition has critical exponents that do not match standard nematic or polar universality classes. They further show that in the locked phase a pair of same-sign $+1/2$ nematic defects can be connected by a string of constant tension, so the pair is confined at a length set by balancing string tension against defect repulsion. If correct, the model identifies a new universality class for coupled orientational order and a new mechanism for stabilizing topological defects in soft and active matter.

What carries the argument

The central object is the locking relation between the two order parameters: minimizing the free energy with respect to $p$ at fixed $Q$ gives $Q=p^2+A_P/(2A_{NP})$ in the ordered state, so nematic magnitude is slaved to polar magnitude, and the deviation from this relation defines the locking order parameter $\sigma=Q-p^2-A_P/(2A_{NP})$, which vanishes in the nematopolar phase. Near the nematopolar-nematic boundary, the coupling term favors locking while the polar disorder term $A_P|\mathbf{p}|^2$ opposes it, and this frustration localizes in a one-dimensional string of width $w=\sqrt{K_p/A_P}$ connecting two $+1/2$ nematic defects; the string can be viewed as an elongated core of a $+1$ polar defect. The string tension $T=(A_{NP}-A^*_{NP})\sqrt{K_p/A_P}$ is fixed by requiring the string energy to vanish at the critical coupling $A^*_{NP}$, and balancing it against the nematic elastic energy $\frac{\pi}{2}K_Q\ln(L/\ell)$ gives the equilibrium length $\ell\sim \frac{\pi}{2}\sqrt{K_Q^2 A_P/K_p}\,/(A_{NP}-A^*_{NP})$.

What would settle it

Directly measuring the polar order parameter $p$ and its susceptibility in a two-dimensional numerical simulation along the nematopolar-nematic boundary would settle the universality-class claim: if the measured $\beta$ and $\gamma$ differ from $(1/2,1)$ (or from $(1/4,1/2)$ at the tricritical point), the mean-field class is not the physical one. A separate analytic check is to recompute the specific-heat exponent and the scaling relation so that the two are mutually consistent, since the reported values satisfy $\alpha+2\beta+\gamma=2$ rather than $1$.

Watch

Extended reading notes

Core claim

The paper's central claim is that minimal coupling between an isotropic-tendency polar field $\mathbf{p}$ and an ordered nematic tensor $\mathbf{Q}$, with free energy $F = A_{NP}|\mathbf{Q}-\mathbf{P}|^2 + A_P|\mathbf{p}|^2 + A_Q(1-|\mathbf{Q}|^2)^2$ plus gradient terms and $\mathbf{P}=\mathbf{p}\mathbf{p}-\tfrac12|\mathbf{p}|^2\mathbf{1}$, produces a phase diagram with isotropic, nematic, and nematopolar phases, including a triple point at $(A_P,A_{NP})=(\sqrt{2},2)$ and a tricritical point at $(8/\sqrt{27},4/3)$. The continuous nematopolar-nematic transition is asserted to have its own universality class: for the polar order parameter the mean-field exponents are $(\alpha,\beta,\gamma)=(0,1/2,1)$, shifting to $(1,1/4,1/2)$ at the tricritical point, while the locking order parameter $\sigma=Q-p^2-A_P/(2A_{NP})$ has exponents $(0,1,0)$. In the strongly coupled phase, a pair of same-sign $+1/2$ nematic defects is predicted to be confined by a string with tension $T=(A_{NP}-A^*_{NP})\sqrt{K_p/A_P}$ and equilibrium length $\ell\sim 1/(A_{NP}-A^*_{NP})$; the authors report quantitative agreement with numerical minimization of the free energy.

Load-bearing premise

The load-bearing premise is that the simplified mean-field calculation gives the true critical exponents in two dimensions, with fluctuations neglected; the derivation also uses the scaling relation $\alpha+2\beta+\gamma=1$, while the listed exponents satisfy $\alpha+2\beta+\gamma=2$.

Editorial extensions

If this is right

  • The continuous nematopolar-nematic transition is claimed to be a distinct mean-field universality class: for the polar order parameter the exponents are $(0,1/2,1)$ away from the tricritical point and $(1,1/4,1/2)$ at it, while for the locking order parameter $\sigma$ they are $(0,1,0)$.
  • Same-sign $+1/2$ nematic defect pairs in the locked phase are predicted to be stably confined, with equilibrium separation $\ell$ that diverges as the coupling approaches the critical value $A^*_{NP}$.
  • The same mechanism makes neutral (opposite-sign) $+1/2$ pairs unstable: both the Coulomb-like repulsion and the string tension are attractive for them, so only same-sign pairs can balance to a finite length.
  • A sufficiently strong external field coupled to the polar order drives all topological defects out of the bulk, providing a way to clear defects from a region.
  • The predicted phase diagram includes exact locations for the triple point, the tricritical point, and the continuous boundary $A_P=2A_{NP}\sqrt{1-A_{NP}/2}$, all matched by numerical simulations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Our inference: because the derivation quotes the scaling relation $\alpha+2\beta+\gamma=1$ while the listed exponents satisfy $\alpha+2\beta+\gamma=2$, the claimed new universality class should be checked by direct numerical measurement of the exponents in two dimensions, where fluctuations may shift them.
  • Our inference: the same locking-frustration mechanism could produce confining strings in other coupled orientational systems, such as hexatic-nematic or ferroelectric-nematic materials, whenever one field's topological charges can split into smaller charges of the other field.
  • Our inference: the string-length formula $\ell\sim 1/(A_{NP}-A^*_{NP})$ is a tunable prediction; in an experimental realization, varying the effective coupling through temperature, concentration, or activity should produce a visibly growing defect separation as the transition is approached.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript introduces a Landau free energy (Eqs. 1-3) coupling a polar order parameter p and a nematic order parameter Q, where the polar sector prefers disorder and the nematic sector prefers order. From this free energy the authors derive a phase diagram with isotropic, nematic, and nematopolar phases, compute mean-field critical exponents for the N-I and N-NP transitions (Table I), and propose that in the locked nematopolar phase a pair of same-sign +1/2 nematic defects is confined by a string whose equilibrium length diverges as the inverse distance to the N-NP boundary, Eq. (8). They also study the effect of external fields on phases and defects. Relaxational simulations are used to map phases and to illustrate defect confinement and expulsion.

Significance. If the universality-class and string-confinement claims were established, the paper would be a valuable contribution to soft-matter physics, since coupled polar-nematic order is relevant to ferroelectric nematics, active matter, and biological tissues. The Landau phase-diagram calculation is competently executed and appears to agree with the simulated phase map in Fig. 2. The proposed string mechanism, with repulsive Coulombic defect interaction balanced by a constant tension, is physically appealing and clearly presented. However, the advertised central novelty is not supported: the exponents are zero-temperature mean-field Landau exponents, presented without any fluctuation analysis or numerical measurement, and the derivation of the string tension is an ansatz rather than a calculation from the model. The paper's strengths are therefore concentrated in the phase-diagram part, while the two headline claims require substantial additional support.

major comments (3)
  1. [Critical exponents and a new universality class; Table I; Eq. (1)] The central claim of a new universality class is not established. The exponents in Table I are derived entirely from a zero-temperature Landau expansion (Eq. (3) and the Supplemental Material), and the paper contains no renormalization-group analysis, fluctuation calculation, or numerical measurement of critical exponents. In d=2, the gradient terms in Eq. (1) place the model in a coupled O(2)/RP(1) class whose physical critical behavior is generally fluctuation-dominated rather than Gaussian or mean-field, so the mean-field exponents cannot be identified with the universality class without further argument. Moreover, the generic p exponents (alpha=0, beta=1/2, gamma=1) in Table I are identical to the N-I Q exponents in the same table, and the tricritical values are standard mean-field tricritical values; the assertion that these constitutes a scaling 'distinct from standard nematic or polar universality classes' is therefore not supported by the table itself.
  2. [Supplemental Material, 'Critical exponents', Eqs. (S.34)-(S.41)] The derivation of the gamma values is internally inconsistent. The Supplement invokes 'the hyperscaling relation alpha + 2 beta + gamma = 1' three times, but the correct scaling relation is alpha + 2 beta + gamma = 2, and the reported exponents satisfy the latter (0+1+1=2, 1+1/2+1/2=2, and 0+2+0=2). With the correct relation, the stated gamma values are not implied by alpha and beta. The gamma column of Table I is thus asserted rather than derived; a direct computation of the susceptibility from the h-dependent free energy is required.
  3. [Confining strings, Eqs. (5)-(8)] The inverse-length scaling in Eq. (8) is built into the ansatz rather than derived. After writing the nematic elastic energy as F_Q ~ (pi/2) K_Q ln(L/l), the paper states 'Near the N-NP curve, the string tension must vanish, requiring f(A_NP) ~ A_NP - A*_NP' and then defines the tension in Eq. (6). No evaluation of the polar and coupling free-energy cost of an elongated +1 defect profile is provided, so the linear vanishing of the tension near the boundary is an input assumption, not a prediction of Eq. (1). Consequently Eq. (8) follows from that assumed linear tension. The numerical test in Fig. 3(f) checks only the resulting power law for two values of A_P; it does not independently measure the tension T.
minor comments (6)
  1. [Fig. 1 caption] The caption refers to 'the white point' but the tricritical point is marked as a purple filled circle and the triple point as a black filled circle; please clarify which point is meant.
  2. [Table I and main text] The locking order parameter sigma is used in Table I but is not defined before the table; it is introduced only in the following paragraph. Define sigma in the table caption or in the main text before the table appears.
  3. [Limiting case (i), main text] The sentence 'For small AP (AP < 4 sqrt(6)/9), nematic ordering dominates... defining the nematopolar phase' is misleading: 4 sqrt(6)/9 is the large-ANP asymptote of the first-order boundary, and the condition should be stated as an asymptotic statement for ANP >> 1.
  4. [Supplemental Material, Eq. (S.40)] The chain-rule expression for d^2 F / d A_P^2 appears incomplete: a full second derivative of F with respect to A_P should include additional terms beyond (dQ/dA_P) d^2 F/dQ^2. Please verify the expansion.
  5. [Fig. 3(f)] The fit line is described only by a slope of -1 on a log-log plot; state the fitted exponent and whether the fitted prefactor is consistent with Eq. (8) or only the power law is tested.
  6. [Throughout] There are several typographical and grammatical slips, e.g., 'reachedd boundary' in the Supplementary Material and 'p has inherent tendency' in the main text; a careful proofread is needed.

Circularity Check

1 steps flagged · score 3.0 of 10

One string-length prediction is partly circular (assumed linear tension plus fitted power law); the phase diagram and exponents otherwise derive from the stated free energy, while the 'new universality class' support has internal inconsistencies that are correctness risks, not circularity.

  1. fitted input called prediction [Confining strings, Eqs. (5)-(8) and Fig. 3(f); main text: "Near the N-NP curve, the string tension must vanish..."]
    "Near the N-NP curve, the string tension must vanish, requiring f (AN P) ∼ AN P− A∗ N P(AP ) ... Using w = p Kp/AP , the string tension is T ≡ (AN P− A∗ N P) q Kp AP . ... minimizing ... gives ℓ ∼ π 2 q K 2 QAP Kp 1 AN P−A∗ N P . [Fig. 3 caption:] red dashed line is fit of Eq. (8)."

    The only divergence in Eq. (8) is inherited from the inserted linear-vanishing ansatz T ∝ (ANP−A*NP); 'must vanish' fixes the zero of the tension but not its scaling power. Minimizing (π/2)KQ ln(L/ℓ)+Tℓ then mechanically yields ℓ ∝ 1/(ANP−A*NP), so the inverse-power divergence is equivalent to the assumed linear tension by construction. Moreover the 'quantitative confirmation' is obtained by fitting Eq. (8) to the same simulation data, as the figure caption states, with a fitted prefactor rather than a comparison of the theoretically predicted constant; the functional form is therefore not an independent test of the divergence. This is a partial circularity in the string-length prediction, not in the phase-diagram derivation.

full rationale

The phase diagram is derived self-contained from the Landau free energy of Eq. (2) / Eq. (S.1): the phase boundaries, tricritical point (8/√27, 4/3), triple point (√2, 2), and the continuous boundary Eq. (4) all follow from minimizing F and equating free energies. I find no load-bearing self-citation: references [35] and [36] are background context, and no uniqueness theorem from prior work is invoked to forbid alternatives. The critical exponents are computed from the same free energy via Landau expansions (Supplement, 'Critical exponents'), so the derivation chain does not reduce to a fitted parameter; the claim that the N-NP transition defines a new universality class is weakened, however, by the absence of fluctuation or RG analysis showing that these mean-field exponents are the physical 2D exponents. I also flag an internal inconsistency: the Supplement writes the hyperscaling relation as α+2β+γ=1 and uses it to fix γ, while all three reported triples satisfy α+2β+γ=2; with the standard relation the same γ values follow, so this is an algebraic inconsistency and a correctness risk, not circularity. The one genuine partial circularity is the confining-string length scaling: the tension is assumed to vanish linearly at the N-NP boundary, and the resulting 1/(ANP−A*NP) divergence is then fitted to simulations rather than independently predicted. Because the central phase-diagram and exponent claims are nevertheless derived from the stated free energy rather than from their own outputs, the overall circularity score is 3.

Assumptions & free parameters 1 free parameters · 4 assumptions · 1 invented entities

The central results rest on (i) the specific Landau free energy, (ii) the mean-field uniform treatment of order parameters, (iii) an assumed linear vanishing of string tension near the phase boundary, and (iv) a misstated hyperscaling relation. The only numerically fitted constant is the prefactor 6.5 in the string-length fit; the model control parameters A_P, A_NP, K_p, and K_Q are inputs, not fitted. No new fundamental entity is postulated beyond the defect string, which is a predicted composite structure with no independent experimental evidence.

free parameters (1)
  • string-length prefactor c = 6.5
    Fig. 3(f) fits the predicted scaling as l sqrt(A_P) = 6.5 / (A_NP - A*_NP). The prefactor predicted from K_Q, K_p, and A_P is not independently evaluated, so the quantitative confirmation includes a fitted constant.
assumptions (4)
  • domain assumption The Landau free energy in Eq. (2) is the minimal model for coexisting nematic and polar order.
    The free energy is posited, not derived from a microscopic or active-matter theory; all subsequent results inherit this choice.
  • domain assumption Uniform mean-field order parameters p and Q capture the phase diagram.
    The minimization in the Supplemental Material treats p and Q as homogeneous and drops gradient terms, so fluctuations and defect cores are not included in the phase boundary derivations.
  • ad hoc to paper String tension vanishes linearly as f ~ A_NP - A*_NP near the N-NP boundary.
    The main text states the tension 'must vanish, requiring' this form; it is an ansatz, not derived from the free energy, and it directly produces the 1/(A_NP - A*_NP) length scaling.
  • ad hoc to paper The paper's version of hyperscaling, alpha + 2 beta + gamma = 1, is used to fix gamma.
    Invoked in the Supplemental critical exponents section without derivation or citation; as written it contradicts the reported exponents, which satisfy alpha + 2 beta + gamma = 2.
invented entities (1)
  • Confining string (elongated +1 polar defect core)
    purpose: Explains confinement of same-sign +1/2 nematic defect pairs in the locked nematopolar phase.
    The string is a prediction of the coupled Landau model with an assumed tension ansatz; it is supported only by simulations of the same free energy, not by independent experimental data.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Phase diagram, confining strings, and a new universality class in nematopolar matter." pith.science (2026). https://pith.science/paper/SK6UUYBY

@misc{pith2026250104769,
  author       = {Pith},
  title        = {Pith review of: Phase diagram, confining strings, and a new universality class in nematopolar matter},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SK6UUYBY}},
  note         = {Machine review of arXiv:2501.04769}
}
abstract

We study a minimal model of a system with coexisting nematic and polar orientational orders, where one field tends to order and the other prefers isotropy. For strong coupling, the ordered field aligns the isotropic one, locking their orientations. The phase diagram reveals three distinct phases--nematopolar (aligned orders), nematic (independent orders), and isotropic (vanishing orders)--separated by continuous and discontinuous transitions, including a triple and a tricritical point. We find unique critical scaling for the nematopolar-nematic transition, distinct from standard nematic or polar universality classes. Additionally, in the locked nematopolar phase, we show nematic $+1/2$ topological defect pairs are connected and confined by strings with constant tension. These strings arise from frustration in locking the orientational orders and can be interpreted as elongated cores of $+1$ polar topological defects. When a sufficiently strong background field couples to the polar order, all topological defects are expelled from the region. Analytical predictions are quantitatively confirmed by numerical simulations.

Figures

Figures reproduced from arXiv: 2501.04769 by the authors.

Figure 1
Figure 1. FIG. 1. Phase diagram of a system with coexisting ne [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Phase diagram obtained from numerical simulations. [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Confining strings. (a) Sketch of elongated +1 polar defect of length [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Schematic phase diagrams in the presence of external ←→←→ [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Self-propulsive active nematics

    cond-mat.soft 2025-09 conditional novelty 5.0 of 10

    Self-propulsion in an active nematic model produces a non-monotonic ordering effect, with order, defect anti-hyperuniformity, and long-range vorticity correlations all peaking at an intermediate self-propulsion speed.

  2. The Interplay of Polar and Nematic Order in Active Matter: Implications for Non-Equilibrium Physics and Biology

    cond-mat.soft 2025-06 conditional novelty 2.0 of 10

    A review argues that polar and nematic order frequently coexist in active biological matter, so unified mixed-symmetry models are needed to describe it.

Reference graph

Works this paper leans on

44 extracted references · 40 canonical work pages · cited by 2 Pith papers

  1. [1]

    P. G. de Gennes and J. Prost, The physics of liquid crys- tals (Clarendon Press, Oxford, 1993)

    work page 1993

  2. [2]

    (S.23) and (S.28), respectively

    The lower and upper segments of the solid phase boundary satisfy Eqs. (S.23) and (S.28), respectively. arXiv:2501.04769v1 [cond-mat.soft] 8 Jan 2025 2 30], as well as polar and nematic-like interactions of a single field [31–35]. In this work, we propose a minimal theory that couples both the phase and magnitude of nematic and polar order parameters and c...

    arXiv 2025

  3. [3]

    Yang and S.-T

    D.-K. Yang and S.-T. Wu, Fundamentals of liquid crystal devices (John Wiley & Sons, 2014)

    work page 2014

  4. [4]

    Doostmohammadi and B

    A. Doostmohammadi and B. Ladoux, Physics of liquid crystals in cell biology, Trends in cell biology (2021)

    work page 2021

  5. [5]

    M. C. Marchetti, J.-F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Hydrody- namics of soft active matter, Reviews of modern physics 85, 1143 (2013)

    2013

  6. [6]

    V. M. Kaganer, H. M¨ ohwald, and P. Dutta, Structure and phase transitions in langmuir monolayers, Reviews of Modern Physics 71, 779 (1999)

    work page 1999

  7. [7]

    Schenning, G

    A. Schenning, G. P. Crawford, and D. J. Broer, Liquid crystal sensors (CRC Press, 2017)

    work page 2017

  8. [8]

    S. W. Ula, N. A. Traugutt, R. H. Volpe, R. R. Patel, K. Yu, and C. M. Yakacki, Liquid crystal elastomers: an 6 introduction and review of emerging technologies, Liquid Crystals Reviews 6, 78 (2018)

    work page 2018

Show all 44 references
  1. [9]

    Volfson, S

    D. Volfson, S. Cookson, J. Hasty, and L. S. Tsimring, Biomechanical ordering of dense cell populations, Pro- ceedings of the National Academy of Sciences 105, 15346 (2008)

    2008

  2. [10]

    Doostmohammadi, S

    A. Doostmohammadi, S. P. Thampi, and J. M. Yeomans, Defect-mediated morphologies in growing cell colonies, Physical review letters 117, 048102 (2016)

    2016

  3. [11]

    Kruse, J.-F

    K. Kruse, J.-F. Joanny, F. J¨ ulicher, J. Prost, and K. Seki- moto, Generic theory of active polar gels: a paradigm for cytoskeletal dynamics, The European Physical Journal E 16, 5 (2005)

    2005

  4. [12]

    T. B. Saw, A. Doostmohammadi, V. Nier, L. Kocgozlu, S. Thampi, Y. Toyama, P. Marcq, C. T. Lim, J. M. Yeo- mans, and B. Ladoux, Topological defects in epithelia govern cell death and extrusion, Nature 544, 212–216 (2017)

    2017

  5. [13]

    X. Chen, E. Korblova, D. Dong, X. Wei, R. Shao, L. Radzihovsky, M. A. Glaser, J. E. Maclennan, D. Bedrov, D. M. Walba, et al. , First-principles ex- perimental demonstration of ferroelectricity in a ther- motropic nematic liquid crystal: Polar domains and strik- ing electro-op...

    2020

  6. [14]

    O. D. Lavrentovich, Ferroelectric nematic liquid crys- tal, a century in waiting, Proceedings of the National Academy of Sciences 117, 14629 (2020)

    2020

  7. [15]

    Sebasti´ an, M.ˇCopiˇ c, and A

    N. Sebasti´ an, M.ˇCopiˇ c, and A. Mertelj, Ferroelectric ne- matic liquid-crystalline phases, Physical Review E 106, 021001 (2022)

    2022

  8. [16]

    Eckert, B

    J. Eckert, B. Ladoux, R.-M. M` ege, L. Giomi, and T. Schmidt, Hexanematic crossover in epithelial mono- layers depends on cell adhesion and cell density, Nature Communications 14, 5762 (2023)

    2023

  9. [17]

    Armengol-Collado, L

    J.-M. Armengol-Collado, L. N. Carenza, J. Eckert, D. Krommydas, and L. Giomi, Epithelia are multiscale active liquid crystals, Nature Physics 19, 1773 (2023)

    2023

  10. [18]

    Armengol-Collado, L

    J.-M. Armengol-Collado, L. N. Carenza, and L. Giomi, Hydrodynamics and multiscale orderin confluent epithe- lia, Elife 13, e86400 (2024)

    2024

  11. [19]

    Bruinsma and G

    R. Bruinsma and G. Aeppli, Hexatic order and herring- bone packing in liquid crystals, Physical Review Letters 48, 1625 (1982)

    1982

  12. [20]

    Dierker and R

    S. Dierker and R. Pindak, Dynamics of thin tilted hex- atic liquid crystal films, Physical review letters 59, 1002 (1987)

    1987

  13. [21]

    Sprunt and J

    S. Sprunt and J. Litster, Light-scattering study of bond orientational order in a tilted hexatic liquid-crystal film, Physical review letters 59, 2682 (1987)

    1987

  14. [22]

    J. V. Selinger and D. R. Nelson, Theory of transitions among tilted hexatic phases in liquid crystals, Physical Review A 39, 3135 (1989)

    1989

  15. [23]

    J. V. Selinger, Dynamics of tilted hexatic phases in liquid- crystal films, Journal de Physique II 1, 1363 (1991)

    1991

  16. [24]

    Drouin-Touchette, P

    V. Drouin-Touchette, P. P. Orth, P. Coleman, P. Chan- dra, and T. C. Lubensky, Emergent potts order in a cou- pled hexatic-nematic xy model, Physical Review X 12, 011043 (2022)

    2022

  17. [25]

    Anisimov, E

    M. Anisimov, E. Gorodetski ˘ ı, and V. Zaprudski ˘ ı, Phase transitions with coupled order parameters, Soviet Physics Uspekhi 24, 57 (1981)

    1981

  18. [26]

    R. B. Meyer, Piezoelectric effects in liquid crystals, Phys- ical Review Letters 22, 918 (1969)

    1969

  19. [27]

    S. M. Shamid, S. Dhakal, and J. V. Selinger, Statis- tical mechanics of bend flexoelectricity and the twist- bend phase in bent-core liquid crystals, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics 87, 052503 (2013)

    2013

  20. [28]

    Mertelj, L

    A. Mertelj, L. Cmok, N. Sebasti´ an, R. J. Mandle, R. R. Parker, A. C. Whitwood, J. W. Goodby, and M. ˇCopiˇ c, Splay nematic phase, Physical Review X 8, 041025 (2018)

    2018

  21. [29]

    Sebasti´ an, L

    N. Sebasti´ an, L. Cmok, R. J. Mandle, M. R. de la Fuente, I. Drevenˇ sek Olenik, M.ˇCopiˇ c, and A. Mertelj, Ferroelectric-ferroelastic phase transition in a nematic liquid crystal, Physical review letters124, 037801 (2020)

    2020

  22. [30]

    M. P. Rosseto and J. V. Selinger, Theory of the splay ne- matic phase: single versus double splay, Physical Review E 101, 052707 (2020)

    2020

  23. [31]

    Paik and J

    L. Paik and J. V. Selinger, Flexoelectricity versus elec- trostatics in polar nematic liquid crystals, arXiv preprint arXiv:2408.10347 (2024)

    2024 arXiv

  24. [32]

    S. E. Korshunov, Possible splitting of a phase transition in a 2 d xymodel, JETP Lett. 41, 216 (1985)

    1985

  25. [33]

    Lee and G

    D. Lee and G. Grinstein, Strings in two-dimensional clas- sical xy models, Physical review letters 55, 541 (1985)

    1985

  26. [34]

    Carpenter and J

    D. Carpenter and J. Chalker, The phase diagram of a gen- eralised xy model, Journal of Physics: Condensed Matter 1, 4907 (1989)

    1989

  27. [35]

    Y. Shi, A. Lamacraft, and P. Fendley, Boson pairing and unusual criticality in a generalized xy model, Physical review letters 107, 240601 (2011)

    2011

  28. [36]

    Amiri, R

    A. Amiri, R. Mueller, and A. Doostmohammadi, Uni- fying polar and nematic active matter: emergence and co-existence of half-integer and full-integer topological defects, Journal of Physics A: Mathematical and The- oretical 55, 094002 (2022)

    2022

  29. [37]

    Ruider, K

    I. Ruider, K. Thijssen, D. R. Vannier, V. Paloschi, A. Sciortino, A. Doostmohammadi, and A. R. Bausch, Topological excitations gov- ern ordering kinetics in endothelial cell lay- ers, bioRxiv 10.1101/2024.09.26.615134 (2024), https://www.biorxiv.org/content/early/2024/09/26/20...

    2024 doi

  30. [38]

    W. E. Schiesser, The numerical method of lines: integra- tion of partial differential equations (Elsevier, 2012)

    2012

  31. [39]

    Press, B

    W. Press, B. Flannery, S. Teukolsky, and W. Vetterling, Multistep, multivalue, and predictor-corrector methods, Numerical Recipes in FORTRAN: The Art of Scientific Computing , 740 (1992). 7 Supplementary Material Phase diagram Here we minimize F , given by F = AN P(Q − p2)2 + ...

    1992

  32. [40]

    p = 0 and Q = 0 for AN P≥ 2

  33. [41]

    3 r − AP 8 + i p |D| # , 2 Re

    p = 0 and Q = q 1 − AN P 2 for AN P≤ 2. We now look for non-zero p solution. We have ∂F ∂p = 0 = ⇒ −2AN P(Q − p2) + AP = 0 (S.6) ∂F ∂Q = 0 = ⇒ AN P(Q − p2) + 2Q(Q2 − 1) = 0 (S.7) From the first equation, we learn that Q = p2 + AP 2AN P (S.8) Thus the nematic order Q is locked ...

  34. [42]

    Isotropic phase: p = 0 and Q = 0 for AN P≥ 2

  35. [43]

    Nematic phase: p = 0 and Q = p 1 − AN P/2 for AN P≤ 2

  36. [44]

    Now that we have identified the phases, we find the boundaries

    Nematopolar phase: p = p Q − AP /(2AN P) and Q = 2 Re 3 q − AP 8 + i p |D| , with D = −1/27 + A2 P /64, provided that AP ≤ 8/ √ 27 and Q ≥ AP /(2AN P). Now that we have identified the phases, we find the boundaries. We first determine the boundary between the nematic and nemat...

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.