Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Phonon self energy

  1. Jan 21, 2010 #1
    Hello guys... Can anyone tell me where I can find a detailed calculation of the phonon self energy for an electron phonon interaction in graphene... The expression I need to get is in T. Ando paper:

    Journal of the Physical Society of Japan
    Vol. 75, No. 12, December, 2006, 124701

    Equation (3.4)....

    or even just hiow to get the self energy (phonon) fr an electron phonon interaction using perturbation theory....

    Thanks very much
  2. jcsd
  3. Jan 22, 2010 #2


    User Avatar
    Science Advisor

    Maybe you can give us the formula, as I do not have access to that journal.
    Generally in lowest order in the coupling constant, the self energy is due to the phonon getting anihilated in a process where an electron is excited. At a later time, the electron hops back to the ground state and the phonon re-emerges.
    Hence the self energy of an (eventually virtual) phonon of energy omega and momentum p is something like (up to constants) [tex] g^2 \int dk d \epsilon G(\omega+\epsilon, k+p)G(-\epsilon,-k)[/tex]
    where g is the electron phonon coupling constant and G is the electronic Greens function (which should be of Dirac form near the conical points.)
  4. Jan 25, 2010 #3
    Thanks very much for the reply. I am trying to get the self energy for graphene (monolayer). I will give you the formula...I just cant get my greens functions right.... The result they get in paper is

    [tex]\Pi[/tex](q,[tex]\omega[/tex])=C[tex]~\sum_{s,s_{1}}[/tex] [tex]\int[/tex][tex]\frac{dk}{(2\pi)^2}[/tex][tex]\frac{f[\epsilon(s,k)]-f[\epsilon(s_{1},k)]}{\hbar\omega-\epsilon(s,k)+\epsilon(s_{1},k)+i0}[/tex]

    where [tex]f[\epsilon(s,k)][/tex] and [tex]f[\epsilon(s_{1},k)][/tex]are the fermi dirac statistics for electrons in a state k in bands s and s´(conduction and valence bands) and C are some constants which i managed to calculate. The e-phonon interaction Hamiltonian is:

    H=[tex]\sum c^{+}_{k+q}c_{k}(b_{q}+b^{+}_{-q})[/tex]
    where c are electronic operators and b are phonon operators...The self energy is coming from retarded Green functions but i am trying to get matsubara Green function and I obtain strange results and I do not get energies within different bands...Can anyone please help..Thanksss
    Last edited: Jan 25, 2010
  5. Jan 27, 2010 #4


    User Avatar
    Science Advisor

    I managed to have a look at that paper. However, I have never done any Matsubara calculation myself, so that I can offer only very limited help.
    First we need the electronic Greens functions. They should be of the form [tex](\epsilon-\epsilon_s)^{-1} [/tex] and [tex](\epsilon-\epsilon_s'+\hbar \omega)^{-1} [/tex] (very schematically and dropping all infinitesimal imaginary constants). We need the product of the two GF's and have to integrate over epsilon.
    To do this, the product has to be written as a sum of two partial fractions. When doing this partial fraction decomposition, I get a term [tex](\hbar\omega-\epsilon(s',k)+\epsilon(s,k))^{-1}[/tex] as appears in the denominator of 3.2. The integrals over the partial fractions yield probably the Fermi occupation probabilities (although I would have expected an [tex]\hbar \omega [/tex] to appear in one of these probabilities), so that 3.2 does not appear implausible. I hope that helps somehow.
  6. Jan 27, 2010 #5
    Thanks very much for your help. I managed to get the Green functions and drew the Feynmann diagrams and discarded the disconnected ones and I got a product of two Green functions of the form

    G([tex]k lambda [/tex], [tex]\tau[/tex]-[tex]\tau1[/tex])*G([tex]k+q \lambda1[/tex], [tex]\tau1[/tex]).

    where (k [tex] \lambda [/tex]) is a state of electron with momentum k and in valence band [tex]\lambda[/tex], and the other one is an a electron in state k+q and conduction band [tex]\lambda1[/tex],. of course [tex]\tau[/tex] and [tex]\tau1[/tex] are the variables of the Fourier Transform.

    Now I need to Fourier transform this to the Matsubara frequency. I cannot manage to do this. I know we need the convolution theorem but how do we appoly it here because we have got two times [tex]\tau[/tex] and [tex]\tau1[/tex]. Then I guess we can apply the partial fraction method and use sums over frequencies or Residue theorem. I am happy to get the Fourier transform now as I can continue the rest of the calculations. If you have any suggestions or a method to get this Fourier transform, please let me know..

    Thanks in advance again

    Last edited: Jan 27, 2010
  7. Jan 27, 2010 #6


    User Avatar
    Science Advisor

    Why don't you work direktly with the Greensfunktions and rules in frequency space?
    But anyhow, what we want to describe is this: a phonon anihilates into an electron and a hole at time [tex]\tau_1[/tex] and resurrects at time [tex]\tau_2[/tex] hence both Greens functions depend on the same time difference [tex] \tau_2-\tau_1[/tex]. Using this time-difference as a new variable, the Fourier transform gives you a convolution of the two Greensfunctions in frequency space, that is, the equivalent of what is written in my post #2.
  8. Jan 28, 2010 #7


    User Avatar
    Science Advisor

    Last evening I had a look at Fetter, Walecka, "Quantum Theory of Many Body Systems", Dover. It contains the calculation of the self energy Pi_0. The Matsubara Greens functions result principally from the ones of my post #4 upon replacing [tex]\epsilon [/tex] with [tex] i \nu [/tex] where [tex] \nu [/tex] are the Matsubara frequencies (the Fourier transform of tau). The sum over the Matsubara frequencies can be evaluated after partial fraction decomposition, but one has to introduce some exponential convergence factor for the sums to be well defined. They become the occupation probabilities as said above.
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook