Polaron

The energy spectrum of an electron moving in a periodical potential of rigid crystal lattice is called the Bloch spectrum, which consists of allowed bands and forbidden bands. An electron with energy inside an allowed band moves as a free electron but has an effective mass that differs from the electron mass in vacuum. However, a crystal lattice is deformable and displacements of atoms (ions) from their equilibrium positions are described in terms of phonons. Electrons interact with these displacements, and this interaction is known as electron-phonon coupling. One possible scenario was proposed in the seminal 1933 paper by Lev Landau, which includes the production of a lattice defect such as an F-center and a trapping of the electron by this defect. A different scenario was proposed by Solomon Pekar that envisions dressing the electron with lattice polarization (a cloud of virtual polar phonons). Such an electron with the accompanying deformation moves freely across the crystal, but with increased effective mass.[3] Pekar coined for this charge carrier the term polaron.

Landau[4] and Pekar[5] formed the basis of polaron theory. A charge placed in a polarizable medium will be screened. Dielectric theory describes the phenomenon by the induction of a polarization around the charge carrier. The induced polarization will follow the charge carrier when it is moving through the medium. The carrier together with the induced polarization is considered as one entity, which is called a polaron (see Fig. 1).

While polaron theory was originally developed for electrons as dressed charges in a crystal field, there is no fundamental reason held against any other charged particle which might interact with phonons. Therefore, also other charged particles such as (electron) holes and ions should generally follow the polaron theory. For example, the proton polaron was identified experimentally in 2017[6] and on ceramic electrolytes after hypothesis[7] of its existence.

Fig. 1: Artist view of a polaron.[8] A conduction electron in an ionic crystal or a polar semiconductor repels the negative ions and attracts the positive ions. A self-induced potential arises, which acts back on the electron and modifies its physical properties.

Usually, in covalent semiconductors the couplings of electrons with lattice deformation is weak and the formation of polarons does not occur. In polar semiconductors the electrostatic interaction with induced polarization is strong and polarons are formed at low temperature, provided that the concentration of polarons is not large and the screening is not efficient. Another class of materials where polarons are observed is molecular crystals, where the interaction with molecular vibrations may be strong. In the case of polar semiconductors the interaction with polar phonons is described by the Fröhlich Hamiltonian. On the other hand, the interaction of electrons with molecular phonons is described by the Holstein Hamiltonian. Usually the models describing polarons may be divided to two classes. The first class represents continuum models where the discreteness of the crystal lattice is neglected. In that case polarons are weakly coupled or strongly coupled depending on whether the polaron binding energy is small or large compared to the phonon frequency. The second class of systems commonly considered are lattice models of polarons. In this case there may be small or large polarons, depending on the relation of the polaron radius and the lattice constant a.

A conduction electron in an ionic crystal or a polar semiconductor is the prototype of a polaron. Herbert Fröhlich proposed a model Hamiltonian for this polaron through which its dynamics are treated quantum mechanically (Fröhlich Hamiltonian).[10][11] The strength of electron phonon interaction is determined by the dimensionless coupling constant ${\displaystyle \alpha =(e^{2}/\kappa )(m/2\hbar ^{3}\omega )^{1/2}}$. Here ${\displaystyle m}$ is electron mass, ${\displaystyle \omega }$ is the phonon frequency and ${\displaystyle \kappa ^{-1}={\epsilon }_{\infty }^{-1}-{\epsilon }_{0}^{-1}}$, ${\displaystyle {\epsilon }_{0}}$, ${\displaystyle {\epsilon }_{\infty }}$ are static and high frequency dielectric constants. In table 1 the Fröhlich coupling constant is given for a few solids. The Fröhlich Hamiltonian for a single electron in a crystal using second quantization notation is:

${\displaystyle H=H_{\rm {e}}+H_{\rm {ph}}+H_{\rm {e-ph}}}$

${\displaystyle H_{\rm {e}}=\sum _{k,s}\xi (k,s)c_{k,s}^{\dagger }c_{k,s}}$

