Nuclear & Fusion PhysicsResearch AnalysisWorking Draft · 20 Oct 2026

Logarchéon and Nuclear Physics: Evaluation and Funding Routes

Logarchéon's AI designs draw on statistical physics, hyperbolic geometry and causal inference, the same mathematics used in nuclear physics, fusion and gravitation. Some of the links are exact identities; others are analogies that still require testing. Each link is graded, followed by the national-security frame and the funding routes.

Terms

  • Critical-entropy attention (CEAS): sets the "temperature" of each attention unit so its output stays near a phase transition, where it is most sensitive.
  • Ψ causal engine: answers cause-and-effect questions from data, and states when the data cannot decide (abstention).
  • GRAIL and MIA: attention computed from distances in hyperbolic space, so results are unchanged under that space's symmetries, including Lorentz boosts.
  • Two-phase scheme: methods learned on small models are applied to a large model by changing only a small set of control parameters.
  • Twin map; DL ≠ LD: a change of parameters that leaves a model's results unchanged; training does not commute with it ("DL ≠ LD"), so applying it changes how training proceeds.
  • Trace store; gate: a store that cross-checks spectra through trace identities, and a statistical test that accepts a change only when its benefit is bounded with stated confidence.

Standard of evidence

A link is listed only when a Logarchéon mechanism has a named counterpart in physics. Results proved or measured in Logarchéon's tests count for more than analogies. The work stays within open, unclassified science: nuclear structure and reactions, nuclear data, emulators, collision data and nonproliferation analytics.

Summary of findings

The overlap with nuclear physics is real but narrower than "almost everything". Several of these tools came from nuclear physics originally: random-matrix theory, for example, was first used for nuclear energy levels. The most promising contributions are those already proved or measured:

  1. DL ≠ LD. Symmetry-generated transformations and preconditioning, with a proved bound on error amplification.
  2. Certified trace checks. The trace store as a check of semiclassical shell structure.
  3. Identification with abstention. The causal engine for nuclear data evaluation.
  4. Gated, traceable updates. The gate and hash chain as verification and provenance.
  5. Exact Lorentz invariance. GRAIL/MIA for collision data.
  6. Compression. Large models on one processor, on site.

Where the designs map onto nuclear physics

Rows are ordered by how likely each is to be a contribution rather than a re-use.

Logarchéon pieceNuclear counterpartStatus
DL ≠ LD: twin jumps, commutator η(L − L−⊤)∇ℓ, boost test eφTransforming the Hamiltonian before an approximate solve (similarity renormalisation group, in-medium, transcorrelated, Lee–Suzuki); bosonic Bogoliubov transformations; preconditioning of variational wavefunctions; covariant fission-path searchesMost promising. New: a proved bound on error growth for non-unitary transformations, and a method for choosing the transformation
Trace store, Selberg trace formula, heat traces Σe−tλPeriodic-orbit theory of shell structure (Gutzwiller, Balian–Bloch, Strutinsky); the heat trace is the nuclear partition function and gives the level densityStrong. A certified check of semiclassical shell corrections against exact spectra
Ψ: identification with abstention; inverse design with infeasibility certificatesNuclear data evaluation, including optical-model ambiguities; attributionDirectly applicable
Gate, hash chain, Lean certificatesVerification, validation and provenance of evaluated data librariesAddresses the verification and traceability that evaluated-data programs require
GRAIL/MIA: Lorentz group, hyperboloidRelativistic kinematics: rapidity is hyperbolic distance, the mass shell is H³; Electron-Ion Collider and heavy-ion dataApplicable. Existing Lorentz-equivariant networks (LorentzNet, PELICAN) set the benchmark to meet
Two phases: small → frozen large through a control setEmulators (eigenvector continuation, reduced basis); the light-to-heavy-nuclei gapUntested; a direct test is defined
Langlands: Hecke operators on GRAIL's arithmetic group; joint Hecke–Laplace eigenfunctionsGlobal spectral training of a GRAIL model of nuclear data (structure, reactions, levels) instead of gradient steps, combined with the two-phase schemeUntested on nuclear data. Hecke operators exist only for arithmetic groups. On a 32-billion-parameter model, Hecke attention ran at a cost of about 0.2% before any training
CEAS: Gibbs rows, freezing point, finite-size criticality, Lee–Yang zerosPhase transitions in finite nuclei (pairing, shape, liquid–gas)The mathematics applies. The model prediction it rests on matched only 1–2 of 12 trained attention units, so any nuclear claim must be tested on nuclear data
Spectral face: GOE/Poisson spacing ratios 0.531 and 0.386, Marchenko–Pastur, BBPRandom matrix theory, which began with nuclear resonancesIncremental: automated, statistically gated analysis of level data
Compression: 32B model on one T4On-site analysis where data cannot leaveRuns a 32-billion-parameter model on one standard GPU (T4), for analysis where data cannot leave the site

DL ≠ LD: tunnelling in the Logarchéon sense

In Logarchéon, tunnelling is the twin jump and its non-commutation with training. Let D be a gradient step and L a twin map, an invertible map with ℓ(Lθ) = ℓ(θ). The following identity is proved:

D(Lθ) − L·D(θ) = η (L − L−⊤) ∇ℓ(θ)

The commutator vanishes exactly when L is orthogonal; for a Lorentz boost of rapidity φ its size is 2 sinh φ. A twin jump reaches a point of equal loss at no cost, and training continues from there under the preconditioner (LᵀL)⁻¹. Mapped back, a boost amplifies error by at most eφ.

Three nuclear counterparts:

  1. Transforming before an approximate solve. Ab initio theory applies transformations that keep exact energies but speed the truncated calculation. Unitary ones are the orthogonal case; non-unitary ones are the boosts, and the eφ bound is a guarantee they lack. Bosonic Bogoliubov transformations have the boost's form, u = cosh r, v = sinh r.
  2. Variational wavefunctions for light nuclei. Neural-network wavefunctions need a good preconditioner (stochastic reconfiguration). A twin jump is a preconditioner change obtained from a symmetry.
  3. Fission paths. The action S = ∫√(2M(q)(V − E)) ds does not depend on the shape coordinates, but gradient-based path searches do unless they use the inertia tensor M as metric. The covariant transport rule is that requirement; M plays the role of LᵀL.

Status. The commutator is classical: it changes how descent moves, not barrier penetration with amplitude e−S. On small character models, twin jumps changed the final loss by under 0.001. The claim to test is symmetry-generated preconditioning or transformation with a proved bound.

K9 faces when the model learns nuclear physics

K9 is the network of Logarchéon's nine coupled components. When the tokens are nucleons, levels or reactions instead of words, the faces take a physical meaning.

