Polaron Effective Mass, Band Distortion, and Self-Trapping in the Holstein Molecular Crystal Model

We present polaron effective masses and selected polaron band structures of the Holstein molecular crystal model in 1-D as computed by the Global-Local variational method over a wide range of parameters. These results are augmented and supported by leading orders of both weak- and strong-coupling perturbation theory. The description of the polaron effective mass and polaron band distortion that emerges from this work is comprehensive, spanning weak, intermediate, and strong electron-phonon coupling, and non-adiabatic, weakly adiabatic, and strongly adiabatic regimes. Using the effective mass as the primary criterion, the self-trapping transition is precisely defined and located. Using related band-shape criteria at the Brillouin zone edge, the onset of band narrowing is also precisely defined and located. These two lines divide the polaron parameter space into three regimes of distinct polaron structure, essentially constituting a polaron phase diagram. Though the self-trapping transition is thusly shown to be a broad and smooth phenomenon at finite parameter values, consistency with notion of self-trapping as a critical phenomenon in the adiabatic limit is demonstrated. Generalizations to higher dimensions are considered, and resolutions of apparent conflicts with well-known expectations of adiabatic theory are suggested.Comment: 28 pages, 15 figure

Effects of dimensionality and anisotropy on the Holstein polaron

We apply weak-coupling perturbation theory and strong-coupling perturbation theory to the Holstein molecular crystal model in order to elucidate the effects of anisotropy on polaron properties in D dimensions. The ground state energy is considered as a primary criterion through which to study the effects of anisotropy on the self-trapping transition, the self-trapping line associated with this transition, and the adiabatic critical point. The effects of dimensionality and anisotropy on electron-phonon correlations and polaronic mass enhancement are studied, with particular attention given to the polaron radius and the characteristics of quasi-1D and quasi-2D structures. Perturbative results are confirmed by selected comparisons with variational calculations and quantum Monte Carlo data

Calculation of excited polaron states in the Holstein model

An exact diagonalization technique is used to investigate the low-lying excited polaron states in the Holstein model for the infinite one-dimensional lattice. For moderate values of the adiabatic ratio, a new and comprehensive picture, involving three excited (coherent) polaron bands below the phonon threshold, is obtained. The coherent contribution of the excited states to both the single-electron spectral density and the optical conductivity is evaluated and, due to the invariance of the Hamiltonian under the space inversion, the two are shown to contain complementary information about the single-electron system at zero temperature. The chosen method reveals the connection between the excited bands and the renormalized local phonon excitations of the adiabatic theory, as well as the regime of parameters for which the electron self-energy has notable non-local contributions. Finally, it is shown that the hybridization of two polaron states allows a simple description of the ground and first excited state in the crossover regime.Comment: 12 pages, 9 figures, submitted to PR

Polaron features of the one-dimensional Holstein Molecular Crystal Model

The polaron features of the one-dimensional Holstein Molecular Crystal Model are investigated by improving a variational method introduced recently and based on a linear superposition of Bloch states that describe large and small polaron wave functions. The mean number of phonons, the polaron kinetic energy, the electron-phonon local correlation function, and the ground state spectral weight are calculated and discussed. A crossover regime between large and small polaron for any value of the adiabatic parameter $\omega_0/t$ is found and a polaron phase diagram is proposed.Comment: 12 pages, 2 figure

Quantum Monte Carlo and variational approaches to the Holstein model

Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principal component representation of the phonon degrees of freedom allows to sample completely uncorrelated phonon configurations. The combination of these elements enables us to perform efficient simulations for a wide range of temperature, phonon frequency and electron-phonon coupling on clusters large enough to avoid finite-size effects. The algorithm is tested in one dimension and the data are compared with exact-diagonalization results and with existing work. Moreover, the ideas presented here can also be applied to the many-electron case. In the one-electron case considered here, the physics of the Holstein model can be described by a simple variational approach.Comment: 18 pages, 11 Figures, v2: one typo correcte

Polaron formation for a non-local electron-phonon coupling: A variational wave-function study

We introduce a variational wave-function to study the polaron formation when the electronic transfer integral depends on the relative displacement between nearest-neighbor sites giving rise to a non-local electron-phonon coupling with optical phonon modes. We analyze the ground state properties such as the energy, the electron-lattice correlation function, the phonon number and the spectral weight. Variational results are found in good agreement with analytic weak-coupling perturbative calculations and exact numerical diagonalization of small clusters. We determine the polaronic phase diagram and we find that the tendency towards strong localization is hindered from the pathological sign change of the effective next-nearest-neighbor hopping.Comment: 11 page

Mass Renormalization in the Su-Schrieffer-Heeger Model

This study of the one dimensional Su-Schrieffer-Heeger model in a weak coupling perturbative regime points out the effective mass behavior as a function of the adiabatic parameter $\omega_{\pi}/J$, $\omega_{\pi}$ is the zone boundary phonon energy and $J$ is the electron band hopping integral. Computation of low order diagrams shows that two phonons scattering processes become appreciable in the intermediate regime in which zone boundary phonons energetically compete with band electrons. Consistently, in the intermediate (and also moderately antiadiabatic) range the relevant mass renormalization signals the onset of a polaronic crossover whereas the electrons are essentially undressed in the fully adiabatic and antiadiabatic systems. The effective mass is roughly twice as much the bare band value in the intermediate regime while an abrupt increase (mainly related to the peculiar 1D dispersion relations) is obtained at $\omega_{\pi}\sim \sqrt{2}J$.Comment: To be published in Phys.Rev.B - 3 figure

Lattice dynamics effects on small polaron properties

This study details the conditions under which strong-coupling perturbation theory can be applied to the molecular crystal model, a fundamental theoretical tool for analysis of the polaron properties. I show that lattice dimensionality and intermolecular forces play a key role in imposing constraints on the applicability of the perturbative approach. The polaron effective mass has been computed in different regimes ranging from the fully antiadiabatic to the fully adiabatic. The polaron masses become essentially dimension independent for sufficiently strong intermolecular coupling strengths and converge to much lower values than those tradition-ally obtained in small-polaron theory. I find evidence for a self-trapping transition in a moderately adiabatic regime at an electron-phonon coupling value of .3. Our results point to a substantial independence of the self-trapping event on dimensionality.Comment: 8 pages, 5 figure

The Holstein Polaron

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 1d polaron energy at intermediate coupling). We compute ground and low-lying excited state properties of the model at continuous values of the wavevector $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 DMRG, Global Local and exact diagonalization. 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.Comment: 24 pages, 11 figures submitted to PR

Path integrals approach to resisitivity anomalies in anharmonic systems

Different classes of physical systems with sizeable electron-phonon coupling and lattice distortions present anomalous resistivity behaviors versus temperature. We study a molecular lattice Hamiltonian in which polaronic charge carriers interact with non linear potentials provided by local atomic fluctuations between two equilibrium sites. We study a molecular lattice Hamiltonian in which polaronic charge carriers interact with non linear potentials provided by local atomic fluctuations between two equilibrium sites. A path integral model is developed to select the class of atomic oscillations which mainly contributes to the partition function and the electrical resistivity is computed in a number of representative cases. We argue that the common origin of the observed resistivity anomalies lies in the time retarded nature of the polaronic interactions in the local structural instabilities.Comment: 4 figures, to appear in Phys.Rev.B, May 1st (2001