arxiv: 2605.05619 · v1 · submitted 2026-05-07 · 🧮 math.NA · cs.NA

Recognition: unknown

A semi-generating function approach to the stability of implicit-explicit multistep methods for nonlinear parabolic equations

Hong-lin Liao , Chaoyu Quan , Tao Tang , Tao Zhou

Authors on Pith no claims yet

Pith reviewed 2026-05-08 07:02 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords implicit-explicit multistepunconditional stabilitysemi-generating functionGrenander-Szegő theoremnonlinear paraboliccontrollability intensityToeplitz matrixtime stepping

0 comments

The pith

Semi-generating function approach establishes unconditional stability for IEMS methods on nonlinear parabolic equations when implicit kernels dominate.

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

This paper presents a semi-generating function approach to analyze stability of implicit-explicit multistep methods for nonlinear parabolic equations. The method represents discrete coefficients as complex rational polynomials on the unit circle and uses the Grenander-Szegő theorem to bound convolution kernels. It provides conditions for unconditional stability based on the minimum eigenvalue of implicit kernels being large and the spectral norm of explicit kernels being small. An implicit-explicit controllability intensity is defined to compare methods and guide the design of new high-order schemes.

Core claim

The authors introduce a semi-generating function approach combined with global discrete energy analysis for the stability of general IEMS methods. Inspired by the Grenander-Szegő theorem, they represent the discrete coefficients via three complex rational polynomials on the unit circle. This yields a unified framework establishing unconditional stability provided the minimum eigenvalue of the implicit composite convolution kernels is sufficiently large and the spectral norm bound of the explicit ones is sufficiently small. They also define an implicit-explicit controllability intensity to evaluate and optimize such methods up to eighth order.

What carries the argument

Semi-generating function approach using three complex rational polynomials on the unit circle to apply the Grenander-Szegő theorem to the Toeplitz matrices arising from the implicit and explicit parts of IEMS methods.

Load-bearing premise

The discrete coefficients of IEMS methods admit representation as three complex rational polynomials evaluated on the unit circle so that the Grenander-Szegő theorem yields the required kernel bounds.

What would settle it

Observing instability or error growth in a simulation of a nonlinear parabolic PDE using an IEMS method where the implicit kernel eigenvalue condition and explicit norm condition are satisfied would falsify the unconditional stability result.

read the original abstract

The rigorous stability analysis of high-order implicit-explicit multistep (IEMS) methods for nonlinear parabolic equations by using discrete energy arguments is a long standing open issue due to their non-A-stable property. A novel semi-generating function approach combined with the global discrete energy analysis is suggested to the stability and convergence analysis of general IEMS methods for nonlinear parabolic equations. Inspired from the Grenander-Szeg\"{o} theorem for the Toeplitz matrix, the semi-generating function approach is used to handle the three groups of discrete coefficients via three complex rational polynomials on the unit circle. A unified theoretical framework is then presented to establish the unconditional stability of IEMS methods if the minimum eigenvalue of composite convolution kernels for the implicit part is properly large and the spectral norm bound of composite convolution kernels for the explicit part is properly small. An indicator, called implicit-explicit controllability intensity, is then introduced to evaluate the degree of controllability of implicit part over explicit part. Some of existing IEMS methods, up to the fifth-order time accuracy, are revisited and compared by computing the associated implicit-explicit controllability intensities such that one can choose certain IEMS method or proper parameter to maintain the unconditional stability for a specific nonlinear parabolic model. We also propose a new parameterized class of IEMS methods, up to the eighth-order time accuracy, which satisfy the priori settings of our theory and have a large value of the implicit-explicit controllability intensity by choosing proper parameter so that they would be well suited for a wide class of nonlinear parabolic problems.

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.

Desk Editor's Note private letter to a colleague

The paper recasts stability of high-order IEMS methods for nonlinear parabolic equations as eigenvalue and norm bounds on composite kernels derived from rational symbols, then uses a controllability intensity to build and tune an eighth-order family.

read the letter

