The ablation of the bounded factor deserves particular attention because it isolates the one genuinely novel ingredient of the proposal. Removing the bounded factor and retaining only the acceleration-dependent response degrades the mean velocity residual by roughly half a kilometre per second, a small but non-zero margin whose bootstrap interval excludes zero. In galaxy-scale fitting, where typical observational uncertainties and mass-to-light systematics are of comparable or larger size, an improvement of this magnitude is best described as a measurable but modest contribution. The authors are explicit that this improvement cannot be taken as evidence for the coherence hypothesis itself, because the factor is constructed endogenously from the model-predicted velocities rather than from any independently measured property of the stellar or gas motions. This is a subtle but important epistemic point: a term that is derived from the model’s own output will, by construction, correlate with the model’s successes, and only an external, independently measured variable can break that circularity. The proposed resolution, testing whether independently measured phase-space order predicts the response on held-out systems, is precisely the kind of falsifiable follow-up that separates a motivating hypothesis from a fitted fudge factor.

The cluster result is arguably the most instructive part of the paper precisely because it is a negative one. Strong gravitational lensing probes the projected mass distribution of clusters at radii and in regimes where disk-galaxy intuition does not automatically transfer, and the baryonic mass budget of a cluster is dominated by hot intracluster gas rather than stars, which changes the character of any baryon-based prescription. Calibrating a single effective amplitude on forty-two strong-lensing systems and then applying it without refitting to seventeen CLASH clusters yields a median predicted-to-observed ratio well above unity, meaning the prescription systematically underpredicts the lensing mass on average, and the discrepancy varies with radius rather than behaving as a simple overall offset. A radius-dependent bias of this kind cannot be repaired by rescaling the amplitude; it indicates that the functional form itself, including the assumption of a single nominal path length for every cluster, is missing physics that differs from the disk case. The authors’ decision to present this as a falsifiable limitation rather than to introduce cluster-specific parameters is methodologically commendable, because adding per-cluster freedom would have converted a genuine transfer failure into an apparent success while telling us nothing new about nature.

The auxiliary QUMOND construction provides a theoretical bridge that clarifies what kind of underlying structure the empirical response could have. QUMOND, the quasi-linear formulation of modified Newtonian dynamics, replaces the Poisson equation with a modified one in which the Newtonian potential feeds back through an algebraic function of its own gradient magnitude. The finding that an action of this type can reproduce the Σ-gravity response when the parameter B is prescribed as a spatially constant quantity, but not when it is generated endogenously from the galaxy model, tells us where the theoretical difficulty lies. Constant-B limits are the easy case for any modified-Poisson framework, since they reduce to a fixed algebraic modification of the field equation. A spatially and temporally varying B that depends on the dynamical state of the baryons is a far stronger requirement, and no existing covariant formulation supplies it without additional structure. This gap is not a fatal objection at the phenomenological stage, but it defines the boundary between what has been demonstrated and what remains aspirational.

It is also worth situating the acceleration scale used in the response function. The quantity g† is built from the speed of light and the Hubble constant, echoing a long-standing observation that the characteristic acceleration below which disk galaxies deviate from Newtonian expectations is numerically close to cH0 divided by a factor of order ten. Any successful theory of this regime must ultimately explain why a cosmological parameter appears in galaxy-scale dynamics, whether through a genuine physical coupling or through a coincidence of structure formation. The baryonic Tully–Fisher relation, which connects the flat asymptotic speed of a disk to the fourth root of its baryonic mass with remarkably small scatter, constrains any candidate response function in its deep-acceleration limit, and the square-root structure of the transition function in Σ-gravity is chosen to respect that constraint. The authors acknowledge that the precise shape of the transition between high and low acceleration is a design choice rather than a derivation, which is the correct posture for an empirical proposal at this stage.

For readers weighing these results against the broader landscape, several contextual points are useful. First, the standard cosmological model rests on evidence far beyond galaxy rotation curves, including the acoustic peak structure of the cosmic microwave background, the growth rate of large-scale structure, and the inferred mass-to-light ratios of galaxy clusters through independent channels. Any alternative phenomenology that aspires to more than galaxy-scale bookkeeping must eventually engage that full body of evidence, and the authors correctly decline to claim otherwise. Second, the residual mass discrepancy in galaxy clusters is a known and persistent challenge for MOND as well, typically requiring either additional unseen matter or cluster-specific modifications, so the transfer failure documented here is consistent with a difficulty that predates this proposal. Third, the exclusion of bulge-dominated radial measurements from the SPARC analysis, while defensible on the grounds that the disk prescription is not expected to describe bulge kinematics, narrows the tested domain in ways that should be kept in mind when comparing headline statistics across the literature.

