{"id":"9c479b9b-3937-480d-aa6c-93f7cc43a65e","arxiv_id":"2501.18611","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The equilibrium entropy of a small subsystem is ln d_A minus a 1/N^2 correction set by c^2 tr(X_A^2), where c and X_A are built from Hamiltonian moments and the initial energy spread.","lead":"This paper derives a formula for the entropy of a small piece of a quantum system that has reached equilibrium, assuming it thermalizes. The formula connects that entropy to the microscopic energy statistics of the Hamiltonian and the initial state, and it roughly matches exact diagonalization data for a spin chain.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The entropy formula stands or falls with Eq. (4), imported from Ref. [7]; if the ETH correction to the reduced density matrix is not exactly I_A/d_A + c X_A with O(1/N^2) trace-norm error, the leading finite-size entropy shift in Eq. (14) is not established.","rationale":"The paper's derivation from Eq. (4) onward is clean and the numerical trend in Fig. 1 is suggestive, but the central claim is conditional on an external ETH refinement. The weakest point is exactly Eq. (4), which is imported from Ref. [7] without re-derivation and is not directly verified in this paper. This is the same load-bearing concern identified by the Reader. The concrete test proposed above would settle the issue by checking the trace-norm distance to I_A/d_A + c X_A as a function of N, and would also test whether the final entropy prediction is independent of the initial state. Since the Reader already conditioned the verdict on this premise, no change of verdict is needed; the condition remains.","tokens_in":3494,"tokens_out":10393,"duration_ms":95732,"concrete_test":"Directly test Eq. (4) by exact diagonalization. For the N-site chain (9) at h_x = -1.05, h_z = 0.5, choose a zero-energy product initial state with exponential correlation decay (the same class used in the paper's comparison with Ref. [6]), compute the dephased reduced density matrix psi_inf_A for A = sites 1-3 (d_A = 8), and evaluate Delta_N = ||psi_inf_A - (I_A/8 + c X_A)||_1 for N = 8, 10, 12, 14, 16. Verify that Delta_N decreases as O(1/N^2), or at least faster than O(1/N); if Delta_N is O(1/N) or has an O(1) plateau, Eq. (4) is the failure point and the prediction in Eq. (14) is not supported. As a secondary check, repeat the fit in Eq. (13) with two different zero-energy initial states: the extracted alpha must match Eq. (14) independently of the initial state, since X_A and tilde_v are state-independent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (4) is the load-bearing input. It asserts that the infinite-time averaged reduced density matrix of an O(1) subsystem is I_A/d_A + c X_A up to O(1/N^2) in trace norm, with X_A fixed by H through Eq. (6). If the true ETH correction contains an additional operator not captured by X_A, or if the trace-norm error is O(1/N) rather than O(1/N^2), then the eigenvalues of psi_inf_A are not 1/d_A + c lambda_j + O(1/N^2). Since the entropy shift is quadratic in the eigenvalue deviations, an O(1/N) density-matrix error produces an O(1) contamination of the leading O(1/N^2) entropy correction via the cross term c*epsilon, and Eq. (8) -- and its specialization Eq. (14) -- fail. The later derivation is internally consistent once Eq. (4) is assumed: with epsilon_j = O(1/N^2), the cross term contributes only O(1/N^3), exactly as claimed. Thus the entire result inherits an unproven, imported premise. The comparison with Ref. [6] is only a single-parameter semi-quantitative check of the final entropy, not a direct verification of Eq. (4), so it cannot distinguish a missing operator in X_A from a coincidental agreement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the equilibrium von Neumann entropy of an O(1) subsystem in a local, translation-non-invariant Hamiltonian on N sites, assuming the eigenstate thermalization hypothesis (ETH). For an initial state with exponential correlation decay and zero total energy, it derives S(psi_inf_A) = ln dA - (c^2 dA tr(X_A^2))/2 + O(1/N^3), where c and X_A are expressed through Hamiltonian moments and the initial energy variance (Eqs. (5), (6), (8)). The derivation expands the entropy around the infinite-temperature state, relying on Eq. (4), which states that the time-averaged reduced density matrix is I_A/d_A + c X_A up to O(1/N^2) in trace norm. The paper applies this result to a spin-1/2 chain and compares the predicted quadratic coefficient with exact diagonalization data from Maceira and Läuchli, reporting semi-quantitative agreement.","tokens_in":3792,"tokens_out":4489,"duration_ms":41614,"significance":"If the central premise Eq. (4) is correct, the paper provides a parameter-free, analytical prediction for the leading finite-size correction to subsystem entropy in thermalizing systems, expressed directly in terms of the Hamiltonian and initial state. The derivation from Eq. (4) to Eq. (8) is transparent and internally consistent, and the explicit example calculation is valuable as a concrete test. The numerical comparison, though semi-quantitative, addresses a recent exact diagonalization study and the prediction in Eq. (14) is falsifiable. However, the main result is entirely conditional on an imported ETH expansion from Ref. [7], and the numerical support is not yet conclusive. The paper is concise and mostly clear, but its central claim requires either a self-contained derivation of Eq. (4) or a much stronger numerical verification.","major_comments":[{"comment":"The entire derivation of Eq. (8) rests on Eq. (4), which is imported from Eq. (15) of Ref. [7] without derivation or a precise statement of the assumptions under which the O(1/N^2) trace-norm bound holds. Since Ref. [7] is a companion paper by the same author, the present manuscript should make the central premise self-contained: either provide a derivation of Eq. (4) in an appendix, state it as an explicit assumption, or quote the precise theorem from Ref. [7] with its hypotheses. As written, the main result is a conditional statement whose key hypothesis is not established in this paper.","section":"Section 2, Eq. (4)"},{"comment":"The numerical validation is only a single-parameter, semi-quantitative test. The dots in Fig. 1 do not clearly converge to the dashed theoretical line by N=30, and the paper explicitly admits that 'one can not conclude whether the dots approach the dashed line as N → ∞'. This does not support the abstract's claim of providing a theoretical explanation for Ref. [6]'s numerical findings. A stronger test would involve multiple parameter sets, a finite-size scaling analysis with estimated O(1/N^3) corrections, or a direct check of Eq. (4) itself rather than only the final entropy coefficient.","section":"Section 3, Fig. 1 and Eq. (13)"},{"comment":"The fitted quantity α in Eq. (13) is extracted from finite-N exact diagonalization data, but the theoretical prediction in Eq. (14) is the coefficient in the thermodynamic limit. The paper does not discuss how the O(1/N^3) error term in Eq. (12) (or the finite-size corrections hidden in 'O(1/N^3)') affects the extracted α at N=20–30. Without such a discussion, the comparison between finite-N fits and the infinite-N prediction is uncontrolled.","section":"Section 3, Eq. (13) versus Eq. (12)"}],"minor_comments":[{"comment":"The expansion ln(1/d_A + μ_j) = ln(1/d_A) + d_A μ_j - (d_A^2 μ_j^2)/2 + ... is valid only if |d_A μ_j| < 1. For finite N this could fail for small systems; the paper should state explicitly that N is taken large enough so that this holds uniformly for all j.","section":"Section 2, Eq. (6) and text after Eq. (7)"},{"comment":"The expression for d_A tr(X_A^2) in Eq. (11) is lengthy and provided without derivation. Including a short derivation or a supplementary note would improve verifiability.","section":"Section 3, Eq. (10) and Eq. (11)"},{"comment":"The Hamiltonian in Eq. (9) uses periodic boundary conditions implicitly (σ_i^z σ_{i+1}^z for i=N should wrap around), but this is not stated. Please clarify the boundary condition.","section":"Section 3, Eq. (9)"},{"comment":"The figure would be clearer with error bars or a discussion of the uncertainty in the fitted α values, and with a statement of the initial state used in the exact diagonalization data of Ref. [6].","section":"Section 3, Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The paper's central result depends almost entirely on Eq. (4) imported from the author's own previous work (Ref. [7]). While self-citation is not itself a problem, the lack of a derivation or precise statement of conditions makes the paper more of an application of Ref. [7] than a standalone derivation. The numerical evidence is also weak for the claimed 'theoretical explanation'. I would encourage the editor to ask for either a self-contained derivation or a more rigorous numerical test before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nShort take: this is a legitimate small step forward. The new output, Eq. (8) for the leading O(1/N^2) finite-size correction to the entropy of an O(1) subsystem, and the explicit rational prediction alpha ≈ 0.818 for the Maceira–Läuchli chain, are genuinely new and parameter-free. Given the input assumption, the algebra from Eq. (4) to Eq. (8) is straightforward and internally consistent. I checked the eigenvalue expansion: because the coefficient c is O(1/N), a trace-norm error of O(1/N^2) in the reduced state contributes only at O(1/N^3) in the entropy, so the formula is not hiding a cross-term contamination.\n\nThe paper also deserves credit for honesty. The numerical comparison is semi-quantitative, and the author explicitly says one cannot yet conclude whether the data approach the dashed line. That is the right attitude for a paper making a finite-size prediction.\n\nWhere the soft spots are: the load-bearing premise is Eq. (4), imported without derivation from the author's own Ref. [7]. The whole result stands or falls on whether the ETH correction to the reduced density matrix really is I_A/d_A + c X_A with an O(1/N^2) trace-norm error. If there is an additional operator, or if the bound is only O(1/N), then the entropy shift is contaminated at leading order and Eq. (14) fails. Self-citation is not itself a flaw, but here the dependency is real and the paper should either re-derive Eq. (4) or state the precise theorem and assumptions from Ref. [7]. The numerical check is a one-parameter fit to the final entropy, not a direct test of Eq. (4), so it cannot distinguish a missing operator from a coincidence. Also, the initial state used in the numerics is not described in this paper, and the fit has no error bars.\n\nThese are addressable issues, not fatal ones. The paper is a short, checkable corollary of an existing formal program, and it makes a concrete falsifiable prediction. I would send it to peer review. The referee should ask for a self-contained derivation of Eq. (4) or a precise citation, and a fuller numerical comparison with specified initial state and uncertainties. This is exactly the kind of paper that reading-group members can take apart in an afternoon, which is a good sign.\n\nMy bottom line: useful, honest, worth engaging with, but read Ref. [7] first.","headline":"A clean, honestly presented corollary: the entropy formula is new and parameter-free, but it inherits its central premise from the author's companion paper, so Referee and reader both need to check Eq. (4) before trusting Eq. (14).","tokens_in":4306,"tokens_out":2651,"would_cite":true,"duration_ms":26392,"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":"For thermalizing quantum systems, the equilibrium von Neumann entropy of a small subsystem is $\\ln d_A$ minus a computable correction fixed by the Hamiltonian and the initial state.","keywords":["eigenstate thermalization hypothesis","subsystem entropy","von Neumann entropy","finite-size corrections","local Hamiltonians","exact diagonalization","thermalization","reduced density matrix"],"falsifier":"For the same spin chain at a different value of the transverse fields, compute the time-averaged reduced density matrix by exact diagonalization for $N=20$ to $34$, fit the entropy to $\\ln 8 - \\alpha((v-\\tilde v)/(2N))^2$, and check whether $\\alpha$ tends to $d_A\\operatorname{tr}(X_A^2)/(2\\tilde v^4)$. If the fitted prefactor moves away from the predicted value as $N$ increases, the $O(1/N^2)$ ansatz is wrong.","tokens_in":3285,"feed_emoji":"⚛️","tokens_out":7687,"duration_ms":63535,"temperature":0.7,"pith_summary":"Assuming the eigenstate thermalization hypothesis, the paper derives a closed-form expression for the equilibrium von Neumann entropy of a subsystem whose size stays fixed as the system grows. The formula expresses the entropy as $\\ln d_A$ minus a positive correction built from traces of powers of the local Hamiltonian and the energy variance of the initial state, up to $O(1/N^3)$ errors. Because every quantity in the correction is computable from the Hamiltonian and the initial state, the result turns a microscopic finite-size effect into a parameter-free prediction. The paper also shows that this prediction tracks, at least semi-quantitatively, the trend seen in recent exact-diagonalization data for a spin chain.","feed_headline":"Small-subsystem entropy gets a parameter-free formula","feed_subtitle":"Under ETH, the finite-size correction is fixed by Hamiltonian moments and the initial state's energy variance.","key_machinery":"The central object is the ETH ansatz for the equilibrated reduced density matrix, Eq. (4): for a subsystem of $O(1)$ size, $\\psi_\\infty^A = I_A/d_A + c X_A + O(1/N^2)$ in trace norm. The scalar $c=(\\tilde v-v)/(2N\\tilde v^2)$ measures how far the initial state's energy variance lies from the infinite-temperature value, and the operator $X_A$ is assembled from traces of $H$, $H^2$, and $H^3$ over the complement. Expanding $-(1/d_A+\\mu)\\ln(1/d_A+\\mu)$ and using $\\sum_j \\mu_j=0$ turns this ansatz into the entropy formula.","core_discovery":"The paper establishes that, under the eigenstate thermalization hypothesis, the von Neumann entropy of a fixed-size subsystem after equilibration is $\\ln d_A - \\frac{c^2 d_A}{2}\\operatorname{tr}(X_A^2) + O(1/N^3)$, where $c$ encodes the difference between the initial state's energy variance and the infinite-temperature variance, and $X_A$ is an explicit operator built from $\\operatorname{tr}(H^3)$, $\\operatorname{tr}_{\\bar A}(H)$, $\\operatorname{tr}(H^2)$, and $\\operatorname{tr}_{\\bar A}(H^2)$. This is the leading finite-size correction to the maximally mixed entropy, and it contains no fitting parameters.","pith_inferences":["A natural extension, not pursued in the paper, is that the same expansion should hold for R\\'enyi entropies with an $\\alpha$-dependent prefactor replacing $1/2$.","The formula could serve as a diagnostic for numerical thermalization simulations: the fitted prefactor should settle at the predicted value as $N$ grows.","A direct test of the underlying ansatz would be to diagonalize the time-averaged reduced density matrix at finite $N$ and check that its leading correction lies in the operator subspace spanned by $X_A$."],"forward_implications":["The equilibrium entropy of a small subsystem depends on the initial state's energy variance through $c$, not just on the Hamiltonian alone.","For any local traceless Hamiltonian, the leading correction to $\\ln d_A$ is negative, so equilibration slightly suppresses subsystem entropy below its maximal value.","The formula provides an analytic explanation of the finite-size scaling seen in exact-diagonalization studies of thermalizing spin chains.","The prediction applies to non-translation-invariant local Hamiltonians and to any initial state with exponential correlation decay and zero total energy."],"supporting_citations":[{"why":"Supplies the equilibrium reduced-density-matrix ansatz, Eq. (15), from which Eq. (4) is imported.","marker":"[7]"},{"why":"Provides the exact-diagonalization data for the quadratic factor that the prediction is compared against.","marker":"[6]"},{"why":"Foundational statement of eigenstate thermalization, which the derivation assumes.","marker":"[1]"},{"why":"Foundational formulation of eigenstate thermalization used to frame the argument.","marker":"[2]"}],"fun_headline_variants":["Subsystem entropy gets a parameter-free ETH formula","No fitting parameters: entropy of small subsystems under ETH","ETH yields exact finite-size entropy correction for subsystems","Parameter-free entropy formula for small subsystems under ETH","Small-subsystem entropy: no free parameters under ETH"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation depends on the assumed form of the equilibrated reduced density matrix of a small subsystem, namely the maximally mixed state plus a specific operator correction with an error that scales as $1/N^2$; if the true correction contains additional operators at that order, the entropy formula fails.","fun_headline_variants_meta":{"raw":{"variants":["Subsystem entropy gets a parameter-free ETH formula","No fitting parameters: entropy of small subsystems under ETH","ETH yields exact finite-size entropy correction for subsystems","Parameter-free entropy formula for small subsystems under ETH","Small-subsystem entropy: no free parameters under ETH"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000266,"raw_usage":{"total_tokens":1514,"prompt_tokens":750,"completion_tokens":764,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":366,"completion_tokens_details":{"reasoning_tokens":690}},"tokens_in":366,"tokens_out":764,"duration_ms":6305,"temperature":1.0,"reasoning_tokens":690,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:14:32.231488+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the same spin chain at a different value of the transverse fields, compute the time-averaged reduced density matrix by exact diagonalization for $N=20$ to $34$, fit the entropy to $\\ln 8 - \\alpha((v-\\tilde v)/(2N))^2$, and check whether $\\alpha$ tends to $d_A\\operatorname{tr}(X_A^2)/(2\\tilde v^4)$. If the fitted prefactor moves away from the predicted value as $N$ increases, the $O(1/N^2)$ ansatz is wrong.","supporting_citations":[{"cited_title":"High-precision simulation of finite-size thermalizing systems at long times","cited_arxiv_id":"2406.05399","evidence_quote":"Supplies the equilibrium reduced-density-matrix ansatz, Eq. (15), from which Eq. (4) is imported."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundational statement of eigenstate thermalization, which the derivation assumes."},{"cited_title":"Srednicki","cited_arxiv_id":null,"evidence_quote":"Foundational formulation of eigenstate thermalization used to frame the argument."}],"review_version":1}