{"id":"429e997d-c916-4995-8c63-fc291e282c83","arxiv_id":"1908.03488","paper_version":3,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"The paper derives an averaged backreaction term with p = -ε from point-mass black hole fluctuations in an FLRW dust universe and claims their masses grow with cosmic time.","lead":"An astrophysicist argues that the combined gravitational ripples from many small black holes in an expanding universe act like a fluid with negative pressure, the same behavior as dark energy. If the mathematics were right, it would offer a classical explanation for cosmic acceleration and for how black holes grow.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (68), the paper's central p=-ε result, is asserted with no derivation and is preceded by fluctuation averages whose N-scaling and dimensions are inconsistent; a first-principles recomputation of δG^(2) is needed.","rationale":"I read the paper as an attempt to show that the quadratic backreaction of linearized spherically symmetric perturbations around N equal point masses in an FLRW dust background gives an effective energy-momentum tensor with p=-ε, thereby explaining cosmic acceleration, and that the same setup gives a mass-growth law m∼t. For the central claim, everything hinges on Eq. (68). The paper's own text supplies no derivation of that equation: it lists two averages, (64) and (67), asserts two more vanish in (66), and then states the result. This is not a matter of style; the second-order Einstein tensor depends on derivative correlations that are not given, and the N-scaling of (64) is opposite to what the linear superposition (61) and the averaging rule (62) imply. The reader's weakest_assumption identifies the same gap (the unexplained jump to δT and inconsistent N/dimensions in the finite averages), so I agree with the rejection. I would not adjust the reader's verdict. A secondary concern I noticed is the Conclusion's claim that the dust-stage growth (58) gives m∼t: with a(η)∝η² in the dust era, η²∝t^{2/3}, so the printed m∼t appears to be a separate error. I kept the primary concern on Eq. (68) because the p=-ε result is the paper's central claim.","tokens_in":9939,"tokens_out":14096,"duration_ms":151551,"concrete_test":"Recompute Eq. (68) from the stated ingredients: insert the metric (46) with h_ik from (61) and the profile ξ_a from (55) [or (59) on the dust branch] into the Einstein tensor, keep all terms quadratic in ξ and its derivatives, apply the averaging rule (62) for N identical sources uniformly distributed in V=4πr0³/3, and perform the renormalization ξ_a→ξ_a−3μ/5r0 described in §2.2. Evaluate symbolically (e.g., xTensor/Sage) and compare each component of ⟨δG^(2)i_k⟩ with -8π times the right-hand side of Eq. (68). As a cross-check, Monte Carlo the averages (64)-(67) for N=10,100,1000 to test whether ⟨ξ²⟩ scales as N or 1/N. If the tensor components do not match Eq. (68), the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Eq. (68). The entire dark-energy conclusion follows from the claim that the averaged second-order Einstein tensor equals a constant times δ^i_k. But the paper never displays δG^(2)i_k. In §2.3 it simply states \"Calculating now the averages ... let us find\" and jumps from the two averages (64) and (67) to (68). A second-order Einstein tensor for the metric (46)/(61) is a quadratic functional of ξ and its first derivatives; its average depends on ⟨ξ²⟩, ⟨(∂η ξ)²⟩, ⟨(∇ξ)²⟩ and cross terms, none of which are given except (64) and (67), and (66) is asserted without proof. Moreover, Eq. (64) has the wrong N dependence: with N uncorrelated identical sources h_ik=-a²δ_ikΣ_a ξ_a and the averaging rule (62), the mean square of the renormalized field should scale as N, whereas (64) has 1/N. Eq. (67) also has a dimension/typo problem in the printed factor. Since Eq. (68) scales as N, this is not harmless. Finally, the renormalization g^(0)→g (Eq. 20) plus the constant shift in (63) is assumed to leave physics unchanged; if the shift alters the macroscopic metric, the effective cosmological constant may be an artifact of that shift rather than of genuine backreaction. Therefore the p=-ε result is not supported by the displayed calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a self-consistent-field formalism for a macroscopic FLRW universe filled with dust and a population of point masses interpreted as black holes. Starting from linearized spherically symmetric perturbations, the author derives mass evolution laws for the ultrarelativistic and dust stages, then averages the perturbations using a prescribed coordinate-averaging rule. The central claim is that the averaged second-order corrections to the Einstein tensor produce an effective energy-momentum tensor δT^i_k = (9N/4r0^2)(2μ0/r0)^2 δ^i_k, which has the vacuum equation of state p = -ε and therefore acts like a cosmological constant. The paper also claims that the black-hole mass grows as m ∼ t during the matter-dominated stage.","tokens_in":10331,"tokens_out":4707,"duration_ms":54648,"significance":"If the central calculation were correct, the paper would offer a purely classical backreaction mechanism for cosmic acceleration and an attendant growth law for primordial black holes, with no new physics beyond general relativity. The formal setup, which adapts the self-consistent-field method to gravitational systems, is a reasonable starting point, and the explicit solutions for the mass evolution in Eqs. (54) and (58) are straightforward and clearly presented. However, the paper's main result, Eq. (68), is not derived: it is stated after two unevaluated averages and contains inconsistencies in N-scaling and dimensions in the preceding equations. The strength of the claimed conclusion is therefore not matched by the displayed computation.","major_comments":[{"comment":"The central result δT^i_k = (9N/4r0^2)(2μ0/r0)^2 δ^i_k is asserted without derivation. The text says 'Calculating now the averages from corrections to Einstein tensor' but never displays δG^(2)i_k. A second-order Einstein tensor built from the metric (46)/(61) is a quadratic functional of ξ and its first derivatives, so its average depends on ⟨ξ²⟩, ⟨(∂ηξ)²⟩, ⟨(∇ξ)²⟩ and cross terms; only (64) and (67) are given, and (66) is stated without proof. The jump from those averages to the isotropic form (68) is therefore unsupported.","section":"§2.3, Eq. (68)"},{"comment":"The N-scaling in Eq. (64) is inconsistent with the superposition (61). With hik = -a²δik Σ_a ξ_a and statistically independent sources, the renormalized average of (Σ_a ξ_a)² is N times the single-source variance, not 1/N times it. Eq. (64) gives ξ² ∝ 1/N, while Eq. (65) and the final result (68) carry factors of N. Since the claimed effect grows with the number of sources, this is not a harmless normalization convention; it indicates a missing or incorrect combinatorial factor in the averaging step.","section":"§2.2, Eqs. (61) and (64)"},{"comment":"Eq. (67) has a dimensional inconsistency: the left side ∂αξ∂βξ has dimension L^-2, whereas the right side (6πN/r0^2)(2μ0/r0^2) has dimension N L^-3 (taking ξ dimensionless and μ0/r0 dimensionless). Even allowing for a typo in the factor, the expression is not dimensionally homogeneous. Because this equation is one of the two inputs to the claimed δG^(2), the inconsistency propagates directly into the central result.","section":"§2.2, Eq. (67)"},{"comment":"The renormalization g(0) → gbar, and especially the constant shift in (63), is assumed to leave the physical macroscopic metric unchanged. If this shift is not merely a gauge or coordinate redefinition, the effective cosmological constant obtained in (68) could be an artifact of the renormalization rather than a genuine backreaction effect. The paper does not establish that the averaging procedure is gauge-invariant or that the shift in (63) is physically equivalent to a scale transformation of the background. This is a load-bearing assumption for the p = -ε conclusion.","section":"§1.2–§2.2, Eqs. (20) and (63)"}],"minor_comments":[{"comment":"The manuscript has numerous typographical and grammatical errors that impede reading, such as 'thepory', 'th e', 'av eraging', and the title's 'a black holes'. A careful English-language edit is needed.","section":"General"},{"comment":"The constants in Eq. (53) are not consistently named: the text introduces µ0 and µ1 but then writes m0 = µ1 + µ2, and the subsequent discussion assumes µ1 = 0 while referring to µ0. This should be clarified.","section":"Eq. (53)"},{"comment":"The Isaacson references are cited as Phys. Rev. 66, 1263 (1966) and 1272 (1966); the standard citation is Phys. Rev. 166, 1263 and 1272 (1968). Please verify and correct.","section":"References [20, 21]"}],"recommendation":"reject","confidential_remarks":"The paper's main claim, Eq. (68), cannot be reproduced from the displayed equations: the averaging calculation is not shown, the N-scaling is inconsistent between Eqs. (64) and (68), and Eq. (67) is dimensionally wrong. These are not local presentation issues but errors in the derivation of the central result. The manuscript also leans heavily on the author's prior work without enough detail for the reader to verify the crucial averaging steps. I therefore recommend rejection, although the underlying idea of averaging self-consistent gravitational perturbations is not without interest."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Ignat'ev is after something genuinely worth asking: can the averaged backreaction from many localized mass concentrations in a dust FLRW universe produce a term with p = -ε, and might that plausibly account for late-time acceleration? The setup — self-consistent field method, linearized spherically symmetric perturbations from earlier papers, averaging over uncorrelated source positions — is coherent, and the paper is honest about building on his own previous work. The mass-growth formula (58) tying black-hole mass to cosmological time is also an interesting consequence, though it does not follow from the paper's own equations in the way claimed.\n\nThe problem is that the main result is not actually derived. Equation (68), the effective stress tensor, appears after a one-line 'Calculating now...' with no expression for δG^(2) in terms of the averaged bilinears. The preceding averages do not provide enough information: (64) gives ⟨ξ²⟩ with a 1/N, while (68) is proportional to N. For N uncorrelated sources, the mean square of the summed field should scale as N, so (64) and (68) cannot both be right. Equation (67) also has a dimension problem: the right side as printed does not have units of a derivative-squared average. Then the jump to (68) is an unstated step, not a calculation.\n\nThere are softer worries too: the averaging rule and the renormalization shift (20)/(63) are not gauge-invariant, and the paper does not address that; Assumption 2 is plausible but strong; and the conclusion that m ~ t in the nonrelativistic stage is inconsistent with m ∝ η² in a dust universe where η ~ t^(1/3). These are not minor quibbles. The p = -ε result is the entire payload, and it is unsupported by anything shown.\n\nI don't see this as salvageable by light editing. It needs a genuine re-derivation of the averaged second-order Einstein tensor, with a clear calculation and consistent N-counting. For a reader interested in backreaction models of dark energy, the paper's premise is worth knowing about, but as a reference it is not usable in its current form.\n\nFor peer review: I would send it out once, not desk-reject, because the question is significant and the framework is not a fringe attempt. A competent referee will likely reject on the grounds above, but the exercise is worthwhile. I would not cite it myself until the calculation is fixed.","headline":"A physically motivated attempt to derive a cosmological constant from averaged black-hole perturbations, but the central calculation is missing and the few displayed averages contain inconsistent N-scaling and dimensional errors.","tokens_in":10757,"tokens_out":7345,"would_cite":false,"duration_ms":72787,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C05","83C55","83C57","83F05"],"pacs":["04.20.Cv","98.80.Cq","96.50.S","52.27.Ny"],"model":"deepseek-v4-flash","headline":"Averaged black-hole metric fluctuations generate a vacuum equation of state, so black-hole backreaction can act as a cosmological constant.","keywords":["macroscopic cosmology","black holes","equation of state","self-consistent field method","metric fluctuations","cosmological constant","Friedmann universe","backreaction"],"falsifier":"Compute the second-order backreaction of a population of point masses in a gauge-invariant cosmological perturbation scheme and compare the isotropic pressure correction to the energy-density correction: if the averaged correction is not exactly $\\delta p = -\\delta\\varepsilon$, the claim fails. Observationally, if the dark-energy equation of state measured from distance indicators deviates from $w = -1$ in a way that tracks the black-hole abundance, the predicted vacuum term would be excluded.","tokens_in":9745,"feed_emoji":"🕳️","tokens_out":7641,"duration_ms":69594,"temperature":0.7,"pith_summary":"This paper tries to show that a universe filled with a dust fluid and many identical black holes can be described macroscopically by the standard Einstein equations plus a correction term that behaves exactly like a cosmological constant. The argument treats the black holes as localized spherically symmetric metric perturbations, averages their combined effect over positions, and finds that the averaged quadratic fluctuations contribute an effective energy-momentum tensor with pressure equal to minus energy density, $p = -\\varepsilon$. If this is right, the backreaction of black holes is a classical mechanism for late-time accelerated expansion, requiring no new physics beyond general relativity. The paper also derives a growth law for the black-hole masses, $m \\sim t$, on the non-relativistic stage, meaning the holes can become very massive today.","feed_headline":"Black-hole backreaction mimics a cosmological constant","feed_subtitle":"Averaged metric fluctuations from many black holes add a vacuum term to the fluid's equation of state.","key_machinery":"The load-bearing mechanism is the self-consistent field method applied to gravity: one particle moves in the average field of all others, and the macroscopic metric is obtained by averaging over the particle-coordinate distribution. The specific technical object is the averaged bilinear combination of metric fluctuations, computed with the averaging rule $\\overline{\\phi} = \\prod_a \\frac{1}{V_a}\\int d^3 r_a \\, \\phi(r|x_1,\\dots,x_N)$, together with the assumptions that differentiation commutes with averaging and that the averaged Einstein tensor inherits the Friedmann symmetry. From the finite averages $\\overline{\\xi^2} = \\frac{108}{175N}\\left(\\frac{2\\mu_0}{r_0}\\right)^2$ and the renormalized derivative averages, the correction to the Einstein tensor is evaluated and converted into the effective stress tensor $\\delta T_i^k$. The mass-growth law follows from the separable linearized equations for the local perturbation and the point-mass mass function.","core_discovery":"The central claim is that the self-consistent-field averaging of linearized spherically symmetric perturbations produced by $N$ equal point masses in a Friedmann dust background yields an effective correction to the energy-momentum tensor $\\delta T_i^k = \\frac{9N}{4r_0^2}\\left(\\frac{2\\mu_0}{r_0}\\right)^2 \\delta_i^k$. Because this correction is proportional to $\\delta_i^k$, it has the vacuum equation of state $\\varepsilon_g + p_g = 0$; it adds a constant negative pressure to the fluid and acts as a cosmological constant. The paper further asserts that the growing mode of the mass function gives $m(\\eta) \\sim \\eta^2$ in the non-relativistic stage and hence $m \\sim t$, so the masses of the black holes grow to very large values at the modern stage. The macroscopic Einstein equations themselves are unchanged; what changes is only the effective equation of state of the fluid.","pith_inferences":["If this mechanism is real, the abundance and mass function of black holes could be constrained by the measured value of dark energy, since the correction scales as $N\\mu_0^2$.","A natural test is to repeat the averaging in a fully gauge-invariant way; the paper's result would be robust only if the vacuum equation of state survives that check.","The same averaging technique could apply to other localized sources, such as compact objects or topological defects, potentially generating similar effective negative-pressure fluids.","One could look for the predicted mass growth in the mass distribution of supermassive black holes across cosmic time; a growth law $m \\sim t$ would be a distinctive signature."],"forward_implications":["The effective stress-energy correction is isotropic and has $\\varepsilon_g + p_g = 0$, so it adds a constant negative pressure independent of the fluid density.","This term is mathematically equivalent to a cosmological constant, so it would drive accelerated expansion at late stages without introducing a scalar field.","The mass of a black hole grows as $m \\sim t$ in the non-relativistic stage, so black-hole masses can reach very large values today.","The correction does not alter the macroscopic Einstein equations themselves; it only changes the effective equation of state of the fluid."],"supporting_citations":[{"why":"Introduces the statistical theory and the averaging over particle coordinates that underlies the macroscopic Einstein equations.","marker":"[1, 2]"},{"why":"Presents the English version and the monograph containing the derivation of the averaged Einstein equations and the perturbation solutions used here.","marker":"[3, 4]"},{"why":"Derives the kinetic equation for massless particles in a Friedmann world with local spherically symmetric fluctuations generated by point masses.","marker":"[5, 6]"},{"why":"Provides the solution for the evolution of spherical local fluctuations in the ultrarelativistic case, giving the growing mode for the mass.","marker":"[7]"},{"why":"Extends the spherical-fluctuation evolution to arbitrary equations of state, including the non-relativistic case used for the $m \\sim t$ growth law.","marker":"[8, 9, 10]"},{"why":"Provides the WKB averaging of microscopic metric fluctuations that justifies replacing the microscopic Einstein equations by macroscopic ones with an effective stress.","marker":"[20, 21]"},{"why":"Gives the standard Lifshitz linear perturbation solutions for the Friedmann universe used to identify the growing and decreasing modes.","marker":"[26, 25]"}],"fun_headline_variants":["Black hole ensemble mimics vacuum energy","Averaged black holes yield a cosmic constant","Black hole backreaction adds a vacuum term","Many black holes impersonate dark energy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The main load-bearing premise is that the averaging rule over black-hole positions is a physically meaningful, gauge-invariant operation that commutes with derivatives and preserves the symmetry of the background; if that fails, the claimed vacuum equation of state does not follow.","fun_headline_variants_meta":{"raw":{"variants":["Black hole ensemble mimics vacuum energy","Averaged black holes yield a cosmic constant","Black hole backreaction adds a vacuum term","Many black holes impersonate dark energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000321,"raw_usage":{"total_tokens":1756,"prompt_tokens":840,"completion_tokens":916,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":456,"completion_tokens_details":{"reasoning_tokens":862}},"tokens_in":456,"tokens_out":916,"duration_ms":10058,"temperature":1.0,"reasoning_tokens":862,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:12:09.637357+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the second-order backreaction of a population of point masses in a gauge-invariant cosmological perturbation scheme and compare the isotropic pressure correction to the energy-density correction: if the averaged correction is not exactly $\\delta p = -\\delta\\varepsilon$, the claim fails. Observationally, if the dark-energy equation of state measured from distance indicators deviates from $w = -1$ in a way that tracks the black-hole abundance, the predicted vacuum term would be excluded.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the solution for the evolution of spherical local fluctuations in the ultrarelativistic case, giving the growing mode for the mass."}],"review_version":1}