Electron relaxation in metals: Theory and exact analytical solutions
Abstract
The non-equilibrium dynamics of electrons is of a great experimental and theoretical value providing important microscopic parameters of the Coulomb and electron-phonon interactions in metals and other cold plasmas. Because of the mathematical complexity of collision integrals theories of electron relaxation often rely on the assumption that electrons are in a ”quasi-equilibrium” (QE) with a time-dependent temperature, or on the numerical integration of the time-dependent Boltzmann equation. We transform the integral Boltzmann equation to a partial differential Schrödinger-like equation with imaginary time in a one-dimensional ”coordinate” space reciprocal to energy which allows for exact analytical solutions in both cases of electron-electron and electron-phonon relaxation. The exact relaxation rates are compared with the QE relaxation rates at high and low temperatures.
pacs:
71.38.-k, 74.40.+k, 72.15.Jf, 74.72.-h, 74.25.FyI Introduction
In recent years investigations of photo-response functions in advanced materials have gone through a vigorous revival. In particular, laser ”pump-probe” techniques, where a second probe pulse is delayed in time with respect to the pump pulse, provide unique information on the strength of electron-electron (e-e) and electron-phonon (e-ph) interactions in metals and doped insulators if an adequate theory is in place.
At present detailed experimental data on relaxation processes is collected for metals (1); (2); (3); (4); (5) and high-temperature superconductors (6); (7); (8); (9); (10); (11); (12). The pump-probe experiments are routinely analyzed in the framework of the so-called two temperature model (TTM) (13); (14). The model is based on the assumption that electrons and phonons are in a thermal quasi-equilibrium (QE) with two different time-dependent temperatures and , respectively. The comprehensive analysis of experimental data collected at room temperature (1) allowed for a determination of the electron-phonon (e-ph) coupling constant of many metals and low-temperature superconductors in the framework of TTM.
Similar experiments and their analysis were performed on high-temperature superconductors. The femtosecond time-resolved measurements on the high-T superconductors TlBaCaCuO (6) and YBaCuO (7); (8) found a relaxation process below T, which is distinct from the equilibration of hot carriers in the normal state. A relatively strong e-ph coupling, (8), a rapid decrease of the photo-response decay rate with decreasing temperature (16) were found in YBaCuO, and the phonon bottleneck (15); (12); (18); (17) or a biparticle recombination (16); (19) were observed below T. More recently a time-resolved photoemission spectroscopy (20) and the standard pump-probe optical measurements (21) have been performed on BiSrCaCuO. Their TTM analysis has led to a rather weak e-ph coupling, .
The pump-probe techniques have a potential to resolve a controversial issue on weather the e-ph interaction is crucial (22) or weak and inessential (23) for the mechanism of high-temperature superconductivity. The pioneering work by Kaganov, Lifshitz and Tanatarov (13) and subsequent TTM studies are based on the assumption that electrons are in the thermal QE-state because the e-e relaxation time is supposed to be much shorter than the e-ph relaxation time. This assumption is of course incorrect on a femtosecond scale comparable with the e-e scattering time of highly excited electrons, but the expectation has been that deviation from QE may not in fact have much influence on the electron energy relaxation on a larger time scale (14) (for discussions of TTM with respect to some experiments see for example (24); (25)).
However later on it has been realized that nonthermal effects are essential even on a picosecond scale, comparable with the e-ph relaxation time, when conditions of low laser excitation power and relatively low temperature are chosen (26). Under these conditions, the e-e collision rate becomes strongly suppressed as a result of the Pauli exclusion principle. Numerically integrating the Boltzmann equation with e-e and e-ph collision integrals Groeneveld, Sprik, and Lagendijk (26) have shown that the electron gas cannot attain a thermal distribution by e-e collisions on the time scale of the e-ph energy relaxation. A departure from QE leads to an increase of the e-ph energy relaxation time with respect to the QE expectation. As a consequence of this departure one might underestimate the e-ph coupling using TTM.
While numerical integrations of the Boltzmann equation can describe the time evolution of the electron distribution function on any time scale, they require a number of input parameters, which might be unknown a priori. Here an analytical approach to this long-standing problem is developed. We reduce the integral Boltzmann equation to a differential Schrödinger-like equation using an auxiliary space reciprocal to energy and find exact analytical time-dependent distributions of electrons in both cases of electron-electron and electron-phonon relaxation. We also derive long-time relaxation rates of response functions and compare them with TTM.
Ii Electron-electron relaxation
Let us first consider a nonthermal relaxation of the electron distribution function caused by electron-electron collisions, which is described by the following Boltzmann equation,
(1) | |||||
(2) |
Here , is the matrix element of the electron-electron scattering (pseudo)potential, is the electron energy with respect to the equilibrium chemical potential. For transparency we drop the time argument in the distribution function. If the distribution function depends only on energy and time, , one can average this equation over the angles of as
(3) |
where is the density of states (DOS) per spin, with the following result
(4) | |||
(5) |
where
(7) | |||||
We restrict our theory to relaxations involving non-equilibrium electron-hole excitations with energies much less than the equilibrium Fermi energy, . Since the kernel has variation on a scale of the Fermi energy, one can approximate it by a constant, . This constant is related to the Coulomb pseudo-potential , important in the theory of superconductivity, . Assuming a low laser excitation power we linearize Eq.(5) by introducing a small non-equilibrium correction, , to the equilibrium distribution, ,
(8) |
where . Keeping terms linear in and measuring energies in units of , which is the only relevant energy scale of the problem, one obtains
(10) | |||||
(11) |
Performing simple integrations in linearized Eq.(11) yields
(12) | |||
(13) |
where
(14) |
is the familiar lifetime of electron-hole excitations in the Fermi liquid. Here we have used the integral with .
The second term on the right-hand side of Eq.(13) describes a source of quasi-particles due to inelastic electron-electron collisions. Collisions in cold degenerate plasmas differ essentially from quasi-elastic collisions in classical (hot) plasmas. In the latter the energy transfer is small compared with the electron energy due to a long-range character of the Coulomb potential, so that one can approximate the Boltzmann collision integral by the differential Landau-Fokker-Plank (LFP) equation (see, for example Ref.(27)). As one can see from Eq.(13) it is not the case in metals. The collision energy transfer in metals is about the same as the excitation energy itself, which makes the differential LFP approximation unacceptable here.
Remarkably the electron-electron collision integral acquires a differential form in a reciprocal auxiliary-time space introduced via the Fourier transform of Eq.(13), rather than in the energy space as in the LFP case. Let us consider non-equilibrium states conserving the electron-hole symmetry, so that the non-equilibrium part of the distribution is an odd function of energy, . If one determines a function,
(15) |
then the Boltzmann equation is simplified as
(16) |
We shall see below that a ”bound state” of the effective ”Schrödinger” equation for the Fourier transform of corresponds to the stationary quasi-equilibrium distribution. Taking the Fourier transform of Eq.(16), we arrive at an exact differential counterpart of the Boltzmann equation,
(17) |
where
(18) |
and
(19) |
Here another integral has been used.
The solution of Eq.(17) is found as a superposition of normalized eigenstates, , of a textbook Hamiltonian (28),
(20) |
where coefficients are determined by the initial non-equilibrium distribution function at ,
(21) |
The eigenstates, , and the eigenvalues, , are found from the Schrödinger equation
(22) |
with . This equation has a finite number of discrete bound states with real ’s (28) and continuum extended states with imaginary ’s, ( is real ). If is an integer there are bound states with and both bound and unbound eigenstates can be expressed in terms of elementary functions (29),
(23) |
where is the normalizing amplitude and . In our case () there are two bound states, the even ground state with () and the odd excited state with (). For relaxations conserving the electron-hole symmetry the ground state contribution to the superposition, Eq.(20), is integrated to zero because the initial non-equilibrium distribution is odd. On the contrary, the excited odd state with is the only state, which survives in Eq.(20) at , so that
(24) |
and (using Eq.(19) and Eq.(15))
(25) |
which is precisely the result of the QE approximation. Indeed expanding the QE distribution function, in powers of , we obtain
(26) |
with the same non-equilibrium correction, , as in Eq.(25).
The exact solution, Eq.(20), allows us to trace the relaxation at any energy and at any time scale. In particular Fig.1 represents the time evolution of the total number of non-equilibrium excitations (electrons plus holes), , for initial non-equilibrium distributions of the shape, , where the distribution width is varied, but the total energy, , is unchanged. For wide sources with large the number of excitations increases with time conserving the total energy in the process of their cooling. On the other hand, when most excitations at are created with the energy less than (i.e. ), their number decreases with time since their individual energies increase due to collisions with equilibrium electrons. One can trace the relaxation for any initial distribution including the case when a single photon initially creates a single electron hole pair, which cascades eventually into a large number of pairs until they get lost in the background thermal distribution.
The asymptotic behavior of response functions can be readily obtained from Eq.(20) by taking into account that only the excited bound state and the extended states with small contribute to the sum in Eq.(20), when . Substituting the extended eigenfunctions, into Eq.(20) with at small and integrating over yield
(27) |
in the saddle-point approximation. Performing the Fourier transform of Eq.(27) with respect to we find
(28) |
The same result is obtained by using a -approximation for the Boltzmann equation (13),
(29) |
which has the following solution
(30) |
Hence one can use the -approximation, Eq.(30), on the time scale much longer than the characteristic collision time, . However this approximation is inaccurate on a shorter time scale because in contrast with the exact solution, Eq.(20), it does not conserve the total energy.
Integrating Eq.(28) yields a universal time-asymptotic of the total number of electron-hole excitations,
(31) |
as also seen from Fig.1.
Importantly the characteristic e-e relaxation time is quite long due to the Pauli exclusion principle. Using realistic and eV we estimate ps at the room temperature K (see also Ref.(26)), which increases further as with cooling.
Iii Electron-phonon relaxation
Now let us consider the electron-phonon relaxation described by the e-ph collision integral,
(33) | |||||
(35) | |||||
(36) |
where is the matrix element of the deformation potential and is the distribution function of phonons with the frequency .
As in the former case of the e-e collisions we average this equation over the momentum angles using Eq.(3) and conventional units:
(38) | |||||
(39) |
Here
(40) |
is the e-ph spectral function (14), which has variation on a scale of the maximum phonon frequency but and variation only on a much larger energy scale of the order of .
Linearizing Eq.(39) with the help of Eq.(8) yields
(41) | |||
(42) | |||
Characteristic electron energies in Eq.(LABEL:collision3) are much less than the Fermi energy, so that
(44) |
is the familiar Eliashberg function. We also assume that phonons are in the thermal equilibrium, , due to their fast thermalization caused by anharmonic interactions (i.e. phonon-phonon collisions) and/or due to a small size of the sample and the pump-laser spot allowing for a fast escape of non-equilibrium phonons. If this condition is not satisfied, one has to solve an equation for the non-equilibrium phonon distribution coupled with Eq.(LABEL:collision3), which is outside the scope of this paper. Under these assumptions Eq.(LABEL:collision3) is transformed into a form similar to the e-e collision integral Eq.(13),
(45) | |||||
(46) |
where
(47) | |||||
and energies are now measured in units of . The Eliashberg function is quite complicated in real metallic compounds because of their complex lattice structures. This complexity can be avoided in a high-temperature regime, and in an opposite low-temperature regime, .
iii.1 High-temperature electron-phonon relaxation
As shown by Allen (14) the energy relaxation in TTM has a particularly simple form in terms of the moments of ,
(49) |
where the coupling constant , which determines the critical temperature in the BCS superconductors, is
(50) |
Here we derive a high-temperature LFP-type equation for the non-equilibrium part of the distribution function to compare our exact approach with the TTM results (13); (14).
At high temperatures the Eliashberg function is a narrow function on the temperature scale, so that one can apply a quasi-elastic approximation expanding the e-ph collision integral, Eq.(LABEL:collision3) or Eq.(46), in powers of the phonon energy, . The zero-order elastic terms are canceled out because the distribution function depends on energy only, while the next order terms yield the LFP-type differential equation,
(51) |
where . Apart from a numerical coefficient of the order of the characteristic e-ph relaxation rate is about the same as the TTM energy relaxation rate (14), , at high temperatures. Indeed multiplying Eq.(51) by and integrating over all energies yield the rate of excitation energy relaxation,
(52) |
where . If we replace in this equation by its argument assuming that has its characteristic energy width of the order of , then Hence the excitation energy relaxes as almost independent on a particular shape of the non-equilibrium distribution.
To verify the numerical coefficient, one can substitute the TTM distribution (Eq.(26)) into Eq.(52) to convert this into the temperature relaxation rate:
(53) |
Eq.(53) is precisely the same as the TTM temperature rate (14), , since the ratio of two integrals in Eq.(53) is .
According to Eq.(52) deviation of from quasi-equilibrium population does not have much influence on the energy relaxation. Hence TTM (13); (14) is the adequate approximation at high temperatures, which agrees well with experimental observations in conventional metals where Debye temperatures are rather low (1).
iii.2 Low-temperature electron-phonon relaxation in poor metals
Characteristic phonon frequencies are exceptionally high in many advanced materials like copper oxides, K, so that the low-temperature regime, is of great importance. Since all dimensionless energies in Eq.(46) are of the order of unity one can apply a low-frequency asymptotic of in this regime. The exponent depends on impurities, disorder, and sample dimensions: in clean bulk crystals while in disordered metals due to a phonon damping (30); (31) and in metallic films (32). Then Eq.(46) can be Fourier-transformed into the Schrödinger equation using the Fourier transform Eq.(19) of .
In the poor-metal case the equation for is almost the same as in the e-e case, Eq.(16), apart from a numerical coefficient in front of the integral term,
(54) | |||
(55) |
where the integrals and have been used. The difference in the numerical coefficients in front of the integral terms originates in different statistics of scatterers, which are bosons in the e-ph case and fermions in the e-e case. At low temperatures the e-ph relaxation time has the same energy and temperature dependence as the e-e relaxation time Eq.(14),
(56) |
We also notice that the temperature dependence of the e-ph relaxation rate at low temperatures, is qualitatively different from its temperature dependence at high temperatures .
The Fourier transform of Eq.(55) yields the Schrödinger-like equation
(57) |
where
(58) |
Different from the e-e case [Eq.(22) with ] the steady-state Schrödinger equation
(59) |
has only one bound (ground) state, , and itinerant states with ,
(60) |
in the e-ph case (). Only itinerant states contribute to the superposition Eq.(20) and determine the time relaxation of the distribution function because the contribution of the even ground state is integrated to zero and there is no excited odd state here. As the result the non-equilibrium part of the distribution function and the number of excitations relax with characteristic time to zero rather than to any quasi-equilibrium state as shown in Fig.3 and Fig.4 by lower curves (). Their time asymptotic is found using the saddle-point approximation as in the case of the e-e collisions,
(61) |
and
(62) |
The time evolution of is widely independent of the width of the initial distribution function at as one can see comparing the lowest curves in Fig.3 and Fig.4.
iii.3 Low-temperature electron-phonon relaxation in clean metals
The low-frequency Eliashberg function is quadratic as a function of frequency, in clean crystalline metals, which makes an analytical expression for the Fourier transform of the Boltzmann equation (46) unavailable in terms of elementary functions. However we can approximate all relevant integrals numerically as
and
where is shown in Fig.2. Then the corresponding Schrödinger-type equation for the Fourier transform of becomes
(63) |
where now
(64) |
The effective ”potential” energy differs only marginally from the poor metal case, Fig.2. At large corresponding to large in Eq.(63) the forth derivative of the low-energy extended eigenstates is small. Hence the asymptotic behavior of response functions in clean metals is qualitatively about the same as in poor metals,
(65) |
but the temperature dependence of the e-ph relaxation time is more pronounced, . In principle, the clean-metal ”potential”, Fig.2, could have ”resonances”, states that are in the continuum but take a long time to leak out resulting in some quantitative differences with the poor-metal relaxation.
Iv Low-temperature electron-phonon relaxation combined with electron-electron relaxation
Finally let us combine both collision integrals into one Boltzmann equation. Performing its Fourier transformation as described above in Sections (II, III) yields the following Schrödinger-like equation in the poor-metal case:
(66) |
where
(67) |
and
(68) |
There are two bound states of the corresponding steady-state Schrödinger-like equation
(69) |
because is larger than but smaller than (28). The excited odd state has the eigenfunction
(70) |
and the energy
(71) |
which determines the asymptotic behavior of all linear response functions. In particular the number of excitations decays at large as
(72) |
When both relaxations are involved the time evolution of calculated using Eq.(66) with different initial distributions differs qualitatively from TTM relaxation as shown in Fig.3 and Fig.4. In fact electrons cannot attain the thermal quasi-equilibrium at any less than in agreement with the numerical results of Ref.(26). Moreover the exact relaxation rate depends on the ratio of the electron-electron relaxation time, Eq.(18), and the electron-phonon relaxation time, Eq.(58),
(73) |
Using Eq.(72) we find
(74) |
where
(75) |
This coefficient changes from at up to at .
The TTM relaxation rate at low temperatures is readily obtained with the Eliashberg function using Eqs. (4,10) of Ref.(14); (33). Linearizing Eq.(10) of Ref. (14) with respect to the temperature difference yields
(76) |
The ratio of our exact relaxation rate to the TTM rate is
(77) |
If e-e collisions are much faster than e-ph collisions (), this ratio is , justifying the TTM approximation also at low temperatures in the limit . However at low temperatures is not necessarily small as assumed in TTM even at small because the Fermi energy in Eq.(73) is often much larger than the phonon energy. Just the opposite limit is feasible at a sizable . In this limit the exact relaxation rate is slower than the low-temperature TTM rate, , so that one may underestimate the electron-phonon coupling constant by about two times using TTM. Also an illegitimate fitting of experimental rates measured at temperatures below with the theoretical high-temperature TTM rate (14) (see section III.1) may underestimate by about times in poor metals and much more in clean metals.
V Conclusions
In conclusion, using the auxiliary Fourier transform we have mapped the linearized Botzmann equation with the electron-electron collision integral onto a Schrödinger-like equation with imaginary time allowing for a simple analytical solution. A similar mapping is also found for the electron-phonon collision integral at low-temperatures both in poor and clean metals. We have analytically traced the time and energy evolution of the non-equilibrium electron distribution function on any time scale and found its asymptotic relaxation rate at .
A low-temperature relaxation rate strongly depends on the temperature: and in poor and clean metals, respectively. The Pauli exclusion principle slows down e-e relaxation, so that e-e and e-ph collisions are strongly entangled at low temperatures. We have shown that electron gas cannot attain a thermal quasi-equilibrium distribution by e-e collisions, Figs.3,4, in agreement with earlier numerical integrations of the Boltzmann equation (26). The rate of return to the equilibrium is not governed solely by electron-phonon processes, but also involves the electron-electron relaxation time, , via the coefficient which depends on the ratio of the e-e collision time to the e-ph collision time. The exact relaxation rate recovers its quasi-equilibrium TTM value only in the limit of the negligible e-ph coupling, . In poor metals the physically realistic ratio is large at low temperatures and the exact relaxation rate is slower than the TTM rate, .
At high temperatures, , we have reduced the e-ph collision integral to the differential Landau-Fokker-Plank form. Using this form we have shown that the deviation of the electron distribution from quasi-equilibrium population does not have much influence on the energy relaxation, so that TTM (13); (14) is a reliable approximation at high temperatures.
Our theory opens up a perspective of determinations of both important microscopic parameters and using single-parameter fitting of response functions in pump-probe experiments at low temperatures, Figs. 3,4. It also allows for an analytical approach to the integral Boltzmann equation for the case of a steady-state source of excitations as in a current-carrying state. In the latter case relaxation times could be different because the current carrying state does not have a distribution function that depends only on energy, as assumed here. The theory could be further extended beyond the assumption that phonons remain in equilibrium by including a linearised Boltzmann equation for the non-equilibrium phonon distribution function.
We are grateful to Alexander Veselov and Rajmund Krivec for illuminating discussions of Calogero-Moser-Sutherland equations and to Dragan Mihailovic for sharing with us his insight into non-equilibrium phenomena in metals and high-temperature superconductors. The work was supported by the Slovenian Research Agency (ARRS) (grant no. 430-66/2007-17) and by EPSRC (UK) (grant no. EP/D035589/1).
References
- S. D. Brorson, A. Kazeroonian, J. S. Moodera, D. W. Face, T. K. Cheng, E. P. Ippen, M. S. Dresselhaus, and G. Dresselhaus, Phys. Rev. Lett. 64, 2172 (1990).
- R. W. Schoenlein, W. Z. Lin, J. G. Fujimoto, and G. L. Eesley, Phys. Rev. Lett. 58, 1680 (1987).
- H. E. Elsayed-Ali, T. B. Norris, M. A. Pessot, and G. A. Mourou, Phys. Rev. Lett 58, 1212 (1987).
- R. H. M. Groeneveld, R. Sprik, and A. Lagendijk, Phys. Rev. Lett. 64, 784 (1990).
- S. D. Brorson, J. G. Fujimoto, and E. P. Ippen, Phys. Rev. Lett. 59, 1962 (1987).
- G. L. Eesley, J. Heremans, M. S. Meyer, G. L. Doll, and S.H. Liou, Phys. Rev. Lett. 65, 3445, (1990).
- S. G. Han, Z. V. Vardeny, K. S. Wong, O. G. Symko, and G. Koren, Phys. Rev. Lett. 65, 2708 (1990).
- S. V. Chekalin, V. M. Farztdinov, V. V. Golovlyov, V. S. Letokhov, Yu. E. Lozovik, Yu. A. Matveets, and A. G. Stepanov Phys. Rev. Lett. 67, 3860 (1991).
- W. Albrecht, Th. Kruse, and H. Kurz, Phys. Rev. Lett. 69, 1451 (1992).
- C. J. Stevens, D. Smith, C. Chen, J. F. Ryan, B. Podobnik, D. Mihailovic, G. A. Wagner, and J. E. Evetts, Phys. Rev. Lett. 78, 2212 (1997).
- J. Demsar, B. Podobnik, V. V. Kabanov, T. Wolf, and D. Mihailovic, Phys. Rev. Lett. 82, 4918 (1999).
- V. V. Kabanov, J. Demsar, B. Podobnik, and D. Mihailovic, Phys. Rev. B 59, 1497 (1999).
- M. I. Kaganov, I. M. Lifshits, and L. B. Tanatarov, Zh. Eksp. Teor. Fiz., 31, 232, (1956) [Sov.Phys. JETP 4, 173 (1957)].
- P. B. Allen, Phys. Rev. Lett. 59, 1460 (1987).
- A. Rothwarf and B. N. Taylor, Phys. Rev. Lett., 19, 27 (1967).
- G. P. Segre, N. Gedik, J. Orenstein, D. A. Bonn, R. Liang, and W. N. Hardy, Phys. Rev. Lett. 88, 137001 (2002); N. Gedik, P. Blake, R. C. Spitzer, J. Orenstein, R. Liang, D. A. Bonn, and W. N. Hardy, Phys. Rev. B 70, 014504 (2004).
- J. Demsar, R. D. Averitt, V. V. Kabanov, and D. Mihailovic, Phys. Rev. Lett. 91, 169701 (2003).
- V. V. Kabanov, J. Demsar, D. Mihailovic Phys. Rev. Lett. 95, 147002 (2005).
- R. A. Kaindl, M. A. Carnahan, D. S. Chemla, S. Oh, and J. N. Eckstein, Phys. Rev. B 72, 060510(R) (2005).
- L. Perfetti, P. A. Loukakos, M. Lisowski, U. Bovensiepen, H. Eisaki, and M. Wolf, Phys. Rev. Lett. 99, 197001 (2007).
- J. X. Zhu, E. E. M. Chia, T. Tamegai, H. Eisaki, K. H. Oh, S. I. Lee, and A. J. Taylor, arXiv:0806.266.
- A. S. Alexandrov and N. F. Mott, Rep. Prog. Phys. , 1197 (1994) and references therein.
- P. W. Anderson, Science 317, 1705 (2007) and references therein.
- J. Demsar, R. D. Averitt, K. H. Ahn, M. J. Graf, S. A. Trugman, V. V. Kabanov, J. L. Sarrao, and A. J. Taylor, Phys. Rev. Lett. 91, 027401 (2003).
- K. H. Ahn, M. J. Graf, S. A. Trugman, J. Demsar, R. D. Averitt, J. L. Sarrao, and A. J. Taylor, Phys. Rev. B 69, 045114 (2004).
- R. H. M. Groeneveld, R. Sprik, and Ad Lagendijk, Phys. Rev. B 51, 11433 (1995).
- V. I. Karas and I. E. Potapenko, Zh. Vychislitel’noi Matematiki i Matematicheskoi Fiziki 46, 307 (2006) [Computational Mathematics and Mathematical Physics 46, 294 (2006)].
- L. D. Landau and E. M. Lifshitz, Quantum Mechanics (Pergamon Press, Oxford 1977).
- We thank A. P. Veselov for a hint, which have led us to this expression.
- D. Belitz, Phys. Rev. B 36, 47, 2513 (1987).
- D. Belitz and M. N. Wybourne, Phys. Rev. B 51, 689 (1995)
- B. I. Belevtsev, Yu. F. Komnik, and E. Yu. Beliayev, Phys. Rev. B 58, 8079 (1998).
- We have noticed that a multiplier is missing in the second term of Eq.(4) in Ref. (14).