Black Holes and Kullback-Leibler Divergence, Decomposing Path-Dependent Processes

1. Introduction

General Relativity achieved remarkable success in predicting a wide range of astronomical observations. However, there remain narrow regimes where Einstein’s theory falls short. In these domains, an Information Geometry approach offers new techniques to probe the underlying deficits. In this paper, we demonstrate that the principle of stationary action is equivalent to minimizing the Kullback-Leibler (KL) divergence in stochastic systems. This equivalence allows us to decompose complex, path-dependent distributions into symmetric sub-components. By leveraging the additivity property of divergence, we isolate an “error” distribution that captures the information where General Relativity fails to match observation.

Building on established work that connects maximum entropy with stationary action to recover dynamical laws [1,2,3], our formulation expresses the KL divergence in terms of the path-dependent distribution $P\left[\mathrm{\Gamma }\right]$, preserving the Euler-Lagrange equations of classical mechanics. Through computational and analytical methods it is shown that $\Delta D_{KL}=0$, an entropy maximum, corresponds to an ideal state where theory and observation align. Specifically, one can decompose a complex distribution into distributions that have specific properties that are invariant under specific transformations. These invariance are closely linked to conserved quantities and informational symmetries akin to Noether’s Theorem. This model-to-observation deficit is leveraged to analyze the "missing" dynamics of General Relativity through deriving a measure for multi-path dependence. Finally, we derive a novel relationship between entropy, energy, and the expansion scalar:

$$\gamma {\displaystyle \frac{TdS+SdT}{dE}}={\displaystyle \frac{dV}{d\tau }}$$

(1)

which offers fresh insight and a more intuitive understanding into black hole relaxation and viscosity.

2. Methods and Established Theory

2.1. Review of Stochastic Action Principle (SAC)

The Stochastic Action Principle (SAP) provides us with a framework whereby the principle of stationary action is equivalent to maximized entropy (Wang, 2009). First hypothesized by Jaynes (1957)[4] Principle of Maximum Entropy, this idea has been built up by various researchers over the years. However, for the benefit of brevity, we will concentrate on formalisms and experimentation provided by Dr. Wang’s Stocastic Action Principle derived from Shannon entropy, and Dr. Yahalom’s (2022) work starting from the Fisher Information Perspective. To start we will have path information encoded as Shannon entropy in relation to the path probability distribution. Specifically, rather than looking at a singular path we instead look at the average path taken by a large number of particles. For entropy S, Path $\mathrm{\Gamma }$ and distribution P:

$$S\left(\mathrm{\Gamma }\right)=-{\Sigma }_{\gamma }P\left[\mathrm{\Gamma }\right]lnP\left[\mathrm{\Gamma }\right]$$

(2)

The Action (A) Value measured via the Action Functional:

$$A\left[\mathrm{\Gamma }\right]={\int }_{t_a}^{t_b}L(x\left(t\right),\dot{x}\left(t\right),t)dt$$

(3)

then determine the path ($\mathrm{\Gamma }$) by maximizing entropy subject to system knowledge constraints, average action:

$$<A>={\Sigma }_{\mathrm{\Gamma }}P\left[\mathrm{\Gamma }\right]A\left[\mathrm{\Gamma }\right]$$

(4)

The entropy maximization framework is seen via normalization, with fixed <A>:

$$\delta (S-\eta {\Sigma }_{\mathrm{\Gamma }}P\left[\mathrm{\Gamma }\right]A\left[\mathrm{\Gamma }\right]-\alpha {\Sigma }_{\mathrm{\Gamma }}P\left[\mathrm{\Gamma }\right])=0$$

(5)

The Lagrangian for variation, functional J extremized with bounding constants $\alpha ,\eta$:

$$J=-S+\alpha ({\Sigma }_{\mathrm{\Gamma }}P\left[\mathrm{\Gamma }\right]-1)+\eta ({\Sigma }_{\mathrm{\Gamma }}P\left[\mathrm{\Gamma }\right]A\left[\mathrm{\Gamma }\right]-<A>)$$

(6)

Take variation with respect $P\left[\mathrm{\Gamma }\right]$ for small variation.

$${\displaystyle \frac{\partial J}{\partial P\left[\mathrm{\Gamma }\right]}}=0=>-(lnP\left[\mathrm{\Gamma }\right]+1+\alpha +\eta A\left[\mathrm{\Gamma }\right])=0$$

(7)

or

$$lnP\left[\mathrm{\Gamma }\right]=-(1+\alpha )-\eta A\left[\mathrm{\Gamma }\right]$$

(8)

and $1+\alpha$ is constant so exponentiation both sides gets.

$$P\left[\mathrm{\Gamma }\right]={\displaystyle \frac{1}{Q}}{e}^{-\eta A\left[\mathrm{\Gamma }\right]}\propto e^{-\eta A\left[\mathrm{\Gamma }\right]}$$

(9)

This is the final path probability distribution, more precisely, the probability of a path decays exponentially with its action deviation. In discrete notation for paths k from a to b:

$$p_{ab}\left(k\right)={\displaystyle \frac{1}{Q}}{e}^{-\eta A_{ab}\left(k\right)}$$

(10)

With extremal condition:

$$lnP\left[\mathrm{\Gamma }\right]+\eta A\left[\mathrm{\Gamma }\right]=const$$

(11)

This shows the balance of maximizing entropy and minimizing action cost term i.e. "Maximum Path Entropy Principle". These results were showing in a computational experiment by Wang (2009) showing a Gaussian cloud of dust particles travel through paths that obey this principle. They also find Shannon entropy:

$$S_{ab}=-\int D\left(r\right)p_{ab}ln{p}_{ab}$$

(12)

which mimics thermodynamics $dQ=dU+dW$, and $U=\overline{E}={\Sigma }_i{p}_i{E}_i$. Finally, we find:

$${\partial }^2{p}_{ab}=-\eta p_{ab}{\partial }^{2}A,\eta ={\displaystyle \frac{1}{2md}}>0$$

(13)

then $\partial p_{ab}=0$ so $\partial A=0$ so ${\partial }^{2}A>0$ or in other words, the least action path is the most probable path [5]. When we have vanishing noise, stochastic dynamics simply tend towards regular dynamics:

$$S_{ab}=lnZ+\eta \overline{A}$$

(14)

implying $\partial A=0$, our stationary action principle, a perspective supported by maximum entropy methods in statistical physics [6].

2.2. Recovering Classical and Quantum Mechanics from SAC and Fisher Information

To recover Classical Mechanics is quite trivial. In the limit where the parameter $\eta$ is large (small fluctuations or low temperature in statistical mechanics analogy), the weight $e^{-}\eta A\left[\mathrm{\Gamma }\right]$ becomes extremely small unless $A\left[\mathrm{\Gamma }\right]$ is at its minimum. This is analogous to the Feynman Path integral construction which for classical systems yields the Euler-Lagrange equations and therefore all of Classical Mechanics.

