pith. sign in

arxiv: 2604.16777 · v1 · submitted 2026-04-18 · 🧮 math.NA · cs.NA

Generalized Scalar Auxiliary Variable Exponential Integrator for A Modified Landau-de Gennes Theory for Smectic Liquid Crystals

Wenshuai Hu , Guanghua Ji , Xiao Li This is my paper

Pith reviewed 2026-05-10 07:41 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords smectic liquid crystalsmodified Landau-de Gennes modelenergy stable schemeexponential integratornumerical analysisdefect dynamics
0
0 comments X

The pith

A reformulated exponential integrator scheme ensures unconditional energy stability for simulations of smectic liquid crystals.

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

The paper develops a numerical method for the modified Landau-de Gennes model of smectic-A liquid crystals, where a tensor order parameter for molecular orientation couples to a scalar for positional order. By combining a generalized scalar auxiliary variable approach with exponential integrators and a relaxed correction, the scheme is made unconditionally stable in a modified energy. This reformulation into a quasi-implicit form removes restrictive time step conditions and allows proofs of bounded solutions and optimal accuracy. Readers interested in computational physics would care because it supports reliable long-time modeling of defect structures without artificial constraints on simulation parameters.

Core claim

The central discovery is that reformulating the exponential time differencing discretization of the GSAV-EI scheme into an equivalent quasi-implicit backward Euler-type structure eliminates CFL mesh-ratio conditions, enabling rigorous fully discrete error analysis while proving unconditional energy stability with respect to a modified discrete energy and uniform boundedness of the tensor order parameter Q.

What carries the argument

The generalized scalar auxiliary variable-exponential integrator (GSAV-EI) with relaxed correction strategy, which reformulates the time stepping to support stability and error analysis for the coupled Q-u system.

If this is right

  • The numerical scheme preserves a modified discrete energy for any time step size.
  • The computed tensor order parameter Q remains uniformly bounded in the simulations.
  • Optimal convergence rates in both time and space are achieved by the method.
  • Complex topological defect dynamics in smectic phases can be captured accurately and efficiently.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar reformulations could extend to other exponential integrator methods for stiff PDE systems in materials science.
  • The approach may enable parameter studies of phase transitions in liquid crystals over extended time scales.
  • Applications to related models in soft matter, such as nematic or cholesteric phases, might benefit from the same stability properties.

Load-bearing premise

The equivalence between the exponential time differencing discretization and the quasi-implicit backward Euler structure holds without introducing additional errors or instabilities.

What would settle it

Numerical experiments demonstrating that the discrete energy increases over time for sufficiently large time steps, or that the error does not converge optimally, would contradict the stability and accuracy claims.

Figures

Figures reproduced from arXiv: 2604.16777 by Guanghua Ji, Wenshuai Hu, Xiao Li.

