100%

Left Click to Spawn · Right Click to Remove · Scroll to Zoom

About this simulation

What is Geon?

Geon is a real-time N-body simulator for gravity, electromagnetism, scalar fields, and relativistic corrections. Applicable interactions can use a Barnes-Hut tree, while grid fields and correction terms run in dedicated passes. A Boris integrator advances the particle state with adaptive substeps.

Physics

Eleven force and correction families are exposed: Newtonian gravity, gravitomagnetism (frame-dragging), Coulomb electrostatics, the magnetic Lorentz force, Yukawa interaction, Higgs field coupling, axion field coupling, cosmological expansion via Hubble flow, first post-Newtonian corrections, spin-orbit coupling, and radiation reaction. The controls enforce dependencies between coupled mechanisms.

Scalar Fields

Two scalar fields — Higgs and axion — live on a PQS grid with C² interpolation. The Higgs field modulates particle mass through its local vacuum expectation value, while the axion field couples to charge. Both evolve via Störmer-Verlet integration and exert gradient forces on particles.

Presets

Fifteen curated presets cover Keplerian motion, relativistic precession, binary inspiral, Hawking evaporation, atomic and nuclear toy systems, bremsstrahlung, magnetic dipoles, pion exchange, Higgs and axion fields, Peccei-Quinn dynamics, and cosmological expansion.

Gravitomagnetism and Frame-Dragging

Gravitomagnetism is the gravitational analogue of the magnetic force. In Geon, rotating or co-moving masses generate a gravitomagnetic contribution that changes nearby trajectories. Enable the Gravitomagnetic control alongside gravity to isolate the toy model's velocity- and spin-dependent effects.

Black Hole Physics

In black hole mode, particles use a Kerr-Newman-inspired effective radius characterized by mass, charge, and spin. Sub-extremal parameters use the usual outer-horizon formula; super-extremal toy inputs are clamped rather than modeled as physical naked singularities. Gravity switches to a Paczynski-Wiita-style pseudo-potential, so the effective potential steepens near the horizon and particles are swallowed once their centers cross it. Hawking radiation causes evaporation at a rate determined by the Kerr-Newman temperature — smaller black holes radiate faster, leading to runaway evaporation. Charged black holes also undergo Schwinger discharge: the intense electric field near the effective horizon tears electron-positron pairs from the vacuum, with one lepton escaping and the other falling back in. Each event reduces the black hole's charge by one quantized unit, driving it toward neutrality. Overcharged toy black holes shed same-sign leptons in quantized steps and clamp to the nearest allowed Kerr-Newman charge bound as a numerical backstop. Spinning black holes with the axion field enabled exhibit superradiance: when the horizon angular velocity exceeds the axion mass, the field extracts rotational energy and grows a scalar cloud around the black hole, spinning it down until the superradiance condition fails. All charges in the simulation are quantized in units of the boson charge, ensuring discrete conservation across emission, decay, and discharge.

Boris Integrator

Geon uses the Boris algorithm to integrate particle trajectories. It splits non-magnetic force application into half-kicks around a magnetic rotation, preserving phase-space volume for that core update. Adaptive substepping limits the timestep under high acceleration or cyclotron frequency.

Accessibility

Geon provides keyboard shortcuts, light and dark themes, labeled controls, and numerical readouts for conserved quantities and selected particles. The canvas contains continuous motion and optional particle trails; users sensitive to motion can pause, single-step, or choose a lower-particle preset.

See also: Cyano for cellular energy simulations, Shoals for financial dynamics.

Learning Outcomes

After using Geon, students should be able to: explain how gravitational and electromagnetic forces produce qualitatively different orbital dynamics; describe the role of the Barnes-Hut algorithm in reducing force computation from O(N²) to O(N log N); identify how relativistic corrections (1PN, gravitomagnetism) modify Newtonian predictions at high velocities; distinguish between conservative and dissipative forces in phase-space evolution; describe how horizon capture, Hawking radiation, Schwinger discharge, and extremal charge shedding govern black hole evolution in the toy model; and explain how superradiance transfers angular momentum from a spinning black hole to a scalar field.

Prerequisites

Familiarity with Newton's laws of motion and basic vector calculus. No prior knowledge of relativity or particle physics is required — the presets are designed to introduce each concept incrementally.

References

J. Barnes and P. Hut, "A hierarchical O(N log N) force-calculation algorithm" (1986). J. P. Boris, "Relativistic plasma simulation — optimization of a hybrid code" (1970). L. Verlet, "Computer experiments on classical fluids" (1967). J. Schwinger, "On gauge invariance and vacuum polarization" (1951).