Probability Distribution for Vacuum Energy Flux Fluctuations in Two Spacetime Dimensions

Christopher J. Fewster; L. H. Ford

Probability Distribution for Vacuum Energy Flux Fluctuations in Two Spacetime Dimensions

T0 review · 0 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The vacuum energy flux in two spacetime dimensions is governed by a symmetric probability distribution given by a modified Bessel function, with exponential tails and, for typical averaging functions, an integrable singularity at zero flux.

desk verdict A correct and clean derivation of the exact vacuum energy flux distribution in 2D CFT, with minor presentation issues that do not affect the main result. read the letter →

arxiv 2412.02917 v2 pith:FTDKILVW submitted 2024-12-04 hep-th gr-qcmath-phmath.MP

Christopher J. Fewster , L. H. Ford This is my paper · ORCID

classification hep-thgr-qcmath-phmath.MP MSC 81T4060E05

keywords vacuumenergyfluxprobabilitydistributiontwo-dimensionalconformalfieldtheoryshiftedGammavariance-gammamodifiedBesselfunctionstresstensorfluctuationsquantuminequalities

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper constructs the probability distribution for vacuum fluctuations of the energy flux in two-dimensional conformal field theory. Building on prior exact results that the time-averaged null stress tensor components follow shifted Gamma distributions, the authors show that their difference, the averaged flux, is distributed as a symmetric variance-gamma law: a modified Bessel function of the second kind. The distribution is even in the flux, decays exponentially for large flux, and for many sampling functions has an integrable singularity at zero flux. The authors also derive the joint distribution of flux and energy density, showing the flux is more centrally concentrated than the density, and obtain the flux distribution conditioned on negative energy density. These are among the few exact non-Gaussian probability laws known for quantum stress tensor fluctuations, and they provide concrete input for four-dimensional models, numerical simulations, and fluid analog experiments.

What carries the argument

The central object is the difference of two independent, identically distributed shifted Gamma random variables, $\omega_R - \omega_L$. The derivation is a convolution integral evaluated through the integral representation of the modified Bessel function $K_\nu$, which reduces the flux density to the symmetric variance-gamma distribution. Decoupling of the left- and right-moving components is proved in the appendix: the joint moment generating function factorizes by translation invariance and the cluster property. The joint density $$P(\rho,F) = \frac12 P_R\!\left(\tfrac12(\rho+F)\right) P_L\!\left(\tfrac12(\rho-F)\right)$$ for flux and density is the second load-bearing mechanism; it yields the marginals, the central-concentration probability, and the conditional flux distribution given negative energy density.

What would settle it

Simulate vacuum fluctuations of a massless scalar field in 1+1 dimensions on a lattice, time-average the flux operator with a Gaussian sampling function, and compare the histogram of outcomes to Eq. (3.4) with $\alpha = c/24$, $\beta = \pi$. A statistically significant deviation near $F=0$, where the integrable singularity exponent $2\alpha-1$ is predicted, or any asymmetry in the histogram would settle the claim.

Watch

Extended reading notes

Core claim

The central claim is that when the right- and left-moving null averaged stress tensor components $\omega_R$ and $\omega_L$ each have the shifted Gamma distribution (2.3), the time-averaged vacuum energy flux $F = \omega_R - \omega_L$ has the probability density $$P_F(F) = \frac{\$\beta$}{\sqrt{\pi}\,\Gamma(\$\alpha$)} \left(\frac{\$\beta$ |F|}{2}\right)^{\$\alpha$-1/2} K_{\$\alpha$-1/2}(\$\beta$ |F|),$$ a symmetric variance-gamma distribution independent of the shift parameter $\omega_0$. This distribution is even, has exponential tails with the same decay rate as the energy density, and exhibits an integrable singularity at $F=0$ when $\alpha \le 1/2$. The paper further claims that the joint distribution of flux and density is supported on $|F| \le \rho + 2\omega_0$, that $P(|F| < |\rho|) > 1/2$ (flux is more centrally concentrated than density), and that conditioning on negative energy density yields a flux distribution with compact support on $[-2\omega_0, 2\omega_0]$ and the same singularity structure.

Load-bearing premise

The entire construction assumes that each time-averaged null stress tensor component $\omega_{L,R}$ has the exact shifted Gamma distribution (2.3); that premise has been established only for particular classes of sampling functions, so if a physically relevant smearing function (for example, a compactly supported one) produces different marginals, the Bessel flux formula (3.4) will not hold for it.

Editorial extensions

If this is right