The main new piece is the semi-generating function technique that turns the three groups of multistep coefficients into complex rational polynomials evaluated on the unit circle, then invokes the Grenander-Szegő theorem to extract the minimum eigenvalue of the implicit composite kernel and the spectral norm of the explicit one. From those two quantities they define an implicit-explicit controllability intensity and show that existing methods up to order five can be ranked by it, while a one-parameter family can be pushed to order eight with intensity values large enough to meet the stability thresholds they derive. The global discrete energy argument then closes under those kernel conditions, giving unconditional stability for the nonlinear problem when the intensity is favorable. This is a concrete way to compare and design schemes without having to redo the energy estimate from scratch each time. The framework is presented cleanly and the new family is constructed explicitly to satisfy the a-priori settings. The soft spot is whether the symbol bounds survive the nonlinear term once it is moved to the explicit side. The energy inner product must absorb that term into the explicit kernel contribution without leftover cross terms or hidden time-step dependence, and the paper does not spell out the precise estimates that guarantee this absorption for arbitrary nonlinearities. If the kernels remain strictly Toeplitz after discretization and the nonlinearity is Lipschitz in the right norm, the argument should hold, but that step is asserted rather than derived in detail. The work is aimed at people who actually pick or invent time-stepping schemes for nonlinear parabolic models and want a systematic way to check or improve their stability properties. A reader who needs to decide between existing IEMS methods or to build higher-order ones with guaranteed unconditional stability will get usable guidance from the intensity calculations and the parameterized family. It deserves a serious referee because the technique is new for this setting, the claims are checkable, and the open issue it targets is real.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a semi-generating function approach, leveraging the Grenander-Szegő theorem, to analyze the stability of implicit-explicit multistep (IEMS) methods for nonlinear parabolic equations. It claims to provide a unified framework for unconditional stability based on conditions on the minimum eigenvalue of implicit composite convolution kernels and the spectral norm of explicit ones, introduces an 'implicit-explicit controllability intensity' indicator, evaluates existing methods up to order 5, and proposes new parameterized IEMS methods up to order 8 that satisfy the stability conditions for a wide class of problems.

Significance. If the central claims hold, this work would offer a valuable tool for selecting and designing high-order time-stepping methods that are unconditionally stable for nonlinear parabolic PDEs, resolving a long-standing challenge in the field. The controllability intensity provides a quantitative criterion for method choice, and the new family of methods could be practically useful. The creative use of generating functions for discrete coefficients is a strength, though the extension to nonlinear cases requires rigorous verification to ensure the bounds close the energy estimates.

major comments (2)

The unified framework claims unconditional stability once the implicit composite kernel has sufficiently large minimum eigenvalue and the explicit one has sufficiently small spectral norm, with these quantities obtained by evaluating three rational polynomials on the unit circle. The global energy analysis must absorb the nonlinear term (treated explicitly) into the explicit kernel contribution. It is unclear whether the symbol bounds from the Grenander-Szegő theorem provide the necessary coercivity independent of Δt when non-quadratic cross terms from nonlinearity are present in the discrete energy inner products.
The new parameterized class of IEMS methods (up to eighth-order) is stated to satisfy the a priori settings of the theory and achieve large implicit-explicit controllability intensity by choice of parameter. However, without explicit derivations of the composite kernel eigenvalues/norms from the rational symbols or numerical verification on a nonlinear parabolic test problem, it is not shown that the controllability intensity enforces the required dominance for general nonlinear problems.

minor comments (2)

The abstract uses 'priori settings' which should read 'a priori settings'.
The three groups of discrete coefficients are represented as complex rational polynomials; explicit definitions and the precise mapping to the composite convolution kernels would improve clarity and reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, explaining our approach and indicating the revisions we will make to improve clarity and completeness.

read point-by-point responses

Referee: The unified framework claims unconditional stability once the implicit composite kernel has sufficiently large minimum eigenvalue and the explicit one has sufficiently small spectral norm, with these quantities obtained by evaluating three rational polynomials on the unit circle. The global energy analysis must absorb the nonlinear term (treated explicitly) into the explicit kernel contribution. It is unclear whether the symbol bounds from the Grenander-Szegő theorem provide the necessary coercivity independent of Δt when non-quadratic cross terms from nonlinearity are present in the discrete energy inner products.

