Cubic-Quartic Equations Overview

• Cubic-quartic equations are a diverse family defined by the interplay between third- and fourth-order structures, covering classical polynomials, nonlinear PDEs, and optimization models.
• They employ geometric, algebraic, and iterative methods to analyze discriminants, root configurations, and convergence properties in distinct contexts.
• The framework is pivotal in bridging theory and application in areas like ultracold gases, nonlinear Schrödinger dynamics, and quartically regularized optimization.

“Cubic-quartic equation” does not denote a single universally fixed object. In the cited literature it refers, depending on context, to classical cubic and quartic polynomial equations, to coefficient-based classification problems for real roots, to solvability criteria by radicals or by iteration, to nonlinear Schrödinger and Gross–Pitaevskii equations containing cubic and quartic terms, and to optimization subproblems formed by a cubic model with quartic regularization (Prodanov, 2022, 2002.04001, Zhu et al., 2023). The unifying feature is the interaction between third- and fourth-order structure, but the mathematical meaning changes substantially across algebra, analysis, and numerical optimization.

1. Terminological scope and principal meanings

In classical algebra, the relevant objects are the monic cubic

$x^3+ax^2+bx+c$

and the monic quartic

$x^4+ax^3+bx^2+cx+d,$

together with their depressed forms, discriminants, resolvents, and root configurations (Prodanov, 2022). In this setting, “cubic-quartic” commonly means the joint study of cubic and quartic equations rather than a single mixed-degree scalar equation.

In nonlinear-wave and condensate theory, the phrase denotes evolution equations whose nonlinear response contains both cubic and quartic terms. A representative example is the trapped cubic–quartic nonlinear Schrödinger equation

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$

where the cubic term is the mean-field Gross–Pitaevskii contribution and the quartic term is a beyond-mean-field correction associated with Lee–Huang–Yang-type physics (2002.04001).

In optimization, the expression refers neither to a scalar cubic nor to a quartic polynomial equation in the classical sense. It refers to the optimality system for minimizing a cubic multivariate model with quartic regularization, and the computational core becomes a secular equation in an auxiliary parameter $\lambda$ (Zhu et al., 2023). A recurrent misconception is therefore that every “cubic-quartic equation” is a one-variable algebraic equation. The cited literature shows that this is false in both PDE and optimization contexts.

2. Algebraic equations, geometry, and coefficient structure

A geometric treatment of cubic and quartic equations associates real roots of a cubic with the Siebeck–Marden–Northshield equilateral triangle and real roots of a quartic with a regular tetrahedron (Prodanov, 2022). For the cubic

$p_3(x)=x^3+ax^2+bx+c,$

the centroid projects to $-a/3$, the critical points are

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$

and the condition $a^2-3b>0$ is the condition that the triangle exists. The paper gives a necessary and sufficient condition for three real roots: $a^2-3b>0$ together with

$$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$

Under this parametrization, variation of the free term rotates the triangle and yields Viète-type trigonometric formulas for the three real roots.

For the quartic

$x^4+ax^3+bx^2+cx+d,$0

the tetrahedral picture places the centroid at $x^4+ax^3+bx^2+cx+d,$1, and the inflection points are the roots of

$x^4+ax^3+bx^2+cx+d,$2

namely

$x^4+ax^3+bx^2+cx+d,$3

The quartic discriminant is cubic in $x^4+ax^3+bx^2+cx+d,$4, so varying the free term again becomes structurally important. The condition for four distinct real roots is expressed as

$x^4+ax^3+bx^2+cx+d,$5

together with a further condition placing $x^4+ax^3+bx^2+cx+d,$6 between the critical values $x^4+ax^3+bx^2+cx+d,$7, or between $x^4+ax^3+bx^2+cx+d,$8 and $x^4+ax^3+bx^2+cx+d,$9 when $$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$0 (Prodanov, 2022).

The same paper derives sharp interval-length and root-bound statements for quartics with four real roots. The maximum interval length containing the roots is

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$1

the minimum is

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$2

and no root can lie farther than

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$3

from $$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$4. These bounds are geometric consequences of the tetrahedral model.

A complementary coefficient-based program studies the quartic by the intersection of the “sub-quartic”

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$5

with the line

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$6

In this approach the first auxiliary cubic

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$7

gives the stationary points of the quartic, while the second auxiliary cubic

$$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$8

governs the $$i\frac{\partial\psi}{\partial t} = -\mathcal{A}\frac{\partial^2\psi}{\partial x^2} +\mathcal{B}|\psi|^2\psi +\mathcal{C}|\psi|^3\psi +\mathcal{D}(x)\psi,$$9 member of the family. Several subsidiary quadratic equations then provide exact threshold values, interval bounds, and isolation information for the real roots (Prodanov, 2020). This makes the quartic problem explicitly cubic-dependent: the derivative cubic and a second cubic organize the root geometry of the quartic.

