---
title: 'AdaCubic: Adaptive Cubic Optimizer'
url: https://www.emergentmind.com/topics/adacubic
type: topic
---

# AdaCubic: Adaptive Cubic Optimizer

Searching arXiv for AdaCubic and related optimizer/context papers.
arXiv search query: AdaCubic optimizer cubic regularization deep learning
AdaCubic is a novel regularization technique that adapts the weight of the cubic term and an adaptive, cubic-regularization–based second-order optimizer specifically engineered to be practical for deep learning [2604.09437]. It starts from cubic-regularized Newton, replaces the fixed cubic weight by an automatically learned one, and uses a Hutchinson-based diagonal Hessian approximation to reduce computational cost and memory. The method is positioned for non-convex empirical risk minimization in large neural networks, with reported evaluations in Computer Vision, Natural Language Processing, Language Modeling, and Signal Processing. A broader, noncanonical use of the label is also suggested by explanatory treatments of cubic equations and cubic Diophantine problems, where “AdaCubic” denotes an adaptive or geometry-aware cubic procedure rather than the 2026 optimizer [1401.5148, 2105.03352].

## 1. Historical placement and naming

In its primary arXiv usage, AdaCubic denotes a deep-learning optimizer built around cubic regularization [2604.09437]. Its conceptual starting point is the classical cubic-regularized Newton model of Nesterov & Polyak,
\[
m_M(s) = f(x_k) + \nabla f(x_k)^\top s + \tfrac12 s^\top H_k s + \tfrac{M}{6} \|s\|^3,
\]
where the cubic term regularizes the local second-order model and supports convergence to an $(\epsilon_g,\epsilon_H)$-stationary point with iteration complexity \(O(\epsilon^{-3/2})\). The motivation for AdaCubic is that direct deployment of cubic-regularized Newton in deep learning is obstructed by the prohibitive cost of computing and storing \(H_k\), and by the need to choose a robust value of \(M\).

AdaCubic addresses both issues simultaneously. It uses Hutchinson’s method to approximate only the diagonal of the Hessian from Hessian–vector products, giving an \(O(d)\)-memory second-order method, and it introduces an auxiliary constrained problem that turns the cubic term into a constraint and recovers the cubic weight \(M\) as a dual variable \(\nu_k\). The optimizer is therefore not a fixed-\(M\) cubic method, but a trust-region–style cubic method in which the regularization weight is adapted per iteration [2604.09437].

A recurring source of terminological ambiguity is that “AdaCubic” has also been used informally in explanatory discussions outside optimization. Kalantari’s work on solving cubic equations by the quadratic formula uses critical points, Voronoi geometry, and the basic family \(B_m\) to obtain a root of a cubic from quadratic-level computations [1401.5148]. A separate line of work on \(x^3 + D = 2^c y\) organizes solutions through 2-adic valuation trees and has been described from an “AdaCubic” perspective because of its adaptive tree traversal and congruence-structured solution generation [2105.03352]. These latter usages are interpretive rather than canonical titles.

## 2. Core mathematical formulation

AdaCubic targets minimization of a non-convex empirical risk
\[
f(x) = \frac{1}{n}\sum_{\ell=1}^n f_\ell(x), \quad x \in \mathbb{R}^d.
\]
Second-order stationarity is measured by
\[
\|\nabla f(x)\|_2 \le \epsilon_g, \qquad \lambda_{\min}(\nabla^2 f(x)) \ge -\epsilon_H.
\]

The central reformulation replaces the unconstrained cubic model by the constrained problem
\[
\begin{aligned}
\min_{s\in\mathbb{R}^d}\quad & \hat m(s) := f(x_k) + \nabla f(x_k)^\top s + \tfrac12 s^\top H_k s \\
\text{s.t.}\quad & \tfrac16(\|s\|_2^3 - \xi)\le 0,
\end{aligned}
\]
with Lagrangian
\[
\mathcal{L}_\xi(s,\nu)=f(x_k)+\nabla f(x_k)^\top s+\tfrac12 s^\top H_k s+\tfrac{\nu}{6}(\|s\|_2^3-\xi), \qquad \nu\ge 0.
\]
The KKT system yields three conditions: a first-order equation,
\[
\Bigl(H_k+\frac{\nu^*}{2}\|s^*\|_2 I\Bigr)s^*=-\nabla f(x_k),
\]
a positive-semidefinite condition,
\[
H_k+\frac{\nu^*}{2}\|s^*\|_2 I \succeq 0,
\]
and complementary slackness,
\[
\nu^*(\|s^*\|_2^3-\xi)=0.
\]
The dual variable \(\nu^*\) therefore plays exactly the role of the cubic regularization weight [2604.09437].

