<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
Well, there are several ways to skin the cat.<br>
1. The most correct way is to calculate "aph2F(w)". How you do this:
you have to numerically get double integral over the calculated Fermi
surface for the function <k|M(k,k')|k'>*ds*ds'. You will need
partial <k,etha| for each point at FS, where etha is relativistic
quantum number. Many years ago I have written the code for RAPW wave
functions basis. Took me 2 years to write matrix element and code for
transition metals, and to get it working. But computers were slow then,
and it was RAPW (non-linearized!). RLMTO and LAPW were not born yet.
Phonon frequency appears through |k-k'|=q in the matrix element.<br>
2. If you don't have heavy metals in your system, you can get away
without relativistic effects. It will make your life easier, but you
still need double integral over the real FS, and real phonon spectra.
General ideas may be found in Ziman's "Electrons and Phonons". Oldest
but goldest.<br>
3. If you don't have substantial p-, f- and d-electron state density at
Ef, (s-electrons->spherical FS), you can use pseudopotential matrix
element, and double integral converts into single integral over dq3 for
expanded zone. You can pull it out using tetrahedron, or whatever BZ
integration method. Much easier, I think it is in Ziman too.<br>
4. You may voluntary ignore the fact that "alpha2" really depends on
omega. It does, but most people don't even think about writing so
involved code. The matter is decided by balance of time at hand,
ambition, braveness and laziness of specific individual. Anyway, if you
believe that your "alpha2" is practically constant, you can use
Hopfield-McMillan approach. It lets you to derive lambda and
lambda(k,k'), but there is some risk (at very extreme this may cost you
20-25% error for lambda). All you will need in this case will be
partial state densities at Ef. You simply circumvent all phonon
dependence of "apha2", and bravely use <w> and <w2>.<br>
5. Finally, if you are extremely lucky, you may find experimental
tunnel spectra or point-contact spectra for your compound. But it is
really a rarity.<br>
Now, where to read more. First of all, look at Ziman. Going beyond that
: it is not a self-promotion, it is just easier for me to do it this
way. Go to <a class="moz-txt-link-freetext" href="http://www.azgre.com/PERSONAL/zhtpersonal.htm">http://www.azgre.com/PERSONAL/zhtpersonal.htm</a>, and from the
list of publications look at 4., 8., 9., 13., 14., 57, and 59 (this one
is hard to find, it's on "pseudizm" in EPI). Look at references in
these articles.<br>
Good luck, and please don't hesitate to ask more questions.<br>
W/regards<br>
Dr. A.Zhako-Titarenko <br>
<blockquote
cite="mid:20090412182253.46206.qmail@f4mail-234-240.rediffmail.com"
type="cite">
<p> <br>
Dear PWSCF users,<br>
<br>
1. I want to calculate the electron-phonon coefficient "lamda" and the<br>
function "square of alpha" F(w) plot, how i do?<br>
<br>
2. I also want to calculate IR spectra and Raman spectra?<br>
<br>
Thanks in advance!!<br>
<br>
Kind Regards<br>
Sanjeev Gupta
</p>
S. K. Gupta
<br>
PhD Candidate,
<br>
Dept. of Physics,
<br>
Bhavnagar University.
<br>
Bhavnagar, 364 002., Gujarat.
<br>
<br>
<br>
<table
style="font-family: Verdana; font-size: 11px; line-height: 15px;"
border="0" cellpadding="0" cellspacing="0" height="57" width="644">
<tbody>
<tr>
<td><a moz-do-not-send="true"
href="http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?"
target="_blank"><FONT size="1">[AD]</FONT></a></td>
</tr>
</tbody>
</table>
<pre wrap="">
<hr size="4" width="90%">
_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
</blockquote>
<br>
</body>
</html>