The matched counter-rotation diagnostic and the MaNGA integral-field comparison serve as secondary checks that broaden the empirical footprint beyond the SPARC sample, and the structural sensitivity test on the assumed photometric scale lengths probes how much of the agreement depends on the input baryon models rather than the response function itself. That the scale-length test does not outperform simple mass-to-light assignment controls is a sobering reminder that in disk-galaxy phenomenology, uncertainties in the stellar population synthesis inputs can rival the differences between competing prescriptions. This is not a reason for despair but a calibration of expectations: distinguishing between closely related algebraic responses requires either tighter photometric constraints or samples selected to maximize the discriminating power of the comparison.

Looking forward, the most consequential test the paper identifies is whether an independently measured coherence variable, drawn from the observed phase-space structure of stellar and gas motions, can predict the response amplitude on galaxies not used in any calibration. If such a variable exists and carries predictive power beyond the acceleration alone, the coherence hypothesis graduates from motivation to measurement. If it does not, the bounded factor joins the long list of empirically motivated corrections that fit well but explain little. Either outcome would be informative, and the framing of the paper ensures that the community can adjudicate it with data rather than rhetoric. In the meantime, the study stands as a careful, self-limiting contribution: a fixed-form galaxy response whose aggregate performance is statistically indistinguishable from the established MOND prescription under identical assumptions, a cluster amplitude that demonstrably fails to transfer, and a clearly marked boundary between what has been tested and what remains hypothesis.

Subject of Research: Σ-Gravity: a coherence-motivated empirical response tested in galaxies and clusters

Article Title: Σ-Gravity: a coherence-motivated empirical response tested in galaxies and clusters

Article References: Speiser, L. (2026). Σ-Gravity: a coherence-motivated empirical response tested in galaxies and clusters. Frontiers in Astronomy and Space Sciences, 13, Article 1866133. https://doi.org/10.3389/fspas.2026.1866133

Image Credits: AI Generated

DOI: 10.3389/fspas.2026.1866133

Keywords: Gravity, coherence-motivated, empirical, response, tested, galaxies, clusters, scientific research

The framework, developed by Gustav Stålhammar of St. Erik Eye Hospital and the Karolinska Institutet in Stockholm, rests on one deceptively simple hypothesis. Ordinary matter is modeled as a collection of localized emitters, each radiating waves of an additional scalar field with random, mutually uncorrelated phases. In neutral matter, a bound atom acts as one emitter; in fully ionized plasma, each free charged particle counts as an independent source. No individual emitter is observable, and the frequencies involved are arbitrary. What matters is only the statistics of the summed field, which by the central limit theorem behaves as a Gaussian random field regardless of the underlying spectrum.

When the waves from roughly 10^57 emitters in a star like the Sun are added together at any point in space and the intensity of the total field is averaged over time, three kinds of terms appear, and each maps onto one aspect of gravitation. The first is the set of self-terms: each wave multiplied by itself, always positive, adding incoherently like the masses of the sources. The framework identifies ordinary Newtonian gravity with this stationary part, taking Newton’s constant from measurement and leaving the tensor sector of gravity exactly that of general relativity. The hypothesized waves are not gravitational waves; they belong to an additional scalar degree of freedom of the kind long studied in scalar-tensor theories.

The second and third sectors arise from the cross-terms in the double sum. Waves arriving from all matter within the causal horizon combine into a fluctuating background whose variance sets a characteristic acceleration scale. Using the cosmic mean baryon density multiplied by the Hubble length, and a geometric factor of 3 over 8π, the framework derives an acceleration scale g_w of about 7.9 × 10^-11 meters per second squared, fixed entirely by measured cosmological parameters rather than fitted. This is tantalizingly close to the acceleration scale near cH_0, the speed of light times the Hubble constant, at which galaxy dynamics are known to depart from Newtonian expectations, a coincidence that modified Newtonian dynamics, or MOND, has long treated as fundamental but never explained.

