{"id":"e88e6578-fdfb-410a-8201-acbab587fc4f","arxiv_id":"2501.11121","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"For the Kitaev-Hubbard chain, the finite-temperature decay rate of the Majorana edge mode follows an Arrhenius law with an effective gap systematically larger than the many-body gap.","lead":"In an interacting electron chain, the authors show that a Majorana edge mode survives for exponentially long times as the temperature is lowered, with a protection scale larger than the ordinary energy gap. This is a step toward understanding whether topological quantum information can remain stable at nonzero temperature.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Δ_eff fit may be corrupted by Lanczos truncation precisely in the large-β regime that defines the exponential slope.","rationale":"The reader's verdict is CONDITIONAL, and I agree that the EDOS ansatz is a modeling assumption worth testing. However, the more load-bearing weakness is the interplay between the Lanczos truncation point and the plateau/staggering decay in the large-β regime. The paper itself acknowledges (Appendix F1) that meaningful results require n_unstable to lie deep inside the plateau and after the staggered component has decayed. The figures suggest this condition is not met for the low-temperature points: for wβ = 2.35, b_n appears to keep growing up to n = 200, so no plateau is observed; and higher β shifts n* upward while pushing n_unstable downward. If the series is truncated while the SSR-type staggering is still present, the artificial chain attached to the homogeneous lead is effectively still topological, which artificially narrows the central EDOS peak and reduces γ(β). Because this bias is monotonic in β, it directly biases the fitted Δ_eff upward, undermining the headline inequality. The ED-based Γ_ϵ check (Sec. IVC) is a valuable independent signal that a low-energy ASZM exists for one parameter point, but it does not by itself certify the temperature-dependent extraction of Δ_eff for all cuts. A concrete numerical test—recomputing the highest-β γ values with larger bond dimension or with an independent time-evolution method—can settle whether the reported Δ_eff survives. Until then, the conditional verdict is appropriate; the paper should be asked to provide this convergence evidence before the quantitative claim is accepted.","tokens_in":25905,"tokens_out":5993,"duration_ms":62725,"concrete_test":"For each parameter point and β used in Fig. 6(a), especially the largest β, compute n* (the smallest n beyond which |h̃_n| stays below 5% of its maximum) and n_unstable for the largest feasible χ. If n_unstable < n* or b_n has not plateaued by the truncation point N, the homogeneous-lead continuation is invalid. Then, using a larger χ (or a more stable Lanczos variant) so that n_unstable > n*, re-extract γ(β); if Δ_eff changes by more than the statistical uncertainty, the reported Δ_eff is an artifact. Independently, for L = 16 compute the ACF with TDVP/MPO up to wt ~ 1/γ at wβ = 2.0 and 2.35 and compare the decay time with the Lanczos prediction; disagreement would confirm the truncation bias.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim Δ_eff > Δ rests on fitting γ(β) to exp(−Δ_eff β) in the large-β regime. In that regime, two requirements of the method are in tension. Section IVA and Fig. 2 show that increasing β shifts n* (where the staggered component h̃_n decays) to larger n, while Appendix F1 shows that the Lanczos instability sets in at smaller n_unstable for larger β. Moreover, for wβ = 2.0 and 2.35, Fig. 2 shows b_n still growing through n = 200, i.e., no plateau is reached within the computed range. The EDOS reconstruction (Appendix D) assumes a homogeneous semi-infinite lead with hopping b_N after truncation; if no plateau has been reached, this continuation is unjustified. If the series is truncated before the staggering has decayed or before b_n saturates, the artificial chain remains effectively in the topological regime at the truncation boundary, producing an artificially narrow Lorentzian and an underestimated γ(β). Since this bias grows with β, it steepens the log γ versus β slope, inflating Δ_eff. Thus Δ_eff > Δ may be partly a truncation artifact rather than a spectral property. The independent ED check (Sec. IVC) supports an ASZM below ~Δ_eff for one parameter point, but it does not validate the γ(β) thermometer used to define Δ_eff across the phase diagram.","agreement_with_reader":"partial"},"referee_report":null,"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper does something genuinely new—it runs the Lanczos algorithm on operators at finite temperature using MPOs, and uses it to extract the temperature-dependent decay rate of a Majorana edge mode in the Kitaev-Hubbard chain. The central observation is that the rate follows τ(β) ~ exp(Δ_eff β) with Δ_eff larger than the many-body gap. I think the qualitative claim is probably right, but the quantitative Δ_eff is softer than it looks.\n\nWhat's good: The method is a real step beyond the infinite-temperature Lanczos studies of Yates et al., and the implementation appears careful. They check convergence with bond dimension, compare against TDVP at infinite temperature, and include a distinct exact-diagonalization probe: the low-energy projected operator Γ_ε has a commutator that decays exponentially for ε below ~Δ_eff and a crossover near Δ_eff. That ED check is the strongest evidence in the paper, and it makes the qualitative story credible. They also ship data and code.\n\nThe soft spot is the extraction of γ(β). The EDOS is fitted to a Lorentzian-plus-background ansatz, and then γ(β) is fitted to an exponential. The paper does not give error bars on either fit, and there are only a handful of temperature points. More importantly, the large-β regime—exactly the one that defines the slope—is where the method is least controlled. Figure 2 shows that for wβ = 2.0 and 2.35 the Lanczos coefficients are still growing at n = 200 and have not reached a plateau, while Appendix F1 states that the Lanczos instability sets in at smaller n_unstable for larger β. Truncating before the plateau means the homogeneous-lead continuation is unjustified; if the truncated chain still looks topological at the boundary, the Lorentzian width is underestimated, and since the bias grows with β, it steepens the log γ vs β slope. That would inflate Δ_eff. The authors' own K-window test (Fig. 12) is only shown at wβ = 0.4, not in the problematic regime, so it does not address the concern.\n\nThe ED check partially rescues them: it supports the existence of an ASZM below ~Δ_eff at one parameter point. But it doesn't validate the γ(β) thermometer across the phase diagram, and Fig. 6(b) shows Δ_eff/Δ for many parameter points without uncertainties.\n\nOverall: worth publishing after revision, but a referee should push for error bars on the fits and a demonstration that the γ(β) extraction converges in the large-β regime, or at least an honest discussion of how the truncation bias affects the exponent. I'd take it to peer review.","headline":"Solid finite-temperature extension of the operator-Lanczos approach, with a credible but not yet airtight claim that Δ_eff > Δ; the quantitative exponent relies on fits in the large-β regime where the method is most fragile.","tokens_in":26775,"tokens_out":4099,"would_cite":true,"duration_ms":40577,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"At finite temperature, the Majorana edge mode of the interacting Kitaev-Hubbard chain decays with lifetime $\\exp(\\Delta_{\\rm eff}\\beta)$, where $\\Delta_{\\rm eff}$ exceeds the many-body gap.","keywords":["Almost strong zero modes","Majorana zero modes","finite temperature","Kitaev-Hubbard chain","Lanczos recursion","tensor network operators","edge density of states","topological protection"],"falsifier":"Compute the edge density of states at a fixed low temperature with a much higher bond dimension (or exact diagonalization for $L\\approx 16$ to 20) and fit no assumed lineshape: if the zero-frequency feature is not Lorentzian, or if its width does not follow $\\gamma=\\gamma_0 e^{-\\Delta_{\\rm eff}\\beta}$ over a range of $\\beta$ spanning at least a decade in $\\gamma$, the central claim fails. A direct check is to compare the autocorrelation function from the four-parameter fit with an independent long-time simulation: if the late-time decay deviates from the fitted exponential by more than the stated error, the effective gap is not a property of the model.","tokens_in":25709,"feed_emoji":"🧊","tokens_out":9292,"duration_ms":96622,"temperature":0.7,"pith_summary":"At zero temperature, Majorana zero modes at the ends of a gapped topological chain live forever; at infinite temperature, the same edge operators survive as long-lived \"almost strong zero modes.\" This paper fills in the gap between those extremes for the interacting Kitaev-Hubbard chain. It argues that at finite temperature the edge mode decays exponentially slowly, with lifetime $\\tau(\\beta)=\\exp(\\Delta_{\\rm eff}\\beta)$, where the effective activation scale $\\Delta_{\\rm eff}$ is systematically larger than the thermodynamic many-body gap $\\Delta$. The paper reads this as evidence that the degeneracy protecting the edge mode is not confined to the ground-state doublet but extends through a low-energy band of the spectrum. If true, this sharpens when cooling actually buys protection for topological qubits: the useful regime is set by $\\Delta_{\\rm eff}$, not by $\\Delta$.","feed_headline":"Cooling stretches Majorana edge-mode lifetimes exponentially","feed_subtitle":"In an interacting chain the edge-mode decay time is exp(Δ_eff β), with Δ_eff above the many-body gap.","key_machinery":"The machinery is the Lanczos recursion for Heisenberg time evolution, lifted to finite temperature and implemented with matrix product operators. The commutation superoperator $[H,\\cdot]$ turns the spreading of the seed operator $\\gamma_{1,a}$ into a single-particle hopping problem on a semi-infinite chain; the Lanczos coefficients $b_n(\\beta)$ are the hopping amplitudes, and the edge density of states $\\nu^E_\\beta(\\omega)$ of this artificial chain is the Fourier dual of the autocorrelation function. The paper approximates the EDOS by a narrow Lorentzian of width $\\gamma(\\beta)$ plus a gapped incoherent background (Eq. (11)), and justifies this by mapping the Lanczos chain to a dimerized single-particle chain in its topological regime attached to a homogeneous lead (Eq. (12)), so that the width extracted from a four-parameter fit is the inverse lifetime. A tensor network ansatz with bond dimension up to $\\chi=2000$ represents the operators and the thermal density matrix, making system sizes $L=22$ and arbitrarily long times accessible.","core_discovery":"Framed as the authors state it: the damping rate $\\gamma(\\beta)$ of the edge-Majorana autocorrelation function obeys $\\gamma(\\beta)=\\gamma_0 e^{-\\Delta_{\\rm eff}\\beta}$ in the temperature range studied, and the effective gap extracted from the slope is always above the many-body gap. At the representative point $\\mu/w=1.2$, $U/w=0.1$ they get $\\Delta_{\\rm eff}/w\\approx 2.7$ against $\\Delta/w\\approx 1.09$. Along three cuts through the topological phase the ratio $\\Delta_{\\rm eff}/\\Delta$ exceeds one, grows as $U\\to 0$ where an exact strong zero mode exists, and turns upward again for larger $U$. Exact diagonalization supports the picture: an operator $\\Gamma_\\epsilon$ built from opposite-parity eigenstates below an energy cut $\\epsilon$ has a commutator with $H$ that vanishes exponentially with system size for $\\epsilon$ below about $\\Delta_{\\rm eff}$, and retains a size-independent overlap with the edge operator $\\gamma_{1,a}$. The conclusion is that the low-energy sector below $\\Delta_{\\rm eff}$ hosts an approximate strong zero mode, so the topological protection extends over a finite energy window rather than only in the ground-state manifold.","pith_inferences":["Beyond the paper, the predicted $\\Delta_{\\rm eff}$ scale could be probed directly in edge-spin coherence experiments: the plateau height of the autocorrelation function and the subsequent decay time should show an activated dependence on temperature with slope $\\Delta_{\\rm eff}$.","The ratio $\\Delta_{\\rm eff}/\\Delta$ turns upward at large $U$, which the authors flag as open; one testable extension is whether this tracks incipient localization or a second topological regime in the phase diagram.","A natural next calculation is to identify which low-energy states form the degenerate pairs below $\\Delta_{\\rm eff}$, for instance by parity-resolved spectroscopy or entanglement diagnostics, since the paper shows their existence but not their microscopic nature.","The method's efficiency suggests extending the same operator-Lanczos tensor network pipeline to two-dimensional ladders or dissipative settings, where the finite-temperature scalar product would need modification."],"forward_implications":["At any finite temperature the Majorana edge operator no longer has infinite lifetime, but the decay is exponentially slow with activation energy $\\Delta_{\\rm eff}$ rather than the thermodynamic gap $\\Delta$.","The many-body spectrum below $\\Delta_{\\rm eff}$ is effectively doubly degenerate in opposite-parity pairs, so an edge operator projected to that window behaves like a strong zero mode.","Cooling a topological chain buys protection continuously: $\\tau(\\beta)=\\exp(\\Delta_{\\rm eff}\\beta)$ connects the zero-temperature infinite lifetime to the infinite-temperature almost strong zero mode.","As the interaction $U$ is tuned toward the integrable limit, $\\Delta_{\\rm eff}/\\Delta$ grows, recovering the exactly protected strong zero mode; the same analysis applies to parafermion chains, Floquet circuits, and particle-conserving Majorana ladders."],"supporting_citations":[{"why":"Supplies the definition of Majorana zero modes and the non-interacting Kitaev chain whose edge mode is the object studied.","marker":"[1]"},{"why":"Establishes the prethermal almost strong zero mode picture and the finite infinite-temperature lifetime that the paper extends to finite temperature.","marker":"[13]"},{"why":"Initiates the Lanczos and EDOS extraction of almost strong edge-mode lifetimes at infinite temperature, the method this paper generalizes.","marker":"[16]"},{"why":"Provides the continuum dimerized-chain-plus-lead approximation that motivates the four-parameter Lorentzian-plus-gap EDOS ansatz.","marker":"[17]"},{"why":"Constructs exact strong zero modes in integrable spin chains, giving the $U=0$ and $\\mu=0$ reference where $\\Delta_{\\rm eff}$ should diverge.","marker":"[21]"},{"why":"Supplies the recursion method and finite-temperature scalar products underlying the temperature-dependent Lanczos series.","marker":"[25]"},{"why":"Gives the universal operator growth hypothesis used to justify the plateau truncation and the attachment of a semi-infinite lead.","marker":"[26]"},{"why":"Shows how to represent the thermal density matrix as a tensor network, enabling finite temperatures within the matrix product operator approach.","marker":"[27]"},{"why":"States the standard exponential relaxation law with the many-body gap against which the paper's larger $\\Delta_{\\rm eff}$ is benchmarked.","marker":"[31]"},{"why":"Gives the spectral construction of many-body Majorana operators that the exact-diagonalization check $\\Gamma_\\epsilon$ is based on.","marker":"[51]"}],"fun_headline_variants":["Edge modes live exponentially longer as temperature drops","Effective gap Δ_eff sets exponential cooling law for edge modes","Warm chains still protect edge modes, but with Δ_eff gap","Finite-T edge-mode decay rate scales as exp(-Δ_eff β)","Better than thermodynamic gap: edge modes survive warm temps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the numerically computed edge density of states really is a narrow Lorentzian peak sitting on a gapped incoherent background, so that fitting its width $\\gamma(\\beta)$ to that shape yields the true decay rate; if the peak is non-Lorentzian or cannot be separated from the sidebands, the exponential scale $\\Delta_{\\rm eff}$ is an artifact of the fit.","fun_headline_variants_meta":{"raw":{"variants":["Edge modes live exponentially longer as temperature drops","Effective gap Δ_eff sets exponential cooling law for edge modes","Warm chains still protect edge modes, but with Δ_eff gap","Finite-T edge-mode decay rate scales as exp(-Δ_eff β)","Better than thermodynamic gap: edge modes survive warm temps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001028,"raw_usage":{"total_tokens":4357,"prompt_tokens":992,"completion_tokens":3365,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":608,"completion_tokens_details":{"reasoning_tokens":3290}},"tokens_in":608,"tokens_out":3365,"duration_ms":27205,"temperature":1.0,"reasoning_tokens":3290,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:37:42.857685+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the edge density of states at a fixed low temperature with a much higher bond dimension (or exact diagonalization for $L\\approx 16$ to 20) and fit no assumed lineshape: if the zero-frequency feature is not Lorentzian, or if its width does not follow $\\gamma=\\gamma_0 e^{-\\Delta_{\\rm eff}\\beta}$ over a range of $\\beta$ spanning at least a decade in $\\gamma$, the central claim fails. A direct check is to compare the autocorrelation function from the four-parameter fit with an independent long-time simulation: if the late-time decay deviates from the fitted exponential by more than the stated error, the effective gap is not a property of the model.","supporting_citations":[{"cited_title":"In the following we discuss how these two approximations influence the results","cited_arxiv_id":null,"evidence_quote":"Supplies the definition of Majorana zero modes and the non-interacting Kitaev chain whose edge mode is the object studied."},{"cited_title":"Jevtic and R","cited_arxiv_id":null,"evidence_quote":"Initiates the Lanczos and EDOS extraction of almost strong edge-mode lifetimes at infinite temperature, the method this paper generalizes."},{"cited_title":"Mahyaeh and E","cited_arxiv_id":null,"evidence_quote":"Provides the continuum dimerized-chain-plus-lead approximation that motivates the four-parameter Lorentzian-plus-gap EDOS ansatz."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Constructs exact strong zero modes in integrable spin chains, giving the $U=0$ and $\\mu=0$ reference where $\\Delta_{\\rm eff}$ should diverge."},{"cited_title":"Boundary Strong Zero Modes","cited_arxiv_id":"2305.16382","evidence_quote":"Supplies the recursion method and finite-temperature scalar products underlying the temperature-dependent Lanczos series."},{"cited_title":"Viswanath and G","cited_arxiv_id":null,"evidence_quote":"States the standard exponential relaxation law with the many-body gap against which the paper's larger $\\Delta_{\\rm eff}$ is benchmarked."},{"cited_title":"Michen, T","cited_arxiv_id":null,"evidence_quote":"Gives the spectral construction of many-body Majorana operators that the exact-diagonalization check $\\Gamma_\\epsilon$ is based on."}],"review_version":1}