{"id":"b71c3afa-b70d-4490-a1ac-bd46d877216a","arxiv_id":"2608.04816","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Ghost-RISB is generalized to electron-phonon systems and reproduces DMFT benchmarks on the Hubbard-Holstein model, showing that superconducting order is suppressed by a phonon-overlap (Franck-Condon) factor in the bipolaronic regime.","lead":"This paper extends a fast variational method, ghost-RISB, to models where electrons interact with lattice vibrations, and tests it on the Hubbard-Holstein model. It reports that with a few auxiliary ghost orbitals the method matches slower DMFT calculations for quasiparticle weights and phonon distributions, and then uses the speedup to analyze superconductivity in the bipolaronic regime.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central new claim—the Franck-Condon suppression of superconductivity—rests on Figs. 5-6, which lack B and phonon-cutoff convergence checks and are not cross-checked against DMFT; the observed drop in Δ and I_dh could be a truncation artifact.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing premise: the finite-ghost/finite-phonon embedding must be converged and accurate in the superconducting regime. I agree. The normal-state benchmarks, while encouraging, do not cover the strong-coupling adiabatic parameters used for the mechanism claim, and the paper itself flags the need for further symmetry-broken-phase studies. The most direct route to settle the concern is a convergence study in B and N_ph, plus an optional DMFT cross-check. Given the absence of such checks, the correct verdict is CONDITIONAL: the method and mechanism are plausible, but the central new physical claim is not yet fully supported.","tokens_in":21937,"tokens_out":7095,"duration_ms":69659,"concrete_test":"Recompute the superconducting calculations of Fig. 5 (U=0, ω0=1.0, λ=0 to 2) and Fig. 6 (U=1, ω0=0.5, λ=0 to 2.5) for B=5 and phonon truncations N_ph=60 and 80, and overlay Δ(λ) and I_dh(λ) on the existing B=1/B=3 curves. If the drop in Δ and I_dh shifts by more than a few percent or vanishes, the mechanism is a truncation artifact. Optionally, benchmark Δ(λ) against DMFT (e.g., EDIpack) for the Holstein model (U=0) at the same parameters to validate the order-parameter magnitude.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the convergence of the superconducting-state calculations in Section III.B. The paper's central mechanism—that the drop of the order parameter Δ tracks the doublon-holon overlap I_dh (Eq. 38) due to bipolaron formation—is based on Figs. 5 and 6, but these figures use only B=1 and B=3 ghost orbitals, with no B=5 data, and the phononic Hilbert space is truncated to 40 states ('We typically truncate the phononic Hilbert space to 40 states', Section III). No convergence study versus B or N_ph is reported for these parameters (U=0, ω0=1.0 and U=1, ω0=0.5, λ up to 2-2.5). The authors themselves state, after Fig. 5, that 'further investigations of various symmetry breakings are in order to validate and/or disprove this statement' (the claim that B=1 suffices in broken-symmetry phases). The normal-state benchmarks in Figs. 1-4 cover different parameter regimes and do not establish accuracy in the strong-coupling adiabatic superconducting regime. If the drop in Δ and I_dh were to shift or disappear upon increasing B to 5 or the phonon cutoff to 80, the Franck-Condon mechanism would be an artifact of the finite bath and phonon truncation. Since the DMFT benchmarks in the superconducting phase are absent, there is no independent check that the variational embedding accurately captures the bipolaronic suppression.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper generalizes the ghost-RISB (ghost-Gutzwiller) formalism to include local dispersionless phonons coupled to arbitrary on-site electronic operators, encoding the phonons in the slave-boson amplitudes. It applies the method to the Hubbard-Holstein model on the Bethe lattice, benchmarking electron quasiparticle weights, mass renormalization, and phonon probability distribution functions against published DMFT data (Figs. 1-4). It then studies s-wave superconductivity in the strong-coupling adiabatic regime (Figs. 5-6), reporting a drop of the order parameter Δ that is interpreted via the inequality |Δ| ≤ I_dh, where I_dh is the overlap of the charge-resolved phonon distributions in the holon and doublon sectors (Eq. (38), Appendix B). The paper identifies a Franck-Condon-like suppression of doublon-holon coherence as the mechanism for this drop.","tokens_in":22269,"tokens_out":6180,"duration_ms":71034,"significance":"If the convergence of the superconducting-state calculations can be established, the approach would provide an efficient, variational, and nonperturbative tool for correlated electron-phonon systems, and Eq. (38) would be a clean, parameter-free relation connecting a phononic observable to the pairing amplitude. The manuscript's strengths include the absence of parameters fitted to the DMFT benchmarks, the use of external DMFT data in the normal state, and the explicit inequality derived in Appendix B. The main open questions are whether the finite ghost-orbital and finite-phonon-basis calculations in the bipolaronic superconducting regime are converged and whether the analytical self-energy expression imported from Ref. [60] remains valid when phonons are encoded in the slave-boson fields.","major_comments":[{"comment":"The central new result, the Franck-Condon mechanism for the suppression of Δ, is based on B=1 and B=3 calculations with the phonon Hilbert space truncated to 40 states; no convergence study with respect to B or the phonon cutoff is reported for these parameters. The authors' own caveat after Fig. 5 ('further investigations of various symmetry breakings are in order') indicates this is a known gap. Please add B=5 (and if possible B=7) and larger phonon cutoffs for the parameter sets of Figs. 5 and 6, and ideally an independent DMFT or other benchmark in the bipolaronic regime. Without this, the drop in Δ and I_dh could be a truncation artifact.","section":"Section III.B, Figs. 5 and 6"},{"comment":"The analytical self-energy expression is imported from Ref. [60], which is a purely electronic ghost-RISB result. Its use here assumes that the form remains valid when phonons are encoded in the slave-boson fields and enter the embedding Hamiltonian. Since the quasiparticle-weight benchmarks in Figs. 1 and 2 rely on Eq. (33), this assumption is load-bearing. Please provide a derivation or explicit justification of Eq. (33) in the electron-phonon setting.","section":"Section II and Section III.A, Eq. (33)"},{"comment":"Appendix B proves only the inequality |Δ| ≤ I_dh; the claim that the suppression of Δ is caused by the drop of I_dh requires more than a coincident decrease of an upper bound. The argument would be conclusive if the ratio |Δ|/I_dh were shown to remain nearly constant before the bipolaronic crossover and to drop only there, or if an additional wavefunction diagnostic directly isolated the doublon-holon coherence. As written, the tracking is suggestive but does not uniquely identify the mechanism.","section":"Appendix B and Section III.B"}],"minor_comments":[{"comment":"The caption of Fig. 6 states ω0/D=1.0, while the text in Section III.B states ω0/D=0.5 for the same calculation; please reconcile this discrepancy.","section":"Figure 6"},{"comment":"The paper states that the phonon Hilbert space is 'typically' truncated to 40 states, but it does not state the cutoff used in each figure or report a convergence test in the normal state; please specify the values and any convergence checks.","section":"Section III"},{"comment":"The abstract and conclusions claim a 'fraction of the computational cost' of DMFT, but no runtime or scaling data are reported; a quantitative comparison would make the practical claim concrete.","section":"Abstract and Conclusions"},{"comment":"The formula 'M=M/2(B−1)' reuses the symbol M on both sides and is difficult to parse; please use a distinct notation for the two quantities.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies on Ref. [60], an arXiv preprint by the same first author, for Eq. (33); if this result is not yet peer-reviewed, the editor may wish to have its correctness independently verified. For a methods paper, releasing data or code for the benchmarks would substantially strengthen reproducibility, especially because several DMFT comparisons are digitized from earlier publications."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know this paper is the first ghost-RISB formulation with local phonons, and the normal-state benchmarks are genuinely good. The authors extend the ghost-orbital construction to include bosonic degrees of freedom in the slave-boson fields, and they show that B=5 matches published DMFT results for quasiparticle weights and mass renormalization across a range of U, filling, and phonon frequency (Figs. 1-2), and B=3 matches phonon PDFs and charge-resolved PDFs (Figs. 3-4). That is a real quantitative advance over earlier Gutzwiller phonon schemes, which had no ghost orbitals and no frequency-dependent self-energy. The upper bound |Δ| ≤ I_dh, derived in Appendix B, is a clean parameter-free result and is probably the most durable piece of the paper.\n\nThe soft spots are in the superconducting section. The mechanism claim—that the drop of Δ tracks the Franck-Condon overlap I_dh—rests on Figs. 5-6, which use only B=1 and B=3, with no B=5, and a phonon truncation of 40 states, with no convergence study. There is no DMFT cross-check in that regime. The authors themselves note that further investigation of broken-symmetry phases is needed. So the suppression of superconductivity is plausible but not yet established; it could still be a truncation artifact, and the paper does not rule that out.\n\nTwo smaller issues: the self-energy formula Eq. (33) is imported from Ref. [60], an unreviewed arXiv preprint by the same first author, and its validity with phonons encoded in the slave-boson fields is assumed rather than demonstrated. And the Figure 6 caption says ω0/D=1.0 while the text says ω0/D=0.5; the text's claim of a 'more adiabatic' regime depends on which is right. No code or data is shipped, and the DMFT comparisons use digitized literature curves, which limits precision but isn't disqualifying.\n\nNone of this undermines the core methodological advance. The paper is internally consistent, the derivations are careful, and the normal-state benchmarks are the right kind of evidence. The SC mechanism is the one place where the evidence is thinner than the claims.\n\nI'd send it to a serious referee. The method deserves exposure, and the convergence question is answerable. With a convergence study and a fixed caption, this could be a strong PRB-type paper.","headline":"A real method extension with convincing normal-state benchmarks, but the headline superconducting mechanism is not yet backed by convergence checks.","tokens_in":22821,"tokens_out":4334,"would_cite":true,"duration_ms":40687,"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":"Ghost-RISB matches DMFT for electrons plus phonons at a fraction of the cost.","keywords":["ghost-RISB","Hubbard-Holstein model","electron-phonon coupling","slave-boson methods","dynamical mean-field theory","quasiparticle weight","polaronic crossover","bipolaronic superconductivity"],"falsifier":"Perform a numerically exact calculation (for example, DMFT with a substantially larger bath than $B=5$ and more than 40 phonon states) of the superconducting order parameter and of the charge-resolved phonon PDFs for the Hubbard-Holstein model at $U/D=1$, $\\omega_0/D=0.5$, $n=0.5$ for $\\lambda/D$ between 1.0 and 2.5; if the exact $|\\Delta|$ does not fall when $I_{dh}$ falls, or if ghost-RISB with $B=3$ to $5$ disagrees with the exact values, the central claim fails.","tokens_in":21707,"feed_emoji":"⚛️","tokens_out":12464,"duration_ms":131686,"temperature":0.7,"pith_summary":"This paper extends the ghost-RISB variational method—a slave-boson approach that adds auxiliary ghost orbitals to mimic dynamical self-energy effects—to systems where electrons are coupled to local phonons. Benchmarking on the Hubbard-Holstein model, it shows that with three to five ghost orbitals the method reproduces DMFT quasiparticle weights, mass renormalization, and phonon displacement distributions across a wide range of couplings, fillings, and phonon frequencies. It then uses the low-cost method to study strong-coupling, adiabatic-phonon superconductivity, where it finds that the superconducting order parameter is bounded by the overlap of the phonon wavefunctions attached to empty and doubly occupied sites, and that this overlap collapses as bipolarons form. The paper presents this Franck-Condon-like reduction of doublon-holon coherence as the mechanism behind the suppression of superconductivity in the bipolaronic regime.","feed_headline":"Ghost orbitals make slave-boson method match DMFT for phonons","feed_subtitle":"A cheaper variational method matches DMFT on quasiparticles and exposes why bipolarons suppress superconductivity.","key_machinery":"The central object is the electron-phonon ghost-RISB embedding: a local impurity problem whose slave-boson amplitudes carry both electronic and phononic quantum numbers, so the phonon mode enters the variational state directly rather than through a perturbative or static renormalization. The identity that carries the mechanism argument is Eq. (38), $|\\Delta| \\le I_{dh} = \\int dx \\sqrt{P_0(x)} \\sqrt{P_2(x)}$, an upper bound on the s-wave superconducting order parameter expressed as the overlap (Bhattacharyya coefficient) of the square roots of the charge-resolved phonon displacement distributions for empty and doubly occupied sites. This bound converts an abstract many-body cancellation into a directly computable phonon quantity: when bipolarons form and $P_0$ and $P_2$ separate, $I_{dh}$ drops and drags the order parameter down with it. The ghost orbitals themselves do the work of reproducing the dynamical self-energy through the analytical expression Eq. (33), which is what makes the quasiparticle-weight benchmarks pass.","core_discovery":"On the paper's own terms, the central discovery is that ghost-RISB remains reliable when local phonons are encoded directly into the slave-boson fields: the embedding Hamiltonian solved by exact diagonalization contains the full local electron-electron, electron-phonon, and phonon terms, and the ghost orbitals supply the frequency structure of the self-energy that static slave-boson methods miss. Benchmark comparisons show that $B=3$ to $5$ ghost orbitals already give DMFT-level results for quasiparticle weights and phonon probability distributions, including the polaronic crossover, while $B=1$ (standard RISB) fails in the adiabatic regime. In the superconducting phase, the paper proves the bound $|\\Delta| \\le I_{dh} = \\int dx \\sqrt{P_0(x)} \\sqrt{P_2(x)}$, where $P_0$ and $P_2$ are the phonon displacement distributions on empty and doubly occupied sites, and demonstrates numerically that $I_{dh}$ tracks the order parameter and falls when the phonon distribution becomes bimodal. The authors interpret this as a phonon-induced orthogonality between holon and doublon sectors—a Franck-Condon effect—that suppresses superconducting coherence even while local pairing remains large.","pith_inferences":["The bound $|\\Delta| \\le I_{dh}$ is derived from the embedding ground state and does not rely on half-filling or on $U=0$, so it should hold in any local-pairing model; computing $I_{dh}$ from charge-resolved phonon PDFs is therefore a cheap general diagnostic for bipolaronic suppression of pairing coherence.","If the Franck-Condon picture is right, increasing the phonon frequency (or otherwise stiffening the lattice) should keep $P_0$ and $P_2$ overlapping and restore superconducting coherence at fixed local pairing—a specific prediction the paper does not make.","The same encoding of phonons inside the slave-boson fields generalizes to several local phonon modes and multiorbital impurities, so the method should be directly applicable to realistic materials where a full dynamical impurity solver is too costly; this is a natural but unproven next step.","The overlap bound might also constrain other off-diagonal orders, such as charge-density-wave or excitonic condensates, whenever the order parameter connects two phonon-dressed charge sectors; whether it is tight in those cases is unexplored."],"forward_implications":["With $B=5$ ghost orbitals, ghost-RISB reproduces DMFT quasiparticle weights and mass renormalizations for the Hubbard-Holstein model, so correlated electron-phonon phase diagrams can be scanned at a fraction of DMFT's cost.","The $B=1$ static slave-boson limit misses the polaronic crossover in the presence of electron-electron repulsion, so a truly dynamical treatment of phonons—here supplied by ghost orbitals—is required in the adiabatic regime.","Because the self-consistency uses static observables, the method avoids DMFT's expensive dynamical self-consistency and can reach strong-coupling, low-phonon-frequency parameters where DMFT becomes unfeasible.","Equation (38) makes the suppression of superconductivity in the bipolaronic regime a quantitative statement about phonon wavefunction overlap, not just about mass renormalization or the effective pairing interaction.","In the doped, more adiabatic case the bound drops faster, producing a sharper suppression of the order parameter even though local double occupations remain numerous."],"supporting_citations":[{"why":"Supplies the DMFT quasiparticle-weight data $Z(\\lambda)/Z(0)$ that ghost-RISB is benchmarked against in Fig. 1.","marker":"[29]"},{"why":"Supplies the DMFT mass-renormalization data as a function of phonon frequency and filling used in Fig. 2.","marker":"[30]"},{"why":"Provides the DMFT phonon displacement PDFs and the polaron-to-bipolaron crossover interpretation used in Figs. 3-4 and for the superconducting mechanism.","marker":"[69]"},{"why":"Additional DMFT phonon PDF data extracted for the polaron crossover benchmark in Fig. 3.","marker":"[72]"},{"why":"Defines the ghost-RISB operatorial formalism that this work extends to include local phonons.","marker":"[42]"},{"why":"Gives the analytical self-energy expression Eq. (33) and the ghost-RISB-to-DMFT correspondence used to justify quasiparticle-weight accuracy.","marker":"[60]"},{"why":"Reports within DMFT the decay of the superconducting order parameter in the polaronic regime, the behavior ghost-RISB reproduces and explains.","marker":"[77]"},{"why":"Documents the accuracy of ghost-RISB as a function of impurity bath size, supporting the claim of DMFT-level accuracy with a smaller bath.","marker":"[45]"}],"fun_headline_variants":["Ghost-RISB brings DMFT accuracy to electron-phonon models","Slave-boson with ghosts outmatches static theory for phonons","New variational method matches DMFT on electron-phonon competition","Ghost orbitals explain superconductivity drop via Franck-Condon effect","Efficient ghost-RISB pinpoints phonon-induced pairing loss"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method's accuracy rests on the assumption that a few ghost orbitals and a phonon Hilbert space truncated to about 40 states capture the full dynamical effect of the lattice on the electrons, especially in the strong-coupling adiabatic and superconducting regimes where the paper does not show a direct comparison with a more exact method.","fun_headline_variants_meta":{"raw":{"variants":["Ghost-RISB brings DMFT accuracy to electron-phonon models","Slave-boson with ghosts outmatches static theory for phonons","New variational method matches DMFT on electron-phonon competition","Ghost orbitals explain superconductivity drop via Franck-Condon effect","Efficient ghost-RISB pinpoints phonon-induced pairing loss"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000758,"raw_usage":{"total_tokens":3430,"prompt_tokens":1072,"completion_tokens":2358,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":688,"completion_tokens_details":{"reasoning_tokens":2267}},"tokens_in":688,"tokens_out":2358,"duration_ms":18360,"temperature":1.0,"reasoning_tokens":2267,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:46:57.063152+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform a numerically exact calculation (for example, DMFT with a substantially larger bath than $B=5$ and more than 40 phonon states) of the superconducting order parameter and of the charge-resolved phonon PDFs for the Hubbard-Holstein model at $U/D=1$, $\\omega_0/D=0.5$, $n=0.5$ for $\\lambda/D$ between 1.0 and 2.5; if the exact $|\\Delta|$ does not fall when $I_{dh}$ falls, or if ghost-RISB with $B=3$ to $5$ disagrees with the exact values, the central claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the DMFT quasiparticle-weight data $Z(\\lambda)/Z(0)$ that ghost-RISB is benchmarked against in Fig. 1."},{"cited_title":"Barone, R","cited_arxiv_id":null,"evidence_quote":"Provides the DMFT phonon displacement PDFs and the polaron-to-bipolaron crossover interpretation used in Figs. 3-4 and for the superconducting mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Additional DMFT phonon PDF data extracted for the polaron crossover benchmark in Fig. 3."},{"cited_title":"Coulter and A","cited_arxiv_id":null,"evidence_quote":"Defines the ghost-RISB operatorial formalism that this work extends to include local phonons."},{"cited_title":"Mejuto-Zaera and M","cited_arxiv_id":null,"evidence_quote":"Gives the analytical self-energy expression Eq. (33) and the ghost-RISB-to-DMFT correspondence used to justify quasiparticle-weight accuracy."},{"cited_title":"Lechermann, A","cited_arxiv_id":null,"evidence_quote":"Documents the accuracy of ghost-RISB as a function of impurity bath size, supporting the claim of DMFT-level accuracy with a smaller bath."}],"review_version":1}