For Quantum Mechanics this is less straightforward, as the path integral weighs each path by a complex exponential:

$$A\left[\mathrm{\Gamma }\right]\propto e^{{\displaystyle \frac{i}{\hslash }}A\left[\mathrm{\Gamma }\right]}$$

(15)

For this we must allow a non-vanishing "spread" of paths, by modifying maximum entropy approach to work with complex probability amplitudes rather than probabilities in an entropic dynamics framework [7]. By incorporating interference between paths, the resulting theory reproduces the quantum dynamics governed by the Schrodinger equation. In effect, classical principle of least action is replaced by a principle where all paths contribute, weighted by $e^{{\displaystyle \frac{i}{\hslash }}A\left[\mathrm{\Gamma }\right]}$

An alternative construction has been shown to lead to the quantum Hamilton-Jacobi equation and ultimately the Schrodinger equation reproducing quantum dynamics as shown in [8] and then Dirac Equation [3].

2.3. Extending to Fisher Information Geometry

Now we will briefly review some concepts with Fisher Information Geometry which will be required for our later analysis. These are, that information geometry has a dualistic structure. In particular, two conjugate torsion-free affine connections coupled with the metric. Analogous geometric ideas also emerge in nonlinear sigma models, where metric flows are interpreted via renormalization group methods (see Friedan, 1985 [9]). Dual parallel transport that is metric-compatable, $2^k$ types of k-gons. Compared with Riemannian manifolds which exhibit opposite holonomies (net rotation of vector) i.e. Einstein Field Symmetry/Geodesic congruence. These features are not necessarily true for information manifolds when: duel affine connections exist such as (e-connect, m-connect), connections with torsion, non-metric connections (Weyl), or non-Abelian gauge field or spin structures i.e. transporting a spinor.

This paper argues that these unique properties of an information manifold actually allow us to better model spacetime dynamics, particularly Kullback-Liebler Divergence which exhibits $D_{KL}\left(P\right|\left|Q\right)\ne D_{KL}\left(Q\right|\left|P\right)$ or more precisely, not necessarily equal.

2.4. Kullback-Liebler Divergence

Divergence between two distributions P and Q is:

$$D_{KL}\left(P\right|\left|Q\right)={\Sigma }_{x}P\left(x\right)ln{\displaystyle \frac{P\left(x\right)}{Q\left(x\right)}}$$

(16)

This canonical definition and its fundamental properties are standard in information theory (see, e.g., Cover and Thomas, 2006 [10]). For appropriate distributions, we can show that available work exists P relative to Q:

$$F_{neq}-Feq=kT{D}_{KL}\left(P\right|\left|Q\right)$$

(17)

In other words the available work exists relative to equilibrium [11]. Moreover, this phase-space perspective on dissipation, which links thermodynamic work to the asymmetry between forward and reverse path probabilities, is thoroughly discussed in Kawai, Parrondo, and Van den Broeck (2007) [11]. Here $\left(Q\right|\left|P\right)$ asks how inefficient it is to use P to model Q. We can easily see that the directional imbalance of distributions go away when:

$$\Delta D=\left|D\right(P\left|\right|Q)-D(Q\left|\right|P\left)\right|=0$$

(18)

It should be intuitive to see this is a maximal entropy state, there are more microstates available when $D\left(P\right|\left|Q\right)=D\left(Q\right|\left|P\right)$, this is already well established as a consequence of:

$$D\left(P\right|\left|u\right)={\Sigma }_{xϵX}P\left(x\right)log{\displaystyle \frac{P\left(x\right)}{u\left(x\right)}}=log\left|X\right|-H\left(P\right)$$

(19)

This is analogous to Ficks Normal Diffusion (proportional to negative concentration gradient). For any process that distorts distributions P or Q, where P and Q are dependent distributions then there is a tendency to maximize entropy (i.e. moving from high-density to low-density states). This emerges from the balanced Log Ration Condition:

$$ln{\displaystyle \frac{P\left(x\right)}{Q\left(x\right)}}so\int (p\left(x\right)-Q\left(x\right))ln{\displaystyle \frac{P\left(x\right)}{Q\left(x\right)}}=0$$

(20)

2.5. Computing Divergence in identifiable Distributions

There is a rich body of literature which shows various distributions, have this property. For example, Gaussian distributions exhibit rotational invariance—reflecting the conservation of angular momentum (see, e.g., [10])—and location shift families demonstrate translation invariance, which underpins momentum conservation; similarly, the Cauchy distribution’s reflection and scale invariance, along with its invariance under Möbius transformations, mirrors conformal symmetries as discussed in [9], while nonequilibrium thermodynamic approaches have linked path-dependent entropy production to these symmetry properties ([11,12]); furthermore, the axial symmetry of von Mises–Fisher distributions ([13]) and the dual affine structure inherent to Bregman divergences ([14]) further corroborate how such statistical symmetries are intimately connected to fundamental conservation laws. These distributional properties can all be verified by a fairly trivial computational experiment. The resulting distributions which produce $\Delta D=0$ all have symmetrical properties which can be directly related to invarience (conservation under translation). From a range of distributions analyzed, below are some notable distributions whereby $\Delta D=0$ (and properties to note):

• Gaussian distributions (same covariance, frame of reference invariance, related to time symmetry),
• location shift families (integrated depends on ${\theta }_1-{\theta }_2$, related to translation symmetry),
• Canchy distributions of same family (same scale, Reflection, Scale Invariant, some Mobius transformations),
• Antipodal or Permuted Discrete Distributions (mirror image or group rotation, rotational symmetry),
• Von Mises-Fisher (preferred direction spin up/down, symmetric under mean direction but not concentration),
• Bregman Divergence (exponential family distributions are symmetric when "distance" d is essentially Euclidean, implying no extra gauge or affine structure, also implies connection used is usual (Levi-Civita) (see, e.g., Amari and Nagaoka, 2000 [?])).

• This should not be too surprising, the distributions on an information manifold can be viewed essentially as field states.

2.6. Decomposing Observed Distributions

The obvious next step is to examine an observed distribution signal, which may appear random or chaotic can actually be represented by a series of decomposable symmetric distributions. Given this, we look towards methods to linearly separate a mixed distribution into their component distribution parts. This again is fairly straightforward using Fourier analysis and signals analysis techniques. In particular, we can find the marginal likelihood for a given model $M_i$ to be a fit within a mixed solution from:

$$p\left(y\right|M_i)=\int p\left(y\right|{\theta }_i,M_i)p\left({\theta }_i\right|M_i)d{\theta }_i$$

