Order-eight splitting methods for quartic potentials from a certified computer search, and why a smaller ℓ1 norm is not enough

We revisit time-symmetric splitting methods of order eight for y'' = g(y), g = -grad V, when V is a polynomial of degree at most four. For this class Escorihuela-Tomàs recently constructed the 19-stage method QA19, whose coefficient norm is l1 = 3.8225. A computer search by a team of AI agents produced twenty certificates, each accepted by two independently written interval-arithmetic checkers. They fall into twelve groups of numerically close 19-stage methods of order eight, all with l1 below that of QA19; the smallest is l1 = 3.0513898849, 20.2% below QA19. High-precision polishing gives numerical evidence that the four best groups contain strict local minima of l1. An independent numerical study shows, however, that a small l1 does not by itself make a method efficient. On three quartic test problems, the two methods with the smallest l1 use 7-20% more force evaluations than QA19 to reach a given energy accuracy. A third method, with l1 = 3.1264 and the smallest leading-error norm of the methods tested, uses 25-27% fewer force evaluations than QA19 on an FPU-beta chain and on a quartic oscillator, but 1-2% more on a phi^4 chain. We also show that, for quartic potentials in dimension d >= 2, the degree-nine component of the relevant Lie algebra has dimension 17, not 18 as tabulated, because modified potentials commute. Files: d1_note.pdf (the paper); d1_deposit.zip (the twenty certificates, the frozen checkers K1 and K2, a one-command re-check script reverify_d1.py, the catalogue of methods with 80-digit coefficients, the numerically polished local-minimum candidates, the efficiency study and the dimension computation, with SHA256SUMS). Use of AI and the author's role: the search was carried out by a multi-agent AI research system (Claude Opus 5.5, Anthropic) operated by the author. Claude also wrote checker K1 and the search and polishing code, re-proved the certificates and drafted the text. GPT-6 Astra (OpenAI) wrote checker K2 blind, carried out the efficiency study and the dimension computation, audited the literature, and fact-checked the manuscript. The author chose the project and its criteria, required cross-family review, chose the stopping rules and decided when to end the search, and re-ran both checkers on all certificates on the author's own computer (every check accepted). The written arguments were drafted with AI and have not yet been checked line by line by the author or by a mathematician. Anyone can check every certified method with one command (python reverify_d1.py). No AI system is an author; the author takes full responsibility for the content. Details are in the section "Use of artificial intelligence" of the paper.

Authors

Institutions

Publication Details

Journal
Zenodo (CERN European Organization for Nuclear Research)
Published
2026-11-04
DOI
https://doi.org/10.5281/zenodo.23934734
Primary Topic
Numerical methods for differential equations
Type
preprint
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
preprint

Order-eight splitting methods for quartic potentials from a certified computer search, and why a smaller ℓ1 norm is not enough

Aoran Zhang
Zenodo (CERN European Organization for Nuclear Research)
Numerical methods for differential equations
preprint

Order-eight splitting methods for quartic potentials from a certified computer search, and why a smaller ℓ1 norm is not enough

Aoran Zhang
preprint en

Abstract

We revisit time-symmetric splitting methods of order eight for y'' = g(y), g = -grad V, when V is a polynomial of degree at most four. For this class Escorihuela-Tomàs recently constructed the 19-stage method QA19, whose coefficient norm is l1 = 3.8225. A computer search by a team of AI agents produced twenty certificates, each accepted by two independently written interval-arithmetic checkers. They fall into twelve groups of numerically close 19-stage methods of order eight, all with l1 below that of QA19; the smallest is l1 = 3.0513898849, 20.2% below QA19. High-precision polishing gives numerical evidence that the four best groups contain strict local minima of l1. An independent numerical study shows, however, that a small l1 does not by itself make a method efficient. On three quartic test problems, the two methods with the smallest l1 use 7-20% more force evaluations than QA19 to reach a given energy accuracy. A third method, with l1 = 3.1264 and the smallest leading-error norm of the methods tested, uses 25-27% fewer force evaluations than QA19 on an FPU-beta chain and on a quartic oscillator, but 1-2% more on a phi^4 chain. We also show that, for quartic potentials in dimension d >= 2, the degree-nine component of the relevant Lie algebra has dimension 17, not 18 as tabulated, because modified potentials commute. Files: d1_note.pdf (the paper); d1_deposit.zip (the twenty certificates, the frozen checkers K1 and K2, a one-command re-check script reverify_d1.py, the catalogue of methods with 80-digit coefficients, the numerically polished local-minimum candidates, the efficiency study and the dimension computation, with SHA256SUMS). Use of AI and the author's role: the search was carried out by a multi-agent AI research system (Claude Opus 5.5, Anthropic) operated by the author. Claude also wrote checker K1 and the search and polishing code, re-proved the certificates and drafted the text. GPT-6 Astra (OpenAI) wrote checker K2 blind, carried out the efficiency study and the dimension computation, audited the literature, and fact-checked the manuscript. The author chose the project and its criteria, required cross-family review, chose the stopping rules and decided when to end the search, and re-ran both checkers on all certificates on the author's own computer (every check accepted). The written arguments were drafted with AI and have not yet been checked line by line by the author or by a mathematician. Anyone can check every certified method with one command (python reverify_d1.py). No AI system is an author; the author takes full responsibility for the content. Details are in the section "Use of artificial intelligence" of the paper.

Zenodo (CERN European Organization for Nuclear Research)
The Ohio State University (US)
Numerical methods for differential equations
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.

Laika access coming soon