${\displaystyle H_{\rm {ph}}=\sum _{q,v}\omega _{q,v}a_{q,v}^{\dagger }a_{q,v}}$

${\displaystyle H_{\rm {e-ph}}={\frac {1}{\sqrt {2N}}}\sum _{k,s,q,v}\gamma (\alpha ,q,k,v)\omega _{qv}(c_{k,s}^{\dagger }c_{k-q,s}a_{q,v}+c_{k-q,s}^{\dagger }c_{k,s}a_{q,v}^{\dagger })}$

The exact form of γ depends on the material and the type of phonon being used in the model. In the case of a single polar mode ${\displaystyle \gamma (q)=i\hbar \omega ({\frac {4\pi \alpha }{V_{0}}}({\frac {\hbar }{m\omega }})^{1/2})^{1/2}{\frac {1}{q}}}$, here ${\displaystyle V_{0}}$ is the volume of the unit cell. In the case of molecular crystal γ is usually momentum independent constant. A detailed advanced discussion of the variations of the Fröhlich Hamiltonian can be found in J. T. Devreese and A. S. Alexandrov.[12] The terms Fröhlich polaron and large polaron are sometimes used synonymously, since the Fröhlich Hamiltonian includes the continuum approximation and long range forces. There is no known exact solution for the Fröhlich Hamiltonian with longitudinal optical (LO) phonons and linear ${\displaystyle \gamma }$ (the most commonly considered variant of the Fröhlich polaron) despite extensive investigations.[5][9][10][11][13][14][15][16][17][18]

Despite the lack of an exact solution, some approximations of the polaron properties are known.

The physical properties of a polaron differ from those of a band-carrier. A polaron is characterized by its self-energy ${\displaystyle \Delta E}$, an effective mass ${\displaystyle m^{*}}$ and by its characteristic response to external electric and magnetic fields (e. g. dc mobility and optical absorption coefficient).

When the coupling is weak (${\displaystyle \alpha }$ small), the self-energy of the polaron can be approximated as:[19]

${\displaystyle {\frac {\Delta E}{\hbar \omega }}\approx -\alpha -0.015919622\alpha ^{2},\qquad \qquad \qquad (1)\,}$

and the polaron mass ${\displaystyle m*}$, which can be measured by cyclotron resonance experiments, is larger than the band mass ${\displaystyle m}$ of the charge carrier without self-induced polarization:[20]

${\displaystyle {\frac {m^{*}}{m}}\approx 1+{\frac {\alpha }{6}}+0.0236\alpha ^{2}.\qquad \qquad \qquad (2)}$

When the coupling is strong (α large), a variational approach due to Landau and Pekar indicates that the self-energy is proportional to α² and the polaron mass scales as α⁴. The Landau–Pekar variational calculation[5] yields an upper bound to the polaron self-energy ${\displaystyle E<-C_{PL}\alpha ^{2}}$, valid for all α, where ${\displaystyle C_{PL}}$ is a constant determined by solving an integro-differential equation. It was an open question for many years whether this expression was asymptotically exact as α tends to infinity. Finally, Donsker and Varadhan,[21] applying large deviation theory to Feynman's path integral formulation for the self-energy, showed the large α exactitude of this Landau–Pekar formula. Later, Lieb and Thomas[22] gave a shorter proof using more conventional methods, and with explicit bounds on the lower order corrections to the Landau–Pekar formula.

Feynman[23] introduced the variational principle for path integrals to study the polaron. He simulated the interaction between the electron and the polarization modes by a harmonic interaction between a hypothetical particle and the electron. The analysis of an exactly solvable ("symmetrical") 1D-polaron model,[24][25] Monte Carlo schemes[26][27] and other numerical schemes[28] demonstrate the remarkable accuracy of Feynman's path-integral approach to the polaron ground-state energy. Experimentally more directly accessible properties of the polaron, such as its mobility and optical absorption, have been investigated subsequently.

