A Novel Mathematical Formalism for Modeling Physical Phenomena

1. Introduction

This paper is about physics, an essential part of which is the activity of knowledge exchange among physicists, at different locations, different eras, different orientations etc. The Lagrangian formalism, also referred to as the (extremum-) action principle, is one of several equivalent tools designed to achieve the common ground necessary for such social activity. Physicist A (mathematically) representing a studied system by $\phi$, is guaranteed that $\phi$ could (in principle) appear in physicist B’s notes if they use a common action to generate the set of their permissible $\phi$’s. To facilitate communication between any two physicists there must also exist a consistent set of dictionaries, translating ${\phi }_{\mathrm{A}}\leftrightarrow {\phi }_{\mathrm{B}}$, which is elegantly provided by the symmetry group of the action. The central role of the Poincaré group in this regard stems from the fact that it provides a necessary and sufficient set of such dictionaries for the vast majority of physicists ever registered. In this paper we ask: Why not extend our Poincaré community of physicists to include also physicists of arbitrary scale?

A cynical reviewer might at this point recommend that this manuscript be resubmitted to a scaled journal, for which there are two good replies. First, to this very day the social activity of physicists is limited to the firm ground of planet earth and at small relative velocities. Yet the mere act of imagining the existence of physicists anywhere else and at large velocities relative to us is what brought as so far, rescuing physics multiple times from long periods of stagnation. Rejecting even the possibility that physicists could have a size other than that of a physics professor, would be a blatant repetition of the original sin of anthropocentrism. So why not imagine a giant observer for whom our galaxy, or even the entire universe is a mere speck of dust? Or miniature ones, experiencing the creation and subsequent annihilation of a short-lived subatomic particle over multiple generations?—which leads to the second reply: We currently don’t know how to identify a scaled physicist. The cynic might object that we do: Just take all the dynamical fields, $\phi$, and scale them according to

$$\phi \left(x\right)\mapsto S_{\lambda }\phi ={\lambda }^{\alpha }\phi \left(\lambda x\right)$$

(1)

for some $\phi$-specific scaling exponent $\alpha$ (its inverse ‘length dimension’); that’s the only way to preserve the multiplicative group property

$$S_{{\lambda }_2}{S}_{{\lambda }_1}\phi =S_{{\lambda }_2{\lambda }_1}\phi$$

(2)

he would argue. This implies that scaled physicists comprise scaled hydrogen atoms, which are seen nowhere. Moreover, scaling (1) without alteration of the constants of nature requires physics to be scale covariant, which it isn’t according to our best understanding.

A detailed model addressing the cynic’s concerns has previously been proposed by the author. It accepts the premise (1) hence also the conventional action principle which was never meant to be invariant under (1). There is only a handful of nontrivial scale invariant actions, none of which come close to being realistic. Attempting realism therefore required a very unorthodox application of the action principle manifesting in various technical subtleties, which is never a good sign. However, this is not the main motivation for the current paper. Rather, the form (1) of a scale transformation is too simplistic for a couple of reasons: First, it is only one part of what occurs when, e.g., zooming out of a picture, the other part being a coarsening/smoothing operator. Such coarsening is familiar from the Renormalization Group formalism where $\phi$ is assumed to be a scale-dependent effective representation of some fundamental underlying reality. However, since we can’t allow such assignment of ontological privilege to any particular scale, $\phi (x,\lambda )$ are equally fundamental irrespective of their $\lambda$. Second, even standalone, (1) presupposes too much about $\phi$. The Hubble expansion, for example, can formally be viewed as a scale transformation satisfying (2), with the cosmological time playing the role of (the log of-) $\lambda$, but not conforming with (1), in which different structures scale differently depending on their ‘$\phi$’. Moreover, (1) doesn’t admit a generally covariant extension, which is a prerequisite for any realistic theory. As (1) results from integrating infinitesimal ‘naive’ scale transformations

$$\lambda {\partial }_{\lambda }\phi =x·{\partial }_x\phi +\alpha \phi ,$$

a more flexible rescaling would ensue from substituting $x·{\partial }_x\phi \mapsto Z·{\partial }_x\phi$ (generalizing to a Lie derivative along Z for tensors) where Z is a $\phi$-dependent scaling field determined on consistency grounds. The $\alpha$-term could likewise locally depend on Z. Combined with some $\lambda$-independent, local coarsening operator $\widehat{C}$, (1) is replaced with

$$\lambda {\partial }_{\lambda }\phi =\widehat{S}\phi ,\mathrm{with}\widehat{S}\phi :=Z·{\partial }_x\phi +\alpha \phi +\widehat{C}\phi$$

(3)

That prescription (3) for scale transformations respects the group property (2) is easily seen by changing the scale variable to $s=ln\lambda$, $s\in \left(-\infty ,\infty \right)$, in which case (3) becomes

$${\partial }_s\phi =\widehat{S}\phi$$

(4)

and (2) becomes the group property of a flow. Above and throughout the paper $\phi$ stands for $\phi (e^s,x)$ whenever the logarithmic scale s is involved, which should be clear from the context. Crucially, while the effect of $\widehat{C}$, as that of its RG counterpart, is to smooth $\phi$ (equivalently, attenuate its high frequencies) it must not result in a projection, or else scale-flow would be possible only in the forward, viz. $+s$ direction, with some $\phi (x,s_0)$ as initial condition, implicitly privileging the scale $s_0$.1 Nonetheless, since $\widehat{C}$is a coarsener, flowing backwards in scale typically leads to a singularity, often at finite-s, as the high frequencies grow without bound. It further turns out that compatibility with Lorentz transformations creates a similar problem also in the $+s$ direction. In the proposed formalism, the tiny subset of solutions of (3) which are well-behaved at any scale, $s=\pm \infty$ included, and any x, denoted by $\mathcal{S}$, plays the role of the set of all extrema of an action.

The set $\mathcal{S}$ is determined solely by the form of $\widehat{S}$ in (3). In other words, given a definition of what a scale transformation is, the mere requirement of consistency between the descriptions of physical phenomena at any scale is what defines the laws of physics. Thus each member of $\mathcal{S}$, referred to as a scale orbit, consists of infinitely many fixed-scale sections, each corresponding to distinct yet compatible representations of spacetime phenomena at different scales. It should therefore not come as a surprise that analyzing an individual section not in the context of its full orbit could lead to ‘bizarre physics’. Critically, unlike in the RG formalism, the scale, s, is not a resolution parameter an experimenter can always control by changing the equipment with which he observes a system, but rather his native scale, encapsulating the totality of instruments and materials he uses when arriving at measurement results. The native scale is an identifier facilitating the dictionary between distinct-scale physicists. A dwarf and a giant can determine whether they are studying the same system, i.e., whether their sections are taken from a common orbit, by propagating in scale their sections at their native scale, to the native scale of the other. Note that by s-translation invariance of (4) one’s native scale is only defined up to a communal constant, i.e. only relative scales matter. Thus a spacetime phenomenon (section) which we, humans, regard as being in the realm of condensed matter physics, a dwarf might label “astrophysical", and so would be his attitude towards us, when slicing the orbit on which we reside at ${\lambda }_{\mathrm{dwarf}}$. However, when slicing his orbit at ${\lambda }_{\mathrm{dwarf}}$ he must arrive at a self-representation which is isomorphic to ours, i.e., ${\phi }_{\mathrm{dwarf}}(x,s_{\mathrm{dwarf}})\sim {\phi }_{\mathrm{human}}(x,s_{\mathrm{human}})$, or else we would not belong to the same community of physicists (Such distinction between the representations of one’s self and of others exists also in action based theories, perhaps the most radical example being a physicist boosted to near light-speed becoming nearly two dimensional). The s-translation invariance of (4) then implies ${\phi }_{\mathrm{dwarf}}(x,s)\sim {\phi }_{\mathrm{human}}(x,s-\mathsf{\Delta }s)$ with $\mathsf{\Delta }s=s_{\mathrm{dwarf}}-s_{\mathrm{human}}$.

Irrespective of its philosophical foundations, the proposed formalism could be used as a phenomenological tool for modeling physical phenomena, no more arbitrary than PDE’s and the like, with multi-scale phenomena, e.g., turbulent flow, being the most natural candidates. However in this paper the main focus is on a model leaning towards the ‘fundamental’. Although special emphasis is put on gravity, there is no escape from touching other fields, pertaining to the nature of matter, as the proposed model cannot ‘borrow’ the representation of (gravitating) matter from standard theories; it must be able to represent matter internally (and hopefully more consistently). Obviously, no single paper and no single brain can cover the full breadth of relevant knowledge. Moreover, the high degree of novelty involved in the proposed formalism meant that, compromises in mathematical rigor in some of the proofs/arguments had to be made. Nonetheless, consistent equations for observables are eventually derived, and solved in simple cases. With such panoramic view of physics through a new lens, albeit of low resolution, multiple aspects of the proposed formalism are explored, enabling novel approaches to longstanding open problems in physics, which are not even expressible in standard language. Finally, to remain at a reasonable length the paper begins with a simple toy model introducing some of the central concepts of the formalism. With that background the reader should be able to fill in certain gaps in the subsequent derivations.

2. Exactly Solvable Linear Toy Model

In order for the notion of native scale to be fully meaningful, $\phi$ must be rich enough to be able to describe: the system being observed; the observer—his equipment included; electromagnetic phenomena involved in most observations etc. This ambitious task is deferred to Section 3. A gentle introduction to the jargon and techniques used in that section is provided by the flow (4) of a time-independent scalar field in Euclidean D-dimensional space. In choosing the generator of coarsening, $\widehat{C}$, the following properties should be included:

• Averaging. If $\phi \left(x_{\mathrm{m}}\right)$ is a local maximum (minimum) then $\widehat{C}\phi \left(x_{\mathrm{m}}\right)\le 0$ ($\ge 0$ resp.)
• Locality. $\widehat{C}$ is second order and does not contain higher order derivatives or powers higher than the first of the second derivative.
• Equivariance. $\widehat{C}$ must commute with translations, rotations and reflections in ${\mathbb{R}}^D$.

The simplest $\widehat{C}$ satisfying the above is the (D-dimensional) Laplacian, corresponding to the (weighted) arithmetic average of $\phi$ in the neighborhood of a point. There are, of course, other choices corresponding to different averages2 such as: $\nabla ·\left(\nabla f\left(\phi \right)\right)$ for a monotonically increasing f, or ${\nabla }^2\phi +{|\nabla g\left(\phi \right)|}^2$ for any g and combinations thereof (∇ and · are both D-dimensional). The locality clause is a corollary of averaging. Indeed, in 1-dimension for simplicity, a ${\partial }_x^4$ term added to the Laplacian would increase a local maximum of $\phi =-x^2+a{x}^4$ at $x=0$ for $a>2\times 4!$, as would an added ${\left({\partial }^2\phi \right)}^2$ do to $\phi =-a{x}^2$ for $a>\frac{1}{2}$. However, when $\phi$ takes values in spaces lacking a clear-cut definition of local extremum, locality becomes an independent clause, defining the local, infinitesimal neighborhood of a point. Clause 3 is obviously needed due to the arbitrariness in positioning and orienting one’s coordinate system. Equivalently, it is what defines the community of physicists in Euclidean space.

Sticking with the Laplacian, and using the simplest scaling field, $Z^i=x^i$, (4) becomes

$${\partial }_s\phi ={\ell }_0^2{\nabla }^2\phi +x·\nabla \phi +\alpha \phi$$

(5)

with ${\ell }_0$ some parameter. It is tempting to attribute a ‘physical dimension of length’ to ${\ell }_0$, balancing the double derivative it multiplies. However, being a description of physics on all scales, the proposed formalism is inexorably an attempted ‘theory of everything’ and as such ought to be able to represent any measurement process. And since the result of any measurement is ultimately a dimensionless number, e.g., the number pointed at by a pointer, or the minimal number of standard-length rods exactly fitting a line segment, the notion of physical dimensions should ultimately be abolished. Moreover, since (5) describes a flow in scale, endowing ${\ell }_0$ with a dimension of length may lead to the wrong expectation that it too would flow in scale. Nonetheless, the developmental stage of the proposed theory is currently insufficient to internally represent any measurement. To make contact with empirical data associated with sections at our native scale, arbitrarily assigned the value $s=0$ or $\lambda =1$, dimensions will occasionally appear in this paper. Unless stated otherwise, ${\ell }_0=1$ is assumed, i.e., the coordinate x at $s=0$ is measured in multiples of ${\ell }_0$. Note that even this innocuous statement relies on the existence of an affine structure of space whose physical validation requires an affine structure of space! Thus without doing away with this circularity via a general covariant extension of (5) (Sec. 3), our proposal cannot even pretend to be a fundamental physical theory.

2.1. The Particle Basis of $\phi$

Of special interest are fixed-point solutions of the flow (5), i.e. scale-invariant $\phi$, of which fixed-points which are further global or local attractors stand out. To find the latter we note that, if $\phi$ is integrable at $s=0$—the case of a non-integrable $\phi$ is dealt with later—its zeroth moment, $m_0$, satisfies ${\partial }_s{m}_0=(\alpha -D)m_0$ and $m_0$ explodes for $\left|s\right|\to \infty$, implying $\phi \notin \mathcal{S}$, unless $\alpha =D$ and $m_0\left(s\right)\equiv 1$ without loss of generality by the linearity of (5), or else $m_0\left(0\right)=0$. Assuming the former for now, it is helpful to represent $\phi$ by its cumulants. Taking the Fourier transform of (5) and dividing by the Fourier transform $\phi :=\mathcal{F}\phi$, which is also the generating function of its moments assumed all to exist, leads to the following equation for the generating function of the cumulants, $\mathcal{Z}(k,s):=ln\phi$

$${\partial }_s\mathcal{Z}=-k^2-k·\nabla \mathcal{Z},(v·w\equiv v_i{w}_i,v^2\equiv v·v)$$

(6)

Equation (6) is an infinite set of uncoupled o.d.e.’s for the coefficients of multinomials $k_1^{n_1}\dots k_D^{n_D}$, whose solutions $\propto exp\left(-s{\sum }_i{n}_i\right)$ all vanish for $s\to \infty$ except those of $k_1^2,\dots ,k_D^2$ which approach $-\frac{1}{2}$. Solving back, $\phi ={\mathcal{F}}^{-1}exp\mathcal{Z}$, we get the Gaussian

$${\phi }_{\mathrm{G}}={\left(2\pi \right)}^{D/2}{e}^{-\frac{1}{2}{\sum }_i{x}_i^2}$$

c flow (5) for (moment-determinate) functions with well defined moments to any order at $s=0$. However, not all moment determinate $\phi (x,s=0)$ are sections of orbits in $\mathcal{S}$. A large subset that is, consists of linear combinations, discrete or continuous, of shifted fixed-point Gaussians. This is so because any such sum has well defined moments to any order and the effect of the flow on individual Gaussians is a trivial shift

$${\phi }_{\mathrm{G}}(x-x_0){|}_{s=0}\mapsto {\phi }_{\mathrm{G}}\left(x-e^{-s}{x}_0\right)$$

(7)

We strongly suspect that the converse is also true, namely, that any $\phi \in \mathcal{S}$ can be decomposed into a linear combination of shifting Gaussians. A formal proof will not be attempted here but the intuition must be clear: The Laplacian increases (decreases) and ‘sharpens’ $\phi$ at its local maximum (resp. minimum) when flowing in the $-s$ direction of (5), and if the scaling part does not sufficiently compensate for this, $\phi$ wildly diverges at a local extremum. Any section of $\phi \in \mathcal{S}$ must therefore be ‘round’ enough on the scale set by ${\ell }_0$($=1$) and if so it should be decomposable into Gaussians of width ${\ell }_0$. A clear illustration of the fate of a function not sufficiently round is provided by a Gaussian of width less than 1 at $s=0$, i.e., $\mathcal{Z}\left(0\right)=-(\frac{1}{2}-ϵ)k^2$. A finite-s singularity is reached at $s=\frac{1}{2}ln\left(2ϵ\right)$ where $\mathcal{Z}$ vanishes, corresponding to $\phi$ which is a ‘delta-function Gaussian’.

Closure under continuous sums is what distinguishes the shifted Gaussians basis of $\mathcal{S}$, referred to as the ‘particle basis’, from the scaled Fourier basis,

$${\lambda }^{\alpha }{e}^{i\lambda k·x-\frac{1}{2}{\lambda }^2{k}^2};k\in {\mathbb{R}}^D,\lambda \equiv e^s$$

(8)

Although individually in $\mathcal{S}$, infinite sums thereof may still lie outside $\mathcal{S}$, as the above narrow-Gaussian example demonstrates. Nonetheless, the pseudo basis (8) is not entirely useless. It clearly shows the rapid decay of waves when their wave-vector is contracted beyond the cutoff frequency (${\ell }_0^{-1}$) and will serve us in the sequel.

Returning to the case of non-integrable, or integrable but zero $m_0$, the corresponding fixed-points of (5) are $f_l\left(r\right)Y_l^m\left(\mathsf{\Omega }\right)$ in $D=3$ for a suitable $f_l$, $l\ge 1$, vanishing at $r=0$, which happen to have a $r^{-3}$ asymptotic tail hence are non-integrable. It can be shown that the basin of attraction of each consists only of itself, rendering it uninteresting from our perspective.

Full justification for the name “particle" attached to Gaussians of width ${\ell }_0$ and $0^{\mathrm{th}}$-moment equal to 1 (or any other normalization) will have to await Section 3 but some can already be given at this stage. If physicists of different native scales are to have isomorphic self-representations—their laboratory etc. included—and if physicists exist at arbitrarily large or small scales (but not necessarily at any scale) then they must all consist of the same particles and their ‘oppositely charged’ antiparticles—Gaussians of width ${\ell }_0$ and $0^{\mathrm{th}}$-moment $-1$. Otherwise particle-antiparticle pairs would not get fully annihilated when flowing via (7) in the $+s$ direction. As a result, the ‘vacuum’ would get increasingly contaminated with particles of arbitrary charge. Conversely, the vacuum could only acquire content when zooming into an empty patch of it, if particle-antiparticle pairs are created out of it. Note that the decomposablility of $\phi$ into a discrete sum of particles is required by the fact that $\phi$ would otherwise trivialize to a uniform $\phi \equiv 0$ for $s\to -\infty$. Thus our model requires for its consistency both particles and the quantization of their charge. Moreover, in the point-particle limit, ${\ell }_0\to 0$, it would take $s\to \infty$ for any pair to annihilate even approximately, contradicting the existence of scaled physicists. Point-particles, which are the source of all evil in mathematical physics, are excluded from the outset.

2.2. Adding Time Dependence

How should the flow (5) be generalized for a time-dependent $\phi$? Guided by the equivariance clause with the Poincaré group replacing the isometry group of Euclidean space, the unique generalization reads

$${\partial }_s\phi ={\nabla }^2\phi -{\partial }_{tt}\phi +x·\nabla \phi +t{\partial }_t\phi +\alpha \phi$$

(9)

Multiplying the ${\partial }_{tt}\phi$ term is a ${({\ell }_0/c)}^2$ coefficient, assumed to equal 1 unless stated otherwise. In covariant notations (9) reads

$${\partial }_s\phi ={\square }^2\phi +x^{\mu }{\partial }_{\mu }\phi +\alpha \phi$$

(10)

