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

Journal of Physics Condensed Matter

Trotter-number scaling in Monte Carlo studies of the small polaron

P. E. Kornilovitch

Published Dec 1, 1997• 9 citations
Publisher

Abstract

A path-integral quantum Monte Carlo method is applied to the two-dimensional single-polaron Holstein model. Simulation data are shown to exhibit -scaling with the number of imaginary-time slices M used in the Trotter decomposition. Numerical extrapolation to yields very accurate estimation of polaron energetic characteristics. Kinetic, potential and total energy are calculated as functions of the electron - phonon coupling strength for different phonon frequencies. The small-polaron regime is found in the strong-coupling limit for all of the frequencies studied. The transition to the polaron state is sharp in the adiabatic regime but broadens as the frequency increases. In the adiabatic limit the polaron forms at .

DOI
10.1088/0953-8984/9/48/011

Related papers

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

Quantum Monte Carlo approach to the Holstein polaron problem

Martin Hohenadler, Hans Gerd Evertz, Wolfgang von der Linden

Abstract A recently developed quantum Monte Carlo approach to the Holstein model with one electron [Phys. Rev. B 69, 024301 (2004)] is extended in several aspects. In addition to a straight‐forward generalization to higher dimensions, a checkerboard breakup for the hopping part of the Hamiltonian is introduced, which allows significantly faster simulations. Moreover, results are extrapolated towards the limit of continuous imaginary time to remove the error due to the Trotter–Suzuki approximation. The performance and accuracy of the current approach are compared to existing quantum Monte Carlo methods. To demonstrate the applicability of our method for almost arbitrary parameters, we study the small‐polaron crossover in one to three dimensions by calculating the electronic kinetic energy. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

arXiv:cond-mat/0401634•physica status solidi (b)•Apr 7, 2005•8 citations

Effect of electron-phonon interaction range on lattice polaron dynamics: A continuous-time quantum Monte Carlo study

Paul Spencer, John H. Samson, P. E. Kornilovitch +1 more

We present the numerically exact ground-state energy, effective mass, and isotope exponents of a one-dimensional lattice polaron, valid for any range of electron-phonon interaction, applying a continuous-time quantum Monte Carlo (QMC) technique in a wide range of coupling strength and adiabatic ratio. The QMC method is free from any systematic finite-size and finite-time-step errors. We compare our numerically exact results with analytical weak-coupling theory and with the strong-coupling $1∕\ensuremath{\lambda}$ expansion. We show that the exact results agree well with the canonical Fr\"ohlich and Holstein-Lang-Firsov theories in the weak and strong coupling limits, respectively, for any range of interaction. We find a strong dependence of the polaron dynamics on the range of interaction. An increased range of interaction has a similar effect to an increased (less adiabatic) phonon frequency: specifically, a reduction in the effective mass.

arXiv:cond-mat/0407250•Physical Review B•May 31, 2005•49 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

Photoemission spectra of many-polaron systems

Martin Hohenadler, Danilo R. Neuber, Wolfgang von der Linden +3 more

The crossover from low to high carrier densities in a many-polaron system is studied in the framework of the one-dimensional spinless Holstein model, using unbiased numerical methods. Combining a novel quantum Monte Carlo approach and exact diagonalization, accurate results for the single-particle spectrum and the electronic kinetic energy on fairly large systems are obtained. A detailed investigation of the quality of the Monte Carlo data is presented. In the physically most important adiabatic intermediate electron-phonon coupling regime, for which no analytical results are available, we observe a dissociation of polarons with increasing band filling, leading to normal metallic behavior, while for parameters favoring small polarons, no such density-driven changes occur. The present work points toward the inadequacy of single-polaron theories for a number of polaronic materials such as the manganites.

arXiv:cond-mat/0412010•Physical Review B•Jun 23, 2005•48 citations

Phonon spectral function of the Holstein polaron

J Loos, M Hohenadler, A Alvermann +1 more