In the strong coupling limit, ${\displaystyle \alpha \gg 1}$, the spectrum of excited states of a polaron begins with polaron-phonon bound states with energies less than ${\displaystyle \hbar \omega _{0}}$, where ${\displaystyle \omega _{0}}$ is the frequency of optical phonons.[29]

In the lattice models the main parameter is the polaron binding energy: ${\displaystyle E_{p}={\frac {1}{2N}}\sum _{q}|\gamma (q)|^{2}/\hbar \omega }$,[30] here summation is taken over the Brillouin zone. Note that this binding energy is purely adiabatic, i.e. does not depend on the ionic masses. For polar crystals the value of the polaron binding energy is strictly determined by the dielectric constants ${\displaystyle \epsilon _{0}}$,${\displaystyle \epsilon _{\infty }}$, and is of the order of 0.3-0.8 eV. If polaron binding energy ${\displaystyle E_{p}}$ is smaller than the hopping integral t the large polaron is formed for some type of electron-phonon interactions. In the case when ${\displaystyle E_{p}>t}$ the small polaron is formed. There are two limiting cases in the lattice polaron theory. In the physically important adiabatic limit ${\displaystyle t\gg \hbar \omega }$ all terms which involve ionic masses are cancelled and formation of polaron is described by nonlinear Schrödinger equation with nonadiabatic correction describing phonon frequency renormalization and polaron tunneling.[18][31] [32] In the opposite limit ${\displaystyle t\ll \hbar \omega }$ the theory represents the expansion in ${\displaystyle t/\hbar \omega }$.[18]

The expression for the magnetooptical absorption of a polaron is:[33]

${\displaystyle \Gamma (\Omega )\propto -{\frac {\operatorname {Im} \Sigma (\Omega )}{\left[\Omega -\omega _{\mathrm {c} }-\operatorname {Re} \Sigma (\Omega )\right]^{2}+\left[\operatorname {Im} \Sigma (\Omega )\right]^{2}}}.\qquad \qquad \qquad (3)}$

Here, ${\displaystyle \omega _{c}}$ is the cyclotron frequency for a rigid-band electron. The magnetooptical absorption Γ(Ω) at the frequency Ω takes the form Σ(Ω) is the so-called "memory function", which describes the dynamics of the polaron. Σ(Ω) depends also on α, β and ${\displaystyle \omega _{c}}$.

In the absence of an external magnetic field (${\displaystyle \omega _{c}=0}$) the optical absorption spectrum (3) of the polaron at weak coupling is determined by the absorption of radiation energy, which is reemitted in the form of LO phonons. At larger coupling, ${\displaystyle \alpha \geq 5.9}$, the polaron can undergo transitions toward a relatively stable internal excited state called the "relaxed excited state" (RES) (see Fig. 2). The RES peak in the spectrum also has a phonon sideband, which is related to a Franck–Condon-type transition.

Fig.2. Optical absorption of a polaron at ${\displaystyle \alpha =5}$ and 6. The RES peak is very intense compared with the Franck–Condon (FC) peak.[15][34]

A comparison of the DSG results [34] with the optical conductivity spectra given by approximation-free numerical [35] and approximate analytical approaches is given in ref.[36]

Calculations of the optical conductivity for the Fröhlich polaron performed within the Diagrammatic Quantum Monte Carlo method,[35] see Fig. 3, fully confirm the results of the path-integral variational approach[34] at ${\displaystyle \alpha \lesssim 3.}$ In the intermediate coupling regime ${\displaystyle 3<\alpha <6,}$ the low-energy behavior and the position of the maximum of the optical conductivity spectrum of ref.[35] follow well the prediction of Devreese.[34] There are the following qualitative differences between the two approaches in the intermediate and strong coupling regime: in ref.,[35] the dominant peak broadens and the second peak does not develop, giving instead rise to a flat shoulder in the optical conductivity spectrum at ${\displaystyle \alpha =6}$. This behavior can be attributed to the optical processes with participation of two[37] or more phonons. The nature of the excited states of a polaron needs further study.

Fig. 3: Optical conductivity spectra calculated within the Diagrammatic Quantum Monte Carlo method (open circles) compared to the DSG calculations (solid lines).[34][35]

