Abstract A general method for studying the properties of free and bound polarons at arbitrary strength of electron‐phonon coupling is given from a unified point of view by first performing a suitable canonical transformation of the electron and phonon field operators and subsequently applying the double‐time Green's function formalism to the transformed polaron‐phonon Hamiltonian. Calculating the one‐polaron Green's function in the lowest approximation and employing the variational principle for the free energy, a pair of coupled differential equations is derived for the effective wave function and a trial function as introduced by the canonical transformation. From these equations, the weakly excited polaron energy levels at low temperatures can be found in the entire range of coupling constants when external fields are present. In this paper, detailed calculations are confined to the optical‐polaron problem in the zero‐field case. Intermediate‐ and strong‐coupling formulas for temperature‐dependent polaron self‐energy and effective mass are given. The strong‐coupling regime appears to increase with temperature.
Abstract A general method for studying the properties of free and bound polarons at arbitrary strength of electron‐phonon coupling is given from a unified point of view by first performing a suitable canonical transformation of the electron and phonon field operators and subsequently applying the double‐time Green's function formalism to the transformed polaron‐phonon Hamiltonian. Calculating the one‐polaron Green's function in the lowest approximation and employing the variational principle for the free energy, a pair of coupled differential equations is derived for the effective wave function and a trial function as introduced by the canonical transformation. From these equations, the weakly excited polaron energy levels at low temperatures can be found in the entire range of coupling constants when external fields are present. In this paper, detailed calculations are confined to the optical‐polaron problem in the zero‐field case. Intermediate‐ and strong‐coupling formulas for temperature‐dependent polaron self‐energy and effective mass are given. The strong‐coupling regime appears to increase with temperature.
The ground-state characterization of the polaron problem is retrieved within the framework of a variational scheme proposed previously by Devreese et al for the bound polaron. The formulation is based on the standard canonical transformation of the strong coupling ansatz and consists of a variationally determined perturbative extension serving for the theory to interpolate in the overall range of the coupling constant. Specializing our considerations to the bulk and strict two-dimensional polaron models we see that the theory yields significantly improved energy upper bounds in the strong coupling regime and, moreover, extrapolates itself successfully towards the well-established weak coupling limits for all polaron quantities of general interest.
We have developed a new equation of motion Green's function theory for the large polaron, that is completely nonperturbative in origin, and therefore gives a unified treatment of the polaron ground state energy for all interaction strengths. The theory also incorporates translational invariance rigorously, a priori. We find the following new results. (a) Previous theories that are valid for all coupling strengths are based on a nonperturbative summation of infinite subsets of entire diagrams (equations of motion), rather than on recoil parts of all diagrams. We can thus recover lowest order results of Feynman for example, via a translationally noninvariant form of the Fock approximation, which resembles that reported in a previous work. (b) The total ground state energy we obtain in the strong coupling limit is E0 = −0.426α2 to lowest order in the coupling constant α, a result substantially below that of previous calculations, E0 = −0.106α2 (using gaussian trial functions in both cases). Our result is also below the lower bound of Lieb and Yamazaki, although we show in a following paper that our result is an upper bound. (c) The origin of the difference of our results compared to that of other strong coupling theories is several fold, and will be presented in detail in subsequent papers. However, the essence of the difference is that the phonons in the field, via their interaction with the electron, assume the role of independent degrees of freedom, in the limit of an infinite number of phonons. We then have a field theoretic problem with an infinite number of degrees of freedom which can exhibit a totally different physical behavior compared to one with a finite number of degrees of freedom. This aspect of the problem has not been given explicit recognition in all other theories, including that of Lieb and Yamazaki.
Ab initio electron propagator methods are developed to study electronic properties of molecular systems with strong electron-electron and electron-phonon interactions. For the calculation of electron Green's functions we apply a canonical small polaron transformation that intrinsically contains strong electron-phonon effects. In the transformed Hamiltonian, the energy levels for the noninteracting particles are shifted down by the relaxation (solvation) energies. The Coulomb integrals are also renormalized by the electron-phonon interaction. For certain values of the electron-phonon coupling constants, the renormalized Coulomb integrals can be negative which implies the attraction between two electrons. Within the small polaron transformation we develop a diagrammatic technique for the calculation of electron Green's function in which the electron-phonon interaction is already included into the multiple phonon correlation functions. Since the decoupling of the phonon correlation functions is impossible, and therefore, a Wick's theorem for such correlation functions is invalid, there is no Dyson equation for the electron Green's function. To find the electron Green's function, we use different approximations. One of them is a link-cluster approximation that includes diagonal transitions for the renormalized zeroth Green's function. In the linked-cluster approach the Dyson equation is derived in the most general case, where the self-energy operator is an arbitrary functional (not only in the Hartree-Fock approximation). It is shown that even a Hartree-Fock electron (hole) is not a particle any longer. It is a quasiparticle with a finite lifetime that depends on energy of particle and hole states in different ways. As a consequence of this, a standard description of a Hartree-Fock approximation in terms of wave functions becomes inappropriate in this problem. To challenge the linked-cluster approximation we develop a different approach: a sequential propagation approximation where scattering events occur only for sequential transitions. A self-consistent Hartree-Fock equation for a four-index Green's function matrix is derived. In conclusion, the proposed schemes can be considered for future method developments for quantum chemical calculations for large molecules with strong nonadiabatic effects, e-e correlated electron transfer reactions, and electron transport in molecular transport junctions.
It is well known that two qualitatively different pictures arise in the description of the motion of an electron in an ionic crystal. In the first case, when the electron interacts weakly with the ionic lattice of the crystal, its motion is the same as the motion of a free band electron whose energy is shifted downward with respect to the bottom of the conduction band, and the effective mass is replaced by the renormalized mass (weak-coupling polaron). In the other limiting case, when the electron interacts strongly with the ionic crystal, an entire series of different self-consistent states of the electron and lattice, each of which has its own effective mass and radius (strong-coupling polaron), arises. Although the motion of an electron in an ionic crystal is described in both limiting cases with the help of the Pekar-Frohlich Hamiltonian, each case requires its own special methods of analysis. Thus, for example, in the strong-coupling theory a special form of adiabatic perturbation theory has been developed in which the trans-lational degeneracy is removed before the perturbation expansion is made. It has been possible to relate both descriptions within the framework of Feynman's varia-tional approach, in which a smooth upper bound was obtained for the energy of the ground state of the polaron for all values of the electron-phonon coupling constant a (see Fig. 1). We call attention in this connection to Ref. 6, where the thermodynamic properties are calculated with the help of a special method for performing the calculations and for finding averages of T-products. The approach developed there is more convenient for performing specific calculations than is the use of path integrals. It seems to us that it has definite advantages in the construction of general proofs concerning the choice of the approximating Hamiltonians (this is illustrated in Ref. 6 for the example of the linear Bogolyubov model), in the investigation of the characteristics of the behavior of a polaron in external fields, etc.
New results for the free energy of a polaron in a strong magnetic field are presented. The method used here is the extended version of Feynman's pathintegral formalism in which the most general quadratic form containing an infinite number of variational parameters is employed to simulate the variational trial action. A set of nonlinear integral equations is established to determine the upper bound for the free energy in the entire range of electron-phonon coupling constants, temperatures and magnetic fields. Analytical solutions are provided for some limiting cases of interest, and are compared critically with existing calculations of the ground-state energy at zero temperature. The free energy thus obtained is exact to the lowest order in the coupling constant and provides the good interpolation up to the strong coupling regime.
A method previously derived for the approximate construction of the nucleon Green's function is here applied to that of the polaron. After making an arbitrary translation of the phonon variables, a linear integral equation for the Green's function is derived by means of a symmetrical treatment and a noncorrelation assumption, in complete analogy to the nucleon problem. This equation is solved through the introduction of a spectral representation in the special case of zero total momentum. The lowest energy state of the system is calculated in terms of the arbitrary translations and then minimized with respect to them. Using the simplest nontrivial cutoff procedure to obtain the variational equation and its solution, results are obtained for values of the coupling parameter $\ensuremath{\alpha}<~3$, and are compared with those of Feynman.