Home / Quantum Gravity / Entropy, action and dynamics

Local rules of change on a sprinkling

Status
Proved: Theorems 1, 3, 8, 9, 10, Corollary 2, Propositions 4, 5, 6, 7, 11, 12; Theorem 1 assumes one rule for all sizes and reversibility. Computed: one rule's error in two dimensions, pairs related by the dominance map, exact sums over small layered classes, variances of smeared actions. Not decided: suppression at Kleitman–Rothschild proportions at depth one or more.
Excludes
Excluded under the stated hypotheses: exactly wave-consistent completion rules on finitely many sprinkled elements in generic position (dimension two or higher), and separate-statistics laws on orders with a symmetry moving a point to one with a different future. Not decided: suppression of Kleitman–Rothschild orders at depth one or more; rules with vanishing error.
Revised

A sprinkling of Minkowski space has no coordinates and no neighbours; it is a partial order. We ask which rules of change can be written on it using finitely many elements at a time. In two dimensions, on the lattice of all pairs of null ranks, an exactly wave-consistent rewriting rule does exist, and the rules whose result does not depend on how a rectangle is cut are the heap operations of groups (Theorem 1). On a Poisson sprinkling, where no two points are null-separated with probability one, no rule that predicts the value of a field at a sprinkled element from its values at finitely many other sprinkled elements is exactly consistent with the massless wave equation, in dimension two or higher (Theorem 3). For path sums, consider orders with three layers of elements in which every element of the first layer precedes every element of the third (section 6). For such orders with outer layers of fixed size and an action built from finitely many interval counts, the sum divided by the number of orders in the class is exponentially small in the size of the middle layer whenever the phase the action gives to one link is not a multiple of 2π2\pi, and the bound on the rate depends on the action only through that phase (Theorems 8 and 9).

A second kind of rule revises the order itself. A law that decides whether xx precedes yy from statistics of xx and of yy taken separately, such as the sizes of their pasts and futures, cannot leave fixed an order that has a symmetry moving a point to a point with a different future (Proposition 5). A Minkowski diamond with a spatial reflection is such an order. A sprinkling with many points is expected to have no such symmetry, so the proposition does not apply to it. For one such law, the dominance map that relates an element to every element with more ancestors and fewer descendants, an open set of directions around every interior point of a smooth density is related to incomparable points (Proposition 6). On sprinklings the computed fraction of incomparable pairs that the map relates is about 0.67 in two dimensions and 0.55 to 0.62 in four.

The sum over orders is not suppressed when the phase the action gives to a link is a multiple of 2π2\pi: exactly so for two layers (Theorem 8), and in the limit of a large middle layer for three (Proposition 11). Almost all orders have three layers of sizes close to n/4,n/2,n/4n/4,n/2,n/4. At these proportions the uniform bound of Theorem 10 is exponentially small in n2n^2 for the link action when ∣cos⁡(βC0/2)∣<1/3\lvert\cos(\beta C_0/2)\rvert<1/\sqrt3, with βC0\beta C_0 the phase per link. For every action of depth at least one, that is with weights on intervals containing one or more elements, the bound exceeds 11 for large nn, so suppression at these proportions is not decided.

A decoherence functional that is a product of one amplitude per history, as in complex transitive percolation, has rank one and does not decohere a binary coarse-graining unless one side has zero measure (Proposition 7). Smearing the interval-count action over a nonlocality scale is a thinning of the sprinkling, and it bounds the variance of the smeared action by the retained fraction times the variance of the unsmeared action at the thinned density (Proposition 12).

Sorkin [Sorkin] argued that Lorentz-invariant discreteness forces non-locality. The d’Alembertians of Benincasa and Dowker [BD] and of Aslanbeigi, Saravani and Sorkin [ASS] are built from interval counts over mesoscopic regions, and Johnston’s propagator [Johnston] sums over chains of the order. The rewriting rule on the null lattice fills the fourth corner of a rectangle from the other three. Its algebra is that of a heap, treated as pure algebra in [ZAB], and in the abelian case it is the parallelogram identity satisfied by solutions of the one-dimensional wave equation [KRud]. Theorem 1 shows that, under two assumptions, no other rule has the cut-independence property. For path sums, Kleitman and Rothschild [KR75] showed that almost all partial orders have three layers. Loomis and Carlip [LC] showed that the Benincasa–Dowker action suppresses bilayer orders, Mathur, Singh and Surya [MSS] that the link action suppresses the Kleitman–Rothschild orders over a range of parameters, and Carlip, Carlip and Surya [CCS1] that these orders are strongly suppressed in the gravitational path integral. Carlip, Carlip and Surya [CCS2] showed that for a typical layered order the action reduces to the link action to leading order. Theorem 9 adds a transfer-matrix bound on the sum over a three-layer class that holds for every finite-depth action and isolates the part of the rate that depends on the action. A map of the quantum-gravity families examined here places these results among other approaches.

1. Setting and three kinds of rule

A sprinkling is a Poisson point process of density ρ\rho in a region of dd-dimensional Minkowski space, with signature mostly plus and the induced causal order. A field theory on it needs a way to propagate values, and a dynamics of the order needs a way to change relations or to weight them.

A completion rule gives, for a target element tt and a finite set SS of other elements, a function that predicts φ(t)\varphi(t) from φ∣S\varphi\vert_S. A pair law takes a finite order and returns a new relation, deciding for each pair x,yx,y whether xx precedes yy in the new relation. The laws of section 4 read only numbers attached to xx and to yy separately. A path-sum weight assigns a phase to each order and leaves the orders themselves unchanged; the dynamics is the sum over orders (section 6).

2. Rectangle rules on the null lattice

In 1+11+1 dimensions use null coordinates u,vu,v. The massless wave equation is ∂u∂vφ=0\partial_u\partial_v\varphi=0, its solutions are φ=f(u)+g(v)\varphi=f(u)+g(v), and on every rectangle with corners b=(u1,v1)b=(u_1,v_1), l=(u1,v2)l=(u_1,v_2), r=(u2,v1)r=(u_2,v_1), t=(u2,v2)t=(u_2,v_2) they satisfy φ(t)=φ(l)+φ(r)−φ(b)\varphi(t)=\varphi(l)+\varphi(r)-\varphi(b). The identity holds for rectangles of every size, so the value at tt can be computed from the other three corners directly, or after cutting the rectangle in two and computing the value at the new corner first. That the two routes agree is the property studied here. We ask which rules t=F(l,b,r)t=F(l,b,r), for values in an arbitrary set XX, have it.

Let Λ={0,…,N}2\Lambda=\{0,\dots,N\}^2 be the null lattice. A rectangle is a quadruple b=(a,c)b=(a,c), l=(a,j)l=(a,j), r=(i,c)r=(i,c), t=(i,j)t=(i,j) with a<ia<i and c<jc<j, of any size. A field is a map from a subset of Λ\Lambda to XX. Given a map F:X3→XF:X^3\to X, a rewriting step picks a rectangle whose corners b,l,rb,l,r carry values and whose corner tt does not, and sets φ(t)=F(φ(l),φ(b),φ(r))\varphi(t)=F(\varphi(l),\varphi(b),\varphi(r)). Two assumptions are made on FF.

  • (H1) The same FF is used for every rectangle, whatever its size.
  • (H2) FF is reversible in the two outer arguments: for fixed bb and rr the map l↦F(l,b,r)l\mapsto F(l,b,r) is a bijection of XX, and for fixed ll and bb the map r↦F(l,b,r)r\mapsto F(l,b,r) is a bijection.

Goursat data are the values of φ\varphi on the two lines {0}×{0,…,N}\{0\}\times\{0,\dots,N\} and {0,…,N}×{0}\{0,\dots,N\}\times\{0\}, agreeing at the origin. A history is a sequence of rewriting steps that fills Λ\Lambda starting from Goursat data. Histories exist, since a point t=(i,j)t=(i,j) with i,j≥1i,j\ge1 and i+ji+j minimal among the empty points has all three other corners of the unit rectangle filled. The rule FF is confluent if for every Goursat data every history ends with the same field on Λ\Lambda.

A heap operation on a group GG is the map F(l,b,r)=l b−1rF(l,b,r)=l\,b^{-1}r. Writing the two ways of cutting a rectangle as identities gives, for all values in XX,

F(F(l,b,m),m,r)=F(l,b,r),(1)F\big(F(l,b,m),m,r\big)=F(l,b,r),\tag{1} F(l,l1,F(l1,b,r))=F(l,b,r).(2)F\big(l,l_1,F(l_1,b,r)\big)=F(l,b,r).\tag{2}

Identity (1) says that cutting a rectangle with bottom edge b→rb\to r along a vertical line through an intermediate bottom point mm and filling the left piece first gives the same value at the top right corner as filling the big rectangle at once. Identity (2) says the same for a cut along a horizontal line through an intermediate point l1l_1 of the left edge.

Theorem 1 (Confluent rectangle rules are heaps). Let F:X3→XF:X^3\to X satisfy (H1) and (H2). The following are equivalent.

(a) FF is confluent.

(b) FF satisfies (1) and (2).

(c) There is a group structure on XX with F(l,b,r)=l b−1rF(l,b,r)=l\,b^{-1}r.

In that case, for Goursat data φ(0,j)=ψL(j)\varphi(0,j)=\psi_L(j) and φ(i,0)=ψB(i)\varphi(i,0)=\psi_B(i), every history produces φ(i,j)=A(j) B(i)\varphi(i,j)=A(j)\,B(i) with A(j)=ψL(j)ψL(0)−1A(j)=\psi_L(j)\psi_L(0)^{-1} and B(i)=ψB(i)B(i)=\psi_B(i). If XX is a manifold and FF is smooth, the group is a Lie group.

Identities (1) and (2) are the two ways of cutting a rectangle, and reversibility turns them into the associativity law of the product x⋅y=F(x,e,y)x\cdot y=F(x,e,y) for a fixed element ee. Both identities are needed, since each alone admits operations that are not heaps (see the counts after the proof).

Proof. (a)⇒\Rightarrow(b). Take N=2N=2 and Goursat values φ(0,0)=b\varphi(0,0)=b, φ(0,1)=l\varphi(0,1)=l, φ(1,0)=m\varphi(1,0)=m, φ(2,0)=r\varphi(2,0)=r, with arbitrary elements of XX also at (0,2)(0,2). One history fills (1,1)(1,1) from the unit rectangle at the origin and then (2,1)(2,1) from the rectangle with bottom edge from (1,0)(1,0) to (2,0)(2,0). It gives F(F(l,b,m),m,r)F(F(l,b,m),m,r) at (2,1)(2,1). Another history fills (2,1)(2,1) at once from the rectangle with bottom corner (0,0)(0,0) and gives F(l,b,r)F(l,b,r). Both extend to histories that fill Λ\Lambda, so confluence forces (1). The transposed array, with data φ(0,0)=b\varphi(0,0)=b, φ(0,1)=l1\varphi(0,1)=l_1, φ(0,2)=l\varphi(0,2)=l, φ(1,0)=r\varphi(1,0)=r and an arbitrary element at (2,0)(2,0), fills (1,2)(1,2) either through (1,1)(1,1) or directly, and gives (2).

