Induced Gravity on the Selection–Stitch Lattice

Induced Gravity on the Selection–Stitch Lattice:
Covariance from Vacuum Fluctuations, Strength from
Entanglement
Raghu Kulkarni
SSMTheory Group, IDrive Inc., Calabasas, CA 91302, USA
raghu@idrive.com
Abstract
In the Selection–Stitch Model (SSM), space is a lattice of entangled bonds. The vacuum’s
quantum fluctuations, the ocean, induce the Regge action on this lattice if they are covariant.
We ask why the ocean is covariant and what fixes Newton’s constant, and answer both with two
forms of one principle. First, the ocean is maximally mixed. At fixed density, the maximally
mixed state of its positions is a Poisson process. If its couplings also keep no memory of an
earlier geometry, its vacuum energy depends on the metric only through the volume, and the
graviton stays massless. A sprinkled ocean responds to a shear at 1
.
001
±
0
.
002 of the covariant
value. The same operator on the lattice gives 1
.
23, and a memory
µ
of an old geometry gives
exactly 1
−µ
. The ocean has no density halo around a defect. Second, maximal entanglement. If
each bond’s entanglement follows its share of area, entanglement equilibrium is the Hamiltonian
constraint. A solved defect slice satisfies it to within 5% once a ball spans two lattice spacings. If
the bond density maximizes the entanglement per unit area, that maximum is 1
/
4
G
. A two-qubit
bond model realizes this law. Since the bonds are frozen ocean, the two results describe one
substance.
Keywords: induced gravity; covariance; maximally mixed state; entanglement equilibrium; Regge
calculus; Newton’s constant
1 Introduction
Sakharov proposed that gravity is induced by the vacuum’s quantum fluctuations [
1
,
2
]. In this
picture, the action of gravity is the energy cost of bending the space the fluctuations live on. The
Selection–Stitch Model (SSM) gives this idea a concrete setting. Its space is a lattice of entangled
bonds. Its curvature sits on the lattice’s hinges as Regge deficits [
3
]. Section 2 describes the model.
Ref. [
4
] showed three things. If the vacuum’s fluctuations are covariant, they induce exactly
the Regge action on the lattice. The fluctuations cannot be the lattice’s own fields, because those
give the graviton a mass. A defect’s verification cost acts as a rest mass. As there, we call the
fluctuations the ocean. Two questions were left open. One is why the fluctuations are covariant.
The other is what fixes the strength of gravity, Newton’s constant G.
This paper answers both questions with two forms of one principle. The principle is maximal
entropy in quantum form. The ocean’s state has the largest von Neumann entropy its density allows,
and the bonds carry the largest entanglement entropy per unit area.
1.
The form of gravity comes from a maximally mixed ocean. The ocean is in the maximally
mixed state at fixed density, and it keeps no memory of earlier geometries. Then its vacuum
1
energy depends only on the volume, which is covariance (Sections 3–6). The same ocean has
no density halo around a defect (Section 7).
2.
The strength of gravity comes from maximal entanglement. The lattice’s bonds carry the vac-
uum’s entanglement. Their equilibrium with matter is the Hamiltonian constraint (Section 8).
Maximizing their entanglement per unit area fixes 1/4G (Section 9).
In the SSM the bonds are made from the ocean, as ice is made from water. The two answers
therefore describe one substance (Section 10). Section 11 lists what remains open and grades each
claim (Table 5).
2 The SSM framework
This section summarizes the parts of the SSM that this paper uses. Each part rests on a published
or archived paper. This paper does not test those papers’ claims.
The vacuum assembles into a lattice. The basic object of the SSM is a bond, a Bell pair that
joins two nodes at a fixed distance [
5
]. The network grows by two operations (Fig. 1a). A stitch
adds a node in the plane of an existing bond. A lift adds a node above a triangle. Stitches build flat
sheets with six neighbors per node. Rare lifts build a foam of regular tetrahedra. That foam cannot
fill space. Five tetrahedra around an edge leave a gap of 7
.
4
◦
. To relieve this strain, the foam turns
into stacked sheets. The result is the face-centered-cubic (FCC) lattice, with twelve neighbors per
node (Fig. 1b,c). Twelve is the most that can touch one sphere in three dimensions, and FCC is
the densest packing of equal spheres [
6
]. A random simulation of the two operations supports this
picture [5].
The lattice is a quantum code. Ref. [
7
] puts one qubit on each edge of the FCC lattice. It
defines two kinds of check (Fig. 2a,b). A
Z
-check acts on the twelve bonds of a node. An
X
-check
acts on the twelve edges of an octahedral void. The two kinds commute, so they define a CSS
code [
8
]. On the smallest periodic box the code is [[192
,
130
,
3]], and its rate tends to 2
/
3. In the
SSM the vacuum is the code space of this code. A state that satisfies every check is a vacuum state.
In the SSM’s reading, a bond persists because the code keeps verifying it.
Matter is a trapped node. If crystallization is incomplete, a node of the foam can be left at
the center of a tetrahedral void [
5
]. It is bonded to the four corners of the void by four short bonds,
which we call spokes (Fig. 3c shows a two-dimensional cartoon). The node cannot leave the void.
To pass through a face, it would have to come closer to the three corners than the lattice allows.
Ref. [
5
] reads this as confinement and identifies the defect with a baryon. It also derives charges
and colors from the defect’s geometry. This paper uses only the geometry of the defect.
Mass is the cost of verification. Ref. [
9
] proposes that a defect’s mass is the cost of checking it.
The cost is a count, the number of qubits the defect touches times the number of checks that can
detect it. Ref. [
10
] states the step from a count to an energy as a postulate for standing information.
Information that is held in place and never erased has energy. Each bit carries a fixed amount
ε
b
, so
a defect’s energy is its count times
ε
b
. Landauer’s principle prices erasure and does not apply here.
Only ratios are used, so mass ratios are ratios of integers. The five costs that survive its selection
rules match the electron, muon, pion, proton, and neutron masses to within 0
.
12%. This paper does
2
stitch: new node
2 bonds, in plane
lift: new node, 3 bonds
out of plane,
h
=
√
2/3
L
0
(a) the two assembly operators
(b)
K
= 12 = 6 + 3 + 3
(FCC, ABC stacking)
6 in plane (stitch)
3 above (lift)
3 below (lift)
Bell pair
K
= 1
stitched sheet
K
= 6
tetrahedral foam
K
= 4
, frustrated
FCC crystal
K
= 12
, Kepler bound
stitch
lift,
P
lift
=
e
−
3
deficit
δ
= 0.128
rad relieved
a node left in a tetrahedral void = matter
(c) the assembly cascade
Figure 1: Assembly in the SSM [
5
]. (a) The two operations. A stitch adds a node in the plane, and
a lift adds a node above a triangle. (b) The final coordination, with six neighbors in the node’s own
sheet and three in each sheet above and below. (c) The cascade from one bond to the FCC crystal.
A node left behind in a tetrahedral void is the SSM’s matter.
(a)
Z
-check on a node (b)
X
-check on an octahedral void
(c) the 13-site cluster
Figure 2: The SSM code [
7
,
9
]. (a) A
Z
-check on the twelve bonds of one node. (b) An
X
-check on
the twelve edges of one octahedral void. (c) The 13-site cluster of a node and its neighbors, whose
counts set the verification costs of Ref. [9].
not use these values. It uses only one result of Ref. [
4
]. A verification cost couples to geometry as a
rest mass, if verification ticks measure proper time.
Gravity is Regge calculus on the lattice. In spacetime the SSM lattice is the
D
4
lattice,
whose slice at fixed time is FCC [
11
]. Its geometry is fixed by its edge lengths alone. Curvature sits
on the hinges as deficit angles, and the Regge action is
P
h
A
h
δ
h
/
8
πG
over the hinges [
3
]. On
D
4
this action reproduces the Einstein–Hilbert action exactly at second and third order. Its static field
gives Newton’s law and the post-Newtonian parameters
β
=
γ
= 1
.
000
±
0
.
002 [
12
]. Ref. [
11
] sets
the lattice scale. It matches one bit of entanglement per bond to the Bekenstein–Hawking entropy,
which puts the bond length at the Planck scale. We use the calibration of Ref. [
4
],
G
=
L
2
f
/
4
ln
2,
with
L
f
the
D
4
nearest-neighbor distance, so
L
f
= 1
.
665
ℓ
P
. Ref. [
11
] also treats black holes, as
regions in which the lattice has broken down.
3
(a) the ocean (b) frozen ocean: the lattice (c) a trapped node: the bonds bend
Figure 3: The SSM’s two phases, drawn in two dimensions. (a) The ocean, made of degrees of
freedom with no orientation. (b) The frozen phase, a lattice of bonds with the ocean around it. (c)
A trapped node (red) in a void, bonded by spokes. Nearby bonds change length, and this is the
curvature. The ocean is excluded only from the node’s own volume (Section 7).
The ocean induces the Regge action. The vacuum’s quantum fluctuations are measured
directly today [
13
]. In the SSM they are older than the lattice. The bonds and defects were
stitched from them, and they persist around the finished lattice [
4
]. They are the ocean of Section 1
(Fig. 3). Ref. [
4
] shows that a covariant ocean induces the Einstein–Hilbert action. On the lattice’s
piecewise-flat geometry this action is exactly the Regge action [
14
]. The check is done on the conical
hinges themselves, from the exact spectrum of a cone. Fields that live on the lattice itself fail.
Their vacuum energy depends on the shape of the lattice cells, and that gives the graviton a mass.
Covariance is the postulate this paper derives.
3 The covariance test
Covariance has a simple test on a lattice. Strain the lattice uniformly. The result is still flat space,
written in other coordinates, and every Regge deficit stays zero. A covariant vacuum energy may
therefore change with the volume, but not with the shape.
We place the ocean’s degrees of freedom at points
X
i
in a periodic four-dimensional box. The
points move with the lattice. We give the box a uniform metric
g
. The ocean’s fluctuation operator
is a graph Laplacian L(g), and its one-loop vacuum term per degree of freedom is
f(g) =
1
N
ln det
′
L(g), (1)
where
det
′
omits the zero mode. The couplings are Gaussian in proper distance, with a smooth cut:
w
ij
= e
−d
2
ij
/2σ
2
1 − d
2
ij
/r
2
c
3
(d
ij
< r
c
), d
2
ij
= ∆x
T
ij
g ∆x
ij
. (2)
Lengths are in half-cells of the D
4
lattice, and there is one point per D
4
lattice site on average.
Let
κ
1
=
1
4
df/dϵ
be the volume response, with
g
= (1 +
ϵ
)
I
. Let
K
s
=
d
2
f/dt
2
be the shear
response, with
g
=
I
+
tH
, where
H
is traceless and
tr H
2
= 1. If
f
depends only on
det g
, then
K
s
= −κ
1
. We report
R =
K
s
−κ
1
, (3)
4
which is 1 for a covariant medium. If
R
= 1, the vacuum energy depends on the shape of flat space.
In the continuum this is a mass term for the graviton [
15
,
16
]. Gravitational waves bound the
graviton mass tightly [17], so R = 1 must hold almost exactly [4].
4 Placement from a maximally mixed ocean
Suppose we know only the mean density n of the ocean’s degrees of freedom. The placement with
the most entropy is then the Poisson process of intensity
n
, in which the points are independent
and uniform [
18
,
19
]. Any correlation, preferred spacing, or preferred direction would be extra
information. It would lower the entropy. In quantum terms, this placement is the maximally mixed
state of the ocean’s positions at fixed density, the state of largest von Neumann entropy with that
constraint. Its position statistics are those of a Poisson process. The state need not be prepared.
Almost every pure state of a large system looks maximally mixed on a small part of it [20].
A Poisson process has the symmetry that covariance needs. Let
A
be a linear map with
det A
= 1. Moving every point by
A
turns a Poisson process into another Poisson process with
the same intensity [
19
]. The couplings depend only on proper distances. Moving the points by
A
and changing the metric to
A
−T
gA
−1
therefore leaves the operator unchanged. Averaged over the
process,
¯
f(A
T
gA) =
¯
f(g) for every A with det A = 1. (4)
Any metric can be brought to a multiple of the identity by such a map. As a result,
¯
f
depends
on
g
only through
det g
, and
R
= 1. This is exact for the average. A single sprinkling fluctuates
around it, and the fluctuation averages out over many points. The density is the only scale the
placement has, so it sets the ocean’s cutoff. In spacetime, the same symmetry is the theorem that a
Poisson sprinkling picks out no frame [
21
]. Causal sets rest on this property [
22
]. Random lattices
were introduced in flat space for a related reason, to restore rotation and translation symmetry on
average [23]. Ref. [24] built induced gravity from matter on a lattice of randomly placed points.
A regular lattice lacks the symmetry. A volume-preserving map moves its points onto a different
lattice. Table 1 compares the two with the same operator. The sprinkled points pass the test. The
lattice points fail it, and they fail more as the kernel narrows and the operator sees the lattice’s
arrangement more closely.
Table 1: Shear response relative to the covariant value, Eq.
(3)
, for the same memoryless operator
on sprinkled points and on
D
4
lattice points. Box side 8 (2048 points). The sprinkled value is an
average over 10 sprinklings with random shear directions. The two lattice values are for a diagonal
and an off-diagonal shear.
points kernel width σ cut r
c
R
Poisson sprinkling 0.7 2.24 1.0009 ± 0.0015
D
4
lattice 0.7 2.24 1.233, 1.234
D
4
lattice 1.0 3.2 0.994, 0.995
5 Couplings without memory
The argument above assumed that the couplings depend only on the current metric. A medium with
no structure of its own has nothing else to depend on, except a second geometry, a reference metric
5
g
0
. A medium that formed in one geometry, and kept its couplings when the geometry changed,
carries such a g
0
. We call this memory.
Memory allows shape terms. With a reference metric, Eq.
(4)
holds only when
g
and
g
0
are
moved together. The vacuum term can then depend on the eigenvalues of
g
−1
0
g
. Terms such as
tr
(
g
−1
0
g
) change under a shear at fixed volume, and nothing forbids them. These are the potential
terms of massive gravity, which needs a reference metric for this reason [
16
,
25
]. A medium with
memory could still be covariant, but only if every shape term happened to cancel. Without memory,
covariance is forced by symmetry.
An exact example. Let every coupling see the blended metric
g
′
= (1
− µ
)
g
+
µg
0
, with
g
0
=
I
and memory
µ
. If the memoryless medium is covariant, then
f
(
g
) =
ϕ
(
det g
′
) for some function
ϕ
.
A compression changes
det g
′
at first order with weight 1
− µ
. A shear changes it at second order
with weight (1 − µ)
2
. The ratio is then
R(µ) = 1 − µ. (5)
Only
µ
= 0 and
µ
= 1 are covariant. At
µ
= 1 the medium ignores the geometry completely. Then
K
s
and
κ
1
both vanish, and the medium induces no gravity. Any value in between gives the graviton
a mass. The shape excess is
K
s
+
κ
1
=
µ
(1
− µ
)
κ
1,0
, where
κ
1,0
<
0 is the memoryless volume
response.
Memory of neighbors. A milder memory keeps the coupling strengths current but remembers
which points are neighbors. In Eq.
(2)
, let the cut factor use the old distance
d
0
for a random
fraction
ϕ
of the pairs. The other pairs use the current distance. The Gaussian factor always uses
the current distance. Table 2 and Fig. 4a show both kinds of memory on sprinkled points. The
blend follows Eq.
(5)
within sampling error. Remembering neighbors also breaks covariance, and
the break grows with the fraction remembered. When every pair remembers,
R
= 0
.
23. A fixed
triangulation is another memory of neighbors. We built the Delaunay triangulation of a sprinkling,
held it fixed, and computed the finite-element Laplacian from the current edge lengths. It gives
R
= 1
.
96 (box side 5, 625 points at twice that density, two sprinklings). Memory can push
R
either
way. Its sign and size depend on how the old structure enters the couplings.
Table 2: Shear response relative to the covariant value for two kinds of memory on a Poisson
sprinkling (box side 6, 648 points, 20 sprinklings). Blend:
σ
= 1,
r
c
= 3
.
2. Remembered neighbors:
σ = 1, r
c
= 2.
memory µ R, blend fraction ϕ R, remembered neighbors
0 0.9995 ± 0.0009 0 0.999 ± 0.003
0.25 0.7497 ± 0.0007 0.25 0.942 ± 0.003
0.5 0.4998 ± 0.0005 0.5 0.845 ± 0.003
0.75 0.2499 ± 0.0002 0.75 0.675 ± 0.002
0.9 0.1000 ± 0.0001 1 0.232 ± 0.001
6
0.0 0.2 0.4 0.6 0.8 1.0
memory
μ
or
ϕ
0.0
0.2
0.4
0.6
0.8
1.0
shear response / covariant value
(a) memory breaks covariance
1
−
μ
blend: metric
(1
−
μ
)
g
+
μ
g
0
remembered neighbors, fraction
ϕ
0.0 0.5 1.0 1.5 2.0 2.5
distance from the defect (half-cells)
0.0
0.2
0.4
0.6
0.8
1.0
1.2
density / far density
defect
(b) no halo around a defect
maximally mixed ocean
a
1/
r
halo, for comparison
Figure 4: (a) Shear response relative to the covariant value against memory, for the blend of Eq.
(5)
(circles, with the line 1
−µ
) and remembered neighbors (squares). (b) Density of a maximally mixed
ocean around an excluded ball of radius 0
.
6, relative to the far density, averaged over 400 sprinklings.
A 1/r profile is shown for comparison.
6 Memory is verification
In the SSM the lattice persists because the code verifies it (Section 2). The ocean is the part of
the vacuum that no check verifies. Its fluctuations arise and decay without being recorded. We
identify these two kinds of persistence with the two kinds of memory. A verified structure keeps its
connections and its formation geometry, so fields that live on it carry a reference geometry. The
unverified ocean keeps no record, so its couplings can depend only on the geometry as it is now,
which is µ = 0.
This identification explains the two main results of Ref. [
4
] as one result. Fields on the lattice
live on remembered structure. The symmetry of Section 4 does not protect them, and every coupling
rule tested in Ref. [
4
] gave a shape-dependent vacuum energy. The ocean carries no reference
geometry. Eq.
(4)
holds for it, and its induced action is covariant, so by Ref. [
4
] it is the Regge
action. The identification is a modeling step. It fits a model that defines persistence by verification,
but this paper does not derive it from the code’s dynamics.
7 The ocean around a defect
A defect occupies a small volume
v
that the ocean cannot enter. Take the maximally mixed
placement again, with that volume excluded. The points are absent inside it and uniform outside
it. Their density outside is
N/
(
V − v
). It is higher than without the defect by the fraction
v/V
,
and the same at every distance. There is no gradient. A Poisson process has no interactions, so
no point can sense that the defect is near. Fig. 4b shows the profile. It is flat within the sampling
noise, right up to the edge of the excluded ball.
A density halo would need an energy cost for gradients of the density,
κ
2
|∇n|
2
. A maximally
mixed medium has no such cost. Even with one, and with a compressibility term
1
2
f
′′
δn
2
, the
density would relax as
e
−r/ξ
/r
with
ξ
=
p
κ/f
′′
. That is short-ranged. A 1
/r
tail would need
f
′′
= 0. Two consequences follow. First, the 1
/r
field of a mass is not a thinning of the ocean. It
is the curvature of the geometry, which the induced Regge action governs (Fig. 3c). Second, the
7
ocean’s density, and so its cutoff, is the same near matter as far away. A Newton constant set by
that cutoff therefore does not vary near masses, as tests of the equivalence principle require.
8 The strength of gravity from entanglement equilibrium
Covariance fixes the form of the induced action. It does not fix the coefficient, 1
/G
. For the
coefficient we use the entanglement that the lattice carries. Each bond is a Bell pair, and the
vacuum’s entanglement across a surface is carried by the bonds the surface cuts. That Einstein’s
equation is a condition on the vacuum’s entanglement is Jacobson’s idea. He first derived it as an
equation of state of local horizons [
26
], and later as the equilibrium of the entanglement in small
balls [
27
]. In holography the linearized Einstein equation follows in the same way from the first law
of entanglement [
28
,
29
]. Others have sought geometry directly from the entanglement of abstract
degrees of freedom [
30
,
31
]. Here the degrees of freedom are the bonds. Two postulates define the
law.
Postulate A: entanglement follows area. The entanglement across any surface Σ in a slice is
S
(Σ) =
A
(Σ)
/
4
G
, where
A
is the area in the lattice’s geometry. In the flat vacuum, each bond that
Σ cuts carries entropy s
0
and a share L
2
f
of its area. Then
G =
L
2
f
4s
0
. (6)
With
s
0
=
ln
2 this is the calibration of Section 2. When the lattice strains, a bond’s share of area
changes. Its entanglement must change with it, as
s
e
=
s
0
a
e
/a
(0)
e
. Here
a
e
is the area of the bond’s
dual face.
Postulate B: the vacuum is in entanglement equilibrium. For every small ball at fixed volume,
the change of entanglement from the geometry balances the change from the matter inside:
δA|
V
4G
+ 2π δ⟨K⟩ = 0, K =
Z
B
R
2
B
− r
2
2R
B
T
00
dV. (7)
Here
K
is the modular Hamiltonian of a ball of radius
R
B
in flat space [
27
]. On a time-symmetric
slice, to first order in the curvature,
δA|
V
=
−
R
B
(
R
2
B
− r
2
)
/
(4
R
B
)
3
R dV
. Equilibrium for every
ball then requires
3
R
= 16
πGρ
at every point. This is the Hamiltonian constraint of general
relativity. Two checks fix the normalization. For a point mass, Schwarzschild initial data give
δA|
V
=
−
4
πGMR
B
, and 2
πδ⟨K⟩
=
πMR
B
. For a straight cosmic string of mass
λ
per unit length,
the two sides balance when its deficit angle is 8πGλ.
Three properties follow.
G
is fixed by the area law, with no sum over particle species. A uniform
vacuum energy changes no entanglement, so it does not enter, and the cosmological constant is only
a constant of integration [
32
,
27
]. With boosts, equilibrium of every small causal diamond gives the
full Einstein equation [
26
,
27
]. We use only the time-symmetric slice. The source is the defect’s rest
mass, which by Ref. [4] is its verification cost.
A test on the lattice. On a piecewise-flat slice,
3
R dV
becomes 2
P
e
ℓ
e
δ
e
on the edges. Equilib-
rium then becomes one equation per node:
X
e∋v
ℓ
e
δ
e
= 16πG m
v
, (8)
the vertex form of the Hamiltonian constraint in 3+1 Regge calculus [
33
]. Here
m
v
is the rest mass
assigned to node
v
. We place the trapped node of Section 2 in the FCC lattice (Fig. 5a). Its rest
8
A
B
C
D
c
(a) trapped node: spokes (red), cage (blue)
0.6 1 2 3 4 6
distance of edge midpoint from void center (half-cells)
10
−
2
10
−
1
largest
|
δ
e
| /
δ
spoke
r
−
3
(b) curvature off the spokes (
q
= 0.1
)
Figure 5: (a) The trapped node
c
at the center of a tetrahedral void, joined by four spokes (red)
to the corners
A
–
D
. The void’s edges are blue. (b) Curvature off the spokes in the solved slice at
q
= 0
.
1 and box side 16. Points show the largest
|δ
e
|
among edges in each distance shell, in units of
the spoke deficit. The dotted line shows r
−3
. The wall at r = 7.6 steepens the last points.
mass is spread over the four spokes, with mass
λ
per unit length, and each spoke gives half its mass
to each end. We write
q
= 8
πGλ
. The slice is conformally flat,
ℓ
e
=
ψ
u
ψ
v
ℓ
(0)
e
. This is an assumption
about the initial data. It is exact for Schwarzschild initial data and at first order. We impose Eq.
(8)
at every node inside a wall of radius
R
w
. A wall is needed, because a closed three-torus admits no
metric whose scalar curvature is nonnegative and somewhere positive [34, 35].
We solve the system by least squares on a periodic FCC box. At the trapped node only the
spokes meet. By the symmetry of the void, each spoke carries the deficit
δ
s
=
q
= 8
πGλ
exactly.
This is the deficit of a cosmic string. Away from the defect, the conformal factor is
ψ
= 1 +
GM/
2
r
.
This is Schwarzschild initial data, with
M
equal to the source mass to 0
.
36% at box side 16. Fig. 5b
shows how the curvature falls off with distance.
We then test Eq.
(7)
directly on the solved slice, at box side 16, wall radius 7
.
6, and
q
= 0
.
01.
Each tetrahedron carries the constant metric fixed by its six edge lengths. Through these metrics
we compute the area of a sphere of radius
R
B
and the volume inside it. We compare the area deficit
at fixed volume,
δA|
V
, with the matter side,
−
8
πG δ⟨K⟩
, from the spoke source. This comparison
checks consistency and does not test the postulates independently. The slice was built to satisfy
Eq.
(8)
, and at first order equilibrium of every round ball is equivalent to it. The check tests whether
this equivalence survives on a lattice whose curvature sits on edges. Fig. 6 shows that it does.
The ratio of the two sides is 0
.
84, 0
.
96, 0
.
98, 0
.
99, and 0
.
99 at
R
B
= 1
.
5, 2
.
5, 3
.
5, 4
.
5, and 5
.
5. It
approaches one once the ball spans about two lattice spacings. Balls with no source inside are in
equilibrium to about 0
.
5%. Regions built from whole tetrahedra fail even without a source, because
their jagged boundaries pick up single edge strains. Equilibrium is a statement about round balls.
9
2 3 4 5
ball radius
R
B
(half-cells)
0.02
0.04
0.06
0.08
magnitude
(a) the two sides of the equilibrium condition
matter side,
8
π
G
δ
⟨
K
⟩
geometry side,
−
δ
A
|
V
2 3 4 5
ball radius
R
B
(half-cells)
0.80
0.85
0.90
0.95
1.00
ratio
two lattice
spacings
(b) geometry side / matter side
Figure 6: Entanglement equilibrium, Eq.
(7)
, on the solved trapped-node slice (box side 16, wall
radius 7
.
6,
q
= 0
.
01, lengths in half-cells, lattice spacing
√
2
). (a) The geometry side and the matter
side. (b) Their ratio.
9 The bonds’ microscopic law from maximal entanglement
Postulate A asks each bond’s entanglement to follow its share of area. This section asks what can
supply that rule.
Counting bonds is not enough. One could take the entanglement across a surface to be
s
0
times the number of bonds it cuts. But when the edge lengths change, that number stays the
same while the area changes. A count measures geometry only if curvature comes from adding or
removing bonds at fixed length. In the solved slices, curvature comes from changed edge lengths.
The code alone cannot supply the rule. In any stabilizer state, the entanglement across a cut
is a whole number of ebits. The checks that straddle the cut fix this number [
36
]. No check involves
the edge lengths. The code, as defined, therefore supplies at most the counting rule. The area rule
needs a bond whose entanglement changes with its length. For a pair of qubits this is the state
cos θ
e
|00⟩ + sin θ
e
|11⟩, with entropy
s(θ
e
) = −cos
2
θ
e
ln cos
2
θ
e
− sin
2
θ
e
ln sin
2
θ
e
, (9)
which rises to ln 2 as θ
e
→ π/4.
A minimal bond model. With
G
fixed by Eq.
(6)
, the area rule requires
d ln s
e
/d ln a
e
= 1 at
the vacuum. Give each bond’s two qubits the Hamiltonian
H
e
= −J(a
e
) (X
1
X
2
+ Z
1
Z
2
) − h (Z
1
+ Z
2
), J(a
e
) = J
0
a
e
/a
(0)
e
p
, (10)
so that the coupling grows with the bond’s dual area. The ground state has the form of Eq.
(9)
.
Table 3 shows its entanglement and response for
p
= 1. Stretching a bond raises its entanglement,
as required. But the slope is one only at h/J
0
= 0.69.
10
Table 3: Bond model of Eq.
(10)
with
p
= 1. Vacuum entanglement and its response to the dual
area.
h/J
0
0.10 0.30 0.50 0.69 1.00 1.50 2.00
s
0
/ ln 2 0.97 0.80 0.60 0.45 0.30 0.17 0.11
d ln s
e
/d ln a
e
0.06 0.39 0.75 1.00 1.25 1.45 1.54
Maximal entanglement fixes the slope. A second use of maximal entanglement selects this
point. It varies something different from Eq.
(7)
. Jacobson’s condition varies the geometry of a
ball, holds the vacuum’s microscopic state fixed, and takes 1
/
4
G
as given [
27
]. Here we hold the
geometry fixed and vary the microscopic state. The ball and its boundary area
A
stay the same,
but the density of bonds crossing the boundary changes. The boundary is crossed by
A/a
bonds,
where
a
is the vacuum’s dual area per bond. The entanglement across it is then
S
= (
A/a
)
s
(
a
).
We postulate that the vacuum’s bond density maximizes this entanglement. Then
s/a
is stationary
at the vacuum:
d
da
s(a)
a
a
0
= 0 ⇐⇒
d ln s
d ln a
a
0
= 1 ⇐⇒
∂S
e
∂a
e
0
=
s
0
a
0
=
1
4G
. (11)
This holds for any bond model in which
s/a
has a maximum. It makes 1
/
4
G
the largest entanglement
per unit area that the bonds can carry (Fig. 7a). Ref. [
37
] conjectured that the vacuum’s entanglement
per unit area equals 1
/
4
G
. Here the equality follows from the maximization of Eq.
(11)
, which we
postulate. In the bond model it fixes the parameters. For
p
= 1 the maximum is at
h/J
0
= 0
.
691,
with
s
0
= 0
.
45
ln
2 (Table 4). A coupling that follows the bond length instead,
p
= 1
/
2, has no
maximum.
This variation is allowed because the number of bonds is not conserved in the SSM. Bonds form
by selection and stitching [
5
], so a vacuum with a different bond density is an allowed state of the
same geometry. The variation compares uniform vacua. Changing the density in one region requires
restitching. This postulate is an addition to Jacobson’s condition, which does not imply it.
Table 4: Bond model at the maximum of
s/a
, Eq.
(11)
, for three exponents
p
. The slope is one
in each case.
G
loc
/G
is the local Newton constant at a 10% and a 49% stretch of the dual area,
under assumptions (i)–(iii) of the text. The last column is the stretch at which
G
loc
/G
= 1
.
08. For
p = 1/2, s/a has no maximum.
p h/J
0
s
0
/ ln 2 G
loc
/G (10%) G
loc
/G (49%) stretch at 1.08
1 0.691 0.45 1.08 1.54 9.8%
1.5 0.450 0.65 1.20 2.62 4.4%
2 0.358 0.74 1.33 4.33 2.8%
Calibration. With
G
fixed, Eq.
(6)
ties the bond length to
s
0
, as
L
f
∝ s
1/2
0
. The bond model’s
s
0
= 0
.
45
ln
2 gives
L
f
= 1
.
121
ℓ
P
instead of 1
.
665
ℓ
P
. Every result stated in lattice units is
unchanged. This includes the deficits, the equilibrium ratios, and the shear tests.
Strong fields. At a maximum,
s/a
curves downward. In any such model, a bond therefore follows
its area exactly at first order and falls below it at second order. Reading this shortfall as a change
of
G
takes three assumptions. (i) The bonds hold the stretch, with no restitching. (ii) The stretch
11
0 1 2 3 4
dual area
a
/
a
0
0.4
0.5
0.6
0.7
0.8
0.9
1.0
entanglement per area,
(
s
/
a
)/(
s
0
/
a
0
)
(a) the vacuum sits at the maximum of
s
/
a
0 10 20 30 40 50
stretch of the dual area
(
%
)
1.0
1.2
1.4
1.6
1.8
2.0
G
loc
/
G
neutron-star bound
(b) local Newton constant if bonds hold the stretch
p
= 1
p
= 1.5
p
= 2
Figure 7: The bond model. (a) Entanglement per unit area against the dual area, for
p
= 1 at
h/J
0
= 0
.
691. The vacuum sits at the maximum, where
d ln s/d ln a
= 1. (b) The local Newton
constant if bonds hold a stretch without restitching, for three exponents
p
. The dashed line is the
bound from two-solar-mass neutron stars [38].
is uniform over a ball small enough for Eq.
(7)
to apply, so that
G
loc
=
G s
′
(
a
0
)
/s
′
(
a
). (iii) The
stretch is measured from the vacuum edge lengths, which we estimate by
ψ
4
on a conformally flat
slice. Under these assumptions, for
p
= 1,
G
loc
grows by 8% at a stretch of 10%. At a neutron star’s
surface (1
.
4 solar masses within 12 km) the stretch is about 49%, and
G
loc
grows by 54% (Fig. 7b).
Two-solar-mass neutron stars bound
G
inside them to at most 8% above its laboratory value at
95% confidence, for a uniform change of
G
[
38
]. Applied to
G
loc
, this bound means that bonds do
not hold such stretches. The lattice must restitch once the stretch reaches a few percent to about
10%, depending on the model. How fully restitching relieves the stretch is not computed here. It is
where this law could depart from general relativity near compact objects.
10 One substance in two phases
The two halves of this paper describe one medium. In the SSM the bonds and defects are made
from the ocean, as ice is made from water [
5
,
4
]. Stitching fixes the ocean’s correlations into Bell
pairs on the bonds. A trapped node is a frozen remnant of the same process (Fig. 3). The liquid
phase is covariant and has no orientation, and it gives gravity its form. The frozen phase carries the
entanglement whose equilibrium is the field equation, and it gives gravity its strength. The ocean’s
entanglement across a surface and the bonds’ entanglement are therefore one quantity, counted once.
The SSM supports this reading, but no calculation here tests it.
This reading turns the induced Newton constant of Ref. [
4
] into a consistency condition. The
fields that froze into the bonds must together carry the bonds’ area law. Their induced coefficient,
1
/G
ind
=
C
Λ
2
/
12
π
, must equal 4
s
0
/L
2
f
. Here
C
is a weighted count of the field species. With one
ocean mode per
D
4
site, the cutoff is Λ
L
f
= 4
.
216 (heat-kernel count) or 5
.
013 (Debye count) [
4
].
The condition reads
C =
48πs
0
(ΛL
f
)
2
=
(
5.88 or 4.16 (s
0
= ln 2),
2.66 or 1.88 (s
0
= 0.45 ln 2).
(12)
In the proper-time scheme, the SSM’s field content gives
C
= 4 [
4
]. This is within a factor of 1
.
5 of
either requirement at
s
0
=
ln
2, and within a factor of 2
.
1 at
s
0
= 0
.
45
ln
2. But the weight of each
12
species depends on the regulator, not only on covariance. We also computed the weights with an
explicit covariant regulator, a sprinkled ocean with couplings in geodesic distance and path-ordered
spin transport. Scalars and gauge fields keep their proper-time signs, enlarged by factors of 2 to 4.
Each Weyl fermion gets a negative weight, from
−
0
.
6 to
−
12 in units of a minimal scalar, depending
on the fermion operator. With 48 Weyl fermions,
C
is then negative for this field content. Whether
Eq.
(12)
holds therefore depends on how the ocean’s fermions couple to curvature at the cutoff. In
this paper
G
is fixed by the bonds, so the answer changes no result above. It decides only whether
the ocean’s own fields reproduce that value.
11 What remains open
The bonds’ law. Two inputs of the bond model are assumed, not derived from the code. One is
how the coupling depends on the dual area,
J ∝ a
p
with
p
= 1. The other is the maximization of
Eq. (11). How fully restitching relieves strong-field stretches is not computed.
Memory and verification. That verified structure carries a reference geometry, and the unverified
ocean does not, is a modeling step (Section 6). Deriving it needs an account of how verification
records geometry.
The ocean’s own coefficient. Whether the ocean’s fields reproduce the bonds’ area law, Eq.
(12)
,
depends on how fermions are regulated at the cutoff. Covariance does not fix this.
Lorentzian signature and higher orders. The tests are Euclidean, and equilibrium is applied
on a time-symmetric slice. In Lorentzian signature a Poisson sprinkling is a causal set, which
picks out no frame [
21
]. But couplings that depend only on proper intervals are nonlocal there.
Entanglement equilibrium has no action, so it gives no corrections beyond Einstein’s equation. The
cosmological constant remains a constant of integration whose value is not fixed.
13
Table 5: Status of the claims.
claim status
SSM framework (Section 2) from Refs. [5, 7, 9, 10, 11, 12, 4]
maximally mixed placement at fixed density is a Poisson process standard [18, 19]
Poisson placement with memoryless couplings gives
R
= 1, Eq.
(4) exact (ensemble average)
sprinkled ocean passes, D
4
lattice points fail (Table 1) computed
a reference geometry allows shape terms known (massive gravity) [16, 25]
blend memory gives R = 1 −µ; remembered neighbors also fail derived and computed
verified structure remembers, the unverified ocean does not identification (SSM)
no density halo around a defect derived and computed
bonds’ entanglement follows area; the vacuum is in equilibrium postulates A and B
equilibrium = Hamiltonian constraint;
G
=
L
2
f
/
4
s
0
; Λ a constant of
integration
derived [27]
spoke deficit 8πGλ exactly; Schwarzschild far field computed
solved slice satisfies equilibrium ball by ball (Fig. 6) computed (consistency check)
stabilizer code gives whole ebits only known [36]
bond density maximizes entanglement, Eq. (11) postulate
bond model: slope one at the maximum; s
0
= 0.45 ln 2 for p = 1 computed
restitching above a model-dependent stretch estimate
bonds are frozen ocean; one quantity counted once interpretation (SSM)
the ocean’s own coefficient reproduces the area law open
12 Conclusion
Induced gravity on the SSM lattice has a form and a strength. Both follow from two forms of
one principle. The ocean is maximally mixed at fixed density, and it keeps no memory of earlier
geometries. Its vacuum energy therefore depends on the metric only through the volume. It sees the
lattice’s geometry but not its directions, the graviton stays massless, and its induced action is the
Regge action. A remembered geometry would allow the shape terms of massive gravity. This is why
fields on the verified lattice fail the test and the unverified ocean passes it. The ocean has no halo
around a defect, so the field of a mass lives in the geometry. The bonds carry as much entanglement
per unit area as they can. Their equilibrium with matter is the Hamiltonian constraint, which the
solved defect slice satisfies ball by ball. Maximizing their entanglement per unit area makes 1
/
4
G
its largest value and fixes the bonds’ microscopic law. Since the bonds are frozen ocean, the stiffness
of the liquid and the entanglement of the solid are one quantity. The remaining question is whether
the ocean’s own fields reproduce it.
Data availability
All results reproduce from the Python scripts archived at https:
//github.com/raghu91302/ssmtheory/blob/main/ssmtheory_form_strength_scripts.zip.
The archive’s README lists the script and the command for each table and figure. The scripts
need only NumPy, SciPy, and Matplotlib. All random numbers come from fixed seeds.
Declarations
Funding: This research received no external funding.
Conflict of interest: The author declares no conflict of interest.
14
Use of artificial intelligence: The author used Claude (Anthropic) to help write and run the
numerical scripts and to draft and edit the manuscript. The author reviewed all results and text
and takes full responsibility for the content.
References
[1]
A. D. Sakharov, Vacuum quantum fluctuations in curved space and the theory of gravitation,
Dokl. Akad. Nauk SSSR 177, 70 (1967) [Sov. Phys. Dokl. 12, 1040 (1968)]; reprinted in Gen.
Relativ. Gravit. 32, 365 (2000). doi:10.1023/A:1001947813563
[2]
M. Visser, Sakharov’s induced gravity: a modern perspective, Mod. Phys. Lett. A 17, 977
(2002). arXiv:gr-qc/0204062
[3]
T. Regge, General relativity without coordinates, Nuovo Cimento 19, 558 (1961).
doi:10.1007/BF02733251
[4]
R. Kulkarni, Vacuum fluctuations as the stiffness of spacetime: induced Regge gravity on the
Selection–Stitch lattice, Quantum Rep. 8, 109 (2026). doi:10.3390/quantum8040109
[5]
R. Kulkarni, Matter as incomplete crystallization: quark charges, color confinement, and the
proton mass from a single extra node in the vacuum lattice, Phys. Open 27, 100423 (2026).
doi:10.1016/j.physo.2026.100423
[6]
T. C. Hales, A proof of the Kepler conjecture, Ann. Math. 162, 1065 (2005).
doi:10.4007/annals.2005.162.1065
[7]
R. Kulkarni, A 67%-rate CSS code on the FCC lattice: [[192
,
130
,
3]] from weight-12 stabilizers
(2026). arXiv:2603.20294
[8]
A. R. Calderbank and P. W. Shor, Good quantum error-correcting codes exist, Phys. Rev. A
54, 1098 (1996). doi:10.1103/PhysRevA.54.1098
[9]
R. Kulkarni, The mass–energy–information equivalence: a bottom-up identification of
the particle spectrum via FCC lattice error correction, Phys. Open 27, 100414 (2026).
doi:10.1016/j.physo.2026.100414
[10]
R. Kulkarni, Dark matter as incomplete crystallization: a geometric construction on the octahe-
dral void of the FCC vacuum lattice, Quantum Rep. 8, 89 (2026). doi:10.3390/quantum8030089
[11]
R. Kulkarni, Black holes in the FCC selection–stitch model: Bekenstein–Hawking entropy,
geometric evaporation, and primordial-black-hole signatures, Eur. Phys. J. Plus 141, 916 (2026).
doi:10.1140/epjp/s13360-026-08148-9
[12]
R. Kulkarni, Regge calculus on the
D
4
lattice through second order: exact identities and the
post-Newtonian static field (2026). doi:10.5281/zenodo.20952333
[13]
C. Riek, D. V. Seletskiy, A. S. Moskalenko, J. F. Schmidt, P. Krauspe, S. Eckart, S. Eggert, G.
Burkard, and A. Leitenstorfer, Direct sampling of electric-field vacuum fluctuations, Science
350, 420 (2015). doi:10.1126/science.aac9788
[14]
J. Cheeger, W. M¨uller, and R. Schrader, On the curvature of piecewise flat spaces, Commun.
Math. Phys. 92, 405 (1984). doi:10.1007/BF01210729
15
[15]
M. Fierz and W. Pauli, On relativistic wave equations for particles of arbitrary spin in an
electromagnetic field, Proc. R. Soc. Lond. A 173, 211 (1939). doi:10.1098/rspa.1939.0140
[16]
K. Hinterbichler, Theoretical aspects of massive gravity, Rev. Mod. Phys. 84, 671 (2012).
doi:10.1103/RevModPhys.84.671
[17]
R. Abbott et al. (LIGO Scientific, Virgo, and KAGRA Collaborations), Tests of general
relativity with GWTC-3, Phys. Rev. D 112, 084080 (2025). doi:10.1103/PhysRevD.112.084080
[18]
E. T. Jaynes, Information theory and statistical mechanics, Phys. Rev. 106, 620 (1957).
doi:10.1103/PhysRev.106.620
[19] J. F. C. Kingman, Poisson Processes (Oxford University Press, Oxford, 1993).
[20]
S. Popescu, A. J. Short, and A. Winter, Entanglement and the foundations of statistical
mechanics, Nat. Phys. 2, 754 (2006). doi:10.1038/nphys444
[21]
L. Bombelli, J. Henson, and R. D. Sorkin, Discreteness without symmetry breaking: a theorem,
Mod. Phys. Lett. A 24, 2579 (2009). doi:10.1142/S0217732309031958
[22]
L. Bombelli, J. Lee, D. Meyer, and R. D. Sorkin, Space-time as a causal set, Phys. Rev. Lett.
59, 521 (1987). doi:10.1103/PhysRevLett.59.521
[23]
N. H. Christ, R. Friedberg, and T. D. Lee, Random lattice field theory: general formulation,
Nucl. Phys. B 202, 89 (1982).
[24]
Y. Shamir, Dynamical-space regular-time lattice and induced gravity (1994). arXiv:hep-
th/9402109
[25] C. de Rham, Massive gravity, Living Rev. Relativ. 17, 7 (2014). doi:10.12942/lrr-2014-7
[26]
T. Jacobson, Thermodynamics of spacetime: the Einstein equation of state, Phys. Rev. Lett.
75, 1260 (1995). doi:10.1103/PhysRevLett.75.1260
[27]
T. Jacobson, Entanglement equilibrium and the Einstein equation, Phys. Rev. Lett. 116, 201101
(2016). doi:10.1103/PhysRevLett.116.201101
[28]
N. Lashkari, M. B. McDermott, and M. Van Raamsdonk, Gravitational dynamics from entangle-
ment “thermodynamics”, J. High Energy Phys. 04, 195 (2014). doi:10.1007/JHEP04(2014)195
[29]
T. Faulkner, M. Guica, T. Hartman, R. C. Myers, and M. Van Raamsdonk, Gravi-
tation from entanglement in holographic CFTs, J. High Energy Phys. 03, 051 (2014).
doi:10.1007/JHEP03(2014)051
[30]
C. Cao, S. M. Carroll, and S. Michalakis, Space from Hilbert space: recovering geometry from
bulk entanglement, Phys. Rev. D 95, 024031 (2017). doi:10.1103/PhysRevD.95.024031
[31]
C. Cao and S. M. Carroll, Bulk entanglement gravity without a boundary: to-
wards finding Einstein’s equation in Hilbert space, Phys. Rev. D 97, 086003 (2018).
doi:10.1103/PhysRevD.97.086003
[32]
M. Henneaux and C. Teitelboim, The cosmological constant and general covariance, Phys. Lett.
B 222, 195 (1989). doi:10.1016/0370-2693(89)91251-3
16
[33]
T. Piran and R. M. Williams, Three-plus-one formulation of Regge calculus, Phys. Rev. D 33,
1622 (1986). doi:10.1103/PhysRevD.33.1622
[34]
R. Schoen and S.-T. Yau, On the structure of manifolds with positive scalar curvature,
Manuscripta Math. 28, 159 (1979). doi:10.1007/BF01647970
[35]
M. Gromov and H. B. Lawson, Spin and scalar curvature in the presence of a fundamental
group. I, Ann. Math. 111, 209 (1980). doi:10.2307/1971198
[36]
D. Fattal, T. S. Cubitt, Y. Yamamoto, S. Bravyi, and I. L. Chuang, Entanglement in the
stabilizer formalism (2004). arXiv:quant-ph/0406168
[37]
D. V. Fursaev, Entanglement entropy in critical phenomena and analogue models of quantum
gravity, Phys. Rev. D 73, 124025 (2006). doi:10.1103/PhysRevD.73.124025
[38]
A. Dobado, F. J. Llanes-Estrada, and J. A. Oller, Existence of two-solar-mass neutron
star constrains gravitational constant
G
N
at strong field, Phys. Rev. C 85, 012801 (2012).
doi:10.1103/PhysRevC.85.012801
17