FaceMeaning for nuclear data
SpectralModel spectra compared with measured level statistics
Number theory / SymbolicCounting shell-model configurations; level density through the trace store
CEASTemperature of the nuclear ensemble
MIA/GRAILExact Lorentz invariance for collision data; twins as transformations that keep observables and speed the approximate solve (DL ≠ LD)
ΨEvaluation of reaction and structure models; attribution
RSI with the gateModel refinement accepted only on statistical evidence; a twin jump is one proposed move
ArithmeticProvenance of the data library (hash chain, fingerprints)
LanguagePrograms of moves reused from light to heavy nuclei (the two phases)
LanglandsGlobal spectral access to the model's own geometry: with an arithmetic group, Hecke and Laplace eigenfunctions give a spectral basis in which training can be solved globally instead of by gradient steps. This applies to any data GRAIL models, nuclear data included. Untested on nuclear data.

From language to nuclei: what each stage means

In a language model the units are words. Here the units are nuclei: each token is a nucleus (proton number Z, neutron number N), one of its energy levels, or a reaction. The stages then read as follows.

  1. Data. Measured masses, levels and cross sections become the sequence the model reads, in place of text.
  2. Geometry (GRAIL). Every token and every parameter is a point in hyperbolic space, and the model's scores are distances between points. Training moves points, so the loss is a landscape over that geometry. With a discrete group Γ, it is a landscape on the quotient Γ\ℍⁿ, a finite space with the symmetries built in.
  3. Twin maps (Lorentz group). Moving all points by the same symmetry leaves every prediction unchanged. This is like changing the basis or coordinates in a nuclear calculation: the physics is the same, but how the calculation converges is not (DL ≠ LD).
  4. Laplace eigenfunctions: the normal modes of the landscape. Any function on the landscape, the loss included, can be written as a sum of the geometry's normal modes, as a vibration is written as a sum of harmonics. The eigenvalue says how fast a mode varies. Low modes carry broad structure across the whole dataset; high modes carry local detail. For nuclear data, the hypothesis is that smooth trends across the chart of nuclides fall in the low modes and shell effects in higher ones. This is untested.
  5. Hecke operators: extra conserved labels. A Hecke operator T_p averages each point with p + 1 related points of the geometry. It commutes with the Laplacian. In a nucleus, angular momentum and parity commute with the energy, so states of equal energy are told apart by those quantum numbers. In the same way, Hecke eigenvalues tell apart modes that the Laplacian alone cannot. The joint Hecke–Laplace eigenfunctions are the landscape's modes with all their labels: its "quantum numbers".
  6. Langlands: global access to the spectrum. By the Langlands correspondence and the trace formula, these labels and frequencies match arithmetic data: counts of closed geodesics, traces of matrices, values at primes. The whole spectrum of the landscape can therefore be read from counting, without exploring the landscape point by point. Training can then fit nuclear data by choosing the coefficient of each labelled mode, solved globally, instead of by many small gradient steps.
  7. Two phases. Modes and programs found on light nuclei, where exact calculations exist, are applied to heavy nuclei by changing only a small set of control parameters.

Status on a 32-billion-parameter model. Stages 1–3 are built and tested on language data. On an open-weight language model of 32.8 billion parameters, held at 2.6 bits per weight in 11.8 GB of memory:

  • all 4,096 attention heads were lifted into hyperbolic geometry with no retraining, at a cost of 0.65% or less in perplexity up to moderate curvature (16.84 → 16.79–16.95);
  • at the strongest curvature tested, perplexity rose to 40.17, and 60 steps of training on the control set brought it back to 18.14; the statistical gate has not yet accepted this, and nearly all of the remaining cost lies in the early layers;
  • Hecke attention ran on lifted heads at a cost of about 0.2% before any training;
  • the CEAS temperature program, accepted by the gate, recovered about 5% of the loss from compression, and a fact update through an adapter passed checks on unrelated text and on 200 general-knowledge questions.

Stages 4–7 on nuclear data are untested. Hecke operators exist only when Γ is arithmetic.

From language to nuclei: what each stage means

In a language model the units are words. Here the units are nuclei: each token is a nucleus (proton number Z, neutron number N), one of its energy levels, or a reaction. The stages then read as follows.

  1. Data. Measured masses, levels and cross sections become the sequence the model reads, in place of text.
  2. Geometry (GRAIL). Every token and every parameter is a point in hyperbolic space, and the model's scores are distances between points. Training moves points, so the loss is a landscape over that geometry. With a discrete group Γ, it is a landscape on the quotient Γ\ℍⁿ, a finite space with the symmetries built in.
  3. Twin maps (Lorentz group). Moving all points by the same symmetry leaves every prediction unchanged. This is like changing the basis or coordinates in a nuclear calculation: the physics is the same, but how the calculation converges is not (DL ≠ LD).
  4. Laplace eigenfunctions: the normal modes of the landscape. Any function on the landscape, the loss included, can be written as a sum of the geometry's normal modes, as a vibration is written as a sum of harmonics. The eigenvalue says how fast a mode varies. Low modes carry broad structure across the whole dataset; high modes carry local detail. For nuclear data, the hypothesis is that smooth trends across the chart of nuclides fall in the low modes and shell effects in higher ones. This is untested.
  5. Hecke operators: extra conserved labels. A Hecke operator T_p averages each point with p + 1 related points of the geometry. It commutes with the Laplacian. In a nucleus, angular momentum and parity commute with the energy, so states of equal energy are told apart by those quantum numbers. In the same way, Hecke eigenvalues tell apart modes that the Laplacian alone cannot. The joint Hecke–Laplace eigenfunctions are the landscape's modes with all their labels: its "quantum numbers".
  6. Langlands: global access to the spectrum. By the Langlands correspondence and the trace formula, these labels and frequencies match arithmetic data: counts of closed geodesics, traces of matrices, values at primes. The whole spectrum of the landscape can therefore be read from counting, without exploring the landscape point by point. Training can then fit nuclear data by choosing the coefficient of each labelled mode, solved globally, instead of by many small gradient steps.
  7. Two phases. Modes and programs found on light nuclei, where exact calculations exist, are applied to heavy nuclei by changing only a small set of control parameters.

Status. Stages 1–3 are built and tested on language data. Stages 4–7 on nuclear data are untested. Hecke operators exist only when Γ is arithmetic.

Status on a 32-billion-parameter model (20 October 2026)

The stages were run on an open-weight language model of 32.8 billion parameters, compressed to 2.6 bits per weight, with each change accepted only when a statistical bound passed on held-out text. Perplexity measures prediction error; lower is better.