3. Root configuration, multiplicity, and positivity regions

A Sturm-sequence treatment of real monic cubics and quartics gives complete coefficient criteria for real and complex root multiplicities and for the order of repeated real roots (Gonzalez et al., 2015). For the cubic

$\lambda$0

the discriminant

$\lambda$1

classifies the simple-root cases: $\lambda$2 When $\lambda$3, the condition $\lambda$4 gives one double real root and one simple real root, while $\lambda$5 gives a triple root. The double root is

$\lambda$6

and the sign of

$\lambda$7

determines whether the simple root lies to its left or right.

For the quartic

$\lambda$8

the discriminant $\lambda$9 is supplemented by

$p_3(x)=x^3+ax^2+bx+c,$0

and

$p_3(x)=x^3+ax^2+bx+c,$1

The classification includes: $p_3(x)=x^3+ax^2+bx+c,$2

$p_3(x)=x^3+ax^2+bx+c,$3

$p_3(x)=x^3+ax^2+bx+c,$4

$p_3(x)=x^3+ax^2+bx+c,$5

$p_3(x)=x^3+ax^2+bx+c,$6

$p_3(x)=x^3+ax^2+bx+c,$7

and

$p_3(x)=x^3+ax^2+bx+c,$8

A different but related line of work studies positivity instead of mere root count. For the cubic

$p_3(x)=x^3+ax^2+bx+c,$9

strict positivity on $-a/3$0 is equivalent to either

$-a/3$1

or

$-a/3$2

where

$-a/3$3

The boundary case $-a/3$4 with either $-a/3$5 or $-a/3$6 gives nonnegativity but not strict positivity, and the unique positive zero is

$-a/3$7

(Qi et al., 2020).

For the normalized quartic

$-a/3$8

nonnegativity on all real numbers is equivalent to

$-a/3$9

and either

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$0

or

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$1

Strict positivity allows either the same inequalities with $\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$2, or the discriminant-zero “appendix curve”

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$3

A notable consequence is that $\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$4 need not indicate a boundary point of the quartic nonnegativity region; on the appendix curve it lies inside the interior (Qi et al., 2020).

4. Solvability, radicals, and iterative cubic computation

Over $\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$5, 1-solvability is stricter than ordinary solvability by radicals. For the depressed cubic

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$6

an irreducible cubic is 1-solvable if and only if it has a rational root or

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$7

with

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$8

For the irreducible depressed quartic

$\mu_{1,2}=-\frac a3\pm \frac13\sqrt{a^2-3b},$9

1-solvability is equivalent to the existence of a rational root $a^2-3b>0$0 of the cubic resolvent

$a^2-3b>0$1

such that $a^2-3b>0$2 and

$a^2-3b>0$3

is a square of a rational number (Akhtyamov et al., 2014). This places the quartic criterion explicitly under the control of an auxiliary cubic.

A structurally different derivation of cubic radicals uses Harrison’s center theory. For

$a^2-3b>0$4

the associated binary cubic can be diagonalized into a sum of two cubes of linear forms, yielding a completed-cube decomposition and, in the depressed case, Cardano’s formula. The same paper extends the “completing powers” philosophy to special higher-degree equations, but explicitly states that even for quartics the method works only for very special equations (Huang et al., 2023).

The literature also contains iterative alternatives that avoid using Cardano as the computational mechanism. One geometric-iterative method starts from a monic cubic complex polynomial with distinct roots and distinct critical points, proves that at least one critical point has the Voronoi property, and then applies Kalantari’s basic family

$a^2-3b>0$5

to converge to a root. The quadratic formula enters only through the computation of the critical points of $a^2-3b>0$6, so the method is “solving a cubic by the quadratic formula” in an iterative rather than radical sense (Kalantari, 2014).

Another numerical program, inspired by Sharaf al-Din Tusi and Smale’s point-estimation theory, reduces any real cubic to one of four canonical forms and then chooses a certified Newton seed from

$a^2-3b>0$7

where $a^2-3b>0$8 approximates $a^2-3b>0$9 within five percent relative error. For suitable canonical forms, after $a^2-3b>0$0 Newton iterations the error satisfies

$a^2-3b>0$1

providing a high-precision route to a real root without Cardano’s formula (Kalantari, 2023).

5. Cubic-quartic nonlinear Schrödinger and Gross–Pitaevskii equations