The dual formulation introduces
\[
s(\nu,r)= -\bigl(H_k + \tfrac{\nu r}{2} I\bigr)^{-1}\nabla f(x_k),
\]
and a scalar function \(\mathscr{L}_\xi(\nu,r)\). The maximizing pair \((\nu^*,r^*)\) satisfies \(r^*=\sqrt[3]{\xi}\), and \(\nu^*\) solves the scalar equation
\[
\phi(\nu^*,r^*) := \frac{1}{\|s(\nu^*,r^*)\|_2} - \frac{1}{\sqrt[3]{\xi}} = 0.
\]
The paper’s equivalence theorem states that the constrained formulation and the cubic-regularized step coincide when the cubic weight \(M\) is set to the optimal Lagrange multiplier \(\nu^*\). In this sense, AdaCubic replaces heuristic tuning of \(M\) by solving for an iteration-dependent dual optimum [2604.09437].

## 3. Algorithmic mechanism and trust-region updates

In practice, AdaCubic uses a batch-based approximate model,
\[
\mathfrak{m}_M(s) = F(x_k) + g_k^\top s + \tfrac12 s^\top B_k s + \tfrac{M}{6}\|s\|^3,
\]
where \(F(x_k)\) and \(g_k\) are minibatch estimates and \(B_k\) is a diagonal Hessian approximation. The inner solver, denoted RootFinder, computes \((s_{k+1},\nu_{k+1})\) as the minimizer of the constrained problem in this approximate model.

The outer iteration follows a trust-region logic in the cubic norm. At each step, the method computes a trial step, evaluates the ratio
\[
\rho_k =
\frac{F(x_k) - F(x_k + s_{k+1})}
     {F(x_k) - \mathfrak{m}_{\nu_{k+1}}(s_{k+1})},
\]
and accepts the step when \(\rho_k \ge \eta_1\). The quantity \(\xi_k\) acts as a trust-region radius parameter through the constraint \(\|s\|^3 \le \xi_k\). The update policy is:

- if \(\rho_k \ge \eta_2\), the iteration is very successful and
  \[
  \xi_{k+1} = \max\{\alpha_1 \|s_{k+1}\|^3,\,\xi_k\};
  \]
- if \(\eta_1 \le \rho_k < \eta_2\), the iteration is successful and \(\xi_{k+1}=\xi_k\);
- if \(\rho_k < \eta_1\), the iteration is unsuccessful and
  \[
  \xi_{k+1} = \max\{\alpha_2 \|s_{k+1}\|^3,\;\epsilon_m\}.
  \]

The inner RootFinder sets \(r=\sqrt[3]{\xi}\), forms
\[
\tilde H(\nu,r)=B_k+\tfrac12 \nu r I,
\]
and solves for \(\nu\) so that the corresponding step
\[
s(\nu,r)= -\tilde H(\nu,r)^{-1} g_k
\]
satisfies \(\|s(\nu,r)\|_2 \approx r\). Because \(B_k + \tfrac{\nu r}{2}I\) is diagonal, the linear algebra reduces to cheap elementwise operations. The one-dimensional equation \(\phi(\nu,r)=0\) is handled by safeguarded Newton iterations, exploiting the monotonicity and concavity properties reported for \(\phi\) [2604.09437].

The hyperparameters are taken directly from trust-region literature and fixed across all experiments:
\[
\eta_1 = 0.05,\quad \eta_2 = 0.75,\quad \alpha_1 = 2.5,\quad \alpha_2 = 0.25,\quad \kappa_{\text{easy}} = 0.01,\quad \epsilon_m = 10^{-6}.
\]
This fixed configuration is central to the paper’s positioning of AdaCubic as an optimizer for settings where fine-tuning is infeasible [2604.09437].

## 4. Hessian approximation, complexity, and convergence theory

A defining scalability feature is that AdaCubic estimates only the diagonal of the Hessian. If \(\mathcal{H}_k(v)=H_kv\) denotes a Hessian–vector product, then Hutchinson’s method forms
\[
b_k = \frac{1}{\mathcal S}\sum_{s=1}^{\mathcal S}\mathcal H_k(v_s)\odot v_s \in \mathbb{R}^d,
\]
using Rademacher probe vectors \(v_s\), and sets
\[
B_k=\mathrm{Diag}(b_k)\in\mathbb{R}^{d\times d}.
\]
The estimator is unbiased for \(\mathrm{diag}(H_k)\), requires one Hessian–vector product per probe, and yields \(O(d)\) memory rather than the storage cost of a dense \(d\times d\) Hessian [2604.09437].