StageWhat was doneResult
Geometry (GRAIL)All 4,096 attention heads lifted to hyperbolic geometry, no retrainingAt moderate curvature, perplexity changes by 0.65% or less (16.84 → 16.79–16.95); at strong curvature it rises to 40.17
Geometry (GRAIL)Recovery at strong curvature: 60 steps of training on per-head temperatures and a small adapter, the stored weights unchangedPerplexity back from 40.17 to 18.14, about 91% of the loss recovered; not yet within the acceptance bound
Geometry (GRAIL)Strong curvature applied to one third of the layers at a timeAlmost all of the cost lies in the early layers: +0.446 nats per token, against +0.011 (middle) and +0.027 (late)
Twin mapsBoosts and Hecke matrices applied to lifted pointsDistances preserved to 7 × 10⁻¹⁵
Hecke operatorsOrbit terms and Hecke operators T₂ and T₃ on six lifted heads, no trainingCost of 0.2% or less (perplexity 16.845–16.868 against 16.839)
Normal modes and LanglandsGlobal spectral trainingNot yet run on this model
Two phases (control set)CEAS temperature program, accepted by the gateRecovers about 5% of the loss from compression, the stored weights unchanged
Two phases (control set)Fact update through a small adapter (0.1% of parameters)Corrected answer preferred on 5 of 5 unseen questions (2 of 5 before); no measurable harm on unrelated web text or on 200 general-knowledge questions

Next. Curvature set by band of layers: strong in the middle and late layers, small in the early ones.

Fusion

Magnetic-confinement fusion gives two exact links, not analogies: Greene's residue and the poloidal-angle twin.

Inertial-confinement fusion overlaps with stockpile stewardship, so work stays on magnetic confinement and open fusion-energy data.

Logarchéon pieceFusion counterpartStatus
SL(2,ℝ) classification by trace; trace store counting periodic pointsThe field-line map of a torus is area-preserving; its linearisation at a periodic field line lies in SL(2,ℝ). Greene's residue R = (2 − Tr M)/4 tells an island O-point (0 < R < 1) from an X-point (R < 0, or R > 1 for an X-point with reflection), and is used to control islands in stellarator designExact identity. The trace store's periodic-point count is the count of periodic field lines
DL ≠ LD: twin mapsChanging the poloidal angle of the Fourier series of a flux surface leaves the surface and the loss unchanged, yet gradient optimisers (DESC, SIMSOPT) converge differentlyExact twin. Choosing the angle by an orbit-optimal twin through the gate is new and testable
Ψ: identification with abstentionEquilibrium reconstruction from external magnetic sensors is non-unique: different internal profiles fit the same dataDirectly applicable: a certified "not identifiable" answer
Two phases: small → large through a control setDisruption predictors trained on small tokamaks transfer poorly to ITER or SPARCStrong and testable
Gate, control-barrier functionsLearned plasma controllers (run on TCV); updates must be safeAddresses the safety requirement for controller updates
BBP spikes over the Marchenko–Pastur edgeCoherent MHD and Alfvén modes above noise in magnetic-sensor arraysModerate
CEAS near criticalityProfiles held near a critical gradient, with avalanchesAnalogy until tested
Hyperboloid geometryRelativistic runaway electrons during disruptionsModerate
Ψ on nuclear dataNeutron cross sections and activation of fusion materialsFits the joint fusion/nuclear-data call
Langlands and Hecke operatorsGlobal spectral training of a GRAIL model of plasma data, as for nuclear dataUntested

First tests. (1) The poloidal-angle twin in an open-source stellarator optimiser: one laptop or free GPU. (2) Greene's residue through the trace store on a computed field-line map: almost pure computation.

Funding. Fusion Energy Sciences lists "AI/ML for fusion and plasma science" in the Office of Science open call (FY26 window closed 30 September 2026; successor not yet confirmed). Fusion-materials data were the main topic of the nuclear-data interagency call. INFUSE vouchers go to fusion companies, so they need a fusion-company partner.

Gravitation (gr-qc)

Gravity in 2+1 dimensions (AdS₃, whose Euclidean form is hyperbolic 3-space) gives exact links, including the Hecke side that has no nuclear use. Logarchéon works with the Lorentz group O(n−1,1) in any dimension, so the image-sum and trace-formula links extend to higher-dimensional AdS quotients and to 3+1 cosmology; the Hecke links are specific to AdS₃.

Logarchéon pieceGravitation counterpartStatus
GRAIL orbit sum Σ_γ e−βd(q,γk); ping-pong admission of Schottky generatorsBTZ black holes and multi-boundary wormholes are quotients of hyperbolic space by discrete (Schottky) groups (Brill; Krasnov). Fields on them are sums over images, G(x, y) = Σ_γ G₀(x, γy)Exact identity: the kernel is an image-sum propagator on the quotient
Selberg trace formula and zeta; resonances by the noisy Ulam methodBTZ quasinormal modes are zeros of a Selberg zeta function (Perry–Williams)Exact. Logarchéon's trace and resonance pipeline applies to black-hole spectra
Patterson–Sullivan critical exponent δ (measured 0.32127 for a rank-3 Schottky group)δ controls convergence of the image sums and which handlebody geometry dominatesExact
Hecke operators; Poincaré series over SL(2,ℤ)Pure-gravity partition function in AdS₃ as a sum over SL(2,ℤ) images (Maloney–Witten); symmetric-orbifold partition functions of the D1–D5 system generated by Hecke operators (Dijkgraaf–Moore–Verlinde–Verlinde)Exact: Langlands-side tools with a gravity meaning
Ψ: identification with abstentionGravitational-wave parameter inference, with known degeneracies such as distance and inclinationModerate to strong
Two phasesWaveform surrogates trained on numerical-relativity runs and extended to new mass ratiosModerate
DL ≠ LDGauge choice in numerical relativity: physics is gauge-independent, discretised evolution is notAnalogy: the gauge freedom there is infinite-dimensional
CEAS, spectral faceBlack-hole chaos: random-matrix spacing statistics and the JT-gravity matrix model (Saad–Shenker–Stanford)Moderate
Logarchéon piece in general dimensionCounterpart beyond 2+1Status
Orbit sums over discrete groups of O(n−1,1)Black holes in AdS of any dimension made by identifying points under a discrete group (Aminneborg et al.); their Euclidean forms are quotients of hyperbolic spaceImage sums extend exactly; a quasinormal-mode–zeta identity beyond 2+1 is an open question
Trace store and Selberg trace formula on quotients of hyperbolic 3-space3+1 cosmology with compact hyperbolic spatial sections: Laplace eigenmodes and CMB fluctuations (Aurich, Lustig, Steiner; eigenmode computation)Exact mathematics in 3+1

Not carried. The wormhole-nucleation conjecture, Lee–Yang wormholes and the gravitational-wave "sensing" proposals, which remain speculative and, by sensitivity estimates, in places undetectable; and Coleman–De Luccia or Hawking–Page tunnelling, which Logarchéon does not compute.

