{"id":"d17e5097-682a-4191-b489-26fa7fc7e248","arxiv_id":"2412.02733","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The phonon-traced ground state of a spin-Holstein cluster model is a bit-flip-corrupted cluster state whose von Neumann and Rényi-2 conditional mutual information detect SPT-to-trivial transitions at different couplings.","lead":"This paper studies a 2D cluster state (a quantum resource for measurement-based computing) coupled to local vibrational modes, and shows that tracing out the vibrations maps its ground state to a bit-flip-corrupted cluster state. It then shows two information-theoretic markers, the von Neumann and Rényi-2 conditional mutual information, each signal a symmetry-protected-topological to trivial phase transition, but at different coupling strengths.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The bit-flip channel Eq. (8) rests on the unverified product ansatz Eq. (6); residual phonon correlations in the exact ground state would renormalize or invalidate the channel and shift both predicted g/ω transition points.","rationale":"The reader's weakest_assumption correctly identifies Eq. (6) as the load-bearing premise, and my independent reading points to the same spot. The concern is not manufactured: the paper itself labels the product state an ansatz in the Fig. 1 caption, and Appendix C optimizes only within a coherent-state family, so α = 0 does not prove the absence of all residual phonon correlations. The algebra from Eq. (7) to Eq. (9) is internally consistent, and the Rényi-2 CMI numerics are consistent with the known bit-flip-channel results, which gives the paper independent support for the channel model itself. The absence of direct von Neumann CMI numerics is a secondary concern because the paper explicitly relies on an analytic mapping there. A direct numerical check of the ground-state factorization and of the resulting reduced state would settle the main risk. Since the reader already assigned CONDITIONAL on essentially these grounds, my stress-test does not move the verdict; it reinforces the condition. I would keep the verdict conditional, with the factorization check and the von Neumann CMI evidence as the conditions for full acceptance.","tokens_in":26490,"tokens_out":4193,"duration_ms":48621,"concrete_test":"Exact-diagonalize Eq. (1) on a small Lieb lattice (e.g., 2×2 or 3×3 vertex plaquettes with edge sites) with each phonon mode truncated to N_ph = 6–10 boson levels, for g/ω ∈ [0.2, 0.6]. Compute (i) the fidelity F = |(⟨Ψ̃_GS| ⊗ ⟨0|) U_p |Ψ_GS⟩|^2 and (ii) the trace distance D(ρ2_exact, ρ2_ansatz) between the exact traced spin state and the channel state Eq. (8) with p from Eq. (9). If F stays ≈1 and D approaches zero with increasing system size, the product ansatz is validated; if D is non-negligible (above the numeric noise floor) in the transition window g/ω ≈ 0.35–0.47, the predicted critical points are not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Eq. (6): |Ψ̃_tot⟩ = |Ψ̃_GS⟩ ⊗ |0⟩, asserted as the ground state of the polaron-rotated Hamiltonian. From this product form alone, the derivation in Appendix B gives the local bit-flip channel and p = 1/2[1 − e^{−2(g/ω)^2}]. But Eq. (3) shows that the polaron transform does not remove spin-phonon coupling: the operators Ĉ_i and Ŝ_i still contain bosonic creation and annihilation operators, so virtual phonon processes dress the spin couplings. Projecting onto the phonon vacuum yields the effective spin Hamiltonian Eq. (4), but the true ground state of Eq. (3) will contain admixtures of excited phonon states. Consequently, the exact reduced spin state after tracing out phonons is not exactly Eq. (8). Appendix C's variational calculation only minimizes within the restricted product family |Ψ̃_GS⟩ ⊗ |α⟩; the result α = 0 does not exclude non-product correlations such as squeezed phonon states or phonon-mediated many-body entanglement. Because both critical couplings (g/ω ≈ 0.352 and g/ω ≈ 0.4694) are obtained from the single error parameter p through Eq. (9), any correction to the channel changes both predicted transition points. This is a correctness risk for the paper's central claim that the microscopic spin-Holstein model realizes these mixed-state transitions, not merely a question of numerical accuracy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a two-dimensional cluster-state (symmetry-protected topological) Hamiltonian on a Lieb lattice locally coupled to bosonic modes. It first applies a polaron transformation and projects onto the phonon vacuum, obtaining an effective spin Hamiltonian whose ground state remains the cluster state at all couplings; the authors conclude that this pure-state approach does not detect the expected SPT-to-trivial transition. The paper then assumes that the ground state in the polaron-rotated frame factorizes as |Ψ̃_GS>⊗|0> (Eq. 6), traces out the phonons, and derives a local bit-flip channel with error rate p = (1/2)[1−e^{−2(g/ω)^2}] (Eqs. 7–9). Using the von Neumann conditional mutual information (CMI) and the Rényi-2 CMI, it predicts mixed-state transitions at g/ω ≈ 0.352 and g/ω ≈ 0.4694, respectively, and presents tensor-network data for the Rényi-2 CMI scaling.","tokens_in":26879,"tokens_out":7795,"duration_ms":86941,"significance":"If the central derivation holds, the paper offers a useful bridge between microscopic spin-phonon models and the recently developed mixed-state phase-transition framework: the phonon bath is shown to act as a local X-noise channel, and the two CMI diagnostics give different transition points, illustrating that von Neumann and Rényi-2 notions of mixed-state phases can differ. The algebraic derivation of the channel from the product ansatz is clean, and the tensor-network method avoids Monte Carlo sampling of noise configurations. However, the key factorization assumption Eq. (6) is not established, and the Rényi-2 scaling analysis fixes the critical point from prior work, so the quantitative predictions are less secure than the presentation suggests.","major_comments":[{"comment":"The bit-flip channel and both predicted critical couplings rely on the factorization |Ψ̃_tot> = |Ψ̃_GS>⊗|0> stated in Eq. (6) and treated as exact at the start of Sec. IV. The polaron-rotated Hamiltonian in Eq. (3) still contains the bosonic operators Ĉ_i and Ŝ_i inside the spin terms, so the exact ground state of that Hamiltonian will generally contain excited phonon components; the projection in Eq. (2) is a low-energy approximation. The variational calculation in Appendix C only minimizes within the product family |Ψ̃_GS>⊗|α>, so the result α=0 does not exclude squeezed or entangled phonon correlations. Because Eq. (9) maps p to g/ω monotonically, any correction to the channel shifts both transition points (g/ω≈0.352 and ≈0.4694). I ask the authors to either prove the factorization for this specific model or benchmark the exact reduced spin state against Eq. (8) using, for example, exact diagonalization or DMRG on small systems.","section":"Sec. IV (Eqs. 6–9)"},{"comment":"The von Neumann CMI transition is identified with p_vN≈0.11 by importing the critical point of the 2D ±J random-bond Ising model from Ref. [96]. The mapping from the decohered cluster state to the RBIM is plausible, but the paper does not derive the relation between the disorder probability in that RBIM and the bit-flip parameter p of Eq. (8), nor does it compute p_vN directly for the model at hand. Since this is one of the two central quantitative predictions, the authors should provide the derivation or determine p_vN numerically with the same tensor-network approach used for the Rényi-2 CMI and show consistency.","section":"Sec. IV A (after Eq. 23)"},{"comment":"The finite-size scaling analysis fixes p_c=0.178, taken from Ref. [35], and then reports that the data validate this value; this is partly circular. The scaling collapse should treat p_c as a free parameter and report its fitted value with an uncertainty, so that the only independent evidence for the critical point is the 'position of the peak approaches this value' statement in Fig. 4(a). In addition, the fitted exponent α changes from 2 (8×8 region, r′=2r) to 3.2 (4×4 region, r′=2r), so the scaling form of Eq. (24) is not universal; the authors should explain this geometry dependence or weaken the universality claim.","section":"Sec. IV B (Eq. 24 and Fig. 4)"},{"comment":"The numerical method is described as computing the Rényi-2 CMI exactly, but the boundary contraction keeps only 'significant' Schmidt values after an SVD truncation. No bond dimension, truncation error, or convergence data are reported. For the quantitative claims p_c=0.178 and the scaling exponents, the authors should report the numerical precision and demonstrate convergence with bond dimension.","section":"Appendix D (boundary contraction)"}],"minor_comments":[{"comment":"The notation is not fully consistent: 'Renyi-2' and 'Rényi-2' appear interchangeably, and 'von Newmann' appears in Sec. IV B where 'von Neumann' is meant.","section":"Throughout"},{"comment":"In the expression for κ_v, the fourth factor is written as ⟨0_e3|cos(B_e3)|0_e4⟩; the bra and ket labels should both be e3 (or e4, consistently).","section":"Appendix A (Eq. A7)"},{"comment":"Figure 2 is referenced parenthetically as 'see Fig. 2b' and 'see Fig. 2(a)' with inconsistent punctuation; please standardize the caption references and make the vertex/edge sublattice labels in the figure match the text.","section":"Sec. II and Fig. 2"},{"comment":"The probability P(l)=p^{|l|}(1−p)^{V_M−|l|} is stated without specifying that the string configuration l must be compatible with the region M; adding this clarification would improve readability.","section":"Sec. IV A (Eq. 18)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of cond-mat.str-el and addresses a topical question. The main concern is the unverified factorization ansatz in Eq. (6); if the authors can add a numerical comparison with exact ground states, the central claim would be substantially strengthened. I also note that the critical point p_c=0.178 is taken from the authors' own prior work and is then used as a fixed input in the scaling analysis; this is acceptable if transparent, but it should not be presented as independent validation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new thing here is the derivation of a local bit-flip channel from a microscopic spin-phonon Hamiltonian: p = 1/2(1 - e^{-2(g/ω)^2}) for the cluster state. The algebra in Appendix B is clean, and the tensor-network computation of the Rényi-2 CMI is exactly contracted and appears carefully done. The observation that the standard polaron pure-state approach cannot see the transition is simple but worth having on record.\n\nThe soft spot is the load-bearing factorization Eq. (6): the ground state in the polaron frame is asserted to be |Ψ̃_GS> ⊗ |0>. The polaron transform does not remove the spin-phonon coupling—Eq. (3) still contains bosonic operators inside C_i and S_i—so the exact ground state will have virtual phonon admixtures. Appendix C optimizes only over coherent displacements and finds α=0, which does not rule out squeezed or correlated phonon states. If those corrections are non-negligible, the channel becomes approximate and both predicted critical g/ω values shift. The paper flags this as an ansatz in places, but later states it as exact, which overstates the status.\n\nThe von Neumann CMI result is also weaker than it first appears: the transition at p_vN ≈ 0.11 is imported from the RBIM literature and mapped through Eq. (9), not computed for this model. The argument is plausible, but the paper does not provide direct evidence that the von Neumann Markov length diverges here. The Rényi-2 numerics are the real evidence, and they are consistent with the known Ising transition—though the scaling collapse fixes p_c = 0.178 from Ref. [35], which shares two authors, so the validation is partly circular. The fits themselves look fine, and the gauge-sector check landing at p_c = 0.5 is a good sanity check.\n\nMinor: no code or data are included, which would help for a numerical paper.\n\nOverall: the qualitative message—that mixed-state CMI diagnostics can see transitions the pure-state approach misses, and that the two CMI variants pick different points—survives my reading. The specific numbers should be treated as approximate until the factorization assumption is checked. This deserves serious refereeing. A good referee should ask for either a justification that phonon correlations beyond coherent displacements are controlled, or a clear statement that the channel is an approximate effective description.\n\nRecommendation: send to peer review with request for more direct von Neumann CMI evidence and an honest treatment of the approximation.","headline":"Clean derivation of a microscopic bit-flip channel and a careful Rényi-2 CMI study, but the key quantitative predictions rest on an unverified factorization ansatz and the von Neumann result is imported rather than computed.","tokens_in":27443,"tokens_out":2707,"would_cite":true,"duration_ms":29293,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.65.Ud","03.67.Mn","71.38.-k"],"model":"deepseek-v4-flash","headline":"In a spin-Holstein cluster model, tracing out phonons gives exactly a bit-flip channel, and two conditional-mutual-information diagnostics put the SPT-to-trivial transition at $g/\\omega \\approx 0.352$ and $\\approx 0.4694$, while the…","keywords":["mixed-state phases","spin-Holstein model","symmetry-protected topological order","conditional mutual information","polaron transformation","bit-flip channel","cluster state","decoherence"],"falsifier":"Compute the exact ground state of Eq. (1) on a finite Lieb lattice with tensor networks, trace out the phonons, and test whether the resulting spin density matrix equals the bit-flip channel Eq. (8); any deviation, such as correlated multi-qubit errors or unequal off-diagonal coherences, would falsify the channel derivation, and the von Neumann and Rényi-2 Markov lengths could then be computed from the actual reduced state to locate the true transitions.","tokens_in":26278,"feed_emoji":"⚛️","tokens_out":17623,"duration_ms":140795,"temperature":0.7,"pith_summary":"The paper studies a 2D cluster Hamiltonian locally coupled to phonon modes and asks whether the symmetry-protected topological order of the cluster state survives arbitrarily strong spin-phonon coupling. It shows that the standard pure-state route, which decouples spins and phonons with a polaron transformation and projects onto the phonon vacuum, always reproduces the same cluster ground state and therefore misses the expected transition. The authors then trace out the phonons entirely and find the spin state becomes the cluster state acted on by a product of local bit-flip channels, with error rate $p = \\frac{1}{2}\\left[1 - e^{-2(g/\\omega)^2}\\right]$ fixed by the microscopic coupling. Using two conditional-mutual-information diagnostics, they argue the mixed state undergoes SPT-to-trivial transitions at $g/\\omega \\approx 0.352$ (von Neumann Markov length) and $g/\\omega \\approx 0.4694$ (Rényi-2 Markov length), reflecting slightly different definitions of mixed-state phase.","feed_headline":"Phonon coupling turns cluster state into bit-flip errors","feed_subtitle":"Tracing out phonons reveals two topological-to-trivial transitions at g/ω ≈ 0.35 and ≈ 0.47.","key_machinery":"The load-bearing object is the polaron disentangling unitary $\\hat{U}_p = \\prod_i \\exp\\left[\\frac{g}{\\omega} X_i (\\hat{b}_i^\\dagger - \\hat{b}_i)\\right]$, together with the factorization assumption that the ground state in the rotated frame is $|\\tilde{\\Psi}_{GS}\\rangle \\otimes |0\\rangle$. Acting with this transformation and tracing out the phonons turns the pure cluster state $\\rho_1$ into $\\rho_2$ through a product of local bit-flip channels, with error rate $p = \\frac{1}{2}\\left[1 - e^{-2(g/\\omega)^2}\\right]$. The phase diagnostics are the von Neumann and Rényi-2 conditional mutual informations, whose exponential decay length (the Markov length) diverges at the transition; analytically these map to a random-bond Ising model and to a single-flavor Ising model, respectively.","core_discovery":"The central discovery is that the reduced spin state obtained by tracing out the phonons in the spin-Holstein cluster model is exactly the pure cluster state passed through a site-wise bit-flip channel: $\\rho_2 = \\prod_i \\left[(1-p)\\rho_1 + p X_i \\rho_1 X_i\\right]$ with $p = \\frac{1}{2}\\left[1 - e^{-2(g/\\omega)^2}\\right]$. This identification converts the microscopic spin-phonon coupling problem into a decoherence problem with a known statistical-mechanics mapping. The von Neumann conditional mutual information then maps to the free energy of a two-dimensional random-bond Ising model, whose critical point predicts a transition at $p_{vN} \\approx 0.11$ ($g/\\omega \\approx 0.352$), while a numerically-exact tensor-network computation of the Rényi-2 CMI locates the transition at the 2D Ising critical point $p_c \\approx 0.178$ ($g/\\omega \\approx 0.4694$). The gauge-sector contribution shows no meaningful transition, since its fitted critical point is $p = 0.5$, the infinite-temperature limit.","pith_inferences":["This channel derivation suggests a general recipe: any spin-phonon model whose coupling is mediated by the symmetry generator of an SPT state will trace out to a convex combination of the state and symmetry-flipped copies, with the channel structure determined by how the symmetry generator enters the coupling.","If the exact ground state contains phonon correlations beyond the single coherent displacement assumed in Eq. (6), the effective noise should acquire spatial correlations, turning the statistical mechanics mapping into a random-bond model with correlated disorder and shifting both critical couplings.","The gap between $p_{vN} \\approx 0.11$ and $p_c \\approx 0.178$ may be the same physics as strong-to-weak spontaneous symmetry breaking, a connection the paper leaves open; if that connection holds, the two Markov-length divergences should coincide with the SW-SSB boundary.","One could test the channel formula directly in a simulator by preparing a small cluster state, coupling it to engineered bosonic modes, and measuring the reduced spin density matrix to check whether it equals Eq. (8)."],"forward_implications":["The phonon bath acts as a local $X$-noise channel whose error rate $p = \\frac{1}{2}\\left[1 - e^{-2(g/\\omega)^2}\\right]$ is fixed by the microscopic coupling, so increasing $g/\\omega$ is exactly equivalent to increasing bit-flip decoherence.","Because the pure-state polaron analysis always returns the original cluster Hamiltonian with positive coefficients, it cannot detect the destruction of SPT order in this model, and mixed-state diagnostics are necessary.","The von Neumann and Rényi-2 Markov lengths identify distinct critical couplings ($g/\\omega \\approx 0.352$ and $\\approx 0.4694$), showing that 'mixed-state phase transition' is defined by the chosen equivalence relation.","The numerical Rényi-2 CMI can be computed exactly for systems up to $92 \\times 92$ with tensor networks, making the divergence of the Markov length a practical diagnostic for decohered SPT states.","Both critical couplings are of the same order as the phonon frequency $\\omega$, so the predicted SPT-to-trivial transition is expected in parameter regimes accessible to strongly coupled spin-phonon systems."],"supporting_citations":[{"why":"It supplies the mapping of bit-flip-decohered cluster states to the 2D Ising model and the critical rate $p_c \\approx 0.178$ used for the Rényi-2 prediction.","marker":"[35]"},{"why":"It provides the finite Markov length criterion for mixed-state phases and the tensor network method that the numerical section extends to exact contractions up to $92 \\times 92$.","marker":"[70]"},{"why":"It gives the critical value $p_{vN} \\approx 0.11$ of the random-bond Ising model used for the von Neumann CMI transition.","marker":"[96]"},{"why":"It defines the Rényi-2 conditional mutual information and its interpretation as a phase transition in the doubled Hilbert space.","marker":"[64]"},{"why":"It is the recovery map connecting a finite von Neumann Markov length to two-way channel connectivity, the equivalence relation behind the von Neumann diagnostic.","marker":"[72]"},{"why":"It is the standard polaron-plus-low-energy-projection approach whose failure to detect the transition the pure-state section demonstrates.","marker":"[94]"}],"fun_headline_variants":["Phonons induce bit-flip errors, revealing two transitions","Two mixed-state transitions emerge from spin-phonon coupling","Pure-state approach fails; mixed states reveal spin-phonon transitions","Spin-phonon coupling maps to bit-flip decoherence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation of the simple bit-flip channel assumes that after the polaron rotation the ground state factorizes exactly into the cluster spin state and the phonon vacuum; any residual phonon correlations beyond this single coherent displacement would invalidate the channel form and the predicted critical couplings.","fun_headline_variants_meta":{"raw":{"variants":["Phonons induce bit-flip errors, revealing two transitions","Two mixed-state transitions emerge from spin-phonon coupling","Pure-state approach fails; mixed states reveal spin-phonon transitions","Spin-phonon coupling maps to bit-flip decoherence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000739,"raw_usage":{"total_tokens":3339,"prompt_tokens":1022,"completion_tokens":2317,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":2245}},"tokens_in":638,"tokens_out":2317,"duration_ms":17914,"temperature":1.0,"reasoning_tokens":2245,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:13:55.250201+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact ground state of Eq. (1) on a finite Lieb lattice with tensor networks, trace out the phonons, and test whether the resulting spin density matrix equals the bit-flip channel Eq. (8); any deviation, such as correlated multi-qubit errors or unequal off-diagonal coherences, would falsify the channel derivation, and the von Neumann and Rényi-2 Markov lengths could then be computed from the actual reduced state to locate the true transitions.","supporting_citations":[{"cited_title":"Spin Liquid and Superconductivity emerging from Steady States and Measurements","cited_arxiv_id":"2408.07125","evidence_quote":"It provides the finite Markov length criterion for mixed-state phases and the tensor network method that the numerical section extends to exact contractions up to $92 \\times 92$."},{"cited_title":"Deterministic quantum teleportation of atomic qubits,","cited_arxiv_id":null,"evidence_quote":"It gives the critical value $p_{vN} \\approx 0.11$ of the random-bond Ising model used for the von Neumann CMI transition."},{"cited_title":"Topological order at nonzero tempera- ture,","cited_arxiv_id":null,"evidence_quote":"It defines the Rényi-2 conditional mutual information and its interpretation as a phase transition in the doubled Hilbert space."},{"cited_title":"Local uni- tary transformation, long-range quantum entanglement, wave function renormalization, and topological order,","cited_arxiv_id":null,"evidence_quote":"It is the recovery map connecting a finite von Neumann Markov length to two-way channel connectivity, the equivalence relation behind the von Neumann diagnostic."},{"cited_title":"Maximal entangle- ment versus entropy for mixed quantum states,","cited_arxiv_id":null,"evidence_quote":"It is the standard polaron-plus-low-energy-projection approach whose failure to detect the transition the pure-state section demonstrates."}],"review_version":1}