(b)⇒\Rightarrow(c). Put r=mr=m in (1). By (H2) the element x=F(l,b,m)x=F(l,b,m) runs over all of XX as ll does, so

F(x,m,m)=xfor all x,m.(3)F(x,m,m)=x\qquad\text{for all }x,m.\tag{3}

Put l=l1l=l_1 in (2). As rr varies, y=F(l1,b,r)y=F(l_1,b,r) runs over XX, so F(l1,l1,y)=yF(l_1,l_1,y)=y for all l1,yl_1,y. Fix e∈Xe\in X and define x⋅y=F(x,e,y)x\cdot y=F(x,e,y). The two identities with the middle argument equal to ee give x⋅e=xx\cdot e=x and e⋅y=ye\cdot y=y. By (H2), for every mm the map l↦l⋅ml\mapsto l\cdot m is a bijection, with inverse x↦x/mx\mapsto x/m, and r↦m⋅rr\mapsto m\cdot r is a bijection with inverse y↦m\yy\mapsto m\backslash y. Now put b=eb=e in (1): F(l⋅m,m,r)=l⋅rF(l\cdot m,m,r)=l\cdot r. Writing x=l⋅mx=l\cdot m this reads

F(x,m,r)=(x/m)⋅r.(4)F(x,m,r)=(x/m)\cdot r.\tag{4}

Put b=eb=e in (2): F(l,l1,l1⋅r)=l⋅rF(l,l_1,l_1\cdot r)=l\cdot r. Writing y=l1⋅ry=l_1\cdot r this reads

F(l,l1,y)=l⋅(l1\y).(5)F(l,l_1,y)=l\cdot(l_1\backslash y).\tag{5}

Comparing (4) and (5) for the same arguments gives (x/m)⋅y=x⋅(m\y)(x/m)\cdot y=x\cdot(m\backslash y). Choose x=a⋅mx=a\cdot m and y=m⋅cy=m\cdot c. Then x/m=ax/m=a and m\y=cm\backslash y=c, so a⋅(m⋅c)=(a⋅m)⋅ca\cdot(m\cdot c)=(a\cdot m)\cdot c for all a,m,ca,m,c. Thus (X,⋅)(X,\cdot) is an associative quasigroup with identity ee, that is a group, with x/m=x m−1x/m=x\,m^{-1}, and (4) is F(x,m,r)=x m−1rF(x,m,r)=x\,m^{-1}r.

(c)⇒\Rightarrow(a). Let FF be the heap operation of a group. Along any history every filled point (i,j)(i,j) carries A(j)B(i)A(j)B(i): this holds on the Goursat lines, since A(0)=eA(0)=e and B(0)=ψB(0)=ψL(0)B(0)=\psi_B(0)=\psi_L(0), and a step with corners A(j)B(a)A(j)B(a), A(c)B(a)A(c)B(a), A(c)B(i)A(c)B(i) produces A(j)B(a) B(a)−1A(c)−1 A(c)B(i)=A(j)B(i)A(j)B(a)\,B(a)^{-1}A(c)^{-1}\,A(c)B(i)=A(j)B(i). The final field is therefore independent of the history.

For the last statement, x⋅y=F(x,e,y)x\cdot y=F(x,e,y) is smooth, and so is the inverse x−1=F(e,x,e)x^{-1}=F(e,x,e). □\square

On sets of n=2,3,4n=2,3,4 elements, an exhaustive search over the maps reversible in the sense of (H2) finds 1,1,41,1,4 that satisfy both (1) and (2), all heap operations, and 1,2,241,2,24 that satisfy (1) alone and 1,2,241,2,24 that satisfy (2) alone, of which only 1,1,41,1,4 are heaps. The number of heap operations on an nn-element set is ∑G(n−1)!/∣Aut G∣\sum_G (n-1)!/\lvert\mathrm{Aut}\,G\rvert over isomorphism classes of groups of order nn.

In a group, a field φ(i,j)=A(j)B(i)\varphi(i,j)=A(j)B(i) has ∂vφ φ−1\partial_v\varphi\,\varphi^{-1} independent of uu, so in the continuum limit the rule reproduces the equation

∂u(∂vφ  φ−1)=0,(6)\partial_u\big(\partial_v\varphi\;\varphi^{-1}\big)=0,\tag{6}

whose solutions are exactly the factorised fields φ=A(v)B(u)\varphi=A(v)B(u). For G=(R,+)G=(\mathbb R,+) it is the free wave equation. For non-abelian GG it is not the principal chiral equation ∂u(φ−1∂vφ)+∂v(φ−1∂uφ)=0\partial_u(\varphi^{-1}\partial_v\varphi)+\partial_v(\varphi^{-1}\partial_u\varphi)=0, which in null coordinates reads 2∂u∂vφ=∂uφ φ−1∂vφ+∂vφ φ−1∂uφ2\partial_u\partial_v\varphi=\partial_u\varphi\,\varphi^{-1}\partial_v\varphi+\partial_v\varphi\,\varphi^{-1}\partial_u\varphi whereas (6) reads ∂u∂vφ=∂vφ φ−1∂uφ\partial_u\partial_v\varphi=\partial_v\varphi\,\varphi^{-1}\partial_u\varphi; the two coincide when the currents ∂uφ φ−1\partial_u\varphi\,\varphi^{-1} and ∂vφ φ−1\partial_v\varphi\,\varphi^{-1} commute, in particular for abelian GG. If X=RX=\mathbb R and FF is smooth, the group is a connected, simply connected one-dimensional Lie group, hence isomorphic to (R,+)(\mathbb R,+) by a diffeomorphism hh, so F=h−1(h(l)−h(b)+h(r))F=h^{-1}(h(l)-h(b)+h(r)): a free wave equation for the redefined field h∘φh\circ\varphi.

Corollary 2 (Linear confluent rules). Let X=RX=\mathbb R and F(l,b,r)=cl l+cb b+cr rF(l,b,r)=c_l\,l+c_b\,b+c_r\,r with real cl,cr≠0c_l,c_r\neq0 and cbc_b. Then FF is confluent if and only if F(l,b,r)=l−b+rF(l,b,r)=l-b+r.

Proof. Expanding (1) gives cl2l+clcbb+(clcr+cb)m+crr=cll+cbb+crrc_l^2l+c_lc_bb+(c_lc_r+c_b)m+c_rr=c_ll+c_bb+c_rr, so cl2=clc_l^2=c_l and clcr+cb=0c_lc_r+c_b=0; with cl≠0c_l\neq0 this gives cl=1c_l=1 and cb=−crc_b=-c_r. Expanding (2) gives cll+(cb+crcl)l1+crcbb+cr2r=cll+cbb+crrc_ll+(c_b+c_rc_l)l_1+c_rc_bb+c_r^2r=c_ll+c_bb+c_rr, so cr2=crc_r^2=c_r, hence cr=1c_r=1 and cb=−1c_b=-1. Conversely l−b+rl-b+r is the heap operation of (R,+)(\mathbb R,+). □\square

A mass term cannot be added to the rule. For (∂u∂v+m2)φ=0(\partial_u\partial_v+m^2)\varphi=0 the combination φ(t)−φ(l)−φ(r)+φ(b)\varphi(t)-\varphi(l)-\varphi(r)+\varphi(b) equals −m2-m^2 times the integral of φ\varphi over the rectangle, and two solutions can share the values at the three corners b,l,rb,l,r and differ in that integral, because φ(t)\varphi(t) depends on the data along the characteristic segments from bb to ll and from bb to rr and not only on their end points. No rule that sees only the corners reproduces it.

Theorem 1 concerns the lattice Λ\Lambda of all pairs of ranks. A two-dimensional sprinkling of nn points is a two-order: each element has integer null ranks (i(x),j(x))∈{1,…,n}2(i(x),j(x))\in\{1,\dots,n\}^2 obtained by ordering the points along each null direction, and x≺yx\prec y exactly when both ranks increase. Every rank occurs once, so no rectangle of the rank lattice has all four corners among the sprinkled elements, since two corners of a rectangle share a rank. The sprinkled elements occupy a fraction 1/n1/n of the lattice sites, so rewriting on the lattice is rewriting on a regular square lattice and not on the sprinkling.

3. Completion rules on sprinkled elements

Let Σ⊂R1,d−1\Sigma\subset\mathbb R^{1,d-1} be a finite set of points, d≥2d\ge2, and let t∈Σt\in\Sigma and V⊂Σ∖{t}V\subset\Sigma\setminus\{t\} be such that tt is not null-separated from any point of VV. A completion rule is any function H:RV→RH:\mathbb R^V\to\mathbb R that may depend on the positions of all points and on the order, and no locality is assumed on the choice of VV. It is exactly consistent with □φ=0\Box\varphi=0 if H(φ∣V)=φ(t)H(\varphi\vert_V)=\varphi(t) for every smooth solution φ\varphi of the wave equation on Minkowski space. The same definition applies to fields with values in a group GG that contains a non-trivial one-parameter subgroup exp⁡(sξ)\exp(s\xi), with the principal chiral equation ∂μ(g−1∂μg)=0\partial^\mu(g^{-1}\partial_\mu g)=0 in place of the wave equation, and with exact consistency meaning H(g∣V)=g(t)H(g\vert_V)=g(t) for all smooth solutions.

Plane waves along a null direction give the counterexample. A function χ(n⋅x)\chi(n\cdot x) with nn null satisfies the wave equation and is constant on each null hyperplane n⋅x=constn\cdot x=\text{const}. The finite set VV meets only finitely many such hyperplanes. If one can choose nn so that the hyperplane through tt contains no point of VV, then χ\chi can vanish on all the hyperplanes through VV and equal 11 on the one through tt.

Theorem 3 (No exact completion rule on sprinkled elements). Let d≥2d\ge2, and let VV and tt be as above. There are smooth solutions φ1,φ2\varphi_1,\varphi_2 of □φ=0\Box\varphi=0 with φ1=0\varphi_1=0, φ2∣V=0\varphi_2\vert_V=0 and φ2(t)=1\varphi_2(t)=1. Consequently no completion rule on (V,t)(V,t) is exactly consistent with the wave equation. The same holds for the principal chiral equation with values in a group GG that contains a non-trivial one-parameter subgroup. For a Poisson sprinkling these hypotheses hold almost surely, simultaneously for all VV and tt.

Proof. Put Δs=xt−xs\Delta_s=x_t-x_s for s∈Vs\in V; each Δs\Delta_s is non-zero and not null by assumption. The null directions are the vectors n=(1,n^)n=(1,\hat n) with n^\hat n on the unit sphere of Rd−1\mathbb R^{d-1} (two points when d=2d=2). For such nn, n⋅Δs=−Δs0+n^⋅Δ⃗sn\cdot\Delta_s=-\Delta_s^0+\hat n\cdot\vec\Delta_s. This vanishes for no n^\hat n if Δs\Delta_s is timelike, because ∣n^⋅Δ⃗s∣≤∣Δ⃗s∣<∣Δs0∣\lvert\hat n\cdot\vec\Delta_s\rvert\le\lvert\vec\Delta_s\rvert<\lvert \Delta_s^0\rvert; and if Δs\Delta_s is spacelike it vanishes on a set of directions that is a proper closed subset of the sphere without interior (empty for d=2d=2 since Δs\Delta_s is not null). A finite union of such sets does not cover the sphere, so there is a null nn with n⋅xt∉{n⋅xs:s∈V}n\cdot x_t\notin\{n\cdot x_s:s\in V\}. Choose a smooth χ:R→R\chi:\mathbb R\to\mathbb R that vanishes at the finitely many numbers n⋅xsn\cdot x_s and equals 11 at n⋅xtn\cdot x_t. Then φ2(x)=χ(n⋅x)\varphi_2(x)=\chi(n\cdot x) satisfies □φ2=nμnμχ′′=0\Box\varphi_2=n^\mu n_\mu \chi''=0, vanishes on VV and equals 11 at tt. If HH were exactly consistent, H(0)=φ1(t)=0H(0)=\varphi_1(t)=0 and H(0)=φ2(t)=1H(0)=\varphi_2(t)=1, a contradiction.

For the group-valued case put g1=eg_1=e and g2(x)=exp⁡(χ(n⋅x)ξ)g_2(x)=\exp\big(\chi(n\cdot x)\xi\big) with χ\chi chosen as before except that χ(n⋅xt)=s0\chi(n\cdot x_t)=s_0 with exp⁡(s0ξ)≠e\exp(s_0\xi)\neq e. Since exp⁡(χξ)\exp(\chi\xi) commutes with its derivative, g2−1∂μg2=χ′nμξg_2^{-1}\partial_\mu g_2=\chi'n_\mu \xi, whose divergence is nμnμχ′′ξ=0n^\mu n_\mu \chi''\xi=0. The two solutions agree on VV and differ at tt.

For the Poisson process, almost surely no two points coincide or are null-separated, and this property holds for all finite subsets at once. □\square

The dimension plays no role, and neither does the choice of VV: the rule may use any finite set of sprinkled elements, however it is selected, and may use the positions as well as the order. The proof uses only generic position. The rectangle rule of Theorem 1 evades it because in d=2d=2 the only null directions are the two lattice directions. In a rectangle the corner rr has the same uu as tt and the corner ll has the same vv, so for each of the two null directions the null hyperplane through tt contains a point of V={b,l,r}V=\{b,l,r\}. The general solution f(u)+g(v)f(u)+g(v) is then determined by the three corners, and the proof breaks down. For sprinkled elements, null separation has probability zero. In d=2d=2 the group-valued solutions exp⁡(f(u)ξ)\exp(f(u)\xi) and exp⁡(f(v)ξ)\exp(f(v)\xi) are factorised fields A(v)B(u)A(v)B(u), so they also solve (6) and the same conclusion holds for rules consistent with (6).

The theorem says nothing about rules whose error tends to zero as the density grows. A rule with a defect that vanishes in the continuum limit, or one that uses a growing number of elements, is not covered. A scaling identity shows what such a rule would have to do.

Proposition 4 (Scaling of the error of an order-only rule). Let a completion rule be defined for each marked finite order (C,≺,t)(C,\prec,t) by a set VV and an affine-equivariant function HH, meaning H(λ φV+λ0)=λ H(φV)+λ0H(\lambda\,\varphi_V+\lambda_0)=\lambda\,H(\varphi_V)+\lambda_0 for real λ,λ0\lambda,\lambda_0, and let both depend on the isomorphism class of (C,≺,t)(C,\prec,t) only. For the linear field φ(x)=a⋅x\varphi(x)=a\cdot x, with a fixed vector aa, on a Poisson sprinkling of density ρ\rho in a causal interval of proper height TT, let eρ(T)e_\rho(T) be the mean of ∣H(φ∣V)−φ(t)∣\lvert H(\varphi\vert_V)-\varphi(t)\rvert over targets. Then

eρ(T)=ρ−1/d e1(ρ1/dT).(7)e_\rho(T)=\rho^{-1/d}\,e_1\big(\rho^{1/d}T\big).\tag{7}

Proof. The dilation x↦ρ1/dxx\mapsto\rho^{1/d}x maps a sprinkling of density ρ\rho in an interval of height TT onto a sprinkling of density 11 in an interval of height ρ1/dT\rho^{1/d}T, preserves the order, hence the selection of VV and the function HH, and turns the data a⋅xa\cdot x into (ρ−1/da)⋅x′(\rho^{-1/d}a)\cdot x'. By equivariance the error is multiplied by ρ−1/d\rho^{-1/d}. □\square

Write θ=ρ1/dT\theta=\rho^{1/d}T for the height in units of the discreteness length. Equation (7) says that the error of an order-only rule at fixed TT tends to zero as ρ→∞\rho\to\infty if and only if e1(θ)=o(θ)e_1(\theta)=o(\theta): the error, measured in discreteness lengths, must grow more slowly than the size of the region in the same units. If the error scales as eρ∝ργe_\rho\propto\rho^{\gamma}, then e1(θ)∝θdγ+1e_1(\theta)\propto\theta^{d\gamma+1}, and a rule whose error stays bounded in discreteness lengths has γ=−1/d\gamma=-1/d.

The rectangle rule has a transcription to order data in two dimensions, which we call the join-heap rule. For a target tt, take two maximal elements l,rl,r of its past (its past links), so that tt is minimal in J+(l)∩J+(r)J^+(l)\cap J^+(r); take bb maximal in J−(l)∩J−(r)J^-(l)\cap J^-(r); and predict φ(t)=φ(l)+φ(r)−φ(b)\varphi(t)=\varphi(l)+\varphi(r)-\varphi(b). We evaluated it on the wave-equation solution φ=u+v\varphi=u+v in null coordinates on the unit square, for targets with u,v>0.5u,v>0.5, averaging over sprinklings, targets and admissible triples (l,r,b)(l,r,b).

points NN 250 500 1000 2000 4000 8000
mean absolute error of the prediction 0.174 0.169 0.156 0.151 0.142 0.134
spread of the prediction over admissible triples for one target 0.168 0.179 0.173 0.175 0.173 0.171

By (7) with d=2d=2, a rule with error bounded in discreteness lengths would give an error falling like N−1/2N^{-1/2}, a factor 5.75.7 over this range; the observed fall is a factor 1.31.3. A least-squares fit of ln⁡e\ln e on ln⁡N\ln N gives γ≈−0.08\gamma\approx-0.08, that is e1(θ)∝θ0.85e_1(\theta)\propto\theta^{0.85}.

The selected elements are not confined to a neighbourhood of tt. The longer null side of the box from bb to tt has mean 0.45,0.39,0.350.45,0.39,0.35 for N=500,2000,8000N=500,2000,8000, while the mean spacing N−1/2N^{-1/2} is 0.045,0.022,0.0110.045,0.022,0.011. The mean area of the box is small (0.077,0.050,0.0390.077,0.050,0.039), because tt is minimal in the common future of ll and rr, but the box is long in one null direction, so the elements involved are far from tt on the scale of the region. These are computed values for one rule in two dimensions with one way of selecting ll and rr.

4. Order-revising laws that read separate statistics

Suppose the relations of an order are to be re-decided by a law that looks, for each pair x,yx,y, at numbers attached to xx and to yy separately, such as the sizes of the past and future cones, the height, or the number of chains through the element, and not at how xx and yy relate to each other. A manifold-like order should be a fixed point of such a law. We show that this is impossible when the order has a symmetry that moves a point to a point with a different future, and then study a concrete law of this kind.

Let (Q,<)(Q,<) be a strict partial order, possibly carrying a measure (as a Minkowski diamond carries Lebesgue measure). An automorphism is an order isomorphism that preserves the measure, if there is one. A statistic is a function s:Q→Γs:Q\to\Gamma with values in a set Γ\Gamma that is invariant under the automorphisms: s(gx)=s(x)s(gx)=s(x) for every automorphism gg. A separate-statistics law is given by ss and a relation Ψ⊆Γ×Γ\Psi\subseteq\Gamma\times\Gamma, and returns the relation x R′ y  ⟺  Ψ(s(x),s(y))x\,R'\,y\iff\Psi(s(x),s(y)). The law fixes (Q,<)(Q,<) if R′R' equals <<.

Proposition 5 (A fixed point has orbit-invariant futures). If a separate-statistics law fixes (Q,<)(Q,<), then for every automorphism gg and all x,y∈Qx,y\in Q, x<yx<y implies gx<ygx<y. Consequently an order that has an automorphism gg, and elements x<yx<y with gx≮ygx\not<y, is not fixed by any separate-statistics law.

Proof. If x<yx<y then Ψ(s(x),s(y))\Psi(s(x),s(y)) holds. Since s(gx)=s(x)s(gx)=s(x), also Ψ(s(gx),s(y))\Psi(s(gx),s(y)), that is gx R′ ygx\,R'\,y, and R′=<R'=< gives gx<ygx<y. □\square

The hypotheses are met by the configurations one wants to be fixed. The Minkowski diamond, with Lebesgue measure and cone volumes as the statistics, is mapped to itself by a spatial reflection gg through a hyperplane containing the time axis. Take an interior point x=(t,x⃗)x=(t,\vec x) with x⃗≠gx⃗\vec x\neq g\vec x, so that gx=(t,gx⃗)gx=(t,g\vec x) is spacelike to xx, and let y=x+ς(1,e^)y=x+\varsigma(1,\hat e) with e^\hat e the unit vector along x⃗−gx⃗\vec x-g\vec x and ς>0\varsigma>0 small. Then yy is in the future of xx, while y−gx=(ς,x⃗−gx⃗+ςe^)y-gx=(\varsigma,\vec x-g\vec x+\varsigma\hat e) is spacelike, so x<yx<y and gx≮ygx\not<y. The same holds in every dimension d≥2d\ge2. In the two-dimensional null lattice [N]2[N]^2 the exchange of the two coordinates is such an automorphism. A randomised law whose output distribution depends only on the pair of statistics has the same limitation in expectation, since Pr⁡(gx R′ y)=Pr⁡(x R′ y)\Pr(gx\,R'\,y)=\Pr(x\,R'\,y).

A sprinkling with many points is expected to have a trivial automorphism group (this is not proved), and the proposition then says nothing about it. What a law does to a sprinkling is a separate question, studied next for one law.

