Home / Quantum Gravity / Discrete Lorentz-invariant substrates

Lorentz-invariant quantum sprinklings are condensates

Status
Proved: Theorems 2 and 4 (bosonic Weyl algebra, (H1)–(H3)), 6 (CAR algebra, finite density on a ball), 7 (determinantal processes) and 10 (complex processes of bounded variation); Proposition 11 (de Sitter kernels). Computed: void probabilities of random lattices in two dimensions. Heuristic: the lattice void law and count tail. Sketch: the Remark after Corollary 5 (strongly continuous translations assumed).
Excludes
Under (H1)–(H3): bosonic states that are not mixtures of coherent condensates, including sub-Poisson, number-squeezed and Fock-number states of non-zero density. Fermionic CAR states other than the vacuum that are invariant under translations and boosts with finite density on a ball. Non-empty determinantal processes invariant under translations and one boost subgroup.
Revised

A sprinkling of Minkowski space is a random discrete set of points whose law is invariant under the Poincaré group; the standard one is the Poisson process, whose number variance in a region equals its mean. We ask whether a quantum field on spacetime, with the quanta as the points, can be in an invariant state of finite density and count more regularly than Poisson, with a number variance below the mean or a fixed number of quanta in a region. For a bosonic field with wave functions on spacetime the answer is no, if the state is locally normal (each bounded region has a density matrix) and the mean number of quanta in every bounded region is finite (section 1.4). Every such state is a mixture of coherent states, each with one constant complex wave function zz. Its counts are Poisson with random intensity ∣z∣2\lvert z\rvert ^2, the same at every point, so the variance is at least the mean (Theorem 4).

Poisson sprinklings are the standard model of a causal set [Surya], [Sorkin]. Bombelli, Henson and Sorkin [BHS] prove that no Lorentz-equivariant choice of a preferred direction at a point of a Poisson sprinkling exists. Dowker, Henson and Sorkin [DHS] explain why a boosted regular lattice is not a sprinkling with the same law as the original, and Saravani and Aslanbeigi [SA] show that Poisson sprinkling gives the best number-volume correspondence in arbitrarily small regions, while two-dimensional Lorentzian lattices do better in large ones. Hudson and Parthasarathy [HP] realise Poisson processes inside Fock space. Theorem 4 classifies all invariant, locally normal states of finite density in that setting, including the phase of the condensate.

The same boost argument leaves only the vacuum for fermions (Theorem 6) and only the empty process for determinantal point processes, the classical model of a fermionic sprinkling (Theorem 7). The classical analogue of the bosonic statement fails: Haar-random lattices are invariant, ergodic and not Poisson (Proposition 9), and invariant point processes with complex weights of bounded total variation are combinations of classical ergodic processes (Theorem 10). The counts of a condensate are covered by Theorem 1 of Integer counts and finite-dimensional algebras in a positive-energy vacuum: such a state is not the vacuum of a positive-energy translation group, with the identification of the counts with the hypotheses of that theorem a sketch (Remark after Corollary 5). The consequences of a random intensity for the causal-set d’Alembertian are worked out in Mixed-Poisson and signed substrates and the four-dimensional d’Alembertian.

1. Sprinklings, counts and the setting

1.1 Classical counts

Let M=R1,3M=\mathbb R^{1,3} with Lebesgue measure dVdV, let V(D)V(D) be the volume of a bounded Borel set DD, and let PP be the proper orthochronous Poincaré group. The same arguments work in R1,d−1\mathbb R^{1,d-1} for every d≥2d\ge 2. The signature is (−,+,+,+)(-,+,+,+), so the Minkowski pairing of a momentum pp and a point xx is p⋅x=−p0x0+p⃗⋅x⃗p\cdot x=-p^0x^0+\vec p\cdot\vec x. A sprinkling is a random locally finite subset of MM, and a causal set is obtained by keeping the points and the causal order between them. Everything below depends only on the counts N(D)N(D), the number of points in DD. The Poisson process of intensity λ\lambda has N(D)N(D) Poisson with mean λV(D)\lambda V(D) and independent counts in disjoint sets. Its mean and variance are both λV(D)\lambda V(D), and its law is invariant under PP because Lebesgue measure is.

1.2 The quantum version

Take the one-particle space to be L2(M,dV)L^2(M,dV): a one-particle state is a wave function on spacetime, and the points of the sprinkling are the positions of the quanta. Let Γ(h)\Gamma(h) be the bosonic Fock space over a Hilbert space hh, with annihilation operators a(f)a(f) (antilinear in ff), creation operators a∗(f)a^*(f) and [a(f),a∗(g)]=⟨f,g⟩[a(f),a^*(g)]=\langle f,g\rangle, where the inner product is antilinear in the first slot. The Weyl operators W(f)=exp⁡(a∗(f)−a(f))W(f)=\exp(a^*(f)-a(f)) satisfy W(f)W(g)=e−i Im⟨f,g⟩W(f+g)W(f)W(g)=e^{-i\,\mathrm{Im}\langle f,g\rangle}W(f+g). For a bounded Borel set DD of positive volume let L2(D)⊂L2(M)L^2(D)\subset L^2(M) be the functions vanishing outside DD, let A(D)\mathfrak A(D) be the C∗^*-algebra generated by W(f)W(f) with f∈L2(D)f\in L^2(D), and let A\mathfrak A be the algebra generated by all of them. The group PP acts by automorphisms αgW(f)=W(f∘g−1)\alpha_gW(f)=W(f\circ g^{-1}). The count of DD is the number operator of Γ(L2(D))\Gamma(L^2(D)),

N(D)=dΓ(1D)=∑ka∗(ek)a(ek),(1)N(D)=d\Gamma(1_D)=\sum_k a^*(e_k)a(e_k), \tag{1}

for any orthonormal basis (ek)(e_k) of L2(D)L^2(D). It has spectrum {0,1,2,… }\{0,1,2,\dots\}, and for disjoint D1,D2D_1,D_2 the Fock space of D1∪D2D_1\cup D_2 is the tensor product of those of D1D_1 and D2D_2, so counts in disjoint sets commute. This is a sharp count: it assigns an integer to every region and is not smeared against a profile.

1.3 Coherent states are Poisson processes

For z∈Cz\in\mathbb C let ωz\omega_z be the state on A(D)\mathfrak A(D) given by the normalised exponential vector e−∣z∣2V(D)/2e(z1D)e^{-\lvert z\rvert ^2V(D)/2}e(z1_D) of Γ(L2(D))\Gamma(L^2(D)). Then

ωz(W(f))=exp⁡ ⁣(−12∥f∥2+2i Im(zˉ ⁣∫ ⁣f)),f∈L2(D),(2)\omega_z(W(f))=\exp\!\big(-\tfrac12\lVert f\rVert ^2+2i\,\mathrm{Im}(\bar z\!\int\! f)\big),\qquad f\in L^2(D), \tag{2}

which is consistent in DD and extends to every f∈L1∩L2f\in L^1\cap L^2. The operator eitN(D)e^{itN(D)} maps the exponential vector e(g)e(g) to e(eitg)e(e^{it}g), so ωz(eitN(D))=exp⁡(∣z∣2V(D)(eit−1))\omega_z(e^{itN(D)})=\exp\big(\lvert z\rvert ^2V(D)(e^{it}-1)\big), the generating function of a Poisson variable of mean ∣z∣2V(D)\lvert z\rvert ^2V(D). The exponential vector also factorises over a partition of DD, with a coherent state in each factor, so counts of disjoint sets are independent. Hence ωz\omega_z is a Poisson sprinkling of intensity ∣z∣2\lvert z\rvert ^2. It is Poincaré invariant, since ∥f∘g−1∥=∥f∥\lVert f\circ g^{-1}\rVert =\lVert f\rVert and ∫f∘g−1=∫f\int f\circ g^{-1}=\int f. It is not a vector in the Fock space of L2(M)L^2(M) when z≠0z\ne0, because the constant z 1Mz\,1_M is not square integrable. For this reason we work with the local algebras A(D)\mathfrak A(D) and not with the Weyl algebra over L2(M)L^2(M), on which ∫f\int f is undefined. Working over L1∩L2L^1\cap L^2 is possible, but then the state on functions of unbounded support is an additional datum that the counts do not see.

1.4 Hypotheses

A state ω\omega on A\mathfrak A is a quantum sprinkling if it satisfies:

(H1) Invariance. ω∘αg=ω\omega\circ\alpha_g=\omega for all g∈Pg\in P.

(H2) Local normality. For every bounded DD there is a density matrix ϱD\varrho_D on Γ(L2(D))\Gamma(L^2(D)) with ω(A)=tr(ϱDA)\omega(A)=\mathrm{tr}(\varrho_D A) for A∈A(D)A\in\mathfrak A(D).

(H3) Finite density. tr(ϱDN(D))<∞\mathrm{tr}(\varrho_D N(D))<\infty for every bounded DD. We write ω(N(D))\omega(N(D)) for this number.

Local normality makes t↦ω(W(tf))t\mapsto\omega(W(tf)) continuous, so the field operators a(f)a(f) exist. The hypothesis (H3) is used in Theorem 2 to make a certain measure finite; section 8 discusses what happens without it. The one-body density is the operator TT on L2(M)L^2(M) with

⟨f,Tg⟩=ω(a∗(g) a(f)),f,g∈Lc2(M),(3)\langle f,Tg\rangle=\omega\big(a^*(g)\,a(f)\big),\qquad f,g\in L^2_c(M), \tag{3}

where Lc2L^2_c is the set of square-integrable functions with bounded support. The right side is finite: (H3) puts the range of ϱD1/2\varrho_D^{1/2} in the domain of N(D)1/2N(D)^{1/2}, and ∥a(f)ψ∥≤∥f∥ ∥N(D)1/2ψ∥\lVert a(f)\psi\rVert \le\lVert f\rVert \,\lVert N(D)^{1/2}\psi\rVert. For a coherent state, (3) reads ⟨f,Tg⟩=∣z∣2 ∫f‾∫g\langle f,Tg\rangle=\lvert z\rvert ^2\,\overline{\textstyle\int f}\int g, a rank-one operator with constant kernel. Theorem 2 shows that this is the only possibility.

2. The one-body density is a constant

Lemma 1 (wandering sets). Let μ\mu be a finite positive Borel measure on Rn\mathbb R^n, and let BB be a linear map of determinant ±1\pm1 with B∗μ=μB_*\mu=\mu. If YY is a Borel set such that the sets BkYB^kY, k∈Zk\in\mathbb Z, are pairwise disjoint, then μ(Y)=0\mu(Y)=0. The same holds for a function G∈L1(Rn)G\in L^1(\mathbb R^n) with G∘B=GG\circ B=G almost everywhere, with ∫Y∣G∣=0\int_Y\lvert G\rvert =0 in place of μ(Y)=0\mu(Y)=0, and then G=0G=0 almost everywhere on ⋃kBkY\bigcup_kB^kY.

Proof. The sets BkYB^kY all have the same measure, and they are disjoint, so ∑kμ(BkY)≤μ(Rn)<∞\sum_k\mu(B^kY)\le\mu(\mathbb R^n)<\infty forces μ(Y)=0\mu(Y)=0. For GG the same argument applies to ∫BkY∣G∣=∫Y∣G∘Bk∣=∫Y∣G∣\int_{B^kY}\lvert G\rvert =\int_Y\lvert G\circ B^k\rvert =\int_Y\lvert G\rvert, which uses that BB preserves Lebesgue measure. □\square

A boost with rapidity θ\theta along x1x^1 acts on momenta in the null coordinates p±=p0±p1p^\pm=p^0\pm p^1 as (p+,p−)↦(eθp+,e−θp−)(p^+,p^-)\mapsto(e^\theta p^+,e^{-\theta}p^-) and fixes p2,p3p^2,p^3. The sets Y={1≤∣p+∣<eθ}Y=\{1\le\lvert p^+\rvert <e^{\theta}\} and Y′={1≤∣p−∣<eθ}Y'=\{1\le\lvert p^-\rvert <e^\theta\} are wandering, and their images cover {p+≠0}\{p^+\ne0\} and {p−≠0}\{p^-\ne0\}. A finite boost-invariant measure is therefore carried by {p+=p−=0}\{p^+=p^-=0\}, and this is what makes the one-body density constant. The same wandering-set argument, for a group containing a hyperbolic element, is Lemma 1 of Lorentz-invariant point sets: what they can and cannot carry.

Theorem 2 (constant one-body density). Under (H1)–(H3) there is a number ρ≥0\rho\ge0 with

ω(a∗(g)a(f))=ρ ∫f‾ ∫gfor all f,g∈Lc2(M),(4)\omega\big(a^*(g)a(f)\big)=\rho\,\overline{\textstyle\int f}\,\textstyle\int g\qquad\text{for all } f,g\in L^2_c(M), \tag{4}

and ω(N(D))=ρ V(D)\omega(N(D))=\rho\,V(D) for every bounded DD.

Translation invariance writes the density as a positive measure μ\mu on momentum space, finite density makes the total mass of μ\mu finite, and a finite boost-invariant measure sits at the origin.

Proof. Step 1: a spectral measure. The form τ(f,g)=⟨f,Tg⟩\tau(f,g)=\langle f,Tg\rangle is positive and translation invariant. For f,g∈L2(D)f,g\in L^2(D) it is bounded: a∗(f)a(f)≤∥f∥2N(D)a^*(f)a(f)\le\lVert f\rVert ^2N(D), because a∗(e)a(e)a^*(e)a(e) with ∥e∥=1\lVert e\rVert =1 is one term of (1) in a basis containing ee, so Cauchy–Schwarz gives ∣τ(f,g)∣≤∥f∥∥g∥ ω(N(D))\lvert \tau(f,g)\rvert \le\lVert f\rVert \lVert g\rVert \,\omega(N(D)). By the Schwartz kernel theorem and the Bochner–Schwartz theorem [RS2] there is a unique positive tempered measure μ\mu on R4\mathbb R^4 with

τ(f,g)=∫f^(p)‾ g^(p) dμ(p),f^(p)=∫e−ip⋅xf(x) dx,f,g∈Cc∞(M),(5)\tau(f,g)=\int\overline{\hat f(p)}\,\hat g(p)\,d\mu(p),\qquad \hat f(p)=\int e^{-ip\cdot x}f(x)\,dx,\quad f,g\in C_c^\infty(M), \tag{5}

where p⋅xp\cdot x is the Minkowski pairing. The products f^‾g^\overline{\hat f}\hat g are the Fourier transforms of the convolutions f~∗g∈Cc∞\tilde f*g\in C_c^\infty, with f~(x)=f(−x)‾\tilde f(x)=\overline{f(-x)}, and these are dense, so τ\tau determines μ\mu.

Step 2: finite density is finite mass. Let DD be an open ball and let (ek)(e_k) be an orthonormal basis of L2(D)L^2(D) made of functions in Cc∞(D)C_c^\infty(D). By (1), (3), (5) and Tonelli,

ω(N(D))=∑kτ(ek,ek)=∫∑k∣e^k(p)∣2 dμ(p)=V(D) μ(R4).\omega(N(D))=\sum_k\tau(e_k,e_k)=\int\sum_k\lvert \hat e_k(p)\rvert ^2\,d\mu(p)=V(D)\,\mu(\mathbb R^4).

The last step is Parseval: e^k(p)=⟨eip⋅1D,ek⟩\hat e_k(p)=\langle e^{ip\cdot}1_D,e_k\rangle up to conjugation, so ∑k∣e^k(p)∣2=∥eip⋅1D∥2=V(D)\sum_k\lvert \hat e_k(p)\rvert ^2=\lVert e^{ip\cdot}1_D\rVert ^2=V(D) for every pp. Hence (H3) gives μ(R4)<∞\mu(\mathbb R^4)<\infty.