Authors: We appreciate the referee's focus on the treatment of nonlinearity in the energy analysis. The nonlinear term is treated explicitly and folded into the explicit kernel contribution. The semi-generating function approach, via the Grenander-Szegő theorem applied to the three rational symbols on the unit circle, yields uniform bounds on the minimum eigenvalue of the implicit composite kernel and the spectral norm of the explicit composite kernel. These bounds ensure that the implicit coercivity term dominates in the global discrete energy estimate, producing a positive definite contribution that is independent of Δt. The non-quadratic cross terms arising from the nonlinearity are controlled using the assumed Lipschitz continuity of the nonlinearity together with standard inequalities (e.g., Young's inequality) that absorb them into the implicit coercivity term, provided the implicit-explicit controllability intensity exceeds a suitable threshold. We will revise the manuscript to include a more detailed, step-by-step derivation of this absorption in the energy estimate section. revision: partial
Referee: The new parameterized class of IEMS methods (up to eighth-order) is stated to satisfy the a priori settings of the theory and achieve large implicit-explicit controllability intensity by choice of parameter. However, without explicit derivations of the composite kernel eigenvalues/norms from the rational symbols or numerical verification on a nonlinear parabolic test problem, it is not shown that the controllability intensity enforces the required dominance for general nonlinear problems.

Authors: We agree that additional explicit derivations and verification would strengthen the presentation. The controllability intensity is obtained directly from the minimum eigenvalue and spectral norm of the composite kernels, which are computed by evaluating the three rational symbols on the unit circle. We will add an appendix that derives these quantities explicitly for the parameterized family and shows how the intensity can be made arbitrarily large by suitable parameter choice while preserving the a priori stability conditions. In addition, we will include numerical experiments on a representative nonlinear parabolic problem (e.g., the Allen-Cahn equation) to confirm that the proposed methods remain unconditionally stable for the chosen parameters, thereby verifying that the intensity enforces the required dominance in practice. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the semi-generating function stability framework

full rationale

The derivation applies the external Grenander-Szegő theorem to three rational polynomial symbols obtained directly from the IEMS coefficient groups, yielding explicit min-eigenvalue and spectral-norm bounds on the composite kernels. Unconditional stability is then proved conditionally on those bounds being sufficiently large (implicit) and small (explicit). The controllability intensity is introduced afterward as a scalar indicator computed from the same bounds to compare methods and select parameters; this is an application of the conditional result rather than a definitional reduction. Parameterized new methods are constructed to satisfy the a-priori kernel conditions by direct evaluation of the symbols, with no fitted quantities renamed as predictions and no load-bearing self-citations. The chain is therefore self-contained against the external theorem and the explicit coefficient-to-symbol mapping.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The framework rests on the standard Grenander-Szegő theorem for Toeplitz matrices and introduces one new scalar metric whose value is tuned via a free parameter in the proposed methods.

free parameters (1)

parameter in new IEMS family
Chosen so that the resulting implicit-explicit controllability intensity is large enough to satisfy the stability conditions of the framework.

axioms (1)

standard math Grenander-Szegő theorem applies to the three groups of discrete coefficients represented as complex rational polynomials on the unit circle.
Invoked to obtain bounds on the composite convolution kernels.

invented entities (1)

implicit-explicit controllability intensity no independent evidence
purpose: Scalar indicator that quantifies how strongly the implicit part controls the explicit part.
New quantity introduced to compare existing methods and to guide construction of the new parameterized family.

pith-pipeline@v0.9.0 · 5588 in / 1540 out tokens · 70956 ms · 2026-05-08T07:02:10.905366+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

37 extracted references · 2 canonical work pages

[1]

Akrivis , Implicit-explicit multistep methods for nonlinear parabolic equations , Math

G. Akrivis , Implicit-explicit multistep methods for nonlinear parabolic equations , Math. Comput. , 82 (2013), pp. 45-68

2013
[2]