The crucial new ingredient is the interference between the locally generated field and this cosmic background. The signed average of the cross-term is zero, so nothing would happen if the medium responded to the instantaneous signed product. The framework’s central assumption, labeled A2, is that the material medium responds instead to the product of the two components’ envelope amplitudes, their absolute magnitudes, independent of phase. Under this rectified response, the mean of the product for two independent Gaussian envelopes is exactly π/2 times the product of their standard deviations, a result verified by Monte Carlo simulation to within 0.12 percent. The result is an additional attractive acceleration proportional to the geometric mean of the Newtonian acceleration and the cosmic scale, directed along the local field, precisely the square-root behavior that characterizes dark-matter phenomenology in the weak-field regime.

Three phenomenological functions modulate this term. A screening factor suppresses the effect where a system’s own column density exceeds the cosmic source column, which is what keeps the Solar System clean. A transport saturation factor describes how the rectified response accumulates over a density-dependent length anchored at the local Galactic disk value. And an ionized-gas enhancement reflects the emitter-counting rule: ionization multiplies the number of independent emitters per unit mass, so plasma environments should feel a stronger effect. Three constants are fixed in advance from cosmology and local Galactic measurements; only one coefficient, the plasma coefficient γ_p, is calibrated on data, using the X-COP sample of twelve massive clusters, and then frozen for every other test.

The performance on galaxy rotation curves is the headline result. Applied with no galaxy-fitted parameters whatsoever to 3,187 weak-field points in 175 galaxies from the SPARC database, the framework achieves a root-mean-square error of 0.228 dex in log-acceleration, with a 95 percent confidence interval of 0.204 to 0.253 dex. For comparison, general relativity with baryons alone yields 0.398 dex, MOND at its canonical acceleration constant yields 0.270 dex, and MOND with the constant re-fitted to the same selection yields 0.213 dex. The framework thus outperforms canonical MOND and baryons-only gravity while landing between the two MOND variants, all without a single parameter tuned to galaxies. Notably, its prediction shifts by only 0.005 dex when the stellar mass-to-light assumptions are changed, whereas the fitted MOND constant moves by a factor of 2.4 under the same variation.

On the cluster side, the story is more mixed but structurally interesting. Leave-one-cluster-out validation on X-COP yields an error of just 0.072 dex, and the calibrated plasma coefficient of 0.69 lies within a factor of 1.6 of an independent, parameter-free ionization estimate of 1.08. What distinguishes the framework from MOND here is not fit quality but internal consistency: MOND requires an acceleration constant 5.5 times larger to fit clusters than to fit galaxies, whereas the same kernel with the same fixed constants spans both regimes. However, applied unchanged to 22 clusters from the eROSITA sample, the framework roughly halves the central shortfall of baryons alone but retains a radially structured under-prediction of 0.435 dex overall, comparable in magnitude to MOND’s own central failure. Four of five dwarf spheroidal galaxies are also under-predicted, bracketing the model’s domain of validity at low baryonic density.

The framework survives the obvious local tests through its own screening structure. Predicted fractional anomalies in the Solar System range from 8.5 × 10^-14 at Mercury to 4.0 × 10^-8 at Neptune, between three and more than eight orders of magnitude below the Cassini sensitivity limit, with a steep radial increase that future ephemeris measurements at Neptune-to-Kuiper distances could probe directly. Wide-binary predictions depend on a stated transport convention for point-pair systems; under the ambient-disk reading, the predicted velocity excesses remain below current Gaia sensitivity, but measurements beyond 100 kAU would discriminate sharply between the two readings. As a conditional cosmological check, a matter-like scalar background reproduces the baryon-acoustic-oscillation sound horizon to within 0.04 percent under a matched drag-epoch prescription, though the full cosmological sector remains an open problem.

The author is explicit about what remains unproven. The rectified response of Assumption A2 has no microscopic derivation, and the Hubble-length causal column is an imposed simplification. The screening and transport functions are phenomenological, with demonstrated functional-form freedom at the 0.01 to 0.06 dex level. Yet the framework makes concrete, falsifiable predictions: the plasma coefficient must re-calibrate to the same value on independent cluster samples, ultra-wide-binary velocities beyond 100 kAU must match the stated curve, matched-radius lensing and dynamical masses must agree, and the scalar background’s energy density must evolve like non-relativistic matter through the drag epoch. All numerical results were produced with a single released analysis pipeline, making every modeling choice auditable. Whether the universe really generates its missing gravity from the interference statistics of ordinary matter is now a question that measurement, not philosophy, will decide.

Subject of Research: An emergent scalar-field model that derives Newtonian gravity and dark-matter phenomenology from random-phase waves emitted by ordinary matter

