{"id":"e6ee5299-e69e-40b4-ad24-31b59d1a6d1c","arxiv_id":"2607.06850","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"A squeezed-polaron variational ansatz with generic Gaussian boson-boson correlations captures the correct critical scaling, critical exponents, and scattering phase shifts of the sub-Ohmic spin-boson localization transition.","lead":"The paper introduces a variational wavefunction for the spin-boson model that includes Gaussian boson-boson correlations, improving the description of the localization transition in the sub-Ohmic regime. It yields better ground-state energies, critical couplings, and critical exponents than standard coherent-state polaron ansatzes, and predicts boson scattering phase shifts.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Non-mean-field exponents for s>0.5 are derived on a non-ground-state branch, making them physically decoupled from the actual transition","rationale":"The reader's identification of the Gaussian restriction as the weakest assumption is correct, and the CONDITIONAL verdict is appropriate. The concern I raise is a sharpening of the reader's: it's not just that Gaussian states miss cat-like fluctuations in general, but specifically that this causes the non-mean-field exponents to be derived on a non-ground-state branch, making them physically decoupled from the actual transition for s>0.5. The paper's strongest results — the critical scaling laws at b=b_c [Eq. (9)], the critical couplings for s≲0.4, and the scattering phase shifts — are well-supported and represent genuine contributions. The analytical derivation of the integral equation for m(y) and the extraction of scaling laws via Appendix E and G are technically sound. The numerical benchmarks (Fig. 3(a), Fig. S2) show good agreement for the critical couplings, particularly for s≲0.4. The scattering analysis (Fig. 4, Appendices H-I) provides a falsifiable prediction (R_c = sin(πs/2)²) that is independently supported by CFT arguments. The CONDITIONAL verdict correctly reflects that the central claim is partially supported: the deep sub-Ohmic results are quantitatively strong, while the shallow sub-Ohmic non-mean-field exponents are only qualitatively suggestive and physically unrealized on the actual ground-state branch. No verdict adjustment is needed.","tokens_in":27796,"tokens_out":4719,"duration_ms":192900,"concrete_test":"For s=0.8, compute the energy difference E_loc−E_deloc as a function of α near the delocalized branch's critical point (where m₀→0, α_∆≈1.15). If this energy difference is finite and positive at the delocalized critical point — as Fig. 2(d) suggests — then the derived ν−1=1−s does not govern any physical transition, confirming the exponents are branch artifacts. Conversely, if the energy gap closes continuously as m₀→0, the exponents would have physical relevance.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the Gaussian restriction as the root concern, but the specific load-bearing issue is more precise. For s≳0.4, the ansatz produces a spurious first-order transition: the delocalized branch (θ=0) and localized branch (θ≠0) cross at a finite energy gap, and the system jumps discontinuously between them. The non-mean-field exponent ν−1=1−s for s>0.5 is derived (Appendix G) from the continuous vanishing of m₀ on the *delocalized* branch at b=b_c. But for s>0.5, this b=b_c point lies on the higher-energy branch — the localized branch has already become the ground state well before m₀ reaches zero (visible in Fig. 2 for s=0.8). Consequently, the power-law scaling ω*∼|α−α_c|^ν is computed at a point that is not the actual phase transition. The exponents are artifacts of the ansatz's internal structure on a non-physical branch, not predictions of the model's critical behavior. The authors acknowledge this ('in a regime where the localized branch is energetically favored'), but the strongest_claim still states the ansatz 'displays non-mean-field critical exponents in the shallow sub-Ohmic regime' without flagging that these exponents are physically unrealized. The deep sub-Ohmic results (s≲0.4, where b_c coincides with the actual transition) and the scattering phase shifts are solid and well-supported.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"The manuscript introduces a 'squeezed polaron' variational wavefunction for the sub-Ohmic spin-boson model that extends the standard coherent-state polaron ansatz by allowing the bosonic bath to be in a generic Gaussian state. This permits the variational principle to select optimal boson-boson correlations induced by the spin impurity. The central technical result is a self-consistent integral equation [Eq. (5)] for a correlation function m(y), from which the authors derive the correct critical power-law scaling of bath observables [Eq. (9)], critical exponents that transition from mean-field (ν−1=s for s<0.5) to non-mean-field (ν−1=1−s for s>0.5), and critical couplings α_c that compare favorably with numerically exact methods (Fig. 3, Table S2). The paper also derives boson scattering phase shifts from the variational background and shows they encode an emergent energy scale ω* that vanishes at the transition. The variational derivation is internally consistent: the energy functional Eq. (3) is minimized systematically via resolvent methods (Appendix A), the critical scaling is verified both analytically (Appendix E) and numerically (Fig. S3), and the correlation-length exponent is derived from a controlled perturbation analysis of the scaling function (Appendix G).","tokens_in":28501,"tokens_out":1172,"duration_ms":155808,"significance":"The paper makes a genuine contribution to the analytical theory of the spin-boson model. Standard coherent-state polaron ansatzes (Silbey-Harris, Chin) are known to be unable to capture the boson-boson correlations that develop at low frequencies near the localization transition. This work provides the first analytical variational treatment that incorporates generic Gaussian bath correlations and demonstrates that they are sufficient to recover the correct critical scaling laws previously accessible only through diagrammatic or numerical methods. The derivation of scattering phase shifts and their connection to the emergent scale ω* via both variational and CFT methods (Appendix I) is a nice addition that provides experimentally falsifiable predictions. The critical coupling values (Table S2) represent a systematic improvement over the coherent-state polaron and are within a few percent of QMC and VMPS benchmarks across the full sub-Ohmic regime. These are concrete, verifiable results.","major_comments":[{"comment":"Abstract and main-text framing of non-mean-field exponents for s>0.5: The abstract states the ansatz 'displays non-mean-field critical exponents in the shallow sub-Ohmic regime,' and the strongest claims of the paper echo this. However, as the authors themselves acknowledge in the 'Delocalized branch near b_c' section, for s≳0.4 the delocalized branch at b=b_c (where m₀→0 and the exponent ν−1=1−s is derived, Appendix G) is NOT the ground state — the localized branch has already crossed to lower energy. The power-law ω*∼|α−α_c|^ν is thus computed on a non-physical branch. The authors note this ('in a regime where the localized branch is energetically favored'), but the abstract and main-text presentation do not adequately flag that these exponents are physically unrealized at the actual phase transition. This is load-bearing because it is one of the two headline claims about critical exop","section":null}],"minor_comments":[{"comment":"Eq. (3): The overbrace notation for D is introduced inline but the definition is somewhat buried. Making the definition of D more prominent (e.g., as a separate numbered equation) would improve readability.","section":null},{"comment":"Fig. 2: The panel labels (a)-(d) are small and the distinction between solid/dashed blue lines is hard to read in the printed version. Consider enlarging or simplifying.","section":null},{"comment":"Fig. 3(d): The caption mentions 'two loop epsilon expansion [37]' but the reference appears to be Fisher-Ma-Nickel (Ref. 37), which is a one-loop result. Please clarify.","section":null},{"comment":"The notation α_Δ is introduced in Eq. (6) but used earlier in the text. Define it explicitly before its first use.","section":null},{"comment":"In the 'Scattering' section, the effective Hamiltonian Eq. (11) is introduced without much derivation. A brief pointer to Appendix H would help the reader navigate to the justification for easier reference.","section":null},{"comment":"Appendix C: The iteration scheme for solving the integral equation is described, but convergence criteria (relative error threshold, number of iterations) are only briefly mentioned. Stating these more precisely would help reproducibility.","section":null},{"comment":"The phrase 'Taken together as a whole, these equations define D, δE, θ, and α_Δ as parametric functions of b' appears after Eq. (8). The reference to Eq. (8) in the surrounding text is somewhat redundant with the equation itself.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The stress-test concern about the non-mean-field exponents being derived on a non-ground-state branch is valid and is the most important issue to address. However, I assess this as a presentation/framing issue rather than a load-bearing error: the deep sub-Ohmic results (s≲0.4) where b_c coincides with the actual transition are solid, and the authors are transparent about the branch crossing in the body text. The fix is to adjust the abstract and framing to avoid implying the s>0.5 exponents are physical predictions of the transition. The Gaussian restriction is a known limitation of variational approaches and the authors acknowledge it appropriately. I recommend minor revision."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"The core idea here is genuinely new: replace the coherent-state polaron's product-state bath with a generic Gaussian state, which introduces a self-consistent function m(y) governing boson-boson correlations. The Silbey-Harris and Chin ansatzes drop out as the m(y)→1 limit. This is a real extension of the variational polaron program, not a parameter tweak. The deep sub-Ohmic results (s≲0.4) are the strongest part. The critical scaling laws in Eq. (9) match the known frequency and s-dependence from diagrammatic arguments, the critical couplings agree with NRG and QMC benchmarks to within a few percent, and the mean-field exponent ν⁻¹=s is recovered correctly. The scattering phase-shift analysis — including the CFT cross-check giving S(E)=e^{iπs} in the critical regime — is a nice bonus that connects the variational ground state to experimentally accessible probes. The derivation through resolvent methods in the appendices is clean and reproducible. The integral equation is solved numerically with a sensible iteration scheme, and the scaling collapse in Fig. S3 checks out. This is solid technical work. The soft spot is real but bounded. For s≳0.4, the delocalized and localized branches cross at a finite energy gap, producing a spurious first-order jump. The non-mean-field exponent ν⁻¹=1−s for s>0.5 is derived from the continuous vanishing of m₀ on the delocalized branch at b=b_c — but by that point the localized branch has already become the ground state. So these exponents are computed on a non-physical branch and are not predictions of the model's actual critical behavior. The authors acknowledge this honestly in the text, but the abstract and strongest-claim framing still present the non-mean-field exponents as a feature of the ansatz without flagging that they're physically unrealized. The Gaussian restriction also misses cat-like fluctuations near criticality, which is the root cause of the spurious first-order jump. None of this invalidates the deep sub-Ohmic results or the scattering predictions. The paper is for researchers working on variational methods for quantum impurity models and open quantum systems. It deserves a serious referee who can check the integral-equation derivation and the scaling arguments — the technical core is sound even where the physical interpretation oversells. I'd recommend acceptance contingent on toning down the non-mean-field exponent claims to match what the authors actually show in the text.","headline":"Gaussian polaron ansatz captures deep sub-Ohmic critical scaling and scattering physics; non-mean-field exponents for s>0.5 are derived on a non-physical branch","tokens_in":28767,"tokens_out":600,"would_cite":true,"duration_ms":136831,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Squeezed polaron captures bath correlations near quantum localization","keywords":[],"falsifier":"If the true ground state near the transition requires non-Gaussian bath correlations (beyond two-point functions) to correctly capture the critical behavior — which is suggested by the spurious first-order jump and the acknowledged need for cat-state superpositions — then the squeezed polaron's critical exponents for s > 0.5 are artifacts of the Gaussian restriction rather than genuine physical predictions.","tokens_in":27954,"feed_emoji":"🌊","tokens_out":1191,"duration_ms":240456,"temperature":0.7,"pith_summary":"The spin-boson model — a two-level quantum system coupled to a bath of harmonic oscillators — undergoes a localization phase transition when the coupling strength exceeds a critical value. Standard variational wavefunctions dress the spin with coherent-state displacements of the bath but cannot capture the boson-boson correlations that proliferate among low-frequency modes near this transition. This paper introduces a squeezed polaron ansatz that replaces the bare bosonic vacuum with a generic Gaussian state, allowing arbitrary two-point correlations among bath modes to emerge from the variational optimization. The central result is that this single modification recovers the correct critical power-law scaling of bath observables, yields critical exponents that transition from mean-field to non-mean-field character as the bath spectral exponent s crosses 1/2, and produces critical coupling values within a few percent of numerically exact methods across the entire sub-Ohmic regime. The authors further derive boson scattering phase shifts from the variational ground state, showing that a constant phase shift pi*s emerges in the critical frequency window and encodes the emergent energy scale that vanishes at the transition.","feed_headline":"Squeezed polaron captures bath correlations near quantum localization","feed_subtitle":"Gaussian boson-boson correlations in a variational ansatz recover correct critical scaling and predict measurable scattering phase shifts at","key_machinery":"The self-consistent integral equation for m(y) [Eq. (5)], which replaces the trivial m=1 of coherent-state polarons and encodes all bath-bath correlations induced by the spin impurity. The parameter b = alpha_D (cos theta)^{3-s} serves as the control parameter; its critical value b_c where m_0 vanishes determines the transition point.","core_discovery":"The squeezed polaron's key object is the function m(y), determined self-consistently by an integral equation [Eq. (5)], which encodes how the spin impurity modifies low-frequency bath correlations beyond what coherent-state displacements can achieve. When m(y) departs from its trivial value of 1 at low frequencies, it generates the correct power-law scaling of displacement and squeezing amplitudes [Eq. (9)] that standard polaron ansatzes miss entirely. The vanishing of m at zero frequency (m_0 -> 0) signals the localization transition, and the scaling function r(omega/omega*) that interpolates between critical and non-critical regimes yields the correlation-length exponent nu^{-1} = min(s, 1","pith_inferences":["The Gaussian restriction on the bath state means the ansatz captures all two-point but no higher-order correlations. Since the true ground state near criticality develops non-Gaussian cat-like structure, the squeezed polaron may systematically underestimate the critical coupling and distort exponents for s > 0.5 — exactly where the spurious first-order jump appears. A non-Gaussian extension (e.g.,","The connection between the variational phase shift phi = pi*s and the CFT result S(E) = e^{i*pi*s} suggests that the squeezed polaron, despite being a variational approximation, captures exact conformal data in the critical regime for s < 0.5 where the underlying CFT is Gaussian. For s > 0.5, multiparticle production corrections would modify the plateau value of the reflection coefficient.","The integral equation for m(y) has the structure of a self-consistent screening equation, analogous to Dyson equations in many-body theory. This suggests the squeezed polaron can be understood as a variational realization of a self-energy approximation, where m(y) plays the role of a frequency-dependent screening function for the spin-bath interaction."],"forward_implications":["The scattering phase shift phi = pi*s in the critical regime is a directly measurable signature of the localization transition, testable in waveguide-QED or circuit-QED experiments where a bosonic mode scatters off an impurity spin.","Adding Gaussian fluctuations to multi-polaron ansatzes could remove the spurious first-order jump between localized and delocalized branches that appears for s > 0.5, potentially yielding quantitatively accurate critical exponents in the non-mean-field regime.","The squeezed-polaron framework extends naturally to the Ohmic case (s=1) and to finite temperature, where bath correlations are expected to play an equally important role.","The emergence of an energy scale omega* << D that governs low-frequency physics suggests a separation of scales that could be exploited in perturbative or renormalization-group treatments beyond the variational ansatz."],"fun_headline_variants":["Squeezed polaron fixes critical scaling near localization transition","Squeezed polaron ansatz predicts boson scattering phase shifts","Squeezed polaron recovers correct scaling in sub-Ohmic spin-boson","Variational squeezed polaron maps criticality in sub-Ohmic model","Squeezed polaron signals localization via vanishing bath correlations"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The ansatz restricts the bosonic bath to a Gaussian state, meaning all bath correlations are fully determined by two-point functions. Near the localization transition, the true ground state develops non-Gaussian cat-like fluctuations that this restriction cannot represent, which likely causes the spurious first-order jump between branches for s > 0.5 and may bias the critical exponents in the non-mean-field regime.","fun_headline_variants_meta":{"raw":{"variants":["Squeezed polaron fixes critical scaling near localization transition","Squeezed polaron ansatz predicts boson scattering phase shifts","Squeezed polaron recovers correct scaling in sub-Ohmic spin-boson","Variational squeezed polaron maps criticality in sub-Ohmic model","Squeezed polaron signals localization via vanishing bath correlations"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":2186,"prompt_tokens":505,"completion_tokens":1681,"prompt_tokens_details":null},"tokens_in":505,"tokens_out":1681,"duration_ms":70826,"temperature":1.0,"reasoning_tokens":1710,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T00:13:20.159451+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the true ground state near the transition requires non-Gaussian bath correlations (beyond two-point functions) to correctly capture the critical behavior — which is suggested by the spurious first-order jump and the acknowledged need for cat-state superpositions — then the squeezed polaron's critical exponents for s > 0.5 are artifacts of the Gaussian restriction rather than genuine physical predictions.","supporting_citations":[],"review_version":1}