Step 3: boosts. A boost BB acts on test functions by f↦f∘B−1f\mapsto f\circ B^{-1} and on transforms by f^↦f^∘B−1\hat f\mapsto\hat f\circ B^{-1}, since the pairing is invariant. Invariance (H1) of τ\tau therefore says that μ\mu is invariant under BB. Lemma 1 for the boost along x1x^1 with the wandering sets above gives μ(p+≠0)=μ(p−≠0)=0\mu(p^+\ne0)=\mu(p^-\ne0)=0, so μ\mu lives on {p0=p1=0}\{p^0=p^1=0\}. The boosts along x2x^2 and x3x^3 give {p0=p2=0}\{p^0=p^2=0\} and {p0=p3=0}\{p^0=p^3=0\}, and the intersection is the origin. Thus μ=ρ δ0\mu=\rho\,\delta_0 with ρ=μ(R4)\rho=\mu(\mathbb R^4), and (5) becomes (4) because f^(0)=∫f\hat f(0)=\int f. Since μ\mu is finite, dominated convergence extends (5) from Cc∞C_c^\infty to f,g∈Lc2f,g\in L^2_c (approximate in L2L^2 and in L1L^1 simultaneously, so that f^n→f^\hat f_n\to\hat f boundedly), which gives (4) as stated. The Parseval computation of Step 2 with an orthonormal basis of L2(D)L^2(D) then gives ω(N(D))=ρV(D)\omega(N(D))=\rho V(D) for every bounded DD. □\square

The proof uses translations and boosts only, not rotations and not positivity of the energy. With translations alone any finite measure μ\mu is possible (section 8).

3. Reduction to a single mode

Fix a bounded set DD of positive volume and let uD=1D/V(D)u_D=1_D/\sqrt{V(D)} be the normalised constant on DD. By (4), a(f)a(f) has zero second moment in ω\omega whenever ∫f=0\int f=0. Proposition 3 turns this into a statement about the state: all quanta of DD are in the single mode uDu_D, and the single-mode states of nested regions are related by loss. The relation (6) is what lets us pass from small regions to large ones in section 4.

Proposition 3 (single-mode reduction and attenuation). Assume (H1)–(H3). Then:

(a) For every bounded DD of positive volume, Γ(L2(D))=Γ(CuD)⊗Γ(uD⊥)\Gamma(L^2(D))=\Gamma(\mathbb Cu_D)\otimes\Gamma(u_D^\perp) and ϱD=σD⊗∣Ω⟩⟨Ω∣\varrho_D=\sigma_D\otimes\lvert \Omega\rangle\langle\Omega\rvert, where σD\sigma_D is a density matrix on the one-mode space Γ(CuD)≅ℓ2(N0)\Gamma(\mathbb Cu_D)\cong\ell^2(\mathbb N_0) and Ω\Omega is the vacuum of Γ(uD⊥)\Gamma(u_D^\perp). In particular N(D)=a∗(uD)a(uD)N(D)=a^*(u_D)a(u_D) on the support of ϱD\varrho_D.

