1. Introduction
The vacuum of a quantum field theory is not empty, it is a seething sea of virtual particle-antiparticle pairs whose quantum states are strongly entangled over short distances [1]. In QFT, the vacuum two-point correlation functions and the entanglement entropy of spatial subregions obey an area law, implying that entanglement is predominantly local [2]. More distant regions share progressively less entanglement. Leonard Susskind has emphasized that “the story seems to be that entanglement stitches space together. It really is the stuff that holds space together” [3]. This idea traces back to the ER = EPR conjecture of Maldacena and Susskind, in which Einstein-Rosen bridges (wormholes) are identified with quantum entanglement between distant regions [3].
We take these insights one step further and propose that gravitational attraction is the macroscopic consequence of this entanglement network. Two massive bodies are gravitationally attracted precisely because they occupy regions of spacetime that, on average, share a greater density of entangled vacuum fluctuations than more distant pairs. The collective effect of countless such microscopic entanglements produces an emergent force that scales inversely with the square of the separation, precisely the Newtonian limit of gravity. In the language of graph theory, the nodes represent localized bundles of QFT degrees of freedom (or discretized spacetime points), and the edge weights are proportional to the mutual information between their vacuum states. Proximity in the embedding space therefore correlates with stronger entanglement edges.
Time itself emerges from the unitary evolution of the quantum state. Each infinitesimal time step corresponds to the application of a unitary operator on the entangled network.
Remark 1. The dual-cover property of unitaries in three real dimensions naturally selects a three-dimensional spatial geometry for the emergent manifold. This observation is consistent with the well-known double-cover relation
[4] [5]. We leave a rigorous treatment of this connection for future work.
The graph deforms continuously under these unitary updates, and the resulting curvature is identified with the Einstein tensor. In this way we obtain a microscopic, entanglement-driven mechanism for both the geometry of space and the dynamics of gravity, all within the framework of ordinary QFT.
Relation to Prior Work
This construction is related to, but distinct from, several recent approaches to emergent gravity. Jacobson [6] derived the Einstein equation as a thermodynamic equation of state, while Verlinde [7] proposed gravity as an entropic force arising from an information-theoretic holographic screen. Van Raamsdonk [8] demonstrated that entanglement entropy directly controls the connectedness of spacetime in AdS/CFT. Our approach differs from all three in working directly with the graph-theoretic Ollivier-Ricci curvature on an entanglement network, without assuming a holographic duality or a thermodynamic ensemble, and in deriving rather than postulating the gravitational coupling constant via an explicit Newtonian limit calculation.
2. Graph Theory and Entanglement Weights
Let the underlying structure be an undirected weighted graph
whose vertex set
consists of nodes located at coordinate positions
. Each node
carries a mass-energy
that will later be identified with the local energy density of the QFT vacuum.
The entanglement strength between neighboring nodes
and
is defined as the mutual information of their vacuum reduced density matrices [2] [9] [10]:
(1)
where
is the von Neumann entropy. In a relativistic QFT vacuum this quantity decays with distance according to the area law:
(2)
where
is the scaling dimension of the relevant field operator and
is a constant fixed by the QFT. The Euclidean distance between nodes is
(3)
3. Graph Structure and Geometry
We use metric signature
throughout. Lowercase Latin indices
run over spatial components 1, 2, 3; Greek indices
run over spacetime components 0, 1, 2, 3. We retain explicit factors of
where needed for the Newtonian limit.
Assumption 1 (Spanning). For each node
, the set of unit vectors
spans
.
We work in a locally Euclidean embedding. Define the entanglement-weighted metric at node
:
(4)
where
and
. This metric is symmetric and positive definite under the Spanning Assumption, and reduces to
in flat space with uniform entanglement.
Define a local curvature scalar between neighboring nodes
and
via the normalized Wasserstein-1 distance
between the mass measures supported on their respective neighborhoods [11]:
(5)
Positive, negative, and vanishing
correspond to positive, negative, and zero curvature, respectively:
(6)
Following Ollivier [11], the weighted transport measure from
is:
(7)
The discrete Ricci curvature at
is
(8)
The components of the discrete Ricci tensor are
(9)
The Ricci scalar follows by tracing with the induced metric:
(10)
The Einstein tensor on the graph is then
(11)
4. Continuum Limit of the Discrete Einstein Tensor
Consider a sequence of graphs
obtained as
-nets of a smooth Riemannian manifold
: a maximal set of points with pairwise distances
, connected whenever
. As
the vertex density grows as
for an
-dimensional manifold.
By the measured Gromov-Hausdorff convergence theorem [12] [13], optimal transport distances converge, and Ollivier’s expansion (Theorem 7 of [11], sharpened for irregular graphs by [15]) gives, for the unit vector
.
(12)
Substituting into
and replacing the discrete sum over neighbors by an integral over directions on
:
(13)
Applying the spherical averaging identity
(14)
and contracting with the symmetric tensor
yields, after appropriate normalization by
:
(15)
Consequently,
(16)
distributionally, recovering the continuum Einstein tensor. The convergence holds pointwise for smooth metrics and in
for metrics with bounded curvature [14].
5. Coupling to Stress Tensor
Identifying the node masses
with the local vacuum energy density, the smoothed density
(17)
(with
a normalized Gaussian kernel) supplies the
component. The full stress-energy tensor on the graph is the perfect-fluid form [16]:
(18)
where the four-velocity
is the mass-weighted average of node velocities,
, satisfying
in the continuum limit and the pressure
receives contributions from kinetic velocity fluctuations among neighbors and from the gradient of the entanglement weights
.
Writing out the components in the slow-motion limit:
(19)
(20)
and
(21)
The tensor is manifestly symmetric (
) and its trace is
, recovering
for pressureless dust and 0 for radiation (
) as expected.
In the continuum limit the graph Einstein tensor satisfies
(22)
with the coupling constant fixed by matching to the Newtonian limit (see Section 7). The right-hand side is sourced entirely by the same mass-energy distribution that defines the entanglement network, so that gravity is the geometric manifestation of vacuum entanglement rather than an independent fundamental force.
6. Discrete Contracted Bianchi Identity
Theorem 1 (Discrete Contracted Bianchi Identity). Let
be a weighted graph satisfying the Spanning Assumption, with Ollivier-Ricci curvature
and discrete Einstein tensor
. Then
for all
. Consequently, the field equation
implies energy-momentum conservation
.
Proof. Because the Wasserstein transport distance is symmetric,
, the edge curvatures satisfy
. Computing the discrete divergence of
and using
:
(23)
The discrete gradient of the Ricci scalar is
(24)
Hence
and therefore
(25)
The deeper reason this identity holds is that the graph possesses a relabeling symmetry: permuting vertex labels leaves all physical quantities invariant. This discrete symmetry plays the role of diffeomorphism invariance, and the Bianchi identity is its Noether consequence [17] [18]. For the analogous result in combinatorial Ricci curvature settings see [19].
Energy-momentum conservation
is therefore not an independent assumption but a theorem forced by the geometry of the entanglement network, recovering the discrete work-energy theorem and the discrete Euler equation.
7. Newtonian Limit and the Coupling Constant
We impose the three standard conditions [20] [21]: 1) weak field:
with
; 2) slow motion:
; and 3) static:
.
(26)
Under these conditions the dominant component of the Einstein tensor reduces to
(27)
where
identifies
as the Newtonian gravitational potential. The stress-energy component is
. The graph Laplacian
in the continuum limit [22]. Inserting into the field equation
:
(28)
Matching to Poisson’s equation
fixes
(29)
recovering the standard Einstein coupling constant. The graph Poisson equation is thus which in the continuum limit is the Newtonian gravitational field equation.
8. Conclusions
We have shown that a graph-theoretic model whose edge weights encode the entanglement structure of a QFT vacuum naturally produces the Einstein tensor and couples it to the stress-energy tensor. The derivation relies only on 1) the locality of vacuum entanglement, 2) the Wasserstein geometry of mass transport on neighborhoods, and 3) unitary time evolution. In the continuum limit the discrete equations become the Einstein field equations of general relativity, with the coupling constant derived from the Newtonian limit rather than inserted by hand, and energy-momentum conservation following as a theorem from the discrete Bianchi identity.
This framework supplies a purely quantum-information-theoretic origin for gravity within ordinary quantum field theory. It reproduces known results on emergent geometry from Ollivier-Ricci curvature on random graphs and opens new avenues for computing gravitational phenomena directly from entanglement entropies. Future work will include explicit lattice-QFT simulations to extract the effective Newton constant from vacuum mutual information and extensions to curved backgrounds and quantum cosmology.