{"id":"f03830a0-2683-4aac-a9e0-668750cf38b0","arxiv_id":"2509.02701","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A single gravitational wave template combining inflaton and evaporating primordial black hole sources, with and without memory burden, yields new frequency-ratio signatures for testing black hole evaporation.","lead":"This paper maps the full gravitational wave spectrum from an early universe in which tiny primordial black holes dominate and then evaporate to reheat the cosmos. It shows how the memory burden effect, a proposed slowing of black hole evaporation, would move and reshape the spectrum and add a telltale double peak.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (92), advertised as giving 'a direct access to n' from a double peak, fails the n=0, q=1 semiclassical limit and is not obtained by dividing Eq. (90) by Eq. (85), so the central burden diagnostic is not established.","rationale":"I read the paper as a phenomenological synthesis: the unburdened spectral template obtained by combining the inflationary plateau, PBH-domination tilt, density-fluctuation peak, and PBH-evaporation peak is coherent, and the text is transparent about approximations such as gamma and the cutoff. The reader's conditional verdict is therefore appropriate. However, the reader located the load-bearing weakness in the physical uncertainty of the memory-burden ansatz (q,n free). I find a sharper, model-independent problem: the diagnostic ratio Eq. (92) that is supposed to make n observable is not consistent with the two frequency formulas from which it is derived. The n=0 limit alone is decisive, and the exponent comparison shows the displayed formula cannot be a typo in a single coefficient. Because this concern invalidates the strongest advertised observable consequence of the memory-burden extension, the conditional verdict is kept, but for this additional reason. The concrete check is analytic and takes minutes; until it is resolved, the double-peak signature should be presented as an illustration of the ansatz rather than as a means to measure n.","tokens_in":35878,"tokens_out":11597,"duration_ms":98264,"concrete_test":"Divide Eq. (90) by Eq. (85) by hand and compare the resulting q- and (M_in/M_P)-exponents with Eq. (92); then evaluate the ratio at n=0, q=1, where it must be 1 because both peaks reduce to Eq. (82). If the exponents do not match or the n=0 limit is not 1, Eq. (92) must be corrected and the claimed direct extraction of n from an observed double peak should be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Even granting the memory-burden ansatz Eq. (12) and the acknowledged freedom in q and n, the paper's headline double-peak diagnostic is internally inconsistent. In the semiclassical limit n=0, q=1, Eqs. (85) and (90) both reduce to Eq. (82), so fev2_0/fev1_0 must equal 1. Substituting n=0, q=1 into Eq. (92) gives a factor (2n)^{2/3}=0, hence a vanishing ratio; Eq. (92) cannot be correct in the limit it is supposed to reproduce. Directly dividing Eq. (90) by Eq. (85) gives fev2_0/fev1_0 = 3^{-1/6} 2^{n/3}(3+2n)^{1/3} q^{2n/3} (M_in/M_P)^{4n/3}, whereas Eq. (92) states 1.2 q^{1+4n/3} [2n(3+2n)]^{2/3} (M_in/M_P)^{4n/3-1}: the q and mass exponents and the n=0 prefactor all disagree. Since Sec. VI and the abstract use this ratio as the observational signature of memory burden, the central claim rests on a formula that does not follow from the displayed equations. This is a correctness issue independent of the physical uncertainty in q and n.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the total present-day gravitational-wave background from a post-inflationary phase in which primordial black holes reheat the Universe while the inflaton oscillates in a V(phi) proportional to phi^k potential. It combines four contributions: inflationary tensor modes re-entering during radiation, PBH domination, and inflaton domination (regions A1-A3); gravitons from inflaton scattering (region C); scalar-induced gravitational waves from PBH number-density fluctuations (regions B1/B2); and direct Hawking graviton emission (region D1). A memory-burden extension is added in which evaporation is suppressed by a factor 1/S^n after the PBH mass drops below q M_in, producing the additional features D'_1 and D'_2. The advertised observational signatures are the low-frequency enlargement of the spectrum under memory burden, the frequency ratio f_UV0/f_IR0, and especially the double-peak ratio f_ev2,0/f_ev1,0, which is claimed to give 'a direct access to n'.","tokens_in":36231,"tokens_out":17231,"duration_ms":133480,"significance":"If the central formulas were correct, this would be a useful unified template for gravitational-wave searches in PBH-reheating scenarios. The paper has genuine strengths: the PGW matching between eras is internally consistent (for example, Eq. (41) and Eq. (45) meet at f_BH0); the authors explicitly flag the non-linear cutoff ambiguity in Appendix D.4 and display its impact in Fig. 9; and they are transparent about the absence of a microscopic derivation for gamma and for the burden-onset parameter q. The burdened UV/IR relation Eq. (79) follows from Eqs. (76) and (46). However, the diagnostic that the paper emphasizes most, the two-peak ratio of Eq. (92), is not supported by the displayed algebra, so the central memory-burden conclusion is not yet established.","major_comments":[{"comment":"The ratio f_ev2,0/f_ev1,0 is announced as the direct probe of n, but Eq. (92) does not follow from Eqs. (85) and (90). Dividing Eq. (90) by Eq. (85) as printed gives f_ev2,0/f_ev1,0 = 3^(-5/6) 2^(-7n/6) (3+2n)^(-2/3) q^(2n/3) (M_in/M_P)^(4n/3), whereas Eq. (92) states 1.2 q^(1+4n/3) [2n(3+2n)]^(2/3) (M_in/M_P)^(4n/3-1). The two expressions disagree in the q exponent, the mass exponent, and the n=0 prefactor. In the semiclassical limit n=0, q=1, both Eqs. (85) and (90) are claimed to reduce to Eq. (82), so the ratio must tend to 1; Eq. (92) instead gives 0 through the factor (2n)^(2/3). Since Sec. VI and the abstract use this ratio as the observational signature of memory burden, this is a load-bearing error that must be repaired.","section":"Sec. VI.B, Eq. (92)"},{"comment":"The two defining equations are inconsistent in the limit they are both stated to reproduce. Substituting n=0, q=1 into Eq. (85) gives a prefactor proportional to 3^(7/12), while Eq. (90) gives a prefactor proportional to 3^(-1/4) for the same physical frequency; the two expressions therefore differ by a factor 3^(5/6) and cannot both equal Eq. (82). The origin appears to lie in the orientation of the scale-factor ratio in Eq. (86) and in the q- and 3-dependences of Eq. (85). The double-peak calculation should be re-derived from the definitions of a_ev1, a_RH, and T_RH before the benchmark frequencies in Eqs. (88) and (91) are quoted.","section":"Sec. VI.B, Eqs. (85) and (90)"},{"comment":"The claim that observing the double peak gives 'a direct access to n' should be qualified by the explicit model dependence of Eq. (12). Since q and n are free parameters and the paper itself states that there are no microscopic motivated values for q, the peak positions and their ratio are predictions of the power-law entropy-suppression ansatz, not of memory burden in general. The paper should present at least a concrete cross-check, for example the behavior of f_ev2,0/f_ev1,0 as n tends to 0 at fixed q, and should state what would change for a non-power-law suppression of the evaporation rate.","section":"Sec. VI.B and Sec. II.C"}],"minor_comments":[{"comment":"The sentence 'on needs to solve the Boltzmann equation' should read 'one needs to solve the Boltzmann equation'.","section":"Sec. IV.B"},{"comment":"The word 'wich' should be 'which' in the sentence describing the scalar fluctuations.","section":"Sec. V.A"},{"comment":"The word 'sumarized' should be 'summarized', and the sentence 'The usual workaround for PBH evaporation is to get the generic instantaneous emission before integrating it with respect to the energy in order to get the proper mass evolution' is grammatically incomplete.","section":"Appendix B.2"},{"comment":"The phrase 'For the first time' is a strong claim; given the overlap with Refs. [30, 80, 81], the phrasing should be softened or the precise new element should be stated explicitly.","section":"Abstract"},{"comment":"The parenthetical 'with g_RH = 106.75' is not closed; please close the parenthesis.","section":"Eq. (84)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the paper is a competent assembly of known ingredients, and the main concern in my report is purely algebraic: Eq. (92) is the paper's most prominent new diagnostic and it does not follow from the equations it claims to combine. This is fixable within the scope of the manuscript, so I recommend major revision rather than rejection. I also suggest asking the authors to check the novelty statement against Refs. [30, 80, 81], since much of the scalar-induced and evaporation spectra builds directly on those works. I saw no indication of scholarly misconduct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the unburdened spectrum is a genuinely useful synthesis, and it is internally consistent. The memory-burden double-peak diagnostic advertised in the abstract, however, is algebraically wrong. I checked the stress-test note against the paper's own equations, and it is correct: Eq. (92) does not follow from Eqs. (85) and (90), and in the semiclassical limit n=0, q=1 it gives a ratio of 0 instead of 1.\n\nWhat is actually new: combining the four GW sources (inflationary PGW, inflaton scattering, PBH density-fluctuation GW, PBH evaporation GW) into a single template, with arbitrary w_phi and a memory-burden extension. The unburdened parts are solid. The matching at f_BH0 between Eqs. (41) and (45) works, the w_phi-dependent slopes look right, and the authors explicitly flag the two main approximations: the non-linear cutoff ambiguity (Appendix D.4, Fig. 9) and the Carr-Hawking gamma. The appendices are substantial, and the paper ships closed-form formulas that others can test.\n\nThe soft spots, in order. First, Eq. (92). Dividing Eq. (90) by Eq. (85) gives a different expression: 3^{-1/6} 2^{n/3}(3+2n)^{1/3} q^{2n/3}(M_in/M_P)^{4n/3}, not the 1.2 q^{1+4n/3} [2n(3+2n)]^{2/3}(M_in/M_P)^{4n/3-1} in the text. The n=0 limit of the published formula is 0, not 1. Since the double peak is the paper's main observational claim for memory burden, this is load-bearing. The paper needs to correct this and temper the 'direct access to n' claim. Second, the memory-burden ansatz itself, Eq. (12), has free q and n with no microscopic input, which the authors admit. That makes the burdened predictions model-dependent signatures, not evidence for the memory-burden mechanism. The semiclassical part does not suffer from this.\n\nWho should read it: anyone doing PBH reheating phenomenology or forecasting high-frequency GW detectors. The unburdened template is a useful contribution; the burdened diagnostic needs a fix before I would rely on it. It deserves a serious referee, and my expectation would be major revision, not rejection.","headline":"A genuinely useful unburdened GW template, but the burdened double-peak diagnostic is algebraically wrong and needs a major fix.","tokens_in":36733,"tokens_out":6598,"would_cite":true,"duration_ms":49358,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that a PBH-reheating universe writes a four-part gravitational-wave fossil record, readable band by band, and that a memory-burden double peak would expose the back-reaction.","keywords":["primordial black holes","gravitational wave background","PBH reheating","memory burden","Hawking evaporation","scalar-induced gravitational waves","inflaton scattering","stochastic gravitational wave background"],"falsifier":"Search for the spectrum with a detector sensitive to $\\Omega_{\\rm GW}h^2\\simeq10^{-18}$ across $10^{-5}$--$10^7$ Hz and check for the flat A1 plateau, the $f_0^{-2}$ slope in A2, the $f_0^{2/5}$ slope in A3, the B1/B2 peak at $f_{\\rm UV}^0\\simeq 4.0\\times10^6\\,(M_{\\rm BH}/1\\,{\\rm g})^{-5/6}$ Hz, and the evaporation bump near $f_{{\\rm ev}}^0\\simeq 1.5\\times10^{13}\\,(M_{\\rm BH}/1\\,{\\rm g})^{1/2}$ Hz. For the burden picture, look for a double evaporation peak with $f_{{\\rm ev}2}^0/f_{{\\rm ev}1}^0\\simeq 1.2\\, q^{1+4n/3}\\,[2n(3+2n)]^{2/3}(M_{\\rm in}/M_P)^{4n/3-1}$; a single peak would rule out the two-stage burdened evaporation, while a different ratio would fix $n$ at a value inconsistent with the power-law ansatz.","tokens_in":35662,"feed_emoji":"🌊","tokens_out":10745,"duration_ms":92974,"temperature":0.7,"pith_summary":"This paper argues that in a universe where primordial black holes (PBHs) reheat the cosmos after inflation, the stochastic gravitational-wave background observed today is the sum of four physically distinct sources, each leaving a recognizable imprint in a different frequency band. It derives analytic shapes and characteristic frequencies for each component: inflationary tensor modes that re-enter during radiation, PBH domination, or inflaton oscillations; gravitational waves from inflaton scattering at the highest frequencies; scalar-induced waves from the Poisson density fluctuations of the PBH population; and gravitons emitted directly by Hawking evaporation. It then redoes the calculation under the memory-burden hypothesis, in which evaporation is slowed once the hole has radiated a significant fraction of its mass, and shows that the burdened spectrum is pushed toward lower frequencies and acquires a double evaporation peak. If the picture is right, a future broadband gravitational-wave measurement would read out the reheating temperature, the PBH mass and abundance, the inflaton equation of state, and possibly the memory-burden parameters $n$ and $q$.","feed_headline":"GW spectrum from PBH reheating splits into four bands","feed_subtitle":"Each band fingerprints a different early-Universe era, from inflationary plateau to evaporation peak.","key_machinery":"The load-bearing machinery is the PBH evaporation law. Semiclassically $dM_{\\rm BH}/dt = -\\epsilon M_P^4/M_{\\rm BH}^2$; the memory-burden generalization divides this by $S^n$, where $S=\\frac12 (M_{\\rm BH}/M_P)^2$ is the black-hole entropy, introducing two free parameters: $q$, the mass fraction at which burdening begins, and $n$, the power of the suppression. This single ansatz fixes the reheating temperature and the duration of PBH domination, and through them every characteristic frequency in the spectrum. Together with transfer-function methods for tensor modes produced during inflation, graviton production by inflaton scattering, and scalar (Poisson) density fluctuations, it produces the region-by-region shapes A1 through D$_2$ in Fig. 1.","core_discovery":"The paper's central claim is a decomposition. With a PBH population that dominates and reheats the Universe, the present-day gravitational-wave abundance $\\Omega_{\\rm GW}h^2$ as a function of present frequency $f_0$ is assembled from region A1 (a flat plateau from modes entering during radiation), A2 (falling as $f_0^{-2}$ from modes entering during PBH domination), A3 (rising as $f_0^{(6w_\\phi-2)/(3w_\\phi+1)}$ from modes entering during inflaton oscillations), a sudden cutoff C from modes that never left the horizon, the peak B1/B2 from scalar-induced waves sourced by PBH number-density fluctuations, and a high-frequency evaporation bump D1. The paper verifies this architecture for potentials $V(\\phi)\\propto \\phi^k$ during reheating, with $k=6$ as the benchmark. Under memory burden it claims the same architecture survives but with reheating occurring later and at lower frequency: the gap ratio $f_{\\rm UV}^0/f_{\\rm IR}^0$ grows from about $4200$ to about $3.3\\times 10^6$ for $n=1$, $q=1/2$, and the evaporation bump splits into two peaks D$'_1$ and D$'_2$ whose frequency ratio $f_{{\\rm ev}2}^0/f_{{\\rm ev}1}^0 \\simeq 1.2\\, q^{1+4n/3}\\,[2n(3+2n)]^{2/3}\\,(M_{\\rm in}/M_P)^{4n/3-1}$ gives direct access to $n$.","pith_inferences":["The paper itself notes two simplifying choices: it assumes the burden effect begins during the PBH-dominated era and it neglects accretion. Either assumption would shift the predicted frequencies, so testing them is necessary before precise data comparisons can be drawn.","Beyond the paper: the double peak is the most falsifiable burden prediction, but its exact ratio depends on the power-law form $S^{-n}$; a two-peak signal with a different ratio would indicate a different back-reaction law rather than ruling out memory burden itself.","Beyond the paper: if the high-frequency evaporation peak is ever detected at the predicted amplitude, it could serve as a standard candle for the reheating epoch even if lower-frequency regions remain below detector sensitivity.","Beyond the paper: the same machinery could be applied to non-Poisson or strongly clustered PBH initial conditions, which would change the slope of the B1/B2 spectrum and could mimic or mask the memory-burden signature."],"forward_implications":["The flat A1 plateau amplitude, $\\Omega_{\\rm A1} h^2\\simeq 2.2\\times10^{-18}$, directly measures the inflationary Hubble scale $H_{\\rm end}$, independent of the PBH parameters.","The break frequencies $f_{\\rm RH}^0$ and $f_{\\rm BH}^0$ and the A3 slope encode the sequence of eras and the inflaton equation of state: for $V(\\phi)\\propto\\phi^6$, $\\Omega_{\\rm GW}\\propto f_0^{2/5}$ in region A3.","The density-fluctuation peak B1/B2 scales as $f_{\\rm UV}^0\\simeq 4.0\\times10^6\\,(M_{\\rm BH}/1\\,{\\rm g})^{-5/6}$ Hz and $f_{\\rm IR}^0\\simeq 947\\,(M_{\\rm BH}/1\\,{\\rm g})^{-3/2}$ Hz in the semiclassical case, so its position measures the PBH mass.","The evaporation bump D1 has an amplitude $\\simeq 7.5\\times10^{-7}$ that is independent of $M_{\\rm BH}$, making it a clean marker of a PBH reheating phase even when other features are obscured.","With memory burden, the whole spectrum is enlarged toward lower frequencies, and the double-peak ratio $f_{{\\rm ev}2}^0/f_{{\\rm ev}1}^0$ provides a direct measurement of the burden exponent $n$."],"supporting_citations":[{"why":"supplies the PBH-reheating scenario and the critical abundance $\\beta_c$ above which PBHs dominate the energy budget before evaporating.","marker":"[16]"},{"why":"supplies the memory-burden evaporation ansatz (suppression by a power of the black-hole entropy) that the burdened spectrum is built on.","marker":"[17]"},{"why":"justifies the $q=1/2$ benchmark and the $S^{-n}$ suppression form used for the numerical spectra.","marker":"[19]"},{"why":"provides the computation of gravitational waves from PBH evaporation that region D1 and the evaporation-peak amplitude are based on.","marker":"[30]"},{"why":"gives the information-storage argument that a black hole's own memory back-reacts to resist further information loss, motivating the burden effect.","marker":"[47]"},{"why":"supports treating inflaton-scattering gravitons as the high-frequency counterpart of modes that never left the horizon, which defines region C.","marker":"[63]"},{"why":"supplies the Poisson density-fluctuation formalism for the PBH distribution that sources the scalar-induced B1/B2 gravitational waves.","marker":"[76]"},{"why":"provides the burdened density-fluctuation gravitational-wave spectra that the paper generalizes to an arbitrary initial equation of state.","marker":"[80]"},{"why":"provides the transfer-function and frequency-ratio framework for scalar-induced gravitational waves during PBH reheating that the burdened B-region calculation extends.","marker":"[81]"}],"fun_headline_variants":["Four GW bands from PBH reheating fingerprint each era","PBH memory burden widens GW gap 800-fold","Memory burden splits PBH evaporation peak into two","GW spectrum decodes inflaton k and PBH memory burden","PBH reheating yields four-band GW signal, burden splits peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire burdened half of the paper rests on the ansatz that once a PBH has radiated a fraction $1-q$ of its mass, its evaporation rate is suppressed by an inverse power of its entropy, $S^{-n}$, with $q$ and $n$ free parameters and benchmark $q=1/2$, $n=1$; the paper itself states that there are no real microscopic motivated values for $q$.","fun_headline_variants_meta":{"raw":{"variants":["Four GW bands from PBH reheating fingerprint each era","PBH memory burden widens GW gap 800-fold","Memory burden splits PBH evaporation peak into two","GW spectrum decodes inflaton k and PBH memory burden","PBH reheating yields four-band GW signal, burden splits peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000295,"raw_usage":{"total_tokens":1743,"prompt_tokens":999,"completion_tokens":744,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":662}},"tokens_in":615,"tokens_out":744,"duration_ms":7642,"temperature":1.0,"reasoning_tokens":662,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:36:08.167815+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Search for the spectrum with a detector sensitive to $\\Omega_{\\rm GW}h^2\\simeq10^{-18}$ across $10^{-5}$--$10^7$ Hz and check for the flat A1 plateau, the $f_0^{-2}$ slope in A2, the $f_0^{2/5}$ slope in A3, the B1/B2 peak at $f_{\\rm UV}^0\\simeq 4.0\\times10^6\\,(M_{\\rm BH}/1\\,{\\rm g})^{-5/6}$ Hz, and the evaporation bump near $f_{{\\rm ev}}^0\\simeq 1.5\\times10^{13}\\,(M_{\\rm BH}/1\\,{\\rm g})^{1/2}$ Hz. For the burden picture, look for a double evaporation peak with $f_{{\\rm ev}2}^0/f_{{\\rm ev}1}^0\\simeq 1.2\\, q^{1+4n/3}\\,[2n(3+2n)]^{2/3}(M_{\\rm in}/M_P)^{4n/3-1}$; a single peak would rule out the two-stage burdened evaporation, while a different ratio would fix $n$ at a value inconsistent with the power-law ansatz.","supporting_citations":[],"review_version":2}