(21)

The posterior probability for model $M_i$ from Bayes Rule:

$$p\left(M_i\right|y)={\displaystyle \frac{p\left(y\right|M_i)p\left(M_i\right)}{p\left(y\right)}}$$

(22)

and from Bayes Factor (see, e.g., Kass and Raftery, 1995 [15]):

$$B{F}_{ij}={\displaystyle \frac{p\left(y\right|M_i)}{p\left(y\right|M_j)}}$$

(23)

Alternative checks for likelihood, for one distribution to be the "correct" decomposition can be found using AIC/BIC or Penalized Likelihood/MAP estimation instead or other techniques. The important thing in this analysis is to find a high degree of confidence that a given signal is decomposed correctly into its sub-signals, and if uncertainty exists between configurations this must be factored into analysis. To better avoid this, we can determine boundary conditions based on physical law. Further confidence is increased when evaluating the time evolution of the distribution as each composable signal has a unique time evolution signature when it is perturbed and drifts back to its idealized state.

2.7. Complex Components (Phase Space Projection)

If we have a distribution that exists in phase space, but its distribution can only be observed in $R^3$, then we inevitably get complex structure and unpredictability. The result of this is a large variation in the probability "certainty" of the linear separation. In other words, a single smooth distribution may appear in $R^3$ as having multiple modes. If we are optimizing for $R^3$ space then $D\left(Q\right|\left|P\right)$ again will exhibit mode-seeking behavior, where under time evolution the divergence of $\left(Q\right|\left|P\right)$ will re-center itself around one mode which acts pseudo "orthogonal" to the leftover modes. This mode-seeking behavior is conceptually very similar to the super-position of states. Selecting one distribution automatically causes cascading effects into all remaining distribution patterns, mimicking wave function collapse. Similarly, we can only capture one mode at a time, so the higher the convergence on one mode, the less information and therefore higher uncertainty we have about the remaining, unselected mode. For a complex function

$$f\left(z\right)=u\left(x_{\mu }\right)+iv\left(x_{\mu }\right)$$

(24)

Projected onto the real axis we define projection:

$$g\left(x_{\mu }\right)=\mathrm{Re}\left(f\left(x_{\mu }\right)\right),h\left(x_{\mu }\right)=\mathrm{Im}\left(f\left(x_{\mu }\right)\right)$$

(25)

Complex function is simply a sum of exponentials, real projection exhibit interference patterns:

$$f\left(z\right)={\Sigma }_n{e}^{i{k}_{n}z},\mathrm{Re}\left(f\left(z\right)\right)={\Sigma }_{n}cos\left(k_{n}x\right)$$

(26)

The complex Fourier modes is simply:

$$f\left(z\right)={\Sigma }_n{A}_n{e}^{i{k}_{n}z}$$

(27)

An interesting pattern emerges whereby for any complex distribution, the real distribution must have its most likely decomposition to include multiple sub-distributions. When this pattern is renormalized, the distribution not captured by mode optimization is free for further signal decomposition. The reader may note this is very similar to the way the Standard Model has essentially been constructed:

$${\mathcal{L}}_{SM}={\mathcal{L}}_{gauge}+{\mathcal{L}}_{fermion}+{\mathcal{L}}_{Higgs}+{\mathcal{L}}_{Yukawa}$$

(28)

Where each Lagrangian term represents its own path and action minimization and therefore the fields can be written in terms of a path distribution. Conceptually we can think of this process in terms of entropic maximization tends towards phase alignment

2.8. Path Dependence

Instead of just a simple distribution of point P(x), the probability of being in state x. We can consider the probability of observing a particular path $P\left[\mathrm{\Gamma }\right]$, where $\mathrm{\Gamma }$ can be either discrete steps in a path or a continuous time evolution. Therefore, path-dependent divergence with the path integral measure $\beth \mathrm{\Gamma }$, defined as the limit of the product measure over all discretized time steps of the path with appropriate normalization is:

$$D_{KL}\left(P\right|\left|Q\right)=\int \beth \mathrm{\Gamma }P\left[\mathrm{\Gamma }\right]ln{\displaystyle \frac{P\left[\mathrm{\Gamma }\right]}{Q\left[\mathrm{\Gamma }\right]}}$$

(29)

This tells us how surprised we would be for a path to deviate from assumed dynamics Q, rather than the observed dynamics P. Entropy production relative to the ratio of forwards and backwards path probabilities, total entropy production is given by the KL divergence [16]. Where $\left[\tilde{\mathrm{\Gamma }}\right]$ is the time-reversed trajectory:

$$\Delta S_{tot}\left[\mathrm{\Gamma }\right]=ln{\displaystyle \frac{P\left[\mathrm{\Gamma }\right]}{Q\left[\tilde{\mathrm{\Gamma }}\right]}}$$

(30)

And the average over all trajectories in a stochastic system is:

$$\Delta S_{tot}=D_{KL}\left(P\left[\mathrm{\Gamma }\right]\right||Q\left[\tilde{\mathrm{\Gamma }}\right])$$

(31)

Similar relations for nonequilibrium steady-state entropy production have been derived by Hatano and Sasa (2001) [12], further reinforcing the connection between path-dependent dynamics and irreversible work. This average KL Divergence quantifies total entropy production for traveling a path, giving us the irreversibility of a process/dynamic (cf. Seifert, 2012 [17]). This in "reverse" $D\left(Q\right[\mathrm{\Gamma }\left]\right|\left|P\right[\mathrm{\Gamma }\left]\right)$ quantifies the error in observed model when trying to capture idealized measures. This is akin to saying what is the system inefficiency or signal loss.

2.9. Identifying Multi-Path Dependence

The dual affine structure of Information Manifolds provides the perfect grounding for evaluating processes with multiple separable path-dependent processes. It is simple to find examples, an object placed in front of a fluid flow will experience path-dependency from the geometry of the space (assuming constant curvature), while the displacement of fluid flow due to relative differences in inertia result fluid being forced around the object in turbulence. This is consistent with the thermodynamic metrics approach [18]. In this, not only is the object affected by the flow geometry, but the geometry itself is affected by the presence of the object. If we remove this object, its impacts can still be seen until equilibrium is re-established and normal, unimpeded flow returns to an equilibrium. In this, total divergence can be separated by Leibniz rule or additivity of KL divergence, so we can decompose these two separate but related path-dependent processes into:

$$D_{KL}\left(P\right|\left|Q\right)=D_{KL}\left(P\right||P_1)+D_{KL}(P_1\left|\right|Q)+{\Sigma }_x[P\left(x\right)-P_1\left(x\right)]ln{\displaystyle \frac{P_1\left(x\right)}{Q\left(x\right)}}$$

