← Gallery

Muonic hydrogen & muon-catalyzed fusion — what RealQM can do

A muon (mμ ≈ 207 me) bound to one or two nuclei, in RealQM. Reduced models; the robust result is geometric (the ~150× compression), not the absolute energetics.
▶ Launch ddμ scan D2+ binding curve, rescaled to the muon · WebGPU/Chrome

1. The key fact: point-nucleus muonic H is the electron problem rescaled

A muon is just a heavy lepton: mμ ≈ 207 me. The bound-state Schrödinger equation with a point Coulomb potential is scale-invariant in the lepton mass — substituting r → r/m maps the muon problem exactly onto the electron problem. So for muonic hydrogen the length shrinks by ~1/m and the energy grows by ~m:

This means that if RealQM reproduces electronic hydrogen (it does), it reproduces muonic hydrogen structure by construction — the same dimensionless grid solve, with the axes relabelled. That is a validation, not a discovery. The interesting muon-specific physics lives where this scaling breaks: finite nuclear size, and molecular ions where the nuclei are not infinitely heavier than the muon. This page tackles the second.

2. Muon-catalyzed fusion: the muonic molecular ion ddμ

Put one muon between two deuterons and you have a hydrogen molecular ion D2+ with the electron replaced by the muon. Because the muon is ~207× heavier, the ion is ~207× smaller — the equilibrium internuclear separation drops from ~0.74 Å (the D2 bond) to a few hundred femtometres. That compression pulls the deuterons into quantum-tunnelling range, and they fuse: this is muon-catalyzed (cold) fusion, observed since the 1950s. The muon is released and catalyses again (~102 fusions per muon for dtμ before it decays or sticks).

What was performed

RealQM (molecule.js, real-space grid, no basis set, no fitted force field) computes the D2+ electronic binding curve: one electron shared by two bare deuterons (Z = 0, Znuc = +1 each), nuclei clamped, scanning the internuclear distance R and relaxing the electron density to convergence → E(R). The muonic ion then follows by exact Born–Oppenheimer mass-scaling (clamped nuclei → use mμ, not the reduced mass):

Rμ = Re / 206.8    Eμ = Ee × 206.8

Result

R (a₀)E (Ha)→ Rμ (fm)
1.4−0.505358
1.7−0.525435
2.0−0.538512
2.4−0.531614
3.0−0.489767

The binding curve has its minimum exactly at Re = 2.0 a₀ — the textbook H2+ equilibrium — rising symmetrically on both sides. Rescaled by the muon mass:

So RealQM, handed only the nuclei, reproduces the geometric heart of muon-catalyzed fusion: the muon screens the two deuterons down to ~500 fm, into tunnelling range.

3. Scope — what this does and does not show

4. Relation to exact / first-principles methods

Muonic molecular ions are a classic three-body Coulomb problem, solved to extreme precision decades ago by variational and adiabatic methods (e.g. Hylleraas-type and complex-coordinate calculations), which give ddμ / dtμ binding levels to meV and the weakly-bound resonance level (Vesman mechanism) that underpins muCF rate theory. Those are the gold standard for the energetics.

 Variational 3-body / adiabatic (state of the art)RealQM (this run)
EnergiesmeV-accurate binding, all rovibrational levelsgeometric Re robust; absolute energy grid-limited
Nuclear motionfull non-adiabatic three-bodyBorn–Oppenheimer (clamped)
Methodtailored basis (Hylleraas / Gaussian), problem-specificreal-space grid, no basis set, parameter-free
What it showsthe numbers that feed muCF rate theorythe ~150× compression mechanism, from scratch

RealQM does not compete with variational three-body methods on the muCF energy levels. Its contribution here is conceptual: a single, parameter-free real-space formalism reproduces electronic chemistry and — by nothing more than the lepton mass — the muonic compression that drives cold fusion, making the continuity between the two regimes explicit.

Bottom line. Point-nucleus muonic hydrogen is electronic hydrogen rescaled, so RealQM gets its structure for free. For the muonic molecular ion ddμ, RealQM computes the D2+ binding curve and the muon mass rescales it to Rμ ≈ 512 fm — a ~144× compression of the deuterons into tunnelling range, the mechanism of muon-catalyzed fusion, from first principles. The energetics are grid-limited and BO; the geometry is the real result.

← Back to Gallery