with ${\square }^2\equiv {\eta }^{\mu \nu }{\partial }_{\mu }{\partial }_{\nu }$ and ${\eta }^{\mu \nu }={\eta }_{\mu \nu }:=\mathrm{diag}(\text{-}1,1,1,1)$. Readers experiencing unease from the appearance of Lorentz symmetry out of the blue are referred to [3]; It has been known since the early years of Relativity Theory that Lorentz transformations were only serendipitously discovered within the framework of electrodynamics. This symmetry group (the Galilean group being a limiting case thereof) is an inevitability when the meaning of synchronized clocks is logically analyzed—which is essentially what is being done in this paper with regard to scale transformations. Thus our proposed time-dependent generalization of ${\mathbf{\nabla }}^2\phi (\mathit{x},t)$ inexorably involves values of $\phi$ at times other than t, but not because space and time form a ‘spacetime’ continuum—the Minkowskian/geometric view (if that were the case then ${\partial }_{tt}$ would be expected to have the opposite sign). Instead, space and time are fundamentally distinct and are mixed together on consistency grounds; without such mixing no community of physicists would exist.

Generalizing static particles are—naturally—moving particles and in particular uniformly so, obtained by boosting a static particle solution. In $D=1$ for simplicity, the boost explicitly reads ${\phi }_{\mathrm{G}}(x,s)\mapsto {\phi }_{\mathrm{G}}\left(\gamma (x-vt),s\right)$ with $\gamma =\sqrt{1/(1-v^2)}$, making the Lorentz contraction of the particle in the direction of motion manifest. Uniformly moving particles are therefore all members of $\mathcal{S}$, but are all members of $\mathcal{S}$ such? More accurately: Do they form a basis for $\mathcal{S}$? To answer this question, (9) is first integrated over three-space. Assuming the integral exists results in the following equation for the zeroth moment

$${\partial }_s{m}_0=-{\partial }_{tt}{m}_0+t{\partial }_t{m}_0+(\alpha -D)m_0$$

(11)

Plugging a scaled Fourier ansatz $m_0(\omega ;t,s)=f_{\omega }\left(s\right)exp\left(i{e}^s\omega t\right)$ and continuing with the $\alpha$=D case, we get: ${\displaystyle \frac{\mathrm{d}}{\mathrm{d}s}}{f}_{\omega }={\omega }^2{f}_{\omega }$ implying $f_{\omega }\equiv 0$$\forall \omega \ne 0$ or else it explodes for $s\to \infty$. For $\omega =0$, $f_{\omega }$ is some constant which can be assumed to equal 1 by the linearity of (9). We conclude that $m_0$ which is constant in both time and scale is a necessary condition for the corresponding $\phi$ to lie in $\mathcal{S}$. Next, consider the generalization of equation (6) for the time-dependent cumulants of $\phi$,

$$\lambda {\partial }_{\lambda }\mathcal{Z}=-k^2-{\partial }_{tt}\mathcal{Z}+{\left({\partial }_t\mathcal{Z}\right)}^2+t{\partial }_t\mathcal{Z}-k·\nabla \mathcal{Z}$$

(12)

similarly obtained by Fourier transforming (9) and dividing by $\phi (k,t,\lambda )$. Using (12) we first argue that the instability of the flow (9) in the $-s$ direction, which mandates particles, in and of itself does not further mandates their uniform motion. To show this we plug the following ansatz into (12), continuing with $D=1$ for clarity

$$\mathcal{Z}(k,\lambda ,t)=\sum _{n=1}^{\infty }{c}_n(\lambda t,\lambda )k^n=\sum _{n=1}^{\infty }\left(\sum _{m=n-1}^{\infty }{c}_n^m\left(\lambda t\right){\lambda }^{2m-n}\right)k^n$$

(13)

where $c_n^m$ is a double-indexed function of the scaled time alone. Equating the coefficient of each power of k to zero in increasing powers of $\lambda$ (note that the series is missing the $n=0$ term by our result $c_0\equiv ln{m}_0=0$). Starting with $c_1^{-1}$, dictating the asymptotic scaling form of the center-of-mass, i.e. ${\lambda }^{-1}{c}_1^{-1}\left(\lambda t\right)$, we see that it decouples from all other terms in the limit $\lambda \to 0$ and can be an arbitrary function. Successive terms, $c_1^k$, can then be iteratively computed as the ${\left({\partial }_{t}Z\right)}^2$ term does not contain the first power of k, and the ${\partial }_{tt}$ term pulls out an extra factor of ${\lambda }^2$, e.g., $c_1^1=-\frac{1}{2}{\ddot{c}}_1^{-1},c_1^3=-{\displaystyle \frac{1}{4}}{\ddot{c}}_1^1,\dots$ (all evaluated at the scaled time $\lambda t$). For $c_1^{-1}$ with bounded time derivatives to all orders, the power series of $c_1$ clearly converges for $\lambda <1$. Moving to $n=2$, the leading order of (minus the-) variance reads $c_2^0=-\frac{1}{2}(1-v^2)$ with $v={\dot{c}}_1^{-1}$ again manifesting the Lorentz contraction. Higher order corrections, $c_2^k$, can then be calculated in terms of ${\ddot{c}}_2^{k-2}$ and products ${\dot{c}}_1^p{\dot{c}}_1^q$, $p+q=k-2$ coming from the nonlinear term, e.g., $c_2^2={\displaystyle \frac{1}{4}}\left(-{\ddot{c}}_2^0+2{\dot{c}}_1^{-1}{\dot{c}}_1^1\right)\dots$. The leading order term of the third cumulant, the so called “skewness", reads $\lambda c_3^1=\lambda {\displaystyle \frac{2}{3}}{v}^2\dot{v}$. It is likewise a relativistic effect in which the ‘front’ and ‘back’ of an accelerating (extended) particle experience different Lorentz contractions. Continuing this way, $c_n^m$ can be computed in terms of $c_1^{-1}$ and its first $n-1$ derivatives. The results in a convergent power series for each at $\lambda <1$. As $c_n=O\left({\lambda }^{n-2}\right)$, at small enough scale the shape of a particle is approximated arbitrarily well by that of a uniformly moving particle (appropriately Lorentz contracted in the direction of motion). The leading order correction appear in the form of a particle’s skewness in the direction of acceleration (predicated on relativistic velocities). We conclude that a sufficient condition for $\phi$ to well-behave at small scales is that it approaches a sum of scaling particle solutions, generalizing the static-particle result. Although more difficult to prove, there is no apparent reason why this conclusion should not carry to nonlinear scale-flows in the case of a single particle.

Returning to our original question, of whether $c_1^{-1}$ necessarily describes a uniformly moving particle, we turn to the fate of such a well localized moving particle solution at small $\lambda$, entirely encoded in the single function $c_1^{-1}$, when it flows to large $\lambda$, outside the convergence radius of the each cumulant’s power series. We prove that the answer is positive, viz., unless $c_1^{-1}$ describes a globally freely moving particles, ${\lambda }^{-1}{c}_1^{-1}\left(\lambda t\right)={\zeta }_{\mathrm{free}}:={\lambda }^{-1}{x}_0+vt$ for some $x_0$ and v, its corresponding $\phi$ is not in $\mathcal{S}$. This is due to a new flow instability in the $+\lambda$ direction, introduced by the minus sign of the ${\partial }_{tt}$ term in (9). To prove this, consider the equation for the center of the particle $\zeta \equiv c_1$ (switching notations in order to not overload the upper index)

$$\lambda {\partial }_{\lambda }{\zeta }^i=-{\partial }_{tt}{\zeta }^i+t{\partial }_t{\zeta }^i-{\zeta }^i,i=1,\dots ,D$$

(14)

obtained by equating the coefficients of $k^i$ in (12) to zero (i not a power!). Now plug into (14) the most general solution

$${\zeta }^i(t,\lambda )={\zeta }_{\mathrm{free}}^i(\lambda ,t)+{\lambda }^{-1}\sum _{\omega }{A}_{\omega }^i{e}^{i\lambda \omega t+\frac{1}{2}{\omega }^2{\lambda }^2}$$

(15)

with the sum representing also an integral, and the $A_{\omega }$’s are the Fourier coefficients of ${\zeta }^i(1,t)$. Clearly, unless $A_{\omega }^i\equiv 0$, and insofar as $\phi$ still describes a localized particle, this particle (wildly) moves around unbounded for $\lambda \to \infty$ which in and of itself implies $\phi \notin \mathcal{S}$; any ‘wrinkle’ in a straight path at $\lambda =1$ gets rapidly amplified at large $\lambda$ by the $-{\partial }_{tt}$ term in (14). However, moving to higher order cumulants, which are all morphological attributes of a particle-like $\phi$ hence independent of ${\zeta }^i$, a similar divergence occurs, implying either the divergence $\phi$ or its complete delocalization.

Another way of seeing why only uniformly moving particles appear in $\mathcal{S}$ is by decomposing a non-uniformly moving particle solution at $s=0$ into its space-time Fourier components, and letting them each flow to $s=\infty$. Their evolution in scale is just (8) with $k·x$ meaning $k_{\mu }{x}^{\mu }$ (and $k^2\equiv k_{\mu }{k}^{\mu }$). Waves with $k^2>0$ are strongly attenuated at large $\lambda$, while those with $k^2<0$ blow-up. Now, it is easily verified that any non-uniform, or uniform but superluminal motion at $s=0$, must have some time-like ($k^2<0$) Fourier components in its decomposition, expelling the orbit on which it resides from $\mathcal{S}$. This method is applicable also to possible superpositions of non-uniformly moving particles having $\zeta ={\zeta }_{\mathrm{free}}$ for their joint $\phi$. Rather than resorting to murky causal paradoxes, or to our current inability to accelerate masses beyond the speed of light, the proposed formalism rejects Tachyons on simple mathematical grounds. The dominance of waves with light-like k’s can also be appreciated even before moving to more complicated models.

3. A Realistic Model

The alert reader must have anticipated the main result of the previous section, namely, that $\mathcal{S}$ consists of freely moving particles. By linearity, particles can move through one another uninterrupted and if so, they are noninteracting particles which should better have straight paths. Enabling their mutual interaction therefore requires some form of nonlinearity, either in the coarsener, $\widehat{C}$, or in the scaling part. Further recalling our commitment to general covariance as a precondition for any fundamental physical theory, nonlinearity is inevitably and, in a sense, uniquely forced upon us. A nonlinear model also supports a plurality of particles, having different sizes which are different from the common ${\ell }_0$ in a linear theory. This frees ${\ell }_0$, ultimately estimated at $\sim {10}^{20}$ km, to play a role at astrophysical scales.

Our realistic model involves a spin-1, $\alpha =1$ field: The A-field. A straightforward (and unique up to terms involving the curvature) way to render a differential operator generally covariant is through the minimal-coupling prescription, of ‘dotting the commas’ which can be applied to the ‘Maxwell coarsener’ ${\eta }_{\mu \nu }{\square }^2-{\partial }_{\mu }{\partial }_{\nu }$. The scaling piece of an $\alpha =1$ covariant vector is just its Lie derivative with respect to the scaling field, $Z^{\mu }$, defined below. Combined, the scale flow of $A_{\mu }$ reads

$${\partial }_s{A}_{\mu }={\ell }_0^2\left({\nabla }^{\nu }{\nabla }_{\nu }{A}_{\mu }-{\nabla }^{\nu }{\nabla }_{\mu }{A}_{\nu }\right)+{\left({\mathcal{L}}_{Z}A\right)}_{\mu }$$

(16)

with

$${\left({\mathcal{L}}_{Z}A\right)}_{\mu }\equiv Z^{\nu }{\nabla }_{\nu }{A}_{\mu }+{\nabla }_{\mu }{Z}^{\nu }{A}_{\nu }$$

(17)

the Lie derivative of $A_{\mu }$ with respect to $Z^{\mu }$. In flat spacetime and $Z^{\mu }\left(x\right)=x^{\mu }$, (17) reduces to ‘naive’ $\alpha =1$ scaling, ${\left({\mathcal{L}}_{Z}A\right)}_{\mu }\to x^{\nu }{\partial }_{\nu }{A}_{\mu }+A_{\mu }$. Equation (16) prescribes the generally covariant scale-flow of a vector in one particular coordinate system common to all scales. Thus $\mathcal{S}$ is partitioned into equivalence classes, the members of each are related by some coordinate transformation.

Analyzing an $\alpha =1$ model in $D=3$ space is much more difficult as no moment of the associated particle exists. This is clearly seen already in a linear, naive scaling, flat-space model. A non-spinning, viz., $\mathit{A}\equiv 0$, fixed-point reads

$$A_0\left(r\right):={\displaystyle \frac{1}{r}}\mathrm{erf}\left({\displaystyle \frac{r}{\sqrt{2}{\ell }_0}}\right)$$

(18)

having a non-integrable $r^{-1}$ tail. To facilitate the analysis of such extended particles, we define a auxiliary $\alpha =3$ model for the center(oid) of $A_{\mu }$

$$J_{\mu }:=-{\nabla }^{\nu }{\nabla }_{\nu }{A}_{\mu }+{\nabla }^{\nu }{\nabla }_{\mu }{A}_{\nu }\equiv -{\nabla }^{\nu }{F}_{\nu \mu }\equiv -{\displaystyle \frac{1}{\sqrt{-g}}}{\partial }^{\nu }\left(\sqrt{-g}{F}_{\nu \mu }\right)$$

(19)

Operating with ${\nabla }^{\mu }$ on (19), using the antisymmetry of $F_{\mu \nu }:={\nabla }_{\mu }{A}_{\nu }-{\nabla }_{\nu }{A}_{\mu }$, the commutators of covariant derivatives

$$\left[{\nabla }_{\mu },{\nabla }_{\nu }\right]V^{\alpha }=-R_{\phantom{1}\rho \mu \nu }^{\alpha }{V}^{\rho },\left[{\nabla }_{\mu },{\nabla }_{\nu }\right]T^{\alpha \beta }=-R_{\phantom{1}\rho \mu \nu }^{\alpha }{T}^{\rho \beta }-R_{\phantom{1}\rho \mu \nu }^{\beta }{T}^{\alpha \rho },\mathrm{etc}.$$

(20)

and the symmetries of the Riemann tensor, gives ${\nabla }^{\mu }{J}_{\mu }\equiv 0$, i.e., $J_{\mu }$ is covariantly conserved at any scale, s. Above and throughout the paper Weinberg’s sign convention for the Riemann tensor is used

$${R^{\mu }}_{\alpha \beta \gamma }=-\left({\mathsf{\Gamma }}_{\alpha \gamma ,\beta }^{\mu }-{\mathsf{\Gamma }}_{\alpha \beta ,\gamma }^{\mu }+{\mathsf{\Gamma }}_{\sigma \beta }^{\mu }{\mathsf{\Gamma }}_{\gamma \alpha }^{\sigma }-{\mathsf{\Gamma }}_{\sigma \gamma }^{\mu }{\mathsf{\Gamma }}_{\beta \alpha }^{\sigma }\right)$$

(21)

Operating with ${\widehat{W}}_{\rho }^{\phantom{1}\mu }:=-{\nabla }_{\alpha }{\nabla }^{\alpha }{g}_{\rho }^{\phantom{1}\mu }+{\nabla }^{\mu }{\nabla }_{\rho }$ on (16), the second term of this operator almost annihilates the coarsener by the above remarks (charge conservation), leaving only a curvature term $R_{\nu \mu }{J}^{\mu }$ coming from commutator (20). The scaling piece combines the covariant generalization of the flat spacetime conversion $\alpha =1\mapsto 3$ with a novel nonlinear term (see Sec. 3.2 below). The l.h.s. reads ${\widehat{W}}_{\rho }^{\phantom{1}\mu }{\partial }_s{A}_{\mu }$. We would like to swap the order of $\widehat{W}$ and ${\partial }_s$, which would give ${\partial }_s{J}_{\nu }$ by (19). However, in curved spacetime ${\widehat{W}}_{\rho }^{\phantom{1}\mu }$ could implicitly depend on the scale s through $g_{\mu \nu }$, introducing a commutator correction. The combined result reads

$${\partial }_s{J}_{\nu }+\left[{\widehat{W}}_{\nu }^{\phantom{1}\mu },{\partial }_s\right]A_{\mu }={\ell }_0^2\left({\nabla }_{\alpha }{\nabla }^{\alpha }{J}_{\nu }+R_{\nu \mu }{J}^{\mu }\right)+{\widehat{W}}_{\nu }^{\phantom{1}\mu }{\left({\mathcal{L}}_{Z}A\right)}_{\mu }$$

(22)

more suited for analysis. For example, in the case of flat spacetime and naive scaling, $Z^{\mu }=x^{\mu }$, the particle solution of the J-field associated with (18) is the familiar Gaussian $J_0=e^{-r^2/2{\ell }_0^2}$ from the previous sections. It is emphasized that $J_{\mu }$ and its associated scale flow are merely analytic tools, not to be put on equal footing with $A_{\mu }$ and its flow (As the nonlinear terms arising from scaling and commutator do not involve $J_{\mu }$ but rather $A_{\mu }$, for a given $A_{\mu }$ and $g_{\mu \nu }$, (22) describes the linear but inhomogeneous—in both spacetime and scale—flow of $J_{\mu }$). Moreover, infinitely many other markers of an A-particle’s center exist, involving, e.g., curvature terms, which may turn out to be more suitable in certain domains. However, for the purpose of this introductory paper, focusing on weak gravitational fields/flat-spacetime, $J_{\nu }$ turns out to be a most useful choice.

Relation (19) is formally equivalent to Maxwell’s equations with $J_{\mu }$ sourcing $A_{\mu }$’s wave equation. However, $J_{\mu }$ is not an independent object as in classical electrodynamics but a marker of the locus of privileged points at which the Maxwell coarsener does not annihilate $A_{\mu }$ (distinct $A_{\mu }$’s differing by some ${\partial }_{\mu }\mathsf{\Lambda }$ therefore have identical $J_{\mu }$’s). For $J_{\mu }$ and $A_{\mu }$ to mimic those of classical electrodynamics, $J_{\mu }$ must also be localized along curved worldlines traced by solutions of the Lorentz force equation in $A_{\mu }$ (which as already shown in the scalar case, necessitates a nonlinear scale flow). And just like in the scalar-particle case, where higher order cumulants ($n>2$) are ‘awakened’ by its center’s nonuniform motion, deforming its stationary shape, so does the $A_{\mu }$ “adjunct" (in the jargon of action-at-a-distance electrodynamics) to each such $J_{\mu }$ gets deformed. Due to the extended nature of an A-particle, and unlike in $\alpha =3$ models3, these deformations at $(t,\mathit{x})$ are not encoded in the local motion of its center at time t, but rather on its motion at retarded and advanced times, $t\pm \left|\mathit{x}\right|$ (assuming flat spacetime for simplicity). However, associating such temporal incongruity with ‘radiation’ can be misleading, as it normally implies the freedom to add any homogeneous solution of Maxwell’s equations to $A_{\mu }$ which is clearly nonsensical from our perspective. Consequently, the retarded solution cannot be imposed on $A_{\mu }$ and in general, $A_{\mu }$ contains a mixture of both advanced and retarded parts, which varies across spacetime. The so-called radiation arrow of time manifested in every macroscopic phenomenon must therefore receive an alternative explanation (see Section 3.4.2).

