{"id":"baf23d79-0b3c-496c-a323-95ce59b61ff7","arxiv_id":"2608.06357","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a low-density one-dimensional Holstein metal, phonon spectral weight spreads from the bare phonon energy down to zero frequency and above it, with no Kohn-anomaly kink at 2k_F, and a dressed RPA scheme captures the DMRG results semi-quantitatively.","lead":"This paper calculates how lattice vibrations change in a one-dimensional metal with very few charge carriers, using a numerical method called DMRG. It finds the vibration spectrum is spread out in frequency rather than showing the textbook sharp feature, and offers a cheap approximation that reproduces the numerical results.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The no-2k_F-kink and low-energy-weight claims rest on N=80 OBC/center-site DMRG without finite-size scaling; at x=0.05 E_F≈0.025t is below η=0.05, so a larger-N check is needed.","rationale":"Reading the paper in good faith, its goal is to establish, via DMRG, that low-density spinless Holstein metals show phonon spectral weight far from Ω and no Kohn-anomaly kink at 2k_F. For that claim to hold, the DMRG spectra must faithfully represent the thermodynamic limit in momentum and frequency. The least secure condition is exactly the numerical setup of Sec. III A: N=80, open boundaries, center-site approximation, and η=0.05, with no finite-size scaling and no comparison of center versus bulk responses. At x=0.05 the Fermi energy is below the broadening, so the low-energy portion of the claimed spectrum is not resolved in the usual sense. The reader's weakest-assumption statement identifies the same issue, and I agree with it. I do not see a separate internal inconsistency: the vertex-correction negativity is explicitly acknowledged and is not part of the central DMRG claim. The absence of a 2k_F kink is also not surprising by itself, because E_F/Ω is small for these parameters and the Migdal-limit condition is not met. Consequently, the original CONDITIONAL verdict should stand unchanged; the condition should be a targeted finite-size and resolution check before the quantitative low-energy claims are treated as established. The lack of deposited data and code is secondary but reinforces the need for such a check.","tokens_in":12402,"tokens_out":7015,"duration_ms":95127,"concrete_test":"Run the x=0.05, λ=1 case with N=160 (Nf=8) and N=240 (Nf=12), keeping t=Ω=1, bond dimension m=1000, local phonon dimension 8, and η=0.05, and overlay B(q,ω) on the q grid near 2k_F≈0.314 and for ω<0.1t. Repeat the N=160 run with η=0.02 to separate genuine low-energy weight from Lorentzian tails. If the ω→0 tail and the absence of a 2k_F feature persist at both larger sizes and at the smaller broadening, the central claim survives; if the tail collapses or a kink emerges, the DMRG evidence is not yet sufficient.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of the paper is that phonon spectral weight is transferred down to ω=0 and that no 2k_F feature appears near Ω. This claim is supported only by the DMRG spectra of Sec. III A, computed with N=80 open chains, the center-site approximation, up to 8 phonon states, and Lorentzian broadening η=0.05. No finite-size scaling and no center-versus-bulk comparison are presented; the statement that OBC and center-site artifacts are 'minor' is not quantified. The weakest point is the x=0.05 case: the bare Fermi energy E_F≈0.025t is smaller than η=0.05, so the claimed weight at ω→0 lies below the nominal resolution. Additionally, the low-energy particle-hole spectrum near q=2k_F is discretized on an energy scale of order v_F(2π/N)≈0.02t, comparable to the frequency range where the central qualitative feature is claimed. If the apparent low-energy tail and the absence of a 2k_F kink are boundary, center-site, or finite-size artifacts rather than thermodynamic-limit properties, the abstract's strongest claim is unsupported. This is a convergence gap, not an internal inconsistency; the DMRG calculations are reported at fixed parameters and standard truncation accuracies, but the decisive scaling evidence is missing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the phonon spectral function B(q,ω) of the one-dimensional spinless Holstein model at low carrier concentrations x≤0.15 and couplings λ=0.25–1, using DMRG on N=80 open chains and comparing with one-loop approximations. The authors report that spectral weight is transferred both below the bare phonon energy Ω, down to ω=0, and above Ω, and that no Kohn-anomaly kink appears at q=2k_F near Ω. They show that a 'dressed RPA' in which the electron addition propagator is replaced by the Momentum Average (MA) polaronic propagator reproduces the DMRG results semi-quantitatively, while a first-order vertex correction yields negative spectral weight.","tokens_in":12787,"tokens_out":6461,"duration_ms":78250,"significance":"If the DMRG results are converged, this is the first direct numerical determination of phonon spectral functions for a finite-density polaron metal, and the observed broad transfer of spectral weight away from Ω would constitute a clear qualitative distinction from Migdal-limit phenomenology. The dressed RPA scheme is simple, inexpensive, and potentially extendable to higher dimensions, and the paper honestly discusses its limitations, including the bare removal propagator and the absence of a physical explanation for the scheme's success. The DMRG setup is standard and the quoted parameters (m=1000, M=8, 8 phonon states, η=0.05) are explicit, but the absence of finite-size scaling and truncation-convergence data leaves the central claims not fully supported.","major_comments":[{"comment":"The central DMRG claims—that spectral weight is transferred all the way down to ω=0 and that no 2k_F feature appears near Ω—are not backed by convergence evidence. Only N=80 open chains are used with the center-site approximation, and the statement that OBC/center-site artifacts are 'minor' is unquantified. No finite-size scaling (e.g., N=120, 160) and no center-versus-bulk or comparison with periodic-boundary results are shown. At x=0.05, the bare Fermi energy E_F≈0.025t is smaller than the Lorentzian broadening η=0.05, so the claimed weight at ω→0 lies below the nominal energy resolution; the low-energy tail could be a broadening or finite-size artifact. Please provide a convergence analysis demonstrating that the qualitative features persist with increasing N and decreasing η.","section":"Section III A, Fig. 1"},{"comment":"The statement that 'up to 8 phonon states are sufficient to obtain well converged results' is made without supporting data. The paper should show the dependence of B(q,ω) on the local phonon-state truncation (e.g., 6, 8, and 10 states) at λ=1, where multi-phonon dressing is strongest, and also on the bond dimension used specifically in the Krylov correction-vector step (m=1000 is quoted for the ground state, but the spectral calculation may require larger m). Without this, the possibility that the broad spectral features are truncation artifacts cannot be excluded.","section":"Section III A, local phonon basis"},{"comment":"The dressed RPA relies on extracting the quasiparticle weight A(k) and energy E_A(k) from MA Green's functions by fitting Lorentzians of width η to the polaron peak. For x=0.05, the polaron energy scale is of order E_F≈0.025t, comparable to η=0.05, so the fit may be dominated by the artificial broadening rather than the intrinsic line shape. The sensitivity of the dressed RPA spectra to the fitting width should be tested, and comparison with the full MA propagator (including incoherent weight) would strengthen the semi-quantitative agreement claim, which the authors themselves note worsens with increasing x and λ.","section":"Section III B, Eq. (14)"}],"minor_comments":[{"comment":"The emphasis on 'stark contrast with the Kohn-anomaly phenomenology expected in the Migdal limit' is somewhat misleading because E_F/Ω is at most 0.15 here, far outside the Migdal regime; even the bare RPA at the same parameters does not produce a 2k_F kink near Ω since Ω exceeds the particle-hole continuum scale. The genuinely new result is the significant transfer of weight away from Ω, which could be highlighted more precisely.","section":"Abstract and Section III A"},{"comment":"The figure plots DMRG for q<0 and RPA for q>0, relying on B(q,ω)=B(-q,ω); the caption should state explicitly that this symmetry is used.","section":"Figure 1 caption"},{"comment":"No data availability statement or raw data are provided; given the novelty of the DMRG results, making the spectra and the MA fitting parameters available would aid reproducibility.","section":"General"},{"comment":"Equation (A2) is very lengthy; consider moving the full expression to supplemental material and summarizing the derivation steps in the main text.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is interesting and likely correct in its qualitative conclusions, but the strongest claims in the abstract rest entirely on DMRG spectra computed at fixed parameters with no finite-size scaling. The stress-test concern about η > E_F at x=0.05 is real and should be addressed with larger N and smaller η, at least for a representative set of parameters. If the authors provide the missing convergence analysis, the paper could become a solid contribution. I do not see an internal inconsistency or a fundamental error in the dressed RPA construction; the negative-weight vertex result is plausibly explained as an inconsistent truncation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you care about polaron physics or electron-phonon coupling beyond the Migdal limit. The paper does something genuinely new: it computes phonon spectral functions of the 1D spinless Holstein model at finite low doping (x = 0.05–0.15) with DMRG, a regime that had only been studied at the single-polaron level or at half-filling. The central result—that spectral weight moves away from the bare phonon energy, spreading down to ω = 0 and above Ω, with no obvious 2k_F Kohn anomaly—is interesting and, if true, changes how one thinks about phonon renormalization in weakly doped insulators.\n\nThe dressed-RPA construction is a good idea: use the Momentum-Average propagators in the one-loop polarization, and it captures the DMRG trends semi-quantitatively at negligible cost. The authors are honest about its limitations: the hole-removal propagator is left bare, and the quantitative agreement worsens as x and λ grow. They also show, correctly, that adding a naive first vertex correction produces negative spectral weight, and they sensibly attribute this to inconsistent truncation rather than claiming a physical effect. That is careful, transparent work.\n\nThe soft spots are real but not disqualifying. The main one is the absence of finite-size scaling. The DMRG is done on N = 80 open chains with the center-site approximation, and the paper says OBC artifacts are minor but gives no N-dependence or center-versus-bulk comparison. At x = 0.05, the Fermi energy E_F ≈ 0.025t is smaller than the broadening η = 0.05, so the claim that weight extends all the way down to ω = 0 is only resolved at the width of the Lorentzian. The 2k_F kink absence could also be checked more rigorously with larger systems or a momentum-resolved analysis. These are convergence gaps, not internal contradictions; the DMRG parameters are standard and well documented, and the qualitative story is plausible.\n\nI would send this to peer review. The referees should ask for finite-size scaling (at least one larger N, e.g., 120 or 160) and a direct statement about the η versus E_F limit at the lowest doping. With those additions, the central claim would be on solid ground. This is a useful paper for polaron and DMRG practitioners, and the dressed-RPA scheme is likely to be reused.","headline":"New, credible DMRG phonon spectra for the low-doping Holstein model plus a cheap dressed-RPA scheme; the main caveat is missing finite-size scaling, but the paper deserves a serious referee.","tokens_in":13282,"tokens_out":2222,"would_cite":true,"duration_ms":29127,"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":"Low-density polaron metals erase the 2k_F phonon kink.","keywords":["phonon spectral function","Holstein model","polaron metal","DMRG","Kohn anomaly","random phase approximation","Momentum Average approximation","low carrier concentration"],"falsifier":"Compute the phonon spectral function with DMRG on chains of length N=160 and N=320 at x=0.05, λ=1, using a Lorentzian broadening η=0.01 and a phonon cutoff of 12 states, and check whether the spectral weight still reaches ω=0 continuously and whether any feature develops at q=2k_F near Ω; if the low-frequency weight vanishes or a 2k_F kink appears, the paper's central claim fails.","tokens_in":12189,"feed_emoji":"🔬","tokens_out":8094,"duration_ms":82622,"temperature":0.7,"pith_summary":"Using the density matrix renormalization group, this paper computes the phonon spectral function of a one-dimensional spinless Holstein model at low carrier concentrations x ≤ 0.15 and couplings λ = 0.25, 0.5, and 1.0. It reports that electron-phonon coupling moves phonon spectral weight away from the bare phonon energy Ω, spreading it continuously down to ω = 0 and up above Ω, with no sharp feature at momentum 2k_F. This contradicts the Kohn-anomaly picture expected in the Migdal limit, where weight stays near Ω and shows a kink at 2k_F. The authors show that a simple one-loop polarization with propagators dressed by the Momentum Average polaron approximation reproduces the DMRG spectra semi-quantitatively, while adding a first vertex correction yields negative spectral weight. A sympathetic reader would care because low-density doped insulators are a common experimental regime, and this work gives a cheap way to predict their phonon renormalization.","feed_headline":"Low-density polaron metals erase the 2kF phonon kink","feed_subtitle":"In a 1D Holstein model, phonon weight spreads from zero to above Ω, overturning the Kohn-anomaly picture.","key_machinery":"The central object is the retarded phonon propagator and its spectral function B(q,ω) = −(1/π) Im D^R(q,ω). On the numerical side it is computed with DMRG using the root-M Krylov correction-vector method on 80-site open chains under the center-site approximation. On the analytic side the argument runs through the polarization Π(q,ω): the lowest-order particle-hole loop (RPA) is evaluated with bare propagators, then with electron-addition propagators replaced by the coherent polaron peak from the Momentum Average (MA) approximation for a low-density Fermi sea, giving a closed-form 'dressed RPA' expression; a first vertex correction is then appended and shown to produce negative weights. The MA-dressed loop is what carries the main physical claim: it shows that strong renormalization of the electron and hole polaron dispersions, not a Fermi-surface nesting effect, is what reshapes the phonon spectrum.","core_discovery":"The paper's central discovery is that the phonon spectral function B(q,ω) of a low-density polaron metal is not the Migdal-limit Kohn-anomaly spectrum: instead of a sharp peak near the bare phonon energy Ω with a kink at q=2k_F, the DMRG results show spectral weight transferred continuously both below Ω, all the way to ω=0, and above Ω, with no resolvable signature at 2k_F. The same qualitative picture is obtained from a random-phase-approximation polarization loop once the bare electron propagators are replaced by Momentum Average (MA) polaron propagators, and this 'dressed RPA' is semi-quantitative across the studied couplings and densities. Adding the lowest-order vertex correction to this dressed scheme produces unphysical negative spectral weight, an inconsistency the authors trace to mixing diagrammatic orders, and which suggests a Ward-identity-consistent treatment of the vertex is needed for a fully controlled approximation.","pith_inferences":["If the DMRG result survives at smaller broadening and larger system sizes, phonon renormalization in low-density metals reflects polaron-cloud physics rather than Fermi-surface nesting; a test would be measuring phonon line shapes in a lightly doped insulator.","The negative-weight failure of the first vertex correction suggests that iterating the ladder vertex and dressed polarization to self-consistency might restore positive spectral weight, a direct next step the paper leaves open.","At x=0.05 the Fermi energy is smaller than the broadening used, so a finite-size and resolution study with η below the Fermi energy is needed to confirm that the reported weight truly reaches ω=0.","The spinless model deliberately avoids bipolaron formation; testing the same ideas in a spinful Holstein model with repulsion would show whether the absence of the 2k_F kink is generic."],"forward_implications":["Phonon spectral functions in low-density electron-phonon systems with x ≤ 0.15 should not be interpreted with the Migdal-limit Kohn-anomaly picture; the DMRG spectra show weight extending from ω=0 to above Ω with no 2k_F kink.","The dressed RPA with Momentum Average propagators is a negligible-cost improvement over bare RPA that reproduces the DMRG spectra semi-quantitatively, with the agreement worsening only at the largest x and λ.","Adding the lowest-order vertex correction to the dressed scheme yields negative spectral weight, showing that vertex and propagator dressings must be treated consistently when propagators are dressed nonperturbatively.","Because the Momentum Average approximation becomes more accurate in higher dimensions, the dressed RPA scheme should extend to two- and three-dimensional low-density polaron liquids, where no comparable results exist.","A better approximation for the hole-removal part of the fermion propagator should further improve the agreement of the dressed RPA."],"supporting_citations":[{"why":"Supplies the standard coupled electron-phonon propagator equations and Migdal-limit formalism that define the polarization and self-energy used throughout.","marker":"[11]"},{"why":"Shows single-polaron phonon spectral weight transfers below and above Ω, the x→0 precursor of the finite-density spectra.","marker":"[14]"},{"why":"Identifies polaron-induced states below and above the one-phonon continuum, another single-polaron precedent for the weight redistribution.","marker":"[16]"},{"why":"Provides a recent weak-coupling, Ward-identity-consistent vertex-corrected single-polaron calculation used as the comparison for the vertex-correction issues.","marker":"[18]"},{"why":"Offers the only other finite-density phonon spectral function results, for the half-filled Holstein-Hubbard model, used as a contrasting insulating case.","marker":"[19]"},{"why":"Derives the Momentum Average Green's function for the low-density electron-polaron that dresses the RPA polarization.","marker":"[20]"},{"why":"Validates the MA addition propagators against DMRG in this doping and coupling range, justifying the dressed RPA.","marker":"[21]"},{"why":"Supplies the root-M Krylov correction-vector algorithm used to compute the DMRG phonon spectral functions.","marker":"[24]"}],"fun_headline_variants":["Phonon spectral weight smears out, no 2kF kink in polaron metal","DMRG shows phonons escaping the Migdal 2kF kink","Dressed RPA captures low-density polaron phonon spread","No 2kF kink: phonons span from zero to above Ω"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the DMRG spectra computed for 80-site open chains with moderate truncations and a Lorentzian broadening of 0.05 being faithful to the true infinite-system spectral function; if finite-size, truncation, or broadening artifacts are significant, the claimed transfer of weight and the missing 2k_F kink could be numerical artifacts.","fun_headline_variants_meta":{"raw":{"variants":["Phonon spectral weight smears out, no 2kF kink in polaron metal","DMRG shows phonons escaping the Migdal 2kF kink","Dressed RPA captures low-density polaron phonon spread","No 2kF kink: phonons span from zero to above Ω"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000992,"raw_usage":{"total_tokens":4249,"prompt_tokens":1037,"completion_tokens":3212,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":3125}},"tokens_in":653,"tokens_out":3212,"duration_ms":26709,"temperature":1.0,"reasoning_tokens":3125,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:28:35.922501+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the phonon spectral function with DMRG on chains of length N=160 and N=320 at x=0.05, λ=1, using a Lorentzian broadening η=0.01 and a phonon cutoff of 12 states, and check whether the spectral weight still reaches ω=0 continuously and whether any feature develops at q=2k_F near Ω; if the low-frequency weight vanishes or a 2k_F kink appears, the paper's central claim fails.","supporting_citations":[{"cited_title":"Engelsberg and J","cited_arxiv_id":null,"evidence_quote":"Supplies the standard coupled electron-phonon propagator equations and Migdal-limit formalism that define the polarization and self-energy used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows single-polaron phonon spectral weight transfers below and above Ω, the x→0 precursor of the finite-density spectra."},{"cited_title":"Vidmar, J","cited_arxiv_id":null,"evidence_quote":"Identifies polaron-induced states below and above the one-phonon continuum, another single-polaron precedent for the weight redistribution."},{"cited_title":"Phonon spectral function of Holstein polaron: Investigation of many-body effects with self-energy and vertex correction","cited_arxiv_id":"2607.02136","evidence_quote":"Provides a recent weak-coupling, Ward-identity-consistent vertex-corrected single-polaron calculation used as the comparison for the vertex-correction issues."},{"cited_title":"Weber, F","cited_arxiv_id":null,"evidence_quote":"Offers the only other finite-density phonon spectral function results, for the half-filled Holstein-Hubbard model, used as a contrasting insulating case."},{"cited_title":"Berciu, Polarons in spinless metals: a variational solution, Journal of Physics: Materials5, 044002 (2022)","cited_arxiv_id":null,"evidence_quote":"Derives the Momentum Average Green's function for the low-density electron-polaron that dresses the RPA polarization."},{"cited_title":"Nocera and M","cited_arxiv_id":null,"evidence_quote":"Validates the MA addition propagators against DMRG in this doping and coupling range, justifying the dressed RPA."},{"cited_title":"Nocera and G","cited_arxiv_id":null,"evidence_quote":"Supplies the root-M Krylov correction-vector algorithm used to compute the DMRG phonon spectral functions."}],"review_version":1}