Pith. sign in

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 →

arxiv 2608.06357 v1 pith:OUHFSDCL submitted 2026-08-06 cond-mat.str-el

classification cond-mat.str-el
keywords phononspectralfunctionHolsteinmodelpolaronmetalDMRGKohnanomalyrandomphaseapproximationMomentumAveragelowcarrierconcentration
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

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)
  1. [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 η.
  2. [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.
  3. [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)
  1. [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.
  2. [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.
  3. [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.
  4. [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

0 steps flagged · score 0.0 of 10

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 3 free parameters · 6 assumptions · 0 invented entities

The central claim depends on the DMRG numerical assumptions (finite chain, phonon truncation, broadening) and, for the approximate scheme, on the MA propagators from prior work. No new particles, forces, dimensions, or conserved quantities are introduced.

free parameters (3)
  • Lorentzian broadening eta = 0.05
    Chosen for DMRG correction-vector and RPA spectra. At x=0.05 the Fermi energy is about 0.025t, smaller than eta, so the low-frequency tail is resolution-limited.
  • Local phonon Hilbert space truncation = 8 phonon states per site
    Chosen after observing convergence. This truncation is a numerical cutoff that the central spectral features depend on, especially at stronger coupling.
  • 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
    Eq. (14) replaces the MA addition propagator by a single coherent pole. These parameters are internal to the dressed RPA scheme and do not affect the DMRG claim.
assumptions (6)
  • standard math Dyson equations and spectral representations for G, D, and B
    Starting point of Section II; standard many-body formalism.
  • domain assumption Spinless Holstein model with nearest-neighbor hopping, Einstein phonons, and linear coupling is the relevant model for low-density polaron metals
    Section I; the spinful model is excluded because of bipolaron instability.
  • domain assumption The MA approximation assumes the metal is frozen in a mean-field Fermi sea plus uniform lattice distortion
    Section III B, citing Ref. [20]; underpins the dressed RPA scheme.
  • 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
    Section III B, Eq. (14); authors state they lack a better removal propagator.
  • ad hoc to paper Local phonon Hilbert space can be truncated to 8 states per site
    Section III A; authors report convergence but show no systematic saturation test across all parameters.
  • ad hoc to paper Lowest-order vertex correction can be combined with nonperturbatively dressed propagators
    Section III C; the result is negative spectral weight, which the authors attribute to inconsistent truncation.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2608.06357 by the authors.

Figure 1
Figure 1. Phonon spectral function B(q, ω) for x = 0.05 (left column), x = 0.10 (middle column) and x = 0.15 (right column) at λ = 0.25 (top row), λ = 0.5 (middle row) and λ = 1 (bottom row). In each panel, the left half (q < 0) shows the DMRG results while the right half (q > 0) shows results obtained from RPA, i.e. keeping only the lowest order polarization loop of Eq. (11). In all cases t = 1, Ω = 1, η = 0.05. a q > 2kF an… view at source ↗
Figure 2
Figure 2. Phonon spectral function B(q, ω) for x = 0.05 (left column), x = 0.10 (middle column) and x = 0.15 (right column) at λ = 0.25 (top row), λ = 0.5 (middle row) and λ = 1 (bottom row). In each panel, the left half (q < 0) shows the DMRG results while the right half (q > 0) shows results obtained from dressed RPA of Eq. (15), see text for more details. In all cases t = 1, Ω = 1, η = 0.05. B. Dressed RPA results The sign… view at source ↗
Figure 3
Figure 3. Phonon spectral function B(q, ω) for x = 0.05 (left column), x = 0.10 (middle column) and x = 0.15 (right column) at λ = 0.25 (top row), λ = 0.5 (middle row) and λ = 1 (bottom row). In each panel, the left half (q < 0) shows the DMRG results while the right half (q > 0) shows results obtained from the polarization that also includes the first vertex correction, see text for more details. In all cases t = 1, Ω = 1, η… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

38 extracted references · 33 canonical work pages

  1. [1]

    L. D. Landau, Electron motion in crystal lattices, Phys. Z. Sow- jetunion3, 664 (1933)

  2. [2]

    L. D. Landau and S. I. Pekar, The effective mass of the polaron, Zh. Eksp. Teor. Fiz.18, 419 (1948)

  3. [3]

    Bardeen, L

    J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Microscopic the- ory of superconductivity, Phys. Rev.106, 162 (1957)

  4. [4]

    Bardeen, L

    J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of super- conductivity, Phys. Rev.108, 1175 (1957)

  5. [5]

    Zhang, J

    C. Zhang, J. Sous, D. R. Reichman, M. Berciu, A. J. Millis, N. V . Prokof’ev, and B. V . Svistunov, Bipolaronic high-temperature superconductivity, Phys. Rev. X13, 011010 (2023)

  6. [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)

  7. [7]

    Peierls,Surprises in Theoretical Physics, Princeton Series in Physics (Princeton University Press, Princeton, NJ, USA, 1979)

    R. Peierls,Surprises in Theoretical Physics, Princeton Series in Physics (Princeton University Press, Princeton, NJ, USA, 1979)

  8. [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)

Show all 38 references
  1. [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)

  2. [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)]

  3. [11]

    Engelsberg and J

    S. Engelsberg and J. R. Schrieffer, Coupled electron-phonon system, Phys. Rev.131, 993 (1963)

  4. [12]

    A. S. Alexandrov,Theory of Superconductivity: From Weak to Strong Coupling(IOP Publishing, Bristol, 2003). 10

  5. [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)

  6. [14]

    J. Loos, M. Hohenadler, A. Alvermann, and H. Fehske, Phonon spectral function of the Holstein polaron, Journal of Physics: Condensed Matter18, 7299 (2006)

  7. [15]

    O. S. Bari ˇsi´c, Holstein light quantum polarons on the one- dimensional lattice, Phys. Rev. B73, 214304 (2006)

  8. [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)

  9. [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)

  10. [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]

  11. [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)

  12. [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)

  13. [21]

    Nocera and M

    A. Nocera and M. Berciu, Electron addition spectral functions of low-density polaron liquids, SciPost Phys.15, 110 (2023)

  14. [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)

  15. [23]

    S. R. White, Density matrix formulation for quantum renormal- ization groups, Phys. Rev. Lett.69, 2863 (1992)

  16. [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)

  17. [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)

  18. [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)

  19. [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)

  20. [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)

  21. [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)

  22. [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)

  23. [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)

  24. [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)

  25. [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)

  26. [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)

  27. [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)

  28. [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)

  29. [37]

    G. D. Mahan,Many-Particle Physics, 3rd ed. (Springer, New York, 2000)

  30. [38]

    G. L. Goodvin, M. Berciu, and G. A. Sawatzky, Green’s func- tion of the Holstein polaron, Phys. Rev. B74, 245104 (2006)

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.