The application of a sufficiently strong external magnetic field allows one to satisfy the resonance condition ${\displaystyle \Omega =\omega _{\mathrm {c} }+\operatorname {Re} \Sigma (\Omega )}$, which {(for ${\displaystyle \omega _{c}<\omega }$)} determines the polaron cyclotron resonance frequency. From this condition also the polaron cyclotron mass can be derived. Using the most accurate theoretical polaron models to evaluate ${\displaystyle \Sigma (\Omega )}$, the experimental cyclotron data can be well accounted for.

Evidence for the polaron character of charge carriers in AgBr and AgCl was obtained through high-precision cyclotron resonance experiments in external magnetic fields up to 16 T.[38] The all-coupling magneto-absorption calculated in ref.,[33] leads to the best quantitative agreement between theory and experiment for AgBr and AgCl. This quantitative interpretation of the cyclotron resonance experiment in AgBr and AgCl[38] by the theory of Peeters[33] provided one of the most convincing and clearest demonstrations of Fröhlich polaron features in solids.

Experimental data on the magnetopolaron effect, obtained using far-infrared photoconductivity techniques, have been applied to study the energy spectrum of shallow donors in polar semiconductor layers of CdTe.[39]

The polaron effect well above the LO phonon energy was studied through cyclotron resonance measurements, e. g., in II–VI semiconductors, observed in ultra-high magnetic fields.[40] The resonant polaron effect manifests itself when the cyclotron frequency approaches the LO phonon energy in sufficiently high magnetic fields.

In the lattice models the optical conductivity is given by the formula:[30]

${\displaystyle \sigma (\Omega )=n_{p}e^{2}a^{2}{\frac {\pi ^{\frac {1}{2}}t^{2}\left(1-e^{\frac {\hbar \Omega }{k_{\rm {B}}T}}\right)}{2\hbar ^{2}\Omega (E_{a}k_{\rm {B}}T)^{1/2}}}\exp \left(-{\frac {(\hbar \Omega -4E_{a})^{2}}{16E_{a}k_{\rm {B}}T}}\right)}$

Here ${\displaystyle E_{a}}$ is the activation energy of polaron, which is of the order of polaron binding energy ${\displaystyle E_{p}}$. This formula was derived and extensively discussed in[41][42][43] and was tested experimentally for example in photodoped parent compounds of high temperature superconductors.[44]

The great interest in the study of the two-dimensional electron gas (2DEG) has also resulted in many investigations on the properties of polarons in two dimensions.[45][46][47] A simple model for the 2D polaron system consists of an electron confined to a plane, interacting via the Fröhlich interaction with the LO phonons of a 3D surrounding medium. The self-energy and the mass of such a 2D polaron are no longer described by the expressions valid in 3D; for weak coupling they can be approximated as:[48][49]

${\displaystyle {\frac {\Delta E}{\hbar \omega }}\approx -{\frac {\pi }{2}}\alpha \ -0.06397\alpha ^{2};\qquad \qquad \qquad (4)\,}$

${\displaystyle {\frac {m^{*}}{m}}\approx 1+{\frac {\pi }{8}}\alpha \ +0.1272348\alpha ^{2}.\qquad \qquad \qquad (5)\,}$

It has been shown that simple scaling relations exist, connecting the physical properties of polarons in 2D with those in 3D. An example of such a scaling relation is:[47]

${\displaystyle {\frac {m_{\rm {2D}}^{*}(\alpha )}{m_{\rm {2D}}}}={\frac {m_{\rm {3D}}^{*}({\frac {3}{4}}\pi \alpha )}{m_{\rm {3D}}}},\qquad \qquad \qquad (6)\,}$

where ${\displaystyle m_{\mathrm {2D} }^{*}}$ (${\displaystyle m_{\mathrm {3D} }^{*}}$) and ${\displaystyle m_{\mathrm {2D} }}$ (${\displaystyle m_{\mathrm {3D} }}$) are, respectively, the polaron and the electron-band masses in 2D (3D).