Article Title: Newtonian gravity and dark-matter phenomenology from an emergent scalar field sourced by ordinary matter

Article References: Stålhammar, G. (2026). Newtonian gravity and dark-matter phenomenology from an emergent scalar field sourced by ordinary matter. Frontiers in Astronomy and Space Sciences, 13, Article 1904894. https://doi.org/10.3389/fspas.2026.1904894

Image Credits: AI Generated

DOI: 10.3389/fspas.2026.1904894

Keywords: dark matter, scalar field, modified gravity, MOND, galaxy rotation curves, galaxy clusters, cosmology, SPARC, X-COP, baryonic matter, gravitation, emergent gravity

The study, conducted by Francisco A. Guzmán, Elías Castellanos, and Jorge Mastache, tackles the Gross–Pitaevskii–Poisson (GPP) system, a pair of coupled equations that describes a Bose–Einstein condensate evolving under its own gravity. In this picture, the entire dark matter halo of a galaxy is represented by a single macroscopic wavefunction whose squared amplitude gives the mass density. The first equation is a nonlinear Schrödinger equation that includes the gravitational potential and a short-range repulsive interaction between the bosons, while the second is the Newtonian Poisson equation, which computes the gravitational field generated by the condensate’s own density. Solving these equations simultaneously is notoriously difficult because the problem is nonlinear and nonlocal: the gravitational potential depends on the wavefunction everywhere, not just at a single point.

What sets this analysis apart from much of the existing literature is a deliberate choice of bookkeeping. Many previous studies exploit mathematical scaling symmetries of the equations, solving everything in dimensionless units and only converting back to physical quantities at the end. The Mexican team instead kept three quantities explicit throughout: the boson mass, the scattering length that characterizes the strength of repulsive self-interaction between particles, and the total particle number N. That last parameter is crucial, because it directly fixes the total halo mass through the simple relation that the halo mass equals the boson mass multiplied by the particle count. By treating N as an independent control parameter, the researchers built a direct bridge between microscopic particle physics and macroscopic observables such as halo size, density profile, and rotation speed.

Technically, the team discretized the radial equations on a finite grid extending to 100 kiloparsecs and solved the resulting nonlinear eigenvalue problem with a Newton–Raphson iterative scheme, updating the wavefunction, the gravitational potential, and the chemical potential self-consistently until convergence to tolerances between 10 and 8 and 10 to the minus 10. To guard against artifacts of the starting guess, they initialized the solver with four different smooth, localized profiles: Gaussian, exponential, linear-exponential, and hyperbolic-secant shapes. A striking robustness emerged. Whenever different starting profiles converged to the same solution branch, the resulting stationary configurations were numerically indistinguishable, with relative differences below roughly one part in ten thousand for the wavefunction. The choice of initial guess affected only how often the iteration converged, not the physics of the answer.

Scanning a 40-by-40 grid of boson masses between 10 to the minus 23 and 10 to the minus 20 electronvolts and particle numbers between 10 to the 93 and 10 to the 99, the researchers classified every converged solution into three families. Ground states are nodeless, meaning the wavefunction never crosses zero, and they represent the lowest-energy equilibrium of the condensate. Excited states possess one or more radial nodes, producing oscillatory density shells around a central core, much like the excited orbitals of an atom. Unbound configurations fail a binding criterion, meaning the energy per particle exceeds the depth of the self-generated gravitational well, so no localized condensate can exist. The resulting phase diagrams show the ground-state branch occupying a well-defined diagonal band in the mass–particle-number plane: heavier bosons need fewer particles to bind, while lighter ones require vastly larger numbers.

One of the most consequential findings concerns how repulsive self-interactions reshape this map. As the scattering length increases, the ground-state branch systematically migrates toward larger boson masses and smaller particle numbers, because the repulsive pressure helps support the condensate against gravitational collapse. The excited-state and unbound regions, by contrast, shift far more slowly. This sensitivity matters for phenomenology: the location of the ground-state band, the region where realistic galactic halos live, is exquisitely dependent on the strength of the self-interaction, while the rest of the solution space is comparatively stubborn.