(b) If D⊂D′D\subset D' and η=V(D)/V(D′)\eta=V(D)/V(D'), then σD\sigma_D is the image of σD′\sigma_{D'} under the pure-loss channel of transmissivity η\eta, equivalently

ΦD(c)=ΦD′(η c),ΦD(c):=tr(σD eca†e−cˉa),(6)\Phi_D(c)=\Phi_{D'}(\sqrt\eta\,c),\qquad \Phi_D(c):=\mathrm{tr}\big(\sigma_D\,e^{c a^\dagger}e^{-\bar c a}\big), \tag{6}

where a=a(uD)a=a(u_D) on the left and a=a(uD′)a=a(u_{D'}) on the right.

In (a) a mode orthogonal to the constant carries no quanta. In (b) the constant on DD is a mixture of the constant on D′D', with weight η\sqrt\eta, and a function orthogonal to the constant on D′D', which carries no quanta, so the mode uDu_D sees a fraction η\eta of the quanta of the mode uD′u_{D'}.

Proof. (a) Let f∈L2(D)f\in L^2(D) with ∫f=0\int f=0, that is, f⊥uDf\perp u_D. By (4), tr(ϱDa∗(f)a(f))=0\mathrm{tr}(\varrho_Da^*(f)a(f))=0, which is ∥a(f)ϱD1/2∥HS2=0\lVert a(f)\varrho_D^{1/2}\rVert ^2_{HS}=0, so the support of ϱD\varrho_D lies in ker⁡a(f)\ker a(f). The common kernel of a(f)a(f) over f∈uD⊥f\in u_D^\perp is Γ(CuD)⊗Ω\Gamma(\mathbb Cu_D)\otimes\Omega, the only vectors of Γ(uD⊥)\Gamma(u_D^\perp) annihilated by all a(f)a(f) being multiples of the vacuum.

(b) Write uD=η uD′+wu_D=\sqrt\eta\,u_{D'}+w with w=uD−η uD′∈L2(D′)w=u_D-\sqrt\eta\,u_{D'}\in L^2(D'). Then ⟨uD′,w⟩=0\langle u_{D'},w\rangle=0, ∫w=0\int w=0 and ∥w∥2=1−η\lVert w\rVert ^2=1-\eta. For c∈Cc\in\mathbb C the Weyl operator W(cuD)=W(cη uD′)W(cw)W(cu_D)=W(c\sqrt\eta\,u_{D'})W(cw) is a product of commuting factors, one in the mode uD′u_{D'} and one in uD′⊥u_{D'}^\perp. By (a) for D′D',

ω(W(cuD))=tr(σD′W(cη uD′)) e−(1−η)∣c∣2/2.\omega(W(cu_D))=\mathrm{tr}\big(\sigma_{D'}W(c\sqrt\eta\,u_{D'})\big)\,e^{-(1-\eta)\lvert c\rvert ^2/2}.

By (a) for DD the left side is tr(σDW(cuD))\mathrm{tr}(\sigma_DW(cu_D)). Passing to the normally ordered function with W(cu)=eca†e−cˉae−∣c∣2/2W(cu)=e^{ca^\dagger}e^{-\bar ca}e^{-\lvert c\rvert ^2/2} turns this into (6). □\square

The relation (6) holds for every pair of nested regions, so the state of a mode in a small region is obtained from the state in a large region by loss, for all loss rates at once. Normally ordered moments scale as ⟨a†jak⟩D=η(j+k)/2⟨a†jak⟩D′\langle a^{\dagger j}a^k\rangle_D=\eta^{(j+k)/2}\langle a^{\dagger j}a^k\rangle_{D'} (with the appropriate mode on each side). Mandel’s parameter QM=(⟨a†2a2⟩−⟨a†a⟩2)/⟨a†a⟩=(Var N−⟨N⟩)/⟨N⟩Q_{\rm M}=(\langle a^{\dagger2}a^2\rangle-\langle a^\dagger a\rangle^2)/\langle a^\dagger a\rangle=(\mathrm{Var}\,N-\langle N\rangle)/\langle N\rangle therefore obeys QM(D)=η QM(D′)Q_{\rm M}(D)=\eta\,Q_{\rm M}(D'). Every state has QM≥−1Q_{\rm M}\ge-1, so QM(D)≥−V(D)/V(D′)Q_{\rm M}(D)\ge-V(D)/V(D'), and letting D′D' grow gives QM(D)≥0Q_{\rm M}(D)\ge0. Whenever the second moments are finite, Var N(D)≥⟨N(D)⟩\mathrm{Var}\,N(D)\ge\langle N(D)\rangle follows already at this point.

Example (a number state cannot occur). Suppose the mode of a region has exactly one quantum, σD=∣1⟩⟨1∣\sigma_D=\lvert 1\rangle\langle1\rvert. By (6) with η=1/k\eta=1/k, k>1k>1, the state σD′\sigma_{D'} of the region kk times larger has ΦD′(c)=1−k∣c∣2\Phi_{D'}(c)=1-k\lvert c\rvert ^2. The normally ordered function of the number state ∣n⟩\lvert n\rangle is the Laguerre polynomial Ln(∣c∣2)L_n(\lvert c\rvert ^2), so 1−k∣c∣2=p0L0+p1L11-k\lvert c\rvert ^2=p_0L_0+p_1L_1 forces p1=k>1p_1=k>1 and p0=1−k<0p_0=1-k<0, which is not a state. For k=3/2k=3/2 the weights are (p0,p1)=(−12,32)(p_0,p_1)=(-\tfrac12,\tfrac32), and the pure-loss channel with η=2/3\eta=2/3 maps this operator back to ∣1⟩⟨1∣\lvert 1\rangle\langle 1\rvert. The obstruction is the positivity of the preimage.

4. Classification

4.1 Bosons

Theorem 4 (invariant quantum sprinklings are coherent condensates). Let ω\omega satisfy (H1)–(H3). There is a unique Borel probability measure m\mathsf m on C\mathbb C with ∫∣z∣2dm(z)=ρ<∞\int\lvert z\rvert ^2d\mathsf m(z)=\rho<\infty such that

ω(A)=∫ωz(A) dm(z),A∈A(D), D bounded.(7)\omega(A)=\int\omega_z(A)\,d\mathsf m(z),\qquad A\in\mathfrak A(D),\ D\ \text{bounded}. \tag{7}

Conversely every such mixture satisfies (H1)–(H3). Consequently, for disjoint bounded sets D1,…,DkD_1,\dots,D_k the counts have the joint law

Pr⁡(N(Di)=ni, i≤k)=E∏ie−ΛV(Di)(ΛV(Di))nini!,Λ=∣z∣2,(8)\Pr\big(N(D_i)=n_i,\ i\le k\big)=\mathbb E\prod_{i}e^{-\Lambda V(D_i)}\frac{(\Lambda V(D_i))^{n_i}}{n_i!},\qquad \Lambda=\lvert z\rvert ^2, \tag{8}

so NN is a Cox (mixed Poisson) process with a spacetime-constant random intensity Λ\Lambda of mean ρ\rho, and Var N(D)=ρV(D)+V(D)2 Var Λ ≥ EN(D)\mathrm{Var}\,N(D)=\rho V(D)+V(D)^2\,\mathrm{Var}\,\Lambda\ \ge\ \mathbb E N(D), with equality for all DD if and only if Λ=ρ\Lambda=\rho almost surely.

By (6) the one-mode state of a small region is obtained from the state of every larger region by losing each quantum independently with probability 1−η1-\eta (a pure-loss channel). The Husimi function of a one-mode state is the probability density of the outcomes of a measurement in coherent states, so its Fourier transform is positive definite; it equals the normally ordered characteristic function ΦD\Phi_D of (6) times a Gaussian whose width grows with the volume. In the limit of large regions the Gaussian drops out, and ΦD\Phi_D becomes the Fourier transform of a positive measure on coherent amplitudes. Passing to this limit V(D)→∞V(D)\to\infty uses that (6) relates every region to arbitrarily larger ones.

Proof. Let Ψ(ζ)=ΦD(ζ/V(D))\Psi(\zeta)=\Phi_D(\zeta/\sqrt{V(D)}). By (6) with c=ζ/V(D)c=\zeta/\sqrt{V(D)} and η c=ζ/V(D′)\sqrt\eta\,c=\zeta/\sqrt{V(D')}, Ψ\Psi is the same for nested regions, hence for all bounded regions (any two lie in a third). It is continuous with Ψ(0)=1\Psi(0)=1, since σD\sigma_D is a normal state.

The Husimi function qD(α)=π−1⟨α∣σD∣α⟩q_D(\alpha)=\pi^{-1}\langle\alpha\lvert \sigma_D\rvert \alpha\rangle is a probability density on C\mathbb C. Its Fourier transform is the anti-normally ordered function tr(σDe−ξˉaeξa†)=ΦD(ξ)e−∣ξ∣2\mathrm{tr}(\sigma_De^{-\bar\xi a}e^{\xi a^\dagger})=\Phi_D(\xi)e^{-\lvert \xi\rvert ^2}. A Fourier transform of a probability density is positive definite, so with ζ=V(D) ξ\zeta=\sqrt{V(D)}\,\xi the function Ψ(ζ)e−∣ζ∣2/V(D)\Psi(\zeta)e^{-\lvert \zeta\rvert ^2/V(D)} is positive definite on R2≅C\mathbb R^2\cong\mathbb C for every V(D)∈(0,∞)V(D)\in(0,\infty). Let V(D)→∞V(D)\to\infty. A pointwise limit of positive-definite functions is positive definite, so Ψ\Psi is a continuous positive-definite function with Ψ(0)=1\Psi(0)=1, and Bochner’s theorem gives a probability measure m\mathsf m on C\mathbb C with

Ψ(ζ)=∫ezˉζ−zζˉ dm(z).\Psi(\zeta)=\int e^{\bar z\zeta-z\bar\zeta}\,d\mathsf m(z).

A coherent state ∣β⟩\lvert \beta\rangle has tr(∣β⟩⟨β∣eξa†e−ξˉa)=eξβˉ−ξˉβ\mathrm{tr}(\lvert \beta\rangle\langle\beta\rvert e^{\xi a^\dagger}e^{-\bar\xi a})=e^{\xi\bar\beta-\bar\xi\beta}, and the normally ordered function determines the state, so σD=∫∣zV(D)⟩⟨zV(D)∣ dm(z)\sigma_D=\int\lvert z\sqrt{V(D)}\rangle\langle z\sqrt{V(D)}\rvert \,d\mathsf m(z). By (2) and Proposition 3(a), ∣zV(D)⟩⊗Ω\lvert z\sqrt{V(D)}\rangle\otimes\Omega is the exponential vector defining ωz\omega_z on A(D)\mathfrak A(D), which gives (7). Uniqueness of m\mathsf m follows from that of the Fourier transform. Finite density gives ∫∣z∣2dm=ρ\int\lvert z\rvert ^2d\mathsf m=\rho.

For the converse, (2) shows that each ωz\omega_z, hence the mixture, is invariant; local normality and finite density are clear. The law (8) follows from the independence and Poisson counts under each ωz\omega_z. The mean and variance of a mixed Poisson variable with random mean ΛV\Lambda V are EΛV\mathbb E\Lambda V and EΛV+Var(Λ)V2\mathbb E\Lambda V+\mathrm{Var}(\Lambda)V^2, which is +∞+\infty if EΛ2=∞\mathbb E\Lambda^2=\infty. □\square

The counts determine only the law of Λ=∣z∣2\Lambda=\lvert z\rvert ^2, whereas the state also contains the phase of zz. For m=δz0\mathsf m=\delta_{z_0} the state is a single coherent condensate. For m\mathsf m uniform on the circle ∣z∣2=ρ\lvert z\rvert ^2=\rho the counts are exactly Poisson and the state is the phase average of a condensate, whose counting statistics coincide with those of the classical Poisson sprinkling.

Corollary 5 (what cannot occur). Let ω\omega satisfy (H1)–(H3) with ρ>0\rho>0. Then no region has Var N(D)<EN(D)\mathrm{Var}\,N(D)<\mathbb E N(D) (sub-Poisson or number-squeezed counts). No region of positive volume has a fixed number of quanta (a Fock number state of the mode uDu_D). The variance does not follow an area law, Var N(D)∝\mathrm{Var}\,N(D)\propto the area of ∂D\partial D. The only fluctuation beyond Poisson is a random but spacetime-constant intensity, which contributes a term V2 Var(Λ)V^2\,\mathrm{Var}(\Lambda) growing like the square of the volume.

The case ρ=0\rho=0 is the vacuum, in which N(D)=0N(D)=0 with certainty. The variance is ρV+V2VarΛ\rho V+V^2\mathrm{Var}\Lambda for every region, with no term that depends on the boundary, which excludes an area law.

Remark (positive energy). Let ω\omega be as in Theorem 4 with EΛ2<∞\mathbb E\Lambda^2<\infty. By local normality the counts (1) act in the GNS space of ω\omega (the Hilbert space built from the state, in which ω\omega is the expectation in a cyclic vector) as commuting self-adjoint operators with integer spectrum, additive on disjoint sets and covariant under the unitary translations that implement (H1); the second moment of every count in the GNS vector is finite by the variance formula of Theorem 4. If the translations are strongly continuous in the GNS representation, these data satisfy the hypotheses (C1)–(C4) and (R) of Integer counts and finite-dimensional algebras in a positive-energy vacuum, and its Theorem 1 gives ρ=0\rho=0 whenever the spectrum of the translation generator lies in the closed forward cone. Hence a condensate with ρ>0\rho>0 is not the vacuum of a positive-energy translation group.

4.2 Fermions

Replace the Weyl algebra by the CAR algebra generated by a(f)a(f), f∈Lc2(M)f\in L^2_c(M), with {a(f),a∗(g)}=⟨f,g⟩\{a(f),a^*(g)\}=\langle f,g\rangle. Every state of this algebra satisfies ∣ω(a∗(g)a(f))∣≤∥f∥∥g∥\lvert \omega(a^*(g)a(f))\rvert \le\lVert f\rVert \lVert g\rVert, because ∥a(f)∥=∥f∥\lVert a(f)\rVert =\lVert f\rVert, so TT is a bounded operator with 0≤T≤10\le T\le1, and no normality or regularity is needed. For a state we define ω(N(D))=∑kω(a∗(ek)a(ek))∈[0,∞]\omega(N(D))=\sum_k\omega(a^*(e_k)a(e_k))\in[0,\infty] for an orthonormal basis of L2(D)L^2(D), which does not depend on the basis.

Theorem 6 (no invariant fermionic sprinkling). Let ω\omega be a state on the CAR algebra over Lc2(M)L^2_c(M) that is invariant under translations and boosts and has ω(N(D))<∞\omega(N(D))<\infty for some open ball DD. Then ω\omega is the Fock vacuum: ω(N(D))=0\omega(N(D))=0 for every DD and ω(a∗(f)a(g))=0\omega(a^*(f)a(g))=0.

For fermions the bound T≤1T\le1 rules out the constant density that the boost argument allows: on a region of volume V(D′)V(D') the constant density gives ρV(D′)2\rho V(D')^2 for the smeared number of the function 1D′1_{D'}, while T≤1T\le1 allows at most V(D′)V(D'). Poincaré-invariant quasi-free CAR states with 0≤T≤10\le T\le1 whose single-mode numbers are sub-Poisson do exist, namely those whose TT is a bounded invariant Fourier multiplier, but they have infinite density (section 8).

Proof. Theorem 2 uses nothing about the statistics except the bound on τ\tau, which here comes from ∥a(f)∥=∥f∥\lVert a(f)\rVert =\lVert f\rVert. Steps 1 to 3 give τ(f,g)=ρ ∫f‾∫g\tau(f,g)=\rho\,\overline{\int f}\int g, with finite ρ\rho since ω(N(D))=V(D)μ(R4)<∞\omega(N(D))=V(D)\mu(\mathbb R^4)<\infty (Step 2 applies to the one ball DD for which the density is finite, and the extension of (4) from Cc∞C_c^\infty to Lc2L^2_c is by continuity, as in the proof of Theorem 2). Now take f=1D′f=1_{D'} for any ball D′D': τ(f,f)=ω(a∗(f)a(f))≤∥f∥2\tau(f,f)=\omega(a^*(f)a(f))\le\lVert f\rVert ^2 gives ρV(D′)2≤V(D′)\rho V(D')^2\le V(D'), that is ρ≤1/V(D′)\rho\le1/V(D'). Since D′D' can have arbitrarily large volume, ρ=0\rho=0. Then ω(a∗(f)a(f))=0\omega(a^*(f)a(f))=0 for all f∈Lc2f\in L^2_c, and ω(N(D))=0\omega(N(D))=0 for every bounded DD. For any BB in the algebra, Cauchy–Schwarz gives ∣ω(Ba(f))∣2≤ω(BB∗) ω(a∗(f)a(f))=0\lvert \omega(Ba(f))\rvert ^2\le\omega(BB^*)\,\omega(a^*(f)a(f))=0, and by taking adjoints ω(a∗(f)B)=0\omega(a^*(f)B)=0. Hence ω\omega vanishes on every normally ordered monomial a∗(f1)⋯a∗(fm)a(g1)⋯a(gn)a^*(f_1)\cdots a^*(f_m)a(g_1)\cdots a(g_n) with m+n≥1m+n\ge1. These monomials and the identity span the algebra, so ω\omega is the vacuum state. □\square

5. Classical determinantal sprinklings

The point process of a quasi-free fermionic state is determinantal, and the classical statement is simpler and needs no Fock space. A point process on MM is determinantal with Hermitian kernel K\mathcal K if K\mathcal K is a self-adjoint operator on L2(M)L^2(M) with 0≤K≤10\le\mathcal K\le1 and 1DK1D1_D\mathcal K1_D trace class for bounded DD, and the process has correlation functions ρn(x1,…,xn)=det⁡[K(xi,xj)]\rho_n(x_1,\dots,x_n)=\det[K(x_i,x_j)], where K(x,y)K(x,y) is a measurable kernel of K\mathcal K with ∫DK(x,x)dx=tr(1DK1D)\int_DK(x,x)dx=\mathrm{tr}(1_D\mathcal K1_D) [Soshnikov]. In particular ρ1(x)=K(x,x)\rho_1(x)=K(x,x) and ρ2(x,y)=ρ1(x)ρ1(y)−∣K(x,y)∣2\rho_2(x,y)=\rho_1(x)\rho_1(y)-\lvert K(x,y)\rvert ^2.

Theorem 7 (no invariant determinantal process). Let Ξ\Xi be a determinantal point process on R1,d−1\mathbb R^{1,d-1}, d≥2d\ge2, whose law is invariant under translations and under one boost subgroup. Then Ξ\Xi is empty almost surely.

The point is that invariance makes ∣K(x,y)∣2\lvert K(x,y)\rvert ^2 a boost-invariant function of x−yx-y, and K2≤K\mathcal K^2\le\mathcal K makes it integrable, so Lemma 1 kills it. Integrability is obtained from the Hilbert–Schmidt norm of K1D\mathcal K1_D, which replaces the pointwise inequality (K2)(x,x)≤K(x,x)(\mathcal K^2)(x,x)\le K(x,x) and so avoids a choice of the diagonal values of the kernel.

Proof. Invariance of the law makes ρ1\rho_1 invariant under all translations almost everywhere, hence equal to a constant ρ\rho almost everywhere (by Fubini, a locally integrable function that is invariant a.e. under every translation is a.e. constant). Likewise ρ2(x+a,y+a)=ρ2(x,y)\rho_2(x+a,y+a)=\rho_2(x,y) gives ∣K(x,y)∣2=G(x−y)\lvert K(x,y)\rvert ^2=G(x-y) for a measurable G≥0G\ge0, and invariance under a boost BB gives G∘B=GG\circ B=G almost everywhere.

Next, K2≤K\mathcal K^2\le\mathcal K. For bounded DD the operator K1D\mathcal K1_D is Hilbert–Schmidt, with ∥K1D∥HS2=tr(1DK21D)≤tr(1DK1D)=ρV(D)\lVert \mathcal K1_D\rVert _{HS}^2=\mathrm{tr}(1_D\mathcal K^2 1_D)\le\mathrm{tr}(1_D\mathcal K1_D)=\rho V(D). The left side is ∫ ⁣ ⁣∫∣K(x,y)∣21D(y) dx dy=V(D)∫G\int\!\!\int\lvert K(x,y)\rvert ^21_D(y)\,dx\,dy=V(D)\int G by the form of ∣K∣2\lvert K\rvert ^2. Hence G∈L1(Rd)G\in L^1(\mathbb R^d) with ∫G≤ρ\int G\le\rho. Lemma 1, applied to the wandering sets of the boost BB (the regions ekθ≤∣r0+r1∣<e(k+1)θe^{k\theta}\le\lvert r^0+r^1\rvert <e^{(k+1)\theta} in the difference variable r=x−yr=x-y), gives G=0G=0 almost everywhere. Then K=0K=0 almost everywhere, so K1D\mathcal K1_D has zero Hilbert–Schmidt norm for every DD and K=0\mathcal K=0. Finally ρV(D)=tr(1DK1D)=0\rho V(D)=\mathrm{tr}(1_D\mathcal K1_D)=0, so Ξ=∅\Xi=\emptyset. □\square

The same lemma gives a statement about hyperuniformity, the property that the number variance grows more slowly than the volume, Var N(D)/V(D)→0\mathrm{Var}\,N(D)/V(D)\to0 along a growing sequence of regions.

Corollary 8 (pair correlations in L1L^1). Let Ξ\Xi be a simple point process on R1,d−1\mathbb R^{1,d-1} with locally integrable intensity ρ1=ρ\rho_1=\rho and pair correlation ρ2(x,y)=ρ2g(x−y)\rho_2(x,y)=\rho^2g(x-y), invariant under translations and one boost subgroup. If g−1∈L1(Rd)g-1\in L^1(\mathbb R^d) then g=1g=1 almost everywhere, and Var N(D)=ρV(D)\mathrm{Var}\,N(D)=\rho V(D) for every bounded DD. In particular such a process is not hyperuniform.

Proof. g−1g-1 is an integrable boost-invariant function, so Lemma 1 gives g−1=0g-1=0 almost everywhere. Then Var N(D)=∫Dρ1+∫D ⁣ ⁣∫D(ρ2−ρ1⊗ρ1)=ρV(D)\mathrm{Var}\,N(D)=\int_D\rho_1+\int_D\!\!\int_D(\rho_2-\rho_1\otimes\rho_1)=\rho V(D). □\square

The hypothesis g−1∈L1g-1\in L^1 is what hyperuniformity normally implies, through the sum rule 1+ρ∫(g−1)=01+\rho\int(g-1)=0. A hyperuniform process invariant under translations and one boost subgroup must therefore have pair correlations that are not integrable.

6. Classical lattices and complex weights

Sections 2 to 4 use the Fock structure of the counts and the positivity of the state. Without the Fock structure there are invariant classical processes that are not mixed Poisson (section 6.1). Without positivity, a process with complex weights of bounded total variation is a combination of classical ergodic processes (section 6.2).

6.1 Haar-random lattices are invariant and not Poisson

Let d≥2d\ge2, let X=ASL(d,R)/ASL(d,Z)X=ASL(d,\mathbb R)/ASL(d,\mathbb Z) with its invariant probability measure (finite by Siegel’s theorem [Siegel]), and let L=gZd+ξ\mathcal L=g\mathbb Z^d+\xi be the random affine lattice of density 1 with (g,ξ)(g,\xi) distributed by this measure. Concretely, gg is distributed by the invariant probability measure on SL(d,R)/SL(d,Z)SL(d,\mathbb R)/SL(d,\mathbb Z) and, given gg, ξ\xi is uniform in the torus Rd/gZd\mathbb R^d/g\mathbb Z^d. A typical lattice of this ensemble has no unusually short vector, but the measure gives positive weight to lattices with arbitrarily short ones, which consist of long, widely separated rows of closely spaced points. These produce the heavy tail of the counts seen below.

Proposition 9 (Haar-random lattices). The point process L\mathcal L on R1,d−1\mathbb R^{1,d-1} has the following properties.

(a) Its law is invariant under PP and under all of ASL(d,R)ASL(d,\mathbb R).

(b) For every bounded DD: EN(D)=V(D)\mathbb E N(D)=V(D) and E[N(D)(N(D)−1)]=V(D)2\mathbb E[N(D)(N(D)-1)]=V(D)^2, so Var N(D)=EN(D)\mathrm{Var}\,N(D)=\mathbb E N(D). The pair correlation is ρ2≡1\rho_2\equiv1.

(c) L\mathcal L is not a Cox process directed by any random measure without atoms; in particular it is neither Poisson nor mixed Poisson.

(d) L\mathcal L is ergodic under PP [Moore].

The proof is Theorem 1 of Lorentz-invariant point sets: what they can and cannot carry, which also describes the construction. It rests on three inputs. For (b), averaging over ξ\xi and applying Siegel’s mean value theorem [Siegel] gives E[N(N−1)]=V(D)2\mathbb E[N(N-1)]=V(D)^2. For (c), x≠y∈Lx\ne y\in\mathcal L forces 2y−x∈L∖{x,y}2y-x\in\mathcal L\setminus\{x,y\}, so the triples (x,y,2y−x)(x,y,2y-x) with x,yx,y in a bounded set AA number at least N(A)(N(A)−1)N(A)(N(A)-1), with expectation V(A)2>0V(A)^2>0, while a Cox process with non-atomic directing measure gives zero probability to a third point at 2y−x2y-x. For (d), Moore’s theorem [Moore] applies to SO+(1,d−1)SO^+(1,d-1) acting on SL(d,R)/SL(d,Z)SL(d,\mathbb R)/SL(d,\mathbb Z).

Proposition 9 shows that “Lorentz invariance and finite density force Poisson” is false for classical point processes. In a Fock state the counts in disjoint regions are independent in every coherent component (section 1.3), whereas in a lattice they are rigidly related: the pair x,yx,y forces the point 2y−x2y-x. The second-order statistics of the random lattice coincide with those of a Poisson process (ρ2≡1\rho_2\equiv1), so they cannot detect the difference, and neither can Corollary 8. The Lorentzian lattices of [SA] have smaller fluctuations of the number-volume correspondence in large regions, but are invariant only under a discrete group.

Void probabilities do separate the random lattice from a Poisson process. In two dimensions the law of N(D)N(D) is the same for every parallelogram DD of area VV, because the lattice law is invariant under SL(2,R)SL(2,\mathbb R) and translations, so it suffices to count the points in the square [0,a]2[0,a]^2 with V=a2V=a^2. We drew 1.6×1071.6\times10^{7} lattices (eight batches of 2×1062\times10^6), with the shape distributed by the hyperbolic measure on the modular fundamental domain, a uniform rotation, covolume 1 and ξ\xi uniform. The entries are means over the batches, with the batch-to-batch standard error in the last digits in parentheses.

VV Pr⁡(N=0)\Pr(N=0) Poisson e−Ve^{-V} V⋅Pr⁡(N=0)V\cdot\Pr(N=0)
1 0.30378 (9) 3.7×10−13.7\times10^{-1} 0.304
4 0.05432 (5) 1.8×10−21.8\times10^{-2} 0.217
16 0.01286 (3) 1.1×10−71.1\times10^{-7} 0.206

The void probability decays roughly like 1/V1/V and not exponentially: the product VPr⁡(N=0)V\Pr(N=0) stays near 0.2 for V≥4V\ge4, while e−Ve^{-V} falls by five orders of magnitude between V=4V=4 and V=16V=16. A heuristic accounts for the power law. A lattice avoids the square when its points lie on widely spaced rows, which requires a lattice vector of length at most of order 1/a1/a, and the invariant measure of such lattices is of order a−2=1/Va^{-2}=1/V. The same short vectors make the tail of the counts heavy. The invariant measure of lattices whose shortest vector is shorter than 1/n1/n is of order n−2n^{-2}. Such a lattice has rows spaced by more than nn, so a square of side of order one meets a row with probability of order 1/n1/n, and then contains of order nn points. This heuristic gives Pr⁡(N>n)\Pr(N>n) of order n−3n^{-3}: the variance is finite, as Proposition 9(b) proves, and the third moment would be infinite, which we do not prove. No formula for the constant 0.2 is available.

6.2 Complex weights of bounded total variation

A complex measure on configurations is a generalisation of a sprinkling in which the weights of configurations may be negative or complex; such weights occur in the study of the causal-set d’Alembertian on Mixed-Poisson and signed substrates and the four-dimensional d’Alembertian. A complex point process (CPP) is a family {PD}\{\mathsf P_D\}, indexed by bounded Borel sets, of complex measures on the space ΓD\Gamma_D of finite counting measures on DD, with PD(ΓD)=1\mathsf P_D(\Gamma_D)=1, finite total variation ∥PD∥<∞\lVert \mathsf P_D\rVert <\infty, and consistency: the image of PD′\mathsf P_{D'} under restriction to D⊂D′D\subset D' is PD\mathsf P_D. It is Poincaré invariant if PgD=g∗PD\mathsf P_{gD}=g_*\mathsf P_D. A positive process has ∥PD∥=1\lVert \mathsf P_D\rVert =1. Positivity entered the proof of Theorem 4 through the Husimi function, and a complex process does not have it. Theorem 10 describes the invariant complex processes whose total variation is bounded: they are combinations, with complex phases, of classical ergodic invariant processes, and they are classical when the bound is 1.

Theorem 10 (bounded total variation: classical ergodic components). Let P\mathsf P be a Poincaré-invariant CPP with S=sup⁡D∥PD∥<∞S=\sup_D\lVert \mathsf P_D\rVert <\infty. Then P\mathsf P extends to a complex measure on the space ΓM\Gamma_M of locally finite configurations, and

P=S∫eiφ(s) Rs κ(ds),(9)\mathsf P=S\int e^{i\varphi(s)}\,R_s\,\kappa(ds), \tag{9}

where κ\kappa is a probability measure, each RsR_s is an ergodic Poincaré-invariant classical (positive) point process, and φ\varphi is a real measurable function with S∫eiφdκ=1S\int e^{i\varphi}d\kappa=1. In particular, if S=1S=1 then φ=0\varphi=0 almost everywhere and P\mathsf P is a classical point process.

The total variations of the restrictions define a positive invariant measure ν\nu that dominates P\mathsf P. The density of P\mathsf P with respect to ν\nu is then invariant, hence constant on ergodic components, and its modulus is 1 because the total variation is exhausted by the restrictions.

Proof. Let DnD_n be the ball of radius nn and πD′→D\pi_{D'\to D} the restriction of configurations.

Step 1 (a dominating positive family). For a measurable map π\pi and a complex measure μ\mu, ∣π∗μ∣≤π∗∣μ∣\lvert \pi_*\mu\rvert \le\pi_*\lvert \mu\rvert. Hence ∣PD∣≤(πD′→D)∗∣PD′∣\lvert \mathsf P_D\rvert \le(\pi_{D'\to D})_*\lvert \mathsf P_{D'}\rvert for D⊂D′D\subset D', and the right side is nondecreasing in D′D' and has mass at most SS. The setwise limit νDm\nu_{D_m} of (πDn→Dm)∗∣PDn∣(\pi_{D_n\to D_m})_*\lvert \mathsf P_{D_n}\rvert as n→∞n\to\infty is then a measure of mass S=lim⁡∥PDn∥S=\lim\lVert \mathsf P_{D_n}\rVert, the family (νDm)(\nu_{D_m}) is consistent, and ∣PDm∣≤νDm\lvert \mathsf P_{D_m}\rvert \le\nu_{D_m}. The limit is the same along any increasing exhaustion of MM by bounded sets, so g∗νD=νgDg_*\nu_D=\nu_{gD}.

Step 2 (Kolmogorov). The spaces ΓDn\Gamma_{D_n} are standard Borel, and their projective limit is ΓM\Gamma_M. By Kolmogorov’s extension theorem for projective systems on standard Borel spaces [Parthasarathy] there is a unique finite measure ν\nu on ΓM\Gamma_M with marginals νDm\nu_{D_m}. It is invariant by Step 1, and has mass SS.

Step 3 (martingale). This step constructs the extension; Steps 1, 2 and 4 to 6 apply standard theorems. Since ∣PDn∣≤νDn\lvert \mathsf P_{D_n}\rvert \le\nu_{D_n} the density hn=dPDn/dνDnh_n=d\mathsf P_{D_n}/d\nu_{D_n} exists with ∣hn∣≤1\lvert h_n\rvert \le1. Consistency says that hnh_n is the conditional expectation of hn+1h_{n+1} on the configuration in DnD_n, so (hn)(h_n) is a bounded martingale for the filtration generated by the restrictions to DnD_n and converges ν\nu-almost everywhere and in L1L^1 to a function hh with ∣h∣≤1\lvert h\rvert \le1. The complex measure hνh\nu has the marginals PDn\mathsf P_{D_n}, and is the unique complex measure with these marginals, because cylinder sets form a π\pi-system that generates the σ\sigma-algebra. This is the extension of P\mathsf P. For projective families of complex Radon measures, uniform boundedness of the variation is necessary for such an extension to exist [Albeverio–Mazzucchi], and together with a tightness condition it is sufficient [Thomas]. The argument above does not need a separate tightness condition because the family is dominated by the positive projective family ν\nu, which extends by Kolmogorov’s theorem.

Step 4 (invariance). For g∈Pg\in P, g∗Pg_*\mathsf P and P\mathsf P have the same marginals, so g∗P=Pg_*\mathsf P=\mathsf P, and likewise g∗ν=νg_*\nu=\nu. Therefore h∘g=hh\circ g=h ν\nu-almost everywhere for each gg, and since PP is locally compact and second countable, hh agrees almost everywhere with a strictly invariant function [Zimmer].

Step 5 (ergodic decomposition). By Varadarajan’s theorem on ergodic decompositions [Varadarajan], ν/S=∫Rs κ(ds)\nu/S=\int R_s\,\kappa(ds) with ergodic invariant probability measures RsR_s. A strictly invariant function is constant on RsR_s-almost every point for κ\kappa-almost every ss; call the constant h(s)h(s). Then P=S∫h(s)Rsκ(ds)\mathsf P=S\int h(s)R_s\kappa(ds).

Step 6 (the modulus is 1). We have ∣P∣(ΓM)=∫∣h∣dν≤S\lvert \mathsf P\rvert (\Gamma_M)=\int\lvert h\rvert d\nu\le S. Conversely ∣P∣(ΓM)≥∥PDn∥→S\lvert \mathsf P\rvert (\Gamma_M)\ge\lVert \mathsf P_{D_n}\rVert \to S. Therefore ∫(1−∣h∣)dν=0\int(1-\lvert h\rvert )d\nu=0 with 1−∣h∣≥01-\lvert h\rvert \ge0, so ∣h∣=1\lvert h\rvert =1 ν\nu-almost everywhere, and h=eiφh=e^{i\varphi}. Normalisation P(ΓM)=1\mathsf P(\Gamma_M)=1 gives S∫eiφdκ=1S\int e^{i\varphi}d\kappa=1. If S=1S=1, a unimodular function with average 1 equals 1 almost everywhere. □\square

In (9) a process that is not classical has S>1S>1: its ergodic components are positive processes, and the phases eiφe^{i\varphi} produce cancellations between them. The simplest case is P=(1+ϵ)R1−ϵR2\mathsf P=(1+\epsilon)R_1-\epsilon R_2 for two mutually singular ergodic invariant processes R1,R2R_1,R_2, which has S=1+2ϵS=1+2\epsilon, κ\kappa the probability measure with weights (1+ϵ)/S(1+\epsilon)/S and ϵ/S\epsilon/S, and eiφ=±1e^{i\varphi}=\pm1. The components RsR_s are arbitrary ergodic invariant classical processes, not necessarily Poisson, and Proposition 9 gives non-Poisson examples. Bounded total variation therefore adds no ergodic components other than classical ones. The cancellations can still change the counting statistics. For a signed combination of Poisson processes with mean intensity ρ\rho, the second cumulant of N(D)N(D) computed from the correlation functions is ρV+V2μ2\rho V+V^2\mu_2, where μ2\mu_2 is the formal second central moment of the intensity under the signed mixing law. For the three-atom law of Example 1 on Mixed-Poisson and signed substrates and the four-dimensional d’Alembertian, whose total variation is at most 7/67/6, μ2=−0.47ρ2\mu_2=-0.47\rho^{2}, so the second cumulant is ρV−0.47ρ2V2<ρV\rho V-0.47\rho^2V^2<\rho V. The inequality Var N≥EN\mathrm{Var}\,N\ge\mathbb EN of Theorem 4 therefore fails for this cumulant, and since μ2<0\mu_2<0 some count probabilities are negative (Proposition 5 of the same page; for example Pr⁡(N=11)=−0.0087\Pr(N=11)=-0.0087 at ρV=2.6\rho V=2.6).

Complex processes with prescribed invariant correlation functions of moderate growth are constructed on the same page (Theorem 7), and for mixed Poisson processes with compactly supported complex mixing law it proves that requiring every count probability to lie in [0,1][0,1] makes the mixing law positive (Proposition 5).

7. Invariant kernels on de Sitter space

The proofs above use the non-compactness of the boosts, and a non-compact stabiliser plays the same role for other symmetry groups. On de Sitter space dS4=G/HdS_4=G/H with G=SO0(4,1)G=SO_0(4,1) and H=SO0(3,1)H=SO_0(3,1) the one-body density of an invariant state would be a positive-definite GG-invariant kernel. We do not set up states on de Sitter space and prove only the statement about kernels.

Proposition 11 (invariant positive-definite kernels on de Sitter space). Let k:dS4×dS4→Ck:dS_4\times dS_4\to\mathbb C be a continuous positive-definite kernel with k(gx,gy)=k(x,y)k(gx,gy)=k(x,y) for all g∈Gg\in G. Then kk is constant.

Proof. Let H\mathcal H be the Hilbert space obtained from finitely supported functions on dS4dS_4 with the inner product ⟨δx,δy⟩=k(x,y)\langle\delta_x,\delta_y\rangle=k(x,y). By invariance, π(g)δx=δgx\pi(g)\delta_x=\delta_{gx} defines a unitary representation of GG, and it is strongly continuous because ∥π(g)δx−δx∥2=2k(x,x)−2Re k(gx,x)→0\lVert\pi(g)\delta_x-\delta_x\rVert ^2=2k(x,x)-2\mathrm{Re}\,k(gx,x)\to0 as g→eg\to e. Fix x0x_0 with stabiliser HH and put v=δx0v=\delta_{x_0}, so π(h)v=v\pi(h)v=v for h∈Hh\in H. Write v=v0+v′v=v_0+v' with v0v_0 invariant under all of GG and v′v' in the orthogonal complement, which contains no non-zero GG-invariant vector. Then v′v' is HH-invariant. The group GG is connected, simple and has trivial centre, so by the Howe–Moore theorem [Howe–Moore] the matrix coefficient ⟨π(g)v′,v′⟩\langle\pi(g)v',v'\rangle tends to zero as gg leaves every compact set. For g=hng=h_n in the closed non-compact subgroup HH it equals ∥v′∥2\lVert v'\rVert ^2 for every nn, so v′=0v'=0. Hence vv is GG-invariant, δx=π(g)v=v\delta_x=\pi(g)v=v for all xx, and k(x,y)=∥v∥2k(x,y)=\lVert v\rVert ^2. □\square

8. Counting statistics without boost invariance or finite density

A fixed frame is a legitimate datum, for instance the rest frame of the cosmic microwave background. A quasi-free fermionic state whose TT is the multiplier by the indicator of a Euclidean ball in momentum space is translation invariant, has finite density and is not boost invariant. Steps 1 and 2 of Theorem 2 hold without boosts, with μ\mu any finite measure, so a Fermi-ball state is allowed. Which counting statistics such states give is open.

Without finite density, a bounded Lorentz-invariant multiplier, for example the indicator of {p0>0, 1≤(p0)2−∣p⃗∣2≤4}\{p^0>0,\ 1\le(p^0)^2-\lvert\vec p\rvert^2\le4\}, defines a Poincaré-invariant quasi-free CAR state with 0≤T≤10\le T\le1. For a unit vector h∈Lc2h\in L^2_c the number a∗(h)a(h)a^*(h)a(h) of the single mode hh has variance t(1−t)t(1-t) with mean t=⟨h,Th⟩t=\langle h,Th\rangle, so single-mode numbers are sub-Poisson. Here μ(R4)=∞\mu(\mathbb R^4)=\infty, so ω(N(D))=V(D)μ(R4)=∞\omega(N(D))=V(D)\mu(\mathbb R^4)=\infty for every DD by the computation in Step 2 of Theorem 2: sharp counts have no finite mean, and the conclusion of Theorem 6 fails without finite density. For bosons the mixtures (7) with ∫∣z∣2dm=∞\int\lvert z\rvert ^2d\mathsf m=\infty are invariant and locally normal, with infinite mean counts and Var N(D)=∞\mathrm{Var}\,N(D)=\infty.

Counts other than dΓ(1D)d\Gamma(1_D), for example a positive operator-valued measure on the one-particle space, are outside the setting; Integer counts and finite-dimensional algebras in a positive-energy vacuum treats integer additive counts of a positive-energy field. A state that is not locally normal has no single-mode reduction, and Proposition 3 does not apply to it. Complex processes of unbounded total variation lie outside Theorem 10.

  1. Corollary 8 excludes a hyperuniform invariant classical process when the pair correlation is ρ2g(x−y)\rho^2g(x-y) with g−1g-1 integrable. The orbit lattice of Kernel-scale kinematics of a single realisation is invariant in law and has a pair correlation carried by hyperboloids of fixed proper time, so Corollary 8 does not apply to it. The variance-to-mean ratio computed there for causal diamonds is about 0.48 to 0.49, 0.33 and 0.20 at volumes 4, 16 and 64 (averages over a ball of hyperbolic radius 5 to 6 approximating the invariant measure), and the proved bound Var N≤C2V3/2\mathrm{Var}\,N\le C^2V^{3/2} for a constant CC does not settle whether it tends to 0. Larger volumes, or a bound better than V3/2V^{3/2}, would settle whether a Poincaré-invariant classical point process can be hyperuniform.
  2. Fock states have independent counts in disjoint regions, given the intensity, and the random lattice lacks this property. Whether it characterises the mixed Poisson processes among classical invariant simple point processes with second moments is open.
  3. The mixtures (7) with infinite mean are the known invariant, locally normal states without (H3). Whether Theorem 4 extends to all such states, which would need a version of Theorem 2 with an infinite measure μ\mu, is open.

References

  1. Bombelli, Henson, Sorkin, “Discreteness without symmetry breaking: a theorem,” Mod. Phys. Lett. A 24 (2009) 2579, arXiv:gr-qc/0605006. Used for: no equivariant choice of a preferred direction exists on a Poisson sprinkling.
  2. Dowker, Henson, Sorkin, “Quantum gravity phenomenology, Lorentz invariance and discreteness,” Mod. Phys. Lett. A 19 (2004) 1829, arXiv:gr-qc/0311055. Used for: local Lorentz invariance of the discreteness of a sprinkled causal set, and the non-invariance of boosted regular lattices.
  3. Saravani, Aslanbeigi, “On the causal set-continuum correspondence,” Class. Quantum Grav. 31 (2014) 205013, arXiv:1403.6429. Used for: Poisson sprinkling is optimal for the number-volume correspondence in arbitrarily small regions; two-dimensional Lorentzian lattices do better in large regions.
  4. Hudson, Parthasarathy, “Quantum Itô’s formula and stochastic evolutions,” Commun. Math. Phys. 93 (1984) 301–323. Used for: the realisation of Poisson processes by operators on Fock space.
  5. Surya, “The causal set approach to quantum gravity,” Living Rev. Relativ. 22 (2019) 5, doi:10.1007/s41114-019-0023-1, arXiv:1903.11544. Used for: the causal-set setting and the definition of a sprinkling.
  6. Sorkin, “Causal sets: discrete gravity (notes for the Valdivia summer school),” arXiv:gr-qc/0309009. Used for: an introduction to causal sets and Poisson sprinkling.
  7. Reed, Simon, Methods of Modern Mathematical Physics II: Fourier Analysis, Self-Adjointness, Academic Press, 1975, section IX.2. Used for: the Bochner–Schwartz theorem (a positive-definite tempered distribution is the Fourier transform of a positive tempered measure).
  8. Soshnikov, “Determinantal random point fields,” Russian Math. Surveys 55 (2000) 923–975, doi:10.1070/RM2000v055n05ABEH000321, arXiv:math/0002099. Used for: determinantal processes with locally trace-class Hermitian kernels.
  9. Siegel, “A mean value theorem in geometry of numbers,” Ann. of Math. 46 (1945) 340, doi:10.2307/1969027. Used for: the mean value theorem for lattice sums over the space of unimodular lattices and the finiteness of its invariant measure.
  10. Moore, “Ergodicity of flows on homogeneous spaces,” Amer. J. Math. 88 (1966) 154, doi:10.2307/2373052. Used for: ergodicity of closed non-compact subgroups on the quotient of a simple group by a lattice.
  11. Thomas, “Projective limits of complex measures and martingale convergence,” Probab. Theory Related Fields 119 (2001) 579–588, doi:10.1007/pl00008772. Used for: the extension of a projective family of complex Radon measures of bounded variation under a tightness condition, in the form presented in [Albeverio–Mazzucchi].
  12. Albeverio, Mazzucchi, “A unified approach to infinite-dimensional integration,” Rev. Math. Phys. 28 (2016) 1650005, arXiv:1411.2853. Used for: uniform boundedness of the total variation as a necessary condition for a projective family of complex measures to arise from a measure on the limit.
  13. Parthasarathy, Probability Measures on Metric Spaces, Academic Press, 1967, Theorem V.3.2. Used for: Kolmogorov extension of projective systems of Borel measures on standard Borel spaces.
  14. Varadarajan, “Groups of automorphisms of Borel spaces,” Trans. Amer. Math. Soc. 109 (1963) 191–220, doi:10.1090/s0002-9947-1963-0159923-5. Used for: the ergodic decomposition of an invariant measure.
  15. Zimmer, Ergodic Theory and Semisimple Groups, Birkhäuser, 1984. Used for: an almost-everywhere invariant function for an action of a second-countable locally compact group has a strictly invariant version.
  16. Howe, Moore, “Asymptotic properties of unitary representations,” J. Funct. Anal. 32 (1979) 72–96, doi:10.1016/0022-1236(79)90078-8. Used for: matrix coefficients of a unitary representation of a connected simple Lie group with finite centre, without non-zero invariant vectors, vanish at infinity.