{"id":"e71aba9e-6113-439e-ab7c-6be9c3f3d62b","arxiv_id":"2411.19601","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Quantum imaginary time evolution on small qubit registers reproduces Fermi-Dirac and Bose-Einstein thermal distributions for 1+1 dimensional field theories.","lead":"This proceedings paper tests whether the quantum imaginary time evolution algorithm can prepare thermal states of fermion and scalar field theories on small qubit lattices. The simulated momentum-space occupancies match Fermi-Dirac and Bose-Einstein curves, which suggests the approach works for tiny 1+1 dimensional systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The interacting-fermion analytical comparison is the load-bearing point: f_p0, f_p1, and m̃ are stated without derivation or parameter values, so the claimed agreement in Fig. 2 could reflect curve-fitting rather than independent evidence for QITE.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the interacting-fermion thermal-state factorization and the effective mass m̃ are under-specified and potentially fitted. My read of Sec. 3 confirms this. The free-fermion and scalar results are reasonable checks of the QITE algorithm, but the strongest demonstration of quantum simulation for interacting thermal field theories is Fig. 2, and that figure's analytical comparison is not self-contained. The manuscript points to reference [4] for the full derivation; if that reference supplies the factorization and the parameter values, the concern may be resolved. However, as a standalone proceedings paper, the central claim rests on an unstated derivation and unspecified simulation parameters. This is not an internal inconsistency, but it is a correctness risk in the sense that the reported agreement could be an artifact of the chosen analytical curves. The appropriate verdict remains CONDITIONAL: the claim is not contradicted, but it cannot be fully verified without the missing derivation, parameters, and exact-diagonalisation comparison. No change to the reader's verdict is needed, so I recommend UNCHANGED.","tokens_in":4629,"tokens_out":5373,"duration_ms":54542,"concrete_test":"Build the 5-qubit Hamiltonian from Sec. 3 (four ψ-qubits plus one homogeneous ψ_B qubit) with the bare parameters used in the simulation, and exact-diagonalize e^{−βH} on that Hilbert space. Compute E_Ω and E_0 directly from the Hamiltonian spectrum, form m̃ without fitting, and evaluate ⟨a_p† a_p⟩_β for each momentum mode. Compare these exact-diagonalisation results with both the analytic f_p0/f_p1 curves and the QITE simulation points. If the analytics match ED only after adjusting m̃, the Fig. 2 agreement is not an independent check; if they match with the unfitted spectrum, the conditional concern is resolved. Also report the QITE initial state and the imaginary-time step count used to realize β.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that QITE prepares thermal states whose momentum-space occupancies agree with Fermi-Dirac and Bose-Einstein distributions. For the free-fermion and scalar cases the comparison is plausible, but the interacting-fermion case in Sec. 3 is where the argument is least secure. The formulas for f_p0 and f_p1 depend on two sector partition functions Z_β^0 and Z_β^1 and on an effective mass m̃ = E_0 − E_Ω. The paper does not derive the factorization of the thermal state into these two sectors, does not define how E_0 and E_Ω are computed from the Hamiltonian, and does not state the bare parameters (m_0, M, g, β, a, N) used in the simulation. The summary claims comparison with exact diagonalisation, but no exact diagonalisation data or error bars are shown in Fig. 2. If m̃ is adjusted to match the simulation, then the agreement is not an independent test of QITE. Additionally, QITE as cited prepares a pure state e^{−βH/2}|ψ⟩; matching a Gibbs trace requires a specified initial state and a purification or averaging protocol, which is not described here. These omissions make the strongest numerical claim impossible to verify from the manuscript alone.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper is an ICHEP2024 proceedings contribution reporting quantum simulations of thermal states for 1+1-dimensional field theories on a small digital quantum simulator. The authors describe qubit encodings for Majorana fermions and a discretised scalar field, and they apply the quantum imaginary time evolution (QITE) algorithm to prepare approximate thermal states. They compare the momentum-space occupancies f_p with analytical Fermi-Dirac (free and interacting Majorana fermions) and Bose-Einstein (scalar) distributions in Figs. 1-3, and they claim strong agreement, including with exact diagonalisation methods.","tokens_in":4946,"tokens_out":3854,"duration_ms":31223,"significance":"If the central claim holds, the paper provides a useful small-scale demonstration that QITE can prepare thermal states of lattice field theories, with the free-fermion and scalar comparisons serving as independent benchmarks. The qubit mappings and the QITE implementation are concrete and reproducible in principle, and the free-field checks are not circular because the analytical distributions are externally derived. The main value is as a stepping stone toward studying thermalisation and thermal fixed points in real-time quantum simulation, as the authors state in the abstract and outlook. The significance is limited by the proceedings format: several technical steps, especially the interacting-fermion analytical expressions, are stated without derivation, and the claimed exact-diagonalisation comparison is not shown.","major_comments":[{"comment":"The analytical distributions for the interacting Majorana fermion model are stated without derivation: f_p^0 and f_p^1 are expressed in terms of sector partition functions Z_beta^0 and Z_beta^1 and an effective mass m̃ = E_0 - E_Omega, but the paper does not derive the factorization Z_beta = Z_beta^0 + Z_beta^1, does not define how E_0 and E_Omega are computed from the Hamiltonian, and does not state the bare parameters (m_0, M, g, beta, a, N) used in the simulation. If m̃ is adjusted to match the numerical data, the agreement in Fig. 2 would not be an independent test of QITE. The authors should provide the derivation or an explicit reference, define all quantities, and state the parameter values; they should also show the claimed exact diagonalisation data for this model, since Fig. 2 contains only the simulation and the analytical curves.","section":"Sec. 3, Eqs. defining f_p^0, f_p^1, Z_beta^0, Z_beta^1"},{"comment":"The manuscript defines the thermal expectation value as a Gibbs trace, but it does not specify how the trace is implemented with QITE. QITE as cited prepares the pure state e^{-beta H/2}|psi> (up to normalisation), so the reported f_p values require a specification of the initial state |psi> and the procedure for averaging over a complete basis or for using a purification. Without this protocol, the numerical results in Figs. 1-3 are not reproducible from the manuscript alone, and the comparison to a thermal ensemble is not fully justified.","section":"Secs. 1-2, QITE protocol"},{"comment":"The paper states in the summary that the results are compared with 'exact diagonalisation methods', but no exact diagonalisation data or error bars appear in any figure. Moreover, the simulations in Figs. 1-3 lack a complete list of input parameters (lattice spacing a, number of sites N, bare masses, coupling g or lambda, temperature beta, boson cutoff N_b, and number of qubits per site n_Q). For the scalar case, the approach to the Bose-Einstein distribution with increasing n_Q is presented without stating the truncation parameters or the finite-volume corrections, so it is unclear whether the agreement reflects convergence to the continuum distribution or merely to the digitised approximation. The authors should add a parameter table, show the exact diagonalisation results or remove the claim, and discuss truncation errors in the scalar-field digitisation.","section":"Figs. 1-3 and Sec. 4"}],"minor_comments":[{"comment":"The displayed free-fermion Hamiltonian has a summation structure that is hard to parse, with a sum over n appearing both outside and inside the second line; please rewrite it with clear parentheses and a single summation convention.","section":"Sec. 2, free-fermion Hamiltonian"},{"comment":"The caption states '4 qubits' and mentions two thermal limits, but it does not give the values of beta, the lattice spacing, or the number of momentum modes; please add these so that the reader can reproduce the curves.","section":"Fig. 1 caption"},{"comment":"The caption refers to 'both quasiparticles' but does not label which markers correspond to f_p^0 and f_p^1; please add a legend or explicitly identify the curves.","section":"Fig. 2 caption"},{"comment":"The text introduces N_b as a boson number cutoff but then uses N_phi = 2^{n_Q} as the local Hilbert space dimension and never specifies how N_b is chosen or how it relates to N_phi; please clarify this truncation procedure.","section":"Sec. 4, digitisation"},{"comment":"Reference [5] is listed as 'work in progress'; if the scalar-field results rely on this work, please provide a preprint identifier or report the results directly in the present paper.","section":"References"},{"comment":"The Wilson parameter r is introduced but its value in the simulations is never stated; please give the value used in Figs. 1 and 2.","section":"Sec. 2, Wilson term"}],"recommendation":"major_revision","confidential_remarks":"This is a short proceedings paper, so some lack of detail is expected, but the interacting-fermion analytical comparison is the central new element and it is currently unverifiable. The authors should either derive the factorization and m̃ from the Hamiltonian or point to a published derivation, and they must provide the simulation parameters and the exact diagonalisation data they claim in the summary. The 'work in progress' reference for the scalar section is also a concern for a published contribution. I believe the paper can be brought to an acceptable state with a major revision, but the current version does not allow the reader to assess whether the Fig. 2 agreement is a genuine QITE result or a fitted analytical curve."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a short proceedings, so I read it at that level. The free-field results are the strongest part. The QITE simulations for free fermions match Fermi-Dirac distributions, and the scalar field results approach Bose-Einstein as the digitization improves. The paper also does a good job explaining the qubit encoding, especially the Jordan-Wigner map for fermions and the field truncation plus discrete Fourier transform for scalars. The scalar digitization results in Fig. 3 are new, while the fermionic sections largely repeat the authors' JHEP paper [4].\n\nThe soft spot is Sec. 3. The analytical formulas for the interacting fermion model are stated without derivation: the two-sector factorization, the partition functions Z_β^0 and Z_β^1, and the effective mass m̃ = E_0 − E_Ω just appear. No parameter values (m_0, M, g, β, a, N) are given, and the claim of comparison with exact diagonalisation is unsupported because Fig. 2 shows only the analytical curves and simulated points, no exact diagonalisation data and no error bars. This makes the strongest numerical claim unverifiable from the manuscript alone. The stress-test concern that m̃ could be adjusted to match the simulation is legitimate, though I don't see direct evidence of curve-fitting—the definition is physically motivated but incomplete.\n\nI also note that QITE as described prepares a pure state e^{−βH/2}|ψ⟩, and matching a Gibbs trace requires a specified initial state and averaging or purification protocol. That step is not described, but for a proceedings this is a brevity issue rather than an error.\n\nThe paper is honest that the scalar section is work in progress, and the free-field checks are genuinely useful sanity checks for the QITE algorithm. I would not place much weight on the interacting-fermion agreement until the full derivation appears.\n\nThe intended audience is people working on quantum simulation of field theories, especially thermal-state preparation. They will get a useful snapshot of the program, and the free-field results are worth seeing. The paper deserves serious referee time because the central claim is concrete and testable, but the referee should ask for the derivation of m̃ or a direct citation to a self-contained version, and for simulated data with uncertainties.\n\nMy recommendation: engage with it, cite the JHEP paper for fermions, and treat Sec. 3 as provisional.","headline":"A credible QITE demonstration for free fields, but the interacting-fermion section is under-specified and needs either a derivation or a clear pointer to the full paper.","tokens_in":5432,"tokens_out":2467,"would_cite":false,"duration_ms":22847,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81P68","81T80"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that the quantum imaginary time evolution (QITE) algorithm can prepare the thermal (Gibbs) states of 1+1 dimensional fermionic and scalar field theories on a small digital quantum simulator, with measured momentum-space…","keywords":["quantum simulation","thermal field theory","QITE algorithm","Fermi-Dirac distribution","Bose-Einstein distribution","Majorana fermions","scalar field theory","qubit encoding"],"falsifier":"Run an exact diagonalisation of the interacting Majorana-plus-background Hamiltonian at the same parameters used for Fig. 2 and compute $\\langle a_p^{\\dagger}a_p\\rangle_\\beta$ and $\\langle a_p^{\\prime\\dagger} a_p^\\prime\\rangle_\\beta$ directly from the full thermal state; if the exact occupancies deviate from the paper's two-sector formulas with $\\tilde m=E_0-E_\\Omega$ determined from the spectrum, the central agreement claim is refuted.","tokens_in":4418,"feed_emoji":"⚛️","tokens_out":13992,"duration_ms":104423,"temperature":0.7,"pith_summary":"Four qubits are enough here for a free Majorana fermion field to produce momentum-space occupancies $f_p$ that coincide with the Fermi-Dirac distribution at high and low temperature. Adding a homogeneous Majorana background field introduces a second quasiparticle sector, and the paper reports that the simulated occupations $f_p^0$ and $f_p^1$ match a two-sector analytic formula with effective masses $\\tilde m$ and $\\tilde m+g$. For a scalar $\\phi^4$ theory with digitised field operators, the simulated occupations approach the Bose-Einstein distribution as the number of qubits per lattice site grows. The authors' stated motivation is that preparing such thermal states is the first step toward simulating thermalisation and thermal fixed points in real time, where classical sign-problem methods cannot go.","feed_headline":"Quantum imaginary-time evolution matches Fermi-Dirac curves","feed_subtitle":"Four-qubit quantum circuits reproduce the exact thermal occupancies of 1+1D fermion and scalar fields.","key_machinery":"The engine of the calculation is the QITE algorithm, which approximates the non-unitary operator $e^{-\\beta H}$ by a sequence of unitary gates chosen from local expectation values, so no auxiliary qubits are required. The fermion part of the argument rests on a staggered-lattice discretisation followed by the Jordan-Wigner transformation, which maps $N$ Majorana modes to Pauli operators on $N$ qubits. The scalar part rests on a digitisation scheme: each lattice site gets a finite Hilbert space spanned by field eigenvalues $\\varphi_\\alpha=\\Delta\\varphi(\\alpha-(N_\\varphi-1)/2)$, and the conjugate momentum operator is defined through a discrete Fourier transform, $\\Pi_n=\\bar m\\,\\mathcal{F}_n\\Phi_n\\mathcal{F}_n^{-1}$. These encodings matter because they determine exactly what state the simulator prepares and how the measured momentum-space occupancies are extracted from coordinate-space operators.","core_discovery":"The central claim is that quantum imaginary time evolution (QITE) prepares the Gibbs state $e^{-\\beta H}/Z_\\beta$ of a lattice-discretised quantum field theory, and that observables read from the prepared state reproduce the exact thermal field theory. In the free Majorana case the measured $f_p$ follow Fermi-Dirac curves for both $T\\gg m$ and $m\\gg T$. In the interacting case, the four-fermion interaction is implemented with a homogeneous spectator Majorana field $\\psi_B$, and the paper defines two partition functions $Z_\\beta^0$, $Z_\\beta^1$ and two occupation functions $f_p^0$, $f_p^1$; the simulation, run on four qubits for $\\psi$ and one for $\\psi_B$, agrees with those analytic forms. For the scalar theory, continuous fields are replaced by finite-dimensional operators on a local Hilbert space of dimension $N_\\varphi=2^{n_Q}$, with conjugate momentum built from a discrete Fourier transform; the simulated thermal occupancies converge toward Bose-Einstein as $n_Q$ increases. The paper takes this as evidence that qubit-based simulators can access equilibrium properties of interacting field theories, and as the groundwork for real-time thermalisation studies.","pith_inferences":["If the two-sector effective-mass ansatz is correct, the same factorization should survive at stronger coupling and on larger lattices, with $\\tilde m=E_0-E_\\Omega$ computable from the single-particle spectrum; those runs would either confirm or break the ansatz.","The QITE-prepared thermal state could be used as an initial condition for a subsequent real-time evolution on the same qubits, turning the equilibrium preparation directly into a thermalisation simulation.","The scalar digitisation results suggest a resource estimate the paper does not provide: the number of qubits per site and the circuit depth should scale polynomially with $N_\\varphi$ and $\\beta$, which could be checked numerically before hardware deployment."],"forward_implications":["If the QITE preparation works as demonstrated, finite-temperature observables in small lattice field theories can be obtained on digital quantum hardware without constructing the full Gibbs state classically.","Because the simulations run in coordinate space while the comparison is made in momentum space, the same preparation can be used for interacting theories whose interaction terms are non-local in momentum.","The interacting Majorana model predicts a doubled quasiparticle thermal spectrum, with occupation branches controlled by $\\tilde m$ and $\\tilde m+g$; this is a concrete signature of four-fermion interactions at finite temperature.","For scalar fields, thermal-state accuracy is controlled by the local register size $n_Q$: increasing it drives the prepared state toward the Bose-Einstein limit, giving a controlled digitisation error.","The successful preparation of equilibrium states is the paper's announced precondition for simulating real-time thermalisation and thermal fixed points in quantum field theory."],"supporting_citations":[{"why":"Supplies the QITE algorithm that implements imaginary time evolution as a sequence of unitary gates.","marker":"[3]"},{"why":"Provides the fermionic thermal field theory construction and the analytic two-sector comparison used in Section 3.","marker":"[4]"},{"why":"Provides the scalar-field simulation whose numerical results are shown in Fig. 3.","marker":"[5]"},{"why":"Gives the staggered-lattice discretisation and Jordan-Wigner qubit mapping for fermion fields.","marker":"[6]"},{"why":"Gives the scalar field digitisation and simulation approach that Section 4 extends.","marker":"[7]"}],"fun_headline_variants":["Qubits simulate thermal field theories accurately","Quantum imaginary-time evolution yields Gibbs states","Four-qubit circuits match thermal Fermi-Dirac curves","Simulating thermal quantum fields on a quantum computer","Quantum circuits reproduce thermal occupancies"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For the interacting fermion model, the analytical comparison assumes that the thermal state factorises into two independent sectors described by $Z_\\beta^0$ and $Z_\\beta^1$ with an effective mass $\\tilde m=E_0-E_\\Omega$; the paper does not derive this factorization or specify how $E_0$ and $E_\\Omega$ are computed, so if that assumption is wrong the agreement in Fig. 2 would be an artifact of the analytic curves.","fun_headline_variants_meta":{"raw":{"variants":["Qubits simulate thermal field theories accurately","Quantum imaginary-time evolution yields Gibbs states","Four-qubit circuits match thermal Fermi-Dirac curves","Simulating thermal quantum fields on a quantum computer","Quantum circuits reproduce thermal occupancies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000423,"raw_usage":{"total_tokens":2133,"prompt_tokens":869,"completion_tokens":1264,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":1197}},"tokens_in":485,"tokens_out":1264,"duration_ms":8782,"temperature":1.0,"reasoning_tokens":1197,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T06:01:15.622010+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run an exact diagonalisation of the interacting Majorana-plus-background Hamiltonian at the same parameters used for Fig. 2 and compute $\\langle a_p^{\\dagger}a_p\\rangle_\\beta$ and $\\langle a_p^{\\prime\\dagger} a_p^\\prime\\rangle_\\beta$ directly from the full thermal state; if the exact occupancies deviate from the paper's two-sector formulas with $\\tilde m=E_0-E_\\Omega$ determined from the spectrum, the central agreement claim is refuted.","supporting_citations":[{"cited_title":"Cuntín, W","cited_arxiv_id":null,"evidence_quote":"Provides the scalar-field simulation whose numerical results are shown in Fig. 3."}],"review_version":1}