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

The Journal of Chemical Physics

Converging toward a practical solution of the Holstein molecular crystal model

Aldo Humberto Romero, David W. Brown, Katja Lindenberg

Published Oct 22, 1998• 78 citations
Publisher

Abstract

We present selected results for the Holstein molecular crystal model in one space dimension as determined by the Global–Local variational method, including complete polaron energy bands, ground state energies, and effective masses. We juxtapose our results with specific comparable results of numerous other methodologies of current interest, including quantum Monte Carlo, cluster diagonalization, dynamical mean field theory, density matrix renormalization group, semiclassical analysis, weak-coupling perturbation theory, and strong-coupling perturbation theory. Taken as a whole, these methodologies are mutually confirming and provide a comprehensive and quantitatively accurate description of polaron properties in essentially any regime. In particular, this comparison confirms the Global–Local variational method as being highly accurate over a wide range of the polaron parameter space, from the nonadiabatic limit to the extremes of high adiabaticity, from weak coupling through intermediate coupling to strong coupling.

DOI
10.1063/1.477305

Related papers

Solving the Holstein molecular crystal model

Aldo Humberto Romero, David W. Brown, Katja Lindenberg

We present selected results for the Holstein molecular crystal model in one space dimension as determined by the Global-Local variational method, including complete polaron energy bands, ground state energies, and effective masses. We juxtapose our results with specific comparable results of numerous other methodologies of current interest, including quantum Monte Carlo, cluster diagonalization, dynamical mean field theory, density matrix renormalization group, semiclassical analysis, weak-coupling perturbation theory, and strong-coupling perturbation theory. This comparison confirms the Global-Local variational method as being highly accurate over a wide range of the polaron parameter space, from the non-adiabatic limit to the extremes of high adiabaticity, from weak coupling through intermediate coupling to strong coupling.

arXiv:cond-mat/9710321•Oct 29, 1997•0 citations

Holstein polaron

J. Bonča, S. A. Trugman, I. Batistić

We describe a variational method to solve the Holstein model for an electron coupled to dynamical, quantum phonons on an infinite lattice. The variational space can be systematically expanded to achieve high accuracy with modest computational resources (12-digit accuracy for the one-dimensional polaron energy at intermediate coupling). We compute ground-state and low-lying excited-state properties of the model at continuous values of the wave vector k in essentially all parameter regimes. Our results for the polaron energy band, effective mass, and correlation functions compare favorably with those of other numerical techniques, including the density-matrix renormalization-group technique, the global-local method, and the exact diagonalization technique. We find a phase transition for the first excited state between a bound and unbound system of a polaron and an additional phonon excitation. The phase transition is also treated in strong-coupling perturbation theory.

arXiv:cond-mat/9812252•Physical review. B, Condensed matter•Jul 15, 1999•281 citations

Dimensionality effects on the Holstein polaron

Li-Chung Ku, S. A. Trugman, J. Bonča

Based on a recently developed variational method, we explore the properties of the Holstein polaron on an infinite lattice in D dimensions, where $1<~D<~4.$ The computational method converges as a power law, so that highly accurate results can be achieved with modest resources. We present the most accurate ground state energy (with no small parameter) to date for polaron problems, 21 digits for the one-dimensional (1D) polaron at intermediate coupling. The dimensionality effects on polaron band dispersion, effective mass, and electron-phonon (el-ph) correlation functions are investigated in all coupling regimes. It is found that the crossover to large effective mass of the higher-dimensional polaron is much sharper than the 1D polaron. The correlation length between the electron and phonons decreases significantly as the dimension increases. Our results compare favorably with those of the quantum Monte Carlo, dynamical mean-field theory, density-matrix renormalization-group, and Toyozawa variational methods. We demonstrate that the Toyozawa wave function is qualitatively correct for the ground-state energy and the two-point electron-phonon correlation functions, but fails for the three-point functions. Based on this finding, we propose an improved Toyozawa variational wave function.

arXiv:cond-mat/0109282•Physical review. B, Condensed matter•Apr 29, 2002•127 citations

Self-trapping line of the Holstein molecular crystal model in one dimension

Aldo Humberto Romero, David W. Brown, Katja Lindenberg