(32)

if $P_1\left(x\right)/Q\left(x\right)$ is constant i.e. when $P\left(x\right)\ne P_1\left(x\right)$, $P_1\left(x\right)=0$ or if $P\left(x\right)=P_1\left(x\right)$ then this cross term vanishes [10]), which is the case if an observed distribution can be perfectly modeled by two separate path models.

If a composite distribution can be viewed as a mixture or sum of sub-distributions (A + B = C in some sense), then its total divergence with respect to baseline Q generally decomposes into the sum of the divergences of the sub-distributions, plus additional cross terms. Those cross-terms vanish under specific conditions. In this treatment we assume symmetry in the sense that their density ratios relative to Q are constant (or nearly so) over their support. Under these conditions, the cross-terms in the full decomposition vanish, justifying the simplified additive form: $D_1=D_2$ and $\Delta D_{12}=D_1-D_2$ relating them to baseline distribution Q, final distribution P and an intermediate distribution $P_1$. As such for a distribution $\left(P\right|\left|Q\right)$ comprised only of symmetric distributions $D_1=(P_1\left|\right|Q)$, and $D_2=\left(P\right||P_1)$ respectively:

$$\Delta D_{KL}(P,Q)=\Delta D_{KL}(P_1,Q)+\Delta D_{KL}(P,P_1),\Delta D_1=-\Delta D_2$$

(33)

Given the maximum entropy is found at $\Delta D_{KL}=0$ then if (P,Q) are in a state of maximum entropy, then its two sub-components must either be reciprocals or equal (i.e, symmetric themselves). If they are independent paths, they should tend towards individual net divergence of zero, but small perturbations essentially are guaranteed to be antisymmetric in a sense. This should make sense, the amount of differences between two composite distributions results in the divergence of the underlying total distribution.

This is very powerful, without requiring any information about a missing path process, we can say that we can measure the energy gap between a theoretical model, and our observed distribution. From this we can gain insight or even potentially recover the path information which will give rise to the underlying mechanisms that our theoretical models do not currently capture.

2.10. Quantum Reference Frame Invariance

The reader might pause here because $D_{KL}$ can be related to entropy, and in this time evolution we have $\partial D_1=-\partial D_2$ this seams unphysical, we should not be seeing a negative path derivative of entropy. However, it is already established that entanglement and coherence of subsystems is invariant under ideal Quantum Reference Frame transformations (Cepollaro et al., 2024) [19]:

$$E^A+C^A=E^B+C^B,\Delta E\Delta C\le 0$$

(34)

Additionally, Vopson & Lepadau (2022) found what they coined "the Second Law of Infodynamics"[20], where under repeated observation and experimentation, across vastly different systems they observed an informational entropy $S_{inf}$ with exactly this property, in particular, that $S_{total}=S_{trad}+S_{info}$ or written differently, $\partial S_{trad}=-\partial S_{inf}$. These two results one would argue exactly meets the condition set for $\partial D_1=-\partial D_2$ as an entropic maximum.

Carlo et al.’s Theorem 1 demonstrates that for pure bipartite quantum states, the sum of a particular entanglement measure and its corresponding coherence measure remains invariant under ideal quantum reference frame (QRF) transformations. Although this result is derived for a simple scenario using specific measures, specifically a pure bipartite quantum states, its core message is that a unitary QRF change merely redistributes quantum resources without altering the overall quantum uncertainty of the system.

In our work, we invoke this theorem solely as an illustrative analogy rather than as a foundational assumption. We explicitly acknowledge its limitations—namely, that it applies only to pure bipartite states and may not hold for more complex or mixed systems. Nevertheless, the underlying principle that total “quantumness” (or quantum uncertainty) is conserved under a change of reference frame is well grounded in the unitarity of quantum mechanics. This trade-off has been noticed in other areas of physics such as the "second law of infodynamics" whereby a system’s informational content decreases with an increase in information. In this they see again a total information-centric conservation law arise where the total informational content within the system is preserved and transformed in a way that preserves this relationship.

Furthermore, as discussed in our Appendix A, we extend this intuition by considering the trade-off between local coherence and system–environment entanglement in dissipative settings. In this extended framework, even if local measures (such as coherence) appear to decay, the overall quantum uncertainty remains invariant when accounting for all parts of the system. While this generalization is speculative, it provides a plausible conceptual bridge from the restrictive conditions of Theorem 1 to the broader phenomena observed in our study.

In summary, although Carlo et al.’s result is limited in scope, its central insight—that the total quantum resources are preserved under ideal QRF transformations—supports our view that decaying informational entropy does not imply a loss of overall quantumness, but rather a redistribution that is frame-independent.

3. Results

It should be clear that, $D_{KL}\left(P\right|\left|Q\right)=0$ occurs when our theoretical model matches observation, i.e. when an ideal physical model is known. $D_{KL}(P_1\left|\right|Q)$ is how much the ideal theoretical model is mismatched by our current model. $D_{KL}\left(P\right||P_1)$ tells us how much the amount of work required to bring a universe following our best current model to match our observed universe. In this idealized case, $D_{KL}\left(P\right|\left|Q\right)=0$, the following must result $D(P_1\left|\right|Q)=-D\left(P\right||P_1)$. This can be described as a decomposition of two separate, but additive types of path dependence (i.e. geometric curvature and drag), but ultimately, knowing the specific underlying mechanism is irrelevant. This is because the optimized feature is an energy deficit from our expected distribution from our models and the observed universe. Lastly, it is well known that Ricci Flow follows as a Lie derivative with respect to a perturbed state. The Linearized operator governing evolution of h(t) has only smooth decaying modes, consistent with the below work. We attempt to reconcile this observational deficit by modeling it as a function of Ricci flow as a spacetime reparative process ensuring a smooth manifold.

3.1. Ricci Flow as Repair

The geometric evolution for a Riemannian metric, the unnormalized flow is given by:

$${\displaystyle \frac{\partial }{\partial t}}{g}_{ij}=-2Ri{c}_{ij}$$

(35)

We can examine the rate of "repair" compared to average flow, this is examined later as a function of temperature and density. Setting aside this for now we know the DeTurck Trick, whereby weakly parabolic (freedom to change coords) adds a diffeomorphism term to flow, leading us to a strickly parabolic system:

$${\displaystyle \frac{\partial }{\partial t}}{g}_{ij}=-2Ri{c}_{ij}+{\mathcal{L}}_V{g}_{ij}$$

(36)

In this, the Lie derivative with respect to a perturbed field V, the linearized operator governing the evolution of h(t) only has decaying modes. This ensures that any perturbation dies out over time in a smooth and exponential rate relative to neighbouring states [?]. In many cases the linearized analysis shows that each mode of the perturbation decays independently and exponentially as well. This is to say in the language of this framework that if h(t) can be decomposed into eigenmodes of the linearaized operator, then each mode decays like:

$$e^{-\lambda t}$$

(37)

with $\lambda >0$ as t approaches infinity and the metric returns (up to a diffeomorphism) to $g_0$. It is important to note however that there may exist upper and lower bounds of decay rate. If this is the case this signature should relate to a rate of information transfer (perhaps quantum speed limit) and be extracted from examining signatures upon decomposition.

3.2. Black Holes

In order to investigate the "missing dynamics" between our observations and ideal models, we can first consider small perturbations of a background metric and that this relaxation process is governed by a linearized Ricci flow relational to the linearized Ricci operator $\mathcal{L}$

$${\displaystyle \frac{\partial h_{ij}}{\partial t}}\approx -\mathcal{L}{h}_{ij}$$

This gives a natural exponential decay such that:

$$h\left(t\right)\approx e^{-\lambda t}$$

In this "Ricci repair" refers to the natural relaxation process of spacetime geometry under Ricci flow, whereby perturbations do not immediately revert to the idealized state but instead decay exponentially over time. we know that black holes can be described in three separate ways consistent with the four laws of black hole mechanics [21]. First, in terms of their event horizon; Second, in terms of the singularity; Third, in terms of the gradient curvature of spacetime between regions. This paper argues that solely examining a black hole from only one of these frameworks is an insufficient description. As such we will relate spacetime compression of states to the hypothesized upper rate limit of information transfer (quantum speed limit) as a function of an upper boundary of decay. By doing so we have an emergent boundary relational to the Bousso or Bekenstein Bound for a given observer’s reference frame. Specifically, for an outside observer, you end up with an energy density which approaches infinity at the singularity, but with maximum entropy boundaries across a finite surface area dependent on the radius of the event horizon. Each of these entropic shells have their own energy distribution as given by Helmholtz Free Energy equations. As the radius of a black hole is directly proportional to its mass, which is then proportional to the number of maximally defined states (atomic states) we can understand the reason why smaller black holes evaporate faster if entropy flow dominates energy dissipation. This is given because if $R_s\sim r$ then higher-order multipoles lead to fields dissipating far faster. However, for $R_s\gg r$ it would appear as a point source and the inverse-square law would dominate.

As such we have two reliable boundary conditions for entropy, particularly as it relates to curvature. Specifically, the maximum entropy state, or curvature state that blocks photon escape (i.e. Bousso, Bekenstein Bounds at event horizon). Assuming a stable system, all free energy approaching the event horizon is drawn into the center of the black hole rapidly until maximum entropy density is reached (finite or infinite), if the slope of the gradient of curvature is continuous between these two points, the rate of black hole decay should directly relate to differences in energy states between curvature steps. This is consistent with our expectation that Ricci Flow is the potential description for geometric repair as exponential decay. For a given system we have:

$$S\left(n\right)=\left({\displaystyle \frac{1}{C}}\right)S(n+1)$$

(38)

which can be expressed by the curve:

$$S\left(x\right)=S_0{e}^{-xlnC}$$

(39)

This "Ricci repair" mechanism, whereby a small deviation in the metric decay exponentially over a characteristic time scale, mirrors relaxation processes in thermodynamics and provides a natural explanation for the gradual repair of spacetime geometry. Below we will investigate if this an appropriate explanation for where observation does not meet our current theories.

However, because the bounding values of curvature are predefined (escape curvature of photon and maximum relative curvature between neighbour states), the only varying parameter is the radius of the black hole event horizon, and therefore we can write that the change in field decay is moderated by the scaling factor k (a function itself of r):

$$C^{\prime }={\displaystyle \frac{C}{k}},C^{\prime }k=C$$

(40)

where C represents the maximum possible field decay. Again, we can use the Von Neumann entropy S of a quantum state $\rho$ given by:

$$S\left(\rho \right)=-\mathrm{Tr}(\rho log\rho )$$

(41)

If we consider two states then the amount of entropic work performed is related to the difference in their entropies and respective energies.

$$F\left(\rho \right)=E\left(\rho \right)-TS\left(\rho \right)$$

(42)

What is interesting here is the equivalence of $TS\left(\rho \right)$ to energy, indicating that given the relatively homogeneous temperature within a stable black hole, the change of temperature relative to energy is $dT/dE\propto 1/E^2$; therefore entropy changes are equivalent to changes in energy at a first-order approximation. This should not be too surprising as a commonly defined macroscopic description of temperature as a comparison of two heat transferring states (V, N) is:

$$T:={\left({\displaystyle \frac{dU}{dS}}\right)}_{V,N}$$

(43)

In this, there are some specific consequences particularly as it relates to the thermodynamic limit for considering temperature in an information context. This will be examined in a later paper, so for now we can focus on the Work extracted which is related to the difference in free energies:

$$W=F\left({\rho }_1\right)-F\left({\rho }_2\right)$$

(44)

This can instead be looked at directly from the perspective that a change in entropy is a change in energy state. Specifically, decreasing the entropy increases free energy, which can then propagate away. Along the event horizon, very high entropy inhibits the release of photons, as this entropy dissipates away. In this case, even if $E\left(\rho \right)$ within the event horizon (outside singularity) falls to zero, i.e. in the state where all mass/energy falls towards the center, there still exists a flow energy resulting from the flow of entropy. This should make sense as the entropy flowing away from the singularity is equal to the entropy flowing away from the event horizon. In this, the flow of entropy is what makes black holes evaporate as vacuum entropy density is significantly lower. This requires its own gradient resulting in the fluid-like nature of the surface of a black hole. We can look at the rate of volumetric contraction of the black hole as it relates to the loss rate. The ratio of the rate of entropy loss S as it relates to field decay with the respective volume change defined by the total entropy/energy within the sphere. Under constant temperature:

$${\displaystyle \frac{\frac{TdS+SdT}{dt}}{\frac{dV}{dt}}}={\displaystyle \frac{\frac{TdS}{dt}}{\frac{dV}{dt}}}={\displaystyle \frac{TdS}{dV}}$$

(45)

and this from the above construction is simply:

$${\displaystyle \frac{TdS}{dV}}={\displaystyle \frac{3Cr}{4\pi r^3}}={\displaystyle \frac{3C}{4\pi r^2}}$$

(46)

The Schwarzschild radius is given by Black hole thermodynamics was first formulated by Bekenstein (1973) [22] and quantum evaporation by Hawking (1975) [23]:

$$R_s={\displaystyle \frac{2GM}{c^2}}$$

(47)

Which we can rearrange to convert back into mass:

$${\displaystyle \frac{1}{R_s}}={\displaystyle \frac{c^2}{2GM}}$$

(48)