Now, why should $J_{\mu }$ be confined to the neighborhood of a worldline? As already seen in the linear, time-dependent case, a scale-flow such as (16) suffers from instability in both s-directions due to the action of the coarsener; Local variations of $A_{\mu }$ along space-like/time-like directions are rapidly amplified in the $-s$/$+s$ direction resp. If we therefore examine the scale flow of $A_{\mu }$ inside a ‘lab’ of dimension much smaller than ${\ell }_0$, centered at the origin without loss of generality, then the coarsener would completely dominate the flow of $A_{\mu }$ whose scale of variation is much smaller than ${\ell }_0$ (unlike (18)!) leading to its rapid divergence, unless it is almost annihilated by ${\widehat{W}}_{\rho }^{\phantom{1}\mu }$, i.e., $A_{\mu }$ varies exclusively in light-like directions. An exception to this rule could take place around privileged points, where the scaling-field grows to such large magnitude as to allow for the scaling piece to counter a large coarsener contribution. This is where $J_{\rho }\equiv {\widehat{W}}_{\rho }^{\phantom{1}\mu }{A}_{\mu }$ is focused and, as shown in Section 3.1.2 below, it is only around world-lines marking the center of an extended particle, where $F_{\mu \nu }$ peaks, that $Z^{\mu }$ can indeed grow to such large magnitude. “Almost" is emphasized above because it is precisely the fact that, at distances from $J_{\mu }$ that are much smaller than ${\ell }_0$, the action of the coarsener is only reduced to the order of that of the scaling piece, which gives rise to deviations from classical physics. This will be a recurrent theme in the rest of the paper.

3.1. Determining the Metric and the Scaling Field

The flow (16) of $A_{\mu }$ requires specifying both $g_{\mu \nu }(x,\lambda )$ and $Z^{\mu }(x,\lambda )$ in a generally covariant way. Starting with the former, we seek the scale flow of $g_{\mu \nu }$. It is well known that, in Riemannian geometry, the Ricci tensor is the unique, symmetric generally covariant tensor which can be constructed from the metric tensor and its first two derivatives and does not contain higher power than the first of its second derivative. By our definition of a coarsener, this leaves $\widehat{C}=R_{\mu \nu }+b{g}_{\mu \nu }R$ with $R_{\mu \nu }:=R_{\phantom{1}\mu \alpha \nu }^{\alpha }$ and $R:=R_{\phantom{1}\mu }^{\mu }$, as the only permissible coarsener for some constant b. Now, the flow (16) of $A_{\mu }$ is ‘guided’ by $g_{\mu \nu }$ via the minimal coupling prescription. On consistency grounds the flow of $g_{\mu \nu }$ must also be guided by $A_{\mu }$, or else the gravitational field would not focus around matter. The simplest way to achieve this mutuality is through the use of the (symmetric) canonical energy-momentum tensor

$${\mathsf{\Theta }}_{\mu \nu }=F_{\mu \rho }{F}_{\nu }^{\phantom{1}\rho }-\frac{1}{4}{g}_{\mu \nu }{F}^{\rho \sigma }{F}_{\rho \sigma }$$

(23)

peaking around $J_{\mu }$, which leads to the following scale flow of $g_{\mu \nu }$:

$${\partial }_s{g}_{\mu \nu }={\ell }_g^2\left(R_{\mu \nu }+b{g}_{\mu \nu }R-4\pi G{\mathsf{\Theta }}_{\mu \nu }\right)+{\left({\mathcal{L}}_{Z}g\right)}_{\mu \nu }$$

(24)

with ${\ell }_g$, G and b some constants, and

$${\left({\mathcal{L}}_{Z}g\right)}_{\mu \nu }\equiv Z^{\rho }{\nabla }_{\rho }{g}_{\mu \nu }+g_{\rho \nu }{\nabla }_{\mu }{Z}^{\rho }+g_{\rho \mu }{\nabla }_{\nu }{Z}^{\rho }={\nabla }_{\mu }{Z}_{\nu }+{\nabla }_{\nu }{Z}_{\mu }$$

(25)

The global sign of the coarsener piece reflects Weinberg’s sign convention for the Riemann tensor (21).

Taking the covariant divergence of (24), and using ${\nabla }^{\mu }\left(R_{\mu \nu }-\frac{1}{2}{g}_{\mu \nu }R\right)\equiv 0$ (as a result of the Bianchi identity) implies that

$${\nabla }^{\mu }\left(-4\pi G{\ell }_g^2{\mathsf{\Theta }}_{\mu \nu }+{\left({\mathcal{L}}_{Z}g\right)}_{\mu \nu }-{\partial }_s{g}_{\mu \nu }+\left(b+\frac{1}{2}\right){\ell }_g^2{g}_{\mu \nu }R\right)=0$$

(26)

is a necessary condition for (24) to have a solution. Defining

$$-4\pi G{\ell }_g^2{P}_{\mu \nu }:={\left({\mathcal{L}}_{Z}g\right)}_{\mu \nu }$$

(27)

equation (26) can be rewritten as

$${\nabla }^{\mu }\left({\mathsf{\Theta }}_{\mu \nu }+P_{\mu \nu }\right)=-{\displaystyle \frac{1}{4\pi G{\ell }_g^2}}{\nabla }^{\mu }{\partial }_s{g}_{\mu \nu }+{\displaystyle \frac{b+\frac{1}{2}}{4\pi G}}{\partial }_{\nu }R$$

(28)

Postulating that energy-momentum conservation is recovered in the limit $g_{\mu \nu }\to {\eta }_{\mu \nu }$ mandates $b=-\frac{1}{2}$, nullifying the last term on the r.h.s. of (28), which is assumed henceforth. The possibility that gravity is essentially involved in the structure of elementary matter [4], and consequently $g_{\mu \nu }\to {\eta }_{\mu \nu }$ is nonphysical, has not been explored.

3.1.1. Determining $Z^{\mu }$ in the flat Spacetime Approximation

Assuming $g_{\mu \nu }={\eta }_{\mu \nu }$ leads to significant simplification when determining $Z^{\mu }$, and is therefore considered first. In this approximation covariant derivatives appear as ordinary derivatives, (28) reduces to energy-momentum conservation

$${\partial }^{\mu }{T}_{\mu \nu }\equiv {\partial }^{\mu }\left({\mathsf{\Theta }}_{\mu \nu }+P_{\mu \nu }\right)=0$$

(29)

and ${\left({\mathcal{L}}_{Z}g\right)}_{\mu \nu }={\partial }_{\mu }{Z}_{\nu }+{\partial }_{\nu }{Z}_{\mu }$. By virtue of definition (19) of $J^{\mu }$ (Maxwell’s equations) and the Maxwell-Bianchi identity ${\partial }_{\mu }{F}^{\mu \nu }\equiv 0$ ($F_{\mu \nu }:={ϵ}_{\mu \nu \rho \sigma }{F}^{\rho \sigma }$) alone, Poynting theorem is satisfied identically

$${\partial }^{\mu }{\mathsf{\Theta }}_{\mu \nu }=-F_{\nu \mu }{J}^{\mu }$$

(30)

hence (29) implies

$${\partial }^{\mu }{P}_{\mu \nu }=F_{\nu \mu }{J}^{\mu }$$

(31)

Maxwell’s equations (19), along with (30) and (31), referred to henceforth as the basic tenets of classical electrodynamics, are nowadays taken as the definition of classical electrodynamics, encapsulating its experimental success while avoiding the traditional use of the Lorentz force equation of a point charge, with its infamous, unresolved classical self-force problem.

Equations (27) and (31) result in four second order equations for the four components of the scaling vector $Z^{\mu }$,

$${\square }^2{Z}_{\nu }+{\partial }_{\nu }{\partial }^{\mu }{Z}_{\mu }=-4\pi G{\ell }_g^2{F}_{\nu \mu }{J}^{\mu }$$

(32)

which, together with the boundary condition $Z^{\mu }\to x^{\mu }$ away from matter, define $Z^{\mu }$ up to a solution of the homogeneous equation (32). A fifth equation ((52) below, derived on the basis of charge conservation) apparently removes this remaining freedom.

Summarizing, as a corollary of defining the scale flow of the metric, a definition of the scaling field at each scale $\lambda$ was obtained, further constrained by (52) and continuity in $\lambda$. A generalization of (32) to curved spacetime follows from adding the r.h.s. of (28) (with $b=-\frac{1}{2}$) to the r.h.s. of (32),

$${\nabla }^{\mu }{\nabla }_{\mu }{Z}_{\nu }+{\nabla }^{\mu }{\nabla }_{\nu }{Z}_{\mu }=-4\pi G{\ell }_g^2{F}_{\nu \mu }{J}^{\mu }+{\nabla }^{\mu }{\partial }_s{g}_{\mu \nu }$$

(33)

and using a covariant form for the boundary condition

$$P_{\mu \nu }\to -{\displaystyle \frac{1}{2\pi G{\ell }_g^2}}{g}_{\mu \nu }\mathrm{far}\mathrm{away}\mathrm{from}\mathrm{matter}$$

(34)

3.1.2. Time Independent Fixed-Point $A_{\mu }$

Now that $A_{\mu }$ depends on $Z^{\mu }$ via (16), and $Z^{\mu }$ on $A_{\mu }$ via (32), the nonlinear nature of the flow (16) can be appreciated, as well as the indirect mixing of $A_0$ and $\mathit{A}$ through $Z^{\mu }$, even in the time-independent case. Analyzing and solving the fixed-point solutions of the nonlinear system (16)(32) in the general, time-independent case, deserves a separate paper. Here we analyze the spherically symmetric, non-spinning ($\mathit{A}\equiv 0$) case in order to demonstrate how the nonlinearity introduces a second length scale, governing particle physics, which unlike ${\ell }_0$ is an attribute of the solution rather than a parameter of the model.

For $A_0\equiv \phi \left(r\right)$, $\mathit{A}\equiv 0$, the r.h.s. of (31) is an outwards-pointing radial force $x_{i}f\left(r\right)$ with

$$f\left(r\right):=\left(2{\displaystyle \frac{{\phi }^{\prime 2}}{r^2}}+{\displaystyle \frac{{\phi }^{\prime }{\phi }^{\prime \prime }}{r}}\right)$$

Plugging

$$Z_i=x_i+z_i\equiv x_i+x_{i}z\left(r\right),Z_0=x_0+h\left(r\right)$$

(35)

into (32) with boundary conditions $z(\infty )=0$, $z^{\prime }\left(0\right)=0$ (for the Laplacian to be well-defined at the origin) and $h^{\prime }\left(0\right)=h(\infty )=0$ translates into $h=0$ and a second order ODE for $z\left(r\right)$

$$2z^{\prime \prime }+{\displaystyle \frac{8z^{\prime }}{r}}=-4\pi G{\ell }_g^{2}f,z(\infty )=0,z^{\prime }\left(0\right)=0$$

(36)

Setting ${\partial }_s{A}_0=0$ in (16) gives

$${\ell }_0^2\left({\phi }^{\prime \prime }+{\displaystyle \frac{2{\phi }^{\prime }}{r}}\right)+r(1+z){\phi }^{\prime }+\phi =0,\phi (\infty )=0,{\phi }^{\prime }\left(0\right)=0$$

(37)

The system (36)(37) is symmetric under $A_{\mu }\mapsto -A_{\mu }$, $z^{\mu }\mapsto z^{\mu }$, guaranteeing that fixed-points come in particle-antiparticle pairs, which is true also in the general case. Since $z(\infty )=0\Rightarrow \phi (\infty )=0$, the system is under-determined, i.e., its solutions involve four integration constants satisfying only three independent conditions, therefore specifying a one-parameter family of solutions. Solutions of (36) are

$$z\left(r\right)=-2\pi G{\ell }_g^2{\int }_0^r\mathrm{d}{r}^{\prime }{\displaystyle \frac{{\int }_0^{r^{\prime }}\mathrm{d}{r}^{\prime \prime }{r^{\prime \prime }}^{4}f\left(r^{\prime \prime }\right)}{r^{\prime 4}}}+2\pi G{\ell }_g^2{\int }_0^{\infty }\mathrm{d}{r}^{\prime }{\displaystyle \frac{{\int }_0^{r^{\prime }}\mathrm{d}{r}^{\prime \prime }{r^{\prime \prime }}^{4}f\left(r^{\prime \prime }\right)}{r^{\prime 4}}}$$

(38)

A particle solution is defined by $J_0\approx 0$, hence also $f\left(r\right)\approx 0$, for $r\gtrsim r_{\mathrm{p}}$, where $r_{\mathrm{p}}$ is the particle’s radius (‘matter’). Inside matter $f\left(r\right)>0$, implying $z\left(r\right)>0$.

Using $f\left(0\right)=\frac{1}{3{\ell }_0^4}{\phi }^2\left(0\right)$, obtained from the analytic solution near $r=0$, the system (36)(37) can be numerically integrated from $r=0$ using $\phi \left(0\right)$ as a free parameter, adjusting $z\left(0\right)$ to meet $z(\infty )=0$, with the result that $r_{\mathrm{p}}$ monotonically and unboundedly decreases with increasing $\phi \left(0\right)$.

The Poincaré stress-energy tensor (27) reads

$$P_{ij}=-{\displaystyle \frac{1}{2\pi G{\ell }_g^2}}\left({\eta }_{ij}(1+z)+{\displaystyle \frac{x_i{x}_j}{r}}{z}^{\prime }\right),P_{00}={\displaystyle \frac{1}{2\pi G{\ell }_g^2}},P_{0i}=0$$

(39)

consisting of a z-independent ‘vacuum energy-momentum’ piece, $-{\left(2\pi G{\ell }_g^2\right)}^{-1}{\eta }_{\mu \nu }$, entirely due to the boundary condition (34), and a z-dependent, particle-specific piece, $p_{\mu \nu }$, which is colocalized with $J_0$. The trace of $p_{ij}$ is positive inside matter; Poincaré would have interpreted this as the negative pressure holding the particle against its internal Coulomb expansion, and the nontrivial scaling field, $z^i$, as the displacement vector due to the Coulomb stress. Since the positive ‘Coulomb pressure’ inside a charged particle is on the order of $q^2/r_{\mathrm{p}}^4$, equilibrium with the Poincaré pressure mandates a huge ${\partial }_i{z}^i\sim 2\pi q^{2}G{\ell }_g^2/r_{\mathrm{p}}^4\approx {10}^{32}$ for ${\ell }_g={10}^{20}\mathrm{km}$ and $r_{\mathrm{p}}$ a typical subatomic scale, e.g., the classical electron radius or the proton’s radius; Remarkably, all attributes of a particle involve gravitational parameters, G, ${\ell }_g$ and ${\ell }_0$. Asymptotically, $z\sim r^{-3}\Rightarrow p_{ii}=o\left(r^{-3}\right)$. Note that the $p_{00}$ component vanishes thus the particle’s energy is attributable entirely to its electrostatic self-energy ${\mathsf{\Theta }}_{00}$.

In conclusion of this section, a few final remarks. First, the above particle solution, although involving a nonlinearity, must not be conflated with soliton solutions of nonlinear PDEs, having a long history in modeling of particles. The existence of particles in the proposed formalism does not hinge on the flow being nonlinear (as seen in Section 2), but rather on a unique scaling operation countering the coarsener, therefore requiring a large scaling field inside a small particle. Second, at large r, the scaling piece completely dominates the fixed-point equation, exactly annihilating only the monopole, $\sim r^{-1}$, of $J_{\mu }$. No higher order multipoles are therefore part of a fixed-point solution at $r>{\ell }_0$ (it would take, e.g., $\alpha =2$ for the dipole). However, since electro- and magneto-statics are only testable at scales much smaller than ${\ell }_0$, higher order multipoles, required by a spinning particle, need only extend to distances$\ll {\ell }_0$, where the coarsener completely dominates, exactly annihilating any term in the multipole expansion. Second, charge quantization could be explained by cosmological considerations, of the type discussed at the end of Section 2.1. As in the linear case, and as seen in the spherically symmetric solution, fixed-points depend on a continuous parameter, controlling all their attributes this time, which is ‘spontaneously’ fixed at its observed value by a global consistency condition. Alternatively, charge quantization may arise just from the fixed-point condition for a spinning particle, where the reduced symmetry of the solution eliminates said free parameter. Finally, some/most/all real-world particles are likely represented by time-dependent solutions, which are fixed-points only in the statistical sense, when averaged over sufficiently long yet microscopic time intervals. Since $p_{00}\propto {\partial }_0{z}_0$, the time-averaged energy $\int {\mathrm{d}}^{3}x{p}_{00}$ would vanish nonetheless. The time dependence of such solutions would need to be chaotic, with a scale invariant power spectrum up to some cutoff frequency. Analyzing the properties of such dynamical fixed-points is suited for a statistical theory, complementing the proposed realistic model on such issues, which presumably is quantum mechanics and its generalizations; see Section 3.3

3.1.3. The Gravitational Field of a Matter Lump

Moving one step beyond the flat spacetime approximation, the distortion to ${\eta }_{\mu \nu }$ and $Z^{\mu }=x^{\mu }+z^{\mu }$ caused by a fixed-point solution is calculated. To this end the flat metric is replaced with a (quasi-)static metric

$$g_{00}=-(1-\mathsf{\Phi }),g_{ij}={\eta }_{ij}(1+\mathsf{\Phi }),g_{i0}=g_{0i}=0,\left|\mathsf{\Phi }\right|\ll 1$$

(40)

with non-vanishing Christoffel symbols

$${\mathsf{\Gamma }}_{i0}^0=-\frac{1}{2}{\partial }_i\mathsf{\Phi },{\mathsf{\Gamma }}_{00}^i=-\frac{1}{2}{\partial }^i\mathsf{\Phi },{\mathsf{\Gamma }}_{jk}^i=\frac{1}{2}\left({\partial }_j\mathsf{\Phi }{\eta }_{\phantom{1}k}^i+{\partial }_k\mathsf{\Phi }{\eta }_{\phantom{1}j}^i-{\partial }^i\mathsf{\Phi }{\eta }_{jk}\right)$$

(41)

and we seek a fixed-point solution of the metric flow (24), with $T_{\mu \nu }:={\mathsf{\Theta }}_{\mu \nu }+P_{\mu \nu }$ now calculated using (40) instead of ${\eta }_{\mu \nu }$, and $Z^{\mu }$ solving (33). Starting with the latter, the non gravitational part on the r.h.s. of (33), involving $J_{\mu }$, vanishes outside matter and is ignored. The remaining, gravitational equation, has on its r.h.s. a source, which for a fixed-point reads:

$$\lambda {\partial }^i{\partial }_{\lambda }{\mathsf{\Phi }}_{\mathrm{p}}(\mathit{x}-{\lambda }^{-1}{\mathit{x}}^{\prime })={\lambda }^{-1}{x^{\prime }}^j{\partial }^i{\partial }_j{\mathsf{\Phi }}_{\mathrm{p}}(\mathit{x}-{\lambda }^{-1}{\mathit{x}}^{\prime })$$

for some offset ${\mathit{x}}^{\prime }$, where ${\mathsf{\Phi }}_{\mathrm{p}}$ is the metric perturbation adjunct to particle/fixed-point p. In the above expression it was assumed that during scale flow, naive scaling alone rigidly shifts the center of an offset fixed-point, which turns out below to be a consistent assumption. Using the metric (40) and an ansatz in (33): $Z^i=x^i(1+{\mathsf{\Phi }}_{\mathrm{p}}\left(\mathit{x}\right))+u^i\left(\xi \right)$, $Z^0=x^0$, with $\xi :=\mathit{x}-{\lambda }^{-1}{\mathit{x}}^{\prime }$, translation invariant equations for $u^i$ are obtained,

$${\partial }^j{\partial }_j{u}^i\left(\xi \right)+{\partial }^i{\partial }_j{u}^j\left(\xi \right)=-{\xi }^j{\partial }^i{\partial }_j{\mathsf{\Phi }}_{\mathrm{p}}\left(\xi \right)$$

with solutions $\left|\mathit{u}\right|=\left|\xi \right|O\left(\mathsf{\Phi }\right)$. This $u^i$ component of $Z^i$, adjunct to a body, is the gravitational counterpart of $z^i$ from the previous section. Note that, as with $z^i$, for multiple fixed-points at close proximity, $Z^i=x^i\left(1+{\sum }_{\mathrm{p}}{\mathsf{\Phi }}_{\mathrm{p}}\left(\mathit{x}\right)\right)+{\sum }_{\mathrm{p}}{u}_{\mathrm{p}}^i\left({\xi }_{\mathrm{p}}\right)$ is not a solution (individual ${\mathit{u}}_{\mathrm{p}}$ must get deformed for that). However, such $O\left(\mathsf{\Phi }\right)$ departures of $Z^i$ from naive scaling only introduce $O\left(\mathsf{\Phi }\right)$ corrections to $O\left(1\right)$ quantities, and are neglected henceforth.