Akrivis , Stability of implicit-explicit backward difference formulas for nonlinear parabolic equations , SIAM J

G. Akrivis , Stability of implicit-explicit backward difference formulas for nonlinear parabolic equations , SIAM J. Numer. Anal. , 53 (2015), pp. 464-484

2015
[3]

Akrivis, M.H

G. Akrivis, M.H. Chen, F. Yu and Z. Zhou , The energy technique for the six-step BDF method , SIAM J. Numer. Anal. , 59(5) (2021), pp. 2449-2472

2021
[4]

Akrivis, M.H

G. Akrivis, M.H. Chen and F. Yu , The weighted and shifted seven-step BDF method for parabolic equations , IMA J. Numer. Anal. , 2024, doi: 10.1093/imanum/drae089

work page doi:10.1093/imanum/drae089 2024
[5]

Akrivis, M

G. Akrivis, M. Crouzeix and Ch. Makridakis , Implicit-explicit multistep finite element methods for nonlinear parabolic problems , Math. Comput. , 67 (1998), pp. 457-477

1998
[6]

Akrivis, M

G. Akrivis, M. Crouzeix and Ch. Makridakis , Implicit-explicit multistep methods for quasilinear parabolic equations , Numer. Math. , 82 (1999), pp. 521-541

1999
[7]

Akrivis, E

G. Akrivis, E. Katsoprinakis , Backward difference formulae: new multipliers and stability properties for parabolic equations , Math. Comput. , 85 (2016), pp. 2195-2216

2016
[8]

Akrivis and F

G. Akrivis and F. Karakatsani , Modified implicit-explicit BDF methods for nonlinear parabolic equations, BIT Numer. Math. , 43 (2003), pp. 467-483

2003
[9]

Akrivis and C

G. Akrivis and C. Lubich , Fully implicit, linearly implicit and implicit-explicit backward difference formulae for quasi-linear parabolic equations , Numer. Math. , 131 (2015), pp. 713-735

2015
[10]

Contri, B

A. Contri, B. Kov\' a cs, and A. Massing , Error Analysis of BDF 1-6 time-stepping methods for the transient Stokes problem: velocity and pressure estimates, SIAM J. Numer. Anal. , 63:4 (2025), pp. 1586-1616

2025
[11]

C. F. Curtiss and J. O. Hirschfelder , Integration of stiff equations, Proc. National Academy Sciences , 38 (1952), pp. 235-243

1952
[12]

Dahlquist , G-stability is equivalent to A-stability , BIT , 18 (1978), pp

G. Dahlquist , G-stability is equivalent to A-stability , BIT , 18 (1978), pp. 384-401

1978
[13]

C. W. Gear , Numerical initial value problems in ordinary differential equations , Prentice-Hall, Englewood Cliffs, N.J., 1971

1971
[14]

Grenander and G

U. Grenander and G. Szeg\" o , Toeplitz Forms and Their Applications , 2nd edition, AMS Chelsea, Providence, RI, 2001

2001
[15]

J. L. Guermond and J. Shen , A new class of truly consistent splitting schemes for incompressible flows, J. Comput. Phys. , 192 (1) (2003), pp. 262-276

2003
[16]

Hairer, S

E. Hairer, S. P. N rsett and G. Wanner , Solving Ordinary Differential Equations I: Nonstiff Problems , Volume 8 of Springer Series in Computational Mathematics, Second Edition, Springer-Verlag, 1992

1992
[17]

Hairer and G

E. Hairer and G. Wanner , Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , Springer Series in Computational Mathematics Volume 14, Second Edition, Springer-Verlag, 2002

2002
[18]

Huang and J

F. Huang and J. Shen , Stability and error analysis of a second-order consistent splitting scheme for the Navier-Stokes equations , SIAM J. Numer. Anal. , 61(5) (2023), pp. 2408-2433

2023
[19]

Huang and J

F. Huang and J. Shen , On a new class of multistep and IMEX schemes for parabolic type equations , SIAM J. Numer. Anal. , 62(4) (2024), pp. 1609-1637

2024
[20]

Huang and J