Individual flux measurements in a 2D CFT vacuum are symmetrically distributed about zero, so positive and negative outcomes are equally likely and large fluctuations are exponentially suppressed with the same decay constant as the energy density.
For sampling functions with $\alpha \le 1/2$ (including Gaussian and Lorentzian averaging of a massless scalar field), the flux distribution has an integrable singularity at $F=0$, making tiny outcomes overwhelmingly likely and the cumulative distribution nearly a step function.
The flux is typically more centrally concentrated than the energy density: $P(|F| < |\rho|) = 1 - 2p + 2p^2 > 1/2$, and for Gaussian averaging of a scalar field this probability is about 0.81.
Conditioning on a negative energy density forces the flux distribution onto the compact interval $[-2\omega_0, 2\omega_0]$, with the same near-zero singularity behavior as the unconditional distribution.
The explicit cumulative distribution in terms of Bessel and Struve functions gives a ready algorithm for Monte Carlo simulation of vacuum flux fluctuations, including correlations between different times.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The variance-gamma form means the flux can be represented as a Gaussian with a Gamma-distributed variance, so these vacuum fluctuations connect to standard normal-variance mixture models; one could use that representation to generate correlated multi-time flux samples.
The paper's own asymptotic result for compactly supported smearing functions is a stretched exponential $e^{-\beta |F|^{\alpha_p}}$ rather than the Bessel form, suggesting the exact distribution for compact support would interpolate between the two and giving a concrete test of how strongly the measurement process shapes the tail.
The sharp signature of an integrable singularity with exponent $2\alpha - 1$ at zero flux, together with the compact support of the conditioned distribution, could be looked for in one-dimensional phonon analog systems.
Factorized left/right moments imply vanishing flux-density correlations for null-separated regions; a lattice simulation could check this directly, since it is the cluster-property step that makes the whole joint distribution product form.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: the paper does exactly what it says. It takes the known shifted Gamma marginals for smeared null energy components in 2D CFT, convolves them to get the flux distribution, and lands on a variance-gamma distribution with a Bessel K form. The core derivation is correct, the presentation is transparent, and the joint distribution analysis in Sec IV adds something the earlier energy-density papers did not have.

New here is not the classical difference-of-Gamma identity — that goes back to Kullback — but the physical package: explicit alpha and beta for Gaussian and Lorentzian averaging, the symmetric flux distribution with exponential tails, the integrable singularity at zero for alpha <= 1/2, and the conditioning on negative energy density. For people who simulate stress tensor fluctuations or work on 1D analog models, this is directly usable.

Strengths: the convolution in Eq (3.3) is clean; the normalization and variance check out; the joint distribution via the product of marginals is justified by the factorization argument in Appendix A, which is a nice addition. The paper also says plainly that compactly supported sampling functions lead to different tail behavior, so it is not overclaiming universality.

Soft spots, in order: (1) Section V.A uses a=23/44 as an example even though Eq (5.1) says a is a nonnegative integer. That is an internal contradiction; either a should be restricted or the analytic continuation should be justified. It is illustrative only, so it does not damage Eq (3.4). (2) Eq (3.8) is missing a minus sign in front of the logarithm; as printed it gives a negative density near zero for alpha=1/2. Clearly a typo. (3) The whole construction inherits the shifted Gamma premise from Refs [1-3]; if a sampling function outside those two classes gives different marginals, the Bessel form won't hold. The paper says as much in Sec V.C, so this is a boundary condition, not a hidden flaw.

Bottom line: this is a solid, honest paper with a correct central result and a couple of small presentation errors. The right venue will want the typos fixed, but the referee time is warranted. I'd send it to review.

Referee Report

0 major / 4 minor

Summary. The paper derives the probability distribution for time-averaged vacuum energy flux fluctuations in (1+1)-dimensional conformal field theory. It assumes the right- and left-moving averaged null stress-tensor components ω_R and ω_L are independent and each have the shifted Gamma distribution (2.3), as established in earlier work for Gaussian-like and Lorentzian-like sampling functions. Under this premise, the flux F = ω_R − ω_L is shown in Sec. III to have the symmetric Bessel-K distribution (3.4); the paper derives its large- and small-argument asymptotics, variance (3.10), and cumulative distribution (3.12). Sec. IV constructs the joint distribution of flux and energy density, exhibits support restrictions, computes Prob(|F| < |ρ|), and analyzes the flux conditioned on negative energy density, including the singular behavior at the support endpoints. Sec. V gives explicit parameters for Gaussian-like and Lorentzian-like sampling functions and discusses compactly supported (non-shifted-Gamma) sampling. The paper concludes with applications to four-dimensional models, condensed matter analogs, and numerical simulations.

