Today, we’re introducing Épi, a self-evolving framework for autonomous research. Working in verifiable environments on our sandbox infrastructure, Épi’s agents produced 120 record-setting results in just a few days, including new constructions and stronger bounds. For example, in polyomino enumeration, Épi lowered the upper bound on Klarner’s constant from 4.5238 to 4.29569. It also proved that the Cohn–Elkies method cannot reach the exact sphere-packing density in dimensions 10 and 11, even at its theoretical limit. For Tuza’s conjecture, Épi produced an independently verified, computer-assisted proof of a new case, covering an infinite family of graphs.
These findings show that frontier models already have capabilities that can advance research. We built Épi to turn those capabilities into a sustained research process that carries verified results, useful methods, and unfinished questions forward. Our goal is to make discovery cumulative, with each advance expanding the scope and depth of future research.
Epi: Turning compute into progress
Épi brings together environments that evolve with research, agents that develop and test their own approaches, and infrastructure that keeps experiments running at scale.

Research that keeps going
Each session carries the investigation forward, preserving its best results, unfinished work, and the reasoning behind earlier decisions.
Strategies that evolve
Agents develop new approaches, revisit earlier ideas, and reshape their research strategy as evidence accumulates.
Experience that accumulates
Results, methods, and failed approaches form a growing research memory that future investigations can draw on.
Experiments in parallel
Explore multiple ideas at once, with persistent experiments that keep running as the investigation moves between sessions.
Progress you can verify
Every candidate faces the same fixed standard, making each accepted improvement a result you can check.
Independent review
An independent perspective challenges the agent’s assumptions and research direction, helping promising ideas stand up to scrutiny.
Environments that evolve
Progress starts with a robust verifier: a reliable way to check whether a candidate is valid and improves on the reference result. Within each investigation, Épi freezes the environment’s configuration and verifier, so every candidate faces the same test as the research changes direction.
Across investigations, the frontier moves as verified advances create new starting points and new research demands shape the environment itself. This is the core of Épi’s self-evolution: agents that make progress within environments, and environments that evolve with their discoveries.
Agents that develop their own methods
Because a difficult problem rarely reveals the right method at the outset, Épi’s agents develop their approaches as evidence accumulates, following a promising direction further or finding a different way into the problem. In the process, the investigation can change the method itself.
Research history remains available across sessions, allowing that development to continue, while independent review brings another perspective. An idea set aside earlier can become useful when the problem is understood differently.
Infrastructure for discovery at scale
Breadboard, the execution infrastructure powering Épi, combines high concurrency in isolated sandboxes with the stability that extended research requires. It sustains many investigations in parallel, so research can keep moving while its longest experiments are still running.
As experiments continue across sessions and across different timescales of computation, one investigation can follow a difficult question further while others open new directions.
Mathematical advances with Epi
Épi’s mathematical advances include new constructions, stronger bounds and a new computer-assisted proof. They also establish limits of the Cohn–Elkies method, showing where a stronger mathematical approach is needed.
| Research area | Problem families | Results |
|---|---|---|
| Numerical bounds | 12 | 20 records |
| Explicit constructions | 8 | 100 records |
| Tuza’s conjecture | 1 | 1 theorem extension |
| Total | 21 | 120 records + 1 theorem extension |
A problem family groups related instances. Each final result for a distinct instance or parameter setting counts once. A certified interval counts once, and repeated improvements to the same instance contribute one record.
Which problem families?
Numerical bounds (12): Boolean-function bounds, Cohn–Elkies linear programming, Spencer discrepancy, cubic-lattice dimers, n-queens, Ramsey multiplicity c₄,₅, MAX-4-CUT, polyomino growth, Bₕ[g] coefficients, Feigenbaum dimension, quantum capacity and water-network design.
Constructions (8): circles in a quadrant, balls in four dimensions, pentagons in a triangle, circles in a semicircle, circles in an octagon, vehicle routing, the third autocorrelation inequality and the Zhang–Zagier bound.
Proof (1): Tuza’s conjecture for a larger class of split graphs.
New constructions and stronger bounds
MAX-4-CUT: Épi reformulated the problem and lowered the NP-hardness ratio from 0.9870848895 to 84931/86101 ≈ 0.986411307651.13
Ramsey multiplicity c₄,₅: A construction on 1,024 vertices set a new upper bound of 4427841832181/2⁴⁹ ≈ 0.007865427122378.12
Spencer discrepancy: The advance fits inside an 8 × 8 matrix. All 256 sign assignments were evaluated exactly, giving an exhaustive certificate that raises the record lower bound from 7/√17 to 5/√8.9
Certified numerical results
N-queens constant: The constant α describes how quickly the number of non-attacking arrangements of n queens on an n × n chessboard grows. Nobel, Agrawal, and Boyd bounded α using large convex optimization problems.11 Épi extended their calculation and narrowed the enclosing interval to a width of 3.9367 × 10⁻⁹, approximately 84× narrower than the published interval. A separate verification checked all 67,108,864 projected coordinates of the construction and all 57,340 constraints. Feasibility was verified exactly, and directed rounding kept the objective bounds valid despite numerical error.
Klarner’s constant: This constant describes the growth in the number of shapes formed by joining squares edge to edge on a grid as more squares are added. Mathematicians have studied bounds on this growth for more than half a century. Building on Bui’s recurrence method, Épi lowered the upper bound from 4.5238 to 4.29569 by extending the calculation to 1,186 neighborhood types. An exact certificate covers the full recurrence system.14
Quantum capacity: Épi found 40-qubit input states that support positive capacity at noise parameters 0.1198803 for the independent X–Z channel and 0.1193415 for the Two-Pauli channel. Both exceed the published thresholds. A separate interval calculation confirms that coherent information remains strictly positive after accounting for numerical error.17
Feigenbaum dimension: For the quartic attractor, Épi obtained the certified dimension interval [0.6425750638035370, 0.6425750659339700], more than 3.6 million times narrower than the rigorous interval proved by Burbanks, Osbaldestin, and Thurlby. This comparison concerns the width of the proved enclosure; it does not measure a gain over every prior numerical estimate.16
Bₕ[g] coefficients: Épi improved all five coefficients for h = 3 through 7. These coefficients bound the size of sets with a limited number of representations of an integer as a sum. The search optimized cosine kernels, and interval bounds certified their minima across the entire domain.15
Water-network design: Épi raised the global lower bound for the waterund36 network-design problem to 662.80702839, within approximately 0.00001001 of the known feasible value. The release includes an exact certificate for the bound.20
Limits of the Cohn–Elkies method
The sphere-packing results cross a different kind of threshold: they establish where a leading mathematical method must fall short. The problem asks how densely identical spheres can fill a space. The Cohn–Elkies method puts an upper bound on that density and gives the exact answer in dimensions 8 and 24. Elsewhere, there can be a gap between what the method allows and what any packing can achieve.
Épi used Li’s discrete-reduction method to bound the best answer the Cohn–Elkies program could ever produce.8 It improved these lower bounds in dimensions 9, 10, 11 and 13. In dimensions 10 and 11, they crossed a decisive threshold: they exceeded known upper bounds on the actual packing density.18
This establishes a limitation of the method itself. Even if the Cohn–Elkies program were solved perfectly, its answer would remain above the true packing density by at least approximately 0.35% in dimension 10 and 1.92% in dimension 11. Closing those gaps requires a stronger method.
A new case of Tuza’s conjecture
Épi’s research extends from improving numerical bounds to proving new cases of mathematical conjectures. For Tuza’s conjecture, it produced a computer-assisted proof for graphs with a fully connected eight-vertex core and an independent set of other vertices. Among the outside vertices that form triangles, it allows up to three distinct ways of connecting to the core, extending Zeng’s two-type result.19 Each pattern can repeat arbitrarily many times, so the theorem covers an infinite family of graphs. It guarantees that if at most k triangles can be chosen without sharing an edge, at most 2k edges suffice to meet every triangle. The proof passed independent verification.
Twenty numerical bound records across twelve problem families
The table lists all 20 numerical bound records. They cover 12 problem families: Cohn–Elkies contributes four dimensions, Bₕ[g] contributes five parameter values, quantum capacity contributes two channels, and the other nine families contribute one result each.
The entries measure different kinds of progress: a tighter upper or lower bound, a narrower enclosing interval, or a higher noise threshold at which quantum communication remains possible. The direction of improvement is marked for each row.
| Constant or bound | Published bound | Épi |
|---|---|---|
| Wellens’ D∞7Upper bound ↓ | 4.3935 | 4.3654674227 |
| Cohn–Elkies LP, d = 98Lower bound on LP optimum ↑ | 1/20 | 0.058044329039 |
| Cohn–Elkies LP, d = 108Lower bound on LP optimum ↑ | 1/24 | 0.056433853144 |
| Cohn–Elkies LP, d = 118Lower bound on LP optimum ↑ | 1/24 | 0.057731808202 |
| Cohn–Elkies LP, d = 138Lower bound on LP optimum ↑ | 1/28 | 0.066440082613 |
| Spencer discrepancy9Lower bound ↑ | 7/√17 ≈ 1.697749 | 5/√8 ≈ 1.767767 |
| Cubic-lattice dimer constant10Upper bound ↓ | 0.452130 | 0.452010 |
| N-queens constant α11Enclosing interval | [1.944000752, 1.944001082] | [1.9440010223327, 1.9440010262694] |
| Ramsey multiplicity c₄,₅12Upper bound ↓ | 2129191/2²⁸ ≈ 0.007931854575872 | 4427841832181/2⁴⁹ ≈ 0.007865427122378 |
| MAX-4-CUT hardness ratio13Upper bound ↓ | 195581/198140 ≈ 0.9870848895 | 84931/86101 ≈ 0.986411307651 |
| Polyomino growth constant λ14Upper bound ↓ | 4.5238 | 4.29569 |
| Bₕ[g] coefficient, h = 315Upper-bound coefficient ↓ | 14.3 | 14.287390218082 |
| Bₕ[g] coefficient, h = 415Upper-bound coefficient ↓ | 71.49 | 70.289552594834 |
| Bₕ[g] coefficient, h = 515Upper-bound coefficient ↓ | 413.07 | 408.104453996958 |
| Bₕ[g] coefficient, h = 615Upper-bound coefficient ↓ | 2774.16 | 2749.492485514635 |
| Bₕ[g] coefficient, h = 715Upper-bound coefficient ↓ | 21294.74 | 21149.189776156220 |
| Quartic Feigenbaum attractor dimension16Enclosing interval | [0.6395131468772885, 0.6473156929016112] | [0.6425750638035370, 0.6425750659339700] |
| Quantum capacity, independent X–Z17Positive-capacity noise threshold ↑ | 0.11837113922723974 | 0.1198803 |
| Quantum capacity, Two-Pauli17Positive-capacity noise threshold ↑ | 0.1180672806619063 | 0.1193415 |
| Water-network design, waterund3620Global lower bound ↑ | 655.5441707 (MINLPLib) | 662.80702839 |
Decimal bounds are rounded outward; fractions and radicals give exact values. Bₕ[g] entries are the coefficients inside the h-th root, compared with Timmons’ Table 1.
Sustained research in practice
The evaluations below examine how Épi improves existing results, from successive gains within a single investigation to outcomes across a broader set of problems.
Seventeen improvements in 56 minutes
Épi improved a published construction for packing 34 regular pentagons into an equilateral triangle. Starting from Hyra’s reference arrangement, credited to Haowei Lin on Friedman’s pages in July 2026, it reduced the triangle’s side length from 12.98141 to 12.97564 through 17 successive improvements in a 56-minute session.6
All 17 candidates passed the same fixed verifier, keeping the standard of validity and improvement consistent throughout the search. The total reduction reached 19 times the minimum improvement required for this case.
The first two candidates accounted for 82.4% of the improvement within 9.4 minutes. Épi continued searching after these early gains, producing another 15 verified improvements that reduced the side length from 12.97665 to 12.97564. Together, those later results contributed the remaining 17.6% of the total reduction.
One hundred construction records across eight problem families
Épi set 100 construction records across eight problem families, covering geometric packing, vehicle routing, and analytic inequalities. These include 46 circle packings in a quadrant, 36 ball packings in four dimensions, nine pentagon packings, five circle packings in a semicircle, and one record in each of the remaining four families.
We compared the results against published tables and registries, including Packomania, Erich Friedman’s packing pages, CVRPLIB, EinsteinArena, and the mathematical optimization registry. The constructions, mathematical certificates, and complete result tables are available in Épi Results.
Building on Hyra’s constructions
Épi also built on recent AI-generated constructions, improving nine pentagon packings released by Tencent’s Hyra. In addition, Épi reduced routing overhead across Hyra’s collection of 61 quantum circuits by 4,140 CNOT gates, or 1.60%.5 The constructions are available in Épi Results.
| Construction | Hyra | Épi |
|---|---|---|
| 34 pentagons Triangle side length ↓ | 12.98141 | 12.97564 |
| 61 quantum circuits Added CNOTs ↓ | 258,366 | 254,226 |
Comparison of released constructions. Routing totals are recomputed from the circuit files, with three CNOTs per SWAP.
When the search finds no improvement
In this failure case, Épi investigated flat polynomials for seven hours, exploring 31 approaches, abandoning 27, and submitting 20 candidates for evaluation. None improved on the published construction. Every candidate had a higher objective value than the reference, on a task where lower values are better.
The tree below traces how these approaches developed, showing where the search branched and which directions were abandoned. Four directions remained open when the run ended.
Even when a search finds no improvement, Épi preserves what was tried and which approaches were abandoned. Later investigations can review this history when deciding what to revisit and what to try next.
Research that builds on itself
These 120 records and the Tuza theorem extension are a starting point. Mathematics is our first testing ground, and we are expanding Épi into other domains, including research on recursive self-improvement. As we develop its agents, verifiable environments, and sandbox infrastructure together, our goal is for Épi to become a more effective research system through the work it carries out.
Our next evaluations will examine how efficiently Épi explores different approaches and whether it can combine discoveries to address more complex research questions. We also want to move beyond tasks whose objectives are already defined, toward identifying worthwhile questions and deciding where further investigation is most likely to matter.
In French, an épi is an ear of grain, carrying the seeds of the next season. It captures our hope for Épi: that each discovery leaves behind methods, ideas, and new questions from which further research can grow.
Models provide capability. Environments make progress verifiable. Épi turns compute into improvement.
We work with teams building verifiable research environments, deploying research agents, and running controlled experiments at scale on sandbox infrastructure that provides unified execution and traceable results. Get in touch to work with us.
Related work
Épi builds on a growing body of work on models as research systems. FunSearch uses execution feedback to guide program search.1 AlphaEvolve extends evolutionary coding to scientific and algorithmic problems.2 The AI Scientist-v2 coordinates experiments through agentic tree search.3 TTT-Discover trains the model during its search.4 Hyra has released constructions across mathematics, scientific computing and AI research.5 Our focus is the environment, agent and infrastructure that let this work continue across investigations, with a verifiable record of what has improved. Most of the results reported here were obtained with OpenAI’s GPT-5.6 Sol and GPT-6-Astra.
References
Research agents and systems
- Bernardino Romera-Paredes et al.. Mathematical discoveries from program search with large language models.
- Alexander Novikov et al.. AlphaEvolve: A coding agent for scientific and algorithmic discovery.
- Yutaro Yamada et al.. The AI Scientist-v2: Workshop-Level Automated Scientific Discovery via Agentic Tree Search.
- Mert Yuksekgonul et al.. Learning to Discover at Test Time.
- Hyra Team. Hyra: Hunyuan Research Agent.
Mathematical sources
- Erich Friedman. Pentagons in Triangles.
- Jake Wellens. Relationships between the number of inputs and other complexity measures of Boolean functions.
- Rupert Li. Dual Linear Programming Bounds for Sphere Packing via Discrete Reductions.
- Terence Tao and contributors. A collection of optimization problems in mathematics.
- Qidong He. A new upper bound on the dimer constant of ℤ³.
- Parth Nobel, Akshay Agrawal and Stephen Boyd. Computing tighter bounds on the n-queens constant via Newton’s method.
- Olaf Parczyk, Sebastian Pokutta, Christoph Spiegel and Tibor Szabó. New Ramsey Multiplicity Bounds and Search Heuristics.
- Ansh Nagda, Prabhakar Raghavan and Abhradeep Thakurta. Reinforced Generation of Combinatorial Structures: Hardness of Approximation.
- Vuong Bui. A convolutional approach to bounding the number of polyominoes.
- Craig Timmons. Upper bounds for Bₕ[g]-sets with small h.
- Andrew D. Burbanks, Andrew H. Osbaldestin and Judi A. Thurlby. Rigorous bounds on the Hausdorff dimension of Feigenbaum attractors.
- Avantika Agarwal et al.. Enhanced quantum capacity thresholds from symmetry.
- Henry Cohn, David de Laat and Andrew Salmon. Three-point bounds for sphere packing.
- Zijian Zeng. Tuza’s Conjecture for Split Graphs with an Eight-Vertex Clique Part and Two Neighborhood Types.
- MINLPLib. waterund36: primal and dual bounds.
Cite this work
Epi Team. (2026). Épi: Turning compute into verifiable improvement. Bake AI.
https://bakeai.inc/research/articles/epi/
@misc{epiteam2026epi,
author = {{Epi Team}},
title = {{\'E}pi: Turning compute into verifiable improvement},
year = {2026},
month = sep,
howpublished = {Bake AI},
url = {https://bakeai.inc/research/articles/epi/}
}
Records may change over time. If any result here is out of date, please let us know.

