REVIEW 3 major objections 4 minor 38 references
Phonon spectral functions of low-density polaron metals
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Low-density polaron metals erase the 2k_F phonon kink.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Section III A, Fig. 1] 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 III A, local phonon basis] 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 III B, Eq. (14)] 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 λ.
minor comments (4)
- [Abstract and Section III A] 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.
- [Figure 1 caption] 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.
- [General] 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.
- [Appendix A] Equation (A2) is very lengthy; consider moving the full expression to supplemental material and summarizing the derivation steps in the main text.
Circularity Check
No significant circularity: the central DMRG result is an independent numerical benchmark, and the dressed-RPA inputs come from a different observable and are validated against DMRG in this paper itself.
full rationale
The paper's central claim, the absence of a 2k_F kink and the transfer of phonon spectral weight down to omega=0, comes directly from DMRG calculations on the Holstein Hamiltonian. The DMRG method computes B(q, omega) from the ground state with no free parameters fitted to the target observable; convergence parameters (m=1000, 8 phonon states, eta=0.05) are numerical settings, not fitted inputs. The RPA result of Eq. (11) is a parameter-free one-loop calculation, and the dressed RPA of Eq. (15) uses quasiparticle weights and energies obtained from the MA electron-addition Green's function, a separate observable, which are then inserted into the polarization loop and compared with DMRG in Fig. 2. That comparison is an external test, not a fit to the phonon spectrum. The self-citations to Refs. [20,21] supply the MA method, but the paper does not rely on them to establish the main qualitative result; its own DMRG data are the evidence. The negative spectral weight found with the vertex correction is an independent outcome. Numerical limitations such as N=80 and finite broadening are convergence concerns, not circularity.
Assumptions & free parameters
free parameters (3)
- Lorentzian broadening eta =
0.05
- Local phonon Hilbert space truncation =
8 phonon states per site
- MA coherent-pole weight A(k) and energy E_A(k) =
Z_k and E_P(k) extracted from Lorentzian fits to MA Green's functions, not tabulated
assumptions (6)
- standard math Dyson equations and spectral representations for G, D, and B
- domain assumption Spinless Holstein model with nearest-neighbor hopping, Einstein phonons, and linear coupling is the relevant model for low-density polaron metals
- domain assumption The MA approximation assumes the metal is frozen in a mean-field Fermi sea plus uniform lattice distortion
- ad hoc to paper In dressed RPA, only the coherent polaron peak of the addition propagator is kept and the removal propagator is left bare
- ad hoc to paper Local phonon Hilbert space can be truncated to 8 states per site
- ad hoc to paper Lowest-order vertex correction can be combined with nonperturbatively dressed propagators
Cite this review
Pith. "Pith review of Phonon spectral functions of low-density polaron metals." pith.science (2026). https://pith.science/paper/OUHFSDCL
@misc{pith2026260806357,
author = {Pith},
title = {Pith review of: Phonon spectral functions of low-density polaron metals},
year = {2026},
howpublished = {\url{https://pith.science/paper/OUHFSDCL}},
note = {Machine review of arXiv:2608.06357}
}
abstract
We use the density matrix renormalization group (DMRG) to compute the phonon spectral function of a one-dimensional spinless Holstein model doped with a low but finite carrier concentration, $x \leq 0.15$, as a function of the electron-phonon coupling $\lambda$. To the best of our knowledge, these are the first such results in this regime, complementing extensive prior work at the single-polaron level ($x\to 0$). We find that significant phonon spectral weight is transferred both below, all the way down to $\omega=0$, and above the bare phonon energy $\Omega$, in stark contrast with the Kohn-anomaly phenomenology expected in the Migdal limit, where weight remains centered near $\Omega$ with a kink at $q=2k_F$. No signature of this $2k_F$ kink appears in our results. This behavior is captured qualitatively by the Random Phase Approximation (RPA), and semi-quantitatively, at negligible extra computational cost, by a ``dressed RPA'' scheme in which the electron addition propagator is renormalized using the Momentum Average (MA) approximation for the low-density electron-polaron. By contrast, adding the lowest-order vertex correction to this dressed scheme produces unphysical negative spectral weight, signaling that vertex and propagator dressings must be treated consistently once the propagators are dressed nonperturbatively. Our results provide an efficient approximation for the phonon spectral functions of low-density polaron metals, a regime relevant to weakly doped insulators.
Figures
Reference graph
Works this paper leans on
-
[1]
L. D. Landau, Electron motion in crystal lattices, Phys. Z. Sow- jetunion3, 664 (1933)
work page 1933
-
[2]
L. D. Landau and S. I. Pekar, The effective mass of the polaron, Zh. Eksp. Teor. Fiz.18, 419 (1948)
work page 1948
-
[3]
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Microscopic the- ory of superconductivity, Phys. Rev.106, 162 (1957)
work page 1957
-
[4]
Bardeen, L
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of super- conductivity, Phys. Rev.108, 1175 (1957)
1957
- [5]
-
[6]
J. Sous, C. Zhang, M. Berciu, D. R. Reichman, B. V . Svistunov, N. V . Prokof’ev, and A. J. Millis, Bipolaronic superconductivity out of a Coulomb gas, Phys. Rev. B108, L220502 (2023)
work page 2023
-
[7]
R. Peierls,Surprises in Theoretical Physics, Princeton Series in Physics (Princeton University Press, Princeton, NJ, USA, 1979)
work page 1979
-
[8]
Kohn, Image of the Fermi surface in the vibration spectrum of a metal, Phys
W. Kohn, Image of the Fermi surface in the vibration spectrum of a metal, Phys. Rev. Lett.2, 393 (1959)
work page 1959
Show all 38 references
-
[9]
Fomichev and M
S. Fomichev and M. Berciu, Renormalized phonon spectrum in the Su-Schrieffer-Heeger model, Journal of Physics: Materials 6, 035003 (2023)
2023
-
[10]
A. B. Migdal, Interaction between electrons and lattice vibra- tions in a normal metal, Sov. Phys. JETP7, 996 (1958), [Zh. Eksp. Teor. Fiz. 34, 1438–1446 (1958)]
1958
-
[11]
Engelsberg and J
S. Engelsberg and J. R. Schrieffer, Coupled electron-phonon system, Phys. Rev.131, 993 (1963)
1963
-
[12]
A. S. Alexandrov,Theory of Superconductivity: From Weak to Strong Coupling(IOP Publishing, Bristol, 2003). 10
2003
-
[13]
A. S. Mishchenko, N. V . Prokof’ev, A. Sakamoto, and B. V . Svistunov, Diagrammatic Monte Carlo method for many- polaron systems, Phys. Rev. B62, 6317 (2000)
2000
-
[14]
J. Loos, M. Hohenadler, A. Alvermann, and H. Fehske, Phonon spectral function of the Holstein polaron, Journal of Physics: Condensed Matter18, 7299 (2006)
2006
-
[15]
O. S. Bari ˇsi´c, Holstein light quantum polarons on the one- dimensional lattice, Phys. Rev. B73, 214304 (2006)
2006
-
[16]
Vidmar, J
L. Vidmar, J. Bon ˇca, and S. A. Trugman, Emergence of states in the phonon spectral function of the Holstein polaron below and above the one-phonon continuum, Phys. Rev. B82, 104304 (2010)
2010
-
[17]
Jansen, J
D. Jansen, J. Bon ˇca, and F. Heidrich-Meisner, Finite- temperature density-matrix renormalization group method for electron-phonon systems: Thermodynamics and Holstein- polaron spectral functions, Phys. Rev. B102, 165155 (2020)
2020
-
[18]
Rai and S
L. Rai and S. Pandey, Phonon spectral function of Holstein po- laron: Investigation of many-body effects with self-energy and vertex correction (2026), arXiv:2607.02136 [cond-mat.str-el]
2026 arXiv
-
[19]
Weber, F
M. Weber, F. F. Assaad, and M. Hohenadler, Phonon spec- tral function of the one-dimensional Holstein-Hubbard model, Phys. Rev. B91, 235150 (2015)
2015
-
[20]
Berciu, Polarons in spinless metals: a variational solution, Journal of Physics: Materials5, 044002 (2022)
M. Berciu, Polarons in spinless metals: a variational solution, Journal of Physics: Materials5, 044002 (2022)
2022
-
[21]
Nocera and M
A. Nocera and M. Berciu, Electron addition spectral functions of low-density polaron liquids, SciPost Phys.15, 110 (2023)
2023
-
[22]
Holstein, Studies of polaron motion: Part I
T. Holstein, Studies of polaron motion: Part I. The molecular- crystal model, Annals of Physics8, 325 (1959)
1959
-
[23]
S. R. White, Density matrix formulation for quantum renormal- ization groups, Phys. Rev. Lett.69, 2863 (1992)
1992
-
[24]
Nocera and G
A. Nocera and G. Alvarez, Root-nKrylov-space correction vec- tors for spectral functions with the density matrix renormaliza- tion group, Phys. Rev. B106, 205106 (2022)
2022
-
[25]
Alvarez, The density matrix renormalization group for strongly correlated electron systems: A generic implementa- tion, Computer Physics Communications180, 1572 (2009)
G. Alvarez, The density matrix renormalization group for strongly correlated electron systems: A generic implementa- tion, Computer Physics Communications180, 1572 (2009)
2009
-
[26]
Nocera and G
A. Nocera and G. Alvarez, Spectral functions with the density matrix renormalization group: Krylov-space approach for cor- rection vectors, Phys. Rev. E94, 053308 (2016)
2016
-
[27]
Jeckelmann and S
E. Jeckelmann and S. R. White, Density-matrix renormalization-group study of the polaron problem in the Holstein model, Phys. Rev. B57, 6376 (1998)
1998
-
[28]
Zhang, E
C. Zhang, E. Jeckelmann, and S. R. White, Density matrix ap- proach to local Hilbert space reduction, Phys. Rev. Lett.80, 2661 (1998)
1998
-
[29]
C. Guo, A. Weichselbaum, J. von Delft, and M. V ojta, Criti- cal and strong-coupling phases in one- and two-bath spin-boson models, Phys. Rev. Lett.108, 160401 (2012)
2012
-
[30]
Brockt, F
C. Brockt, F. Dorfner, L. Vidmar, F. Heidrich-Meisner, and E. Jeckelmann, Matrix-product-state method with a dynamical local basis optimization for bosonic systems out of equilibrium, Phys. Rev. B92, 241106 (2015)
2015
-
[31]
Jansen, C
D. Jansen, C. Jooss, and F. Heidrich-Meisner, Charge density wave breakdown in a heterostructure with electron-phonon cou- pling, Phys. Rev. B104, 195116 (2021)
2021
-
[32]
Stolpp, T
J. Stolpp, T. K ¨ohler, S. R. Manmana, E. Jeckelmann, F. Heidrich-Meisner, and S. Paeckel, Comparative study of state-of-the-art matrix-product-state methods for lattice models with large local Hilbert spaces without U(1) symmetry, Com- puter Physics Communications269, 108106 (2021)
2021
-
[33]
Stolpp, J
J. Stolpp, J. Herbrych, F. Dorfner, E. Dagotto, and F. Heidrich- Meisner, Charge-density-wave melting in the one-dimensional Holstein model, Phys. Rev. B101, 035134 (2020)
2020
-
[34]
Jansen, J
D. Jansen, J. Bon ˇca, and F. Heidrich-Meisner, Finite- temperature optical conductivity with density-matrix renormal- ization group methods for the Holstein polaron and bipolaron with dispersive phonons, Phys. Rev. B106, 155129 (2022)
2022
-
[35]
K ¨ohler, J
T. K ¨ohler, J. Stolpp, and S. Paeckel, Efficient and flexible ap- proach to simulate low-dimensional quantum lattice models with large local Hilbert spaces, SciPost Physics10, 058 (2021)
2021
-
[36]
Mardazad, Y
S. Mardazad, Y . Xu, X. Yang, M. Grundner, U. Schollw ¨ock, H. Ma, and S. Paeckel, Quantum dynamics simulation of in- tramolecular singlet fission in covalently linked tetracene dimer, The Journal of Chemical Physics155, 194101 (2021)
2021
-
[37]
G. D. Mahan,Many-Particle Physics, 3rd ed. (Springer, New York, 2000)
2000
-
[38]
G. L. Goodvin, M. Berciu, and G. A. Sawatzky, Green’s func- tion of the Holstein polaron, Phys. Rev. B74, 245104 (2006)
2006
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.