This gives us the rate of entropy change as a ratio of total volume change (curvature change):

$${\displaystyle \frac{TdS}{dV}}={\displaystyle \frac{3C{c}^4}{16\pi G^2{M}^2}}$$

(49)

The total entropy of a black hole is given by the Bekenstein-Hawking formula and is proportional to the area of the event horizon:

$$S={\displaystyle \frac{k_B{c}^3}{\hslash G}}{\displaystyle \frac{A}{4}}$$

(50)

For which we find the mass loss rate derived from the power equation:

$${\displaystyle \frac{dM}{dt}}\sim -{\displaystyle \frac{\hslash c^4}{G^2{M}^2}}$$

(51)

And we can factor out the mass loss:

$${\displaystyle \frac{TdS}{dV}}={\displaystyle \frac{dM}{dt}}\left({\displaystyle \frac{3C}{16\pi \hslash G}}\right)$$

(52)

Since the rate of entropy loss is the same as the rate of mass loss (energy/mass equivalence, see Appendix A) we now work towards a cosmological constant for the maximum field decay rate between two states.

$$C={\displaystyle \frac{16\pi \hslash G}{3}}$$

(53)

where ${\mathbb{C}}_3={\displaystyle \frac{V}{r^3}}={\displaystyle \frac{4\pi }{3}}$ is a generalization for an n-dimensional sphere of 3 spatial dimensions. The generalized equation for the coefficient of an n-dimensional sphere (n-ball) with radius r is:

$$C_{n\left(\mathrm{dimensions}\right)}={\displaystyle \frac{V_n\left(r\right)}{r^n}}={\displaystyle \frac{{\pi }^{\frac{n}{2}}}{\mathrm{\Gamma }\left(\frac{n}{2}+1\right)}}$$

(54)

So, for a spacetime with dimensions D = (1 + 3), we have:

$$C=D{\mathbb{C}}_3\hslash G$$

(55)

This can be factored back into the exponential decay equation to get the mass decay rate due to entropy dilution.

• Similarly, we can replace mass conversion here with energy to get:

$${\displaystyle \frac{TdS}{dV}}={\displaystyle \frac{d{E}_t}{dt}},\mathrm{bounded}\mathrm{by}\text{:}{C}^u\left({\displaystyle \frac{m^3}{s}}\right)$$

(56)

where adjusting C to relativistic SI units:

$$C^u={\displaystyle \frac{C}{c^2}}={\displaystyle \frac{D{\mathbb{C}}_3\hslash G}{c^2}}$$

(57)

This is expressed in terms of SI units $({\mathrm{m}}^3/\mathrm{s})$. Specifically, it appears to be a maximum expansion rate on the order of one cubic Planck length per Planck time which is a startling, if not intuitive result. Additionally, in evaluating time as proper time $\left(\tau \right)$ we can see the rearrangement of the equation to the form of:

$${\displaystyle \frac{dt}{d\tau }}{\displaystyle \frac{TdS}{d{E}_t}}=\gamma {\displaystyle \frac{TdS}{d{E}_t}}={\displaystyle \frac{dV}{d\tau }}$$

(58)

expansion scalar change is a measure of relative relaxation over time is respective to the change of entropy in respect to the energy of the system more concretly:

$${\displaystyle \frac{dV}{d\tau }}=\theta V$$

(59)

with $\left(C_e\right)$ as the maximum energy relaxation relative to the observer i.e.

• Formulated relativistically with the Lorentz factor:

$${\gamma }_i{\displaystyle \frac{TdS}{dE}}={\displaystyle \frac{dV(\gamma i,j,k)}{d\tau }}$$

(60)

or more simply in a truly covariant solution:

$${\displaystyle \frac{TdS}{dE}}={\displaystyle \frac{dV}{d\tau }}$$

(61)

What we get from integration as a result is:

$${\displaystyle \frac{TS}{E_t}}=\tilde{V}\left(\tau \right)={\displaystyle \frac{V\left(\tau \right)}{V_0}}=e^{{\int }_0^{\tau }\theta \left({\tau }^{\prime }\right)d{\tau }^{\prime }}$$

one might be interested, particularly if they can observe free propagating energy that:

$${\displaystyle \frac{TS}{E_t}}=1-{\displaystyle \frac{E_f}{E_t}}$$

which relates to the neighboring relaxation rate that is bounded by $C^{\mu }$

$$S\left(x\right)=S_0{e}^{-xln{C}_e}$$

(62)

It is important to note however that this Lorentz factor does cancel out, given that the scaling factor $\theta$ is measured with its own Lorentz factor, giving us a truly covariant solution. This is interesting because S and E can both be described in terms of state functions which can be expanded upon from the completed work in Section 2.3:

$${\displaystyle \frac{S}{E_t}}={\displaystyle \frac{-Tr\left(\rho log\right(\rho \left)\right)}{Tr\left(\rho \widehat{H}\right)}}$$

(63)

And if $\rho$ is thermal then $\rho =-e^{-\beta H}/Z$: ${\lambda }_i=e^{-\beta E_i}/Z$ where $E_i$ is the energy of the i-th eigenstate of H and:

$$T={\displaystyle \frac{1}{k_B\beta }}$$

(64)

Or in other words, the change in symmetry compression is relative to the observer. Specifically, we can see dE as a change to the states’ eigenvalues and dS as a change to the states’ distribution, particularly as it relates to normalization. This implies that a component of curvature relaxation directly relates to the emergence of symmetries in the state function path. Again, it is obvious how this relates back to signals processing theory and in particular data modulation.

This should make intuitive sense, we can view entropy as a measure of concentration of energy on a manifold. Just like weight on an elastic surface, deformity is a function of the concentration of weight, not simply the total amount of weight on the surface. Evenly distributed weights in this case may not break the threshold where surface tension is exceeded and no deformity exists, conversely, if all the weight was concentrated in a singular point, likely the surface would tear. Temperature is a measure of the mean velocity of these particles, this too makes sense why rapid translations across a surface result in a storing of energy if the geometry is not instantaneously repaired back to idealized state in the path of the particles. It also gives a somewhat intuitive description for why particles with mass cannot pass the speed of light, because their forward deformation of spacetime would require spacetime to deform faster than the maximum limit. By adding momentum (via Lorentz factor) we are pre-deforming spacetime in front of us, through dissipation, to ensure this maximum is not reached locally. For massless particles, they do not exceed the tensile capacity of spacetime so they are free to transmit at maximum speed.

3.3. Black Hole Viscosity

Finally, the slope of the entropy flow determined by the radius of the black hole can be related to fluid dynamics where back pressure is directly related to flow in a classical system and is characterized by the slope of the fluid height profile:

$$P\left(x\right)=\rho g(h_0-Sx)$$