Next, using (40) for the metric in (24), setting ${\partial }_s{g}_{\mu \nu }=0$, i.e., a fixed-point stationed at the origin of scaling for simplicity (offset ${\mathit{x}}^{\prime }=0$) and temporarily dropping the ‘p’ label from $\mathsf{\Phi }$ gives

$${\nabla }^2\mathsf{\Phi }=4\pi G\left({\mathsf{\Theta }}_{00}+p_{00}^{\mathrm{met}}\right)+{\displaystyle \frac{2}{{\ell }_g^2}}$$

(42)

with $p_{00}^{\mathrm{met}}$ being the metric part of $p_{00}$, resulting from ${\partial }_0{Z}_0\mapsto {\nabla }_0{Z}_0$ in (27) (the matter part vanishes; see previous section), and the vacuum piece (34) of $P_{\mu \nu }$, now reading $P_{00}^{\mathrm{vac}}=(1-\mathsf{\Phi })/\left(2\pi G{\ell }_g^2\right)\approx 1/\left(2\pi G{\ell }_g^2\right)$. Since in Newtonian gravity, a uniform background density can have no consistent, non-vanishing effect, the tiny $P_{00}^{\mathrm{vac}}$ is ignored in this section, reinterpreted in the context of a cosmological model, Section 3.4.

The scaling piece coming from the Lie derivative, explicitly appearing in the fixed-point equation (37) for $A_{\mu }$, appears to be missing from (42). This is because it has been absorbed into $p_{00}^{\mathrm{met}}$:

$$p_{00}^{\mathrm{met}}\approx -{\displaystyle \frac{1}{4\pi G{\ell }_g^2}}\left(x^i+z^i\right){\partial }_i\mathsf{\Phi }$$

(43)

Moving $p_{00}^{\mathrm{met}}$ to the l.h.s. of (42), using $z^i\ll x^i$ outside matter from the previous section, we obtain the fixed-point equation for an inhomogeneous flow of a scalar $\mathsf{\Phi }$

$${\ell }_g^2\left({\nabla }^2\mathsf{\Phi }-4\pi G{\mathsf{\Theta }}_{00}\right)+x^i{\partial }_i\mathsf{\Phi }=0$$

(44)

with the $z^i{\partial }_i\mathsf{\Phi }$ term inside matter omitted by our assumption that gravity plays a negligible role in the structure of matter. At distances from the source smaller than ${\ell }_g$, solutions of (44) are approximately the usual solutions of Newtonian gravity, ${\mathsf{\Phi }}_m=-Gm/\left|\mathit{x}\right|$, with $m=\int {\mathrm{d}}^{3}x{\mathsf{\Theta }}_{00}$. A covariant generalization of (42) would be

$$R_{\mu \nu }-\frac{1}{2}{g}_{\mu \nu }R=4\pi G\left({\mathsf{\Theta }}_{\mu \nu }+P_{\mu \nu }\right)$$

(45)

which is just Einstein’s fields for $G\mapsto -\frac{1}{2}G$. The physical meaning of $g_{\mu \nu }$ solving (45) is only revealed through its effect on $A_{\mu }$, and in the following sections the reason for this peculiar coefficient of G becomes apparent. Note that $P_{\mu \nu }$ now involves solutions to the wave equation (33) and, consequently, the uniform vacuum energy in (42) becomes nonuniform and time-dependent (having ‘vacuum ripples’) possibly having an observable effect in strong fields.

In most realistic scenarios, ${\mathsf{\Theta }}_{\mu \nu }$ is time-dependent, focused on the worldline, $\zeta (t,\lambda )$, of a particle. The flow of the adjunct metric perturbation solving (24) then involves a form more general then (40): $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, $\left|h_{\mu \nu }\right|\ll 1$, with $h_{\mu \nu }$ also time-dependent. However, for non relativistic $\zeta (t,\lambda )$ the flow of $\mathsf{\Phi }\equiv h_{00}$ takes the simple form

$$\lambda {\partial }_{\lambda }\mathsf{\Phi }={\ell }_g^2\left({\nabla }^2\mathsf{\Phi }-{\partial }_{tt}\mathsf{\Phi }-4\pi G{\mathsf{\Theta }}_{00}\right)+x^i{\partial }_i\mathsf{\Phi }+t{\partial }_t\mathsf{\Phi }$$

(46)

with the negligible $-{\partial }_{tt}\mathsf{\Phi }$ term above retained for the pedagogical reason: A solution for a uniformly moving $\zeta$ can then be obtained by simply boosting the previous static perturbation as in Section 2.2 (though only to nonrelativistic velocities for consistency). For a non-uniform $\zeta$, we first note the freedom to choose the scaling center at will: If $\mathsf{\Phi }(\mathit{x},t,\lambda )$ solves (46) with ${\mathsf{\Theta }}_{00}$ focused on $\zeta (t,\lambda )$, then so does $\mathsf{\Phi }(\mathit{x}-{\lambda }^{-1}{\mathit{x}}^{\prime },t-{\lambda }^{-1}{t}^{\prime },\lambda )$ with $\zeta \mapsto \zeta (t-{\lambda }^{-1}{t}^{\prime },\lambda )+{\lambda }^{-1}{\mathit{x}}^{\prime }$ for any ${\mathit{x}}^{\prime },t^{\prime }$. Using this freedom we can focus on any point on the worldline traced by $\zeta$, assuming $t=0$. It is then readily verified that at sufficiently small accelerations, $\ddot{\zeta }$ (e.g., for a path $\zeta :={\lambda }^{-1}{\zeta }_{\mathrm{scl}}\left(\lambda t\right)$ in its scaling regime), ${\mathsf{\Phi }}_m(\mathit{x}-\zeta )$ becomes a pointwise, arbitrarily good solution for (46) in a neighborhood,

$$\left|\mathit{x}-\zeta \right|\le c^2{\left|\ddot{\zeta }\right|}^{-1}\ll {\ell }_g$$

(47)

Since outside this neighborhood ${\mathsf{\Phi }}_m\approx 0$, the standard, Newtonian approximation

$$\mathsf{\Phi }(\mathit{x},t,\lambda )\approx \sum _{\mathrm{p}}{\mathsf{\Phi }}_{m_{\mathrm{p}}}\left(\mathit{x}-{\zeta }_k(t,\lambda )\right)$$

(48)

where ‘p’ is a particle label, becomes an arbitrarily good, point-wise, global solution of (46). The time and scale dependence of $\mathsf{\Phi }$ is therefore entirely inherited from that of $\left\{{\zeta }_k\right\}$.

As $\mathsf{\Phi }$ (48) flows to larger $\lambda$ the approximation (47) it involves gradually breaks due to the shrinkage of that neighborhood of $\zeta$ at large accelerations. Moreover, for relativistic $\zeta$, the flow equation for components other than $h_{00}$ are no longer satisfied automatically, and (46) is to be replaced with

$$\lambda {\partial }_{\lambda }{h}_{\mu \nu }={\ell }_0^2\left(G_{\mu \nu }^{\left(1\right)}-4\pi G\left({\mathsf{\Theta }}_{\mu \nu }+P_{\mu \nu }\right)\right)$$

where $G_{\mu \nu }^{\left(1\right)}$ is the usual linearized Einstein tensor. The linearized EFE must then be satisfied in the ${\ell }_0$-neighborhood of $\zeta$ for $h_{\mu \nu }$ to not rapidly diverge at large $\lambda$.

3.2. The Motion of Matter Lumps in a Weak Gravitational Field

Equation (48) prescribes the scale flow of the metric (in the Newtonian approximation), given the a set, $\left\{{\zeta }_k(t,\lambda )\right\}$, of worldlines associated with matter lumps. To determine this set, an equation for each ${\zeta }_k(t,\lambda )$, given $\mathsf{\Phi }(t,\mathit{x},\lambda )$, is obtained in this section. This is done by analyzing the scale flow of the first moment of $J_0$ associated with a general matter lump, using (22). An obstacle to doing so comes from the fact that, $\zeta (t,\lambda )$ now incorporates both gravitational and non gravitational interactions in a convoluted way, as the existence of gravitating matter depends on it being composed of charged matter. In order to isolate the effect of gravity on $\zeta (t,\lambda )$, we first analyze the motion of a body in the absence of gravity, i.e., $g_{\mu \nu }\equiv {\eta }_{\mu \nu }$, $Z^{\mu }=x^{\mu }+z^{\mu }$ with $x^{\mu }$ the Minkowskian coordinates and $z^{\mu }$ the structural component from the previous sections. To this end, a better understanding of the scale flow (22) of $J_{\mu }$ is needed. Using (17) and (19) plus some algebra the scaling piece in (22) reads

$${\widehat{W}}_{\nu }^{\phantom{1}\mu }{\mathcal{L}}_Z{A}_{\mu }=x^{\mu }{\partial }_{\mu }{J}_{\nu }+3J_{\nu }+M_{\nu }$$

(49)

with

$$M_{\nu }:={\partial }_{\mu }{M}_{\phantom{1}\nu }^{\mu }:={\partial }_{\mu }\left(z^{\rho }{\partial }_{\rho }{F_{\nu }}^{\mu }+\left({\partial }^{\mu }{z}^{\rho }\right)F_{\nu \rho }+\left({\partial }_{\nu }{z}^{\rho }\right){F_{\rho }}^{\mu }\right)$$

(50)

(where $M_{\phantom{1}\nu }^{\mu }$ is of course only defined up to a divergence free piece). The first two terms in (49) are the familiar $\alpha =1\mapsto \alpha =3$ conversion, to which a ‘matter vector’, $M_{\nu }$, is added. In the absence of gravity the commutator on the l.h.s. of (22) vanishes. Combined, we get

$${\partial }_s{J}_{\nu }={\ell }_0^2{\square }^2{J}_{\nu }+x^{\mu }{\partial }_{\mu }{J}_{\nu }+3J_{\nu }+M_{\nu }$$

(51)

Since ${\partial }^{\nu }{J}_{\nu }\equiv 0$, taking the divergence of (51) implies

$${\partial }^{\nu }{M}_{\nu }\equiv 2\left({\square }^2{z}^{\rho }\right)J_{\rho }+2\left({\partial }^{\alpha }{z}^{\rho }\right)\left({\partial }_{\alpha }{J}_{\rho }+{\partial }_{\rho }{J}_{\alpha }\right)\equiv 2\left({\square }^2{z}^{\rho }\right)J_{\rho }-2\left({\partial }^{\alpha }{z}^{\rho }\right){\square }^2{F}_{\alpha \rho }=0$$

(52)

where the second equality follows from the Maxwell-Bianchi identity. This constraint, relating $A_{\mu }$ to $z^{\mu }$ inside matter, is the counterpart of (32), derived this time from the flow of $A_{\mu }$ rather than that of $g_{\mu \nu }$, and from local charge conservation rather than local energy-momentum conservation; Note that our particle solution from Section 3.1.2 automatically satisfies (52) by virtue of $z^0=J_i=0$. In the case of a general $g_{\mu \nu }$, operating with ${\nabla }^{\nu }$ on (22), the l.h.s vanishes identically as it is just the covariant divergence of $J_{\nu }$ (19) with $A_{\mu }\mapsto {\partial }_s{A}_{\mu }$ there. Using (20) the r.h.s gives a covariant counterparts of (52) involving curvature terms, but still expressing a section constraint.

Defining $q\left(s\right):=\int {\mathrm{d}}^{3}x{J}_0$ and $M\left(s\right):=\int {\mathrm{d}}^{3}x{M}_0$ the (conserved in time-) electric and ‘matter’ charges resp., and integrating (51) over three-space implies ${\displaystyle \frac{\mathrm{d}}{\mathrm{d}s}}q=M$ i.e., electric charge is conserved in scale if and only if the matter charge vanishes. That the latter is identically true follows from the divergence form (50) of $M_{\nu }$ and $M_{\phantom{1}0}^0\equiv 0$. Whether this result generalizes to curved spacetime is unclear to the author at this point4. Clearly, since charge is conserve in time, ${\displaystyle \frac{\mathrm{d}}{\mathrm{d}s}}q\approx 0$ whenever spacetime is approximately locally flat around some point on the world-line of a charge, but it is possible that a center more complicated than $J_{\nu }$ (19) is required for identical charge conservation in scale. Alternatively, charge may not be identically conserve in scale when strong gravitational fields are involved, or it could emerge as a necessary condition for inclusion in $\mathcal{S}$, as in the linear model.

Next, multiplying (51) by $x^i/q$ and integrating over a ball, B, containing a body of charge q, results in

$$\lambda {\partial }_{\lambda }{\zeta }^i=-{\ell }_0^2{\partial }_{tt}{\zeta }^i+t{\partial }_t{\zeta }^i-{\zeta }^i+q^{-1}{\int }_B{\mathrm{d}}^{3}x{x}^i{M}_0(\mathit{x},t,s)$$

(53)

where ${\zeta }^i(t,s)=q^{-1}{\int }_B{\mathrm{d}}^{3}x{x}^i{J}_0(\mathit{x},t,s)$ is an object’s ‘center-of-charge’ (c.o.c.). Above and in the rest of this section, the charge of a body, assumed nonzero for simplicity, is only used as a convenient tracer of matter. Now, as $J_0$ can be both positive and negative, the c.o.c. is not necessarily confined to the support of $J_0$, as with positive distributions. Nonetheless, since $J_0\left(\mathit{x}\right)\mapsto J_0(\mathit{x}-{\mathit{x}}^{\prime })\Rightarrow \zeta \mapsto \zeta +{\mathit{x}}^{\prime }$, $\zeta$ does follow the particle up to some constant displacement, reflecting to a large extent the arbitrariness in defining the exact position of an extended body, and assumed much smaller than any competing length.

At the (sub-)atomic level, $J_{\mu }$ would be rapidly fluctuating in time, endowing $\zeta$ with a ‘jitter motion’ component. To remove it, (53) is convolved with a normalized symmetric kernel of macroscopic extent T. Defining

$$\overline{\zeta }\left(t\right):=\int \mathrm{d}{t}^{\prime }{K}_T(t-t^{\prime })\zeta \left(t^{\prime }\right)$$

the form of (53) is retained for $\zeta \left(t\right)\mapsto \overline{\zeta }\left(t\right)$ (with $M_0\mapsto {\overline{M}}_0$), except for a correction

$$-q^{-1}\int (t-t^{\prime })K_T(t-t^{\prime }){\partial }_{t^{\prime }}{\zeta }^i\left(t^{\prime }\right)\mathrm{d}{t}^{\prime }$$

coming from the time-scaling term. If ${\partial }_t\overline{\zeta }$ is slowly varying over T, this correction is at most on the order of $T^2{\partial }_{tt}\overline{\zeta }$, negligibly renormalizing the ${\ell }_0^2{\partial }_{tt}\overline{\zeta }$ term for $T\ll {\ell }_0$. For economical notations, then, the ‘bar’ is dropped henceforth from all quantities.

The integral in (53) is the first moment of a distribution, $M_0$, whose zeroth moment vanishes, which is therefore invariant under $M_0\left(\mathit{x}\right)\mapsto M_0(\mathit{x}-{\mathit{x}}^{\prime })$—as is expected of a ‘force term’, competing with the ${\partial }_{tt}{\zeta }^i$, acceleration term. Buried in it are presumably all forms of non gravitational interactions preventing ${\zeta }_{\mathrm{free}}:={\lambda }^{-1}{\mathit{x}}^{\prime }+(t-{\lambda }^{-1}{t}^{\prime })\mathit{v}$ for some constant $(t^{\prime },{\mathit{x}}^{\prime })$ and $\mathit{v}$, from being a solution of (53) in the absence of gravity (In particular, elements of the Lorentz force can explicitly be excavated from it, but this route shall not be explore here because the Lorentz force is implicit in the basic tenets of classical electrodynamics and, moreover, as there is only a single, global A-field, its decomposition into $J_{\mu }$, an adjunct ‘self-field’, and an external field, is far from being obvious). Accordingly, this term is ignored for an isolated body, for which all forces are internal and necessarily integrate to zero (this can be verified explicitly, e.g., for a spherically symmetric $M_{\phantom{1}0}^i$). Thus (53) becomes (14), and as proved in that case, solutions for ${\zeta }^i$ must all be straight, non tachyonic worldlines ${\zeta }_{\mathrm{free}}$.

The effect of gravity on those straight worldlines is derived by including a weak field in the flow of the first-moment projection of (22). Recalling from Section 3.1.2 that weak gravity has no effect on the scaling field, and since gravity is assumed to play a negligible role in the structure of matter, the way this field enters the flat spacetime analysis is by ‘dotting the commas’ in partial derivatives. Before analyzing the first moment projection we first have to show that the identically conserved-in-time charge by virtue of the last identity in (19),

$$\int {\mathrm{d}}^{3}x\sqrt{-g}{J}_0(\mathit{x},t,\lambda )=q\left(s\right)$$

(54)

is also conserved in scale. As previously remarked, this might automatically follow from the previously mentioned covariant counterpart of (52), but since that more complicated constraint had not been explored, and as a preparation for the first-moment analysis, we explicitly show next that (54) hold up to $O\left(\mathsf{\Phi }\right)$ corrections. To this end (22) is evaluated using the Newtonian metric (40), then multiplied by $\sqrt{-g}\approx 1+\mathsf{\Phi }$ and integrated over three-space. A straightforward calculation to first order in $\mathsf{\Phi }$, incorporating ${\nabla }^2\mathsf{\Phi }=0$, gives

$${\square }_{\mathrm{GR}}^2{J}_0\equiv {\nabla }^{\alpha }{\nabla }_{\alpha }{J}_0\approx -(1+\mathsf{\Phi }){\partial }_{tt}{J}_0+(1-\mathsf{\Phi }){\partial }_j\left({\partial }_j{J}_0+{\partial }_j\mathsf{\Phi }{J}_0\right)-{\partial }_j\mathsf{\Phi }{\partial }_t{J}_j$$

(55)

Ignoring $O\left(\mathsf{\Phi }\right)$ corrections to the isotropic coarsener, the net effect of of the potential in the Newtonian approximation is to render the coarsener anisotropic through its gradient, with an added relativistic correction in the form of the last term on the r.h.s. of (55). The reader can then verify that the contribution of (55) vanishes up to $O\left(\mathsf{\Phi }\right)$ corrections (see also below, the corresponding contribution to the first moment projection). The same is true for the contribution of the scaling piece, since $Z^{\mu }\approx x^{\mu }$ for weak fields and ${\nabla }_{\mu }{J}_0\approx {\partial }_{\mu }\left(J_0(1+O\left(\mathsf{\Phi }\right))\right)$. As for the contribution of the non-vanishing commutator on the l.h.s. of (22), using $F_{\mu \nu }\equiv {\partial }_{\mu }{A}_{\nu }-{\partial }_{\nu }{A}_{\mu }$ (viz., ordinary derivatives can replace covariant ones) in the definition (19) of $J_{\mu }$, we have

$$\left[{\widehat{W}}_0^{\phantom{1}\mu },{\partial }_s\right]A_{\mu }={\partial }_s\mathsf{\Phi }{J}_0+{\partial }_s{\partial }^i\mathsf{\Phi }{F}_{i0}$$

(56)