F. Huang and J. Shen , Stability and error analysis of a new class of higher-order consistent splitting schemes for the Navier-Stokes equations , Math. Comput. , 2025, doi: 10.1090/mcom/4132

work page doi:10.1090/mcom/4132 2025
[21]

R. A. Horn and C. R. Johnson , Matrix Analysis , Second edition, Cambridge University Press, Cambridge, UK, 2013

2013
[22]

Ji and H.-L

B. Ji and H.-L. Liao , A unified L2 norm error analysis of SAV-BDF schemes for the incompressible Navier-Stokes equation , J. Sci. Comput. , 2024, 100:5, doi: 10.1007/ s10915-024-02555-9

2024
[23]

B. Jin, B. Li, and Z. Zhou , Correction of high-order BDF convolution quadrature for fractional evolution equations, SIAM J. Sci. Comput. , 39(6) (2017), pp. A3129-A3152

2017
[24]

Johnston and J

H. Johnston and J. G. Liu , Accurate, stable and efficient Navier-Stokes solvers based on explicit treatment of the pressure term, J. Comput. Phys. , 199 (2004), pp. 221-259

2004
[25]

B. Li, K. Wang and Z. Zhou . Long-time accurate symmetrized implicit-explicit BDF methods for a class of parabolic equations with non-self-adjoint operators , SIAM J. Numer. Anal. , 58(1) (2020), pp. 189-210

2020
[26]

Li and J

Q. Li and J. Xie . A linear multistep method for solving stiff ordinary differential equations , Journal of Tsinghua University , 31(6) (1991), pp. 1-11 (in Chinese)

1991
[27]

Liao, and Y

H.-L. Liao, and Y. Kang , L^2 norm error estimates of BDF methods up to fifth-order for the phase field crystal model , IMA J. Numer. Anal. , 44 (2024), pp. 2138-2164

2024
[28]

H.-L. Liao, T. Tang and T. Zhou , Discrete energy technique of the third-order variable-step BDF time-stepping for diffusion equations , J. Comput. Math. , 41 (2023), pp. 325-344

2023
[29]

H.-L. Liao, T. Tang and T. Zhou , A new discrete energy technique for multi-step backward difference formulas , CSIAM Trans. Appl. Math. , 3(2) (2022), pp. 318-334

2022
[30]

H.-L. Liao, T. Tang and T. Zhou , Positive definiteness of real quadratic forms resulting from variable-step L1-type approximations of convolution operators , Sci. China. Math. , 67:2 (2024), pp. 237-252

2024
[31]

Liao and Z

H.-L. Liao and Z. Zhang , Analysis of adaptive BDF2 scheme for diffusion equations , Math. Comput. , 90 (2021), pp. 1207-1226

2021
[32]

J. G. Liu, J.-G. Liu and R.L. Pego , Stability and convergence of efficient Navier-Stokes solvers via a commutator estimate , Commun. Pure Appl. Math. , 60:10 (2007), pp. 1443-1487

2007
[33]

Lubich, D

C. Lubich, D. Mansour, and C. Venkataraman , Backward difference time discretization of parabolic differential equations on evolving surfaces , IMA J. Numer. Anal. , 33 (2013), pp. 1365-1385

2013
[34]

Nevanlinna and F

O. Nevanlinna and F. Odeh , Multiplier techniques for linear multistep methods , Numer. Funct. Anal. Optim. 3 (1981), pp. 377-423

1981
[35]

Rosales, B

R. Rosales, B. Seibold, D. Shirokoff and D. Zhou , Unconditional stability for multistep ImEx schemes: theory , SIAM J. Numer. Anal. , 55 (2017), pp. 2336-2360

2017
[36]

Seibold, D

B. Seibold, D. Shirokoff and D. Zhou , Unconditional stability for multistep ImEx schemes: practice , J. Comput. Phys. , 376 (2019), pp. 295-321

2019
[37]

Wang and S

D. Wang and S. T. Ruuth , Variable step-size implicit-explicit linear multistep methods for time-dependent partial differential equations , J. Comput. Math. , 26(6) (2008), pp. 838-855

2008