In ultracold-gas theory, the cubic–quartic nonlinear Schrödinger equation models the competition between mean-field and beyond-mean-field interactions. In quasi-1D, the stationary reduction of

$a^2-3b>0$2

leads to

$a^2-3b>0$3

With a matched cnoidal trap,

$a^2-3b>0$4

the equation admits exact $a^2-3b>0$5- and $a^2-3b>0$6-families. The $a^2-3b>0$7-branch yields, in the limit $a^2-3b>0$8, the localized solution

$a^2-3b>0$9

while the $$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$0-branch yields the delocalized kink–antikink form

$$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$1

in the sense stated in the source formulas. A central conclusion is the absence of a nontrivial sinusoidal mode: in both branches the $$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$2 limit forces the oscillatory amplitude to zero (2002.04001).

A driven version in a bi-chromatic optical lattice adds a phase-locked source term and admits the exact periodic solution

$$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$3

with

$$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$4

Its density forms a stripe or density-wave pattern, and VK analysis gives stability for $$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$5. The same work extends the model to quadratic-cubic-quartic and quadratic-cubic nonlinearities in order to describe dimensional crossover between quasi-1D and strict 1D (Debnath et al., 2021).

A nonlocal cubic–quartic Gross–Pitaevskii equation in $$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$6,

$$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$7

is treated variationally on the mass shell $$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$8. Because the constrained energy is unbounded below, ground states are constructed as mountain-pass critical points on the $$- \frac{2a^3}{27}+\frac{ab}{3}-\frac{2}{27}\sqrt{(a^2-3b)^3} \le c \le - \frac{2a^3}{27}+\frac{ab}{3}+\frac{2}{27}\sqrt{(a^2-3b)^3}.$$9-sphere. The analysis uses localized Palais–Smale sequences and the nonlocal Pohozaev identity

$x^4+ax^3+bx^2+cx+d,$00

which has no remainder term for the dipolar kernel. For suitable masses, ground states exist, are positive up to a constant phase, and have positive chemical potential (Luo et al., 2018).

Another integrable usage occurs in an extended nonlinear Schrödinger equation with third- and fourth-order corrections,

$x^4+ax^3+bx^2+cx+d,$01

Although the title speaks of cubic and quartic nonlinearity, the full model also contains derivative nonlinear couplings and the quintic term $x^4+ax^3+bx^2+cx+d,$02. Using the generalized Darboux transformation, the paper constructs higher-order smooth positons and breather-positons and derives the exact time-dependent displacement

$x^4+ax^3+bx^2+cx+d,$03

showing direct dependence on the higher-order parameters $x^4+ax^3+bx^2+cx+d,$04 and $x^4+ax^3+bx^2+cx+d,$05 (Monisha et al., 2022). This is another context in which “cubic-quartic” is descriptive but not exhaustive.

6. Cubic-quartic regularization in higher-order optimization

In third-order tensor methods for optimization, the relevant model is the quartically regularized cubic subproblem

$x^4+ax^3+bx^2+cx+d,$06

The proposed CQR framework replaces the exact cubic tensor term by a simpler local model

$x^4+ax^3+bx^2+cx+d,$07

where the cubic term approximates local tensor information and the quartic term regularizes progress (Zhu et al., 2023).

The global minimizer $x^4+ax^3+bx^2+cx+d,$08 is characterized by

$x^4+ax^3+bx^2+cx+d,$09

with

$x^4+ax^3+bx^2+cx+d,$10

Introducing

$x^4+ax^3+bx^2+cx+d,$11

reduces the optimality system to

$x^4+ax^3+bx^2+cx+d,$12

together with the scalar secular equation

$x^4+ax^3+bx^2+cx+d,$13

Accordingly, the “cubic-quartic equation” here is a coupled optimality system plus scalar root-finding reduction, not a classical scalar quartic equation.

The paper proves necessary and sufficient optimality conditions for global minimizers and shows that a CQR variant achieves the optimal-order evaluation complexity

$x^4+ax^3+bx^2+cx+d,$14

for minimizing the quartically regularized cubic subproblem, with improvements in special cases. In the univariate case, the quartically regularized cubic can be solved within a single CQR iteration by minimizing two scalar cubic-quartic models, and numerically the method is reported to be competitive with ARC and QQR, with particularly strong performance on ill-conditioned and badly scaled instances (Zhu et al., 2023).

Across these literatures, the cubic–quartic theme is structurally consistent even when the objects differ. In algebra it organizes discriminants, resolvents, and geometric root models; in analysis it encodes the competition of mean-field and higher-order effects; and in optimization it links higher-order local models to secular equations and regularization. The term is therefore best understood not as a single definition, but as a family of mathematically precise settings in which cubic and quartic mechanisms interact.