Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Phys.

Sec. Statistical and Computational Physics

Volume 13 - 2025 | doi: 10.3389/fphy.2025.1672745

Simulating Gravitational Dynamics via Scalar Field Propagation

Provisionally accepted
  • DIRECTV LLC, El Segundo, United States

The final, formatted version of the article will be published soon.

We present a deterministic scalar-field propagation framework in which structured transport R(x) (symmetric positive-definite) and damping Lambda(x) generate reproducible media analogs of gravitational-like behavior. With the source term off (S = 0), dynamics arise purely from initial conditions under d2Phi/dt2 − div( R(x) * grad Phi ) + Lambda(x) * dPhi/dt = 0. The suite demonstrates deflection and ray bending, transit delay with escape thresholds, collapse/trapping and orbital containment, anisotropy-induced drift, repulsion via curvature inversion, and phase-sensitive interference. Each phenomenon is stated as a primary metric with an explicit falsifier and verified by energy-budget closure (Rayleigh loss + boundary flux) within 1–3% after transients, alongside spectral-safety and robustness checks. The framework yields concrete predictions: bend angle scales with grad R and flips sign under gradient reversal; transit delay grows monotonically with the path integral of Lambda and can prevent exit; bounded orbits satisfy r_min <= r(t) <= r_max over a finite capture band; the radial drift in a radial R profile follows |dr/dt| ∝ r^(−alpha) with alpha ≈ 2; transverse drift sign matches the sign of the off-diagonal/tilt in R; and interference visibility follows a cosine law in relative phase. Null and negative controls (uniform R, Lambda = 0, gradient reversal) are included. Numerically, we use stability-respecting discretizations and record discrete energies, loss, and boundary flux each step. Main figures use 256x256 grids; representative 512x512 confirmations for deflection and containment reproduce primary metrics and are included in the supplement. All code, configs, and outputs needed to regenerate figures are archived (versioned DOIs), supporting end-to-end transparency and reuse. The claims are explicitly operational—transport and loss in structured media—not statements about mass or spacetime curvature. Our contribution is twofold. First, we show that thirteen gravity-like kinematics (e.g., bending, orbital containment, transit delay) arise from a single, unified scalar transport law with spatially varying resistance and damping, d2Phi/dt2 − div( R(x) * grad Phi ) + Lambda(x) * dPhi/dt = 0. Second, we codify each phenomenon as an experiment with observables, acceptance gates, and falsification routes, and we release a complete reproducibility bundle. We emphasize these are operational analogs—not claims about spacetime curvature.

Keywords: analogue gravity, scalar-field propagation, damped anisotropic wave equation, Engineered media, reproducible computational physics

Received: 24 Jul 2025; Accepted: 26 Sep 2025.

Copyright: © 2025 Toupin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

* Correspondence: Brendan Toupin, bbtoupin@gmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.