Venue and first test. AdS₃ work goes to hep-th cross-listed to gr-qc; gravitational-wave inference to gr-qc. First test: BTZ quasinormal modes through the Selberg and resonance pipeline, compared with the known closed form.

A dimensional map: from nuclear origins to new scales

Planck's constant entered physics in 1900 as a fitting parameter in the blackbody law. The relation it sat in, E = hν, fixed its dimension as energy × time, that of action or angular momentum, before its meaning was understood. Angular momentum was a known quantity of classical mechanics. h is not a classical angular momentum, but the shared dimension showed where to look: Bohr quantised the electron's angular momentum in units of h/2π in 1913, and electron spin, ħ/2, followed in 1925. The same reasoning applies here. Each candidate's acquired dimension identifies it with a known quantity, a temperature, a mass, an energy or a quantum number, which guides what to measure and what to build, without claiming that the candidate is that quantity. Each Logarchéon construct is a mathematical method; placed in a nuclear relation, its parameters acquire physical dimensions, and those that are universal across independent data are new physical scales.

Symbols

SymbolMeaningDimension
ħ, kBReduced Planck constant and Boltzmann constant, both set to 1, as is the speed of light—
Ej, Ē, σEnergy of state j; centroid (mean) energy of a model space; width (standard deviation) of its energiesMeV
NIn formulas, the number of states of a model space; in "nucleus (Z, N)", the neutron number1
ZIn formulas, the partition function Σj e−βEj; in "nucleus (Z, N)", the proton number1
β; 1/βInverse temperature; temperatureMeV⁻¹; MeV
S, SmaxEntropy of the ensemble; its maximum, ln N1
CHeat capacity, β² Var(E)1
TGap score, βσ/√(2 ln N); T = 1 is the freezing point1
TfFreezing temperature, σ/√(2 ln N)MeV
Tc, ΔCritical temperature of the pairing transition; pairing gapMeV
A; a; E*; ρMass number; level-density parameter (≈ A/8 MeV⁻¹); excitation energy; level density1; MeV⁻¹; MeV; MeV⁻¹
ħωSpacing of the major nuclear shells, ≈ 41 A−1/3 MeVMeV
sj; q, kAttention score of key j; query and key vectors of an attention headas mapped by ι
cCurvature of the hyperbolic (lifted) attention of a headMeV⁻² when q is a momentum
κDimensionless curvature, c⟨|q|²⟩1
ΛMass scale c−1/2MeV
p, m, v, γThree-momentum; mass; speed as a fraction of the speed of light; Lorentz factor 1/√(1 − v²)MeV; MeV; 1; 1
meElectron mass, 0.511 MeVMeV
φ; eφRapidity of a Lorentz boost; relativistic Doppler factor1
θ; ℓ; ∇ℓModel parameters; loss (a variational energy under ι); its gradient1; MeV; MeV
L; L−⊤Twin map, an invertible linear map of the parameters; its inverse transpose1
η; K; τStep size of gradient descent; number of steps; imaginary time KηMeV⁻¹; 1; MeV⁻¹
H; G; E₁Hamiltonian; quantum geometric tensor; energy of the first excited stateMeV; 1; MeV
ℍⁿ; Γ; Γ\ℍⁿHyperbolic space of dimension n; a discrete group of its isometries; the quotient space—
PSL(2,ℤ)The modular group, the standard arithmetic group acting on ℍ²—
ΔΓ; λ; rLaplacian on Γ\ℍⁿ; its eigenvalue; spectral parameter, λ = 1/4 + r² on ℍ²1
h, gTest function of the trace formula and its Fourier transform—
γ (in sums over Γ); ℓγ; ℓ₁A group element or closed geodesic; its length; the shortest length1
tγ; ε₀Period of the orbit of γ; energy per unit geodesic lengthMeV⁻¹; MeV
Tp; λpHecke operator at the prime p; its eigenvalue1
μ; b, α; Δt; γRGRenormalisation scale; constants of the depth law; renormalisation-group time per layer; anomalous dimensionMeV; 1; 1; 1
ιDictionary from model variables to physical quantities—
S₀; ΔEConfigurational entropy of a protein's misfolded states; width of their energies1; energy
Tg; Tfold; TphysGlass temperature; folding temperature; physical temperature 1/βenergy

Origins: the tools came from nuclear physics

Tool used in LogarchéonNuclear originYear
Gibbs ensemble, partition function, temperatureLevel density of the Fermi gas (Bethe); evaporation (Weisskopf); compound-nucleus decay (Hauser–Feshbach)1936, 1937, 1952
Random matrices, level spacingsWigner's random-matrix model of neutron resonances; Porter–Thomas widths; Nuclear Data Ensemble fits GOE (Bohigas–Haq–Pandey)1951–1956, 1956, 1983
Gaussian state densitiesStatistical spectroscopy of shell-model spaces; two-body random ensembles (French–Wong; Bohigas–Flores)1970–1971
Phase transitions in finite systemsPairing in nuclei (Bohr–Mottelson–Pines); lattice-gas zeros (Lee–Yang)1958, 1952
Energy landscapesFission barrier (Bohr–Wheeler); shell correction (Strutinsky)1939, 1967
TunnellingAlpha decay (Gamow)1928
Periodic-orbit trace formulasShell structure from periodic orbits (Balian–Bloch; Gutzwiller; Strutinsky–Magner)1970–1976
Symmetry groups, quantum numbersWigner SU(4); Elliott SU(3); interacting boson model U(6) (Arima–Iachello)1937, 1958, 1975
Effective field theory, renormalisation groupYukawa meson; chiral EFT (Weinberg); similarity RG for nuclear forces (Bogner–Furnstahl–Perry)1935, 1990, 2007
Lorentz group, rapidityRelativistic heavy-ion kinematics (Bjorken)1983
Inverse problemsR-matrix theory (Wigner–Eisenbud); optical-model ambiguities (Igo)1947, 1958

Definitions

  1. Units. Natural units ħ = c = kB = 1. Every dimension is a power of energy, [X] = MeVd, with lengths converted by ħc = 197.327 MeV·fm.
  2. Construct. A Logarchéon construct is a map with parameters: the Gibbs selector pj(β) = eβsj/Z, the gap score T, the lifted score with curvature c, the Laplacian ΔΓ and Hecke operators Tp on Γ\ℍⁿ, the twin map L and the gradient step D(θ) = θ − η∇ℓ.
  3. Dictionary. A dictionary ι assigns to each model variable a physical quantity: scores to minus energies, sj ↦ −Ej; query and key vectors to four-momenta on a mass shell; the loss to a variational energy.
  4. Acquired dimension. The dimension a parameter must carry for every relation containing it to be dimensionally homogeneous under ι.
  5. Candidate scale. A parameter whose acquired dimension is non-zero.
  6. Physical constant. A candidate scale whose fitted value is universal: equal, within stated uncertainty, across independent systems and observables. Planck's h met this test: one value fitted blackbody spectra (1900) and the photoelectric effect (1905).

