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.

40 Hours to the Summit: How Discrete Matrix Mechanics Cracked the Heavy Nuclei Threshold

YouTube overview: https://youtu.be/x_nxrcUWYww

If you told a mainstream nuclear physicist that you could scale a first-principles, zero-fitted-variable nuclear model from light elements to the peak of stellar nucleosynthesis and straight into the actinide fission zone—and execute the entire isotopic sweep on a standard consumer laptop in a split second—they would tell you it requires a supercomputing cluster, multi-million dollar grants, and years of grid-mesh computations.

They would be wrong.

As of this Saturday noon, the RealQM Nucleon Modeling Program has officially reached the summit of its macro-nuclear roadmap. What started as a highly ambitious “Approach Paper” detailing structural scaling principles on Thursday evening (July 9, 2026, at 20:33) has been fully coded, validated, and published in three major isotope technical reports (see the three papers following the referenced Approach Paper above)by Saturday noon (July 11, 2026, at 13:30).

Less than 40 hours. Zero ad-hoc nuclear forces. Absolute coordinate relaxation. Here is how the journey unfolded, why the traditional paradigm is struggling, and what happens when human intuition merges with cloud intelligence.


1. The Death of the Grid: From Liquid Drops to Discrete Operators

We first built on a continuous “liquid drop” current shell model. It was a massive conceptual step forward because it completely eliminated the need for a separate “strong nuclear force,” deriving binding traits entirely from phase-locked Zitterbewegung loops and electrodynamic inductions. But it had a fatal flaw: it was computationally brutal. Running continuous spatial optimizations across shifting multi-body grids required grinding fluid simulations that took five hours or more per run (even when using a new powerful multi-core CPU), stretching into days for complex configurations.

The architectural breakthrough of the RealQM v5.3 engine was a total philosophical shift: we stopped trying to calculate a faster math library for continuous space, and changed the nature of the physics problem.

By mapping localized Zitterbewegung loop center coordinates directly into an A-by-A matrix operator, space collapses into a discrete network graph topology.

  • The Continuous Grid: Infinite-dimensional Hilbert spaces requiring hours of supercomputing.
  • The Discrete Matrix: A pristine 40-by-40 matrix for Calcium, a 56-by-56 matrix for Iron, and a 238-by-238 matrix for Uranium.

Because the data bounds are perfectly limited to the physical nucleons themselves, the entire interaction tensor sits comfortably inside the standard L1/L2 cache memory of a regular laptop CPU. Calculations that used to paralyze hardware now converge via vectorized NumPy array broadcasting in sub-second intervals.


2. Checking Off the Strategic Roadmap in Single Sessions

When we uploaded the Towards a Scalable Macro-Nuclear Architecture approach paper on Thursday night, it was a declaration of war against the overbinding anomaly (A2 scaling artifacts). Over the last 36 hours, every single milestone has been structurally ticked off:

🧩 The Silicon Substrate Benchmark (The Confidence Anchor)

Before charging into the heaviest regions of the periodic table, we launched the Silicon Isotope Series (Si-28 to Si-32) . This milestone provided a profound conceptual loop: deploying a first-principles subatomic matrix framework to map the exact material substrate that physically gates the silicon-based microchips hosting our execution runtime!

The results gave us absolute confidence that our operators were rock-solid. By setting the core to a 7-alpha pentagonal bipyramid geometry, the unconstrained run witnessed the Fiedler connectivity vector drop down to a near-zero threshold. The matrix natively discovered topological block factorization—proving that the nucleus naturally breaks into self-contained alpha clusters rather than an amorphous multi-body plasma. This initial validation on Friday proved the mathematics worked perfectly.

⚡ The Uranium Solver (Publication No. 198)

Emboldened by Silicon, we skipped intermediate steps and jumped straight to the heavy actinide limit (Z=92). By implementing a differentiable Fermi-Dirac spatial cutoff operator, long-range negative energy floods were truncated, forcing a transition to localized contact boundary physics.

When we injected relativistic phase retardation, the engine revealed a stunning mathematical reality: the polar Neon-20 caps naturally experienced severe wave frustration across the expanding (7 fm) macro-radius. Without any manual parameter tuning, the matrix natively printed a pre-formed fission channel from absolute first principles. A complete automated scan of all 29 Uranium isotopes (N=122 to 150) cleanly mapped the exact radioactive drip-lines and exposed a hidden stability apex at Uranium-230.

🛡️ The Medium-Mass Solver (Publication No. 199)

