How CoPilot Independently Rebuilt My Uranium Solver

Over the past weeks, I’ve been developing RealQM Isotope Solvers — a discrete, matrix‑mechanics alternative to the traditional continuous nuclear models. It’s a framework that treats the nucleus not as a quantum fluid, but as a structured electrodynamic network: alpha blocks, neutron buffers, polar caps, and fringe satellites, all interacting through real-valued, phase‑retarded matrix paths.

It’s a model that scales from light nuclides like Silicon‑28 all the way up to the heavy actinides. And it’s open-source: the full code lives on GitHub, ready for anyone to run.

But here’s the thing: DeepSeek and Gemini had co-developed the code and, therefore, understood both the model as well as how to work with the code. And then I tried to get ChatGPT to explore the physics, math and the code, and it just could not do that. So I tried CoPilot: I just gave it the physics paper and pointed to the GitHub repository, and I then asked it to re-build the code to run the full Uranium isotopic chain (N = 122 to 150), extract the Laplacian spectrum and give me the Fiedler curve. It executed the job flawlessly.

Moreover, it clearly understood what it was doing as it noticed the Fiedler curve did – unlike the Gemini code – not yield the expected stability peak at Uranium-230 and it immediately knew why: because it used a linear neutron‑allocation rule, while the RealQM v5.3 Solver uses a nonlinear buffer/fringe transition tuned to the Lead‑208 core geometry. Here are the two curves:

This wasn’t “talk about RealQM.” This was RealQM in action. As mentioned above, the CoPilot graph is not identical to the Gemini/DeepSeek curve — because Copilot used a linear neutron‑allocation rule — but that doesn’t matter all that much in terms of testing the model: the important part is that the solver reacted to geometry and produced a real connectivity landscape: it shows the expected oscillations, pre-apex stiffening and then the high-mass degradation. And CoPilot knew it. In other words, it thought about what is was doing.

Why This Matters

The RealQM UraniumSolver is not a toy model. It’s a heavy nuclear solver incorporating:

  • 476 angular degrees of freedom
  • 238×238 interaction matrices
  • Fermi-Dirac horizon windows
  • relativistic retardation operators
  • multi-core alpha clustering
  • neutron buffer sheets
  • polar Neon‑20 caps
  • fringe valence satellites

It’s the kind of physics that normally requires HPC clusters. Yet three AI systems (Gemini, DeepSeek and now CoPilot) entered the codebase, understood the architecture, and executed the full isotopic chain on my laptop.

I now think there are two types of LLMs: LLMs that talk, and LLMs that think. That’s the dividing line.

LLMs that talk:

  • summarize
  • paraphrase
  • explain
  • rephrase
  • stay at the surface

LLMs that think:

  • internalize structure
  • generalize logic
  • manipulate code
  • execute physics
  • produce new scientific output

The RealQM test exposes the difference instantly. Give them the paper and the code. See what they do. The ones that think take the code and run with it. The others don’t.

What Comes Next

Now that Copilot has independently validated the solver, we can:

  • refine the geometry rules
  • implement the nonlinear buffer/fringe transition
  • reproduce the exact U‑230 apex
  • map the full actinide series
  • extend RealQM to Thorium, Plutonium, and beyond

And we can write about it — not as another ResearchGate paper, but as a story of how open physics and modern AI can work together.

Closing Thoughts

This wasn’t just a computation. It was a demonstration. A demonstration that RealQM is reproducible. A demonstration that AI systems can genuinely think. A demonstration that open-source physics can scale.

And a demonstration that the future of scientific computing may not be about bigger clusters — but about smarter tools.

Revisiting the Neutron and Deuteron puzzle

My previous note on the proton model utilized radically simplified semi-classical reasoning to recover empirical metrics without introducing free parameters.

This new paper scales that exact framework up into the multi-body nuclear domain, treating the neutron and deuteron not as static configurations bound by unobservable “glue” forces, but as an elegant, non-linear synchronization problem involving coupled electromagnetic phase clocks.

Oddly enough, by shifting the ontology away from isolated particles toward relational, phase-locked coherence, the math naturally operates within realistic nuclear regimes—generating an internal neutron magnetic radius of 0.81-0.93 fm, a finite spatial interaction boundary of about 2 fm, and a near-field locking energy of about 2 MeV. These values all closely match experimentally observed ranges.

We, therefore, think this is quite significant. If anything, it shows, perhaps, that progress sometimes does not come from adding more parameters to describe some ‘black box’, but from acknowledging that stable matter may correspond to highly constrained, coherent oscillatory organization.

Read the paper here: “Relational Stability and Synchronization Geometry in the Neutron–Deuteron System

Post Scriptum (23 May 2026):
A subsequent multi‑stage sanity check, involving adversarial cross‑checking between DeepSeek, ChatGPT, and Google Gemini, resulted in three companion pieces that should be read alongside the main paper (click on ‘public files’ on the above‑referenced RG page).

  1. “On the Factor 2 in the Electron’s Ring‑Current Model: A Clarification of Scales” resolves a long‑standing confusion about the electron’s Compton radius and the equipartition of energy, showing that the model is internally consistent.
  2. “On the Binding Energy of the Deuteron: A Correction and Reinterpretation” corrects a numerical error in the static magnetic dipole‑dipole calculation (the correct value is ~15 keV, not 2.2 MeV) and reinterprets the deuteron binding energy as a non‑linear phase‑locking energy.
  3. “The Fine‑Structure Constant and the Deuteron Binding Energy” (with an appended sanity check by Gemini) completes the arc: from the heuristic proposal ηα(mpc2/2)2.31ηα⋅(mpc2/2)≈2.31 MeV (4% error) to the logically and numerically superior expression (1η)αmpc22.22(1−η)⋅αmpc2≈2.22 MeV (error <0.3%), using only the incoherent neutron deficit (1η)(1−η) and the full proton rest energy. The fine‑structure constant αα enters naturally as the electromagnetic coupling strength.

All notes are available on the ResearchGate page. I thank DeepSeek for its careful analytical assistance and for helping to turn an initial overreach into a refined, honest, and testable hypothesis.