Postulates

  1. Homogeneity. Every argument of exp, log, sinh, cosh and of a softmax is dimensionless; both sides of every equation carry the same dimension.
  2. Fidelity of the dictionary. Under ι, a construct computes the physical quantity it is mapped to; for example the Gibbs selector computes Boltzmann occupation probabilities.
  3. Empirical decision. Dimension analysis fixes dimensions only. Whether a candidate scale is a physical constant is decided by the universality test of Definition 6, never by the analysis itself.

Common notion. Quantities equal to the same quantity carry the same dimension.

Propositions: the dimensions the constructs acquire

Proposition 1 (inverse temperature). Under sj ↦ −Ej, [β] = MeV⁻¹, and 1/β is a nuclear temperature.

Proof. The selector contains e−βEj; by Postulate 1 the exponent βEj is dimensionless, so [β][E] = 1. ∎

Proposition 2 (the gap score is a heat-capacity ratio). Let a nuclear model space of N states have a Gaussian density of states of centroid Ē and width σ, as finite shell-model spaces have to good approximation (French–Wong). Then, while the mean energy lies above the ground-state edge,

Smax − S(β) = C(β)/2, T² = β²σ²/(2 ln N) = C/(2Smax) = 1 − S/Smax

with Smax = ln N and heat capacity C = β² Var(E). Freezing, T = 1, occurs exactly when the mean energy reaches the expected lowest of N independent Gaussian levels to leading order in ln N, Ē − σ√(2 ln N), at the temperature

Tf = σ / √(2 ln N)

Proof. ln Z = ln N − βĒ + β²σ²/2, so ⟨E⟩ = Ē − βσ², Var(E) = σ² and C = β²σ². Then S = ln Z + β⟨E⟩ = ln N − β²σ²/2 = Smax − C/2; divide by Smax = ln N. Setting ⟨E⟩ = Ē − σ√(2 ln N) gives β = √(2 ln N)/σ, that is T = 1. ∎

Numbers. For ²⁸Si in the sd shell, N = 93,710 (M-scheme), ln N = 11.448 and Tf = 0.209 σ. The gap score therefore acquires a meaning in nuclear thermodynamics: T² is half the heat capacity in units of the maximal entropy. It is comparable with the pairing critical temperature Tc = 0.567 Δ ≈ 6.8 A−1/2 MeV (BCS, with Δ ≈ 12 A−1/2 MeV), and with heat capacities measured by the Oslo method (Schiller et al., 2001).

Proposition 3 (curvature acquires a mass). Under the dictionary that maps query and key vectors to spatial momenta, [c] = MeV⁻² and Λ = c−1/2 is a mass. On the mass shell p₀² − |p|² = m², Λ = m, and the dimensionless curvature of the K9 sweep is

κ = c⟨|q|²⟩ = ⟨|p|²⟩/m² = ⟨(γv)²⟩

Proof. The lifted score contains sinh(√c |q|), so [c][q]² = 1 by Postulate 1. The mass shell is a hyperboloid of radius m in Minkowski momentum space, of sectional curvature −1/m²; hence c = 1/m². Finally |p|/m = γv. ∎

Numbers. κ = 0.1, 0.3 and 1 correspond to typical speeds v = 0.302, 0.480 and 0.707 (from v = x/√(1 + x²), x = √κ). The curvature axis of the sweep is the axis from non-relativistic to relativistic kinematics.

Proposition 4 (the boost twin is relativistic). For a boost twin of rapidity φ, the commutator size and the acceptance bound are

||L − L−⊤|| = 2 sinh φ = 2|p|/m = 2γv, eφ = √((1 + v)/(1 − v))

the relativistic Doppler factor.

Proof. sinh φ = γv and eφ = γ(1 + v) = √((1 + v)/(1 − v)) for v = tanh φ. ∎

Numbers. φ = 0.2 gives v = 0.197 and eφ = 1.221.

Proposition 5 (a twin step is imaginary-time evolution). For any invertible linear L, gradient descent with step η in the coordinates φ = Lθ, on ℓ′(φ) = ℓ(L⁻¹φ), equals

θ ↦ θ − η (LᵀL)⁻¹ ∇ℓ(θ)

If ℓ is a variational energy E[θ] (MeV) and θ is dimensionless, then [η] = MeV⁻¹: the step is an interval of imaginary time. With LᵀL equal to the quantum geometric tensor G, the step is stochastic reconfiguration, which is imaginary-time evolution e−ηH projected onto the variational manifold (Sorella; Carleo–Troyer), up to the conventional factor 2 in the force.

Proof. ∇φℓ′ = L−⊤∇ℓ; then φ′ = φ − ηL−⊤∇ℓ and θ′ = L⁻¹φ′ = θ − η(LᵀL)⁻¹∇ℓ. The dimension follows from [θ] = [η][∇ℓ] = [η] MeV. ∎

Numbers. Excited admixtures decay as e−E₁τ after imaginary time τ = Kη, E₁ being the energy of the first excited state of the same quantum numbers. For ⁴He, E₁ = 20.21 MeV (first excited 0⁺); suppression by e−10 needs τ = 0.495 MeV⁻¹, that is 3.26 × 10⁻²² s. A twin map changes the number of steps K needed at fixed accuracy; K is the measured quantity.

Proposition 6 (closed geodesics acquire the dimension of time). In the Selberg trace formula

Σj h(rj) = (Area/4π) ∫ r h(r) tanh(πr) dr + Σγ ℓ(γ₀) g(ℓ(γ)) / (2 sinh(ℓ(γ)/2)), g(u) = (1/2π) ∫ h(r) e−iru dr

the product rℓ is dimensionless. Under r ↦ E/ε₀, a closed geodesic of length ℓγ is an orbit of period tγ = ℓγ/ε₀ and produces an oscillation of the level density with period

ΔEγ = 2π ε₀ / ℓγ

Proof. g enters through e−irℓ, so [r][ℓ] = 1 (Postulate 1). In Gutzwiller's formula an orbit of period T contributes cos(S(E)/ħ − μπ/2) with dS/dE = T, whose period in energy is 2πħ/T. ∎

Numbers. The gross-shell spacing is ħω ≈ 41 A−1/3 MeV. For Γ = PSL(2,ℤ) the shortest closed geodesic has length ℓ₁ = 2 arccosh(3/2) = 1.9248, so matching ΔE₁ = ħω fixes ε₀ = ħω ℓ₁/2π = 12.56 A−1/3 MeV; for ²⁰⁸Pb, ħω = 6.92 MeV and ε₀ = 2.12 MeV. By Huber's theorem (1959) the length spectrum and the Laplace spectrum determine each other, so the periodic-orbit periods determine the levels: the global access of the Langlands side.

