{"id":"8cd0114f-d3cd-4d93-89ac-b4e3fe09012f","arxiv_id":"2507.06683","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Under Lindblad-type environmental coupling, the flux-free sector of the Kitaev honeycomb model relaxes to a maximally mixed state for most noise operators, its energy gap closes, and only the periodic Kitaev chain limit keeps a non-mixed steady state.","lead":"The paper simulates the Kitaev honeycomb spin liquid disturbed by an environment, modeled as noise that makes quantum states lose information. It finds that in most cases the final state of the system is featureless and the protected edge modes dissolve, a warning that substrate coupling must stay weak for Majorana detection experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The flux-free sector truncation is the load-bearing weak point: the physical spin-bath dissipator creates fluxes, so the steady-state and MZM conclusions rest on an unquantified approximation.","rationale":"The reader's weakest-assumption identification is exactly the concern that survives stress-testing: the LME is solved only within the flux-free sector, while the physical environment (STM substrate) couples through local spin operators that generate fluxes. This is not a peripheral technicality. The central claim—that the steady state is mostly maximally mixed, that the bulk gap closes, and that MZMs lose their spectral distinction—is stated as a property of the flux-free dissipative dynamics. If fluxes are dynamically generated, the effective fermionic Hamiltonian in Eq. (2) is no longer the correct generator, and the entire map in Tab. II, the QZE analysis, and the band-structure collapse could change qualitatively. The paper itself flags the uncertainty with 'may be a reasonable approximation,' but it does not supply the missing benchmark. The other issues noted by the reader (e.g., the topological-invariant statement in Sec. IV exceeding what was computed, and the fit-based fidelity extension) are real but less load-bearing: the topological transition claim is explicitly hedged, and the fidelity fit is presented as a description rather than a derivation. The truncation affects every quantitative result and the qualitative steady-state classification. Because the reader already issued a conditional verdict, no verdict change is needed; the concern reinforces CONDITIONAL without displacing it. The paper deserves credit for an explicit and systematic numerical map under an admittedly simplified model, but the central claim is conditional on a validation that has not been performed.","tokens_in":17755,"tokens_out":4533,"duration_ms":56163,"concrete_test":"Benchmark the LME on a small cluster (e.g., 2x3 honeycomb with the same zigzag cylinder geometry, 12 sites) in the full spin Hilbert space (dimension 2^12) with physical dephasing jump operators L_i = σ^z_i (or the local spin component appropriate for the substrate), and compare the exact steady state to the flux-free-sector result under the same parameters (Jx=Jy=Jz=1/3, γ=1 and γ=100). Compute ‖ρ_full^ss − ρ_fluxfree^ss‖_tr and the steady-state purity Tr[(ρ_full^ss)^2]; also monitor the flux-sector populations. If the trace distance is not small (say, >0.1) or purity is not near 2^{-12}, the flux-free truncation is invalid and the central conclusions need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's LME is solved entirely inside the flux-free sector of the cylinder, using quadratic fermionic jump operators (Tab. I). The physical STM substrate, however, couples via local spin operators, which in the Kitaev model create flux pairs on adjacent plaquettes and cause transitions between flux sectors. Section II acknowledges this: 'extending beyond the flux-free sector introduces quartic terms... restricting the analysis to the flux-free sector may be a reasonable approximation, and we, therefore, adopt this simplest approach.' The approximation is load-bearing because every central result—steady state ρmm, gap closing, MZM spectral indistinguishability, and QZE—is a property of the truncated dissipative dynamics. If flux-creating processes are non-negligible at the γ ~ Jα couplings where the QZE and gap closing are fastest, the steady state need not be maximally mixed, the band energies need not decay to zero, and the MZM-vs-bulk distinction could survive. The paper does not quantify the flux-creation rate, benchmark against the full Hilbert space, or provide a microscopic derivation of the bilinear jump operators. This is a correctness risk, not a disagreement with consensus; the truncation is a model assumption, and the paper's own caveat ('may be reasonable') indicates it is the least secure condition needed for the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Lindblad dynamics of the Kitaev honeycomb model on a zigzag cylinder, working entirely within the flux-free sector of the Majorana-fermion representation. For a catalog of quadratic fermionic jump operators, the authors compute the steady states of the Liouvillian, the time evolution of entropy, fidelity, band-resolved energies, and the spectral gap, and they identify a quantum Zeno effect with the fastest relaxation around γ ∼ Jα. The central numerical deliverable is Table II, which classifies which jump-operator sets drive the density matrix to the maximally mixed state, with the Kitaev-chain-with-periodic-boundary-conditions parameter regime as the main exception. The authors conclude that in most cases the energy gap closes in finite time and that the spectral distinction between Majorana zero modes and other bands disappears as the system relaxes.","tokens_in":17821,"tokens_out":3532,"duration_ms":44787,"significance":"If the results are taken at face value, the paper provides a useful reference catalog for how dephasing-type noise acts on the flux-free fermionic representation of the Kitaev honeycomb model, and it gives a concrete, falsifiable prediction (the γ ∼ Jα location of the Zeno minimum) that can be tested in controlled synthetic or analog quantum simulators. The numerical implementation is standard and internally consistent, and Table II systematically covers a broad set of jump operators. However, the significance for the stated STM motivation is limited by an explicit and unquantified model assumption: the dissipator is restricted to the flux-free sector even though the physical spin environment creates fluxes. The paper is therefore best read as a study of a fermionic model with quadratic Lindblad operators; its reach to the spin Kitaev model is conditional on a truncation that is acknowledged but not benchmarked.","major_comments":[{"comment":"The central approximation is the restriction of the Liouvillian to the flux-free sector, and the paper's own text in Sec. II states that this 'may be a reasonable approximation' without quantifying it. In the spin language, the local spin operators that a substrate couples to create flux pairs on neighboring plaquettes, so the dissipator used in the numerics drops exactly the physical processes most relevant for the STM setup invoked in the introduction. Because the steady-state classification, gap closing, and Majorana-mode indistinguishability are all properties of this truncated dissipative dynamics, this is a load-bearing assumption. The authors should either provide a quantitative estimate of the flux-creation rate (for example, by computing the relevant dissipator matrix elements to excited flux sectors), benchmark the truncated dynamics against the full spin Hilbert space on a small cluster, or substantially reframe the claims so that they are explicitly about the fermionic model rather than about a Kitaev spin layer on a surface.","section":"Sec. II (Eq. (3)) and Table I"},{"comment":"The claim that the fidelity result 'extends' the analytical expressions of Refs. [49,64] to the MZM case is not supported as stated. Equation (9) is the analytical expression of those references used as a fitting function with three free parameters (C, a, b) fitted to the authors' own numerics. The agreement is therefore a consistency check, not an independent derivation or a parameter-free prediction. The text should be reworded to make this distinction explicit, and the fitted parameter values and fit quality should be reported so that the reader can judge the transferability of the functional form.","section":"Sec. III B (Eq. (9))"},{"comment":"The statement in Sec. IV that 'a topological phase transition occurs over time' is stronger than what the computed quantity supports. The quantities shown in Figs. 6 and 7 are the energies E(k) of the eigenbands of the closed Hamiltonian weighted by the time-dependent density matrix; their convergence to zero shows gap closing in a spectral sense, but a mixed-state topological invariant is not evaluated. The authors already note this in Sec. IV, but the preceding discussion phrases the result as 'topological features of the band structure vanish.' This should be softened to 'the spectral distinction between the MZM band and the other bands vanishes within the flux-free truncated model,' which is exactly what the numerics show.","section":"Sec. III C and Sec. IV"}],"minor_comments":[{"comment":"There appears to be an unmatched square bracket in the Jy term of the Hamiltonian; the bracketing should be checked for consistency.","section":"Eq. (2)"},{"comment":"Equation (7) is called the Uhlmann fidelity, but the expression Tr[ρ1 ρ2] is the Hilbert-Schmidt overlap, not the Uhlmann fidelity. Please either use the standard Uhlmann definition or rename the quantity as a state overlap.","section":"Eq. (7)"},{"comment":"The heading of Sec. III E contains a typo, 'V ariation', which should be corrected.","section":"Section heading"},{"comment":"The notation in the conditions column of Table II is ambiguous; for example, 'L≥ 2, Jy = 0 ≠ 1' mixes a parameter condition with an inequality for deg(λ0). Please separate the parameter conditions from the degeneracy values.","section":"Table II"},{"comment":"The justification for the flux-free truncation relies on Refs. [47-49], but the text does not explain what those references actually establish. A sentence summarizing the physical or mathematical argument would make the limitation much clearer to the reader.","section":"Sec. II"}],"recommendation":"major_revision","confidential_remarks":"This is a competent numerical study with a clear catalog of results, but the central STM-motivated claim rests on a truncation that the authors themselves flag as an assumption. The revision should either supply a benchmark or visibly narrow the paper's scope to the flux-free fermionic model. If the authors can do that, the paper would be a solid contribution; without it, the abstract and conclusions overreach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you work on open Kitaev systems. The paper computes Lindblad evolution in the flux-free sector of the Kitaev honeycomb on a zigzag cylinder and maps out, for 14 jump-operator sets, which ones drive the steady state to the maximally mixed state and which do not (Tab. II). That table, plus the consistent finding of a quantum Zeno effect in tau, entropy, and fidelity, is genuinely new as far as I can tell, and the numerics look standard and internally consistent.\n\nThe authors are also honest about the main limitation: the whole calculation lives in the flux-free sector, and the jump operators are quadratic fermionic operators after Jordan-Wigner. A real substrate couples via local spin operators, which create flux pairs and move the system out of that sector. They cite Refs. [47-49] and say a restriction 'may be a reasonable approximation,' but they never quantify the flux-creation rate or benchmark against the full spin Hilbert space. That makes the central statements—steady state, gap closing, loss of MZM distinction—conditional on an unproven assumption. For the STM motivation in the abstract, that is a real gap, not a minor caveat.\n\nLesser issues: the fidelity 'extension' in Sec. III B is a fit of the analytical form of Refs. [49,64] to their own numerics (Eq. (9)), so it is a consistency check, not an independent prediction. The 'topological phase transition' language in Sec. IV is stronger than what they computed; they concede no invariant was evaluated. Eq. (2) is typeset in a garbled way that makes it hard to reconstruct the Hamiltonian, and there is no code or data release.\n\nWho this is for: people doing open Kitaev models or Lindblad studies of Majorana edge modes will find Tab. II and the QZE analysis a useful reference. The paper deserves a serious referee—the core numerical map is solid—but it needs revision: a small-system benchmark with flux sectors, a corrected Eq. (2), error bars or convergence checks, and code/data to back the figures. I'd send it out rather than desk reject.","headline":"Useful numerical map of the flux-free Kitaev sector under Lindblad noise, but the flux-free truncation keeps it a model study rather than a direct STM prediction.","tokens_in":18574,"tokens_out":2560,"would_cite":true,"duration_ms":79115,"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":"Environmental dephasing drives the Kitaev honeycomb's Majorana modes out of existence, closing the bulk gap in finite time.","keywords":["Kitaev honeycomb model","Majorana zero modes","Lindblad master equation","quantum Zeno effect","decoherence","flux-free sector","quantum spin liquid","open quantum systems"],"falsifier":"Exact-diagonalize the full spin Hamiltonian of a small honeycomb cluster together with the local $\\sigma^z$ dephasing Lindblad terms without any flux-free truncation, and compare the steady state and the gap-closing time with the flux-free prediction; a steady state different from the maximally mixed state, or a surviving gap at the predicted relaxation time, would falsify the central claim.","tokens_in":17332,"feed_emoji":"🧲","tokens_out":8119,"duration_ms":85147,"temperature":0.7,"pith_summary":"The authors ask what happens to the Majorana zero modes at the edges of the Kitaev honeycomb model when the spin liquid is coupled to a decohering environment, as in an STM measurement setup. Working inside the flux-free sector on a zigzag cylinder and using a Lindblad master equation with local dephasing-type jump operators, they find that for most jump operators the steady state is the maximally mixed state. As the system relaxes, all band energies decay to zero, the bulk energy gap closes, and the Majorana edge modes lose any spectral distinction from the bulk. The only clear exception is the parameter regime that reproduces the Kitaev chain with periodic boundary conditions, where the steady state is not maximally mixed. The authors also find a quantum Zeno effect in essentially every quantity they track: relaxation is fastest when the dissipation strength is comparable to the Kitaev couplings, and slower both for weaker and stronger dissipation.","feed_headline":"Substrate noise erases Kitaev honeycomb edge Majoranas","feed_subtitle":"Open-system calculation: the steady state is maximally mixed, the gap closes, and STM detection must race the relaxation.","key_machinery":"The analysis rests on the fermionized, flux-free-sector Hamiltonian of the Kitaev honeycomb model on a zigzag cylinder, in which the spins are mapped via a Jordan-Wigner transformation to complex fermions $\\alpha_{(l,k)}$ and the ground state lies in the flux-free sector. On top of this quadratic Hamiltonian the authors solve a Lindblad master equation, with jump operators obtained from the local $\\sigma^z$ operators (dephasing) together with related number, pair, and hopping operators. The time evolution is computed by vectorizing the Liouvillian via the Choi-Jamiolkowski isomorphism and exponentiating it, and the band structure is obtained from the Bogoliubov-de Gennes spectrum of the quadratic Hamiltonian. The quantum Zeno effect and the relaxation time are read off from the spectral gap of the Liouvillian, i.e. the slowest nonzero decay rate.","core_discovery":"Within the flux-free sector of the Kitaev honeycomb model on a zigzag cylinder, dephasing-type Lindblad dissipation drives the system to the maximally mixed state $\\rho_{\\mathrm{mm}} = \\mathbb{1}/2^L$ for almost all jump operators built from number operators, pair creation/annihilation, and incoherent hopping. Because this steady state has energy zero by particle-hole symmetry, every band converges to $E=0$, the gap closes, and the Majorana zero modes become indistinguishable from the other bands; in the authors' words, the topological features of the band structure vanish. The relaxation is governed by a quantum Zeno effect: the approach to the steady state is fastest when the dissipation strength $\\gamma$ is of order the Kitaev coupling $J_\\alpha$, and slows down for both smaller and larger $\\gamma$. The exceptions are the parameter sets realizing the Kitaev chain with periodic boundary conditions, where the steady state is not maximally mixed, and certain jump operators that preserve only half the Nambu degrees of freedom and produce slower decay.","pith_inferences":["The authors truncate to the flux-free sector and do not quantify the rate at which the full spin dissipator creates flux excitations; if that rate is comparable to $\\gamma$, the steady state and the gap-closing time could differ from their predictions.","A natural testable extension is to compute a time-dependent topological invariant; the paper shows the gap closes but does not determine whether a dynamical topological phase transition occurs.","The result suggests a general principle for local dephasing on Majorana-carrying systems: any environment that effectively measures the local fermion parity will destroy the spectral signature of zero modes, independent of the microscopic origin of the jump operators.","The fidelity fitting function, previously derived for the Kitaev honeycomb with periodic boundary conditions, is shown here to also describe open-boundary (zigzag) systems with Majorana modes for $\\gamma \\ll J_\\alpha$, extending its range of applicability."],"forward_implications":["If the central claim is correct, STM detection of Kitaev-spin-liquid Majorana zero modes requires the substrate coupling to be weak enough that the system does not relax to the maximally mixed state before the measurement.","The bulk gap closes in finite time under environmental coupling, so any spectroscopic feature that relies on the coexistence of a bulk gap and gapless edge states disappears on the relaxation timescale.","The quantum Zeno effect gives a concrete design rule: minimum robustness occurs when $\\gamma \\sim J_\\alpha$, so operating far from this crossover (either much weaker or much stronger dissipation) prolongs the topological signatures.","The Kitaev chain with periodic boundary conditions is the notable protected case, since its steady state is not maximally mixed and some distinction between bands survives."],"supporting_citations":[{"why":"Supplies the zigzag-cylinder flux-free fermionized Hamiltonian and the Majorana edge modes used throughout.","marker":"[33]"},{"why":"Defines the Kitaev honeycomb model and its exact fermionic solution that underlies the whole calculation.","marker":"[29]"},{"why":"Shows that zigzag edges of the Kitaev honeycomb model host Majorana zero modes, the object whose fate is studied.","marker":"[32]"},{"why":"Cited, together with [48], as the basis for treating the flux-free sector as a reasonable approximation under dissipation.","marker":"[47]"},{"why":"Cited, together with [47], as justification for restricting the dissipator to the flux-free sector.","marker":"[48]"},{"why":"Provides the analytic fidelity result that the paper uses as a fitting function and extends to open boundaries.","marker":"[49]"},{"why":"Gives the fidelity form for dissipative many-body systems that the paper adapts to fit entropy and fidelity data.","marker":"[64]"}],"fun_headline_variants":["Environment wipes out Kitaev honeycomb Majorana modes","Dephasing collapses Kitaev edge modes to mixed state","Kitaev honeycomb edge modes erased by decoherence","Maximally mixed state kills Majorana zero modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire calculation assumes the environment only acts inside the flux-free sector, so that the dissipator never creates flux excitations on the plaquettes; if real spin dephasing does create fluxes at a significant rate, the steady state and the fate of the Majorana modes could be different.","fun_headline_variants_meta":{"raw":{"variants":["Environment wipes out Kitaev honeycomb Majorana modes","Dephasing collapses Kitaev edge modes to mixed state","Kitaev honeycomb edge modes erased by decoherence","Maximally mixed state kills Majorana zero modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000159,"raw_usage":{"total_tokens":1242,"prompt_tokens":969,"completion_tokens":273,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":585,"completion_tokens_details":{"reasoning_tokens":207}},"tokens_in":585,"tokens_out":273,"duration_ms":3686,"temperature":1.0,"reasoning_tokens":207,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:59:13.426073+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Exact-diagonalize the full spin Hamiltonian of a small honeycomb cluster together with the local $\\sigma^z$ dephasing Lindblad terms without any flux-free truncation, and compare the steady state and the gap-closing time with the flux-free prediction; a steady state different from the maximally mixed state, or a surviving gap at the predicted relaxation time, would falsify the central claim.","supporting_citations":[{"cited_title":"Mizoguchi and T","cited_arxiv_id":null,"evidence_quote":"Supplies the zigzag-cylinder flux-free fermionized Hamiltonian and the Majorana edge modes used throughout."},{"cited_title":"Thakurathi, K","cited_arxiv_id":null,"evidence_quote":"Shows that zigzag edges of the Kitaev honeycomb model host Majorana zero modes, the object whose fate is studied."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Cited, together with [48], as the basis for treating the flux-free sector as a reasonable approximation under dissipation."},{"cited_title":"Kanega, T","cited_arxiv_id":null,"evidence_quote":"Cited, together with [47], as justification for restricting the dissipator to the flux-free sector."},{"cited_title":"Roberts, M","cited_arxiv_id":null,"evidence_quote":"Provides the analytic fidelity result that the paper uses as a fitting function and extends to open boundaries."},{"cited_title":"Tonielli, R","cited_arxiv_id":null,"evidence_quote":"Gives the fidelity form for dissipative many-body systems that the paper adapts to fit entropy and fidelity data."}],"review_version":1}