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.
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).
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
| R (a₀) | E (Ha) | → Rμ (fm) |
|---|---|---|
| 1.4 | −0.505 | 358 |
| 1.7 | −0.525 | 435 |
| 2.0 | −0.538 | 512 |
| 2.4 | −0.531 | 614 |
| 3.0 | −0.489 | 767 |
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.
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) | |
|---|---|---|
| Energies | meV-accurate binding, all rovibrational levels | geometric Re robust; absolute energy grid-limited |
| Nuclear motion | full non-adiabatic three-body | Born–Oppenheimer (clamped) |
| Method | tailored basis (Hylleraas / Gaussian), problem-specific | real-space grid, no basis set, parameter-free |
| What it shows | the numbers that feed muCF rate theory | the ~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.