Following integration by parts of the second term, the contribution of the commutator, combined with that of ${\partial }_s{J}_0$, reads ${\partial }_s\left(q(1+O(\mathsf{\Phi })\right)$. Finally, since $R_{\mu \nu }\propto {\nabla }^2\mathsf{\Phi }$, the contribution of the curvature term vanishes, as is that of the covariant $M_0$ term up to $O\left(\mathsf{\Phi }\right)$ corrections.

Moving to the first-moment projection of (22), we multiply it by $\sqrt{-g}{x}^i/q\approx (1+\mathsf{\Phi })x^i/q$, plug expression (55) for the d’Alembertian and integrate over B, assuming ${\partial }_i\mathsf{\Phi }$ is approximately constant over the extent of the body. At non-relativistic velocities and neglecting $O\left(\mathsf{\Phi }\right)$ corrections, the double time-derivative piece contributes $-{\partial }_{tt}{\zeta }^i$, where ${\zeta }^i=q^{-1}\int {\mathrm{d}}^{3}x{x}^i\sqrt{-g}{J}_0$. In the second term on the r.h.s. of (55) the $(1+\mathsf{\Phi })$ cancels (to first order in $\mathsf{\Phi }$) the $(1-\mathsf{\Phi })$ factor multiplying it, which is then integrated by parts. Neglecting $O\left(\mathsf{\Phi }\right)$ corrections the result is $-{\partial }^i\mathsf{\Phi }(\zeta ,t,\lambda )$. Using the continuity equation for $J_{\mu }$, integration by parts of the last term in (55) yields a negligible $O\left(\mathsf{\Phi }\right)$, vanishing at non-relativistic velocities. As with the zeroth-moment project, the modification to the scaling piece (49) only introduces an $O\left(\mathsf{\Phi }\right)$ correction to the $-{\zeta }^i$ term in (53) which is neglected in the Newtonian approximation. Finally, the contribution ${\partial }_s{J}_0$, combined with the commutator (56), gives on the l.h.s.: ${\partial }_s\left({\zeta }^i(1+O\left(\mathsf{\Phi }\right)\right)-q^{-1}\int {\mathrm{d}}^{3}x{\partial }_s\mathsf{\Phi }{F}_{i0}$. This last integral is evaluated by multiplying ${\partial }^j{F}_{j0}\approx -J_0$ by $x^i$ and integrating by parts the l.h.s., which only modifies the $O\left(\mathsf{\Phi }\right)$ in the first piece.

Combining all pieces, the first moment projection of (22) reads

$$\lambda {\partial }_{\lambda }{\zeta }^i={\ell }_0^2\left(-{\partial }_{tt}{\zeta }^i-{\partial }^i\mathsf{\Phi }(\zeta ,t,\lambda )\right)-{\zeta }^i+t{\partial }_t{\zeta }^i$$

(57)

This equation is just (14) with an extra ‘force-term’ on its r.h.s. which could salvage a non uniformly moving solution, ${\zeta }^i$, from the catastrophic fate at $\lambda \to \infty$ suffered by its linear counterpart.

At sufficiently large scales, $\lambda$, when all relevant masses contributing to $\mathsf{\Phi }$ occupy a small ball of radius $\ll {\ell }_0$ centered at the origin of scaling without loss of generality ($t\ll {\ell }_0$ can similarly be assumed) the scaling part on the r.h.s of (57) becomes negligible compared with the force term, rather benefiting from such crowdedness. It follows that each $\left|\zeta \right|$ would grow—extremely rapidly as we show next—with increasing $\lambda$ even when the weak-field approximation is still valid, implying that the underlying $A_{\mu }$ is not in $\mathcal{S}$. The only way to keep the scale evolution of $\zeta$ under control is for the acceleration term to similarly grow, almost canceling the force term but not quite, which is critically important; it is the fact that the sum of these two terms, both originating from the coarsener, remains on the order of the scaling term, which is responsible for a nontrivial, non pure scaling $\zeta$. This means that each worldline converges at large scales to that satisfying Newton’s equation

$${\partial }_{tt}{\zeta }^i=-{\partial }^i\mathsf{\Phi }(\zeta ,t,\lambda )$$

(58)

At small scales the opposite is true. The scaling part dominates and any scaling path, i.e., ${\zeta }^i(t,\lambda )={\lambda }^{-1}{\zeta }_{\mathrm{scl}}^i\left(\lambda t\right)$ is well behaved. Combined: at large scales ${\zeta }^i$ is determined, then simply scaled at small scales, gradually converging locally to a freely moving particle. Finally, (58) must be true also for a loosely bound system, e.g. a wide binary moving in a strong external field, implying that Newton’s law applies also to their relative vector, not merely to their c.o.m.

Deriving a manifestly covariant generalization of (57) is certainly a worthwhile exercise. However, in a weak field the result could only be

$$\lambda {\partial }_{\lambda }{\zeta }^{\mu }={\ell }_0^2\left(-{\partial }_{\tau \tau }{\zeta }^{\mu }-{\overline{\mathsf{\Gamma }}}_{\nu \rho }^{\mu }{\partial }_{\tau }{\zeta }^{\nu }{\partial }_{\tau }{\zeta }^{\rho }\right)+\tau {\partial }_{\tau }{\zeta }^{\mu }-z^{\mu }$$

(59)

with $\tau$ some scalar parameterization of the worldline traced by ${\zeta }^{\mu }$, and $z$ the gravitational part of the scaling field, well-approximated by $x^{\mu }$. Above, ${\overline{\mathsf{\Gamma }}}_{\nu \rho }^{\mu }$ are the Christoffel symbols associated with $g_{\mu \nu }(x,\lambda ;-2G)$, i.e., the analytic continuation of the metric, seen as a function of Newton’s constant, to $-2G$. Recalling from Section 3.1.3 that the fixed-point $g_{\mu \nu }$ is a solution of the standard EFE analytically continued to $G\mapsto -\frac{1}{2}G$, $\overline{\mathsf{\Gamma }}$ in that case is therefore just the Christoffel symbol associated with standard solutions of EFE’s. The previous, Newtonian approximation is a private case of this, where $\mathsf{\Phi }$ contains a factor of G. Note however that the path of a particle in our model is a covariantly defined object irrespective of the analytic properties of $g_{\mu \nu }$. Resorting to analyticity simply provides a constructive tool for finding such paths whenever $g_{\mu \nu }$ is analytic in G. In such cases, the covariant counterpart of (58) becomes the standard geodesic equation of GR which gives great confidence that this is also the case for non-analytic $g_{\mu \nu }$.

The reasons for trusting (59) are the following. It is manifestly scale- and general-covarinat, as is our model; it is $\tau$-shift invariant, i.e., ${\zeta }^{\mu }(\tau -{\lambda }^{-1}{\tau }^{\prime },\lambda )$, parameterizing the same, scale dependent world-line, also solves (59) for any ${\tau }^{\prime }$; At nonrelativistic velocities in a Minkowskian background, ${\zeta }^t\equiv t\approx \tau$ solves (59) which, when substituted into the i-components of (59), recovers (57); The scaling regime ansatz, ${\zeta }^{\mu }(\tau ,\lambda )\sim {\zeta }_{\mathrm{scl}}^{\mu }\left(\lambda \tau ,\lambda \right)$, solves $\lambda {\partial }_{\lambda }{\zeta }_{\mathrm{scl}}^{\mu }=-z^{\mu }({\zeta }_{\mathrm{scl}},\lambda )$, i.e., each point on the world-line traced by ${\zeta }^{\mu }$, indexed by a fixed $\lambda \tau$, flows along integral curves of the scaling field—as must be the case when the coarsener is negligible; It only involves local properties of ${\zeta }^{\mu }$ and $g_{\mu \nu }$, i.e., their first two derivatives, which must also be a property of a covariant derivation, as is elucidated by the non-relativistic case. Thus (59) is the only candidate up to covariant, higher derivatives terms involving ${\zeta }^{\mu }$ and $g_{\mu \nu }$, or nonlinear terms in their first or second derivatives, all becoming negligible in weak fields/ at small accelerations.

In conclusion of this section we wish to relate (59) to the fact that $T_{\mu \nu }$ is not covariantly conserved by virtue of (28). It was well known already to Einstein that the geodesic equation follows from local energy-momentum conservation under reasonable assumptions. Similarly, the Lorentz force equation follows quite generally from the basic tenets of classical electrodynamics (31),(30),(19). Referring to Figure 1, both results are derived by integrating ${\nabla }_{\mu }{T}^{\mu \nu }=0$ over a world-cylinder, C, with the ${\partial }_{\tau \tau }{\zeta }^{\nu }$ term obtained by converting part of the volume integral into a surface integral over the $\Sigma$’s via (a relativistic generalization of-) Stoke’s theorem, leaving the remaining part for the ‘force term’, which gives ${\overline{\mathsf{\Gamma }}}_{\nu \rho }^{\mu }{\partial }_{\tau }{\zeta }^{\nu }{\partial }_{\tau }{\zeta }^{\rho }$ in the case of gravity. The integral over T represents small radiative corrections to the geodesic/Lorentz-force equation which can be ignored in what follows. A crucial point in that derivation is that the result is insensitive to the form of C so long as the $\Sigma$’s contain the support of the particle’s energy-momentum distribution. This insensitivity would not carry to the r.h.s. of (28) should it be transferred sides, whether or not $g_{\mu \nu }$ is to include the self-field of the particle. Thus attempting to generalize the geodesic equation based on (28), in the hope that it would reproduce fixed-$\lambda$ sections of (59) solutions, is bound to fail. Nonetheless, since the conservation-violating r.h.s. of (28) is a tiny $O\left({\ell }_g^{-2}\right)$, it can consistently be ignored for any reasonable choice of C whenever the individual terms in the geodesic equation are much larger, i.e., at large scales; it is only at small scales that each term becomes comparable to the r.h.s. even for reasonable choices of C. And since in the limit ${\ell }_0\to \infty$ paths become simple geodesics at any scale according to (59), the two length parameters, ${\ell }_0$, ${\ell }_g$ are not independent, constrained by consistency of solutions such as the above, and likely additional ones involving the structure of matter, as hinted to by the fixed-point example in Section 3.1.2.

3.2.1. Application: The Rotation Curve of Disc Galaxies

As a simple application of (57), let us calculate the rotation curve, $v\left(r\right)$, of a scale-invariant mass, M, located at the origin, as it appears to an astronomer of native scale $\lambda =1$. Above, r is the distance to the origin of a test mass orbiting M in circles at velocity v. Since $\mathsf{\Phi }$ in (57) is time-independent, the time-dependence of ${\zeta }^i$ can only be through the combination $\mathsf{\Omega }\left(\lambda \right)t$ for some function $\mathsf{\Omega }$. Looking for a circular motion solution in the $x^1-x^2$ plane,

$${\zeta }^1(\lambda ,t)=r\left(\lambda \right)sin\left(\mathsf{\Omega }\left(\lambda \right)t\right),{\zeta }^2(\lambda ,t)=r\left(\lambda \right)cos\left(\mathsf{\Omega }\left(\lambda \right)t\right)$$

and equating coefficients of $sin\left(\mathsf{\Omega }t\right)$ and $cos\left(\mathsf{\Omega }t\right)$ for each component, the system (57) reduces to two, first order ODE’s for $\mathsf{\Omega }\left(\lambda \right)$ and $r\left(\lambda \right)$. The equation for $\mathsf{\Omega }$ readily integrates to $\mathsf{\Omega }=\omega \lambda$ for some integration constant $\omega$, and for r it reads

$$\lambda {\displaystyle \frac{\mathrm{d}r}{\mathrm{d}\lambda }}={\ell }_0^2\left({\lambda }^2{\omega }^{2}r-{\displaystyle \frac{GM}{r^2}}\right)-r$$

(60)

Solutions of (60) with $r\left(1\right)\equiv r_1$ as initial condition, all diverge in magnitude for $\lambda \to \infty$ except for a single value of $\omega$ for which $r\left(\lambda \right)\to 0$ in that limit; for any other ${\omega }^{\prime }\gtrless \omega$, $r\left(\lambda \right)$ rapidly diverges to $\pm \infty$ respectively; the map $r_1\mapsto \omega$ is invertible. We note in advance that, for a mistuned $\omega$ this divergence starts well before the weak field approximation breaks down due to the $r^{-2}$ term, and neglected relativistic and self-force terms become important, and being so rapid, $\sim e^{{\lambda }^2}$, those would not tame a rogue solution. It follows that there is no need to complicate our hitherto simple analysis in order to conclude that ${\omega }^{\prime }=\omega$ is a necessary conditions for r to correspond to an $A_{\mu }\in \mathcal{S}$.

Solutions of (60) which are well-behaved for $\lambda \to \infty$ admit a relatively simple analytic form. Reinstating c and defining $\overline{\lambda }:=\omega {\ell }_0\lambda /c$ the result is

$$r\left(\overline{\lambda }\right):={\left({\displaystyle \frac{GM{\ell }_0^2}{6c^2}}\right)}^{1/3}{\displaystyle \frac{1}{\overline{\lambda }}}{\left[6\overline{\lambda }-\sqrt{6\pi }{e}^{{\displaystyle \frac{3}{2}}{\overline{\lambda }}^2}\left(\mathrm{erf}\left(\sqrt{3/2}\overline{\lambda }\right)-1\right)\right]}^{1/3},$$

(61)

having the following power law asymptotic forms

$$r\left(\overline{\lambda }\right)\sim {\left({\displaystyle \frac{GM{\ell }_0^2}{c^2}}\right)}^{1/3}\times \left\{\begin{array}{cc}{\overline{\lambda }}^{-2/3}\hfill & \mathrm{for}\overline{\lambda }\gg 1,\text{‘}\mathrm{coarsening}\mathrm{regime}\text{’}\hfill \\ {\left({\displaystyle \frac{\pi }{6}}\right)}^{1/6}{\overline{\lambda }}^{-1}\hfill & \mathrm{for}\overline{\lambda }\ll 1,\text{‘}\mathrm{scaling}\mathrm{regime}\text{’}\hfill \end{array}\right.$$

(62)

with the corresponding asymptotic circular velocity, $v=\left(\omega \lambda \right)r$

$$v\left(\overline{\lambda }\right)\sim {\left({\displaystyle \frac{GMc}{{\ell }_0}}\right)}^{1/3}\times \left\{\begin{array}{cc}{\overline{\lambda }}^{1/3}\hfill & \mathrm{for}\overline{\lambda }\gg 1,\hfill \\ {\left({\displaystyle \frac{\pi }{6}}\right)}^{1/6}\hfill & \mathrm{for}\overline{\lambda }\ll 1\hfill \end{array}\right.$$

(63)

With these asymptotic forms the reader can verify that, in the large $\overline{\lambda }$ regime, (58), which in this case takes the form:

$${\displaystyle \frac{v^2}{r}}={\displaystyle \frac{GM}{r^2}}$$

(64)

is indeed satisfied for any $\overline{\lambda }$.

It is conjectured that the above result generalizes as follows to the case, of multiple, gravitationally interacting, scale-independent masses $m_{\mathrm{p}}$ (p a particle index): The smooth, $\lambda$-dependent family of solutions ${\zeta }_{\mathrm{p}}(t,\lambda )$ to (57)(48) for a bound system, approaches

$${\zeta }_{\mathrm{p}}(t,\lambda )\underset{\lambda \gg 1}{\sim }{\lambda }^{-2/3}{\zeta }_{\mathrm{p}}^{\mathrm{N}}\left(\lambda t\right)$$

(65)

at large scale, where $\left\{{\zeta }_{\mathrm{p}}^{\mathrm{N}}\right\}$ are Newtonian paths of interacting point-masses $\left\{m_{\mathrm{p}}\right\}$. The scaling form (65) is an exact symmetry of Newtonian gravity, and since well behaved solutions for the flow (57)(48) must approach at large $\lambda$ a smooth-in-scale family of Newtonian solutions, a solution other than (65) seems impossible in generic cases.

Finally, for a scale-dependent mass in (60), an $r(\lambda ;\omega )$ is obtained by the large-scale regularity condition which is not of the form $r\left(\overline{\lambda }\right)$. This results in a rotation curve $v\left(r_1\right)\equiv \omega \left(r_1\right)r_1$ which is not flat at large $r_1$, and an $r\left(\lambda \right)$ which, depending on the form of $M\left(\lambda \right)$, may not even converge to zero at large $\lambda$.

Moving to a more realistic representation of disc galaxies. For a general gravitational potential, $\mathsf{\Phi }$, sourced via (48) by a planar, axially symmetric mass density $T_{00}^{\mathrm{mat}}\left(\rho \right)$—$\rho$ being the radial distance from the galactic axis in the galactic plane—and for a test mass circularly orbiting the symmetry axis in the galactic plane, the counterpart of (60) reads

$$\lambda {\displaystyle \frac{\mathrm{d}r}{\mathrm{d}\lambda }}={\ell }_0^2\left({\omega }^2{\lambda }^{2}r-{\partial }_{\rho }\mathsf{\Phi }(r,\lambda )\right)-r$$

(66)

The time-independent mass density, $T_{00}^{\mathrm{mat}}$, is approximated by a (sufficiently dense) collection of concentric line-rings, each composed of a (sufficiently dense) collection of particles evenly spaced along the circumference. Next, recall that in the above warm-up exercise, the solution of (66) for each such particle is well behaved only for one, carefully tuned value of $\omega$. This sensitivity results form instability of the o.d.e. (66) in the $+\lambda$ direction, inherited from that of $A_{\mu }$, and is not a peculiarity of the Coulomb potential. To find the rotation curve one needs to simultaneously propagate with (66) each ring—or rather a single representative particle from each ring—$r_{r_1}\left(\lambda \right)$, using an initial guess for $\omega \left(r_1\right)$ (where $r_1$ is now a ring index, labeling the ring whose radius at $\lambda =1$ equals $r_1$, viz., $r_{r_1}\left(1\right)=r_1$). Unlike in the previous case, the (mean-field) Newtonian potential of the disc at scale $\lambda$, $\mathsf{\Phi }(\rho ,\lambda )$, solution of (42), must be re-computed at each $\lambda$. The rotation curve is obtained as that (unique) guess, $\omega \left(r_1\right)$, for which no ring diverges in the limit $\lambda \to \infty$. In so finding the rotation curve the scale dependence of individual particles comprising the disc needs to be specified. If those are fixed-point particles then their mass is scale-independent by definition. Moreover, as mentioned above, a scale-dependent mass leads to manifest contradictions with observations. In light of this, a scale independent mass is assumed modulo some caveats discussed in Section 3.4. Note the implication of the scale-invariant-mass approximation, applicable to any gravitating system: Although ${\zeta }^i(t,s=0)$ is the desired spacetime path, by construction and the s-translation invarinace of (57), ${\zeta }^i(t,s)$ for any s would also be a permissible path at $s=0$. In other words, (57) with the regularity condition at $s\to \infty$, generates a continuous family of spacetime paths which could be observed at any fixed native scale, $s=0$ in particular.

The algorithm described above for finding the rotation curve, although conceptually straightforward, could be numerically challenging and will be attempted elsewhere (see below a more promising algorithm). However, much can be inferred from it without actually running the code. Mass tracers lying at the outskirts of a disc galaxy experience almost the same, $-MG/r$ potential, where M is the galactic mass, independently of $\lambda$. This is clearly so at $\lambda =1$, as higher order multipoles of the disc are negligible far away from the galactic center, but also at larger $\lambda$, as all masses comprising the disc converge towards the center, albeit at different paces. The analytic solution (62) can therefore be used to a good approximation for such traces, implying the following power-law relation between the asymptotic velocity, $v_{\mathrm{f}}$, of a galaxy’s rotation curve and its mass, M,

$$M={\left({\displaystyle \frac{6}{\pi }}\right)}^{1/2}{\displaystyle \frac{{\ell }_0}{Gc}}{v}_{\mathrm{f}}^3$$

(67)

Such an empirical power law, relating M and $v_{\mathrm{f}}$, is known as the Baryonic Tully-Fisher Relation (BTFR), and is the subject of much controversy. There is no concensus regarding the conssistency of observations with a zero intrinsic scatter, nor is there an agreemnet about the value of the slope—3 in our case—when plotting $logM$ vs. $log{v}_{\mathrm{f}}$. Some groups [5] see a slope $\sim 3$ while other [6] insisting it is closer to 4 (both ‘high quality data’ representatives, using primary distance indicaors). While some of the discrepancy in slope estimates can be attributed to selection bias and different methods of estimating the galactic mass, the most important factor is the inclusion of relatively low-mass galaxies in the latter. When restricting the mass to lie above $\sim {10}^{10}{M}_{\odot }$, almost all studies support a slope close to 3. The recent study [7] which includes some new, super heavy galaxies, found a slope $\sim 3.26$ and a $logM/M_{\odot }$-axis intercept of $\sim 3.3$ for the massive part of the graph. Since the optimization method used in finding those two parameters is somewhat arbitrary, imposing a slope of 3 and fitting for the best intercept is not a crime against statistics. By inspection this gives an intercept of $\sim 4.2$, consistent with [5], which by (67) corresponds to ${\ell }_0\sim 4.6\times {10}^{20}$ km.

With an estimate of ${\ell }_0$ at hand, yet another prediction of our model can be put to test, pertaining to the radius at which the rotation curve transitions to its flat part. The form (62) of $r\left(\lambda \right)$ implies that the transition from the scaling to the coarsening regime occurs at $\overline{\lambda }\approx \sqrt{2/3}$. At that scale the radius assumes a value $r_{\mathrm{tr}}\approx {\left(MG{\ell }_0^2/c^2\right)}^{1/3}=v_{\mathrm{f}}{\ell }_0/c$, which is also the radius at which the force term equals the scaling term. Using standard units where velocities are given in km/s and distances in kpc, gives $r_{\mathrm{tr}}\approx {\displaystyle \frac{1}{20}}{v}_{\mathrm{f}}$. Now, in galaxies with a well-localized center—a combination of a massive bulge and (exponential) disc—most of the mass is found within a radius $r_{\mathrm{M}}<r_{\mathrm{tr}}$, lying to the right of the Newtonian curve’s maximum. Approximating the potential at $r\ge r_{\mathrm{M}}$ by $-GM/r$, the transition of the rotation curve from scaling to coarsening, with its signature rise from a flat part seen in Figure 2, is expected to show at $r_{\mathrm{tr}}$, followed by a convergence to the galaxy-specific Newtonian curve. This is corroborated in all cases—e.g. galaxies NGC2841, NGC3198, NGC2903, NGC6503, UGC02953, UGC05721, UGC08490... in Figure 12 of [8]

The above sanity checks indicate that the rotation curve predicted by our model cannot fall too far from that observed, at least for massive galaxies; it is guarantied to coincide with the Newtonian curve near the galactic center, depart from it approximately where observed, eventually flattening at the right value. However, these checks do not apply to diffuse, typically gas dominated galaxies, several orders of magnitude lighter. More urgently, a slope$=3$ is difficult to reconcile with [6] which finds a slope$\approx 4$ when such diffuse galaxies are included in the sample. Below we therefore point to two features of the proposed model possibly explaining said discrepancy. First, our model predicts that $v_{\mathrm{f}}$ attributed in [6] to such galaxies would turn out to be a gross overestimation should their rotation curves be significantly extended beyond the handful of data points of the flat portion. To see why, consider an alternative, reverse scheme for finding a rotation curve, also applicable to any bound system:

1. Start with a guess for the mass distribution of a galaxy at some large enough scale, ${\lambda }_{>}$, such that its rotation curve is fully Newtonian (If our conjecture regarding (65) is true then the flow to even larger $\lambda$ is guaranteed not to diverge for any such initial guess).

2. Let this Newtonian curve flow via (66)(48) to $\lambda =1$—no divergence problem in this, stable direction of the flow—comparing the resultant mass distribution at $\lambda =1$ with the observed distribution

3. Repeat step 1 with an improved guess based on the results of 2, until an agreement is reached.

By construction the solution curve is Newtonian at ${\lambda }_{>}$, having a $v\sim \sqrt{GM/r}$ tail past the maximum, whose rightmost part ultimately evolves into the flat segment at $\lambda =1$. We can draw two main distinctions between the flows to $\lambda =1$ of massive and diffuse galaxies’ rotation curves. First, since the hypothetical Newtonian curve at $\lambda =1$—that which is based on baryonic matter only—is rising/leveling at the point of the outmost velocity tracer in the diffuse galaxies of [6], we can be certain that this tracer was at the rising part/maximum of the ${\lambda }_{>}$ curve, rather than on its $v\sim \sqrt{GMr}$ tail as in massive galaxies. This means that, in massive galaxies, the counterpart of the short, flat segment of a diffuse galaxy’s r.c., is rather the short flat segment near its maximum, seen in most such galaxies near the maximum of the hypothetical Newtonian curve. Second, had tracers further away from the center been measured in diffuse galaxies, the true flat part would have been significantly lower relative to this maximum than in massive galaxies. With some work this can be shown via the inhomogeneous flow of $v\equiv \lambda \omega r$ derived from (66)

$$\overline{\lambda }{\partial }_{\overline{\lambda }}v={\overline{\lambda }}^{2}v-{\partial }_{\rho }\mathsf{\Phi }(r,\overline{\lambda })\overline{\lambda }$$

(68)

where r is a solution of (66) (re- expressed as a function of $\overline{\lambda }$). The gist of the argument is that, solutions of (68) deep in the coarsening regime (i.e., ${\ell }_0^2{c}^{-2}Gm\left(r\right)/r^2\gg r$, with $m\left(r\right)$ the enclosed mass at r) upon flowing to smaller $\overline{\lambda }$, decay approximately as ${\overline{\lambda }}^{1/3}$, whereas in the scaling regime they remain constant (see also (63)). In massive galaxies the entire flow from ${\overline{\lambda }}_{>}$ to $\overline{\lambda }=1$ of a tracer originally at the maximum of the r.c., $v_{\max }$, is in the coarsening regime, while in diffuse galaxies it transitions to a scaling regime.

The second possible explanation for the slopes discrepancy, which could further contribute to an intrinsic scatter around a straight BTFR, involves a hitherto ignored transparent component of the energy-momentum tensor. As emphasized throughout the paper, the A-field away from a non-uniformly moving particle (almost solving Maxwell’s equations in vacuum) necessarily involves both advanced and retarded radiation. Thus even matter at absolute zero constantly ‘radiates’, with advanced fields compensating for (retarded) radiation loss, thereby facilitating zero-point motion of matter. The A-field at spacetime point $(t,\mathit{x})$ away from neutral matter is therefore rapidly fluctuating, contributed by all matter at the intersection of its worldline with the light-cone of $(t,\mathit{x})$. We shall refer to it as the Zero Point Field (ZPF), a name borrowed from Stochastic Electrodynamics although it does not represent the very same object. Being a radiation field, the ZPF envelopes an isolated body with an electromagnetic energy ‘halo’, decaying as the inverse distance squared—which by itself is not integrable!—merging with other halos at large distance. Such ‘isothermal halos’ served as a basis for a ‘transparent matter’ model in a previous work by the author [2] but in the current context its intensity likely needs to be much smaller to fit observations. Space therefore hosts a non-uniform ZPF peaking where matter is concentrated, in a way which is sensitive to both the type of matter and its density. This sensitivity may result both in an intrinsic scatter of the BTFR, and in a systematic departure from ZPF-free slope=3 at lower mass. Indeed, in heavy galaxies, typically having a dominant massive center, the contribution of the halo to the enclosed mass at $r_{\mathrm{tr}}$ is tiny. Beyond $r_{\mathrm{tr}}$ orbiting masses transition to their scaling regime, minimally influenced by additional increase in the enclosed mass at r. The situation is radically different in light, diffuse galaxies, where the ratio of ${\rho }_{\mathrm{ZPF}}/{\rho }_{\mathrm{baryon}}$ is much higher throughout the galaxy, and much more of the non-integrable tail of the halo contributes to the enclosed mass at the point where velocity tracers transition to their scaling regime. This under estimation of the effective galactic mass, increasing with decreasing baryonic mass, would create an illusion of a BTFR slope greater than 3.

3.2.2. Other Probes of ‘Dark Matter’

Disc galaxies are a fortunate case in which the worldline of a body transitions from scaling to coarsening at a common scale along its entire worldline (albeit different scales for different bodies). They are also the only systems in which the velocity vector can be inferred solely from its projection on the line-of-sight. In pressure supported systems, e.g., globular clusters, elliptical galaxies or galaxy clusters, neither is true. Some segments of a worldline could be deep in their scaling regime while others in the coarsening, rendering the analysis of their collective scale flow more difficult. The above reverse scheme leverages the fact that, all the worldlines of a bound system are deep in their coarsening regime at sufficiently large scale, where their fixed-$\lambda$ dynamics is well approximated using Newtonian gravity. Starting with such a Newtonian system at sufficiently large ${\lambda }_{>}$, the integration of (57) to small $\lambda$ is in its stable direction, hence not at risk of exploding for any initial choice of Newtonian paths. If the Newtonian system at ${\lambda }_{>}$ is chosen to be virialized, a ‘catalog’ of solutions of pressure supported systems extending to arbitrarily small $\lambda$ can be generated, and compared with line-of-sight velocity projections of actual systems. As remarked above, the transition from coarsening to scaling generally doesn’t take place at a common scale along the worldline of any single member of the system. However, if we assume that there exists a rough transition scale, ${\lambda }_{\mathrm{tr}}$, for the system as a whole in the statistical sense, which is most reasonable in the case of galaxy clusters, then immediate progress can be made. Since in the scaling regime velocities are unaltered, the observed distribution of the line-of-sight velocity projections should remain approximately constant for $\lambda <{\lambda }_{\mathrm{tr}}$, that of a virialized system, viz., Gaussian of dispersion ${\sigma }_v$. On the other hand, at $\lambda >{\lambda }_{\mathrm{tr}}$ a virialized system of total mass M satisfies

$${\sigma }_v\approx \sqrt{GM/r}$$

(69)

where ${\sigma }_v$ is the velocity dispersion, and r is the radius of the system, which is just (64) with ${\sigma }_v\mapsto v$. On dimensional grounds it then follows that ${\sigma }_v$ would be the counterpart of $v_f$ from (67), implying ${\sigma }_v\propto M^{1/3}$ which is in rough agreement with observations. The proportionality constant can’t be exactly pinned using such huristic arguments, but its observed value is on the same order of magnitude as that implied by (67).

Applying our model to gravitational lensing in the study of dark matter requires better understanding of the nature of radiation. This is murky territory even in conventional physics and in next section initial insight is discussed. To be sure, Maxwell’s equations in vacuum are satisfied away from $J_{\mu }$, although only ‘almost so’, as discussed in Section 3. However, treating them as an initial value problem, following a wave-front from emitter to absorber is meaningless for two reasons. First, tiny, $O\left({\ell }_0^{-1}\right)$ local deviations from Maxwell’s equations could become significant when accumulated over distances on the order of ${\ell }_0$. Second, in the proposed model extended particles ‘bump into one another’ and their centers jolt as a result—some are said to emit radiation and other absorb it, and an initial-value-problem formulation is, in general, ill-suited for describing such process. Nonetheless, incoming light—call it a photon or a light-ray—does posses an empirical direction when detected. In flat spacetime this could only be the spatial component of the null vector connecting emission and absorption events, as it is the only non arbitrary direction. A simple generalization to curved spacetime, involving multiple, freely falling observers, selects a path, ${\zeta }^{\mu }$, everywhere satisfying the light-cone condition${\dot{\zeta }}^2\left(\tau \right)\equiv 0$. Every null geodesic satisfies the light-cone condition, but not the converse. In ordinary GR, the only non arbitrary path connecting emission and absorption events which respects the light-cone condition and locally depends on the metric and its first two derivatives is indeed a null geodesic. In our model, a solution of (59) which is well behaved on all scales, further satisfying the light-cone condition for $\lambda \to \infty$, is an appealing candidate: It must (almost) be a null geodesic in strong fields/at large scales, and when transitioning to weak fields/small scales, the light-cone condition inherited from large scales is preserved by the scaling operation; Indeed, denoting ${\dot{\zeta }}^{\mu }\equiv {\partial }_{\tau }{\zeta }^{\mu }$, applying ${\partial }_{\tau }$ to both sides of

$${\partial }_s{\zeta }^{\mu }=\tau {\dot{\zeta }}^{\mu }-z^{\mu }(\zeta ,s)$$

multiplying the result by $g_{\mu \nu }{\dot{\zeta }}^{\nu }$, one gets

$${\partial }_s{\dot{\zeta }}^2=\tau {\partial }_{\tau }\left({\dot{\zeta }}^2\right)+2{\dot{\zeta }}^2$$

(70)

In arriving at (70) use of (24) has been made in converting ${\partial }_s{g}_{\mu \nu }$ to spacetime derivatives of $g_{\mu \nu }$ and $z^{\mu }$, and a term ${\ell }_g^2{G}_{\mu \nu }{\dot{\zeta }}^{\mu }{\dot{\zeta }}^{\nu }$ has been neglected (which is justified in the context of cosmology, discussed in Sec. 3.4.2). Likewise, all spacetime derivatives of $g_{\mu \nu }$ have been neglected (weak field domain). It then follows that ${\partial }_s{\dot{\zeta }}^2(\tau ,s)\equiv 0$ if ${\dot{\zeta }}^2(\tau ,s)\equiv 0$.

Furthermore, in GR the deflection angle of a light ray due to gravitational lensing, by a compact gravitating system of mass M, is $\varphi =4GM/\left(c^{2}R\right)$, where R is the impact parameter. When ${\zeta }^{\mu }$ is in its scaling regime, our model’s $\varphi$ remains constant. If the system is likewise in its scaling regime, (69) implies that its virial mass, $M_{\mathrm{vir}}={\sigma }_v^{2}r/G$, scales according to $M_{\mathrm{vir}}\mapsto {\lambda }^{-1}{M}_{\mathrm{vir}}$, as does the impact parameter of ${\zeta }^{\mu }$, $R\mapsto {\lambda }^{-1}R$. The conventional mass estimate based on the virial theorem, of this $\lambda$-dependent family of gravitating systems, would then agree with that which is based on (conventional) gravitational lensing, $M_{\mathrm{lens}}=c^{2}R\varphi /4G$, up to a constant, common to all members—recall that this entire family appears in the ‘catalog’ of $\lambda =1$ systems. Extending this family to large $\lambda$, the two estimates will coincide by construction. Thus if the system and ${\zeta }^i$ transition approximately at a common $\lambda$, this proportionality constant can only be close to 1. This is apparently the case in most observations pertaining to galaxy clusters. Nonetheless, the two methods according to our model need not, in general, produce identical results. The degree to which they disagree depends on the exact scale-flow of ${\zeta }^{\mu }$, which is one of those calculations avoided thus far, involving a path non-uniformly (in scale) transitioning from coarsening to scaling.

A final caveat regarding the application of (59) is that, we have no reason to expect that the light-cone condition is satisfied exactly (other than in the $\lambda \to \infty$ limit). Nonetheless, in and of itself this caveat does not invalidate such solutions, even as candidates for exact solutions, so long as no conflict with observations arises.

3.3. Quantum Mechanics as a Statistical Description of the Realistic Model

The basic tenets of classical electrodynamics (19), (30) and (31), which must be satisfied at any scale on consistency grounds (up to neglected curvature terms), strongly constrain also statistical properties of ensembles of members in $\mathcal{S}$, and in particular constant-$\lambda$ sections thereof. In a previous paper by the author [1] it was shown that these constraints could give rise to the familiar wave equations of QM, in which the wave function has no ontological significance, merely encoding certain statistical attributes of the ensemble via the various currents which can be constructed from it. It is through this statistical description that ℏ presumably enters physics, and so does ‘spin’ (see below).

This somewhat non-committal language used to describe the relation between QM wave-equations and the basic tenets is for a reason. Most attempts to provide a realist (hidden variables) explanation of QM follow the path of statistical mechanics, starting with a single-system theory, then postulating a ‘reasonable’ ensemble of single-systems—a reasonable measure on the space of single-system solutions—which reproduces QM statistics. Ignoring the fact that no such endeavor has ever come close to fruition, it is rarely the case that the measure is ‘natural’ in any objective way, effectively defining the statistical theory/measure (uniformity over the impact parameter in an ensemble representing a scattering experiment being an example of an objectively natural attribute of an ensemble). Even the ergodicity postulate, as its name suggests, is a postulate—external input. When sections of members in $\mathcal{S}$ are the single-systems, the very task of defining a measure on such a space, let alone a natural one, becomes hopeless. The alternative approach adopted in [1] is to derive constraints on any statistical theory of single-systems respecting the basic tenets, showing that QM non-trivially satisfies them. QM then, like any measure on the space of single-system solutions, is postulated rather than derived, and as such enjoys a fundamental status, on equal footing with the single-system theory. Nonetheless, the fact that the QM analysis of a system does not require knowledge of the system’s orbit makes it suspicious from our perspective. And since a quantitative QM description of any system but the simplest ones involves no less sorcery than math, that fundamental status is still pending confirmation (refutation?).

Of course, the basic tenets of classical electrodynamics are respected by all (sections of-) members of $\mathcal{S}$, not only those associated Dirac’s and Schrödinger’s equations. The focus in [1] on ‘low energy phenomena’ is only due to the fact that certain simplifying assumptions involving the self-force can be justified in this case. In fact, the current realization of the basic tenets, involving fields only instead of interacting particles, is much closer in nature to the QFT statistical approach than to Schrödinger’s.

3.3.1. The Origin of Quantum Nonlocality

“Multiscale locality", built into the proposed formalism, readily dispels one of QM’s greatest mysteries—its apparent non-local nature. In a nutshell: Any two particles, however far apart at our native scale, are literally in contact at sufficiently large scale.

Two classic examples where this simple observation invalidates conventional objections to local-realist interpretations of QM are the following. The first is a particle’s ability to ‘remotely sense’ the status of the slit through which it does not pass, or the status of the arm of an interferometer not traversed by it (which could be a meter away). To explain both, one only needs to realize that for a giant physicist, a fixed-point particle is scattered from a target not any larger than the particle itself, to which he would attribute some prosaic form-factor; At large enough $\lambda$ the particle literally passes though both arms of the interferometer (and through none!). This global knowledge is necessarily manifested in the paths chosen by it at small $\lambda$. Of course, at even larger $\lambda$ the particle might also pass through two remote towns etc., so one must assume that the cumulative statistical signature of those infinitely larger scales is negligible. A crucial point to note, though, is that the basic tenets, which imply local energy-momentum conservation at laboratory scales, are satisfied at each $\lambda$ separately. For this large-$\lambda$ effect to manifest at $\lambda =1$, local energy-momentum conservation alone must not be enough to determine the particle’s path, which is always the case in experiment manifesting this type of nonlocality. Inside the crystal serving as mirror/beam-splitter in, e.g., a neutron interferometer, the neutron’s classical path (=paths of bulk-motion derived from energy-momentum conservation) is chaotic. Recalling that, what is referred to as a neutron—its electric neutrality notwithstanding—only marks the center of an extended particle, and that the very decomposition of the A-field into particles is an approximation, even the most feeble influence of the A-field awakened by the neutron’s scattering, traveling through the other arm of the interferometer, could get amplified to a macroscopic effect. This also provides an alternative, fixed-scale explanation for said ‘remote sensing’. In the double-slit experiment such amplification is facilitated by the huge distance of the screen from the slits compared with their mutual distance.

The second kind of nonlocality is demonstrated in Bell’s inequality violations. As with the first kind, the conflict with one’s classical intuition can be explained both at a fixed scale, or as a scale-flow effect. Starting with the former, and ever so slightly dumbing down his argument, Bell assumes that physical systems are small machines, with a definite state at any given time, propagating (deterministically or stochastically) according to definite rules. This generalizes classical mechanics, where the state is identified with a point in phase-space and the evolution rule with the Hamiltonian flow. However, even the worldlines of particles in our model, represented by sections of members in $\mathcal{S}$, are not solutions of any (local) differential equation in time. Considering also the finite width of those worldlines, whose space-like slices Bell would regard as possibly encoding their ‘internal state’, it is clear that his modeling of a system is incompatible with our model; particles are not machines, let alone particle physicists. Spacetime ‘trees’ involved in Bell’s experiments—a trunk representing the two interacting particles, branching into two, single particle worldlines—must therefore be viewed as a single whole, with Bell’s inequality being inapplicable to the statistics derived from ‘forests’ of such trees.5 This spacetime-tree view gives rise to a scale-flow argument explaining Bell’s inequality violations: The two branches of the tree shrink in length when moving to larger scale, eventually merging with the trunk and with one another. Thus the two detectors at the endpoints of the branches cannot be assumed to operate independently, as postulated by Bell.

3.3.2. Fractional Spin

Fractional spin is regarded as one of the hallmarks of quantum physics, having no classical analog, but according to [1], much like ℏ, it is yet another parameter—discrete rather than continuous—entering the statistical description of an ensemble. At the end of the day, the output of this statistical description is a mundane statement in ${\mathbb{R}}^3$, e.g., the scattering cross-section in a Stern-Gerlach experiment, which can be rotated with $O\left(3\right)$. Neither Bell’s- nor the Kochen-Specker theorems are therefore relevant in our case as the spin is not an attribute of a particle. For this reason the spin-0 particle from Section 3.1.2 is a legitimate candidate for a fractional-spin particles, such as the proton, for its ‘spin measurement/polarization’ along some axis is by definition a dynamical happening, in which its extended world-current bends and twists, expands and contracts in a way compatible with- but not dictated by the basic tenets. As stressed above, there is no natural measure on the space of such objects, and the appearance of two strips on Stern$\&$Gerlach’s plate rather than one, or three etc. need not have raised their eyebrows. Nonetheless, the proposed model does support spinning solutions, viz. $\mathit{J}\not\equiv 0$ in the rest frame of the particle, and there is a case to be made that those are more likely candidates for particles normally attributed with a spin, integer or fractional.

3.3.3. Photons and Neutrinos (Or Illusion Thereof?)

Einstein invented the ‘photon’ in order to explain the apparent violation of energy conservation occurring when an electron is jolted at a constant energy from an illuminated plate even when the plate is placed far enough from the source, such that the time-integrated Poynting flux across it becomes smaller than the energy of the jolted electron. It is entirely possible that Einstein’s explanation can be realized in the proposed formalism, although the rest-frame analysis of a fixed-point particle from Section 3.1.2 must obviously be modified for massless (neutral) particles which might further require extending $A_{\mu }$ and $Z^{\mu }$ to include distributions. Maxwell’s equations would then act as the photonic counterpart of a massive-particle’s QM wave equation, describing the statistical aspects of ensembles of photons. Indeed, since in a ‘lab’ of dimension $\ll {\ell }_0$ individual photons (almost) satisfy the basic tenets of classical electrodynamics (and (32)) for a chargeless current (i.e., $\int {\mathrm{d}}^3{J}_0=0$), the construction from [1] would result in Maxwell’s equations, with the associated ${\mathsf{\Theta }}_{\mu \nu }$ being the ensemble energy-momentum tensor. However, since the A-field (almost) satisfies Maxwell’s equations regardless of it being a building block of photons, it is highly unlikely that photons exhaust all radiation-related phenomena. For example, is there any reason to think that a radio antenna transmits its signal via radio photons, rather than radio (A-) waves? This suggests an alternative explanation for photon-related phenomena, which does not require actual, massless particles. Its gist is that, underlying the seeming puzzle motivating Einstein’s invention of the photon, is the assumption that an electron’s radiation field is entirely retarded which, as emphasized throughout the paper, cannot be the case for the A-field. Advanced radiation converging on the electron could supply the energy necessary to jolt it, further facilitating violation of Bell’s inequality in entangled ‘photons’ experiments. This proposal, first appearing in [1] and further developed in [2], was, at the time, the only conceivable realist explanation of photon related phenomena. In the proposed model, apparently capable of representing ‘light corpuscles’, it may very well be the wrong explanation. Photons would then be just ephemeral massless particles created in certain structural transitions of matter, then disappearing when detected. Note that these two processes are entirely mundane, merely representing a relatively rapid changes in $A_{\mu }$ and $Z^{\mu }$ at the endpoints of a photon’s (extended) worldline. Such unavoidable transient regions might result in an ever-so-slight smoothing of said distributions, which are otherwise excluded from $\mathcal{S}$.

“God is subtle but not malicious" was Einstein’s response to claims that further repetitions of the Michelson-Morley experiment did show a tiny directional dependence of the speed of light. This attitude is adopted vis-a-vis the neutrino’s mass problem. All direct measurements based on time-of-flight are consistent with the neutrino being massless; the case for a massive neutrino relies entirely on indirect measurements and a speculative extension of the Standard Model. Neutrinos would then be quite similar to photons, probably spinning ($\mathit{J}\ne 0$), whose creation and annihilation involve structural transitions at the subatomic scale. However, as with photons, and even more so due to their elusiveness, neutrinos might not be the full story, or even the real one. The classical model of photons cited above assumes that only $A_{\mu }$ contributes to the radiative $T_{\mu \nu }$ which is therefore identified with ${\mathsf{\Theta }}_{\mu \nu }$. In the proposed model $T_{\mu \nu }$ consists also of

$$P_{\mu \nu }=-{\displaystyle \frac{1}{4\pi G{\ell }_g^2}}\left({\partial }_{\mu }{Z}_{\nu }+{\partial }_{\nu }{Z}_{\mu }\right)$$

with $Z_{\nu }$ satisfying (32) in the flat spacetime approximation, rewritten here

$${\square }^2{Z}_{\nu }+{\partial }_{\nu }{\partial }^{\mu }{Z}_{\mu }=-4\pi G{\ell }_g^2{F}_{\nu \mu }{J}^{\mu }$$

This is a massless wave equation, not too dissimilar to Maxwell’s, therefore expected to participate in radiative, energy-momentum transfer. However, two features set it apart. First, the two terms on the l.h.s. of (32) enter with the ‘wrong’ relative sign, spoiling gauge covariance. As a result an extra longitudinal mode exists, i.e., $Z^{\mu }={ϵ}^{\mu }f\left(k_{\nu }{x}^{\nu }\right)$ with ${ϵ}^{\mu }\propto k^{\mu }$ (which in the Maxwell case is a pure gauge), on top of the two transverse modes, ${ϵ}^{\mu }{k}_{\mu }=0$. Second, unlike ${\mathsf{\Theta }}_{\mu \nu }$, $P_{\mu \nu }$ is only linear in $Z^{\mu }$, an impossibility for a Noether current. Combined, these two features imply that only the longitudinal mode can radiate energy-momentum and only during transient, ‘structural changes’ to the radiating system. Indeed, consider the integral of the energy flux of $P_{\mu \nu }$ over $T=S^2\times \mathbb{R}$ in Figure 1.

$$\int \mathrm{d}t{\int }_{S^2}{P}_{0i}\mathrm{d}{\sigma }^i\equiv \int \mathrm{d}t{\int }_{S^2}\mathit{P}·\mathrm{d}\sigma$$

(71)

where $S^2$ is a large sphere centered at the location of the system and $\mathrm{d}\sigma$ is an outward pointing vector orthogonal to $S^2$ of length $\mathrm{d}\sigma$. Clearly, only the longitudinal mode, whose energy flux at each point on $S^2$ is $\propto k_0\mathit{k}\Vert \mathrm{d}\sigma$, contributes to the integral (71). Moreover, we saw in Section 3.1.1 that outside of fixed-points, $z^{\mu }$ must be negligible. So long as the system qualifies as a fixed-point, as during bulk acceleration, no net flux is being generated by it, and it is therefore only while transitioning between distinct fixed-points that $\mathit{P}$ is involved in energy changes (and even then only its ${\partial }_r{Z}_0$ piece—r being the radial coordinate when $S^2$ and the system are co-centered at the origin—as the ${\partial }_t{Z}_r$ piece integrates to zero over time). The Z-field is therefore a natural candidate for a ‘classical neutrino field’, whose relation to neutrino phenomena parallels that of the A-field to photon phenomena. As with photons, it is a particle’s advanced Z-field converging on it which supplies the energy-momentum necessary to jolt it, conventionally interpreted as the result of being struck by a neutrino. Similarly, hitherto ignored retarded Z-field is allegedly generated in structural changes of a system, e.g., when nuclei undergo $\beta$-decay. As pointed out in Section 3.3.2 above, the (fractional) spin-$\frac{1}{2}$ attributed to the neutrino, as is the spin-1 of the photon, only labels the statistical description of phenomena involving such jolting of charged particles.

3.4. Cosmology

Cosmological models are stories physicists entertain themselves with; they can’t truly know what happened billions of years ago, billions of light-years away, based on the meager data collected by telescopes (which covers $0\%$ of the electromagnetic spectrum). Moreover, in the context of the proposed model, the very ambition implied by the term “cosmology" is at odds with the humility demanded of a physicist, whose entire observable universe could be another physicist’s oven. On the other hand, astronomical observations associated with cosmology, also serve as a laboratory for testing ‘terrestrial’ physical theories, e.g., atomic-, nuclear-, quantum-physics, and this would be particularly true in our case, where the large and the small are so intimately interdependent. When the most compelling cosmological story we can devise, fitting those observations, requires contrived adjustments to terrestrial-physics theories, confidence in those, including GR, should be shaken.

Reluctantly, then, a cosmological model is outlined below. Its purpose at this stage is not to challenge $\mathsf{\Lambda }$CDM in the usual arena of precision measurements, but to demonstrate how the novel ingredients of the proposed formalism could, perhaps, lead to a full-fledged cosmological model free of the aforementioned flaw.

3.4.1. A Newtonian Cosmological Model

As a warm-up exercise, we wish to solve the system (57)(48) for a spherical, uniform, expanding cloud of massive particles originating from the scaling center (without loss of generality). The path of a typical particle is described by

$${\zeta }^i(t,\lambda )=r^{i}a(t,\lambda )$$

(72)

where $r^i$ a constant vector. It is easily verified that the same homogeneous expanding cloud would appear to an observer fixed to any particle, not just the one at the origin. The mass density of the cloud depends on a via $\rho \propto a^{-3}$, retaining its uniformity at any time and scale if creation/annihilation of matter in scale is uniform across space. The gravitational force acting on a particle is given by $f^i=-{\displaystyle \frac{4\pi G}{3}}\rho (t,\lambda )r^{i}a$ (the uniform vacuum energy is ignored as its contribution to the force can only vanish by symmetry) and (57) gives a single, particle-independent equation for a

$$\lambda {\partial }_{\lambda }a={\ell }_0^2\left(-\ddot{a}-{\displaystyle \frac{4\pi G}{3}}\rho a\right)-a+t\dot{a}$$

(73)

with $\dot{a}\equiv {\partial }_{t}a$ etc.

Two types of solutions for (73) which are well behaved at all scales should be distinguished: Bounded and unbounded. In the former $a(t,\lambda )$ is identically zero at $t=0$ and a $\lambda$-dependent ‘big-crunch’ time, $t_{\mathrm{f}}$. By our previous remarks, at large scale the coarsening terms—those multiplied by ${\ell }_0^2$ on the r.h.s. of (73)—dominate the flow and must almost cancel each other or else a would rapidly blow up with increasing $\lambda$. The resulting necessary condition for a regular $a(t,\lambda )$ on all scales is a $\lambda$-dependent o.d.e. in time, which is simply the time derivative of the (first) Friedmann equation for non-relativistic matter

$${\dot{a}}^2+k\left(\lambda \right)={\displaystyle \frac{8\pi G}{3}}\rho (t,\lambda )a^2$$

(74)

The k above disappears as a result of this derivative, meaning that it resurfaces as a second integration constant of any magnitude—not just $k\in \left\{-1,0,1\right\}$. Denoting

$$\underset{t\to 0}{lim}\rho (t,\lambda )a^3(t,\lambda ):={\rho }_0\left(\lambda \right)$$

bounded solutions in which mass is conserved in time are therefore described by some flow in $k-{\rho }_0$ parameter space for which ${max}_{t}a(t,\lambda )$ shrinks to zero for $\lambda \to \infty$. For example, as k in (74) plays the role of minus twice the total energy of the explosion per unit mass, for a scale independent ${\rho }_0$, $k\left(\lambda \right)$ monotonically increases with increasing $\lambda$.

Given a solution of (74) at large enough $\lambda$ one can then integrate (73) in its stable, small $\lambda$ direction, where the scaling piece becomes important, but due to the $a(t=0,\lambda )=a(t_{\mathrm{f}},\lambda )\equiv 0$ constraint, some parts of a solution remain deep in their coarsening regime. The same is true for unbounded solutions, but in this case there is no $a(t,{\lambda }_{\mathrm{large}})$ to start from, rendering the task of finding solutions more difficult; instead of b.c. $a(t,\lambda )=0$ for $t=0,t_{\mathrm{f}}$, we have $a(t_{\mathrm{i}},\lambda )=0$ for some initial time, $t_{\mathrm{i}}$, and the large-t asymptotic $\dot{a}(t,\lambda )\sim v_{\infty }$ for some $v_{\infty }\ge 0$. Note the consistency with $a=v_{\infty }t+{\lambda }^{-1}C$ for some constant C, which is an exact solution for the $\rho$-free (73) (and its only solution not wildly diverging in magnitude at large t). One exception to the hardness of the open-solution (applicable also to closed solutions) is a scaling solution, $a(t,\lambda )={\lambda }^{-1}{a}_{\sigma }\left(\lambda t\right)$, where $a_{\sigma }$ is an exact solution of (74) with ${\rho }_0\propto {\lambda }^{-1}$ and $k\left(\lambda \right)=\mathrm{const}$. Note that the asymptotic b.c. is automatically satisfied for $k\le 0$. Another is the scale invariant solution of (73), integrated backwards from $t=\infty$ to $a=0$, implicitly defining $t_{\mathrm{i}}$ (integrating forward from $a(0,\lambda )=0$ leads to nonphysical solutions).

The Newtonian-cloud model, while mostly pedagogical, nonetheless captures a way—perhaps the only way—cosmology is to be viewed within the proposed framework: It does not pertain to the Universe but rather to a universe—an expanding cloud as perceived by a dwarf amidst it. A relative giant, slicing the cloud’s orbit at a much larger $\lambda$, might classify the corresponding section as, e.g., the expanding phase of a Cepheid/red-giant, or a runaway supernova. An even mightier giant may see a decaying radioactive atom. Of course, matter must disappear in such flow to larger and larger scales—a phenomenon already encountered in the linear case which is further discussed below. The rate (in scale) at which this takes place, ${\partial }_s{\rho }_0\left(s\right)$ in the above models, must be compatible with our analysis of galaxies, where mass was assumed conserved in scale. This would be true for small enough global rate, or if around our native scale, mass annihilation takes place primarily outside galaxies (commencing in a galaxy only after scale flow has compressed it to an object currently not identified as a galaxy).

Suppose now for concreteness that a giant’s section is an expanding star. The dwarf’s entire observable universe would in this case correspond to a small sphere, non-concentrically cut from the star. The hot thermal radiation inside that sphere at ${\lambda }_{\mathrm{giant}}$, after flowing with (16) to ${\lambda }_{\mathrm{dwarf}}$, would be much cooler, much less intense, and much more uniform, except for a small dipole term pointing towards the star’s center, approximately proportional to the star’s temperature gradient at the sphere, multiplied by the sphere’s diameter. Similarly for the matter distribution at ${\lambda }_{\mathrm{dwarf}}$, only in this case the distribution of accumulated matter created during the flow is expected to decrease in uniformity if new matter is created close to existing matter. Thus the distribution of matter at ${\lambda }_{\mathrm{dwarf}}$ is proportional to the density at ${\lambda }_{\mathrm{giant}}$ only when smoothed over a large enough ball, whose radius coresponds to a distance at ${\lambda }_{\mathrm{giant}}$ much larger than the scale of density fluctuations. This would elegantly explain the so-called dipole problem [9]—the near perfect alignment of the CMB dipole with the dipole deduced from matter distribution, but with over $5\sigma$ discrepancy in magnitude; Indeed, the density and temperature inside a star typically have co-linear, inward-pointing gradients, but which differ in magnitude. Note that a uniform cloud ansatz is inconsistent with the existence of such a dipole discrepancy and should therefore be taken as a convenient approximation only, rendering the entire program of precision cosmology futile. The horizon problem of pre inflation cosmology is also trivially explained away by such orbit view of the CMB. Similarly, the tiny but well-resolved deviations from an isotropic CMB (after correcting for the dipole term) might be due to acoustic waves inside the star.

Returning to the scale-flow of $\phi \equiv \left(A_{\mu },g_{\mu \nu }\right)$ interpolating between ‘a universe’ and a star, and recalling that $\phi (x,s)$ stands for a spacetime phenomenon as represented by a physicist of native scale s, a natural question to ask is: What would this physicist’s lab notes be? A primary anchor facilitating this sort of note-sharing among physicists of different scales is a fixed-point particle, setting both length and mass standard gauges. We can only speculate at this stage what those are, but the fact that the mass of macroscopic matter must be approximately scale invariant—or else rotation curves would not flatten asymptotically—makes atomic nuclei, where most of the mass is concentrated, primary candidates. Note that in the proposed formalism the elementarity of a particle is an ill-defined concept, and the entire program of reductionism must be abandoned. For if zooming into a particle were to ‘reveal its structure’, even a fixed-point would comprise infinitely many copies of itself as part of its attraction basin.

If nuclei approximately retain their size under scale-flow to large $\lambda$, while macroscopic molecular matter shrinks, then some aspects of spacetime physics (at a fixed-scale section) must change. Instinctively, one would attribute the change to a RG flow in parameter space of spacetime theories, e.g., the Yukawa couplings of the Standard Model of particle physics, primarily that of the electron. However, this explanation runs counter to the view advocated in this paper, that (spacetime) sections should always be viewed in the context of their (scale) orbit; If the proposed model is valid, then the whole of spacetime physics is, at best, a useful approximation with a limited scope. Moreover, an RG flow in parameter space cannot fully capture the complexity involved in such a flow, where, e.g., matter could annihilate in scale (subject to charge conservation); ‘electrons’ inside matter, which in our model simply designate the A-field in between nuclei—the same A-field peaking at the location of nuclei—‘merging’ with those nuclei (electron capture?); atomic lattices, whose size is governed by the electronic Bohr radius $\sim m_{\mathrm{e}}^{-1}$, might initially scale, but ultimately change structure. At sufficiently large $\lambda$ an entire star or even a galaxy would condense into a fixed-point—perhaps a mundane proton, or some more exotic black-hole-like fixed-point which cannot involve a singularity by definition. Finally, we note that, by definition, the self-representation of that scaled physicist slicing $\phi (x,s)$ at his native scale s, is isomorphic to ours, viz., he reports being made of the same organic molecules as we are made of, which are generically different from those he observes, e.g., in the intergalactic medium. So either actual physicists (as opposed to hypothetical ones, serving as instruments to explain the mathematical flow of $\phi$) do not exist in a continuum of native scales, only at those (infinitely many) scales at which hydrogen atoms come in one and the same size; or else they do, in which case we, human astronomers, should start looking around us for odd-looking spectra, which could easily be mistaken for Doppler/gravitational shifts.

3.4.2. Relativistic Cosmology

In order to generalize the Newtonian-cloud universe to relativistic velocities, while retaining the properties of no privileged location and statistical homogeneity, it is convenient to transfer the expansion from the paths of the particles to a maximally symmetric metric—a procedure facilitated by the general covariance of the proposed formalism. Formally, this corresponds to an ‘infinite cloud’ which is a good approximation whenever the size of the cloud and the distance of the observer from its edge are both much greater than ${\ell }_0$ and ${\ell }_g$. Alternatively, the cosmological principle could be postulated as an axiom. For clarity, the spatially flat ($k=0$), maximally symmetric space, with metric

$$g_{tt}=-1,g_{ij}={\eta }_{ij}{a}^2(t,\lambda ),g_{it}=0$$

(75)

is considered first, for which the only non-vanishing Christoffel symbols are

$${\mathsf{\Gamma }}_{ij}^t=a{\partial }_{t}a{\delta }_{ij},{\mathsf{\Gamma }}_{tj}^i\equiv {\mathsf{\Gamma }}_{jt}^i={\displaystyle \frac{{\partial }_{t}a}{a}}{\delta }_{\phantom{1}j}^i$$

(76)

The gravitational part, $z^{\mu }$, of the scaling field, appropriate for the description of a universe which is electrically neutral on large enough scales, i.e., $〈J_{\mu }〉\approx 0$, is given by solutions of (33) which, for the metric (75), reads

$${\nabla }^{\mu }{\nabla }_{\mu }{z}_{\nu }+{\nabla }^{\mu }{\nabla }_{\nu }{z}_{\mu }=\left\{\begin{array}{cc}-6\left(\lambda {\partial }_{\lambda }a/a\right)\left({\partial }_{t}a/a\right)\hfill & \mathrm{for}\nu =0\hfill \\ 0\hfill & \mathrm{for}\nu =1,2,3\hfill \end{array}\right.$$

(77)

However, the generally covariant boundary condition (34) “far away from matter" is not applicable here. Instead, ${\nabla }_i{z}_j+{\nabla }_j{z}_i$ is required to be compatible with the (maximal) symmetry of space—its Lie derivative along any Killing field of space must vanish.

The general form of $z^{\mu }$ consistent with the metric (75) is

$$z^t=\varphi (t,\lambda ),z^i=0\Rightarrow z_t=-\varphi ,z_i=0$$

(78)

Spatial scaling is taken care of by the metric, hence the vanishing $z^i$. This implies that, in cosmic coordinates the size of a gravitationally bound system whose outmost matter is deep in its scaling regime, e.g. a galaxy with a flat r.c., also scales as a, rather than ${\lambda }^{-1}$ in Minkowskian coordinates.

Inserting (78) into (77) results is a single equation

$${\partial }_{tt}\varphi +3{\displaystyle \frac{\dot{a}}{a}}{\partial }_t\varphi -3{\left({\displaystyle \frac{\dot{a}}{a}}\right)}^2\varphi =-3{\displaystyle \frac{\lambda {\partial }_{\lambda }a}{a}}{\displaystyle \frac{\dot{a}}{a}}$$

(79)

Importantly, $\varphi =\pm t$ is a solution when a is of either scaling forms, $a\sim {\lambda }^{\mp }{a}_{\sigma }\left({\lambda }^{\pm }t\right)$ resp. This exposes the fact that, in the generally covariant setting, the scale-direction of giants could be either $\lambda \uparrow$ (equiv. $s\uparrow$) or $\lambda \downarrow$, depending on the coordinate system and its associated solution for the scaling field. As soon becomes apparent, compatibility with the Newtonian model selects a negative $\varphi$—which at any fixed $\lambda$ contains two, free, $\lambda$ dependent integration constants, referred to below—and a $\lambda \downarrow$ direction of giants.

An important issue which must be addressed before proceeding, concerns the ontological status of the energy-momentum tensor. In GR, sourcing the Einstein tensor is a phenomenological device, equally valid when applied to the hot plasma inside a star, or to the ‘cosmic fluid’. In contrast, ${\mathsf{\Theta }}_{\mu \nu }$ and the scaling field from which $P_{\mu \nu }$ is derived, both enter (24) as fundamental quantities, on equal footing with $g_{\mu \nu }$. To make progress, this fundamental status must be relaxed, and the following way seems reasonable: The fundamental scaling field is written $Z^{\mu }=z^{\mu }+z^{\mu }$, with $z^{\mu }$ the above, coarse grained gravitational part, and $z^{\mu }$ the field inside matter. The space averages of the fundamental ${\mathsf{\Theta }}_{\mu \nu }$ and $p_{\mu \nu }$ (derived from $z^{\mu }$) are written at $〈{\mathsf{\Theta }}_{00}+p_{00}〉=\rho$, $〈{\mathsf{\Theta }}_{ij}+p_{ij}〉=a^{2}p{\delta }_{ij}$. That such coarse grained pseudo tensor, respecting the symmetries of the coarse grained metric (75), has the perfect fluid form, can easily be shown.

Plugging $Z^{\mu }$ thus defined and (75) into the metric flow (24), results in space-space and time-time components given, respectively, by

$$2a\lambda {\partial }_{\lambda }a={\ell }_g^2\left(2a\ddot{a}+{\dot{a}}^2-4a^2\pi Gp\right)+2a\dot{a}\varphi$$

(80)

$$0={\ell }_g^2\left(-3{\displaystyle \frac{{\dot{a}}^2}{a^2}}-4\pi G\rho \right)-2{\partial }_t\varphi$$

(81)

with $\rho$ and p incorporating ${\mathsf{\Theta }}_{\mu \nu }$ and $p_{\mu \nu }$ while the remaining terms are entirely due to $z^{\mu }$. Another equation which can be extracted from those two, or directly from (28), in conjunction with (79), is energy-momentum conservation in time

$${\partial }_{t}p{a}^3-{\partial }_t\left(a^3(\rho +p)\right)=0$$

(82)

Only two of the above three equations are independent due to the Bianchi identity and (79).

Equations (80) and (81) can be combined to

$$\lambda {\partial }_{\lambda }a={\ell }_g^2\left(\ddot{a}-{\displaystyle \frac{2\pi Ga}{3}}\left(\rho +3p\right)\right)-{\displaystyle \frac{1}{3}}a{\partial }_t\varphi +\dot{a}\varphi$$

(83)

Remembering that paths of co-moving masses can be deduced by analytically continuing solutions of (24), $G\mapsto -2G$, and solving (59) in the resultant metric (which for the metric (75) gives: ${\zeta }^i=r^i$, ${\zeta }^t=\tau$ with $r^i$ a constant), we might as well solve (83) directly for $G\mapsto -2G$. In accordance with $\varphi <0$ one should also change $a(t,s)\mapsto a(t,-s)$ (or $\lambda \mapsto {\lambda }^{-1}$), for $s\uparrow$ in (83) to be the direction of giants. The result is an equation which, for $p=0$, is quite similar to (73)

$$\lambda {\partial }_{\lambda }a={\ell }_g^2\left(-\ddot{a}-{\displaystyle \frac{4\pi Ga}{3}}\left(\rho +3p\right)-{\displaystyle \frac{2}{3{\ell }_g^2}}a{\partial }_t\varphi \right)+a{\partial }_t\varphi -\dot{a}\varphi$$

(84)

only with ${\ell }_0\mapsto {\ell }_g$ multiplying the coarsening piece (due to the different scale-flows involved) and a dark energy term resulting from splitting $\frac{1}{3}a{\partial }_t\varphi$, such that for $\varphi =-t$ the scaling piece in (73) is recovered.

With the above modifications in mind, (81) becomes

$${\left({\displaystyle \frac{\dot{a}}{a}}\right)}^2={\displaystyle \frac{8\pi G}{3}}\rho -{\displaystyle \frac{2{\partial }_t\varphi }{3{\ell }_g^2}}$$

(85)

which is the first Friedmann equation with an extra term mimicking dark energy. Since (85) is satisfied (also) at $\lambda =1$, reasonably assuming that $\rho {|}_{t=t_0,\lambda =1}$ is on the order of the current baryonic density based on direct ‘count’, ${\rho }_b\approx 3\times {10}^{-19}\mathrm{kg}/{\mathrm{km}}^3$, most likely a lower bound, and $H_0=\dot{a}/a\approx 70\mathrm{km}/\mathrm{s}/\mathrm{Mpc}$ based on local measurements (validated below), an estimate ${\ell }_g\approx 2\sqrt{-{\partial }_t\varphi \left(t_0\right)}\times {10}^{23}$km is obtained, hence ${\partial }_t\varphi (t_0,\lambda =1)<0$, i.e., the ${\partial }_t\varphi$ term in (85) mimics dark energy which is currently positive.

Let us summarize the computational task of finding a solution for the relativistic cosmological model. The single scale-flow equation is for a, (84), whose solutions must be positive, not wildly diverging at large $\lambda$. Equations (85), and (79) with $a(t,\lambda )\mapsto a(t,{\lambda }^{-1})$, act as constraints, which for a given a and ${\partial }_{\lambda }a$ couple $\varphi$ and $\rho$ at any (fixed-$\lambda$) section. The Propagation of a in scale depends on p which, as in a standard Friedmann model, requires extra physical input regarding the nature of the energy-momentum tensor, e.g., an equation-of-state relating p and $\rho$. Since both $\rho$ and p represent some large-volume average of ${\mathsf{\Theta }}_{\mu \nu }+P_{\mu \nu }$, removed of $P_{\mu \nu }$’s ‘dark’ component, the contribution from inside matter (where $J_0\not\equiv 0$), denoted ${\rho }_{\mathrm{m}}$, can be assumed to be that of non-relativistic (“cold") matter, i.e., $p_{\mathrm{m}}\approx 0$. Outside matter the A-field is nearly a vacuum solution of Maxwell’s equations with an associated traceless ${\mathsf{\Theta }}_{\mu \nu }$ contributing ${\rho }_{\mathrm{r}}$ and $p_{\mathrm{r}}=\frac{1}{3}{\rho }_{\mathrm{r}}$ to the total $\rho$ and p. If we proceed as usual, identifying ${\rho }_{\mathrm{r}}$ with the energy of retarded radiation emitted by matter, observations would then imply ${\rho }_{\mathrm{r}}\ll {\rho }_{\mathrm{m}}$ in the current epoch. However, ${\rho }_{\mathrm{r}}$ incorporates also the ZPF which could potentially even outweigh ${\rho }_{\mathrm{m}}$. The contribution of the ZPF, being an ‘extension’ of matter outside the support of its $J_{\mu }$, although having a distinct $p\left(\rho \right)$ dependence, is not an independent component. Properly modeling the combined matter-ZPF fluid, e.g. as interacting fluids, or using some exotic equation-of-state, will be attempted elsewhere.

Returning to the fixed-$\lambda$ constraint, (85) and (79), we first note that $\varphi$ in the latter describes the motion of a damped, harmonic oscillator with a negative spring coefficient and a force term (whose sign and magnitude depend on ${\partial }_{\lambda }a$). A general negative solution has a single local maxima at $t_{\max }<t_0$ (by ${\partial }_t\varphi (t_0,\lambda =1)<0$). Since at a fixed $\lambda$, (79) is second order in time, only one of its two integration constants is fixed by (85) evaluated at $t_0$. The second one can then be used to further tune ${\partial }_{tt}\varphi \left(t_0\right)$ in matching the observed, current acceleration, via

$${\partial }_t\left({\displaystyle \frac{8\pi G}{3}}\rho +{\displaystyle \frac{2{\partial }_t\varphi }{3{\ell }_g^2}}\right){|}_{t=t_0,\lambda =1}={\partial }_t\left({\displaystyle \frac{8\pi G}{3}}{\rho }_m^{\mathsf{\Lambda }}\mathrm{CDM}\right){|}_{t=t_0}$$

Above, ${\rho }_m^{\mathsf{\Lambda }}\mathrm{CDM}\approx 5{\rho }_b$ is the $\mathsf{\Lambda }$CDM cold matter density estimate based on supernova- and transverse BAO-distance observations, and ${\partial }_t\rho (t_0,\lambda =1)$ is determined by (28) (evaluated at $t=t_0,\lambda =1$). Thus, the two integration constants of (79) can conspire to result in an illusion of both a positive cosmological constant, and a cold dark matter addition to $\rho$ (even if ${\rho }_{\mathrm{ZPF}}\approx 0$).

Moving to the early universe, or star-phase of the explosion, at $t<t_{\max }$, the ${\partial }_t\varphi$ term in (85) switches sign, and rapidly decreases with decreasing t, countering the opposite trend in $\rho$, dramatically slowing the shrinkage of $a\left(t\right)$, likely eliminating the horizon problem plaguing a generic Friedmann model. The precise outcome of such a battle of divergences depends on the details of a solution, but a natural, physically motivated scenario follows from the fact that, $\dot{a}=0$, $\rho ={\mathrm{const}}_1$ and $\dot{\varphi }={\mathrm{const}}_2$ is a solution for the system (79)(85) when the two constants are chosen so that the r.h.s. of (85) vanishes. Namely, the growth of ${\partial }_t\varphi$ eventually ‘catches-up’ with that of $\rho$, meaning that there is no big-bang in the remote past, just a static universe/star. During that epoch, a pertubative analysis of $\rho$, p, $g_{\mu \nu }$ and $z^{\mu }$ can be performed. A notable departure from standard such analysis is the appearance of ‘vacuum waves’, perturbative solutions of (36), with associated $\delta P_{\mu \nu }$ masquerading as dark matter of some sort.

Relating cosmological observations to $a(t,\lambda =1)$ entails extra steps which are different in the proposed formalism, therefore expected to yield different relations. Remarkably, this isn’t so in most cases. Consider, e.g., the redshift. To calculate the redshift of a distant, comoving object at $t=t_1$, two adjacent, time-ordered points along its worldline are to be matched with similar two points for earth at $t_0$. The matching is done by finding two solutions of (59) which are well behaved on all scales, satisfying the light-cone condition, connecting the corresponding points at $\lambda =1$. For the metric (75) and scaling field (78), the equation for ${\zeta }^r$ (denoting $r=\left|\mathit{x}\right|$) and ${\zeta }^t$ of each path becomes

$$\lambda {\partial }_{\lambda }{\zeta }^r={\ell }_0^2\left(-{\partial }_{\tau \tau }{\zeta }^r-2{\displaystyle \frac{\dot{a}}{a}}{|}_{{\zeta }^t}{\partial }_{\tau }{\zeta }^t{\partial }_{\tau }{\zeta }^r\right)+\tau {\partial }_{\tau }{\zeta }^r$$

(86)

$$\lambda {\partial }_{\lambda }{\zeta }^t={\ell }_0^2\left(-{\partial }_{\tau \tau }{\zeta }^t-a\dot{a}{|}_{{\zeta }^t}{\left({\partial }_{\tau }{\zeta }^r\right)}^2\right)+\tau {\partial }_{\tau }{\zeta }^t-\varphi ({\zeta }^t,\lambda )$$

(87)

subject to the light-cone condition. The two adjacent solutions at $\lambda =1$, indexed by $\alpha$ (earlier) and $\beta$, trace trajectories $r_{\alpha }\left(t\right)$, $t\in \left[t_1,t_0\right]$ and $r_{\beta }\left(t\right)$, $t\in \left[t_1+{\delta }_1,t_0+{\delta }_0\right]$, and the redshift is calculate from the equality

$$r_1-r_0={\int }_{t_1}^{t_0}{\dot{r}}_{\alpha }\left(t\right)\underset{\dot{}}{\mathrm{t}}={\int }_{t_1+{\delta }_1}^{t_0+{\delta }_0}{\dot{r}}_{\beta }\left(t\right)\underset{\dot{}}{\mathrm{t}}$$

(88)

as $z={\delta }_0/{\delta }_1-1$. Now, on the two, non overlapping, $O\left({\delta }_0\right)$ parts of their supports, $r_{\alpha /\beta }$ (almost) satisfy the light-cone condition which, for the highly symmetric metric (75) implies ${\dot{r}}_{\alpha }={\dot{r}}_{\beta }=a^{-1}$. Assuming that a changes very little over ${\delta }_0$, the difference between the two integrals in (88) coming from those end parts is

$${\displaystyle \frac{{\delta }_0}{a\left(t_0\right)}}-{\displaystyle \frac{{\delta }_1}{a\left(t_1\right)}}$$

(89)

which would give the standard expression for the redshift in terms of a. However, in addition there is also an $O\left({\delta }_0\right)$ contribution from the overlap

$${\int }_{t_1}^{t_0}{\displaystyle \frac{\mathrm{d}}{\underset{\dot{}}{\mathrm{t}}}}\delta r\left(t\right)\underset{\dot{}}{\mathrm{t}},\delta r:=r_{\alpha }-r_{\beta }$$

(90)

(rounding the boundary points for clarity which is legitimate to leading order), as the two solutions for (86)(87) see slightly different Hubble parameters, on the order of $\dot{H}\delta {\zeta }^t$, and a slightly different scaling field in (87), $\dot{\varphi }\delta {\zeta }^t$. Nonetheless, since $\delta r$ vanishes at its endpoints, contribution (90) also vanishes. Note that only the light-condition entered the above analysis, rather than the explicit, conjectured (86) and (87). Similarly for the angular diameter distance and the luminosity distance, the latter further requiring exact conservation of ${\mathsf{\Theta }}_{\mu \nu }$, which is true also in our model.

Finally, the flatness problem of pre-inflation cosmology is elegantly dismissed as follows. First, generalizing the relativistic model to a curved-space FLRW metric is straightforward, and the Friedmann equation (85) receives a $-k/a^2$ addition to its r.h.s. Denoting the ratio between its k and $\rho$ r.h.s. terms, $\chi (t,\lambda ):=3k/\left(8\pi G\rho a^2\right)$, the flatness problem can be stated as the “unrealistic" fine-tuning of $\chi$ to near mathematical zero at early times, needed to bring its current, observed value to zero within measurement uncertainties. For example, if $\rho (t,\lambda )={\rho }_0\left(\lambda \right)a^{-3}(t,\lambda )$, with ${\rho }_0\left(\lambda \right)$ encoding the creation/annihilation of matter in scale flow, then $\chi \propto a/{\rho }_0$ which, at a fixed $\lambda$, grows by many orders of magnitude over the history of the universe. However, in our formalism the universe is not a machine, propagating in time its state at an earlier time, as previously explained in the context of Bell’s theorem; Friedmann’s equation (85) enters the relativistic model as a constraint, not an evolution rule, and a cosmological solution is just what emerges out of the set of all constraints. Moreover, even when seen as an evolution rule, (85) may lead to the following counter argument: At a fixed time, a reasonable ${\rho }_0\left(\lambda \right)$ interpolating between a star and a ‘universe’ would counter the growth of a in the $\lambda \downarrow$ direction (i.e., the scale-rate of density growth due to matter creation is greater than the third root of its geometric depletion rate). Thus, unless $\chi$ is fantastically large close to the scale, ${\lambda }_{\mathrm{giant}}$, at which a giant’s section corresponds to a star—and why should it be?—$\chi (t_0,{\lambda }_{\mathrm{human}})=0$ to within measurement uncertainties is a perfectly “realistic".