Proposition 7 (the number of global parameters). For the modular surface, of area π/3, the number of Laplace eigenvalues below λ is N(λ) ~ λ/12 (Weyl–Selberg). Global spectral training with modes up to λmax therefore fits about λmax/12 coefficients per channel; 60 modes need λmax ≈ 720.

Proof. Weyl's law N(λ) ~ Area·λ/4π with Area = π/3. ∎

Numbers. A liquid-drop mass formula has 5 parameters and global mass models tens; measured masses number a few thousand. A mode basis of tens to hundreds of coefficients is overdetermined by the data, which is the condition for a stable global fit.

Proposition 8 (Hecke eigenvalues are dimensionless quantum numbers). The Hecke operators Tp, normalised by 1/√p, are dimensionless and satisfy [Tp, ΔΓ] = 0 and [Tp, Tq] = 0. Their joint eigenfunctions with ΔΓ are determined by the eigenvalues (λp) (strong multiplicity one), and |λp| ≤ p7/64 + p−7/64 (Kim–Sarnak).

Proof. An average over p + 1 points carries no units. Commutation and multiplicity one are the cited theorems. ∎

Consequence. Under ι, Hecke eigenvalues correspond to dimensionless quantum numbers, such as parity, signature and isospin projection, which label states of equal energy. They do not correspond to angular momentum, which carries ħ. The simplicity of the Laplace spectrum alone is conjectural for PSL(2,ℤ) (Cartier), so the Hecke labels play exactly the role of the extra commuting observables that make labels unique.

Proposition 9 (transfer rule from light to heavy nuclei). A move recorded as a target gap score T* is realised in a nucleus A by βA = T*√(2 ln NA)/σA.

Proof. Solve T* = βσ/√(2 ln N) for β. ∎

Numbers. With the Fermi-gas density ln ρ ≈ 2√(aE*), a ≈ A/8 MeV⁻¹ and E* = 8 MeV: ln N = 8.0 for ¹⁶O and 28.8 for ²⁰⁸Pb. At equal σ, the same T* needs a β 1.90 times larger in ²⁰⁸Pb than in ¹⁶O.

Proposition 10 (conditional: the depth law is a running exponent). If the layer index ℓ corresponds to renormalisation-group time t = ℓΔt = ln(μ₀/μ), the depth law βℓ = β₀ bαℓ becomes β(μ) = β₀(μ₀/μ)γRG with γRG = α ln b/Δt, a dimensionless anomalous dimension. The depth law itself is an unproved ansatz.

Proof. Substitute ℓ = t/Δt. ∎

When an acquired scale is new physics

Proposition 11 (classification). Every parameter of the constructs falls in exactly one of three classes under ι:

  1. a dimensionless label or indicator (T, κ, φ, λp, γRG);
  2. an algorithmic quantity, fixed by the procedure and not by nature (η, the number of steps K, the tolerance of the gate);
  3. a candidate scale of non-zero dimension (1/β, Tf, Λ, ε₀).

A parameter of class 3, or a universal value of class 1, is new physics if and only if it passes the universality test of Definition 6. Class 2 never is.

Proof. The classes partition the parameters by dimension and by whether the value is set by the data (Propositions 1–10). Postulate 3 makes universality the deciding condition. An algorithmic quantity changes with the procedure for the same data, so it cannot be universal across procedures. ∎

The h of 1900 belonged to class 3: a fitted parameter of non-zero dimension whose value proved the same across independent phenomena.

New: candidate physical scales and universal numbers

ParameterAcquired dimensionPhysical readingPredicted or known valueUniversality test
Tf = σ/√(2 ln N)MeVFreezing temperature of a finite model space0.209 σ for ²⁸Si (sd shell)Compare with Tc = 0.567 Δ across A; a fixed ratio Tf/Tc would be a new regularity
T (gap score)1√(C/2Smax), entropy-deficit fractionFreezing at C = 2SmaxSame T at the heat-capacity anomaly across nuclei measured by the Oslo method
ε₀ = ħω ℓ₁/2πMeVEnergy per unit geodesic length12.56 A−1/3 MeV for PSL(2,ℤ)One coefficient across all shell closures would be a new universal number
Λ = c−1/2MeVMass scale of a lifted collision modelm on the mass shellFitted Λ equal across reactions: a known mass (139.57 MeV, 938.27 MeV) or a new universal scale
λp1Commuting quantum numbers|λp| ≤ p7/64 + p−7/64Same labels split the same near-degenerate levels across nuclei
γRG (conditional)1Anomalous dimension of the depth lawα ln b/ΔtSame γRG across model sizes and nuclei

Known physics and algorithmic quantities

ParameterAcquired dimensionPhysical readingPredicted or known valueUniversality test
1/βMeVNuclear temperatureLevel-density temperaturesKnown physics
κ1⟨(γv)²⟩κ = 1 ↔ v = 0.707Kinematic label
φ, eφ1Rapidity; Doppler factore0.2 = 1.221Known physics
η, KMeV⁻¹; 1Imaginary-time step; steps to accuracyτ = 0.495 MeV⁻¹ for ⁴He at e−10Algorithmic: measure K, not a constant

Take-away. The constructs came largely from nuclear physics and return to it with dimensions. Four acquire non-zero dimension (1/β, Tf, Λ, ε₀), and three dimensionless quantities (T, λp, γRG) can carry universal values. Each has a stated universality test on public data. Passing any one of them would be new physics in the sense in which h was in 1900.

Map of applications: where the candidates appear

Each candidate scale of the dimensional map appears in at least two branches of nuclear science, and every branch receives at least one candidate. The map stays at the level of branches and observables.

Definition (branches). Nuclear science is partitioned into fourteen branches: eight of basic science, following the divisions of the U.S. long-range planning for nuclear physics, and six of application.

Definition (appearance). A candidate appears in a branch when the branch measures or computes the physical quantity the candidate acquires under the dictionary: a temperature or heat capacity for 1/β, T and Tf; a mass or kinematic scale for Λ; a shell or orbit period for ε₀; a discrete quantum number for λp; a running scale for γRG. The algorithmic quantity K (steps of imaginary time) appears wherever a branch solves a many-body or variational problem; it is a tool, not new physics.