From the converged solutions, the team extracted empirical scaling laws for the characteristic halo radius, defined as the sphere enclosing 99 percent of the mass. In the weakly interacting regime, the fitted relation gives a radius proportional to the boson mass raised to the power of minus 2.996 and the particle number to the power of minus 0.998, values astonishingly close to the analytic prediction of minus 3 and minus 1 for non-interacting solitons. Rewritten in terms of halo mass, this recovers the famous mass–radius relation of Schrödinger–Poisson solitons, in which the radius scales inversely with the boson mass squared and the halo mass, providing a clean consistency check of the numerics. When self-interactions dominate, the scaling shifts toward the Thomas–Fermi expectation, with the radius growing as the square root of the scattering length and falling as the boson mass to the power of minus 1.545. Yet the fitted exponents reveal that the solutions occupy an intermediate regime where quantum pressure, self-interaction pressure, and gravity all contribute simultaneously, deviating measurably from the asymptotic Thomas–Fermi limit. A three-parameter fit confirms a clear hierarchy: the boson mass exerts the strongest influence on halo size, the particle number a moderate one, and the scattering length a softer but non-negligible effect.

The astrophysical payoff comes when these equilibrium configurations are confronted with real galaxies. Using ground-state solutions with a boson mass of 10 to the minus 22 electronvolts, the team modeled the rotation curves of two dwarf galaxies from the SPARC catalog, KK98-251 and UGC01281, adding the baryonic contribution from gas and stars to the condensate’s contribution computed from the enclosed mass profile. The results were impressive. For KK98-251, a halo of about 896 million solar masses with a characteristic radius of 5.05 kiloparsecs reproduced the observed kinematics, while UGC01281 required roughly 3.58 billion solar masses spread over 5.88 kiloparsecs. Remarkably, in both cases the solitonic core alone accounted for the dark matter, with no need for the extended interference-supported envelope that cosmological simulations of fuzzy dark matter typically predict around the central condensate.

Here, however, the story collides with cosmological constraints. Observations of the Lyman-alpha forest, the thicket of absorption features imprinted by intergalactic hydrogen on the light of distant quasars, place stringent lower bounds on the boson mass when ultralight bosons constitute all of the dark matter. Depending on the dataset and modeling assumptions, masses below roughly 10 to the minus 21 up to 2 times 10 to the minus 20 electronvolts are disfavored, excluding the canonical fuzzy dark matter scale of 10 to the minus 22 electronvolts. The new study shows that repulsive self-interactions can indeed rescue the equilibrium solutions: for scattering lengths around 10 to the minus 69 meters, ground-state branches extend into the mass range of roughly 2 times 10 to the minus 20 to 10 to the minus 18 electronvolts, comfortably compatible with the Lyman-alpha bounds. But there is a catch. Those high-mass, self-interacting configurations generate circular velocities far too small to match the dwarf-galaxy rotation curves that the lighter bosons fit so well. Within the class of stationary halos examined, the tension between the masses preferred by galactic dynamics and those demanded by the Lyman-alpha forest remains unresolved.

The broader significance of the work lies in its framework rather than any single verdict. By keeping the particle number explicit and mapping the full solution space, the authors have produced phase diagrams and scaling relations that connect fundamental boson properties directly to observable galactic quantities, offering a controlled laboratory for testing Bose–Einstein condensate dark matter models. The authors note that future work should probe the dynamical stability of the ground-state and excited branches, extend the analysis to halos with outer envelopes, and confront a much larger sample of rotation curves. For now, the message is twofold: quantum-wave halos remain a viable and elegant description of dwarf galaxies at the canonical mass scale, but escaping the tightening noose of Lyman-alpha constraints through self-interactions may cost the model the very galaxies it explains best. The universe, it seems, is still auditing every proposal for its missing matter.

Subject of Research: Equilibrium solutions of the Gross–Pitaevskii–Poisson system as a model for self-gravitating Bose–Einstein condensate dark matter halos

Article Title: Equilibrium halo solutions of the Gross–Pitaevskii–Poisson system: the role of the particle number

Article References: Guzmán, F. A., Castellanos, E., & Mastache, J. (2026). Equilibrium halo solutions of the Gross–Pitaevskii–Poisson system: the role of the particle number. The European Physical Journal C, 86(9), Article 1103. https://doi.org/10.1140/epjc/s10052-026-16373-w

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16373-w

Keywords: dark matter, Bose–Einstein condensate, Gross–Pitaevskii–Poisson system, ultralight bosons, fuzzy dark matter, galaxy rotation curves, solitonic core, Thomas–Fermi limit, Lyman-alpha forest, self-interaction, halo scaling relations, computational astrophysics