Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0510828 · Physical Review B

Polaron action for multimode dispersive phonon systems

P. E. Kornilovitch

Published Mar 27, 2006• 11 citations
arXivHTMLPDF

Abstract

The path-integral approach to the tight-binding polaron is extended to multiple optical phonon modes of arbitrary dispersion and polarization. The nonlinear lattice effects are neglected. Only one electron band is considered. The electron-phonon interaction is of the density-displacement type, but can be of arbitrary spatial range and shape. Feynman's analytical integration of ion trajectories is performed by transforming the electron-ion forces to the basis in which the phonon dynamical matrix is diagonal. The resulting polaron action is derived for the periodic and shifted boundary conditions in imaginary time. The former can be used for calculating polaron thermodynamics while the latter for the polaron mass and spectrum. The developed formalism is the analytical basis for numerical analysis of such models by path-integral Monte Carlo methods.

DOI
10.1103/physrevb.73.094305

Related papers

Path-integral approach to lattice polarons

P. E. Kornilovitch

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.

Journal of Physics Condensed Matter•May 30, 2007•6 citations

Extending the Feynman variational principle and analytical methods to lattice polarons

Serghei N. Klimin, Jacques Tempere, Matthew Houtput +6 more

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.

arXiv:2603.09609•Physical review. B./Physical review. B•May 27, 2026•0 citations

Polaron effective mass from Monte Carlo simulations

P. E. Kornilovitch, Edward R. Pike

A Monte Carlo algorithm for calculating polaron effective mass is proposed. It is based on the path-integral representation of a partial partition function with fixed total quasimomentum. Phonon degrees of freedom are integrated out analytically resulting in a single-electron system with retarded self-interaction and open boundary conditions in imaginary time. The effective mass is inversely proportional to the covariance of total energy calculated on an electron trajectory and the square distance between ends of the trajectory. The method has no limitations on values of model parameters and on the size and dimensionality of the system although large statistics are required for stable numerical results. The method is tested on the one-dimensional Holstein model for which simulation results are presented.

arXiv:cond-mat/9702226•Physical review. B, Condensed matter•Apr 1, 1997•53 citations

Lattice polaron formation: Effects of nonscreened electron-phonon interaction

Holger Fehske, Jan Loos, Gerhard Wellein

We explore the quasiparticle properties of lattice polarons on the basis of a quite general electron-phonon Hamiltonian with a long-range displacement-type interaction. To treat the dynamical quantum phonons without significant loss of accuracy we adapt an exact Lanczos diagonalization method and compute various static and dynamical quantities, such as the electron-lattice correlation function, the polaron band dispersion, the effective polaron mass, the kinetic energy, the single-particle spectral function, and the optical conductivity, on finite one-dimensional lattices for a wide range of model parameters. We compare the results with those obtained for the standard Holstein model with short-range electron-phonon interaction only.

arXiv:cond-mat/9911218•Physical review. B, Condensed matter•Mar 15, 2000•96 citations

Method for obtaining polaron mobility using real and imaginary time path-integral quantum Monte Carlo

Suzana Miladić‐Tešić, Nenad Vukmirović

We developed a path-integral quantum Monte Carlo-based methodology for calculation of polaron mobility in systems with electron-phonon interaction. Within the method, the current-current correlation function in both the imaginary and real time is calculated in a numerically exact way. The choice of basis for representation of the path integral enabled us to reduce the sign problem and perform real-time calculations for longer times. The DC polaron mobility was extracted by performing analytic continuation that makes use of both the real and imaginary-time data. The method was applied to the Holstein polaron model in one dimension. We obtained reliable results for the temperature dependence of the Holstein polaron mobility for interactions ranging from weak to strong and temperatures that are not too low.

Physical review. B./Physical review. B•May 31, 2023•10 citations

Ground-state dispersion and density of states from path-integral Monte Carlo: Application to the lattice polaron

P. E. Kornilovitch

A formula is derived that relates the ground-state dispersion of a many-body system with the end-to-end distribution of paths with open boundary conditions in imaginary time. The formula does not involve the energy estimator. It allows direct measurement of the ground-state dispersion by quantum Monte Carlo methods without analytical continuation or auxiliary fitting. The formula is applied to the lattice polaron problem. The exact polaron spectrum and density of states are calculated for several models in one, two, and three dimensions. In the adiabatic regime of the Holstein model, the polaron density of states deviates spectacularly from the free-particle shape.

arXiv:cond-mat/9903011•Physical review. B, Condensed matter•Aug 1, 1999•43 citations
The path-integral approach to the tight-binding polaron is extended to multiple optical phonon modes of arbitrary dispersion and polarization. The nonlinear lattice effects are neglected. Only one electron band is considered. The electron-phonon interaction is of the density-displacement type, but can be of arbitrary spatial range and shape. Feynman's analytical integration of ion trajectories is performed by transforming the electron-ion forces to the basis in which the phonon dynamical matrix is diagonal. The resulting polaron action is derived for the periodic and shifted boundary conditions in imaginary time. The former can be used for calculating polaron thermodynamics while the latter for the polaron mass and spectrum. The developed formalism is the analytical basis for numerical analysis of such models by path-integral Monte Carlo methods.
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.
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.
A Monte Carlo algorithm for calculating polaron effective mass is proposed. It is based on the path-integral representation of a partial partition function with fixed total quasimomentum. Phonon degrees of freedom are integrated out analytically resulting in a single-electron system with retarded self-interaction and open boundary conditions in imaginary time. The effective mass is inversely proportional to the covariance of total energy calculated on an electron trajectory and the square distance between ends of the trajectory. The method has no limitations on values of model parameters and on the size and dimensionality of the system although large statistics are required for stable numerical results. The method is tested on the one-dimensional Holstein model for which simulation results are presented.
We explore the quasiparticle properties of lattice polarons on the basis of a quite general electron-phonon Hamiltonian with a long-range displacement-type interaction. To treat the dynamical quantum phonons without significant loss of accuracy we adapt an exact Lanczos diagonalization method and compute various static and dynamical quantities, such as the electron-lattice correlation function, the polaron band dispersion, the effective polaron mass, the kinetic energy, the single-particle spectral function, and the optical conductivity, on finite one-dimensional lattices for a wide range of model parameters. We compare the results with those obtained for the standard Holstein model with short-range electron-phonon interaction only.
We developed a path-integral quantum Monte Carlo-based methodology for calculation of polaron mobility in systems with electron-phonon interaction. Within the method, the current-current correlation function in both the imaginary and real time is calculated in a numerically exact way. The choice of basis for representation of the path integral enabled us to reduce the sign problem and perform real-time calculations for longer times. The DC polaron mobility was extracted by performing analytic continuation that makes use of both the real and imaginary-time data. The method was applied to the Holstein polaron model in one dimension. We obtained reliable results for the temperature dependence of the Holstein polaron mobility for interactions ranging from weak to strong and temperatures that are not too low.
A formula is derived that relates the ground-state dispersion of a many-body system with the end-to-end distribution of paths with open boundary conditions in imaginary time. The formula does not involve the energy estimator. It allows direct measurement of the ground-state dispersion by quantum Monte Carlo methods without analytical continuation or auxiliary fitting. The formula is applied to the lattice polaron problem. The exact polaron spectrum and density of states are calculated for several models in one, two, and three dimensions. In the adiabatic regime of the Holstein model, the polaron density of states deviates spectacularly from the free-particle shape.