For a finite strict order RR on nn elements let nR−(i)n^-_R(i) and nR+(i)n^+_R(i) be the numbers of elements below and above ii. The cone-count dominance map relates ii to jj when ii has fewer predecessors and more successors than jj:

Φ(R)={(i,j): nR−(i)<nR−(j) and nR+(i)>nR+(j)}.(8)\Phi(R)=\{(i,j):\ n^-_R(i)<n^-_R(j)\ \text{and}\ n^+_R(i)>n^+_R(j)\}.\tag{8}

It is a separate-statistics law with s=(nR−,nR+)s=(n^-_R,n^+_R). Part (a) of the next proposition gives its exact properties on finite orders, and part (b) shows that on a smooth density of points it relates incomparable pairs. The quantities n∓n^\mp become the masses ν∓(x)\nu^\mp(x) of the sets below and above xx for a probability measure ν\nu.

Proposition 6 (Properties of the dominance map). (a) For every strict order RR on nn elements, R⊆Φ(R)R\subseteq\Phi(R), Φ(R)\Phi(R) is a strict order of order dimension at most two (an intersection of two linear orders), the sequence R,Φ(R),Φ2(R),…R,\Phi(R),\Phi^2(R),\dots is increasing and constant after at most n(n−1)/2n(n-1)/2 steps, and R=Φ(R)R=\Phi(R) if and only if (nR−(i)−nR−(j))(nR+(i)−nR+(j))≥0(n^-_R(i)-n^-_R(j))(n^+_R(i)-n^+_R(j))\ge0 for every pair i,ji,j that is incomparable in RR. In particular no order of dimension ≥3\ge3 is a fixed point.

(b) Let ν\nu be a probability measure on the open unit square in null coordinates (u,v)(u,v) with a bounded density hh that is continuous and positive on the open square, with the componentwise order, and let ν−(x)\nu^-(x) and ν+(x)\nu^+(x) be the ν\nu-masses of the sets below and above xx. At an interior point x0=(u0,v0)x_0=(u_0,v_0) put

