{"id":"66392326-eecc-4d4d-b718-f48cb2e251a3","arxiv_id":"2411.17658","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"PSpectCosmo applies a pseudo-spectral method to inflationary lattice simulations and claims to resolve the divergent vacuum energy density by tracking a quantum-to-classical ratio R.","lead":"PSpectCosmo is a new C++ lattice simulation code that uses Fourier-space pseudo-spectral derivatives to evolve scalar fields during inflation and reheating. The paper argues this avoids a wave-vector mismatch of finite difference methods and introduces a quantum-to-classical ratio R to tame the divergent vacuum energy at the start of inflation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper defines a quantumness ratio R in Eq. (3.22) but never states the threshold or mode-selection rule that maps R to a convergent energy density, leaving the abstract's central 'resolving divergence' claim undemonstrated and untested by the Fig. 2 power-spectrum check.","rationale":"The reader's rejection is well supported, and I agree with the identification of the weakest assumption. The strongest claim has two parts: accurate power spectra and convergent energy density. The first part is reasonably supported by Fig. 2, by the standard pseudo-spectral construction, and by the public code; I do not object to the spectral method itself. The second part is the load-bearing gap. Section 3.3 introduces R as a diagnostic but gives no concrete mapping from R to ρ. A threshold-stability test is the right way to settle the matter because it directly probes whether a unique, convergent energy density exists. If the code contains a threshold, varying it reveals sensitivity; if no cutoff exists, the advertised mechanism is absent. The current paper therefore cannot support its central claim, and the verdict should remain as the reader set it.","tokens_in":7234,"tokens_out":8547,"duration_ms":81962,"concrete_test":"Open the public PSpectCosmo repository, locate the energy-density routine that uses R (e.g., search for 'R' or 'classical' in the source), and extract the exact threshold R_c and the mode-selection/filtering rule. Then rerun the ϕ²-model configuration of Fig. 2 with two thresholds, R_c = 0.1 and 0.5, at N=128 and N=256, and output the volume-averaged energy density and the Hubble rate from Eq. (3.3). If ρ changes by more than the numerical resolution error when R_c changes, or if no R-dependent cutoff exists in the code, the 'convergent energy density' claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's headline capability is a 'robust mechanism to compute convergent energy density.' Section 3.3 motivates this by the UV-divergent vacuum energy, then defines R in Eq. (3.22) and asserts that R distinguishes quantum (R ∼ 1) from classical (R ≪ 1) modes. The manuscript never states how R is converted into an energy-density calculation: there is no threshold R_c, no formula for a mode cutoff or weight, and no algorithm tying Eq. (3.22) to Eq. (3.15). It also does not show that the resulting ρ is finite, resolution-independent, or consistent with the Friedmann equation (3.3). Figure 1 is an unlabeled plot with no threshold marked. The only quantitative validation, Fig. 2, tests the scalar power spectrum against the linear theory used to set initial conditions in §3.2; it does not exercise the R-based energy subtraction. Because the prescription is unspecified, the central claim cannot be checked, and there is no way to tell whether the advertised convergence is a physical renormalization or a tuned subtraction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents PSpectCosmo, a publicly available C++ pseudo-spectral lattice code for simulating scalar-field dynamics during inflation and reheating. The two advertised advances are (i) the use of a Fourier-space pseudo-spectral method so that the effective wave vector equals the lattice Fourier wave vector, avoiding known finite-difference artifacts, and (ii) a mechanism, based on a ratio R, to separate quantum from classical modes so that the energy density can be computed without the ultraviolet vacuum divergence. The paper derives the equations of motion, describes the lattice initialization from Bunch-Davies vacuum mode functions, defines R in Eq. (3.22), and shows two figures: one unlabeled plot of R and one power-spectrum comparison for a ϕ² model. The central claims of precise non-linear evolution and convergent energy density are asserted in the abstract and summary, but the manuscript does not specify the operational algorithm for the energy cut-off or provide non-linear validation.","tokens_in":7498,"tokens_out":6800,"duration_ms":63788,"significance":"If the claims were fully established, PSpectCosmo would be a useful open-source tool for studying non-linear effects during inflation and preheating, a topic of current interest for primordial black holes and gravitational waves. The pseudo-spectral approach is well motivated, and the public availability of the code, together with the use of standard output formats, is a strength. However, as written, the paper does not demonstrate the central advertised capabilities: the R-based energy-density prescription is not defined, the only quantitative test is a linear self-consistency check, and no interacting-field or gravitational-wave validation is shown. The paper is better characterized as a brief code description than as a validated methods paper.","major_comments":[{"comment":"The initial conditions for the perturbations are written as the deterministic mode functions f_k/a with no random phases. A realistic lattice simulation of a quantum field requires populating each Fourier mode as f_k e^{iθ_k}, with θ_{-k} = -θ_k for a real field, to represent a stochastic realization of the vacuum state. Without this, the initial field is a single coherent configuration whose power spectrum is |f_k|^2 by construction, so the agreement in Fig. 2 reflects the initialization rather than the accuracy of the integrator. The authors should specify the phase sampling procedure and indicate whether the reported spectra are ensemble-averaged.","section":"§3.2, Eqs. (3.11)–(3.12)"},{"comment":"The abstract's claim of a 'robust mechanism to compute convergent energy density' is not supported because the manuscript never specifies how R is used. There is no threshold value R_c, no formula for excluding or weighting modes in the energy sum of Eqs. (3.15)–(3.21), and no demonstration that the resulting ρ is finite, resolution-independent, or consistent with the Friedmann equation (3.3). The text asserts that R∼1 marks quantum modes and R≪1 marks classical modes, but it does not give the operational rule that maps R onto a mode-cut or weight. The authors should state the algorithm completely and test it by varying the threshold and lattice resolution.","section":"§3.3, Eq. (3.22)"},{"comment":"The power-spectrum validation is a linear self-consistency check: the lattice is initialized with the same Bunch-Davies mode functions whose theoretical power spectrum is used as the 'prediction.' This verifies that the integrator reproduces the linear evolution used to set the initial conditions, but it does not exercise the code's advertised capabilities for interacting scalar fields or gravitational waves. A validation against a known non-linear benchmark, or a comparison with an independent lattice code for a self-interacting model, is needed to support the Introduction and Abstract claims.","section":"§5, Fig. 2"},{"comment":"The paper omits the numerical setting for the reported results: lattice size N, box length, physical parameters of the ϕ² model, the adaptive time-step tolerance, and the values of any thresholds used for the energy-density computation. These details are necessary for the reader to reproduce Fig. 2 and to assess the claimed 'exceptional accuracy.' A convergence study (e.g., power spectrum versus N) would substantiate the pseudo-spectral advantage over finite differences.","section":"§4 and §5"}],"minor_comments":[{"comment":"The axes of Fig. 1 are unlabeled, and no threshold value is marked; the reader cannot determine what is plotted on the x-axis or how R evolves without additional information.","section":"Fig. 1"},{"comment":"The caption reads 'when the perturbation exist horizon'; this should be 'when the perturbations have exited the horizon.'","section":"Fig. 2 caption"},{"comment":"The notation Δf and Δπ is not defined precisely; these should be the standard deviations of the mode operators (or their uncertainties), rather than the operators themselves.","section":"Eq. (3.22)"},{"comment":"The statement that the pseudo-spectral method 'can accommodate both periodic and non-periodic boundary conditions' is misleading for the Fourier method described, which is inherently periodic; non-periodic domains require a different spectral basis such as Chebyshev polynomials.","section":"§2"},{"comment":"The paper states natural units with ℏ=1, but Eq. (3.22) retains ℏ explicitly; the unit convention used for the code's internal variables and for the reported power spectrum should be clarified.","section":"§4"},{"comment":"The R-based method is attributed to Baumann's textbook [23], but no specific section or equation is cited; a precise pointer would help the reader locate the underlying argument.","section":"Reference [23]"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is not ready for publication in its current form: the headline capability of convergent energy density is not actually defined in the text, and the only numerical test is circular. However, the pseudo-spectral core is sound, and the missing pieces—the R-threshold algorithm, stochastic initial conditions, and a non-linear benchmark—are additions rather than fundamental errors. I therefore recommend major revision rather than rejection, with the expectation that the authors provide the operational details and substantive validation. The editor may also wish to consider whether the paper's length and depth are sufficient for JCAP."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jordan, here's my take on Jiang's PSpectCosmo paper. The code is real and public, and the pseudo-spectral machinery is standard and competently explained. The power spectrum test in Fig. 2 matches the linear theory, but that is mostly a consistency check: the initial conditions come from the same Bunch-Davis mode functions. It does not exercise the energy-density claim.\n\nThe load-bearing gap is in Sec. 3.3. Eq. (3.22) defines R as the uncertainty ratio, with R~1 quantum and R<<1 classical, but the paper never says how R is converted into a cutoff or weight for the energy density. There is no threshold value, no formula connecting R to Eq. (3.15), and no description of an algorithm. Fig. 1 is an unlabeled curve with no threshold marked. So the abstract's promise of a 'robust mechanism to compute convergent energy density' is not actually demonstrated. You cannot tell whether the subtraction is physical or tuned.\n\nThere is also no independent validation: no convergence test in lattice size, no check of the energy density against an analytic limit, and no test that the resulting rho self-consistently drives the scale factor via the Friedmann equation. Given that PSpectRe and CosmoLattice already use pseudo-spectral methods, the only genuinely new element is the R-based cutoff, and that is exactly the part that is not pinned down.\n\nCredit where due: the code is on GitHub, the paper is clearly written, and the citations are appropriate, including the textbook source for R. A resourceful referee could help the author specify and test the R prescription. If that happens, this could become a modest but useful code paper. As it stands, the central claim is unsupported.\n\nMy recommendation: I would not accept it as is. If it crossed my desk, I'd send it to a serious referee anyway, because the underlying code is real and the gap might be fixable. But I'd expect heavy revision, and I'd want to see the exact mapping from R to the energy density, plus at least one non-trivial validation, before I'd endorse publication.","headline":"Fresh code, clean writing, but the central energy-density claim is a black box: R is defined, yet no threshold or algorithm is stated, so the headline result is unvalidated.","tokens_in":7954,"tokens_out":2969,"would_cite":false,"duration_ms":28141,"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":"PSpectCosmo claims that a Fourier-space pseudo-spectral lattice code preserves the inflationary power spectrum and yields a finite energy density through a quantum-classical ratio R.","keywords":["pseudo-spectral method","lattice simulation","inflation","reheating","power spectrum","quantum-to-classical transition","energy density","gravitational waves"],"falsifier":"Run one inflationary model with the $R$ cutoff at several distinct values, such as $R = 0.5$, $R = 0.1$, and $R = 0.01$, and compare the total energy density and the resulting scale-factor history; if these quantities change materially between runs, the claimed convergence is a calibration artifact rather than a physical result. A second check is to double the lattice resolution $N$ with the physical box fixed and confirm that the power spectrum and energy density are unchanged.","tokens_in":7038,"feed_emoji":"🌌","tokens_out":10698,"duration_ms":86688,"temperature":0.7,"pith_summary":"PSpectCosmo is a C++ lattice code designed to follow scalar-field dynamics from inflation through reheating. The paper's central claim is that replacing finite-difference spatial derivatives with Fourier-space pseudo-spectral derivatives removes the mismatch between effective and lattice wave vectors, so the simulated inflationary power spectrum agrees with the theoretical prediction all the way to the highest resolvable frequency. It further claims that a mode-by-mode ratio $R = \\hbar/2\\,/\\,(|f||\\pi|)$ distinguishes quantum vacuum fluctuations from classical fluctuations, allowing the code to exclude divergent sub-horizon vacuum energy at the start of inflation and include the energy of modes after horizon exit, yielding a convergent total energy density that can drive the background expansion. If these claims hold, the code offers a way to study non-linear inflationary and reheating processes, such as primordial black hole formation and gravitational wave production, on a lattice without the artifacts of finite-difference methods.","feed_headline":"Pseudo-spectral code makes inflation's energy density finite","feed_subtitle":"Fourier-space derivatives plus a quantum-classical ratio R give finite energy density and a matching power spectrum.","key_machinery":"The paper rests on two mechanisms. The first is the Fourier-space pseudo-spectral derivative: the discrete Fourier transform of the field, multiplication by $(ik)^\\nu$ in Fourier space, and an inverse transform, which gives exponential convergence for smooth fields and, crucially, makes the effective derivative wave vector identical to the lattice wave vector $k_{\\mathrm{Lat}}$. The second is the quantum-classical ratio $R \\equiv \\hbar/2\\,/\\,(|f||\\pi|)$ of Eq. (3.22), built from the mode function $f \\equiv a\\,\\delta\\phi$ and its conjugate momentum $\\pi$; it marks whether a lattice mode is still a quantum vacuum fluctuation ($R \\sim 1$) or has become classical ($R \\ll 1$) and therefore whether its energy belongs in the total energy density. Time integration is handled by an adaptive-step second-order symplectic integrator, and the field and energy are evaluated in Fourier space so that the discrete Fourier identity relating position-space and frequency-space sums of squares gives the volume averages used in the background expansion equation.","core_discovery":"The paper's central claim is that the pseudo-spectral method resolves two obstacles that have kept lattice simulations from covering inflation itself. First, because the derivative is applied as multiplication by $(ik)^\\nu$ in Fourier space, the wave vector used in the dynamics is exactly the lattice wave vector $k_{\\mathrm{Lat}}$; the finite-difference operator's effective wave vector $k_{\\mathrm{eff}}$ no longer enters, so the scalar power spectrum is not distorted and gravitational wave mode-mode coupling is represented faithfully. Second, the energy density is computed in Fourier space with the $k = 0$ background separated, and the ratio $R \\equiv \\hbar/2\\,/\\,(|f||\\pi|)$ is used as a diagnostic of the quantum-to-classical transition: modes with $R \\sim 1$ are quantum vacuum fluctuations whose energy is excluded, while modes with $R \\ll 1$ have exited the horizon and are included. This yields an energy density that is finite at the onset of inflation and converges as inflation proceeds, allowing the scale factor to be evolved self-consistently. The paper validates the method with a $\\phi^2$ model, showing a nearly scale-invariant power spectrum in agreement with theory up to the highest resolvable frequency.","pith_inferences":["The paper does not state a numerical threshold for $R$ at which a mode switches from excluded to included; testing whether the final energy density and power spectrum are insensitive to that threshold would establish whether the convergence claim is a physical result or a calibration choice.","In a genuinely nonlinear simulation, mode coupling transfers power between scales, so the free-field-based $R$ may mislabel modes as classical soon after they leave the horizon; one could compare $R$ against a separate criterion such as $k/(aH)$ or the mode occupancy to see whether the classification is stable.","If the convergent-energy construction works as stated, it would allow lattice codes to compute the backreaction of sub-horizon fluctuations on the background expansion during inflation without imposing linear theory, a regime the paper only begins to explore."],"forward_implications":["In the $\\phi^2$ model the simulated dimensionless power spectrum is nearly scale invariant and matches the theoretical prediction for all wave numbers below the highest resolvable frequency, so the pseudo-spectral derivative does not introduce the finite-difference distortion of the scalar spectrum.","Because the Fourier derivative gives $k_{\\mathrm{eff}} = k_{\\mathrm{Lat}}$, the gravitational-wave source convolution is computed with the true lattice wave vector, removing the ambiguity finite differences create for mode-mode coupling.","The ratio $R = \\hbar/2\\,/\\,(|f||\\pi|)$ lets the code exclude quantum vacuum modes from the energy density at early times and include them after horizon exit, so the total energy density stays finite and can drive the background expansion through inflation into reheating.","The code adopts the output and power-spectrum conventions of an existing lattice package, so its results can be compared directly with other early-universe simulations."],"supporting_citations":[{"why":"Identifies the finite-difference mismatch between effective and lattice wave vectors and shows it corrupts gravitational-wave mode coupling.","marker":"[10]"},{"why":"Supplies the lattice-variable rescaling strategy and the discussion of derivative-operator consistency the code adopts.","marker":"[13]"},{"why":"Provides the earlier pseudo-spectral code whose approach this work extends from preheating to inflation.","marker":"[15]"},{"why":"Documents that the wave-vector mismatch distorts the scalar power spectrum in inflationary lattice simulations.","marker":"[21]"},{"why":"Establishes the exponential-convergence property of spectral differentiation for smooth solutions.","marker":"[22]"},{"why":"Outlines the quantum-to-classical mode criterion used to define the ratio $R$.","marker":"[23]"},{"why":"Defines the output and power-spectrum conventions the code follows for comparability with other simulations.","marker":"[14]"}],"fun_headline_variants":["Pseudo-spectral code tames inflation's energy density","Fourier method yields finite energy in inflation runs","PSpectCosmo: Inflation simulations with convergent energy","New code solves inflation's divergent energy problem"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the ratio $R = \\hbar/2\\,/\\,(|f||\\pi|)$ being a valid and threshold-stable classifier of quantum versus classical lattice modes, but the paper does not state at what value of $R$ a mode switches from excluded to included or how that cutoff is applied, so if nonlinear mode coupling shifts $R$ for partly quantum modes, or the cutoff is tuned to force a desired energy, the convergent energy density claim would fail.","fun_headline_variants_meta":{"raw":{"variants":["Pseudo-spectral code tames inflation's energy density","Fourier method yields finite energy in inflation runs","PSpectCosmo: Inflation simulations with convergent energy","New code solves inflation's divergent energy problem"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000655,"raw_usage":{"total_tokens":3022,"prompt_tokens":990,"completion_tokens":2032,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":606,"completion_tokens_details":{"reasoning_tokens":1970}},"tokens_in":606,"tokens_out":2032,"duration_ms":13112,"temperature":1.0,"reasoning_tokens":1970,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:51:40.883136+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run one inflationary model with the $R$ cutoff at several distinct values, such as $R = 0.5$, $R = 0.1$, and $R = 0.01$, and compare the total energy density and the resulting scale-factor history; if these quantities change materially between runs, the claimed convergence is a calibration artifact rather than a physical result. A second check is to double the lattice resolution $N$ with the physical box fixed and confirm that the power spectrum and energy density are unchanged.","supporting_citations":[],"review_version":1}