(65)

In this sense, it takes more work to change the flow of a system with high back pressure. This relates directly to viscosity along the event horizon. A slower entropy decay due to high back pressure would directly appear as a fluid with high viscosity. This directly matches the observed phenomena of fluid-like black hole boundaries. Further calculation is required to get a purely quantum entropy description for this back pressure.

3.4. Expansion Scalar Examination

The expansion scalar $\theta$ is used to denote fluid-like flow for an observer. It is defined as the covariant divergence of the four-velocity field:

$$\theta ={\nabla }_{\mu }{u}^{\mu }$$

(66)

Where ${\nabla }_{\mu }$ is the covariant derivative associated with the spacetime metric $g_{\mu v}$ and $u^{\mu }$ is the four-velocity field of observers or fluid elements. This expansion scalar is interpreted as expanding or diverging worldlines when $\theta >0$, contracting or convergent worldlines when $\theta <0$, and static when equal to 0. Alternatively, we find:

$${\displaystyle \frac{d\theta }{d\tau }}=u^{\mu }{\nabla }_{\mu }\theta =-{\displaystyle \frac{1}{3}}{\theta }^2-u_{\mu v}{u}^{\mu v}+{\omega }_{\mu v}{\omega }^{\mu v}-R_{\mu v}{u}^{\mu }{u}^v$$

(67)

This scalar is regularly defined as an infinitesimal volume element dV moving with the flow with observers’ proper time (see Raychaudhuri, 1955 [24]; Hawking and Ellis, 1973 [25]). The dramatically simplified version for the reader is such that:

$${\displaystyle \frac{d\theta }{d\tau }}=-{\displaystyle \frac{1}{3}}{\theta }^2+\Lambda$$

(68)

In other words $\theta$ is the fractional rate of volume change per unit of proper time. For homogeneous and isotropic universes, like as described by the Friedman-Lemaitre-Robertson-Walker (FLRW) metric, the expansion scalar is related to the Hubble parameter:

$$\theta =3H\left(t\right)$$

(69)

The behavior of $\theta$ is crucial in the Penrose-Hawking singularity theorems, with highly negative values leading to singularities, this highlights this derived result acting in congruence with the above formulation. Further, for examining tidal forces and its impact on expansion, shear and rotation congruence we can examine the evolution of the expansion scalar along the flow as governed by the Raychaudhuri equation:

$${\displaystyle \frac{d\theta }{d\tau }}=-{\displaystyle \frac{1}{3}}{\theta }^2-{\sigma }_{\mu v}{\sigma }^{\mu v}+w_{\mu v}{w}^{\mu v}-R_{\mu v}{u}^{\mu }{u}^v+a_{\mu }{a}^{\mu }$$

(70)

Implying a deep thermodynamic nature of gravity [26]. This above formulation should give ample evidence of the relationship between entropy, energy, and the expansion scalar. In particular, it should become apparent to the reader how the expansion scalar is directly related to dark energy. In particular how the cosmological constant $\mathrm{\Lambda }$ contributes to the Ricci tensor $R_{\mu v}$ such that:

$$R_{\mu v}=4\pi G{g}_{\mu v}(\rho -3p)+\Lambda g_{\mu v}$$

(71)

From the first Friedman Equation (energy density equation) we can also see:

$$H{\left(t\right)}^2={\displaystyle \frac{8\pi G}{3}}\rho \left(t\right)-{\displaystyle \frac{k}{a{\left(t\right)}^2}}+{\displaystyle \frac{\mathrm{\Lambda }}{3}}$$

(72)

and the second Friendman Equation (acceleration equation) which the $\mathrm{\Lambda }$ term drives observed accelerated expansion:

$${\displaystyle \frac{\frac{d^{2}a\left(t\right)}{d{t}^2}}{a\left(t\right)}}=-{\displaystyle \frac{4\pi G}{3}}(\rho \left(t\right)+3p\left(t\right))+{\displaystyle \frac{\mathrm{\Lambda }}{3}}$$

(73)

Also important to note that dark energy has equation state $w=-1$ which means it has negative pressure, whereby we denote pressure as having three times more effect on dynamics than energy density (hence 3p(t)).

If the overall system’s expansion scalar approaches a constant value, this reflects a state of continuous accelerated expansion otherwise seen, or more accurately previously modeled as de-Sitter expansion.

3.5. Special Conditions

This derivation results in some obvious special conditions, primarily when a system has constant entropy and energy i.e. conservation of energy and entropy. In this case, we get:

$$d\tau \propto dV\mathrm{where}\theta \mathrm{bounded}\mathrm{by}{C}_e$$

(74)

This should make intuitive sense, as time increases, the total volume enclosed in the light cone of all interactions observed within it (system evolution) increases as well. Additionally, we have another special case that is potentially worth investigating:

$$dE\propto dSd\tau \mathrm{where}\theta \mathrm{bounded}\mathrm{by}{C}_e$$

(75)

Where a local change of entropy (perturbation of a field) with discrete time steps. This describes a system of two overlapping light cones that are contained in another larger “older” light cone. As time increases, the entropy of the system increases, as such the change of energy contained within curvature also increases. Spacetime contraction requires some quanta of energy, for photons propagating this would look like a redshift as wavelength increases. It is important to note here that this phenomenon would relate directly to the gradient of neighboring states’ entropy, as a photon’s energy would be transferred into the vacuum state if the two states’ normalization is not congruent (i.e. the scalar deformation) consistent with a unitary interaction. On a universal scale, this should on average appear as a linear relationship with distance traveled, with only local deformations resulting from regions with highly similar entropic states or regions where photons gain energy from vacuum states with higher potential. Unlike other "tired light" mechanisms, as this process is based on the relative energy or entropy of the neighboring states that a particle transits through, rather than the particle’s energy itself, it should not be wavelength-dependent. As such, this process should not exhibit blurring which is a categorical problem with other "light scattering" or "tired light" models.

3.6. Temperature Regulation and Dark Matter

Another consequence of this formulation is that temperature evaluated between two points, $T=U\left(x\right)/S\left(x\right)$ is a macro, statistical description that is valid under certain assumptions. T, the average energy transmission through spacetime (i.e. average momentum of particles) may mediate the flow of the expansion scalar.

In systems with little or no relative free energy, with Von Neuman entropy, on can determine from the eigenvalue populations in the density matrix where:

