cert-machine · report · a new kernel, re-derived at every build
The regularization, measured
Ferreira, Gomes and Üçer prove that stationary first-order mean-field games have solutions by casting the system as a monotone operator on a Banach space and adding a small p-Laplacian to make it coercive; the regularization is then sent to zero. The proof needs the ε. Does the solution? On the one-dimensional game with a quadratic Hamiltonian and a cosine potential this page encloses every regularized equilibrium and the unregularized one, cell by cell across an atlas in the potential's amplitude and the regularization's strength, and decides the distance between them. The unregularized game certifies as readily as any regularized one and further than the strongly regularized ones, because after the density is eliminated the game is a second-order elliptic equation for the value function wherever the density is positive. The regularization was for the proof.
tl;dr
The finding. 34 of 50 cells PROVED, each a unique classical even solution with density bounded away from zero; 16 refused by name. Along the amplitude every row is proved then refused: ε = 0 reaches A = 0.7, ε = 1 only A = 0.4. Where both are proved, ‖u_ε − u_0‖_2 is decided to the two radii: 1.043 ε at ε = 0.01 and A = 0.4, sublinear beyond.
The mechanism. Eliminating m = u + u′²/2 − V turns the transport equation into F_ε(u) = m − (m u′ + ε u′³)′ + ε u³ − 1 = 0, whose linearization is −(c h′)′ + e h with c = m + (1 + 3ε) u′² ≥ m. A radii polynomial in the derivative-weighted ℓ¹ space X_2 (explicit columns to six times the truncation, an analytic tail beyond, the algebra property for the nonlinear term) closes a ball about a Galerkin candidate.
Check it.node instruments/regatlas/battery.js (17 checks, 5 red controls, about 25 s) re-derives the atlas, holds the constant solution at A = 0 against the cubic it must solve, and an independent finite-difference solve against three certified cells.
cells proved
34 of 50
radii from 9.1e-16 to 6.8e-8 in X_2; every refusal is Z1 ≥ 1 or Z1 so near 1 that no radius closes
ε = 0 reaches
A = 0.7
as far as ε = 0.01 and further than ε = 1 (A = 0.4): the certifier needs no regularization; the strong regularization costs it
the density at the frontier
0.320
the certified floor of m at the last proved cell of ε = 0; the refusal beyond is the certificate's (a diagonal tail against a varying coefficient), not the density's
distance at ε = 0.01
0.010
at A = 0.4, decided to 6.1e-14: the regularized solution sits ε away, to first order
at ε = 1
0.358
a third of ε: sublinear, as the cubic ε u³ + u = 1 already shows at A = 0 (u = 0.682)
a new kernel
X_2
the first radii polynomial in this machine for a quasilinear, variable-coefficient problem; its bounds are stated in the code and their tail term is the frontier
§1 · the game and its regularization
Elliptic after all
The paper's Problem 1 on the torus, with the power-growth Hamiltonian H = |p|²/2 − m and a discount, reads
−u − u′²/2 + m + V = 0, m − (m u′)′ − 1 = 0, V = A cos 2πx.
Its regularized operator (3.1) adds, in the transport slot, ε times the γ̄-Laplacian of u and the γ̄-power of u with γ̄ = α(β + 1)/β; for α = 2, β = 1 that is γ̄ = 4, so the added terms are ε(u′³)′ and ε u³ — polynomial. The first equation gives m outright, and the second becomes one scalar equation for u:
F_ε(u) = m − (m u′ + ε u′³)′ + ε u³ − 1 = 0, m = u + u′²/2 − V.
Its linearization at ū is a Sturm–Liouville operator, −(c h′)′ + e h, with c = m̄ + (1 + 3ε)ū′² and e = 1 − ū″ + 3εū². The coefficient c is at least the density. Wherever m > 0 the operator is elliptic — at ε = 0 as at ε > 0 — and that is the whole reason the page can certify the unregularized game with the same instrument: the first-order game, after eliminating its density, is a second-order elliptic equation in disguise. The paper's ε is what makes the abstract operator coercive on the Banach space where existence is proved; it is not what makes the equation solvable.
§2 · the atlas
Fifty cells, each proved or refused by name
For every amplitude A from 0 to 0.9 and every ε in {0, 0.01, 0.1, 0.3, 1}, a Galerkin candidate with N = 24 to 40 cosine modes is found by Newton continuation from the constant solution, and a radii polynomial in the space X_2 = {u : Σ|u_k|(1 + k)²ν^k < ∞}, ν = 1.05, is asked to close a ball about it. The approximate inverse is the dense inverse of the Galerkin Jacobian on the first N modes and the diagonal 1/(c₀(2πk)² + e₀) beyond; Z1 is computed column by column to six times N and bounded analytically past that; Z2 comes from the algebra property of the weighted norm. A zero of F in X_2 is twice differentiable with an analytic Fourier series — a classical solution — and the density is certified positive on the whole ball, so it is a strong solution in the paper's sense with m bounded away from zero.
Figure 1 · The atlas. Solid: PROVED, the shade darkening as the contraction factor κ falls; hatched: REFUSED, with the Z1 that refused it. Every row is proved then refused as A grows, and the proved region shrinks as ε grows.
Every refusal is the same refusal: Z1, the norm of I − A·DF(ū), reaches 1. It is the tail's doing. A diagonal approximate inverse beyond N is right when the principal coefficient c is nearly constant and wrong in proportion to its variation, and c = m + (1 + 3ε)u′² varies more as A grows (m follows −V) and as ε grows (the 4-Laplacian weights u′² three times more). So the frontier moves left with ε: the regularization that gives the paper its coercivity takes contraction from the certificate. The density is nowhere near zero at the frontier — the certified floor at the last proved cell of ε = 0 is 0.320, and the float candidate at the first refused cell still has min m = 0.225. The void of this atlas is the certificate's boundary, not the equation's; the two would only coincide where m reaches zero, which is where the paper's strong solutions stop being classical, and that is beyond the map.
Figure 2 · The contraction factor along the amplitude for three regularization strengths, on the proved cells. Each curve rises toward 1 and stops at its frontier; the strongly regularized one first.
§3 · the regularization measured
ε away, to first order
Where a regularized cell and its unregularized neighbour are both proved, the distance ‖u_ε − u_0‖_2 is decided: the norm of the difference of the two candidates, exactly, plus or minus the two radii. At A = 0 the solutions are constants — u_ε is the real root of εu³ + u = 1 — and the distance is 1 − u_ε, which the battery holds against the cubic. At every amplitude the distance grows with ε, equals about ε at ε = 0.01, and is a third of ε at ε = 1: the regularized game is a first-order perturbation of the unregularized one in this norm, and the perturbation saturates because the cubic term is bounded. The paper passes ε → 0 by compactness; here the limit is watched.
Figure 3 · The decided distance ‖u_ε − u_0‖_2 against ε for three amplitudes, on logarithmic axes; the brackets are narrower than the stroke. The dashed guide is the diagonal.
A
ε = 0.01
ε = 0.1
ε = 0.3
ε = 1
0
[9.7115e-3, 9.7115e-3]
[7.8301e-2, 7.8301e-2]
[1.7095e-1, 1.7095e-1]
[3.1767e-1, 3.1767e-1]
0.1
[9.8152e-3, 9.8152e-3]
[7.9193e-2, 7.9193e-2]
[1.7308e-1, 1.7308e-1]
[3.2231e-1, 3.2231e-1]
0.2
[9.9527e-3, 9.9527e-3]
[8.0393e-2, 8.0393e-2]
[1.7602e-1, 1.7602e-1]
[3.2900e-1, 3.2900e-1]
0.3
[1.0146e-2, 1.0146e-2]
[8.2107e-2, 8.2107e-2]
[1.8034e-1, 1.8034e-1]
[3.3950e-1, 3.3950e-1]
0.4
[1.0434e-2, 1.0434e-2]
[8.4723e-2, 8.4723e-2]
[1.8719e-1, 1.8719e-1]
[3.5779e-1, 3.5779e-1]
0.5
[1.0896e-2, 1.0896e-2]
[8.9049e-2, 8.9049e-2]
[1.9914e-1, 1.9914e-1]
UNDECIDED
0.6
[1.1713e-2, 1.1713e-2]
[9.7017e-2, 9.7017e-2]
UNDECIDED
UNDECIDED
0.7
[1.3368e-2, 1.3368e-2]
UNDECIDED
UNDECIDED
UNDECIDED
0.8
UNDECIDED
UNDECIDED
UNDECIDED
UNDECIDED
0.9
UNDECIDED
UNDECIDED
UNDECIDED
UNDECIDED
Figure 4 · The density at A = 0.6 for the proved regularization strengths. All three are certified to radii below 1e−9; the regularized densities are slightly flatter and, through the cubic term, carry slightly less mass where u is large.
Figure 5 · The value function at A = 0.6, the same cells. The regularized ones sit lower by about ε: the zero-order term εu³ acts like an added discount.
§4 · the honest boundary
What is claimed, and what is not
Proved, cell by cell. A unique zero of F_ε in a ball of X_2 about the candidate, in the even subspace; classical; density positive on the ball. Each certificate carries five falsifiers, and the battery fires them.
Decided. The distances between doubly-proved cells; the ellipticity of the linearization (c₀ > 0, e₀ ≥ 0) on every proved cell; the constant solution's cubic at A = 0.
Not claimed. Uniqueness outside the ball or outside the even subspace; anything on a refused cell — existence there is neither asserted nor denied; the paper's d-dimensional theorems; the singular-congestion and weak-growth cases (Theorems 1.5, 1.6). The one-dimensional quadratic game with a single cosine mode is the instance, and it is the instance because its regularization is polynomial.
The frontier is the certificate's. Every refusal is Z1 ≥ 1 from a diagonal tail against a varying coefficient. A banded tail inverse would move it; the density would not. The page says so rather than drawing a boundary it does not have.
A new kernel. The radii-polynomial bounds for this quasilinear problem are written out in kernel.js and are this machine's, not lifted. They are checked here by a finite-difference solve on a different discretization and by the constant solution; they have not been reviewed by anyone else.
§5 · check it
Half a minute on your machine
node instruments/regatlas/battery.js # 17 checks, 5 red controls: the atlas re-derived, the cubic at A = 0, finite differences, the frontier
node instruments/regatlas/run.js --check # the record, re-derived and compared byte for byte
The kernel is instruments/regatlas/kernel.js, with derive.js; the record is certs/regatlas.json.
references
Sources
R. Ferreira, D. A. Gomes, M. Üçer, Solving mean-field games with monotonicity methods in Banach spaces, arXiv:2506.21212v3 (2026) — Problem 1, Assumption 2.4 (power growth), the regularized operator (3.1), Proposition 3.1 (coercivity), Theorem 1.4.
J. B. van den Berg, J.-P. Lessard, Rigorous numerics in dynamics, Notices AMS 62 (2015); S. Day, J.-P. Lessard, K. Mischaikow, SIAM J. Numer. Anal. 45 (2007) — the radii polynomial; the derivative-weighted ℓ¹ space is the standard device for quasilinear problems.