Significance. If correct, Eq. (3.4) is a clean, exact result for flux fluctuations in 2D CFT, complementary to the shifted Gamma result for the energy density. The derivation is transparent and checks out: the convolution is evaluated correctly, and the normalization, variance, and symmetry of the resulting distribution are mutually consistent. The paper explicitly and appropriately limits the claims to sampling functions whose marginals are shifted Gamma (Sec. V.C), so the non-universality of Eq. (3.4) is disclosed rather than hidden. The joint-distribution construction and the conditional flux distribution are useful extensions. No free parameter is fitted to reach the main result; the input distributions come from previously published independent results. The main caveat is that the premise (2.3) is not universal, so Eq. (3.4) should not be read as a universal law for all sampling functions.

minor comments (4)

[Section V.A, after Eq. (5.13)] The illustrative example 'a = 23/44' contradicts the definition of the Gaussian-like class in Eq. (5.1), which requires a to be a nonnegative integer. Since 23/44 is not an integer, the claimed case α = 1 with c = 1 is not a member of that class. Please replace it with a valid example (for instance c = 8, a = 1, which gives α = 1 in Eq. (5.3)) or explicitly extend the domain of a in Eq. (5.1) with a supporting derivation.
[Section III, Eq. (3.8)] The asymptotic formula for α = 1/2 has a sign error: since K_0(z) ~ -log(z/2) for small z, the right-hand side should be -(β/π) log(β|F|/2), not +(β/π) log(β|F|/2). As written, the displayed distribution would be negative for small |F|.
[Section V.C and References] There are several typographical errors: 'foe' should be 'for' in 'appropriate descriptions foe physical measures', 'was assume to model' should be 'was assumed to model', and Ref. [19] gives the arXiv identifier as '2409:02855', which should be '2409.02855'.
[Section V.B] In the Lorentzian case, the sentence 'Here PF ∝ |ω|^{-17/18}' uses ω where the independent variable has been denoted F throughout; please use F for consistency.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: Eq. (3.4) is an honest convolution of independently established shifted-Gamma marginals, and the self-citations that occur are not load-bearing.

full rationale

Section III takes as its premise the shifted Gamma marginal distributions of Eq. (2.3), which the paper attributes to Refs. [1–3]. The flux F is defined in Eq. (3.1) as the difference of the two null stress components, and Eq. (3.2) is a genuine convolution over the independent marginals; the evaluation to Bessel-K form in Eqs. (3.3)–(3.4) uses standard integral identities (Gradshteyn–Ryzhik 3.387.3, DLMF 10.32.8), with no fitted parameter and no use of the flux distribution to define its inputs. Because the shifted-Gamma inputs are prior analytical results for explicitly listed sampling-function classes, and because Section V.C explicitly states that compactly supported functions lead to different tails, the conditional character of Eq. (3.4) is transparent rather than circular. The self-citations in Section IV (e.g., Ref. [1] for the Hamburger moment condition) are not load-bearing: the product form Eq. (4.3) follows from the explicit independent densities and the factorization argument of Appendix A, and the moment-determinacy remark is not needed to obtain the later joint-distribution formulas. I therefore find no circular reduction. Separately, and outside circularity, the illustrative example in Section V.B sets a = 23/44 in Eqs. (5.1) and (5.3) although Eq. (5.1) states that a is a nonnegative integer, and Eq. (3.8) appears to have a sign error for α = 1/2; these are correctness or typographical issues that do not affect the derivation of Eq. (3.4).

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper adds no free parameters of its own; α and β are analytical inputs from previous work. The main axioms are the shifted Gamma marginal distributions, statistical independence of left/right movers, and standard QFT cluster properties.

assumptions (4)

domain assumption The averaged null stress tensor components PL,R have shifted Gamma distributions (Eq. 2.3) for the sampling functions considered.
Taken as established from Refs. [1-3]; the flux derivation depends entirely on this form.
domain assumption The vacuum is translation invariant and satisfies the cluster property, ensuring factorization of joint moments (Appendix A).
Standard QFT properties; argued in Appendix A of the paper.
domain assumption The right- and left-moving components have the same central charge and are statistically independent.
Standard parity-invariant 2D CFT assumption; the convolution in (2.5) relies on it.
standard math Moment determinacy via Petersen's theorem (Ref. [16]).
Used to justify that the joint distribution is determined by moments; standard mathematical result.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Probability Distribution for Vacuum Energy Flux Fluctuations in Two Spacetime Dimensions." pith.science (2026). https://pith.science/paper/FTDKILVW

@misc{pith2026241202917,
  author       = {Pith},
  title        = {Pith review of: Probability Distribution for Vacuum Energy Flux Fluctuations in Two Spacetime Dimensions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FTDKILVW}},
  note         = {Machine review of arXiv:2412.02917}
}