qu−=∫0v0h(u0,v′) dv′,qv−=∫0u0h(u′,v0) du′,qu+=∫v01h(u0,v′) dv′,qv+=∫u01h(u′,v0) du′.q^-_u=\int_0^{v_0}h(u_0,v')\,dv',\quad q^-_v=\int_0^{u_0}h(u',v_0)\,du',\quad q^+_u=\int_{v_0}^1h(u_0,v')\,dv',\quad q^+_v=\int_{u_0}^1h(u',v_0)\,du'.

For every slope τ\tau with 0<τ<min⁡(qu−/qv−,  qu+/qv+)0<\tau<\min(q^-_u/q^-_v,\;q^+_u/q^+_v) there is δ0>0\delta_0>0 such that for 0<δ<δ00<\delta<\delta_0 the point y=x0+δ(1,−τ)y=x_0+\delta(1,-\tau) is incomparable to x0x_0 and satisfies ν−(x0)<ν−(y)\nu^-(x_0)<\nu^-(y) and ν+(x0)>ν+(y)\nu^+(x_0)>\nu^+(y): the continuum analogue of Φ\Phi relates x0x_0 to yy. For τ>max⁡(qu−/qv−, qu+/qv+)\tau>\max(q^-_u/q^-_v,\,q^+_u/q^+_v) the same holds with the two inequalities reversed, so that yy is related to x0x_0.

Proof. (a) If i<ji<j the down-set of jj contains the down-set of ii and ii itself, so n−(i)<n−(j)n^-(i)<n^-(j); the same argument for up-sets gives n+(i)>n+(j)n^+(i)>n^+(j). Hence R⊆Φ(R)R\subseteq\Phi(R). The relation Φ(R)\Phi(R) is the strict componentwise order on the points (n−(i),−n+(i))∈Z2(n^-(i),-n^+(i))\in\mathbb Z^2, which is irreflexive and transitive, and it is the intersection of two linear orders: order the elements lexicographically by (n−,n+)(n^-,n^+) ascending and by (−n+,−n−)(-n^+,-n^-) ascending, and break exact ties by the index in one order and by the reverse index in the other. So it has dimension at most two. Applying R⊆Φ(R)R\subseteq\Phi(R) to each term shows the sequence is increasing; a strict order has at most n(n−1)/2n(n-1)/2 relations. A fixed point satisfies Φ(R)⊆R\Phi(R)\subseteq R. For a pair i,ji,j, one of (i,j)(i,j) and (j,i)(j,i) lies in Φ(R)\Phi(R) exactly when (n−(i)−n−(j))(n+(i)−n+(j))<0(n^-(i)-n^-(j))(n^+(i)-n^+(j))<0. For a pair comparable in RR the product is negative and the pair is already in RR or its reverse, so only incomparable pairs can violate Φ(R)⊆R\Phi(R)\subseteq R, and they do so exactly when the product is negative. Fixed points have dimension at most two because they equal Φ(R)\Phi(R).

(b) Since hh is bounded and continuous, differentiation under the integral shows that ν−(u,v)=∫0u∫0vh\nu^-(u,v)=\int_0^u\int_0^vh and ν+(u,v)=∫u1∫v1h\nu^+(u,v)=\int_u^1\int_v^1h are C1C^1 on the open square, with ∂uν−=qu−\partial_u\nu^-=q^-_u, ∂vν−=qv−\partial_v\nu^-=q^-_v, ∂uν+=−qu+\partial_u\nu^+=-q^+_u and ∂vν+=−qv+\partial_v\nu^+=-q^+_v at x0x_0. For y=x0+δ(1,−τ)y=x_0+\delta(1,-\tau) the first-order changes are Δν−=δ(qu−−τqv−)+o(δ)\Delta\nu^-=\delta(q^-_u-\tau q^-_v)+o(\delta) and Δν+=δ(−qu++τqv+)+o(δ)\Delta\nu^+=\delta(-q^+_u+\tau q^+_v)+o(\delta), which are positive and negative respectively when τ<qu−/qv−\tau<q^-_u/q^-_v and τ<qu+/qv+\tau<q^+_u/q^+_v. The point yy is incomparable to x0x_0 because its uu-coordinate is larger and its vv-coordinate smaller. The case τ\tau large reverses both signs. □\square

Part (a) leaves open whether the map moves a given order at all, since an order of dimension two with the right statistics could be a fixed point. Part (b) shows that on a smooth density it does move the order, and that the set of promoted directions is open. For the uniform measure on the unit square, ν−=uv\nu^-=uv and ν+=(1−u)(1−v)\nu^+=(1-u)(1-v). Numerical integration of the fraction of incomparable pairs that are related in one orientation or the other gives 0.666670.66667; the identification with 2/32/3 rests on this value alone.

On sprinklings the fraction of incomparable pairs that one application of Φ\Phi relates is, averaged over three samples each,

dimension 2, points nn 250 500 1000 2000
fraction related 0.655 0.664 0.668 0.663
dimension 4, points nn 250 500 1000 1500
fraction related 0.547 0.580 0.614 0.620

Iterating Φ\Phi on sprinklings of n=600n=600 points reaches a fixed point after a sample-dependent number of steps, 3232 in the two-dimensional sample reported here, with the fraction of related pairs rising from 0.5180.518 to 0.9880.988 and the height from 4242 to 257257. In a four-dimensional sample it took 99 steps, with the fraction rising from 0.0980.098 to 0.8740.874 and the height from 99 to 163163.

A law that reads data of the pair itself, such as the number of common predecessors ∣↓x∩↓y∣\lvert\downarrow x\cap\downarrow y\rvert, the number of elements in the interval, or the links between them, is not covered, because a symmetry that preserves the separate statistics of xx and yy need not preserve those.

5. Decoherence for a rank-one functional

Section 4 concerned laws that revise an order. Sequential growth instead builds an order by adding one element at a time, and its quantum version is a candidate dynamics of the order. The question here is whether the decoherence functional of the simplest such model separates histories.

Classical sequential growth builds a causal set one element at a time, with probabilities for the new element’s past [RS]. Complex sequential growth replaces the probabilities by complex amplitudes [SZ]. Complex transitive percolation [DJS] is the case with a single parameter pp, in which a labelled order with N0N_0 links and NrelN_{\mathrm{rel}} relations has amplitude pN0(1−p)(n2)−Nrelp^{N_0}(1-p)^{\binom n2-N_{\mathrm{rel}}}. The amplitude A(E)\mathcal A(E) of an event EE, a set of labelled orders, is the sum of their amplitudes. A decoherence functional D(E,E′)\mathcal D(E,E') measures the interference between two events, and a coarse-graining into disjoint events E,E′E,E' is decoherent when D(E,E′)=0\mathcal D(E,E')=0; the quantum measure of EE is D(E,E)\mathcal D(E,E). Here D\mathcal D is the product D(E,E′)=A(E)‾ A(E′)\mathcal D(E,E')=\overline{\mathcal A(E)}\,\mathcal A(E'), with A(Ω)=1\mathcal A(\Omega)=1 for the whole space Ω\Omega. This product form, and the one-dimensional space of histories it implies, are noted in [DJS], which also discusses the extension of the functional to covariant events together with [SZ]. The consequence for coarse-grainings is the following.

Proposition 7 (A rank-one functional does not decohere a binary coarse-graining). Let D(E,E′)=A(E)‾A(E′)\mathcal D(E,E')=\overline{\mathcal A(E)}\mathcal A(E') with A(Ω)=1\mathcal A(\Omega)=1 and A\mathcal A additive on disjoint events, and let EE be an event with complement EcE^c. If A(E)A(Ec)≠0\mathcal A(E)\mathcal A(E^c)\neq0 then ∣D(E,Ec)∣/D(E,E)D(Ec,Ec)=1\lvert\mathcal D(E,E^c)\rvert/\sqrt{\mathcal D(E,E)\mathcal D(E^c,E^c)}=1. The off-diagonal element D(E,Ec)\mathcal D(E,E^c) vanishes if and only if A(E)∈{0,1}\mathcal A(E)\in\{0,1\}, that is if one of E,EcE,E^c has zero quantum measure.

Proof. ∣D(E,Ec)∣=∣A(E)∣∣A(Ec)∣=D(E,E)D(Ec,Ec)\lvert\mathcal D(E,E^c)\rvert=\lvert\mathcal A(E)\rvert\lvert\mathcal A(E^c)\rvert=\sqrt{\mathcal D(E,E)\mathcal D(E^c,E^c)}; this is the Cauchy–Schwarz inequality with equality in the one-dimensional space C\mathbb C. Since A(Ec)=1−A(E)\mathcal A(E^c)=1-\mathcal A(E), the product A(E)A(Ec)\mathcal A(E)\mathcal A(E^c) vanishes exactly when A(E)=0\mathcal A(E)=0 or A(E)=1\mathcal A(E)=1. □\square

The statement holds for every event, every number of elements and every value of pp, including real pp (where the model is classical percolation seen as an amplitude). Decoherence between events needs a history space of dimension at least two, that is a vector-valued amplitude with at least two independent components. Models of that kind lie outside Proposition 7.

6. Path sums over three-layer orders

A sum over orders, Z=∑OeiS(O)/ℏZ=\sum_{\mathcal O}e^{iS(\mathcal O)/\hbar}, can be dominated by manifold-like orders only if the phases suppress the rest. The difficulty is entropy. By the theorem of Kleitman and Rothschild [KR75], almost all partial orders on nn elements have three layers, of sizes close to n/4n/4, n/2n/2 and n/4n/4 (the Kleitman–Rothschild proportions), and their number is 2n2/4+o(n2)2^{n^2/4+o(n^2)}. Almost all of these orders are very far from a manifold. For a given action the question is whether its phases cancel enough in the sum over the three-layer class that the class does not dominate. For an action built from interval counts the sum over the class can be done exactly.

6.1 The action, the class and two layers

For a finite order, let NkN_k be the number of pairs x≺yx\prec y with exactly kk elements strictly between them, so N0N_0 is the number of links, and let nn be the number of elements. A finite-depth action is

Sℏ=β(ζ n−∑k=0LCkNk),(9)\frac{S}{\hbar}=\beta\Big(\zeta\,n-\sum_{k=0}^{L}C_kN_k\Big),\tag{9}

with real parameters β,ζ,C0,…,CL\beta,\zeta,C_0,\dots,C_L and CL≠0C_L\neq0; the integer LL is the depth of the action. The four-dimensional Benincasa–Dowker action [BD] has L=3L=3, ζ=1\zeta=1 and C=(1,−9,16,−8)C=(1,-9,16,-8), with β\beta proportional to the squared ratio of the discreteness length to the Planck length. The two-dimensional one has L=2L=2, ζ=1\zeta=1 and C=(2,−4,2)C=(2,-4,2). In this convention CkC_k is the coefficient of −Nk-N_k in the action, with metric signature mostly plus and the weight of an order eiS/ℏe^{iS/\hbar}; writing the action as +∑kckNk+\sum_kc_kN_k gives ck=−Ckc_k=-C_k. The CkC_k are the layer coefficients of the d’Alembertian in the Remark of section 5 of horizon counts and the area law, and in the action they multiply −Nk-N_k, as in (9). The two-dimensional values are twice the weights (1,−2,1)(1,-2,1) of the two-dimensional d’Alembertian, the factor being absorbed in the normalisation of β\beta. The coupling β\beta is a free parameter. The link action has only C0≠0C_0\neq0 and depth 00. The link phase is z=e−iβC0z=e^{-i\beta C_0}, the phase the action assigns to one link, and f(k)=e−iβCkf(k)=e^{-i\beta C_k}, with f(k)=1f(k)=1 for k>Lk>L and f(0)=zf(0)=z.

Fix layer sizes (n1,n2,n3)(n_1,n_2,n_3) and call the layers Y1,Y2,Y3Y_1,Y_2,Y_3. An order of the class K(n1,n2,n3)\mathcal K(n_1,n_2,n_3) is given by two bipartite relations P12⊆Y1×Y2P_{12}\subseteq Y_1\times Y_2 and P23⊆Y2×Y3P_{23}\subseteq Y_2\times Y_3 with a≺va\prec v for (a,v)∈P12(a,v)\in P_{12}, v≺cv\prec c for (v,c)∈P23(v,c)\in P_{23}, and a≺ca\prec c for all a∈Y1a\in Y_1, c∈Y3c\in Y_3. It has ∣K∣=2n2(n1+n3)\lvert\mathcal K\rvert=2^{n_2(n_1+n_3)} members, which is 2n2/42^{n^2/4} at proportions (n/4,n/2,n/4)(n/4,n/2,n/4), the leading exponent of the number of all orders. For a∈Y1a\in Y_1 and c∈Y3c\in Y_3 the number of elements strictly between them is kac=∣{v:(a,v)∈P12,(v,c)∈P23}∣k_{ac}=\lvert\{v:(a,v)\in P_{12},(v,c)\in P_{23}\}\rvert. The links are the elements of P12P_{12} and P23P_{23} together with the pairs (a,c)(a,c) with kac=0k_{ac}=0. We write

Z(n1,n2,n3)=∑O∈KeiS(O)/ℏ.Z(n_1,n_2,n_3)=\sum_{\mathcal O\in\mathcal K}e^{iS(\mathcal O)/\hbar}.

With two layers there is nothing between any pair, so only links occur: the bilayer class consists of the orders on Y1∪Y2Y_1\cup Y_2 with a≺va\prec v for (a,v)∈G⊆Y1×Y2(a,v)\in G\subseteq Y_1\times Y_2.

Theorem 8 (Exact bilayer sum). For layers of sizes m1,m2m_1,m_2 and the action (9),

Z=eiβζ(m1+m2) (1+z)m1m2,∣Z∣2m1m2=∣cos⁡βC02∣m1m2.(10)Z=e^{i\beta \zeta(m_1+m_2)}\,(1+z)^{m_1m_2},\qquad\frac{\lvert Z\rvert}{2^{m_1m_2}}=\Big\lvert\cos\frac{\beta C_0}{2}\Big\rvert^{m_1m_2}.\tag{10}

For every finite-depth action the bilayer sum depends on the action only through the link phase βC0\beta C_0 (and a global phase). If βC0∈2πZ\beta C_0\in2\pi\mathbb Z the sum is not suppressed at all.

Proof. In a bilayer every interval is empty, so Nk=0N_k=0 for k≥1k\ge1 and N0=∣G∣N_0=\lvert G\rvert, and eiS/ℏ=eiβζnz∣G∣e^{iS/\hbar}=e^{i\beta\zeta n}z^{\lvert G\rvert}. Summing over all G⊆Y1×Y2G\subseteq Y_1\times Y_2 gives ∑Gz∣G∣=(1+z)m1m2\sum_Gz^{\lvert G\rvert}=(1+z)^{m_1m_2}, and ∣1+z∣=2∣cos⁡(βC0/2)∣\lvert1+z\rvert=2\lvert\cos(\beta C_0/2)\rvert. □\square

6.2 The transfer matrix

In the three-layer class the sum organises itself by the elements of the middle layer. Each v∈Y2v\in Y_2 has a set uv⊆Y1u_v\subseteq Y_1 of lower neighbours and wv⊆Y3w_v\subseteq Y_3 of upper neighbours, chosen independently of the other middle elements, with weight z∣uv∣+∣wv∣z^{\lvert u_v\rvert+\lvert w_v\rvert} for its links. The number kack_{ac} is the number of vv with (a,c)∈uv×wv(a,c)\in u_v\times w_v. The weight of the whole order is therefore

eiS/ℏ=eiβζn  z∣P12∣+∣P23∣∏a∈Y1,c∈Y3f(kac).(11)e^{iS/\hbar}=e^{i\beta\zeta n}\;z^{\lvert P_{12}\rvert+\lvert P_{23}\rvert}\prod_{a\in Y_1,c\in Y_3}f(k_{ac}).\tag{11}

Adding the middle elements one at a time, what is needed to continue is the matrix kˉ=(min⁡(kac,L+1))\bar k=(\min(k_{ac},L+1)) of common-neighbour counts capped at L+1L+1, since f(k)=1f(k)=1 beyond the cap. It takes (L+2)n1n3(L+2)^{n_1n_3} values. A middle element with data (u,w)(u,w) adds 11 to the entries in u×wu\times w, with weight z∣u∣+∣w∣z^{\lvert u\rvert+\lvert w\rvert}. Call MM the matrix of these weights. Its entries are sums of powers of zz. Counts never decrease, so MM is triangular for the componentwise order on states. Its diagonal entry at a state depends only on the set U={(a,c):kˉac≤L}U=\{(a,c):\bar k_{ac}\le L\} of uncapped entries, because a move leaves the state unchanged exactly when (u×w)∩U=∅(u\times w)\cap U=\varnothing. Summing over such moves, and for each uu over the allowed w⊆Y3∖NU(u)w\subseteq Y_3\setminus N_U(u), gives

λU=∑(u,w): (u×w)∩U=∅z∣u∣+∣w∣=∑u⊆Y1z∣u∣(1+z) n3−∣NU(u)∣,NU(u)={c:(a,c)∈U for some a∈u}.(12)\lambda_U=\sum_{(u,w):\,(u\times w)\cap U=\varnothing}z^{\lvert u\rvert+\lvert w\rvert}=\sum_{u\subseteq Y_1}z^{\lvert u\rvert}(1+z)^{\,n_3-\lvert N_U(u)\rvert},\qquad N_U(u)=\{c:(a,c)\in U\text{ for some }a\in u\}.\tag{12}

The final weight of a state is g(kˉ)=∏f(kˉac)g(\bar k)=\prod f(\bar k_{ac}), and Z=eiβζn (Mn2)0⋅ gZ=e^{i\beta\zeta n}\,(M^{n_2})_{0\cdot}\,g, where 00 is the zero state. Let λ∗=max⁡U∣λU∣\lambda_*=\max_U\lvert\lambda_U\rvert, the maximum over all subsets U⊆Y1×Y3U\subseteq Y_1\times Y_3, and κ=λ∗/2n1+n3\kappa=\lambda_*/2^{n_1+n_3}. The quantities λU\lambda_U, λ∗\lambda_* and κ\kappa depend on zz, n1n_1 and n3n_3 only. They do not depend on ζ\zeta, on C1,…,CLC_1,\dots,C_L, or on LL.

6.3 Suppression at fixed outer layers

The weights C1,…,CLC_1,\dots,C_L do not enter the entries of MM, which depend on zz alone; they enter through the cap L+1L+1, which fixes the set of states, and through the final weights gg. The exponential rate in the bound is set by the diagonal entries of MM, which see only the links through zz, and the weights beyond the link enter the polynomial prefactor of the bound.

Theorem 9 (Suppression at fixed outer layers). Fix n1,n3≥1n_1,n_3\ge1 and a finite-depth action (9). Let κ=max⁡U∣λU∣/2n1+n3\kappa=\max_U\lvert\lambda_U\rvert/2^{n_1+n_3} with λU\lambda_U as in (12), J=(L+1)n1n3J=(L+1)n_1n_3, and μ=max⁡(∣cos⁡(βC0/2)∣,12)\mu=\max(\lvert\cos(\beta C_0/2)\rvert,\tfrac12). Then

κ≤μmax⁡(n1,n3),κ≥∣cos⁡βC02∣n1+n3,(13)\kappa\le\mu^{\max(n_1,n_3)},\qquad\kappa\ge\Big\lvert\cos\frac{\beta C_0}{2}\Big\rvert^{n_1+n_3},\tag{13}

and for κ>0\kappa>0 and all n2≥1n_2\ge1

∣Z(n1,n2,n3)∣∣K∣≤(J+1) n2 J κ n2−J.(14)\frac{\lvert Z(n_1,n_2,n_3)\rvert}{\lvert\mathcal K\rvert}\le(J+1)\,n_2^{\,J}\,\kappa^{\,n_2-J}.\tag{14}

If κ=0\kappa=0 then Z=0Z=0 for n2>Jn_2>J. In particular, if ∣cos⁡(βC0/2)∣<1\lvert\cos(\beta C_0/2)\rvert<1 then κ<1\kappa<1 and the sum over the class is exponentially small in n2n_2. The exponential factor in (14) depends on the action only through z=e−iβC0z=e^{-i\beta C_0}. The action enters the prefactor, a polynomial in n2n_2 of degree at most JJ.

The bound (14) is below 11 only once (n2−J)ln⁡(1/κ)(n_2-J)\ln(1/\kappa) exceeds Jln⁡n2+ln⁡(J+1)J\ln n_2+\ln(J+1). For the four-dimensional Benincasa–Dowker action (L=3L=3) with n1=n3=3n_1=n_3=3 and ∣cos⁡(βC0/2)∣=0.4\lvert\cos(\beta C_0/2)\rvert=0.4 one has J=36J=36 and κ=0.0623\kappa=0.0623, and this needs n2n_2 of order 10210^2.

Proof. Triangularity was shown above: moves only increase entries, so Mσσ′≠0M_{\sigma\sigma'}\neq0 implies σ′≥σ\sigma'\ge\sigma componentwise. A term of (Mn2)0kˉ(M^{n_2})_{0\bar k} is a product along a path of states 0=σ0≤σ1≤⋯≤σn20=\sigma_0\le\sigma_1\le\dots\le\sigma_{n_2}. A step that changes the state raises at least one entry by 11, so it raises the sum of the entries by at least 11; that sum of the entries lies between 00 and (L+1)n1n3=J(L+1)n_1n_3=J, hence a path has at most JJ steps that change the state. At the n2−jn_2-j steps that do not change the state the factor is a diagonal entry, of modulus at most λ∗\lambda_*. At a step that changes the state, the entries from the current state to all other states are sums over distinct moves of unimodular terms, so their moduli add up to at most the number 2n1+n32^{n_1+n_3} of moves. There are (n2j)\binom{n_2}{j} choices of the positions of the jj changing steps, and therefore, summing over the destination of each changing step and over the final state,

∑kˉ∣(Mn2)0kˉ∣≤∑j=0min⁡(J,n2)(n2j)2(n1+n3)jλ∗ n2−j.\sum_{\bar k}\big\lvert(M^{n_2})_{0\bar k}\big\rvert\le\sum_{j=0}^{\min(J,n_2)}\binom{n_2}{j}2^{(n_1+n_3)j}\lambda_*^{\,n_2-j}.

Since ∣g∣=1\lvert g\rvert=1, ∣Z∣\lvert Z\rvert is at most this sum. Dividing by ∣K∣=2n2(n1+n3)\lvert\mathcal K\rvert=2^{n_2(n_1+n_3)} turns each term into (n2j)κ n2−j\binom{n_2}{j}\kappa^{\,n_2-j}, and with κ≤1\kappa\le1, (n2j)≤n2 j≤n2 J\binom{n_2}{j}\le n_2^{\,j}\le n_2^{\,J} and κ n2−j≤κ n2−J\kappa^{\,n_2-j}\le\kappa^{\,n_2-J} the sum over the at most J+1J+1 values of jj gives (14). If κ=0\kappa=0, every diagonal entry vanishes, a chain of states has at most J+1J+1 members, and Mn2=0M^{n_2}=0 for n2>Jn_2>J.

For the upper bound in (13), ∣1+z∣=2∣cos⁡(βC0/2)∣≤2μ\lvert1+z\rvert=2\lvert\cos(\beta C_0/2)\rvert\le2\mu, and ∣NU(u)∣≤n3\lvert N_U(u)\rvert\le n_3, so each term of (12) satisfies

∣z∣u∣(1+z)n3−∣NU(u)∣∣2n3=∣cos⁡βC02∣n3−∣NU(u)∣2−∣NU(u)∣≤μn3,\frac{\lvert z^{\lvert u\rvert}(1+z)^{n_3-\lvert N_U(u)\rvert}\rvert}{2^{n_3}}=\Big\lvert\cos\frac{\beta C_0}{2}\Big\rvert^{n_3-\lvert N_U(u)\rvert}2^{-\lvert N_U(u)\rvert}\le\mu^{n_3},

since both ∣cos⁡∣≤μ\lvert\cos\rvert\le\mu and 12≤μ\tfrac12\le\mu. The sum has 2n12^{n_1} terms, so ∣λU∣/2n1+n3≤μn3\lvert\lambda_U\rvert/2^{n_1+n_3}\le\mu^{n_3}. The first form of λU\lambda_U in (12) is symmetric in Y1Y_1 and Y3Y_3, so summing over uu first instead gives μn1\mu^{n_1}, hence κ≤μmax⁡(n1,n3)\kappa\le\mu^{\max(n_1,n_3)}. For the lower bound, U=∅U=\varnothing has NU(u)=∅N_U(u)=\varnothing for all uu, so λ∅=(1+z)n3∑uz∣u∣=(1+z)n3+n1\lambda_\varnothing=(1+z)^{n_3}\sum_uz^{\lvert u\rvert}=(1+z)^{n_3+n_1}. □\square

For ∣cos⁡(βC0/2)∣=0.4\lvert\cos(\beta C_0/2)\rvert=0.4, so μ=1/2\mu=1/2, (12) gives λ∗=1.510, 2.280, 3.986\lambda_*=1.510,\,2.280,\,3.986 for n1=n3=1,2,3n_1=n_3=1,2,3, hence exponential rates −ln⁡κ=0.974, 1.948, 2.776-\ln\kappa=0.974,\,1.948,\,2.776 per middle element. The exact sums, computed by enumerating the middle layer by the multiset of its neighbourhood types, decrease with the following local slopes −Δln⁡(∣Z∣/∣K∣)/Δn2-\Delta\ln(\lvert Z\rvert/\lvert\mathcal K\rvert)/\Delta n_2.

layers action local slopes for n2=8→24n_2=8\to24 (step 4) rate −ln⁡κ-\ln\kappa
(1,n2,1)(1,n_2,1) link 0.974,0.974,0.974,0.9740.974,0.974,0.974,0.974 0.9740.974
(1,n2,1)(1,n_2,1) Benincasa–Dowker, d=4d=4 0.585,0.739,0.802,0.8360.585,0.739,0.802,0.836 0.9740.974
(2,n2,2)(2,n_2,2) link 2.10,2.07,1.77,2.032.10,2.07,1.77,2.03 1.9481.948
(2,n2,2)(2,n_2,2) Benincasa–Dowker, d=4d=4 0.91,1.17,1.35,1.490.91,1.17,1.35,1.49 1.9481.948
(3,n2,3)(3,n_2,3) link 2.95,2.91,2.90,2.832.95,2.91,2.90,2.83 (step 2, n2=8→16n_2=8\to16) 2.7762.776
(3,n2,3)(3,n_2,3) Benincasa–Dowker, d=4d=4 0.80,1.23,1.47,1.710.80,1.23,1.47,1.71 (step 2, n2=8→16n_2=8\to16) 2.7762.776

For the link action the local slopes lie close to the rate. Theorem 9 bounds ∣Z∣/∣K∣\lvert Z\rvert/\lvert\mathcal K\rvert with a polynomial prefactor, so a local slope can fall below the rate while the asymptotic decay is not slower; the entry 1.771.77 in the (2,n2,2)(2,n_2,2) row against 1.9481.948 is of this kind. For the Benincasa–Dowker coefficients the local slopes are below the rate and rising, as a prefactor of large degree produces.

The theorem gives the rate only as an upper bound for ∣Z∣/∣K∣\lvert Z\rvert/\lvert\mathcal K\rvert. Whether it is attained for a given action depends on whether the dominant term has a non-zero coefficient, which the theorem does not decide. For the Benincasa–Dowker rows with layers (2,n2,2)(2,n_2,2) and (3,n2,3)(3,n_2,3) (J=16J=16 and 3636) the bound (14) is not yet below 11 at these sizes, so the table illustrates κ\kappa as a rate and does not test (14).

The rate in (14) is for fixed n1n_1 and n3n_3. At Kleitman–Rothschild proportions all three layer sizes grow together, and the next result covers all sizes at once.

6.4 A uniform bound and its logarithmic gap

Theorem 10 (Uniform bound, and its logarithmic gap). Let BL(ℓ)=∑j=0L(ℓj)B_L(\ell)=\sum_{j=0}^{L}\binom{\ell}{j} and μ\mu as above. For all n1,n2,n3n_1,n_2,n_3 and every finite-depth action,

∣Z∣∣K∣≤μ n2n3(Eℓ∼Bin(n2,1/2)[(1+2BL(ℓ))n3])n1.(15)\frac{\lvert Z\rvert}{\lvert\mathcal K\rvert}\le\mu^{\,n_2n_3}\Big(\mathbb E_{\ell\sim\mathrm{Bin}(n_2,1/2)}\big[(1+2B_L(\ell))^{n_3}\big]\Big)^{n_1}.\tag{15}

(a) For the link action (L=0L=0, B0=1B_0=1) the right-hand side is (μn23n1)n3\big(\mu^{n_2}3^{n_1}\big)^{n_3}. At proportions (n/4,n/2,n/4)(n/4,n/2,n/4) it is exp⁡[n28ln⁡(3 μ)]\exp\big[\tfrac{n^2}{8}\ln(\sqrt3\,\mu)\big], exponentially small in n2n^2 whenever ∣cos⁡(βC0/2)∣<1/3\lvert\cos(\beta C_0/2)\rvert<1/\sqrt3.

(b) Let L≥1L\ge1 and ∣cos⁡(βC0/2)∣<1\lvert\cos(\beta C_0/2)\rvert<1. If n1(Lln⁡(n2+1)+ln⁡3)≤(1−δ) n2ln⁡(1/μ)n_1\big(L\ln(n_2+1)+\ln3\big)\le(1-\delta)\,n_2\ln(1/\mu) for some δ>0\delta>0, the right-hand side of (15) is at most exp⁡(−δ n2n3ln⁡(1/μ))\exp(-\delta\,n_2n_3\ln(1/\mu)). For n1=o(n/ln⁡n)n_1=o(n/\ln n) with n2,n3n_2,n_3 proportional to nn this gives e−cn2e^{-cn^2} for large nn.

(c) Let L≥1L\ge1 and ℓ0=⌈n2/2⌉\ell_0=\lceil n_2/2\rceil. The right-hand side of (15) is at least 2−n1μn2n3(1+2BL(ℓ0))n1n32^{-n_1}\mu^{n_2n_3}(1+2B_L(\ell_0))^{n_1n_3}, and at proportions (n/4,n/2,n/4)(n/4,n/2,n/4) it exceeds 11 for large nn.

The weights of the interval counts k=1,…,Lk=1,\dots,L are unimodular but not summable, so the bound has to control them without using cancellation. The device is to expand the product over a∈Y1a\in Y_1 in (11) as a sum over subsets, so that each term depends on the neighbourhood of the middle layer only through a bounded number of constraints.

Proof. Fix P12P_{12} and write Pa⊆Y2P_a\subseteq Y_2 for the upper neighbourhood of a∈Y1a\in Y_1, ℓa=∣Pa∣\ell_a=\lvert P_a\rvert. In (11) the columns c∈Y3c\in Y_3 are independent given P12P_{12}, so Z=eiβζn∑P12z∣P12∣W(P12)n3Z=e^{i\beta\zeta n}\sum_{P_{12}}z^{\lvert P_{12}\rvert}W(P_{12})^{n_3} with

W(P12)=∑ω⊆Y2z∣ω∣∏a∈Y1f(∣Pa∩ω∣).W(P_{12})=\sum_{\omega\subseteq Y_2}z^{\lvert \omega\rvert}\prod_{a\in Y_1}f(\lvert P_a\cap \omega\rvert).

Here ω\omega is the set of middle elements below a given top element. Write f=1+gf=1+g, where g(k)=0g(k)=0 for k>Lk>L and ∣g∣≤2\lvert g\rvert\le2, and expand ∏a(1+g)=∑I⊆Y1∏a∈Ig(∣Pa∩ω∣)\prod_a(1+g)=\sum_{I\subseteq Y_1}\prod_{a\in I}g(\lvert P_a\cap \omega\rvert). For fixed II the product depends on ω\omega only through ω∩PI\omega\cap P_I, with PI=⋃a∈IPaP_I=\bigcup_{a\in I}P_a. Summing over the free part of ω\omega outside PIP_I gives (1+z)n2−∣PI∣(1+z)^{n_2-\lvert P_I\rvert}. The sum over ω′=ω∩PI\omega'=\omega\cap P_I has non-zero terms only if ∣Pa∩ω′∣≤L\lvert P_a\cap \omega'\rvert\le L for all a∈Ia\in I. The map ω′↦(ω′∩Pa)a∈I\omega'\mapsto(\omega'\cap P_a)_{a\in I} is injective, so at most ∏a∈IBL(ℓa)\prod_{a\in I}B_L(\ell_a) terms survive, each of modulus at most 2∣I∣2^{\lvert I\rvert}. With ∣1+z∣n2−∣PI∣≤(2μ)n2\lvert1+z\rvert^{n_2-\lvert P_I\rvert}\le(2\mu)^{n_2}, because 2μ≥12\mu\ge1, we get

∣W(P12)∣≤(2μ)n2∑I∏a∈I2BL(ℓa)=(2μ)n2∏a∈Y1(1+2BL(ℓa)).\lvert W(P_{12})\rvert\le(2\mu)^{n_2}\sum_I\prod_{a\in I}2B_L(\ell_a)=(2\mu)^{n_2}\prod_{a\in Y_1}\big(1+2B_L(\ell_a)\big).

Dividing by 2n22^{n_2}, raising to the power n3n_3 and averaging over uniform P12P_{12}, for which the row sums ℓa\ell_a are independent and binomial with parameters n2n_2 and 1/21/2, gives (15). For (a), B0=1B_0=1.

For (b), BL(ℓ)≤(ℓ+1)L≤(n2+1)LB_L(\ell)\le(\ell+1)^L\le(n_2+1)^L, so 1+2BL(ℓ)≤3(n2+1)L1+2B_L(\ell)\le3(n_2+1)^L and the right-hand side of (15) is at most exp⁡[ n3(n2ln⁡μ+n1(Lln⁡(n2+1)+ln⁡3)) ]\exp[\,n_3(n_2\ln\mu+n_1(L\ln(n_2+1)+\ln3))\,], which is at most exp⁡(−δn2n3ln⁡(1/μ))\exp(-\delta n_2n_3\ln(1/\mu)) under the stated condition.

For (c), BLB_L is non-decreasing and Pr⁡[Bin(n2,12)≥ℓ0]≥12\Pr[\mathrm{Bin}(n_2,\tfrac12)\ge \ell_0]\ge\tfrac12, so E[(1+2BL(ℓ))n3]≥12(1+2BL(ℓ0))n3\mathbb E[(1+2B_L(\ell))^{n_3}]\ge\tfrac12(1+2B_L(\ell_0))^{n_3}, which gives the stated lower bound. Since μ≥12\mu\ge\tfrac12, its logarithm is at least n3[n1ln⁡(1+2BL(ℓ0))−n2ln⁡2]−n1ln⁡2n_3\big[n_1\ln(1+2B_L(\ell_0))-n_2\ln2\big]-n_1\ln2. At proportions (n/4,n/2,n/4)(n/4,n/2,n/4) we have ℓ0=n/4\ell_0=n/4 and BL(ℓ0)≥(n/4L)≥(n/4L)LB_L(\ell_0)\ge\binom{n/4}{L}\ge(n/4L)^L, so this is at least n216[Lln⁡n4L−2ln⁡2]−n4ln⁡2\tfrac{n^2}{16}\big[L\ln\tfrac n{4L}-2\ln2\big]-\tfrac n4\ln2, which tends to +∞+\infty. □\square

Part (b) states what the bound does give and part (c) what it cannot give: the exponent of the right-hand side of (15) carries a term n1n3Lln⁡n2n_1n_3L\ln n_2 against the term n2n3ln⁡(1/μ)n_2n_3\ln(1/\mu) that produces the decay, which is the logarithmic gap.

The bound (15) is far from tight at the sizes where it can be compared with the exact sum. For the link action at ∣cos⁡∣=0.4\lvert\cos\rvert=0.4 it is 7.4×10−117.4\times10^{-11} against the exact 5.0×10−145.0\times10^{-14} on layers (1,10,4)(1,10,4). For the Benincasa–Dowker coefficients it is 1.5×10−41.5\times10^{-4} against 2.4×10−92.4\times10^{-9} on the same layers, and it exceeds 11 on the other four layer sets examined, (2,12,4)(2,12,4), (2,16,6)(2,16,6), (3,12,3)(3,12,3) and (4,8,4)(4,8,4), and on all five at 0.80.8.

The exact value of ∣Z∣/∣K∣\lvert Z\rvert/\lvert\mathcal K\rvert at (m,2m,m)(m,2m,m) with ∣cos⁡(βC0/2)∣=0.4\lvert\cos(\beta C_0/2)\rvert=0.4 (so βC0=2arccos⁡0.4\beta C_0=2\arccos0.4) is:

n=4mn=4m link action Benincasa–Dowker, d=4d=4 0.4n2/40.4^{n^2/4}
4 2.59×10−12.59\times10^{-1} 1.90×10−11.90\times10^{-1} 2.6×10−22.6\times10^{-2}
8 9.54×10−39.54\times10^{-3} 5.51×10−45.51\times10^{-4} 4.3×10−74.3\times10^{-7}
12 4.16×10−54.16\times10^{-5} 6.48×10−66.48\times10^{-6} 4.7×10−154.7\times10^{-15}
16 1.86×10−81.86\times10^{-8} 2.94×10−82.94\times10^{-8} 3.4×10−263.4\times10^{-26}

The last column is the value (10) would give if each of the n2/4n^2/4 relation slots n2(n1+n3)n_2(n_1+n_3) contributed a factor ∣cos⁡∣\lvert\cos\rvert. All entries are suppressed, much more weakly than that. The data stop at n=16n=16 and do not decide the asymptotic form.

A trivial link phase removes the suppression in the limit of a large middle layer, whatever the other weights are.

Proposition 11 (Trivial link phase). If βC0∈2πZ\beta C_0\in2\pi\mathbb Z, then for all n1,n2,n3n_1,n_2,n_3

∣Z∣∣K∣≥1−2 n1n3 Pr⁡[Bin(n2,14)≤L].\frac{\lvert Z\rvert}{\lvert\mathcal K\rvert}\ge1-2\,n_1n_3\,\Pr\big[\mathrm{Bin}(n_2,\tfrac14)\le L\big].

For fixed n1,n3n_1,n_3 the right-hand side tends to 11 as n2→∞n_2\to\infty, and it also tends to 11 at proportions (n/4,n/2,n/4)(n/4,n/2,n/4).

Proof. With z=1z=1, (11) reduces to eiβζn∏a,cf(kac)e^{i\beta\zeta n}\prod_{a,c}f(k_{ac}), and the product equals 11 unless some kac≤Lk_{ac}\le L. Under the uniform measure on K\mathcal K each kack_{ac} is binomial with parameters n2n_2 and 1/41/4. Then ∣E∏f−1∣≤2Pr⁡[∃(a,c):kac≤L]≤2n1n3Pr⁡[Bin(n2,14)≤L]\lvert\mathbb E\prod f-1\rvert\le2\Pr[\exists(a,c):k_{ac}\le L]\le2n_1n_3\Pr[\mathrm{Bin}(n_2,\tfrac14)\le L], and a Chernoff bound makes the probability e−cn2e^{-cn_2}. In both limits n1n3n_1n_3 grows at most polynomially in n2n_2. □\square

For two layers the statement is exact (Theorem 8). For three layers the ratio is bounded below by a quantity that tends to 11 only as n2n_2 grows. With (C0,C1,C2,C3)=(0,−9,16,−8)(C_0,C_1,C_2,C_3)=(0,-9,16,-8), ζ=1\zeta=1, β=1\beta=1 and layers (2,n2,2)(2,n_2,2), the exact value of ∣Z∣/∣K∣\lvert Z\rvert/\lvert\mathcal K\rvert is 0.097,0.111,0.881,0.99990.097,0.111,0.881,0.9999 for n2=8,16,32,64n_2=8,16,32,64, and the bound of the proposition is negative at 88 and 1616 and equals 0.800.80 at 3232. These weights have a trivial link phase by construction and are not claimed to have a continuum limit.

6.6 Relation to earlier results on layered orders

Theorem 9 shows that at fixed outer layers the weights C1,…,CLC_1,\dots,C_L do not enter the exponential rate bound, which is in line with the finding in [LC] for bilayers and in [CCS2] for typical layered orders that the link term governs the leading suppression. The statement of [CCS2] evaluates the action on a typical member of the class, whereas the sums above weight every order of the class with its phase. The transfer-matrix bounds do not establish the reduction to the link action at Kleitman–Rothschild proportions, and the prefactor that Theorem 9 leaves free is visible in the exact data above. The results of [MSS] and [CCS1], that the link action and the gravitational action suppress the Kleitman–Rothschild orders over a range of parameters, concern a different quantity from the full phase sum of this section.

7. Smeared actions and the variance of the action

The sums of section 6 weight orders by the phase eiS/ℏe^{iS/\hbar}. For manifold-like orders, the average of this phase over sprinklings of one region is, heuristically and for an approximately Gaussian SS, small when SS has a large variance, so the variance of the action bears on whether such orders contribute to a path sum. We ask how it depends on the nonlocality scale. Write F=ζN−∑kCkNk\mathcal F=\zeta N-\sum_kC_kN_k for the raw functional in (9), where NN is the number of elements, so that in d=4d=4 the action is a constant times ρ−1/2F\rho^{-1/2}\mathcal F. For a sprinkling of density ρ\rho and a nonlocality density ρk=ϵρ\rho_k=\epsilon\rho with 0<ϵ≤10<\epsilon\le1, the smeared action of [BD] replaces the weights CkC_k by their binomial averages C(ϵ)(m)=∑kCk(mk)ϵk(1−ϵ)m−kC^{(\epsilon)}(m)=\sum_{k}C_k\binom mk\epsilon^k(1-\epsilon)^{m-k} over the number mm of elements in the interval:

Sϵ=ϵ ζ N−ϵ2∑x≺yC(ϵ)(mxy),(16)S_\epsilon=\epsilon\,\zeta\,N-\epsilon^2\sum_{x\prec y}C^{(\epsilon)}(m_{xy}),\tag{16}

up to the same prefactor as the unsmeared action; mxym_{xy} is the number of elements strictly between xx and yy.

Proposition 12 (Smearing is thinning, and damps the variance). Let η\eta be a Poisson process of density ρ\rho on a bounded region, η′\eta' an independent ϵ\epsilon-thinning of η\eta (each point is kept independently with probability ϵ\epsilon), so that η′\eta' is Poisson of density ϵρ\epsilon\rho, and let F\mathcal F be the raw functional above. Then (i) Sϵ(η)=E[F(η′)∣η]S_\epsilon(\eta)=\mathbb E[\mathcal F(\eta')\mid\eta], and (ii) Var Sϵ≤ϵ Var F(η′)\mathrm{Var}\,S_\epsilon\le\epsilon\,\mathrm{Var}\,\mathcal F(\eta').

The proof of (i) is direct. The proof of (ii) uses the chaos expansion of Last and Penrose [LP], Var G=∑k≥11k!∥E DkG∥2\mathrm{Var}\,\mathcal G=\sum_{k\ge1}\frac1{k!}\lVert\mathbb E\,D^k\mathcal G\rVert^2, where DyG(η)=G(η+δy)−G(η)D_y\mathcal G(\eta)=\mathcal G(\eta+\delta_y)-\mathcal G(\eta) is the add-one-point difference and the norm is that of L2L^2 of the kk-fold product of the intensity measure; it holds for every square-integrable functional.

Proof. (i) A pair of η\eta survives the thinning with probability ϵ2\epsilon^2, and given that, the number of surviving points of its interval is binomial with parameters mxym_{xy} and ϵ\epsilon, so the expected contribution to −∑C-\sum C of the pair is −ϵ2C(ϵ)(mxy)-\epsilon^2C^{(\epsilon)}(m_{xy}); the number of points contributes ϵζN\epsilon\zeta N.

(ii) A point yy added to η\eta survives the thinning with probability ϵ\epsilon, so DyG=ϵ E[DyF(η′)∣η]D_y\mathcal G=\epsilon\,\mathbb E[D_y\mathcal F(\eta')\mid\eta] for G=E[F(η′)∣η]\mathcal G=\mathbb E[\mathcal F(\eta')\mid\eta]. Iterating, with independent survival indicators for the added points, E DkG(y1,…,yk)=ϵk TkF(y1,…,yk)\mathbb E\,D^k\mathcal G(y_1,\dots,y_k)=\epsilon^k\,T_k\mathcal F(y_1,\dots,y_k), where TkF=E DkF(η′)T_k\mathcal F=\mathbb E\,D^k\mathcal F(\eta') is the corresponding quantity for the process η′\eta' of intensity measure ϵρ dx\epsilon\rho\,dx. Then ∥ϵkTkF∥L2((ρdx)k)2=ϵk∥TkF∥L2((ϵρ dx)k)2\lVert\epsilon^kT_k\mathcal F\rVert^2_{L^2((\rho dx)^k)}=\epsilon^k\lVert T_k\mathcal F\rVert^2_{L^2((\epsilon\rho\,dx)^k)}, and summing over k≥1k\ge1 with ϵk≤ϵ\epsilon^k\le\epsilon gives Var G≤ϵ Var F(η′)\mathrm{Var}\,\mathcal G\le\epsilon\,\mathrm{Var}\,\mathcal F(\eta'). Square integrability holds because ∣F∣\lvert \mathcal F\rvert is bounded by a polynomial in the number of points of the bounded region. □\square

So at fixed nonlocality density ρk\rho_k the variance of the smeared action is at most ρk/ρ\rho_k/\rho times the variance of the unsmeared action at density ρk\rho_k, for every finite-depth kernel.

For the four-dimensional Benincasa–Dowker functional F=N−N0+9N1−16N2+8N3\mathcal F=N-N_0+9N_1-16N_2+8N_3 on a causal interval of unit height with mean point number Nˉ=250,500,1000,2000\bar N=250,500,1000,2000, the variance of the raw functional is

Nˉ\bar N 250 500 1000 2000
Var F\mathrm{Var}\,\mathcal F 1.6×1051.6\times10^{5} 6.1×1056.1\times10^{5} 2.4×1062.4\times10^{6} 8.1×1068.1\times10^{6}
relative statistical error of Var F\mathrm{Var}\,\mathcal F 6% 8% 12% 18%

so Var F\mathrm{Var}\,\mathcal F grows like Nˉ1.9\bar N^{1.9} over this range. With S∝ρ−1/2FS\propto\rho^{-1/2}\mathcal F the variance of the action grows like Var F/Nˉ\mathrm{Var}\,\mathcal F/\bar N, a fitted power Nˉ0.89\bar N^{0.89}; the asymptotic exponent is not determined by these data. At Nˉ=1000\bar N=1000 and ϵ=1/4\epsilon=1/4 the sampled variance of SϵS_\epsilon (in units of the raw functional) is of order 10310^{3}, against ϵ Var F(η′)≈4×104\epsilon\,\mathrm{Var}\,\mathcal F(\eta')\approx4\times10^{4}, so the inequality (ii) is far from tight.

8. Rules with vanishing error and suppression of layered orders

Theorem 3 leaves open rules whose error tends to zero as ρ→∞\rho\to\infty, among them the interval-count d’Alembertian [BD], whose nonlocality scale is a free input. Constraints on operators of this kind are derived in the four-dimensional causal-set d’Alembertian has no positive spectral representation, for the four-dimensional operators of the form considered there, and in horizon-count kernels that follow the area law at every gradient order, for the weights of horizon counts. For a rule built from the order alone, Proposition 4 reduces the question to whether e1(θ)=o(θ)e_1(\theta)=o(\theta) on linear data, which is open. A theory whose only objects are the sprinkled elements needs a different argument from Theorem 1, which concerns fields on all pairs of ranks.

What a separate-statistics law does to an actual sprinkling is known only from the computations of section 4, and a proof that typical sprinklings have a trivial automorphism group together with an analytic value for the related fraction would settle the matter for the dominance map. A decoherence functional of rank greater than one needs a history space of dimension at least two; whether a non-commutative sequential-growth model of that kind satisfies the causality conditions is open.

For actions of depth at least one and a non-trivial link phase, the suppression at Kleitman–Rothschild proportions may be of order n2n^2, of lower order, or absent; the uniform bound (15) and the exact data up to n=16n=16 do not decide it. A refinement of (15) restricted to typical P12P_{12} would show whether e−cn2e^{-cn^2} holds for ∣cos⁡(βC0/2)∣\lvert\cos(\beta C_0/2)\rvert in some range, or whether the rate at these proportions equals the link rate. Theorems 9 and 10 cover finite-depth interval actions on three-layer classes with the complete relation Y1≺Y3Y_1\prec Y_3; infinite-depth kernels, actions with an independent nonlocality scale and sums over layer sizes with unrestricted layer assignment need separate treatment. The exact sums at small sizes are the only evidence on whether the rate in (14) is attained for a given action, and the asymptotic growth exponent of Var S\mathrm{Var}\,S for finite-depth actions is determined only over the range of the data in section 7.

References

  1. Sorkin, “Does locality fail at intermediate length-scales?,” arXiv:gr-qc/0703099.
  2. Benincasa, Dowker, “The scalar curvature of a causal set,” Phys. Rev. Lett. 104, 181301 (2010), arXiv:1001.2725.
  3. Aslanbeigi, Saravani, Sorkin, “Generalized causal set d’Alembertians,” JHEP 06 (2014) 024, arXiv:1403.1622.
  4. Johnston, “Feynman propagator for a free scalar field on a causal set,” Phys. Rev. Lett. 103, 180401 (2009), arXiv:0909.0944.
  5. Zapata-Carratalá, Arsiwalla, Beynon, “Heaps of Fish: arrays, generalized associativity and heapoids,” arXiv:2205.05456.
  6. Korzyuk, Rudzko, “Curvilinear parallelogram identity and mean-value property for a semilinear hyperbolic equation of second-order,” Eurasian Math. J. 15 (2024) 61, arXiv:2204.09408.
  7. Kleitman, Rothschild, “Asymptotic enumeration of partial orders on a finite set,” Trans. Amer. Math. Soc. 205 (1975) 205–220.
  8. Loomis, Carlip, “Suppression of non-manifold-like sets in the causal set path integral,” arXiv:1709.00064.
  9. Carlip, Carlip, Surya, “Path integral suppression of badly behaved causal sets,” arXiv:2209.00327.
  10. Carlip, Carlip, Surya, “The Einstein-Hilbert action for entropically dominant causal sets,” arXiv:2311.18238.
  11. Mathur, Singh, Surya, “Entropy and the link action in the causal set path-sum,” arXiv:2009.07623.
  12. Last, Penrose, “Poisson process Fock space representation, chaos expansion and covariance inequalities,” arXiv:0909.3205.
  13. Rideout, Sorkin, “A classical sequential growth dynamics for causal sets,” Phys. Rev. D 61, 024002 (2000), arXiv:gr-qc/9904062.
  14. Dowker, Johnston, Surya, “On extending the quantum measure,” J. Phys. A 43, 505305 (2010), arXiv:1007.2725.
  15. Surya, Zalel, “A criterion for covariance in complex sequential growth models,” arXiv:2003.11311.