{"id":"a63cd5f1-9507-4426-9dbd-9a925a140f5b","arxiv_id":"2602.01794","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A semidefinite program decides whether any local-conservation-preserving Lindblad equation can match the leading-order steady state of a weakly coupled quantum system, yielding no-go results for few-qubit XXZ chains.","lead":"This paper turns the question “can a Lindblad equation describe this quantum steady state?” into a computer-checkable optimization problem. For several small quantum spin chains connected to heat baths, the method finds that no locally-conserving Lindblad description can reproduce the correct populations and coherences.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Numerical SDP results lack solver certificates and duality gaps, so the claimed 'rigorous no-go' verdicts (Eqs. 30, 34) are not established; a solver tolerance near δ_tol=1e-6 could change an 'impossible' conclusion.","rationale":"I read the paper as making two connected claims: (i) the SDP formulation exactly decides, in exact arithmetic, whether an LE satisfying local conservation laws can match the Redfield NESS to leading order; (ii) the numerical results provide rigorous no-go statements for specific models. Claim (i) is mathematically sound and is a genuine contribution. Claim (ii) is where the load-bearing weakness sits. The paper emphasizes duality and rigorous bounds, but the reported numbers come from a standard numerical solver without certificates. A 'rigorous no-go' requires a guaranteed lower bound on the optimal value, not merely a solver's approximate optimum. The figures show large margins, so I do not think the core qualitative conclusions are necessarily wrong; but the strength of the language ('fundamentally unattainable', 'rigorous no-go', 'all parameter regimes') exceeds what the presented evidence supports. The reader's weakest assumption focused on the Redfield equation being the correct leading-order description. That is a legitimate physical concern, especially because Refs. 34 and 35 are from the same group, and I partially agree with it. However, even granting Redfield accuracy, the SDP no-go results are not certified. Conversely, even with certificates, if Redfield were inaccurate, the target would be wrong. I see the numerical certification issue as more immediately load-bearing because it directly undermines the paper's advertised rigor, and it is concretely testable. The Redfield-accuracy question is somewhat mitigated by prior literature (including Ref. 35, which is independent), whereas no such mitigation exists for the missing certificates. Therefore I keep the reader's CONDITIONAL verdict: the SDP method is valuable, but the 'rigorous no-go' conclusions should be presented as numerical evidence supported by solver output, pending certified verification.","tokens_in":22129,"tokens_out":10562,"duration_ms":116946,"concrete_test":"Pick the strongest no-go point, e.g. N_L=N_R=1, g=0.01, β_L=1, β_R=5 (Fig. 2a). Re-run both SDPs (29) and (33) with a rigorous solver: SDPA-GMP in high precision or a verified SDP package, and extract explicit dual certificates / primal-dual gaps. Check whether the certified lower bound on τ_pop_opt and τ_pop_coh_opt exceeds δ_tol=1e-6. Repeat for at least 20 sampled parameter points across Figs. 2–5. If at any point the certified lower bound is ≤ δ_tol, the 'impossible' verdict is not established and the table entries 'Impossible' must be downgraded to 'not certified at that point'.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central decision rule is Eqs. (30) and (34): if τ_pop_opt or τ_pop_coh_opt ≥ δ_tol with δ_tol=10^-6, the desired LE is declared 'impossible', and the abstract/Table I call these 'rigorous no-go results'. The SDP formulation is correct in exact arithmetic, and the paper correctly notes that SDP duality can yield tight bounds. However, all numerical values in Figs. 2–5 and A1–A2 are obtained from CVX with default solvers, and the paper reports neither primal-dual gaps nor any certified lower bounds. A standard SDP solver returns an approximate optimal value within its tolerances; it does not by itself prove τ_opt ≥ δ_tol. If the true optimum at a displayed point is below δ_tol but the solver returns above it (or vice versa), the classification in Eqs. (30)/(34) is wrong. The margins in the figures are often orders of magnitude above δ_tol, so this is unlikely to change the qualitative picture for the plotted points; but the 'rigorous no-go / fundamentally unattainable' language and the extrapolation to 'all parameter regimes' require certified bounds or an analytical argument, neither of which is provided. This is the most load-bearing weakness in the paper's central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper asks whether, for a finite-dimensional non-degenerate system weakly coupled to thermal baths, there exists a completely positive, trace-preserving Lindblad master equation that (i) preserves local conservation laws in the sense of Eq. (22), (ii) reproduces the exact NESS populations at leading order, and/or (iii) reproduces both leading-order populations and coherences. The authors formulate these questions as semidefinite programs. The quantities τ_pop and τ_pop,coh measure the violation of the population and coherence conditions, and the decision rule compares optimized values with a tolerance δ_tol. A lemma lower-bounds the trace distance between the exact zeroth-order NESS and the NESS of any allowed LE. The method is applied to isotropic XXZ chains with one, two, or all boundary qubits coupled to bosonic baths with Ohmic spectral functions. Numerical results are reported: for N_L=N_R=1, no LE satisfying local conservation laws can reproduce even leading-order populations; for N_L=N_R=2, populations are possible only for sufficiently small inter-qubit coupling, while coherences are generally impossible; attaching all qubits in equilibrium admits a consistent LE, but non-equilibrium settings do not. The paper presents these as rigorous no-go results.","tokens_in":22424,"tokens_out":9957,"duration_ms":110518,"significance":"The SDP formulation is a genuinely useful contribution. It recasts a difficult existence question about Lindblad descriptions as a convex optimization problem, provides candidate LEs when feasible, extends the earlier thermalization-only SDP of Ref. 19 to NESS and coherences, and the code is publicly available. If the numerical no-go statements are supported by certified bounds or by appropriately softened language, the results will be valuable for the open quantum systems community. The lower-bound lemma is also an interesting idea, but its proof currently needs correction.","major_comments":[{"comment":"The central claims of 'rigorous no-go' and 'fundamentally unattainable' are not established by the numerical results as presented. The SDP formulation is exact in exact arithmetic, and duality can in principle provide certified bounds, but the paper reports no primal-dual gaps, no dual certificates, and no verified lower bounds. CVX with default solvers returns approximate optima within solver tolerances; it does not prove τ_opt ≥ δ_tol. Since the decision rule (30)/(34) hinges on the comparison with δ_tol=10^-6, the classification at points near the threshold could be wrong. The qualitative picture may survive because many plotted margins are orders of magnitude above δ_tol, but the 'rigorous' language requires, at minimum, certified lower bounds for the no-go points or an analytic argument. Please either provide such certificates or revise the no-go statements to 'no LE found by numeri","section":"§II E/F, Eqs. (30),(34); Figs. 2–5; Table I"},{"comment":"The 'exact' NESS is defined through the second-order Redfield dissipator L2, and the no-go results inherit this benchmark. The paper should state prominently that the claims are conditional on the Redfield equation giving the exact leading-order populations and coherences. This is standard in the weak-coupling regime, but the phrase 'exact NESS' is stronger than what is computed; a reader could mistake the target for the exact microscopic NESS. A concrete check would be to compare the SDP target with exact TCL or numerically converged results for at least one parameter point, or to state the required conditions explicitly (non-degeneracy, small ϵ, no level crossings, etc.).","section":"§II B/C, Eqs. (14)–(15), (19), (40)"},{"comment":"The paper claims impossibility 'across all parameter regimes' in several places, but the numerical study samples only finite ranges of β_L and g. Extrapolation from those plots to 'all parameter regimes' is not justified. Please restrict the conclusions to the ranges actually computed, or provide an analytic or certified argument covering the full parameter space.","section":"§III, Table I, Figs. 2–5"},{"comment":"The proof of Lemma 1 appears to contain a normalization error. Eq. (B12) writes the prefactors as 1/[2(N_M+N_R)] and 1/[2(N_M+N_L)], but Eq. (23) uses Hilbert-space dimensions d_M d_R and d_L d_M in the denominators. If N_M,N_R denote numbers of qubits rather than dimensions, the displayed formula is inconsistent with the model. Please correct this and re-derive the bound. The lemma is not needed for the SDP decision rule itself, but as a main-text mathematical claim it must be correct.","section":"Appendix B, Eq. (B12); Lemma 1"}],"minor_comments":[{"comment":"The norm in Eq. (32) is not defined; specify that it is the Hilbert–Schmidt norm on the operator space (or state the chosen p-norm explicitly).","section":"§II F, Eq. (32)"},{"comment":"The text says 'or get τ_pop,coh [Eq. (29)]'; this should refer to Eq. (33), not Eq. (29).","section":"§III C, text before Eq. (40)"},{"comment":"The derivation uses several inequalities without fully specifying the intermediate steps. Clarify how the bound Σ_k |C_{ka}| ≤ α is obtained, especially the role of d_L/d and the δ_{ka} term.","section":"Appendix B, Eq. (B16)"},{"comment":"Minor inconsistencies: Fig. 2 uses 'orange circle' while Fig. 3 uses 'red circle' for the same β_R=5.0 marker; also Table I says 'energy biases' in the caption, but the main text treats zero-bias except in Appendix C. Please harmonize captions.","section":"Table I and Fig. 3 captions"},{"comment":"The notation for ρ vs \\barρ, L2 vs L'_2, and bracket symbols is sometimes inconsistent (e.g., L2[ρ] vs L2(ρ)). A unified notation summary would improve readability.","section":"General notation"}],"recommendation":"major_revision","confidential_remarks":"The core SDP method is sound and publishable in principle, but the 'rigorous no-go' claims are currently ahead of the numerical evidence. I would not reject: the formulation is useful and the plotted margins suggest the qualitative results are likely correct. The revision must supply certified bounds or revise the language, and the Lemma 1 proof needs correction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does something useful: it turns the question 'does there exist a Lindblad equation that, at weak coupling, reproduces the exact NESS populations and coherences while preserving complete positivity and local conservation laws?' into a semidefinite program. That is a genuine extension of the authors' earlier thermalization-optimization work to non-equilibrium steady states, and the trace-distance lower bound in Lemma 1 is a nice addition. The SDP formulation itself is clean; Eqs. (29) and (33) are linear in the dissipator, and the decision rules in Eqs. (30) and (34) are sound in exact arithmetic.\n\nThe qualitative results are also credible and consistent with prior work: with only one boundary qubit per bath, even correct populations are impossible; with two boundary qubits, populations become feasible at weak coupling but coherences do not; and attaching baths to every qubit restores feasibility at equilibrium but not in the non-equilibrium case. These are the kind of no-go statements that people in the open-quantum-systems community will want to know about, even if they are not surprising in hindsight.\n\nNow the soft spots, in proportion. The main one is the 'rigorous no-go' language. The numerical results come from CVX with default solvers, and the paper reports neither primal-dual gaps nor certified lower bounds. A solver's approximate optimum does not by itself prove that the true optimum exceeds the tolerance delta_tol = 1e-6. In practice the plotted values are usually orders of magnitude above that tolerance, so I doubt the qualitative picture changes, but the word 'rigorous' is not earned for the numbers as presented. The fix is straightforward: report duality gaps, or provide an analytical certificate for at least one representative no-go regime.\n\nSecond, the abstract and Table I say 'impossible' and 'all parameter regimes' when the search is actually over LEs satisfying local conservation laws, with fixed bath parameters and no explicit time dependence. That is a well-defined class, and the authors correctly restrict to the most general local-conservation-preserving form, but it is still a restricted class. The abstract should not imply an unconditional statement about all Markovian embeddings.\n\nThird, the target 'exact' NESS is taken from the Redfield accuracy results cited as Refs. [34] and [35]. One is the authors' own earlier paper, but the other is an independent result by Fleming and Cummings, so this is not a crippling circularity. Still, the no-go conclusions are conditional on Redfield being the correct leading-order steady state, and that caveat should be explicit.\n\nWho is this for? Anyone constructing or benchmarking local Lindblad models for transport in small spin chains, and anyone interested in rigorous methods for open quantum systems. It deserves a serious referee. I would send it to peer review with the expectation of a revision that tightens the claims and adds numerical certificates or an explicit disclaimer about solver precision.","headline":"Solid SDP framework for testing whether local-conservation-preserving Lindblad equations can match Redfield-accurate NESS; numerical no-go results are plausible but the 'rigorous' label is a step too far without certified bounds.","tokens_in":22943,"tokens_out":2510,"would_cite":true,"duration_ms":29960,"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":"A semidefinite program can decide whether any Lindblad master equation can describe a weakly coupled quantum system—and certify when none can.","keywords":["semidefinite programming","Lindblad master equation","complete positivity","Redfield equation","non-equilibrium steady state","local conservation laws","XXZ qubit chain","Markovian approximation"],"falsifier":"Construct an explicit completely positive, trace-preserving Lindblad generator with Lindblad operators supported only on the boundary qubits that solves Eq. (25) [and Eq. (27)] in a parameter regime where the paper reports τ_pop^opt (or τ_pop,coh^opt) above δ_tol—for instance, for the N_L=N_R=1 XXZ chain at g=0.01 and β_L=β_R=1. Since the SDP optimum is a certified lower bound, any such construction directly contradicts the no-go claim. A cheaper check: run the same SDP on a system with a known exact Lindblad description and verify that τ_pop,coh^opt falls below tolerance.","tokens_in":21984,"feed_emoji":"🚫","tokens_out":5015,"duration_ms":48623,"temperature":0.7,"pith_summary":"This paper asks a yes/no question: given a small quantum system weakly coupled to heat baths, is there any Lindblad master equation that is physically consistent, obeys local conservation laws, and reproduces the exact steady-state populations (and possibly coherences) to leading order in the coupling? The authors show this question is a convex optimization problem—a semidefinite program—whose optimal value gives a rigorous answer: below a tolerance, a suitable Lindblad equation may exist and the solver returns one; above it, no such equation exists, so a Markovian description is fundamentally unattainable. Applied to isotropic XXZ qubit chains, the method yields no-go results: with one boundary qubit per bath, even correct populations are impossible; with two, populations become possible at weak inter-qubit coupling but coherences still fail, except near equilibrium. The work turns a vague worry about Markovian approximations into a computable, certified statement.","feed_headline":"A semidefinite program certifies when Lindblad equations fail","feed_subtitle":"For weakly coupled qubit chains, the method proves whether any completely positive master equation can match the exact steady state.","key_machinery":"The central object is the semidefinite program built from the most general form of a Markovian generator that is simultaneously CPTP and preserves local conservation laws: variable Hermitian Lamb-shift Hamiltonians on the left/right boundary and positive semidefinite rate matrices Γ^(L), Γ^(R). The conditions that the generator reproduce the Redfield zeroth-order populations (Eq. 25) or also the coherences (Eq. 27) are linear in these variables; hence minimizing the mismatch metrics τ_pop (Eq. 28) and τ_pop,coh (Eq. 32) is a convex SDP whose duality gives rigorous bounds. This machinery turns existence of a Lindblad description into a numerically certified decision.","core_discovery":"The central claim is that the existence of a completely positive, trace-preserving Lindblad equation that preserves local conservation laws and matches the exact leading-order non-equilibrium steady state can be decided by semidefinite programming. For a finite-dimensional system with a non-degenerate Hamiltonian, the authors define two mismatch metrics: τ_pop measures how far any such Lindblad equation's zeroth-order populations sit from the Redfield-accurate populations, and τ_pop,coh measures the distance between the candidate dissipator and the Redfield dissipator acting on ρ⁰_NESS. Minimizing these metrics over all Hermitian Lamb shifts and positive semidefinite rate matrices—the most g","pith_inferences":["The SDP criterion effectively gives a computable 'Markovianity certificate' for steady states; the same template could be extended to time-dependent dynamics or to other properties (correct currents, correlation functions) as long as they are linear in the generator variables.","The trace-distance lower bound suggests a quantitative grading: even when an LE is impossible, τ_pop^opt measures how badly the populations must be wrong, which could guide the choice between local and global master equations in practice.","Scanning the SDP over parameter space maps out the phase boundary of Markovian describability; near level crossings or large couplings the no-go results may weaken because the Redfield premise itself becomes questionable—an extension one could test with exact small-system simulations.","The method could be combined with exact microscopic simulations (beyond Redfield) to check whether the Redfield-defined 'exact' leading-order steady state is itself the right target in strongly coupled or near-degenerate regimes."],"forward_implications":["For isotropic XXZ chains with one qubit per bath, no completely positive, trace-preserving Lindblad equation preserving local conservation laws can give even correct leading-order populations, either in equilibrium or in a non-equilibrium steady state.","With two qubits per bath, an LE can give correct leading-order populations for small inter-qubit coupling (g ≲ 0.1), but correct coherences together with populations remain impossible except near equilibrium.","With all qubits attached to baths at the same temperature, a consistent LE is attainable; detaching a single site or making any bath temperature differ makes it unattainable.","Whenever the SDP certifies impossibility, the dynamics necessarily has non-Markovian features even at weak system-bath coupling.","The method is not restricted to XXZ chains: any finite-dimensional non-degenerate Hamiltonian with local baths can be plugged into the same SDP formulation."],"fun_headline_variants":["SDP test reveals when Lindblad equations are impossible","Semidefinite programming decides Lindblad feasibility","New method proves Markovian limits for quantum baths","When can Lindblad master equations be trusted? SDP says no","Certifying Lindblad failure via semidefinite programming"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole no-go machinery treats the second-order Redfield dissipator as the exact leading-order dynamics: the target populations and coherences the Lindblad equation is required to match are the ones Redfield predicts, so if Redfield is not actually correct to leading order in the regime studied, the SDP is certifying impossibility against a mis-specified target.","fun_headline_variants_meta":{"raw":{"variants":["SDP test reveals when Lindblad equations are impossible","Semidefinite programming decides Lindblad feasibility","New method proves Markovian limits for quantum baths","When can Lindblad master equations be trusted? SDP says no","Certifying Lindblad failure via semidefinite programming"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000349,"raw_usage":{"total_tokens":1784,"prompt_tokens":822,"completion_tokens":962,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":886}},"tokens_in":566,"tokens_out":962,"duration_ms":6917,"temperature":1.0,"reasoning_tokens":886,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T05:32:33.823860+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct an explicit completely positive, trace-preserving Lindblad generator with Lindblad operators supported only on the boundary qubits that solves Eq. (25) [and Eq. (27)] in a parameter regime where the paper reports τ_pop^opt (or τ_pop,coh^opt) above δ_tol—for instance, for the N_L=N_R=1 XXZ chain at g=0.01 and β_L=β_R=1. Since the SDP optimum is a certified lower bound, any such construction directly contradicts the no-go claim. A cheaper check: run the same SDP on a system with a known exact Lindblad description and verify that τ_pop,coh^opt falls below tolerance.","supporting_citations":[],"review_version":1}