read the original abstract

The probability distribution for vacuum fluctuations of the energy flux in two dimensions will be constructed, along with the joint distribution of energy flux and energy density. Our approach will be based on previous work on probability distributions for the energy density in two dimensional conformal field theory. In both cases, the relevant stress tensor component must be averaged in time, and the results are sensitive to the form of the averaging function. Here we present results for two classes of such functions, which include the Gaussian and Lorentzian functions. The distribution for the energy flux is symmetric, unlike that for the energy density. In both cases, the distribution may possess an integrable singularity. The functional form of the flux distribution function involves a modified Bessel function, and is distinct from the shifted Gamma form for the energy density. By considering the joint distribution of energy flux and energy density, we show that the distribution of energy flux tends to be more centrally concentrated than that of the energy density. We also determine the distribution of energy fluxes, conditioned on the energy density being negative. Some applications of the results will be discussed.

Figures

Figures reproduced from arXiv: 2412.02917 by the authors.

**Figure 1.** Here PF ∝ |F| −11/12 near the origin. The cumulative distribution function P<(F) is plotted for the Gaussian averaged scalar field case with c = 1 in [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗

**Figure 2.** Note that the integrable singularity in PF (F) leads to P<(F) approximating a step function. Numerical evaluation gives P<(10−10) = 0.585 and P<(−10−10) = 0.415, so the average slope in the interval |F| < 10−10 is ∆P< ∆F ≈ 8.5 × 108 . (5.6) which is consistent with the lower panel in [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

**Figure 2.** FIG. 2. The cumulative probability distribution, [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

Figures from the paper (2 more)

**Figure 3.** Figure 3: FIG. 3. The probability distribution of flux, conditioned on the energy density being negative. [PITH_FULL_IMAGE:figures/full_fig_p009_3.png]

**Figure 4.** Figure 4: FIG. 4. The cumulative probability distribution, [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 21 canonical work pages

[1]

C. J. Fewster, L. H. Ford and T. A. Roman, Probability distributions of smeared quantum stress tensors, Phys. Rev. D 81, 121901(R) (2010), arXiv:1004.0179 [quant-ph]. 12

work page Pith review arXiv 2010
[2]

C. J. Fewster and S. Hollands, Probability distributions for the stress tensor in conformal field theories, Lett. Math. Phys. 109, 747 (2019), arXiv:1805.04281

work page Pith review arXiv 2019
[3]

M. C. Antony and C. J. Fewster, Explicit examples of probability distributions for the energy density in two-dimensional conformal field theory, Phys. Rev. D 101, 025010 (2020), arXiv:1908.00393

work page Pith review arXiv 2020
[4]

C. J. Fewster and S. Hollands, Quantum energy inequalities in two-dimensional conformal field theory, Rev. Math. Phys. 17, 577 (2005), arXiv:math-ph/0412028

arXiv 2005
[5]

C. J. Fewster and L. H. Ford, Probability distributions for quantum stress tensors measured in a finite time interval, Phys. Rev. D 92, 105008 (2015), arXiv:1508.02359

work page Pith review arXiv 2015
[6]

C. J. Fewster, L. H. Ford and T. A. Roman, Probability distributions for quantum stress tensors in four dimensions, Phys. Rev. D 85, 125038 (2012), arXiv:1204.3570 [quant-ph]

work page Pith review arXiv 2012
[7]

C. J. Fewster and L. H. Ford, Probability distributions for space and time averaged quantum stress tensors, Phys. Rev. D 101, 025006 (2020), arXiv:1909.07295

work page Pith review arXiv 2020
[8]

E. D. Schiappacasse, C. J. Fewster and L. H. Ford, Vacuum quantum stress tensor fluctuations: A diagonalization approach, Phys. Rev. D 97, 025013 (2018), arXiv:1711.09477 [hep-th]

work page Pith review arXiv 2018

Show all 23 references

[9]

P. Wu, E. D. Schiappacasse and L. H. Ford, Space and time averaged quantum stress tensor fluctuations, Phys. Rev. D 103, 125014 (2021), arXiv:2104.04446 [hep-th]

2021 arXiv
[10]

P. Wu, E. D. Schiappacasse and L. H. Ford, Frequency spectra analysis of space- and time-averaged quantum stress tensor fluctuations, Phys. Rev. D 107, 036013 (2023), arXiv:2211.12001

2023 arXiv
[11]

Carlip, R

S. Carlip, R. A. Mosna and J. P. M. Pitelli, Vacuum fluctuations and the small scale structure of spacetime, Phys. Rev. Lett. 107, 021303 (2011), arXiv:1103.5993

2011 arXiv
[12]

Carlip, Ricardo A

S. Carlip, Ricardo A. Mosna, J. P. M. Pitelli, Quantum fields, geometric fluctuations, and the structure of spacetime, Phys. Rev. D 102, 126018 (2020), arXiv:1809.08265

2020 arXiv
[13]

I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series, and products, 7th ed. Academic Press (2007)

2007
[14]

NIST Digital library of mathematical functions, https://dlmf.nist.gov/, Release 1.2.2 of 2024-09-15. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, eds

2024
[15]

Kullback, The distribution laws of the difference and quotient of variables independently distributed in Pearson type III laws

S. Kullback, The distribution laws of the difference and quotient of variables independently distributed in Pearson type III laws. The Annals of Mathematical Statistics 7, 51–53 (1936)

1936
[16]

Schm¨ udgen, The moment problem,Graduate Texts in Mathematics, vol

K. Schm¨ udgen, The moment problem,Graduate Texts in Mathematics, vol. 277. Springer, Cham (2017)

2017
[17]

Wu and L

P. Wu and L. H. Ford, Large zero point density fluctuations in fluids, Phys. Rev. Research 2, 032028 (2020), arXiv:2005.04266

2020 arXiv
[18]

E. R. Taylor, S. Yencho and L. H. Ford, Numerical simulation of quantum field fluctuations, Phys. Rev. D 109 116010 (2024), arXiv:2312.17155

2024 arXiv
[19]

L. H. Ford, Vacuum radiation pressure fluctuations on electrons, Phys. Rev. D 110, 096002 (2024), arXiv:2409:02855

2024
[20]

L. H. Ford and and T. A. Roman, Energy Density-Flux Correlations in an Unusual Quantum State and in the Vacuum, Phys. Rev. D 76, 064012 (2007), arXiv:0706.1970v2

2007 arXiv
[21]

R. F. Streater, A. S. Wightman, PCT, spin and statistics, and all that. Princeton Landmarks in Physics. Princeton University Press, Princeton, NJ (2000). Corrected third printing of the 1978 edition

2000
[22]

Maison, Eine Bemerkung zu Clustereigenschaften

D. Maison, Eine Bemerkung zu Clustereigenschaften. Comm. Math. Phys. 10, 48–51 (1968)

1968
[23]

Petrov, Asymptotic behaviour of the unitary representations of the Poincar´ e group

A. Petrov, Asymptotic behaviour of the unitary representations of the Poincar´ e group. Rep. Mathematical Phys.4, 183–188 (1973)

1973

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.

[1] [1]

C. J. Fewster, L. H. Ford and T. A. Roman, Probability distributions of smeared quantum stress tensors, Phys. Rev. D 81, 121901(R) (2010), arXiv:1004.0179 [quant-ph]. 12

work page Pith review arXiv 2010

[2] [2]

C. J. Fewster and S. Hollands, Probability distributions for the stress tensor in conformal field theories, Lett. Math. Phys. 109, 747 (2019), arXiv:1805.04281

work page Pith review arXiv 2019

[3] [3]

M. C. Antony and C. J. Fewster, Explicit examples of probability distributions for the energy density in two-dimensional conformal field theory, Phys. Rev. D 101, 025010 (2020), arXiv:1908.00393

work page Pith review arXiv 2020

[4] [4]

C. J. Fewster and S. Hollands, Quantum energy inequalities in two-dimensional conformal field theory, Rev. Math. Phys. 17, 577 (2005), arXiv:math-ph/0412028

arXiv 2005

[5] [5]

C. J. Fewster and L. H. Ford, Probability distributions for quantum stress tensors measured in a finite time interval, Phys. Rev. D 92, 105008 (2015), arXiv:1508.02359

work page Pith review arXiv 2015

[6] [6]

C. J. Fewster, L. H. Ford and T. A. Roman, Probability distributions for quantum stress tensors in four dimensions, Phys. Rev. D 85, 125038 (2012), arXiv:1204.3570 [quant-ph]

work page Pith review arXiv 2012

[7] [7]

C. J. Fewster and L. H. Ford, Probability distributions for space and time averaged quantum stress tensors, Phys. Rev. D 101, 025006 (2020), arXiv:1909.07295

work page Pith review arXiv 2020

[8] [8]

E. D. Schiappacasse, C. J. Fewster and L. H. Ford, Vacuum quantum stress tensor fluctuations: A diagonalization approach, Phys. Rev. D 97, 025013 (2018), arXiv:1711.09477 [hep-th]

work page Pith review arXiv 2018