The effect of the confinement of a Fröhlich polaron is to enhance the effective polaron coupling. However, many-particle effects tend to counterbalance this effect because of screening.[45][50]

Also in 2D systems cyclotron resonance is a convenient tool to study polaron effects. Although several other effects have to be taken into account (nonparabolicity of the electron bands, many-body effects, the nature of the confining potential, etc.), the polaron effect is clearly revealed in the cyclotron mass. An interesting 2D system consists of electrons on films of liquid He.[51][52] In this system the electrons couple to the ripplons of the liquid He, forming "ripplopolarons". The effective coupling can be relatively large and, for some values of the parameters, self-trapping can result. The acoustic nature of the ripplon dispersion at long wavelengths is a key aspect of the trapping.

For GaAs/Al_xGa_1−xAs quantum wells and superlattices, the polaron effect is found to decrease the energy of the shallow donor states at low magnetic fields and leads to a resonant splitting of the energies at high magnetic fields. The energy spectra of such polaronic systems as shallow donors ("bound polarons"), e. g., the D₀ and D⁻ centres, constitute the most complete and detailed polaron spectroscopy realised in the literature.[53]

In GaAs/AlAs quantum wells with sufficiently high electron density, anticrossing of the cyclotron-resonance spectra has been observed near the GaAs transverse optical (TO) phonon frequency rather than near the GaAs LO-phonon frequency.[54] This anticrossing near the TO-phonon frequency was explained in the framework of the polaron theory.[55]

Besides optical properties,[9][17][56] many other physical properties of polarons have been studied, including the possibility of self-trapping, polaron transport,[57][58] magnetophonon resonance, etc.

Significant are also the extensions of the polaron concept: acoustic polaron, piezoelectric polaron, electronic polaron, bound polaron, trapped polaron, spin polaron, molecular polaron, solvated polarons, polaronic exciton, Jahn-Teller polaron, small polaron, bipolarons and many-polaron systems.[9] These extensions of the concept are invoked, e. g., to study the properties of conjugated polymers, colossal magnetoresistance perovskites, high-${\displaystyle T_{c}}$ superconductors, layered MgB₂ superconductors, fullerenes, quasi-1D conductors, semiconductor nanostructures.

The possibility that polarons and bipolarons play a role in high-${\displaystyle T_{c}}$ superconductors has renewed interest in the physical properties of many-polaron systems and, in particular, in their optical properties. Theoretical treatments have been extended from one-polaron to many-polaron systems.[9][59][60]

A new aspect of the polaron concept has been investigated for semiconductor nanostructures: the exciton-phonon states are not factorizable into an adiabatic product Ansatz, so that a non-adiabatic treatment is needed.[61] The non-adiabaticity of the exciton-phonon systems leads to a strong enhancement of the phonon-assisted transition probabilities (as compared to those treated adiabatically) and to multiphonon optical spectra that are considerably different from the Franck–Condon progression even for small values of the electron-phonon coupling constant as is the case for typical semiconductor nanostructures.[61]

In biophysics Davydov soliton is a propagating along the protein α-helix self-trapped amide I excitation that is a solution of the Davydov Hamiltonian. The mathematical techniques that are used to analyze Davydov's soliton are similar to some that have been developed in polaron theory. In this context the Davydov soliton corresponds to a polaron that is (i) large so the continuum limit approximation in justified, (ii) acoustic because the self-localization arises from interactions with acoustic modes of the lattice, and (iii) weakly coupled because the anharmonic energy is small compared with the phonon bandwidth.[62]

It has been shown that the system of an impurity in a Bose–Einstein condensate is also a member of the polaron family.[63] This allows the hitherto inaccessible strong coupling regime to be studied, since the interaction strengths can be externally tuned through the use of a Feshbach resonance. This was recently realized experimentally by two research groups.[64][65] The existence of the polaron in a Bose–Einstein condensate was demonstrated for both attractive and repulsive interactions, including the strong coupling regime and dynamically observed.[66]

• Exciton
• Sigurd Zienau
• TI-polaron