The deterministic theory assumes that each \(f_i(x)\) is twice continuously differentiable and bounded below, and that \(f_i\), \(\nabla f_i\), and \(\nabla^2 f_i\) are Lipschitz continuous on a closed convex set \(\mathcal F\). With the exact diagonal Hessian \(D(x)=\mathrm{Diag}(\nabla^2 f(x))\), the resulting diagonal cubic-regularized iteration satisfies
\[
\min_{1\le i\le k}\|\nabla f(x_i)\|_2 \le \mathcal O(k^{-2/3}),
\]
and therefore reaches \(\min_i\|\nabla f(x_i)\|\le \epsilon\) in
\[
k \ge \mathcal O(\epsilon^{-3/2})
\]
iterations. The paper presents this as the classical cubic-regularization iteration complexity inherited in the diagonal approximation setting [2604.09437].

For the stochastic case, the analysis imposes inexact-model bounds
\[
\|g_k-\nabla f(x_k)\|\le C_g\|s_k\|^2,
\qquad
\|B_k-\mathrm{Diag}(\nabla^2 f(x_k))\|\le C_B\|s_k\|.
\]
Probabilistic deviation bounds are then derived with vector and matrix Bernstein inequalities. The gradient bound holds with probability \(1-\delta\) provided
\[
|\mathcal B_k^g|
\ge
32L_f^2\frac{\ln \tfrac1\delta + 1/4}{C_g^2\|s_k\|^4},
\]
and the Hessian-diagonal bound holds with probability \(1-\delta\) provided
\[
|\mathcal B_k^H| \ge \sqrt{d}L_g
\frac{\ln \tfrac{2d}{\delta}}{\mathcal S C_B\|s_k\|}.
\]
The stated implication is asymptotic: as \(\|s_k\|\to 0\), batch sizes must become large, approaching full-data behavior, and the deterministic cubic-regularization guarantees apply locally [2604.09437].

## 5. Empirical performance across application domains

The experimental program covers Computer Vision, GLUE natural language understanding, Language Modeling, and audio-based Camera Model Identification [2604.09437]. In all experiments AdaCubic uses the same universal hyperparameters, whereas the baselines are tuned.

On CIFAR-10, with ResNet20 and ResNet32, AdaCubic clearly outperforms SGD and Adam and remains close to AdaHessian. For ResNet20, the reported accuracies are AdaHessian \(91.64\% \pm 0.46\), AdaCubic \(91.49\% \pm 0.46\), Adam \(90.26\%\), and SGD \(88.52\%\). For ResNet32, the reported values are AdaHessian \(93.15\% \pm 0.12\), AdaCubic \(92.65\% \pm 0.19\), Adam \(91.24\%\), and SGD \(89.02\%\). On CIFAR-100 with ResNet18, AdaCubic is slightly behind tuned baselines: without spatial averaging, SGD \(72.62\%\), Adam \(72.04\%\), AdaHessian \(72.45\%\), and AdaCubic \(71.81\%\); with spatial averaging, AdaHessian \(72.59\% \pm 0.27\) and AdaCubic \(72.00\% \pm 0.34\).

On GLUE with SqueezeBERT, AdaCubic is reported as either best or second-best across tasks despite the absence of per-task tuning. Across SST-2, QNLI, RTE, and WNLI, average accuracy is SGD \(77.17\%\), AdaCubic \(76.86\%\), and AdaHessian \(75.37\%\). On MRPC, AdaCubic reaches F1 \(0.9042\) and accuracy \(86.76\%\), compared with SGD F1 \(0.9094\), accuracy \(87.25\%\), and AdaHessian F1 \(0.8562\), accuracy \(78.18\%\). On QQP, AdaCubic reports F1 \(0.8723\) and accuracy \(90.40\%\), against SGD F1 \(0.8775\), accuracy \(90.89\%\), and AdaHessian F1 \(0.8742\), accuracy \(90.82\%\). On MNLI, AdaCubic reports \(81.88/81.89\) for matched/mismatched, compared with SGD \(82.45/82.05\) and AdaHessian \(81.65/81.57\).

On language modeling, the optimizer is often much better than AdaHessian and sometimes better than SGD. On WikiText-2, perplexity is reported as \(3.756\) for AdaCubic and \(4.374\) for AdaHessian on RoBERTa, \(5.759\) for AdaCubic and \(16.151\) for AdaHessian on BERT, and \(6.565\) for AdaCubic and \(6.822\) for AdaHessian on DistilBERT. On PTB, the reported perplexities are \(5.145\) for AdaCubic and \(7.582\) for AdaHessian on RoBERTa, \(14.170\) for AdaCubic and \(20.851\) for AdaHessian on BERT, and \(7.334\) for AdaCubic and \(10.182\) for AdaHessian on DistilBERT.

