{"id":"60a9ecc5-09dd-4654-bceb-e295f26e2fcb","arxiv_id":"1908.09245","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A large variational calculation places the dissipation-driven localization transition of two qubits in a common Ohmic bath at alpha_c = 0.316(8), and predicts a first-order transition for strong antiferromagnetic coupling.","lead":"Using a large variational wave function, the authors compute the ground state of two quantum bits coupled to a shared low-frequency environment, locating the localization transition at a coupling about 0.316, far higher than earlier simulations suggested. A generalist might care because this controversy sits at the heart of designing open quantum systems and suggests the transition type changes when the two qubits repel each other.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The §3.1 filter requiring |⟨σz⟩| monotone in α is applied to the same order parameter whose discontinuity defines αc; αc=0.316(8) is therefore not established independently of that filter.","rationale":"The reader correctly identifies the variational convergence and the Section 3.1 exclusion rule as the weakest points; my stress-test sharpens this into a concrete circularity. The central claim is that αc≈0.316, far above previous NRG/QMC values, and the principal evidence is an abrupt jump in ⟨σz⟩ combined with bath observables that are analyzed after the same monotonicity filter has been applied. Because the filter is imposed on the order parameter whose discontinuity is the reported transition signature, the numerical result cannot be considered independent of the method unless the filter is removed or justified by an energy-based selection that reproduces the jump. The agreement with the simple analytic boundary αc=K/4 for strong antiferromagnetic coupling is a real supporting argument, but it applies only to large K and does not validate the K=0 value. The Bethe-ansatz slope mismatch (3.83 versus 2) is acknowledged in the text and further undercuts the assumption of quantitative convergence. I thus maintain the reader's conditional verdict: the paper is worth publishing only if the central number survives a test that removes the selection rule or is confirmed by an independent method. This does not require rejection, since the calculation is legitimate, reproducible in principle, and the claims are falsifiable.","tokens_in":17779,"tokens_out":5017,"duration_ms":55537,"concrete_test":"Re-run the NVM calculation for Δ=0.025, K=ε=0, α∈[0.28,0.35] with no cross-α selection: for each α, keep only the lowest-energy state over all restarts and annealing runs, and plot |⟨σz⟩|, Δr/Δ, and ⟨B↑↑|B↓↓⟩ versus α without imposing monotonicity or jump uniqueness. If the sharp jump still occurs at α≈0.316 and Δr collapses at the same point, the filter is not load-bearing. If the unfiltered data show a smooth crossover or the jump moves, the central αc claim depends on the selection rule. As an auxiliary check, compute the unfiltered Δ/⟨σx⟩ slope and compare with the Bethe-ansatz value of 2; a persistent value near 3.83 would indicate that the variational state is not converged in the relevant regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the data-selection rule stated in Section 3.1: 'some data from the metastates have to be discarded according to the criterion that the absolute value of ⟨σz⟩ monotonically increases with α, and the sharp jump in ⟨σz⟩ is unique.' The transition point is then read off from the discontinuity in ⟨σz⟩ (Fig. 1b), with the bath observables in Figs. 4 and 7 aligned to that same filtered magnetization. This is close to circular: imposing monotonic increase plus a unique jump forces the retained curve to have exactly the qualitative shape used to locate αc. The problem is sharpened by Fig. 3, which shows competing states at α=0.292 differing in energy by only about 3×10⁻⁷ while differing hugely in ⟨σz⟩; energy alone cannot reliably choose the 'true' state at that resolution. The Bethe-ansatz check in Fig. 2a inset also weakens rather than strengthens the convergence claim: the slope of Δ/⟨σx⟩ is 3.83(3), not the expected 2, so the ansatz is not quantitatively converged in the weak-tunneling Ohmic regime where αc is inferred. Unless the selection rule is shown to be harmless, the central value αc=0.316(8) remains vulnerable to the possibility that the filter manufactured the jump that is then interpreted as the transition.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a variational study of the two-impurity spin-boson model with a common Ohmic bath, using a multi-D1 coherent-state ansatz with up to 500 bath modes and 6 coherent-state superpositions (over 12,000 variational parameters). The central claim is a ground-state localization transition at alpha_c = 0.316(8) for tunneling Delta = 0.025 and zero spin-spin coupling, which the authors assign to the Kosterlitz-Thouless universality class, and a first-order transition for strong antiferromagnetic coupling. The paper also computes bath observables such as variances, correlation functions, renormalized tunneling, and coherent-state displacements, extracting critical exponents and a phase diagram in the (alpha, K) plane. The authors argue that previous estimates (NRG 0.18, QMC 0.22, VM2 0.125) are underestimated due to bias, lack of continuum convergence, or trapping in metastable states.","tokens_in":18028,"tokens_out":3804,"duration_ms":42131,"significance":"If the central result holds, the paper would resolve a long-standing controversy by placing the transition at significantly stronger dissipation than most earlier estimates, and it would provide a consistent picture of both spin and bath critical properties across the transition. The work is valuable for introducing a high-dimensional variational approach with a linear discretization grid, which yields lower ground-state energies than logarithmic-grid calculations, and for systematically studying bath correlation functions near a Kosterlitz-Thouless transition. The claimed alpha_c = 0.316(8) and the proposed first-order behavior for large antiferromagnetic coupling are falsifiable predictions that could be checked by independent methods. However, the central value rests on a post-hoc data-exclusion rule and on convergence benchmarks that are not fully quantitative; these issues must be resolved before the main claims can be considered established. The paper does not provide code or raw data, so reproducibility is limited.","major_comments":[{"comment":"The first-order transition claimed for strong antiferromagnetic coupling (K = 3.0) is inferred in Fig. 10(a) from two straight-line fits to E_g with slopes of 0.00 and 2.00, with no quoted uncertainties or statistical test for the discontinuity of dE_g/dalpha. The accompanying analytical estimate alpha_c = K/4 uses the classical displacement f_k = lambda_k/omega_k and neglects renormalization effects and the residual tunneling Delta_r, so it is only a heuristic. To support the first-order claim, the authors should provide a quantitative measure of the discontinuity (e.g., the jump in dE_g/dalpha with error bars), demonstrate that it is robust to the variational convergence, and confirm with an independent observable such as a bimodal order-parameter distribution or a symmetry-breaking indicator in the variational wave function.","section":"Section 3.1, Fig. 1(b)"}],"minor_comments":[{"comment":"The derivation of the second derivative of the ground-state energy and the exponential fit (inset of Fig. 1(a)) would benefit from explicit error bars on the fitted exponent b = 23.2(2) and a justification of the exponential form, given that the paper then uses the absence of a discontinuity to infer the Kosterlitz-Thouless class.","section":"Abstract and conclusions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a well-known controversy and proposes a substantially different critical coupling, so the topic is appropriate for the journal. The main risk is the circularity of the data-selection rule: the same observable that defines the transition is used to filter the data, and the energy differences between retained and discarded states are smaller than the variational accuracy claimed in Fig. 3. If the authors can demonstrate that the filter is harmless by locating the transition with independent observables and providing a systematic convergence study, the paper could become a significant contribution. I would advise the editor to request those additions before accepting."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper you should know about reports a large-scale variational (multi-D1 ansatz, linear discretization) study of the two-impurity spin-boson model in a common Ohmic bath. The headline claim is that the localization transition sits at alpha_c ~ 0.316(8) for Delta = 0.025, well above the NRG 0.18 and QMC 0.22 estimates, and that strong antiferromagnetic coupling (K = 3) makes the transition first order with alpha_c = K/4.\n\nWhat is genuinely new here is the numerical phase diagram and the consistent set of bath observables. The trial wave function is a legitimate extension of the authors' earlier work, with more than 12,000 variational parameters and a linear grid, and it does produce lower variational energies than the log-grid ED and Silbey-Harris-style comparisons. The power-law scaling of bath correlation functions in the delocalized phase and the collapse of the data onto a KT-like picture are reasonable as far as they go. The alpha_c = K/4 line for large antiferromagnetic K is also a nice consistency check, even if it follows from their own energy-balance argument.\n\nThe soft spot the stress-test note puts its finger on is real. Section 3.1 states that data from metastable states are discarded when |<sigma_z>| is not monotonically increasing in alpha and the jump is not unique. That same filtered <sigma_z> is then used to read off alpha_c. At alpha = 0.292 the competing states differ in energy by 3e-7 but in magnetization by about 0.85, so energy alone cannot choose. Without a demonstration that the filtering is harmless, alpha_c = 0.316(8) is not independent of the selection rule. The Bethe-ansatz slope discrepancy (3.83 vs 2) further suggests the ansatz is not quantitatively converged in the weak-tunneling regime. The ED benchmark is small (M = 12, Ntr = 4, log grid) and gives a somewhat different alpha_c ~ 0.26. No code or data is shipped, so the error bars overstate confidence.\n\nI would not call this a rejection-level paper. The method is serious, the controversy is real, and the authors are honest about the metastability problem; they just do not yet resolve it. A serious referee should be engaged, but the letter should ask for an independent method (e.g., VMPS or unbiased QMC) or at least transparent reporting of raw and filtered data plus systematic error estimates. This is a paper for the spin-boson and open-quantum-systems community, and the true critical coupling will only be settled by a method that does not rely on the same order parameter it filters.","headline":"A serious but not-yet-conclusive variational claim for the two-impurity Ohmic spin-boson critical coupling; the central number depends on a data-selection rule that needs an independent check.","tokens_in":18648,"tokens_out":2734,"would_cite":false,"duration_ms":27494,"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":"This paper claims that two spins sharing one Ohmic bosonic bath localize only at dissipation strength $\\alpha_c=0.316(8)$ for tunneling $\\Delta=0.025$, well above earlier estimates, and that the transition is Kosterlitz-Thouless for…","keywords":["two-impurity spin-boson model","Ohmic bath","Kosterlitz-Thouless transition","ground-state phase transition","multi-D1 ansatz","variational method","quantum dissipation","coherent-state expansion"],"falsifier":"An unbiased, continuum-extrapolated calculation of the same model at $\\Delta = 0.025$, $K = 0$, $\\varepsilon = 0$ that finds the magnetization jump at $\\alpha_c$ below $0.25$, or a variational run with more coherent states and modes that moves the jump by more than the quoted uncertainty, would falsify the central value. A simpler partial test is the Bethe-ansatz slope: the paper fits $\\Delta/\\langle\\sigma_x\\rangle = (2\\alpha - 1)\\omega_c$ with slope $3.83$ rather than $2$, so a calculation reproducing slope $2$ while keeping the jump near $0.316$ would settle the issue.","tokens_in":17500,"feed_emoji":"⚛️","tokens_out":7271,"duration_ms":68000,"temperature":0.7,"pith_summary":"The paper argues that the ground state of two spins coupled to a common Ohmic bath undergoes a sharp localization transition at $\\alpha_c = 0.316(8)$ for tunneling $\\Delta = 0.025$, considerably larger than the values $0.18$, $0.22$, and $0.125$ reported by earlier NRG, QMC, and variational studies. The discrepancy is attributed to small bias fields that smooth the magnetization jump and to variational minimization becoming trapped in metastable states. Using a high-dimensional coherent-state trial wave function and a linear discretization of the bath spectrum, the authors find that the transition is of Kosterlitz-Thouless type when the spin-spin coupling is ferromagnetic or absent, but becomes first order for sufficiently strong antiferromagnetic coupling. If correct, the result revises the phase diagram of a basic two-qubit dissipation model and identifies bias and metastability as the causes of the earlier scatter.","feed_headline":"Two qubits in one bath localize only at dissipation 0.316","feed_subtitle":"Critical dissipation rises to 0.316, and strong spin coupling turns the transition first order.","key_machinery":"The load-bearing object is the multi-D1 variational ansatz: a trial ground state built as a superposition of up to six coherent bosonic states for each of the four spin configurations ($\\uparrow\\uparrow$, $\\uparrow\\downarrow$, $\\downarrow\\uparrow$, $\\downarrow\\downarrow$), with independent displacement parameters per bath mode, giving more than 12,000 variational parameters for $M = 500$ bath modes. Minimizing the energy with simulated annealing and a strict convergence criterion is meant to escape the metastable states that the paper argues have contaminated earlier estimates. A second, equally important choice is the linear discretization of the Ohmic spectral density ($\\Lambda_k = k/M$) instead of the logarithmic grids with $\\Lambda > 1$ used in earlier work; the paper claims the linear grid yields lower energies and the correct continuum limit, and it is this combination of ansatz and grid that produces the sharp magnetization jump at $\\alpha_c = 0.316(8)$.","core_discovery":"On the paper's own terms, the central discovery is that the two-impurity spin-boson model with a common Ohmic bath localizes at $\\alpha_c = 0.316(8)$ for $\\Delta = 0.025$, $K = 0$, $\\varepsilon = 0$, with a discontinuous jump of $\\langle\\sigma_z\\rangle$ from $0$ to $-1$ marking the transition. All derivatives of the ground-state energy remain continuous, while bath observables such as the basis correlation $\\langle B_{\\uparrow\\uparrow}|B_{\\downarrow\\downarrow}\\rangle$ and the renormalized tunneling $\\Delta_r$ change sharply, supporting a Kosterlitz-Thouless transition. The same variational machinery yields a phase diagram in which ferromagnetic spin-spin coupling barely moves the boundary, while strong antiferromagnetic coupling pushes it to $\\alpha_c = K/4$ and switches the transition to first order, as signaled by a kink in the ground-state energy. The paper also reports power-law critical behavior of bath observables, such as $2\\Delta X_b - 1 \\sim 1/\\omega_k$ with exponent $1.00(1)$ and $\\mathrm{CorX} \\sim 1/\\sqrt{\\omega_k}$ with exponent $0.500(3)$, connecting the Ohmic criticality to the two-dimensional XY universality class.","pith_inferences":["As an extension the authors do not pursue: re-running unbiased QMC or NRG with continuum extrapolation at $\\Delta = 0.025$ should move $\\alpha_c$ toward $0.32$ if the paper's attribution of all lower values to bias and metastable trapping is right; if it does not, the variational selection rule is the culprit.","The slope mismatch in the Bethe-ansatz check ($3.83$ versus the expected $2$ for $\\Delta/\\langle\\sigma_x\\rangle$) is reported but not resolved; I infer that the ansatz is not yet quantitatively exact near criticality, so an independent calculation is needed before $\\alpha_c = 0.316(8)$ hardens.","The optimal-displacement form $|f_k| = |p_k| = \\lambda_k/(\\omega_k + \\chi)$ with $\\chi \\propto \\Delta_r$ suggests a concrete testable extension: measure $\\chi$ in the single-impurity Ohmic spin-boson model, where the exact Kosterlitz-Thouless boundary is known from the anisotropic Kondo mapping, and check whether the same variational machinery recovers $\\alpha_c \\approx 1$."],"forward_implications":["Earlier NRG and QMC boundaries ($0.18$ and $0.22$) would be underestimates caused by small bias fields smoothing the magnetization jump and by metastable states; the true boundary in the unbiased continuum limit would be near $\\alpha_c = 0.316(8)$.","For two qubits sharing an Ohmic bath, delocalized coherence would survive up to $\\alpha \\approx 0.31$ at weak tunneling, roughly 40% stronger dissipation than the previously quoted QMC value.","At $K \\le 0$ the transition is Kosterlitz-Thouless: no energy-derivative discontinuity, but abrupt order-parameter jumps in $\\langle\\sigma_z\\rangle$ and bath correlations, analogous to the universal jump in the two-dimensional XY model.","For strong antiferromagnetic coupling the boundary follows $\\alpha_c = K/4$ and the transition is first order; the delocalized phase is then the antiparallel spin state with energy $-K/4$ and noninteracting bath modes, while the localized phase is the parallel state with energy $-2\\alpha + K/4$."],"supporting_citations":[{"why":"NRG calculation that gives $\\alpha_c \\approx 0.18$ with a small bias field; it is the main baseline the paper argues is underestimated.","marker":"[36]"},{"why":"QMC simulation giving $\\alpha_c \\approx 0.22$; the comparison baseline the paper claims is unreliable due to convergence and metastability.","marker":"[39]"},{"why":"Variational calculation giving $\\alpha_c = 0.125$ under a bias of $10^{-5}\\omega_c$; the paper attributes the low value to the imposed bias.","marker":"[41]"},{"why":"Earlier Silbey-Harris variational result of $\\alpha_c = 0.5$; it supplies the VM1 comparison and the constrained ansatz the multi-D1 form goes beyond.","marker":"[38]"},{"why":"Bethe-ansatz and Kosterlitz-Thouless expectations for the Ohmic spin-boson model, including the $\\Delta/\\langle\\sigma_x\\rangle = (2\\alpha - 1)\\omega_c$ relation used as a check.","marker":"[16]"},{"why":"Variational treatment of the single-impurity Ohmic spin-boson model with coherent-state expansions; it supplies the method and observables the paper extends to two impurities.","marker":"[27]"},{"why":"Earlier application of the multi-D1 ansatz to spin-boson problems; it is the methodological source of the trial wave function.","marker":"[29]"},{"why":"Study of bath observables near the sub-Ohmic spin-boson transition; it provides the framework for using bath correlations and variances as transition indicators.","marker":"[46]"}],"fun_headline_variants":["Two qubits localize at dissipation 0.316","Critical dissipation 0.316 marks two-qubit localization","Strong antiferromagnetic coupling shifts transition to first order","Dissipation 0.316 localizes two-impurity system","First-order transition emerges at strong spin coupling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire location of the transition rests on trusting that the variational trial wave function with 500 linear-grid modes and six coherent states is fully converged, and that discarding solutions whose magnetization does not rise monotonically with coupling removes metastable states rather than genuine ground states.","fun_headline_variants_meta":{"raw":{"variants":["Two qubits localize at dissipation 0.316","Critical dissipation 0.316 marks two-qubit localization","Strong antiferromagnetic coupling shifts transition to first order","Dissipation 0.316 localizes two-impurity system","First-order transition emerges at strong spin coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000429,"raw_usage":{"total_tokens":2228,"prompt_tokens":1015,"completion_tokens":1213,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":1134}},"tokens_in":631,"tokens_out":1213,"duration_ms":9155,"temperature":1.0,"reasoning_tokens":1134,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:17:55.114472+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An unbiased, continuum-extrapolated calculation of the same model at $\\Delta = 0.025$, $K = 0$, $\\varepsilon = 0$ that finds the magnetization jump at $\\alpha_c$ below $0.25$, or a variational run with more coherent states and modes that moves the jump by more than the quoted uncertainty, would falsify the central value. A simpler partial test is the Bethe-ansatz slope: the paper fits $\\Delta/\\langle\\sigma_x\\rangle = (2\\alpha - 1)\\omega_c$ with slope $3.83$ rather than $2$, so a calculation reproducing slope $2$ while keeping the jump near $0.316$ would settle the issue.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"NRG calculation that gives $\\alpha_c \\approx 0.18$ with a small bias field; it is the main baseline the paper argues is underestimated."},{"cited_title":"Zheng, Z","cited_arxiv_id":null,"evidence_quote":"Variational calculation giving $\\alpha_c = 0.125$ under a bias of $10^{-5}\\omega_c$; the paper attributes the low value to the imposed bias."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier Silbey-Harris variational result of $\\alpha_c = 0.5$; it supplies the VM1 comparison and the constrained ansatz the multi-D1 form goes beyond."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Bethe-ansatz and Kosterlitz-Thouless expectations for the Ohmic spin-boson model, including the $\\Delta/\\langle\\sigma_x\\rangle = (2\\alpha - 1)\\omega_c$ relation used as a check."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Variational treatment of the single-impurity Ohmic spin-boson model with coherent-state expansions; it supplies the method and observables the paper extends to two impurities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier application of the multi-D1 ansatz to spin-boson problems; it is the methodological source of the trial wave function."},{"cited_title":"Blunden-Codd, S","cited_arxiv_id":null,"evidence_quote":"Study of bath observables near the sub-Ohmic spin-boson transition; it provides the framework for using bath correlations and variances as transition indicators."}],"review_version":1}