The ground state of the Holstein molecular crystal model of a single electron in one dimension is studied using the Global-Local variational method, analyzing in particular the total energy, kinetic energy, phonon energy, and interaction energy over a broad region of the polaron parameter space. Through the application of objective criteria, a curve is identified that simply, accurately, and robustly locates the self-trapping transition separating small polaron and large polaron behavior. Particular attention is given to the kinetic energy, which is compared quantitatively with other perturbative and nonperturbative methods.

arXiv:cond-mat/9806031•Physical review. B, Condensed matter•Aug 15, 1999•37 citations

Variational approach for the Holstein molecular-crystal model

V. Cataudella, G. De Filippis, Giuseppe Iadonisi

A variational technique is developed to investigate the polaronic features of the Holstein molecular-crystal model. It is based on a linear superposition of Bloch states that describe large and small polaron wave functions. It is shown that this method provides a very good description of the regime characterized by intermediate values of the electron-phonon coupling constant (the so-called intermediate polaron) for any value of the adiabatic parameter ${\ensuremath{\omega}}_{0}/t.$ The polaron ground state energy in one and two dimensions is calculated and successfully compared with the best estimates available providing a clear physical interpretation of the intermediate polaron. The band structure, the spectral weight of the ground state and the lattice displacement associated to the polaron are also calculated and discussed. The method has the advantage of requiring very little computational effort.

arXiv:cond-mat/9909305•Physical review. B, Condensed matter•Dec 1, 1999•42 citations

Self-trapping problem of electrons or excitons in one dimension

Gerhard Wellein, Holger Fehske

We present a detailed numerical study of the one-dimensional Holstein model with a view to understanding the self-trapping process of electrons or excitons in crystals with short-range particle-lattice interactions. Applying a very efficient variational Lanczos method, we are able to analyze the ground-state properties of the system in the weak- and strong-coupling, adiabatic and nonadiabatic regimes on lattices large enough to eliminate finite-size effects. In particular, we obtain the complete phase diagram and comment on the existence of a critical length for self-trapping in finite (closed) one-dimensional systems. In order to characterize large and small polaron states we calculate self-consistently the lattice distortions and the particle-phonon correlation functions. In the strong-coupling case, two distinct types of small polaron states are shown to be possible according to the relative importance of static displacement field and dynamic polaron effects. Special emphasis is on the intermediate-coupling regime, which we also study by means of direct diagonalization, preserving the full dynamics and quantum nature of phonons. The crossover from large to small polarons shows up in a strong decrease of the kinetic energy accompanied by a substantial change in the optical absorption spectra. We show that our numerical results in all important limiting cases reveal excellent agreement with both analytical perturbation theory predictions and very recent density matrix renormalization group data.