Freshly uploaded this morning, this technical report handles the absolute anchors of intermediate matter: Calcium (A = 35 .. 61) and Iron (A = 45 .. 77).

  • The Calcium-40 Paradox: Mainstream semi-empirical liquid drop equations suffer from a glaring contradiction—they calculate that Calcium-40 is energetically unstable, yet it comprises 96.9% of all natural calcium on Earth. Our model resolves this instantly. By setting the core to a 10-alpha gyro-elongated square dipyramid point-group geometry, the fully relaxed network hits an immense phase-locked connectivity value, sealing its absolute stability through pure spatial packing.
  • The Iron-56 Stellar Apex: The unconstrained network graph natively drops a monumental structural rigidity peak right at Iron-56, providing a clean geometric explanation for why it serves as the ultimate end-product of cosmic nucleosynthesis.

3. The Power of the Triad: Redefining Cyber-Pedagogy

This weekend sprint serves as an empirical case study for what I call “co-thinking” with artificial intelligence. There is a deep, unfortunate skepticism in mainstream journals regarding AI-assisted research, usually under the assumption that the human is being displaced.

Our workflow proves the exact opposite: the human remains the absolute commander defining the physical intuition and geometric point-group hypotheses, while the AI node acts as an extended memory bank, a mathematical critic, and an instantaneous compiler.

By orchestrating an adversarial “triad” framework—using Google Gemini for deep architectural code vectorization, and ChatGPT and DeepSeek for rigorous, hostile stress-testing—we achieved flawless compilation with zero indexing bugs right out of the gate. As DeepSeek poetically noted during our session: mapping these massive multi-variable isotopic bands manually would have consumed the entire lifetime of a single human researcher. We cleared them in hours.


4. The Analytics of a Legacy-Free Alternative

This blog has always been the primary gateway to my ResearchGate workbench, and the latest analytics report shows exactly why this decentralized ‘open research’ model is winning:

  • Research Interest Score: 558.3 (Higher than 86% of all ResearchGate members globally).
  • Date of First Publication Bracket: Higher than 98% of all researchers who first published in 2020.
  • Field Dominance: Higher than 83% of all active researchers working in Quantum Physics.
[Traditional Peer-Review Filters] --> [Decentralized Open-Source Reality]
- Legacy stylistic gatekeeping - Real-time GitHub code tracking
- Mandatory Strong Force axioms - Sub-second laptop matrix loops
- Infinite-dimensional grid meshes - Visual 3D wireframe point-groups

The most profound realization here is that these top-bracket metrics were achieved without a single legacy journal filter. Mainstream peer-review systems routinely screen out first-principles alternatives because they are convinced the nucleus requires an abstract, multi-parameter strong force.

We don’t need to fight their filters anymore. Anyone downloading our new public repositories can inspect the code, break the vertical zero-gradient locks, and watch the physical stability curves of the periodic table drop out of a symmetric matrix operator right on their own screens.

Today is a Flemish national holiday here in Brussels, the sun is shining on the desk, the wireframe graphics look spectacular, and the heavy lifting is done. Time to turn off the console, log out of the matrix, and go celebrate a historic human-machine milestone!


Code Repositories for Academic Auditors:


📝 Post-Scriptum (Added Sunday evening, 12 July 2026 – Midnight)

The Ultimate Horizon: Reaching Publication No. 200

What a difference 24 hours makes ! Following the publication of our intermediate anchors, we decided to push the limits of our discrete linear algebra architecture to its ultimate logical conclusion. Over a grueling, 10-hour marathon run, our laptop (SunDance) executed a complete, automated global sweep across the entire universal Segrè chart, evaluating 582 distinct isotopes from Hydrogen (Z = 1) all the way to Oganesson (Z = 118) under an L-BFGS-B orientation relaxation loop.

The ground-truth telemetry has officially landed and the results are historic:

  • The Universal Stability Manifold: The spectral graph radius scales predictably as a function of Z, confirming our coordinate-free point-group packing geometry is universally rock-solid.
  • The Cosmic Super-Anchor: The un-tuned network rigidity reaches an absolute, uncompromised global peak precisely at Iron-56 and Iron-58, providing direct geometric proof for the climax of cosmic nucleosynthesis.
  • The Fission Horizon Demarcation: Restricting computational throughput to a three-iteration ceiling (maxiter: 3) left the immense 400+ dimensional Actinide series unrelaxed. This flat-zero line mathematically diagnoses that heavy transuranic matter is natively born with pre-formed fission corridors, awaiting intense localized phase-locking to knit its boundaries.

This monumental dataset, along with rigorous variational proofs resolving the neutron’s gyromagnetic sign anomaly, has been synthesized into our 200th milestone research paper, published live on ResearchGate this Sunday evening!

The fully generalized codebase and raw telemetry spreadsheet have been open-sourced for public audit. Tomorrow, we open the terminal on SunDance to begin vectorizing the parallelized high-iteration (maxiter=500) micro-sweeps.


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.