Figure 1
Figure 1. Figure 1: Evolution of the smectic-A liquid crystal system. The snapshots of the [PITH_FULL_IMAGE:figures/full_fig_p037_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Temporal evolution of the supremum norm (left) and the total energy (right) for the liquid crystal system. The [PITH_FULL_IMAGE:figures/full_fig_p038_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Steady-state profiles of the density field [PITH_FULL_IMAGE:figures/full_fig_p039_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Temporal evolution of the supremum norm (left) and the total energy (right) for the smectic-A liquid crystal system [PITH_FULL_IMAGE:figures/full_fig_p040_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Dynamic evolution of the smectic-A liquid crystal system. The top row displays the density field [PITH_FULL_IMAGE:figures/full_fig_p041_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Performance comparison of different time-stepping strategies and numerical schemes. The left panel shows the [PITH_FULL_IMAGE:figures/full_fig_p042_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of numerical stability and energy dissipation between the SAV-enhanced schemes and the standard [PITH_FULL_IMAGE:figures/full_fig_p043_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Profiles of the adaptive time step sizes [PITH_FULL_IMAGE:figures/full_fig_p044_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Three-dimensional dynamic evolution of the smectic-A phase at representative time instances [PITH_FULL_IMAGE:figures/full_fig_p045_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: provides a rigorous quantitative validation of the temporal accuracy and the strict energy 0 20 40 60 80 100 Time 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Supremum norm sESAV-EI with ==1e-2 sESAV-EI with ==1e-4 0 20 40 60 80 100 Time -150 -100 -50 0 50 100 150 Total Energy Modified energy sESAV-EI with ==1e-2 Original energy sESAV-EI with ==1e-2 Modified energy sESAV-EI with ==1e-4 Original energy sESAV-EI with ==1e… view at source ↗
read the original abstract

The Smectic-A (SmA) phase is modeled by a modified Landau-de Gennes (mLdG) model proposed by Xia et al. [Phys. Rev. Lett., 126 (2021), 177801], in which a tensor order parameter $\mathbf{Q}$ for the orientational order is coupled with a real scalar $u$ characterizing the positional order. In this paper, we propose and analyze a novel, highly efficient, and unconditionally energy-stable numerical scheme for this coupled system by combining the generalized scalar auxiliary variable-exponential integrator (GSAV-EI) approach with a relaxed correction strategy. In particular, we reformulate the exponential time differencing time discretization into an equivalent quasi-implicit backward Euler-type structure, a pivotal step that eliminates the restrictive CFL mesh-ratio conditions of the original GSAV-EI method and enables a rigorous fully discrete error analysis. Theoretically, we rigorously establish the unconditional energy stability with respect to a modified discrete energy and the uniform boundedness of the numerical solutions $\mathbf{Q}$, along with optimal error estimates in both time and space. Comprehensive numerical experiments are presented to demonstrate the accuracy, efficiency, and structural preservation of the algorithm, as well as its capability in capturing complex topological defect dynamics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript proposes and analyzes a numerical scheme for the modified Landau-de Gennes (mLdG) model of smectic liquid crystals, which couples a tensor order parameter Q with a scalar positional order parameter u. It combines the generalized scalar auxiliary variable (GSAV) approach with exponential integrators (EI) and a relaxed correction strategy. The central technical step is a reformulation of the exponential time differencing discretization into an equivalent quasi-implicit backward-Euler structure that removes CFL mesh-ratio restrictions. The authors claim to prove unconditional energy stability with respect to a modified discrete energy, uniform boundedness of the numerical solutions Q, and optimal error estimates in both time and space, supported by numerical experiments demonstrating accuracy, efficiency, and defect dynamics.

Significance. If the algebraic equivalence of the reformulation holds rigorously and the stability/error proofs are complete, the work supplies a practical, unconditionally stable, and theoretically grounded integrator for a coupled tensor-scalar system arising in liquid-crystal modeling. Such schemes are valuable for large-scale simulations of topological defects where explicit time-step restrictions would otherwise be prohibitive.

major comments (1)
  1. The pivotal reformulation step (described in the abstract) that converts the GSAV-EI exponential time differencing into an equivalent quasi-implicit backward-Euler structure must be shown to remain algebraically exact once the relaxed correction and the nonlinear Q-u coupling are included. The relaxation parameter introduces a perturbation whose effect on the discrete energy law and on the local truncation error needs to be controlled uniformly with respect to the mesh ratio; without an explicit derivation or bound, both the unconditional stability claim and the optimal error estimates rest on an unverified equivalence.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for identifying the need for greater clarity on the central reformulation step. We address this point directly below.

read point-by-point responses

  1. Referee: The pivotal reformulation step (described in the abstract) that converts the GSAV-EI exponential time differencing into an equivalent quasi-implicit backward-Euler structure must be shown to remain algebraically exact once the relaxed correction and the nonlinear Q-u coupling are included. The relaxation parameter introduces a perturbation whose effect on the discrete energy law and on the local truncation error needs to be controlled uniformly with respect to the mesh ratio; without an explicit derivation or bound, both the unconditional stability claim and the optimal error estimates rest on an unverified equivalence.

    Authors: We thank the referee for this observation. In Section 3.2 we derive the reformulation for the fully coupled system, showing that the GSAV-EI discretization with relaxed correction is algebraically equivalent to a quasi-implicit backward-Euler scheme (see the chain of equalities leading to (3.15)). The equivalence is exact because the auxiliary variable update and the relaxation are applied after the linear exponential integrator step and do not alter the implicit treatment of the nonlinear terms. The perturbation induced by the relaxation parameter is absorbed into the modified discrete energy; Theorem 3.1 establishes unconditional stability for any fixed relaxation parameter in (0,1) without reference to the mesh ratio. For the error analysis, Lemma 4.1 bounds the local truncation error of the relaxed scheme by a term that is O(Δt) uniformly in the mesh ratio, since the exponential integrator treats the stiff linear part exactly. We will add an explicit lemma in the revised manuscript that isolates the algebraic equivalence and the uniform bound on the relaxation perturbation, thereby making the argument self-contained. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained via standard reformulation and analysis

full rationale

The paper's central steps consist of combining the established GSAV-EI framework with a relaxed correction and an algebraic reformulation of the exponential integrator into a quasi-implicit backward-Euler structure. This reformulation is presented as an equivalence that removes CFL restrictions and permits standard energy estimates and error analysis on the modified discrete energy. No load-bearing step reduces by definition to a fitted parameter, self-referential quantity, or prior self-citation chain; the stability and boundedness claims follow from the reformulated scheme's structure rather than tautological inputs. External citations (e.g., the mLdG model) are to independent prior work and do not carry the uniqueness or ansatz burden. The derivation therefore remains non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated. The mLdG model itself is taken from prior literature (Xia et al.).

pith-pipeline@v0.9.0 · 5530 in / 1234 out tokens · 50699 ms · 2026-05-10T07:41:20.111081+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    J. M. Ball and A. Zarnescu, Orientability and energy minimization in liquid crystal models, Archive for rational mechanics and analysis, 202 (2011), pp. 493–535

    work page 2011

  2. [2]

    Baskaran, J

    A. Baskaran, J. S. Lowengrub, C. W ang, and S. M. Wise, Convergence analysis of a second orderconvexsplittingschemeforthemodifiedphasefieldcrystalequation, SIAMJournalonNumerical Analysis, 51 (2013), pp. 2851–2873

    work page 2013

  3. [3]

    Biscari, M

    P. Biscari, M. C. Calderer, and E. M. Terentjev, Landau–degennes theory of isotropic-nematic-smectic liquid crystal transitions, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 75 (2007), p. 051707

    work page 2007

  4. [4]

    Chen and T

    J.-h. Chen and T. Lubensky, Landau-ginzburg mean-field theory for the nematic to smectic-c and nematic to smectic-a phase transitions, Physical Review A, 14 (1976), p. 1202

    work page 1976

  5. [5]

    Cheng and J

    Q. Cheng and J. Shen, Multiple scalar auxiliary variable (msav) approach and its application to the phase-fieldvesiclemembranemodel, SIAM Journal on Scientific Computing, 40 (2018), pp. A3982– A4006

    work page 2018

  6. [6]

    Q. Du, L. Ju, X. Li, and Z. Qiao, Stabilized linear semi-implicit schemes for the nonlocal cahn–hilliardequation, Journal of Computational Physics, 363 (2018), pp. 39–54

    work page 2018

  7. [7]

    , Maximumprinciplepreservingexponentialtimedifferencingschemesforthenonlocal allen–cahn equation, SIAM Journal on Numerical Analysis, 57 (2019), pp. 875–898. 47

    work page 2019

  8. [8]

    Q. Du, L. Ju, X. Li, and Z. Qiao, Maximum bound principles for a class of semilinear parabolic equations and exponential time-differencing schemes, SIAM Review, 63 (2021), pp. 317–359

    work page 2021

  9. [9]

    D. J. Eyre, Unconditionally gradient stable time marching the cahn-hilliard equation, MRS online proceedings library (OPL), 529 (1998), p. 39

    work page 1998

  10. [10]

    M. Fei, W. W ang, P. Zhang, and Z. Zhang, On the isotropic–nematicphase transition for the liquid crystal, Peking Mathematical Journal, 1 (2018), pp. 141–219

    work page 2018

  11. [11]

    X. Feng, T. Tang, and J. Yang, Stabilized crank-nicolson/adams-bashforth schemes for phase field models, East Asian Journal on Applied Mathematics, 3 (2013), pp. 59–80

    work page 2013

  12. [12]

    J. Han, Y. Luo, W. W ang, P. Zhang, and Z. Zhang, Frommicroscopic theory to macroscopic theory: a systematic study on modeling for liquid crystals, Archive for rational mechanics and analysis, 215 (2015), pp. 741–809

    work page 2015

  13. [13]

    Y. Hu, Y. Qu, and P. Zhang, On the disclination lines of nematic liquid crystals, Communications in Computational Physics, 19 (2016), pp. 354–379

    work page 2016

  14. [14]

    Huang and S

    J. Huang and S. Ding, Global well-posedness for the dynamical Q-tensor model of liquid crystals, Science China Mathematics, 58 (2015), pp. 1349–1366

    work page 2015

  15. [15]

    Izzo and M

    D. Izzo and M. J. De Oliveira, Landau theory for isotropic, nematic, smectic-a, and smectic-c phases, Liquid Crystals, 47 (2020), pp. 99–105

    work page 2020

  16. [16]

    Jiang, Z

    M. Jiang, Z. Zhang, and J. Zhao, Improving the accuracy and consistency of the scalar auxiliary variable (sav) method with relaxation, Journal of Computational Physics, 456 (2022), p. 110954

    work page 2022

  17. [17]

    L. Ju, X. Li, and Z. Qiao, Generalized sav-exponential integrator schemes for allen–cahntype gradient flows, SIAM journal on numerical analysis, 60 (2022), pp. 1905–1931

    work page 2022

  18. [18]

    C. Liu, J. Shen, and X. Yang, Dynamics of defect motion in nematic liquid crystal flow: modeling and numerical simulation, Communications in Computational Physics, 2 (2007), pp. 1184–1198

    work page 2007

  19. [19]

    Y. Liu, C. Quan, and D. W ang, On the maximum bound principle and energy dissipation of exponential time differencing methods for the matrix-valued allen–cahnequation, IMA Journal of Numerical Analysis, 45 (2025), pp. 3342–3377

    work page 2025

  20. [20]

    Z. Liu, Y. Zhang, and X. Li, A novelenergy-optimized techniqueof sav-based(eop-sav)approaches for dissipative systems, Journal of Scientific Computing, 101 (2024), p. 38

    work page 2024

  21. [21]

    Majumdar, Equilibrium order parameters of nematic liquid crystals in the landau-de gennes theory, European Journal of Applied Mathematics, 21 (2010), pp

    A. Majumdar, Equilibrium order parameters of nematic liquid crystals in the landau-de gennes theory, European Journal of Applied Mathematics, 21 (2010), pp. 181–203

    work page 2010

  22. [22]

    W. L. McMillan, Simple molecular model for the smectic a phase of liquid crystals, Physical Review A, 4 (1971), p. 1238. 48

    work page 1971

  23. [23]

    Pazy, Semigroups of linear operators and applications to partial differential equations, Springer Science & Business Media, 2012

    A. Pazy, Semigroups of linear operators and applications to partial differential equations, Springer Science & Business Media, 2012

    work page 2012

  24. [24]

    M. Y. Pevnyi, J. V. Selinger, and T. J. Sluckin, Modeling smectic layersin confinedgeometries: Order parameter and defects, Physical Review E, 90 (2014), p. 032507

    work page 2014

  25. [25]

    J. Shen, J. Xu, and J. Yang, The scalar auxiliary variable (sav) approach for gradient flows, Journal of Computational Physics, 353 (2018), pp. 407–416

    work page 2018

  26. [26]

    , A new class of efficientand robust energy stable schemes for gradient flows, SIAM Review, 61 (2019), pp. 474–506

    work page 2019

  27. [27]

    Shen and X

    J. Shen and X. Yang, Numerical approximationsof allen-cahn and cahn-hilliard equations, Discrete Contin. Dyn. Syst, 28 (2010), pp. 1669–1691

    work page 2010

  28. [28]

    B. Shi, Y. Han, C. Ma, A. Majumdar, and L. Zhang, A modified landau–degennes theory for smectic liquid crystals: phase transitions and structural transitions, SIAM Journal on Applied Mathematics, 85 (2025), pp. 821–847

    work page 2025

  29. [29]

    W ang, L

    W. W ang, L. Zhang, and P. Zhang, Modelling and computation of liquid crystals, Acta Numerica, 30 (2021), pp. 765–851

    work page 2021

  30. [30]

    Xia and P

    J. Xia and P. E. F arrell, Variational and numerical analysis of a q-tensor model for smectic-a liquid crystals, ESAIM: Mathematical Modelling and Numerical Analysis, 57 (2023), pp. 693–716

    work page 2023

  31. [31]

    Xia and Y

    J. Xia and Y. Han, Simple tensorial theory of smectic c liquid crystals, Physical Review Research, 6 (2024), p. 033232

    work page 2024

  32. [32]

    J. Xia, S. MacLachlan, T. J. Atherton, and P. E. F arrell, Structural landscapes in geometrically frustrated smectics, Physical review letters, 126 (2021), p. 177801

    work page 2021

  33. [33]

    Xu and T

    C. Xu and T. Tang, Stability analysis of large time-stepping methods for epitaxial growth models, SIAM Journal on Numerical Analysis, 44 (2006), pp. 1759–1779

    work page 2006

  34. [34]

    Z. Xu, X. Yang, H. Zhang, and Z. Xie, Efficientand linear schemes for anisotropic cahn–hilliard model using the stabilized-invariantenergy quadratization (s-ieq) approach, Computer Physics Com- munications, 238 (2019), pp. 36–49

    work page 2019

  35. [35]

    Yang and D

    X. Yang and D. Han, Linearly first-and second-order, unconditionally energy stable schemes for the phase field crystal model, Journal of Computational Physics, 330 (2017), pp. 1116–1134

    work page 2017

  36. [36]

    Yang and G.-D

    X. Yang and G.-D. Zhang, Convergence analysis for the invariant energy quadratization (ieq) schemes for solving the cahn–hilliardand allen–cahnequations with general nonlinear potential, Jour- nal of scientific computing, 82 (2020), p. 55. 49

    work page 2020

  37. [37]

    Zhang, C

    B. Zhang, C. Zhou, and H. Fu, A novel efficient generalized energy-optimized exponential sav schemewithvariable-stepbdfk method forgradientflows, Applied Numerical Mathematics, 210 (2025), pp. 39–63

    work page 2025

  38. [38]

    Zhang and J

    Y. Zhang and J. Shen, A generalized sav approach with relaxation for dissipative systems, Journal of Computational Physics, 464 (2022), p. 111311

    work page 2022

  39. [39]

    J. Zhao, X. Yang, Y. Gong, and Q. W ang, A novel linear second order unconditionally energy stable scheme for a hydrodynamic q-tensor model of liquid crystals, Computer Methods in Applied Mechanics and Engineering, 318 (2017), pp. 803–825. 50

    work page 2017