Including the effect of lattice anharmonicity on electron-phonon interactions has recently garnered attention due to its role as a necessary and significant component in explaining various phenomena, including superconductivity, optical response, and temperature dependence of mobility. This study focuses on analytically treating the effects of anharmonic electron-phonon coupling on the polaron self-energy, combined with numerical Diagrammatic Monte Carlo (DiagMC) data. Specifically, we incorporate a quadratic interaction into the method of squeezed phonon states, which has proven effective for analytically calculating the polaron parameters. Additionally, we extend this method to nonparabolic finite-width conduction bands while maintaining the periodic translation symmetry of the system. Our results are compared with those obtained from Diagrammatic Monte Carlo, partially reported in a recent study [S. Ragni et al., Phys. Rev. B 107, L121109 (2023)], covering a wide range of coupling strengths for the nonlinear interaction. Remarkably, our analytic method predicts the same features as the Diagrammatic Monte Carlo simulation.
Including the effect of lattice anharmonicity on electron-phonon interactions has recently garnered attention due to its role as a necessary and significant component in explaining various phenomena, including superconductivity, optical response, and temperature dependence of mobility. This study focuses on analytically treating the effects of anharmonic electron-phonon coupling on the polaron self-energy, combined with numerical Diagrammatic Monte Carlo (DiagMC) data. Specifically, we incorporate a quadratic interaction into the method of squeezed phonon states, which has proven effective for analytically calculating the polaron parameters. Additionally, we extend this method to nonparabolic finite-width conduction bands while maintaining the periodic translation symmetry of the system. Our results are compared with those obtained from Diagrammatic Monte Carlo, partially reported in a recent study [S. Ragni et al., Phys. Rev. B 107, L121109 (2023)], covering a wide range of coupling strengths for the nonlinear interaction. Remarkably, our analytic method predicts the same features as the Diagrammatic Monte Carlo simulation.
In the present work, the problem of an all-coupling analytic description for the optical conductivity of the Fr\"ohlich polaron is treated, with the goal being to bridge the gap in the validity range that exists between two complementary methods: on the one hand, the memory-function formalism and, on the other hand, the strong-coupling expansion based on the Franck-Condon picture for the polaron response. At intermediate coupling, both methods were found to fail as they do not reproduce diagrammatic quantum Monte Carlo results. To resolve this, we modify the memory-function formalism with respect to the Feynman-Hellwarth-Iddings-Platzman approach in order to take into account a nonquadratic interaction in a model system for the polaron. The strong-coupling expansion is extended beyond the adiabatic approximation by including in the treatment nonadiabatic transitions between excited polaron states. The polaron optical conductivity that we obtain at T=0 by combining the two extended methods agrees well, both qualitatively and quantitatively, with the diagrammatic quantum Monte Carlo results in the whole available range of the electron-phonon coupling strength.
We develop and compare several analytical approximations for the ground state and zero-temperature polaron dispersion in finite-width, nonparabolic conduction bands. The main focus of the work is an extension of the Feynman variational method to a tight-binding lattice, where the effective-mass approximation is no longer applicable. The resulting variational formulation is not restricted to a specific phonon dispersion or electron-phonon interaction and provides a uniform description across weak-, intermediate-, and strong-coupling regimes. In addition, we revisit and generalize other analytical approaches traditionally formulated for continuum polarons, including canonical transformations and self-consistent Wigner-Brillouin-type approximations. For lattice polarons, these methods exhibit qualitative features absent in the continuum case, such as a nontrivial connection between weak- and strong-coupling limits. We show that an improved Wigner-Brillouin scheme yields a momentum-dependent polaron self-energy free of resonances and in good agreement with numerically exact results. All methods are applied to the Holstein model on a tight-binding lattice and are benchmarked against numerically exact calculations, including diagrammatic Monte Carlo (both our calculations and preceding works) and exact diagonalization results. Furthermore, the analytical approaches are extended to polarons with Rashba-type spin-orbit coupling, providing a stringent test of their applicability in systems with nontrivial band structure. Our results demonstrate that the modified Feynman variational method yields ground-state energies and dispersions with accuracy comparable to, and in many cases exceeding, that of other established analytical approaches. The developed framework offers a versatile and reliable analytical description of lattice polarons beyond the continuum approximation.
The basic principles behind a path integral approach to the lattice polaron are reviewed. Analytical integration of phonons reduces the problem to one self-interacting imaginary-time path, which is then simulated by Metropolis Monte Carlo. Projection operators separate states of different symmetry, which provides access to various excited states. Shifted boundary conditions in imaginary time enable calculation of the polaron mass, spectrum and density of states. Other polaron characteristics accessible by the method include the polaron energy, number of excited phonons and isotope exponent on mass. Monte Carlo updates are formulated in continuous imaginary time on infinite lattices and as such provide statistically unbiased results without finite-size and finite time-step errors. Numerical data are presented for models with short-range and long-range electron?phonon interactions.
Polarons—quantum particles dressed by a cloud of bath excitations—arise whenever a particle interacts strongly with its host medium. Electron and hole polarons, formed when carriers bind to a phonon cloud through electron-phonon (e-ph) interactions, govern charge transport in transition metal oxides, organic semiconductors, and halide perovskites. The concept generalizes to exciton-phonon interactions, giving rise to exciton polarons that produce the large Stokes shifts and broad photoluminescence (PL) lineshapes central to light-emitting applications. A predictive treatment of polarons has remained elusive. Numerically exact methods such as diagrammatic Monte Carlo (DMC) have been restricted to simplified model Hamiltonians, while first-principles approaches based on density functional theory provide material-specific e-ph interactions but treat the e-ph correlation approximately. This thesis bridges that gap by developing first-principles diagrammatic Monte Carlo (FPDMC), which stochastically sums all e-ph Feynman diagrams using accurate ab-initio interactions, yielding numerically exact polaron properties in real materials. Two methodological advances make FPDMC possible. A matrix-product formalism sums over electronic band indices exactly, alleviating the multi-band sign problem in realistic band structures. A data-driven compression of e-ph interactions based on singular-value decomposition (SVD) in the Wannier basis then accelerates the method by more than 10³×. The SVD compression also reveals a hidden low-dimensional structure of quantum interactions, maintaining accuracy across charge transport, spin relaxation, band renormalization, and superconductivity calculations. Extending this philosophy to phonon-phonon interactions via a permanent CP (PCP) decomposition yields 10³–10⁴× compression of anharmonic force constants and orders-of-magnitude speedups in thermal conductivity calculations. Applying FPDMC, we study polarons in LiF, SrTiO₃, and TiO₂, finding excellent agreement with experimental mobilities and resolving the contrasting small- and large-polaron transport in the two TiO₂ polymorphs. To make the microscopic picture concrete, we introduce an e-ph correlation function that visualizes lattice distortion and polaron orbital character. Extending FPDMC to exciton polarons, we obtain first-principles exciton-phonon interactions from finite-momentum Bethe-Salpeter equation calculations and derive a Fermi-golden-rule PL theory. Applied to CrBr₃, FPDMC predicts a broad PL spectrum agreeing with experiment in both peak position and lineshape.
We consider two large polaron systems that are described by a Fröhlich type of Hamiltonian, namely the Bose–Einstein condensate (BEC) polaron in the continuum and the acoustic polaron in a solid. We present ground-state energies of these two systems calculated with the Diagrammatic Monte Carlo (DiagMC) method and with a Feynman all-coupling approach. The DiagMC method evaluates up to very high order a diagrammatic series for the polaron's self-energy. The Feynman all-coupling approach is a variational method that has been used for a wide range of polaronic problems. For the acoustic and BEC polaron both methods provide remarkably similar non-renormalized ground-state energies that are obtained after introducing a finite momentum cutoff. For the renormalized ground-state energies of the BEC polaron, there are relatively large discrepancies between the DiagMC and the Feynman predictions. These differences can be attributed to the renormalization procedure for the contact interaction.
Many condensed matter problems involve a particle coupled to its environment. The polaron, originally introduced to describe electrons in a polarizable medium, describes a particle coupled to a bosonic field. The Holstein polaron model, although simple, including only optical Einstein phonons and an interaction that couples them to the electron density, captures almost all of the standard polaronic properties. We herein investigate polarons that differ significantly from this behaviour. We study a model with phonon-modulated hopping, and find a radically different behaviour at strong couplings. We report a sharp transition, not a crossover, with a diverging effective mass at the critical coupling. We also look at a model with acoustic phonons, away from the perturbative limit, and again discover unusual polaron properties. Our work relies on the Bold Diagrammatic Monte Carlo (BDMC) method, which samples Feynman diagrammatic expansions efficiently, even those with weak sign problems. Proposed by Prokof'ev and Svistunov, it is extended to lattice polarons for the first time here. We also use the Momentum Average (MA) approximation, an analytical method proposed by Berciu, and find an excellent agreement with the BDMC results. A novel MA approximation able to treat dispersive phonons is also presented, along with a new exact solution for finite systems, inspired by the same formalism.