In audio-based Camera Model Identification on the VISION dataset with ResNet18, AdaCubic is consistently more accurate and more stable than Adam. The reported mean accuracies over 5 folds are \(94.77\% \pm 1.43\) versus \(93.99\% \pm 1.76\) on native videos, \(93.68\% \pm 2.59\) versus \(91.11\% \pm 2.66\) on WhatsApp, and \(93.51\% \pm 1.01\) versus \(91.89\% \pm 4.98\) on YouTube.

Wall-clock comparisons indicate the expected trade-off between iteration cost and training progress. On ResNet20/CIFAR-10, to reach target loss \(0.15\), the reported figures are 83 epochs and 35.16 minutes for SGD, 81 epochs and 61.85 minutes for AdaHessian, and 55 epochs and 42.40 minutes for AdaCubic [2604.09437].

## 6. Limitations, misconceptions, and related usages beyond deep learning

AdaCubic is a second-order method, but it is not a full-Hessian cubic optimizer. Its curvature model is diagonal, so off-diagonal structure is ignored. The paper identifies this as a plausible explanation for weaker performance on CIFAR-100 and for small gaps versus tuned SGD on some GLUE tasks [2604.09437]. The method also incurs additional cost: Hutchinson estimation requires an extra backward pass, or a few, per iteration. Its practical memory overhead is described as around \(2d\) versus SGD’s \(d\), and below AdaHessian’s \(4d\). Theoretical local convergence in the stochastic setting is asymptotic in the sense that batch sizes must become large as \(\|s_k\|\to 0\).

A further misconception is terminological. Outside the 2026 optimizer, the label has been used loosely for cubic procedures whose adaptivity is geometric or arithmetic rather than optimization-theoretic. In Kalantari’s “Solving Cubic Equations By the Quadratic Formula,” a monic cubic with distinct roots and distinct critical points has at least one critical point with the Voronoi property, meaning that the critical point lies in the Voronoi cell of some root [1401.5148]. For any \(\xi\in V(\theta)\), the basic sequence
\[
B_m(\xi)=\xi-p(\xi)\frac{d_{m-2}}{d_{m-1}}
\]
converges to \(\theta\), where
\[
d_m=p'(\xi)d_{m-1}-0.5\,p(\xi)p''(\xi)d_{m-2}+p(\xi)^2d_{m-3},
\qquad d_0=1,\quad d_{-1}=d_{-2}=0.
\]
Because at least one critical point lies in a root’s Voronoi cell, solving \(p'(z)=0\) by the quadratic formula supplies a guaranteed-convergence starting point for a root of \(p(z)\). The rate of convergence is proportional to the ratio between the distance from the starting point to the nearest root and the distance to the second-closest root. This suggests an adaptive cubic procedure in a geometric sense, but it is numerically and conceptually distinct from the deep-learning optimizer.

An arithmetic analogue appears in work on \(x^3 + D = 2^c y\) with \(y\) odd, where 2-adic valuation trees classify the possible values of \(v_2(x^3+D)\) and the congruence classes of \(x \bmod 2^\ell\) that realize them [2105.03352]. The central cubic theorem states that if \(D = 2^{3k}(2j+1)\), then the valuation tree of \(x^3 + D\) is infinite and \(x^3 + D = 2^c y\) has nontrivial solutions for all \(c\ge 0\) except finitely many values; if \(D \neq 2^{3k}(2j+1)\), only finitely many exponents \(c\) occur. For the example \(D=1\),
\[
\nu_2(n^3+1)=c \quad\Longleftrightarrow\quad n\equiv 2^c-1 \pmod{2^{c+1}}.
\]
This line of work has likewise been described from an “AdaCubic” perspective because the solution space is traversed adaptively through a valuation tree, but it is a 2-adic Diophantine framework rather than an optimizer.

Taken together, these usages identify two distinct meanings. In the strict arXiv sense, AdaCubic is the adaptive cubic-regularization optimizer introduced in 2026 for scalable deep learning [2604.09437]. In a broader interpretive sense, the name also evokes cubic procedures whose adaptivity is governed by Voronoi geometry or valuation-tree structure rather than by a dual trust-region mechanism [1401.5148, 2105.03352].

Source: https://www.emergentmind.com/topics/adacubic