BranchObservables in playCandidates that appear
B1 Nuclear structure and rare isotopesLevel densities, γ-ray strength, shell evolution far from stability, pairing1/β, T, Tf (heat capacity, pairing transition); ε₀ (shell gaps); λp (level labels)
B2 Nuclear astrophysicsReaction rates in stars and explosions, masses along the r-process path, dense matter in neutron stars1/β, T (statistical reaction rates through level densities); ε₀ (shell closures N = 82, 126); Λ (dense-matter scales)
B3 Hot and dense QCD (heavy ions)Freeze-out temperatures, collective flow, rapidity distributions1/β, Tf (freeze-out of a hadron ensemble); Λ and κ (relativistic kinematics)
B4 Hadron structure (Electron-Ion Collider)Parton distributions, three-dimensional imaging, form factorsΛ (kinematic and mass scales); γRG (scale dependence)
B5 Fundamental symmetries and neutrinosSymmetry violation in nuclei, double-beta-decay matrix elementsλp (breaking of a conserved label signals symmetry violation)
B6 Nuclear theory and computationAb initio and lattice calculations, effective field theoryγRG (running in the renormalisation group); ε₀ (semiclassical shell structure)
B7 Nuclear dataEvaluated cross sections, levels, masses, uncertainties1/β, T (level densities in evaluations); λp (spin and parity assignment); ε₀ (masses)
B8 Accelerators and detectorsEvent reconstruction, beam kinematics, spectraΛ, κ (relativistic kinematics); λp (line and level identification)
B9 Fission and nuclear energyFission paths and yields, reactor cross sections, reactor criticality, delayed neutrons and decay heatε₀ (shell effects in fragments); 1/β, T (statistical decay; level densities of fuel nuclei)
B10 Fusion energyPlasma ensembles, periodic field lines, runaway electrons, fusion-materials data, burn and reaction rates in inertial fusion energy1/β, T (plasma ensembles near critical gradients); ε₀ analogue (periods of periodic field lines); Λ = me = 0.511 MeV (runaway electrons)
B11 Isotopes for medicine and industryProduction cross sections and yields1/β, T (statistical production through level densities); λp (level labels of products)
B12 Nonproliferation and safeguards (open science)Spectral signatures, data uncertainty, attributionλp (spectral labels); 1/β (data uncertainty)
B13 Gravitation and multimessenger astrophysicsNeutron-star equation of state, gravitational-wave inference, black-hole spectraΛ (dense-matter scale); ε₀ analogue (quasinormal-mode periods in 2+1 dimensions)
B14 Radiation effects and space nuclear systemsRadiation effects in electronics and materials, the particle environment of spaceΛ, κ (relativistic particle kinematics); 1/β (statistical transport)

Proposition 12 (coverage). Every candidate appears in at least two branches, and every branch receives at least one candidate.

Proof. By inspection of the table: 1/β appears in B1, B2, B3, B7, B9–B12, B14; T in B1, B2, B7, B9–B11; Tf in B1, B3; Λ in B2–B4, B8, B10, B13, B14; ε₀ in B1, B2, B6, B7, B9, with analogues in B10 and B13; λp in B1, B5, B7, B8, B11, B12; γRG in B4, B6. Each row B1–B14 lists at least one candidate. ∎

Implication. A universal value found in one branch is immediately testable in every other branch where the same candidate appears. A freezing ratio Tf/Tc fixed in B1, for example, would predict statistical reaction rates in B2, evaluations in B7 and isotope yields in B11.

Energy: fission power, inertial fusion energy and the stars. Three energy settings draw on the same candidates. In a fission power plant, the level density read from a measured spectrum by the gap score fixes the capture and fission cross sections of fuel nuclei in the unresolved-resonance region; the transfer of Proposition 9 carries measured nuclei to the neutron-rich fragments whose decays supply the delayed neutrons used for reactor control and the decay heat after shutdown; the twin step, with η read as imaginary time, treats penetration of the fission barrier. In inertial fusion energy, as in magnetic fusion, the D–T rate weight exp(−βE − √(EG/E)) is a thermal ensemble multiplied by a tunnelling factor, with EG = 1.182 MeV and a Gamow peak at 30.9 keV for kT = 10 keV. In stars, the same weight with EG = 493 keV places hydrogen burning in the Sun at 6.0 keV; slow neutron capture at kT ≈ 30 keV runs through level densities; rapid neutron capture runs through nuclei no laboratory has yet made, the setting in which the transfer of Proposition 9 is the only route; and dense matter in neutron stars sits at κ = (3/5)(pF/m)² = 0.047 at saturation density and 0.097 at three times saturation.

Other fields that share the mathematics

The mathematics of Logarchéon is not specific to nuclear physics. Five fields share several of its tools at once; eight more share one or two.

Fields that share several tools

  1. Protein folding and biophysics. Folding theory rests on the random energy model (Bryngelson–Wolynes, 1987), the model behind CEAS. Its glass temperature is the freezing temperature of the dimensional map, and the gap score is the ratio Tg/Tphys of glass temperature to physical temperature (Proposition 13). Energy landscapes, inverse design of sequences for a target fold, and causal inference apply as well.
  2. Condensed matter on hyperbolic lattices. Hyperbolic lattices have been built in circuit quantum electrodynamics (Kollár et al., 2019); their band theory uses Fuchsian groups and automorphic functions (Maciejko–Rayan, 2021), the geometry of GRAIL. The Selberg trace formula applies, and for arithmetic lattices so do Hecke operators. Spin glasses, Ising models, Lee–Yang zeros and Anderson localisation are standard tools of the field.
  3. Quantum chaos, black holes and quantum information. Random-matrix spacing statistics and spectral form factors; the SYK model of chaotic many-body physics; the JT-gravity matrix model; modular sums and Hecke operators in AdS₃; and quantum error-correcting codes on hyperbolic surfaces, built as quotients of hyperbolic space by Fuchsian groups (Breuckmann–Terhal, 2016).
  4. Particle physics and QCD. The spectrum of the QCD Dirac operator follows chiral random-matrix theory (Shuryak–Verbaarschot, 1993); collision analysis uses Lorentz-equivariant networks; the sign problem of lattice QCD is attacked by contour deformations, which are twin maps; the renormalisation group runs throughout.
  5. Quantum chemistry. Similarity transformations before an approximate solve (transcorrelated and coupled-cluster methods), which is DL ≠ LD; neural-network wavefunctions for molecules trained with natural-gradient preconditioning; level statistics of molecular spectra.

Proposition 13 (protein folding). In the random-energy model of folding, a protein has N = eS₀ misfolded configurations with Gaussian energies of width ΔE. Its glass temperature is the freezing temperature of Proposition 2, and the gap score is the ratio of the glass temperature to the temperature:

Tg = ΔE/√(2S₀), T = βΔE/√(2S₀) = Tg/Tphys

Proof. Proposition 2 with σ = ΔE and ln N = S₀ gives the freezing temperature σ/√(2 ln N) = ΔE/√(2S₀), which is the glass temperature of the random energy model (Derrida; Bryngelson–Wolynes). Then βσ/√(2 ln N) = Tg/Tphys with Tphys = 1/β. ∎