arXiv:cond-mat/9712238•Physical review. B, Condensed matter•Sep 1, 1998•101 citations
We present selected results for the Holstein molecular crystal model in one space dimension as determined by the Global–Local variational method, including complete polaron energy bands, ground state energies, and effective masses. We juxtapose our results with specific comparable results of numerous other methodologies of current interest, including quantum Monte Carlo, cluster diagonalization, dynamical mean field theory, density matrix renormalization group, semiclassical analysis, weak-coupling perturbation theory, and strong-coupling perturbation theory. Taken as a whole, these methodologies are mutually confirming and provide a comprehensive and quantitatively accurate description of polaron properties in essentially any regime. In particular, this comparison confirms the Global–Local variational method as being highly accurate over a wide range of the polaron parameter space, from the nonadiabatic limit to the extremes of high adiabaticity, from weak coupling through intermediate coupling to strong coupling.
We present selected results for the Holstein molecular crystal model in one space dimension as determined by the Global-Local variational method, including complete polaron energy bands, ground state energies, and effective masses. We juxtapose our results with specific comparable results of numerous other methodologies of current interest, including quantum Monte Carlo, cluster diagonalization, dynamical mean field theory, density matrix renormalization group, semiclassical analysis, weak-coupling perturbation theory, and strong-coupling perturbation theory. This comparison confirms the Global-Local variational method as being highly accurate over a wide range of the polaron parameter space, from the non-adiabatic limit to the extremes of high adiabaticity, from weak coupling through intermediate coupling to strong coupling.
We describe a variational method to solve the Holstein model for an electron coupled to dynamical, quantum phonons on an infinite lattice. The variational space can be systematically expanded to achieve high accuracy with modest computational resources (12-digit accuracy for the one-dimensional polaron energy at intermediate coupling). We compute ground-state and low-lying excited-state properties of the model at continuous values of the wave vector k in essentially all parameter regimes. Our results for the polaron energy band, effective mass, and correlation functions compare favorably with those of other numerical techniques, including the density-matrix renormalization-group technique, the global-local method, and the exact diagonalization technique. We find a phase transition for the first excited state between a bound and unbound system of a polaron and an additional phonon excitation. The phase transition is also treated in strong-coupling perturbation theory.
Based on a recently developed variational method, we explore the properties of the Holstein polaron on an infinite lattice in D dimensions, where 1< D< 4.1<~D<~4.1< D< 4. The computational method converges as a power law, so that highly accurate results can be achieved with modest resources. We present the most accurate ground state energy (with no small parameter) to date for polaron problems, 21 digits for the one-dimensional (1D) polaron at intermediate coupling. The dimensionality effects on polaron band dispersion, effective mass, and electron-phonon (el-ph) correlation functions are investigated in all coupling regimes. It is found that the crossover to large effective mass of the higher-dimensional polaron is much sharper than the 1D polaron. The correlation length between the electron and phonons decreases significantly as the dimension increases. Our results compare favorably with those of the quantum Monte Carlo, dynamical mean-field theory, density-matrix renormalization-group, and Toyozawa variational methods. We demonstrate that the Toyozawa wave function is qualitatively correct for the ground-state energy and the two-point electron-phonon correlation functions, but fails for the three-point functions. Based on this finding, we propose an improved Toyozawa variational wave function.
The ground state of the Holstein molecular crystal model of a single electron in one dimension is studied using the Global-Local variational method, analyzing in particular the total energy, kinetic energy, phonon energy, and interaction energy over a broad region of the polaron parameter space. Through the application of objective criteria, a curve is identified that simply, accurately, and robustly locates the self-trapping transition separating small polaron and large polaron behavior. Particular attention is given to the kinetic energy, which is compared quantitatively with other perturbative and nonperturbative methods.
A variational technique is developed to investigate the polaronic features of the Holstein molecular-crystal model. It is based on a linear superposition of Bloch states that describe large and small polaron wave functions. It is shown that this method provides a very good description of the regime characterized by intermediate values of the electron-phonon coupling constant (the so-called intermediate polaron) for any value of the adiabatic parameter ${\ensuremath{\omega}}_{0}/t.$ The polaron ground state energy in one and two dimensions is calculated and successfully compared with the best estimates available providing a clear physical interpretation of the intermediate polaron. The band structure, the spectral weight of the ground state and the lattice displacement associated to the polaron are also calculated and discussed. The method has the advantage of requiring very little computational effort.
We present a detailed numerical study of the one-dimensional Holstein model with a view to understanding the self-trapping process of electrons or excitons in crystals with short-range particle-lattice interactions. Applying a very efficient variational Lanczos method, we are able to analyze the ground-state properties of the system in the weak- and strong-coupling, adiabatic and nonadiabatic regimes on lattices large enough to eliminate finite-size effects. In particular, we obtain the complete phase diagram and comment on the existence of a critical length for self-trapping in finite (closed) one-dimensional systems. In order to characterize large and small polaron states we calculate self-consistently the lattice distortions and the particle-phonon correlation functions. In the strong-coupling case, two distinct types of small polaron states are shown to be possible according to the relative importance of static displacement field and dynamic polaron effects. Special emphasis is on the intermediate-coupling regime, which we also study by means of direct diagonalization, preserving the full dynamics and quantum nature of phonons. The crossover from large to small polarons shows up in a strong decrease of the kinetic energy accompanied by a substantial change in the optical absorption spectra. We show that our numerical results in all important limiting cases reveal excellent agreement with both analytical perturbation theory predictions and very recent density matrix renormalization group data.