The phonon spectral function of the one-dimensional Holstein model is obtained within weak and strong-coupling approximations based on analytical self-energy calculations. The characteristic excitations found in the limit of small charge-carrier density are related to the known (electronic) spectral properties of Holstein polarons such as the polaron band dispersion. Particular emphasis is laid on the different physics occurring in the adiabatic and anti-adiabatic regimes, respectively. Comparison is made with a cluster approach exploiting exact numerical results on small systems to yield an approximation for the thermodynamic limit. This method, similar to cluster perturbation theory, confirms the analytical findings, and yields accurate results also in the intermediate-coupling regime.

arXiv:cond-mat/0604040•Journal of Physics Condensed Matter•Jul 21, 2006•20 citations
A path-integral quantum Monte Carlo method is applied to the two-dimensional single-polaron Holstein model. Simulation data are shown to exhibit -scaling with the number of imaginary-time slices M used in the Trotter decomposition. Numerical extrapolation to yields very accurate estimation of polaron energetic characteristics. Kinetic, potential and total energy are calculated as functions of the electron - phonon coupling strength for different phonon frequencies. The small-polaron regime is found in the strong-coupling limit for all of the frequencies studied. The transition to the polaron state is sharp in the adiabatic regime but broadens as the frequency increases. In the adiabatic limit the polaron forms at .
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.
Abstract A recently developed quantum Monte Carlo approach to the Holstein model with one electron [Phys. Rev. B 69, 024301 (2004)] is extended in several aspects. In addition to a straight‐forward generalization to higher dimensions, a checkerboard breakup for the hopping part of the Hamiltonian is introduced, which allows significantly faster simulations. Moreover, results are extrapolated towards the limit of continuous imaginary time to remove the error due to the Trotter–Suzuki approximation. The performance and accuracy of the current approach are compared to existing quantum Monte Carlo methods. To demonstrate the applicability of our method for almost arbitrary parameters, we study the small‐polaron crossover in one to three dimensions by calculating the electronic kinetic energy. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
We present the numerically exact ground-state energy, effective mass, and isotope exponents of a one-dimensional lattice polaron, valid for any range of electron-phonon interaction, applying a continuous-time quantum Monte Carlo (QMC) technique in a wide range of coupling strength and adiabatic ratio. The QMC method is free from any systematic finite-size and finite-time-step errors. We compare our numerically exact results with analytical weak-coupling theory and with the strong-coupling $1∕\ensuremath{\lambda}$ expansion. We show that the exact results agree well with the canonical Fr\"ohlich and Holstein-Lang-Firsov theories in the weak and strong coupling limits, respectively, for any range of interaction. We find a strong dependence of the polaron dynamics on the range of interaction. An increased range of interaction has a similar effect to an increased (less adiabatic) phonon frequency: specifically, a reduction in the effective mass.
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.
The crossover from low to high carrier densities in a many-polaron system is studied in the framework of the one-dimensional spinless Holstein model, using unbiased numerical methods. Combining a novel quantum Monte Carlo approach and exact diagonalization, accurate results for the single-particle spectrum and the electronic kinetic energy on fairly large systems are obtained. A detailed investigation of the quality of the Monte Carlo data is presented. In the physically most important adiabatic intermediate electron-phonon coupling regime, for which no analytical results are available, we observe a dissociation of polarons with increasing band filling, leading to normal metallic behavior, while for parameters favoring small polarons, no such density-driven changes occur. The present work points toward the inadequacy of single-polaron theories for a number of polaronic materials such as the manganites.
The phonon spectral function of the one-dimensional Holstein model is obtained within weak and strong-coupling approximations based on analytical self-energy calculations. The characteristic excitations found in the limit of small charge-carrier density are related to the known (electronic) spectral properties of Holstein polarons such as the polaron band dispersion. Particular emphasis is laid on the different physics occurring in the adiabatic and anti-adiabatic regimes, respectively. Comparison is made with a cluster approach exploiting exact numerical results on small systems to yield an approximation for the thermodynamic limit. This method, similar to cluster perturbation theory, confirms the analytical findings, and yields accurate results also in the intermediate-coupling regime.