$$T={\displaystyle \frac{1}{k_B\beta }},{\displaystyle \frac{dt}{d\tau }}{\displaystyle \frac{TdS}{dE}}={\displaystyle \frac{d\theta }{d\tau }}={\displaystyle \frac{dt}{d\tau }}{\displaystyle \frac{Td(-Tr(\rho log\left(\rho \right))}{dTr\left(\rho \widehat{H}\right)}}$$

(76)

and:

$$\rho ={\Sigma }_{i=1}^n{\lambda }_i|v_i{v}_i$$

(77)

so if $\rho$ is thermal, $\rho =e^{-\beta H}/Z$

and ${\lambda }_i=e^{-\beta E_i}/Z$ where $E_i$ is energy of the i-th eigenstate of H.

• In other words, average interaction frequency relates to higher rate of unitary interactions which affect the rate that a systems entropy is updated. Another way of thinking about this is it is the degree that a systems internal movement resists curvature relaxation. Intuitively this makes sense, if we repetitively block flow in a river downstream, the geometry of the flow retains deformation proportional to the frequency of blockage and the rate of flow repair. If we look back to Section 2.2, 2.3 we see that a series of unitary operations occur when path evaluation (curvature evaluation) occurs, in other words, a particle transits through a region. Any transiting particle, such as a photon emitted from a source, would cause interactions resulting in the updating of the states normalization, therefore entropy and the curvature of spacetime. Conversely, in very low-temperature settings, fewer and lower energy particles transit through the medium resulting in less impactful changes to curvature. This may give us insight into why Dark Matter halos are predominantly viewed at the edge of galaxies (see, e.g., Bullock and Boylan-Kolchin, 2017 [27]), as dense regions around star formation tend to have higher temperatures and higher average flux density. As such, We can view the above boundary condition as a measure of the proportion of the free energy and locked entropic energy within a system, in other words:

$$T\propto {\displaystyle \frac{d{E}_f}{d{E}_l}}$$

(78)

In other words, the higher the temperature, the more dominant free energy will be within the system compared to locked energy, i.e. entropically defined energy. A direct examination of this ratio will be conducted in a later paper. Until then the reader can ponder if the velocity of the leading edge of gravitational waves, and thus gravitational energy density can be directly affected by temperature as it acts on gravitational wave flux density.

4. Discussion

The relationship between quantized spacetime and entropy should be clear. Entropic work is performed by the observation or translation of any particle. For simple particles without complex spin symmetries, this time evolution is forced into gradient dispersion as a direct consequence of the projection of a particle its prior state location influencing its current state location. That entropy is further increased requiring more work to be performed when the neighboring state is dissimilar. Since curvature increases with increased energy, this provides a mechanism for the experience of smooth spacetime curvature as spacetime contracts or expands to become more homogeneous with its surrounding states. It is important to note here that a relaxing volume necessitates its neighbor volume to contract as a particle with a given energy is essentially annihilated and then created in the new volume because of the field propagation. If contraction and expansion of a volume occur non-linearly (such as in General Relativity), a macro-observer in a stable system would observe a change of relative volumes dramatically smaller than the lower bound of the non-zero minimum length. This would make observing this mechanism only possible under extreme circumstances where the gravitational field is violently disturbed. Also note that under this formulation gravitational fields are the combinatory field of all energies in all existing fields, I find no reason to assume that this entropic process prefers energy propagation in any particular mode.

This framework gives us a time evolution description for curvature and a relativistic energy density which relates directly to entropy and increasing system symmetries. This description translates directly to expected results for black holes and dark energy. In particular, we expect dark matter density within the halo to fall off proportional to the cube of radius at large distances.

Numerical simulations and observational studies consistently indicate that dark matter halos exhibit a characteristic two-regime density profile. In the inner regions of galaxies, dark matter is expected to follow a shallow cusp with a density scaling roughly as $\rho \left(r\right)\propto 1/r$, as predicted by the Navarro–Frenk–White (NFW) profile [28]. Observational evidence from rotation curve analyses, such as those of M33, supports this central cusp behavior [29], while weak lensing measurements further confirm that the inner mass distribution is well described by such a profile [30]. High-resolution N-body simulations—including results from the Aquarius Project—demonstrate that halos transition from this $1/r$ behavior in their cores to a steeper $\rho \left(r\right)\propto 1/r^3$ decline in the outer regions, a change that is essential for ensuring a finite total halo mass [31,32]. Moreover, empirical models constructed using nonparametric techniques have further validated that the outer halo density follows this steep decline [33]. This steep decline is exactly what we would expect from a Ricci Flow decay mechanism that is based on the ratio of free propagating energy to total energy i.e. $1-(E_f/E_t)=TS/E_t$. Additionally, we can view the expansion of our observable universe to be an expansion of our light cone whereby low entropy states are being added to our relativistic worldview. By using theory which shows that entanglement information is inversely related to coherence information, also seen in information dynamics (Vopson, 2022) you can derive a relationship between the gradient of entropy. Specifically, one can find an entropic relationship between entanglement and coherence that is in covariant and contravariant forms of a new metric tensor. Additionally, we can find that curvature can be described here as worldlines expressing harmonized modes of information where the real and complex components of the eigenvalues are. These scalar multiples relate directly to modes of force that correspond to the 2D plane delta curvature (as compared to Riemann Flat space). This geometric analysis reveals interesting results given the entropically defined time evolution boundary conditions.

4.1. Limitations

Several key limitations must be considered within this paper. First, in deriving the relationship between the change of entropy and energy we assumed stable, non-changing temperature (hence from product rule $SdT=0$. This cannot be assumed for the vast majority of applications, however, for black holes $dT/dE\propto -1/E^2$ and as such should not play a dominant role, particularly on a first-order approximation. Further work should be conducted in fully incorporating dT as simplifications are highly likely. Secondly, we constructed this formulation at extreme boundaries near a singularity, where mass is pulled towards the center. As such we assumed that the entropic component of total energy was dominant for the region in between the singularity and the event horizon. In systems with elevated levels of free energy, this component must be considered for its impact on curvature deformation according to the stress-energy tensor.

More significantly, is the potential for an alternative decomposition which may better "fit" the distribution divergences. While this author attempted to ground the decomposition in sound logic and fundamental principles (least action, geometric interpretation) one cannot ignore the fact that much of these ideas are derived from approximation and assumption from best fit.

5. Conclusions

In conclusion, our investigation challenges the traditional path-independent assumptions of tensor contraction in spacetime, revealing that significant information is lost in this simplification. The evidence points to the necessity of path-dependent analyses to account for the entropic dynamics inherent in spacetime, particularly at quantum scales. This realization opens new avenues for research, particularly in the examination of the Riemann Tensor with an emphasis on preserving information and embracing the entropic nature of spacetime evolution. Future work should explore these entropic path-dependent processes in greater detail, aiming to integrate them with existing theoretical frameworks to enhance our understanding of spacetime dynamics. This approach not only promises to refine our grasp of General Relativity also bridges the gap between classical and quantum physics by acknowledging the fundamental role of information entropy in shaping the universe.