{"id":"1bd1fedc-e600-478d-a84e-0c959c74c32e","arxiv_id":"2412.02917","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Time-averaged vacuum energy flux in 2D CFT follows a symmetric variance-gamma distribution, and the joint flux-energy density distribution is explicitly constructed.","lead":"This paper derives exact probability formulas for vacuum energy flux fluctuations in two-dimensional quantum field theory, including the joint distribution of flux and energy density. The results may aid numerical simulations of quantum noise and analog experiments in one-dimensional condensed matter systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; Eq. (3.4) follows correctly from the stated shifted-Gamma premise, with only minor presentation issues in the illustrative example.","rationale":"The reader's verdict is CONDITIONAL, motivated mainly by the inconsistent use of a = 23/44 in Section V.B. I agree that this inconsistency should be fixed, and that the shifted-Gamma premise is the key input assumption. However, the premise is explicitly scoped by the paper, and the mathematical derivation of Eq. (3.4) from that premise is correct and transparent. I therefore find no load-bearing objection to the central claim; the issues are presentation-level and do not affect the flux distribution result. Keeping the verdict as CONDITIONAL is reasonable pending the trivial fix, but the correctness risk remains low.","tokens_in":13038,"tokens_out":21507,"duration_ms":183376,"concrete_test":"Check the original Antony-Fewster reference [3] for the exact domain of the parameter a in Eq. (5.3); if a = 23/44 is outside that domain, replace the illustrative example with a valid integer-parameter instance (for example, Gaussian sampling with central charge c = 24 and a = 0 gives alpha = 1) and verify that Eq. (3.4) remains unchanged.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation of Eq. (3.4) is mathematically sound: the convolution of two independent shifted Gamma variables is evaluated correctly, and the resulting Bessel-K distribution has the right normalization, variance, and asymptotics. The paper explicitly restricts attention to two classes of sampling functions for which the shifted Gamma marginals (2.3) are known, and it acknowledges in Section V.C that compactly supported sampling functions lead to different tail behavior. Thus the non-universality of Eq. (3.4) is disclosed rather than hidden, so I do not regard the premise about the input distributions as a load-bearing flaw. The one concrete internal inconsistency is in Section V.B, where the example a = 23/44 is used in Eqs. (5.1) and (5.3) even though Eq. (5.1) states that a is a nonnegative integer. This is a problem for that illustrative claim about reaching alpha = 1 with c = 1, but it does not affect the derivation or validity of Eq. (3.4). There is also a typo in Eq. (3.8), where the asymptotic for alpha = 1/2 should carry a minus sign so that the distribution is positive for small |F|. Neither issue undermines the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13294,"tokens_out":7463,"duration_ms":72400,"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.","major_comments":[],"minor_comments":[{"comment":"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":"Section V.A, after Eq. (5.13)"},{"comment":"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":"Section III, Eq. (3.8)"},{"comment":"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":"Section V.C and References"},{"comment":"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.","section":"Section V.B"}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nNew 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.\n\nStrengths: 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.\n\nSoft 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.\n\nBottom 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.","headline":"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.","tokens_in":13776,"tokens_out":2442,"would_cite":true,"duration_ms":61567,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T40","60E05"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["vacuum energy flux","probability distribution","two-dimensional conformal field theory","shifted Gamma distribution","variance-gamma distribution","modified Bessel function","stress tensor fluctuations","quantum inequalities"],"falsifier":"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.","tokens_in":12874,"feed_emoji":"⚛️","tokens_out":8316,"duration_ms":74369,"temperature":0.7,"pith_summary":"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.","feed_headline":"Vacuum flux odds in 2D follow a symmetric Bessel law","feed_subtitle":"In 2D conformal field theory, single flux measurements are centered at zero, exponentially tailed, and often sharply peaked at zero flux.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Establishes the shifted Gamma form for Gaussian averaged energy density, the base case used in Sec. V.","marker":"[1]"},{"why":"Derives shifted Gamma distributions for smeared stress tensor components in general 2D CFTs, the input assumption for the flux construction.","marker":"[2]"},{"why":"Supplies explicit values of alpha and beta for the Gaussian-like and Lorentzian-like sampling function classes used to plot the results.","marker":"[3]"},{"why":"Provides the integral identity (3.387.3) that turns the flux convolution integral into the modified Bessel function.","marker":"[13]"},{"why":"Supplies the DLMF standard formulae for Bessel integrals used in evaluating the flux density and its variance.","marker":"[14]"},{"why":"Gives the classical result that the difference of two identically distributed Gamma variables has the symmetric Bessel (variance-gamma) law.","marker":"[15]"},{"why":"Provides the moment problem theorem used to justify that factorized joint moments uniquely determine the joint distribution of flux and density.","marker":"[16]"}],"fun_headline_variants":["Vacuum flux in 2D: symmetric, Bessel-distributed","2D vacuum flux odds: even, Bessel-tailed, sharp at zero","Symmetric vacuum flux distribution emerges in 2D CFT","Flux fluctuations in 2D vacuum: a symmetric Bessel law","2D vacuum energy flux: symmetric distribution, Bessel core"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Vacuum flux in 2D: symmetric, Bessel-distributed","2D vacuum flux odds: even, Bessel-tailed, sharp at zero","Symmetric vacuum flux distribution emerges in 2D CFT","Flux fluctuations in 2D vacuum: a symmetric Bessel law","2D vacuum energy flux: symmetric distribution, Bessel core"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001318,"raw_usage":{"total_tokens":5395,"prompt_tokens":999,"completion_tokens":4396,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":4303}},"tokens_in":615,"tokens_out":4396,"duration_ms":30851,"temperature":1.0,"reasoning_tokens":4303,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:57:05.388423+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Probability distributions of smeared quantum stress tensors","cited_arxiv_id":"1004.0179","evidence_quote":"Establishes the shifted Gamma form for Gaussian averaged energy density, the base case used in Sec. V."},{"cited_title":"Probability distributions for the stress tensor in conformal field theories","cited_arxiv_id":"1805.04281","evidence_quote":"Derives shifted Gamma distributions for smeared stress tensor components in general 2D CFTs, the input assumption for the flux construction."},{"cited_title":"Explicit examples of probability distributions for the energy density in two-dimensional conformal field theory","cited_arxiv_id":"1908.00393","evidence_quote":"Supplies explicit values of alpha and beta for the Gaussian-like and Lorentzian-like sampling function classes used to plot the results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the DLMF standard formulae for Bessel integrals used in evaluating the flux density and its variance."},{"cited_title":"Kullback, The distribution laws of the difference and quotient of variables independently distributed in Pearson type III laws","cited_arxiv_id":null,"evidence_quote":"Gives the classical result that the difference of two identically distributed Gamma variables has the symmetric Bessel (variance-gamma) law."},{"cited_title":"Schm¨ udgen, The moment problem,Graduate Texts in Mathematics, vol","cited_arxiv_id":null,"evidence_quote":"Provides the moment problem theorem used to justify that factorized joint moments uniquely determine the joint distribution of flux and density."}],"review_version":1}