{"id":"332c1e23-fbd0-4a92-9ac8-3b0a4d8413f2","arxiv_id":"2608.02798","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"In a loop quantum cosmology model with negative cosmological constant, semiclassical states built by integrating over all superselection sectors behave no differently from single-sector states, with differences below numerical precision.","lead":"This paper tests whether an alternative way to build a separable Hilbert space in loop quantum cosmology, averaging over all superselection sectors, changes the evolution of semiclassical universes. On a bouncing and recollapsing model universe, it finds no detectable difference in volume dynamics or decoherence, and confirms a standard numerical library can replace a problem-specific code.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Semiclassicality is probed only through volume mean/variance; higher moments, the Hubble-rate observable, or non-Gaussian states could still distinguish integral from single-sector states.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the representative-ness of the volume mean/variance diagnostics and the Gaussian state family. My stress-test agrees and sharpens it: the integral-state variance formula (22c) contains a cross-sector term, so the paper's check is specifically a check on the first two moments of one observable; it does not certify full phase-space semiclassicality. This is not an internal inconsistency or a numerical error; it is an evidence gap for the breadth of the conclusion. The authors' own Sec. V caveat about model-specificity and continuous spectra reinforces that the claim should not be over-generalized. The two independent numerical methods and the reported eigenvalue agreement to 10^{-14} are genuine supporting evidence, so the result is credible as far as it goes. The appropriate verdict remains CONDITIONAL: the paper should either broaden its diagnostics (higher moments, b-observable, non-Gaussian states) or restrict its claim to volume expectation and variance. Since the reader already assigned CONDITIONAL, the verdict is unchanged.","tokens_in":19830,"tokens_out":11723,"duration_ms":123688,"concrete_test":"Recompute the Sec. IVB comparison using the full volume probability distribution P(v, phi) (or equivalently the third and fourth cumulants of V) and the variance of the Hubble-rate observable b, for the same states of Eq. (20) and additionally for a non-Gaussian spectral profile (e.g., a squeezed Gaussian with a linear phase in the index j). If, at late times (phi = 160 G^{-1/2}) and for any state in this extended grid, the integral-state P(v) or the b-variance deviates from the single-sector ensemble average by more than the numerical tolerance of 10^{-9}, then the volume mean/variance test missed a real sector-induced loss of semiclassicality.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. V: 'Within the precision of the methods used the integral states preserve the properties of the single sector ones') is inferred entirely from the first two moments of a single observable, the volume, evaluated for Gaussian spectral profiles via Eqs. (21)-(22). This is weaker than the claim for two reasons. First, the integral variance (22c) contains a sector-variance term, and the paper only bounds it indirectly through the relative difference delta_sigma_r of Eq. (23); a sector-induced broadening that leaves the volume mean and variance unchanged would not be detected. Second, volume localization does not by itself certify phase-space semiclassicality: the momentum-conjugate observable b, or the shape of the volume probability distribution, could in principle show broadening that the first two moments miss. The parameter scan (Sec. IVA) covers only omega_o in {1000,4000} G^{1/2}, sigma_omega in {50/sqrt(2),200/sqrt(2)} G^{1/2}, and Lambda in {-0.2,-0.1,-0.05,-0.01} l_Pl^{-2}, all with Gaussian profiles. Non-Gaussian or broader states, and higher cumulants of V, are not probed. The authors themselves caution in Sec. V that the result is model-specific and may fail for continuous spectra; the same caution applies to the narrow observable set used to certify 'no increased decoherence'.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the integral superselection-sector Hilbert space construction of Barbero et al. in the concrete LQC model of a flat isotropic universe with a negative cosmological constant and a massless scalar field. It computes the spectrum of the evolution operator Θ_Λ by two independent numerical methods—the purpose-built \"shots\" method of earlier work and a direct application of the Eigen library—and compares the volume expectation values and variances of Gaussian semiclassical states evolved over many bounce–recollapse cycles, in single-sector and integral-sector constructions. The main reported finding is that, within numerical precision, the integral states preserve the volume semiclassicality properties of single-sector states and show no observable increase of decoherence; a secondary finding is that the standard eigensolver reproduces and in most cases improves on the dedicated shots method, at the price of significantly larger computational resources.","tokens_in":20125,"tokens_out":13238,"duration_ms":128033,"significance":"If correct, the central numerical result provides direct evidence that the integral Hilbert space can serve as a practical replacement for single-sector projections in models where the latter are intractable, while the method comparison validates the older shots code and demonstrates that standard numerical libraries are usable for this class of problems. The paper is a numerical case study rather than a theorem; its strengths are the use of two independent implementations, the precision-escalation check (double vs long double), the explicit parameter scan over Λ, ω_o, σ_ω and 64 sectors, and the clear statement of model-specific limitations in Sec. V. The main claims are therefore plausible but rest on a fairly narrow set of observables and states, which is the central point needing attention.","major_comments":[{"comment":"The manuscript never states whether the discrete spectral coefficients obtained by sampling the Gaussian profile (20) are renormalized in each superselection sector. Because the spectra of Θ_Λ differ across ϵ, the sums ∑_j |Ψ̃_ϵ,j|² are not automatically equal, and if the sector states are not unit-normalized the objects in (21) are not expectation values; the comparison of δV_r and δσ_r in (23) would then mix normalization fluctuations with physical sector differences. Please specify the normalization procedure and confirm that every sector state entering the integrals in (22) has unit norm.","section":"Sec. IVA, Eqs. (19)–(22)"},{"comment":"The conclusion that \"there is no noticeable increase of decoherence due to differences between sectors\" is drawn from the first two moments of a single observable, the volume, for Gaussian spectral profiles with ω_o ∈ {1000, 4000} G^{1/2}, σ_ω ∈ {50/√2, 200/√2} G^{1/2}, and Λ ∈ {-0.2, -0.1, -0.05, -0.01} l_Pl^{-2}. This diagnostic cannot exclude sector-induced broadening of the momentum-conjugate observable b, growth of higher cumulants of V, or qualitatively different behavior for non-Gaussian or much wider profiles. Since the abstract phrases the result as a general statement about \"long term semiclassicality properties,\" the claim should either be restricted to the volume-moment criterion used here or supported by additional probes.","section":"Sec. IVB and Sec. V"},{"comment":"The independent validation of the eigenlibrary basis against the shots method is shown for the first 3000 eigenvalues in three sectors, but the most demanding dynamical runs (ω_o = 4000, Λ = -0.05) require about 3400 eigenvalues per sector (Sec. IIID). The accuracy of eigenpairs 3001–3400 is therefore not directly benchmarked, and since truncation and boundary errors typically increase with eigenvalue, this is precisely the range where the \"lowest 10%\" assumption is least safe. Please extend the comparison to the full set of eigenvalues used in the dynamics or provide an explicit error estimate for that range.","section":"Sec. IIIC, Sec. IIID, and Fig. 9"}],"minor_comments":[{"comment":"The text states that the evolution is performed within the range ϕ ∈ [0,160] G^{-1/2}, but all figures and the discussion of the 120th cycle correspond to ϕ up to 60000 G^{-1/2} (e.g., Fig. 3); this inconsistency should be corrected.","section":"Sec. IVB, point 1"},{"comment":"The abstract describes the Eigen-based method as \"more efficient\" than the shots method, but Sec. IIID and the Conclusions state that the new method requires at least an order of magnitude more computational resources; recommend rewording to \"easier to apply\" or \"more accurate in most cases\".","section":"Abstract and Sec. IIID/Conclusions"},{"comment":"The symbol δσ_r is used both for the sector-resolved quantity δσ_r(ϵ,ϕ) in (23) and for the sector-averaged quantity δσr(ϕ) in (24); please use distinct notation to avoid confusion.","section":"Eqs. (23) and (24)"},{"comment":"The word \"bissection\" should be \"bisection\" in the two places where it appears.","section":"Sec. IIIB"},{"comment":"The word \"relatiely\" should be \"relatively\".","section":"Sec. V"},{"comment":"Reference [24] is missing the article number and page range; it should be Phys. Rev. D 77, 124025 (2008).","section":"References"},{"comment":"The captions use \"(+)\" and \"(−)\" without explaining that these denote the positive and negative branches of the plotted quantity; the captions should be self-contained.","section":"Figs. 6 and 7 captions"},{"comment":"The name \"Friedman-Lemaitre-Robertson-Walker\" should be \"Friedmann-Lemaître-Robertson-Walker,\" and the inconsistent use of \"FRLW\" versus \"FLRW\" should be unified.","section":"Sec. II and Introduction"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the authors are appropriately cautious in Sec. V. The main concern is not the correctness of the numerics but the gap between the broad wording of the abstract and the narrow diagnostic set; if the authors tighten the wording or add the requested checks, I would support publication. The abstract's characterization of the Eigen-based method as \"more efficient\" should also be corrected, since the body reports it as more resource-intensive."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jan, here is my read of arXiv:2608.02798. The paper does something useful and concrete: it tests whether the integral Hilbert space construction, which is proposed as a way around non-separable sector spaces in LQC, preserves semiclassicality in a specific isotropic model with negative cosmological constant. This is the first explicit dynamical comparison of integral versus single-sector states in this setting, and it also benchmarks a standard numerical library (Eigen) against the problem-specific 'shots' method. On those terms it does well. The numerics are careful: two independent methods agree on eigenvalues to 1e-14, a double-versus-long-double precision test shows that the slow drift in differences is numerical error, and the averaged variance difference shows no trend over 120 cycles. The paper is also honest about its limits, explicitly noting that the result is model-specific and may not extend to spectra that are already continuous. The soft spots are real but proportionate. The stress-test note is right: the paper infers 'no increased decoherence' from the first two moments of a single observable, the volume, for Gaussian spectral profiles only. Higher moments of V, or the conjugate Hubble-rate observable, could in principle distinguish integral from single-sector states, and the parameter scan is quite thin (two peak frequencies, two widths, four values of Lambda). So the actual result is: within numerical precision, volume expectation and variance of these Gaussian states do not show sector-induced broadening over 120 cycles. That is a solid first test, not the general statement the abstract seems to promise. The authors do hedge in the conclusions, but the abstract and Section V could be more precise. I would also have liked the code and data to be released; for a numerical result, this is not a luxury, and it would let others check the truncation and extrapolation choices. Who is this for? LQC people working on superselection sectors, especially those planning to apply the integral construction to Bianchi I or models where single-sector techniques fail. For that audience it deserves a serious referee. I would send it to peer review and ask for a tightened statement of the claim, a slightly broader observable set or at least an explicit acknowledgment of the limitation, and ideally a public code repository. It is not a transformative paper, but it is a competent and useful one, and the numerics look trustworthy within their stated scope.","headline":"A careful, honest numerical test of the integral Hilbert space construction in one LQC model; the central claim holds for what is actually probed, but the headline conclusion is narrower than the paper's wording suggests.","tokens_in":694,"tokens_out":709,"would_cite":false,"duration_ms":25409,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.60.Pp","98.80.Qc"],"model":"deepseek-v4-flash","headline":"This paper claims that integral superselection-sector states preserve the semiclassical properties of single-sector states, within numerical precision, making the integral Hilbert space a viable replacement in this model.","keywords":["loop quantum cosmology","polymer quantization","superselection sectors","integral Hilbert space","semiclassical states","volume operator","decoherence","numerical spectra"],"falsifier":"Compute the volume trajectories with a basis accurate to better than $10^{-15}$ over the full range of $\\epsilon$ and $\\phi$, using a non-Gaussian spectral profile or an observable beyond the volume; if the sector-averaged variance then grows monotonically with $\\phi$ or the relative difference $\\delta\\sigma_r$ exceeds $10^{-9}$ away from recollapse points, the paper's claim of no added decoherence would be refuted.","tokens_in":19644,"feed_emoji":"🌌","tokens_out":10278,"duration_ms":79589,"temperature":0.7,"pith_summary":"Loop quantum cosmology's quantization of geometry produces Hilbert spaces that are nonseparable, so calculations typically restrict to a single superselection sector. This paper tests an alternative: a separable 'integral' Hilbert space built by averaging over all sectors with the Lebesgue measure, an approach needed for models like Bianchi I where no single sector is viable. Using a quasi-cyclic isotropic universe with negative cosmological constant and a massless scalar field, it builds Gaussian semiclassical states in both approaches and evolves the volume operator over many bounce-recollapse cycles. The central claim is that, within numerical precision, the integral states behave like single-sector states: relative differences stay below $10^{-9}$, and the apparent tendency to decohere disappears when the numerical precision is increased. If right, the construction lets standard Hilbert-space and numerical tools be applied where sector-by-sector analysis fails.","feed_headline":"Sector averages preserve quantum cosmology's semiclassicality","feed_subtitle":"Loop quantum cosmology's integral Hilbert space matches single-sector dynamics with no added decoherence.","key_machinery":"The load-bearing object is the LQC evolution operator $\\Theta_\\Lambda = \\Theta_o - \\pi G \\gamma^2 \\Delta \\Lambda v^2\\,\\mathbb{1}$, a second-order difference operator on each superselection lattice $\\mathcal{L}_\\epsilon = \\epsilon + 4\\mathbb{Z}$. Its non-negative spectrum $\\{\\omega_j^2\\}$ supplies the energy basis $e_{\\epsilon,j}$ used in the Gaussian spectral profiles. The integral Hilbert space construction $\\mathcal{H}_{\\mathrm{int}} = \\frac{1}{4}\\int_0^4 \\mathcal{H}_\\epsilon\\, d\\epsilon$ is what makes the nonseparable kinematical space separable without choosing a sector. The paper's test proceeds by comparing the volume expectation $\\langle V\\rangle$ and variance $\\Delta V$ computed from single-sector states versus sector-integrated states, using two numerical routes to the same spectrum: the problem-specific shots method and a standard tridiagonal eigensolver.","core_discovery":"On the paper's own terms, the discovery is that integrating over superselection sectors does not spoil the semiclassical dynamics of the example universe. The integral Hilbert space $\\mathcal{H}_{\\mathrm{int}} = \\frac{1}{4}\\int_0^4 \\mathcal{H}_\\epsilon\\, d\\epsilon$ is separable and carries a well-defined action of observables through their restrictions to each sector. For the flat isotropic model with $\\Lambda<0$ and a massless scalar, the authors compute the spectra of the evolution operator $\\Theta_\\Lambda$ on 64 sectors using two independent numerical methods, form Gaussian spectral profiles peaked at $\\omega_0 = 10^3$ and $4\\cdot 10^3$ (in $G^{1/2}$) with two variances, for four values of $\\Lambda$, and evolve the volume expectation value and variance through $\\phi \\in [0,160]\\,G^{-1/2}$. They report that relative differences between integral and single-sector trajectories stay below $10^{-9}$ throughout, that the slow growth of the apparent difference falls by several orders of magnitude when the basis is computed in higher precision, and that the averaged variance difference oscillates around zero rather than trending, so there is no detectable extra decoherence. The two methods' spectra agree to relative differences near $10^{-14}$, and the standard eigensolver is at least as accurate while being much easier to apply, at the cost of an order of magnitude more computation.","pith_inferences":["A natural extension would be to probe semiclassicality with observables other than the volume operator, such as the Hubble parameter or bounce-time moments, and with non-Gaussian spectral profiles; if those also show no sector-induced differences, the conclusion would be much stronger.","The paper leaves open whether the integral construction remains well-behaved when the individual sectors have continuous spectra; one concrete test would be an anisotropic model like Bianchi I where this is the actual obstacle.","The recollapse-point spikes in the shots method suggest that tail extrapolation, not the physics, is the limiting error source; replacing the exponential tail extrapolation with the exact decay of the corresponding continuum-geometry eigenfunctions would likely eliminate them and make the dedicated method competitive in accuracy at high $\\omega_0$.","The success of the integral construction here hints that superselection sectors are not physically distinguished: the dynamics averaged over them is the same, so observables may be definable without ever choosing a sector."],"forward_implications":["In this model, the integral Hilbert space can replace single-sector projections: the volume dynamics of semiclassical states is unchanged within numerical precision, so no sector selection is needed.","The standard eigensolver reproduces the dedicated shots method's spectra to about $10^{-14}$ relative agreement, so future studies can use generic numerical libraries rather than problem-specific code.","The apparent small growth of differences between integral and single-sector states over time is numerical, not physical; increasing basis precision reduces it by several orders of magnitude.","Because the eigenvalue curves $\\omega_j(\\epsilon)$ vary continuously with $\\epsilon$, the integral Hilbert space supports a continuous range of energies, whereas each sector alone has a discrete spectrum.","The paper stresses that the result is model-specific: extending the conclusion to models whose sectors already have continuous spectra or non-trivial measure spaces is not automatic."],"supporting_citations":[{"why":"Supplies the model's quantization, the shots method, and the single-sector semiclassical dynamics that this work extends and cross-checks.","marker":"[24]"},{"why":"Introduces the integral Hilbert space construction over superselection sectors, the approach tested here.","marker":"[22]"},{"why":"Defines the mainstream LQC quantization, the area gap $\\Delta$, and the form of the evolution operator $\\Theta_\\Lambda$.","marker":"[27]"},{"why":"Provides the lapse choice $N=2V$ and regularization adjustments that modify the model relative to [24].","marker":"[28]"},{"why":"The standard tridiagonal eigensolver used as the independent numerical method for spectra and eigenfunctions.","marker":"[38]"},{"why":"Establishes non-negativity of $\\Theta_o$, justifying the $\\omega^2$ spectral parametrization.","marker":"[29]"},{"why":"Exhibits the nonseparability of polymer-quantized Hilbert spaces and the integral-space analogy used in the construction.","marker":"[15]"},{"why":"Supplies the effective-dynamics approximation for the bounce value $v_b(\\omega)$ used to set initial data.","marker":"[37]"}],"fun_headline_variants":["Integral Hilbert space matches single-sector dynamics","No added decoherence from sector averages in loop quantum cosmos","Quantum cosmology: integral sectors preserve semiclassicality","Sector-integrated Hilbert space passes semiclassicality test"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison treats semiclassicality as meaning small variance of the volume operator for Gaussian spectral profiles with a handful of parameters; if other observables, higher moments, or non-Gaussian states would reveal sector-induced spreading, the conclusion could fail.","fun_headline_variants_meta":{"raw":{"variants":["Integral Hilbert space matches single-sector dynamics","No added decoherence from sector averages in loop quantum cosmos","Quantum cosmology: integral sectors preserve semiclassicality","Sector-integrated Hilbert space passes semiclassicality test"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000364,"raw_usage":{"total_tokens":1977,"prompt_tokens":980,"completion_tokens":997,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":932}},"tokens_in":596,"tokens_out":997,"duration_ms":8704,"temperature":1.0,"reasoning_tokens":932,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:59:36.236901+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the volume trajectories with a basis accurate to better than $10^{-15}$ over the full range of $\\epsilon$ and $\\phi$, using a non-Gaussian spectral profile or an observable beyond the volume; if the sector-averaged variance then grows monotonically with $\\phi$ or the relative difference $\\delta\\sigma_r$ exceeds $10^{-9}$ away from recollapse points, the paper's claim of no added decoherence would be refuted.","supporting_citations":[{"cited_title":"Bojowald, Living Rev","cited_arxiv_id":null,"evidence_quote":"Supplies the model's quantization, the shots method, and the single-sector semiclassical dynamics that this work extends and cross-checks."},{"cited_title":"Loop quantization of the Gowdy model with local rotational symmetry","cited_arxiv_id":"1706.05673","evidence_quote":"Introduces the integral Hilbert space construction over superselection sectors, the approach tested here."},{"cited_title":"finite box","cited_arxiv_id":null,"evidence_quote":"Defines the mainstream LQC quantization, the area gap $\\Delta$, and the form of the evolution operator $\\Theta_\\Lambda$."},{"cited_title":"Homogeneous Loop Quantum Cosmology","cited_arxiv_id":"gr-qc/0303073","evidence_quote":"Provides the lapse choice $N=2V$ and regularization adjustments that modify the model relative to [24]."},{"cited_title":"Pawłowski and A","cited_arxiv_id":null,"evidence_quote":"The standard tridiagonal eigensolver used as the independent numerical method for spectra and eigenfunctions."}],"review_version":1}