Numbers. Sequences that fold well have a high ratio of folding temperature to glass temperature, estimated at about 1.6 for small proteins (Onuchic et al., 1995). At the folding temperature the gap score is therefore about 1/1.6 ≈ 0.63, below freezing; a gap score of 1 marks the glass transition. Here Tphys is the physical temperature, distinct from the folding temperature Tfold.

Fields that share one or two tools

FieldShared mathematics
Machine-learning theorySpin-glass loss landscapes; random-matrix spectra of trained weights; natural gradient; hyperbolic representations of hierarchies (Nickel–Kiela, 2017)
Quantitative financeCleaning of correlation matrices at the Marchenko–Pastur edge (Laloux et al., 1999); causal inference
Wireless communicationsRandom-matrix capacity of multi-antenna systems (Telatar, 1999)
Network scienceHyperbolic geometry of complex networks and greedy routing (Krioukov et al., 2010); Ramanujan graphs as optimal expanders
NeuroscienceIsing models of neural populations (Schneidman et al., 2006); criticality and avalanches; hyperbolic geometry of hippocampal representations (Zhang et al., 2023)
OptimisationTemperature control in simulated annealing; Boltzmann sampling; spin-glass formulations of combinatorial problems
CosmologyCompact hyperbolic universes and their eigenmodes; random energy landscapes
Number theory and cryptographyHecke operators, zeta functions, elliptic curves; post-quantum constructions

Take-away. Ranked by how many tools a field uses at once and how direct the translation is, protein biophysics comes first, where the CEAS freezing point is the glass temperature; hyperbolic condensed matter and quantum codes come second, where Hecke operators and the Selberg trace formula apply as they stand; quantum chaos and QCD come third. Of the candidates of the dimensional map, Tf and the gap score have their clearest second home in protein folding and spin glasses, and ε₀ and λp on hyperbolic lattices and in quantum codes.

National security within open science

The national-security value lies in open lanes that do not need a clearance:

  • uncertainty quantification of nuclear data, on which reactors and national programmes depend (the nuclear-data working group is interagency);
  • nonproliferation detection analytics and safeguards;
  • nuclear forensics attribution, where identification with abstention fits;
  • trustworthy, verified AI: gated, traceable updates.

Boundary. Open nuclear-structure and reaction data are unrestricted. Technical assistance on reactors or enrichment can fall under 10 CFR 810 and export controls; weapons-related work is classified. Work stays on open data, and a national-laboratory partner advises on any boundary.

Funding routes checked 20 October 2026

The best near-term fit is a DOE SBIR or STTR Phase I proposal, with a national laboratory as partner.

RouteSizeEligibility for LogarchéonTimingFit
DOE Office of Science SBIR/STTRPhase I $200k guideline, $250k maxSmall business: yesNormal year: proposals early DecemberNP topics are now electronics, accelerators, instrumentation; the software/data topic was dropped in FY26 (overlap with ASCR), so AI work fits ASCR or Genesis topics
Genesis Mission SBIR/STTRAbout 40 Phase I awardsSmall business: yesOpened 22 July 2026; a broader Phase I opportunity was expected in the fallFour AI-for-science topics; watch for the broader release
Genesis Mission main callAbout $294M totalTeams of two of: national lab, industry, university; for-profit cost share 20%Deadline 28 April 2026 passed; a new FY27 round planned21 challenges, including nuclear energy faster, safer, cheaper
Office of Science open call$5k–$5M per awardFor-profit allowed with a 200-word basic-science appendixFY26 window closed 30 September 2026; successor not yet confirmedNuclear Physics lists an AI/ML topic
NP AI/ML call$22M total (2024–25)See the callPast; watch for a repeatAccelerators, detectors, NP computations
Nuclear-data interagency call$100k–$1MAll applicant types; basic research2025 cycle past; watch for a repeatUncertainty quantification emphasised

Partner. An STTR research partner can be a federally funded research and development centre such as a national laboratory, not only a university.

Programmes for each branch

The direct entry for a small business is the SBIR/STTR route; other programmes take a university or national-laboratory partner. Status as of 20 October 2026.

BranchProgrammeEntry for Logarchéon
All branchesU.S. Department of Energy SBIR/STTR, including the Genesis Mission topicsSmall business (Phase I guideline $200k, maximum $250k)
All branchesGenesis Mission rounds (teams of two of: national laboratory, industry, university; 20% cost share for a for-profit); next round planned for FY27Lead or team member
B1–B5DOE Office of Science, Nuclear Physics (AI/ML topic of the open call)With a university or laboratory partner
B6DOE Advanced Scientific Computing Research (scientific machine learning, AI for science)With a partner; DOE SBIR/STTR
B7, B12Nuclear Data Interagency Working Group call (Nuclear Physics and Fusion Energy Sciences; uncertainty quantification)All applicant types
B8DOE Nuclear Physics SBIR/STTR topics (electronics, accelerators, instrumentation)Small business
B9Genesis Mission challenge on nuclear energyTeam member
B10DOE Fusion Energy Sciences (AI/ML for fusion and plasma science)With a partner; DOE SBIR/STTR
B11DOE Isotope R&D and ProductionWith a partner; DOE SBIR/STTR
B13DOE High Energy Physics (computational and theoretical research)With a partner; DOE SBIR/STTR
B3, B8, B10, B14; mathematics and AIAir Force Office of Scientific Research open BAA FA9550-26-S-0001: High-Energy Radiation-Matter Systems (nuclear and plasma physics); Computational Mathematics; Dynamical Systems and Control; physics-enhanced machine learning; Quantum Information Sciences; open until supersededIndustry and small businesses
B6, B8, B14; mathematicsArmy Research Office and Army Research Laboratory BAA W911NF-23-S-0001 for foundational research; open to 20 November 2027For-profit organizations, large and small

Take-away. Three entries cover the map: DOE SBIR/STTR, the Office of Science programmes with a partner, and the Genesis Mission rounds; the Air Force and Army open announcements add two continuously open defense entries that also fund the mathematics. The candidate with the widest reach, 1/β with T, appears in nine branches; Λ in seven and λp in six. A first result on one of them opens the most programmes.

First experiments

Each experiment is preregistered, uses public data and is compared with published baselines.

  1. Twin jump in variational training of a light nucleus (for example ⁴He with a neural-network wavefunction). One GPU. Tests DL ≠ LD directly; the strongest preliminary result for a proposal.
  2. Ψ abstention on optical-model fits to public cross sections (EXFOR).
  3. Trace-formula check of shell corrections against exact spectra.
  4. Spacing statistics on evaluated level schemes (ENSDF), with the gate deciding flags.
  5. Lorentz-exact classifier on simulated Electron-Ion Collider events, against LorentzNet and PELICAN.

Results will be posted on arXiv (nucl-th, gr-qc) with code, and used as preliminary results for a DOE Phase I proposal with a national-laboratory co-investigator.

Sources