<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
May or may not help:<br>
<br>
RINDEX = 11 => r**krad for krad > 10 => :XOP 104670.12
(large number)<br>
<br>
Adding a negative sign will change the calculation:<br>
<br>
RINDEX = -11 => 1/r**krad for krad < 10 => :XOP ? (small
number)<br>
<br>
File of reference: $WIENROOT/SRC_lapwdm/radint.f<br>
<br>
On 6/7/2012 7:08 AM, David Tompsett wrote:
<blockquote
cite="mid:CALGV4hME+S4jz+ZvET7wxpwiWTmLkrdn7m=JSbnTvnJ8x7F0zQ@mail.gmail.com"
type="cite">Hi Peter,<br>
<br>
Thanks for the response. The option r-index = 11 is not in the UG,
but is documented in Pavel Novak's technical report on the
Hyperfine Field calculation so other users may strike it in
future.<br>
<br>
Hopefully he will have some input.<br>
<br>
Thanks,<br>
David.<br>
<br>
<div class="gmail_quote">On Thu, Jun 7, 2012 at 7:47 AM, Peter
Blaha <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">You have to
check directly in the code. Without looking in the code I do
not know what r-index = 11 should do ? I guess this is not a
documented option, but Pavel Novak migth have put in something
....<br>
<br>
Am 06.06.2012 14:11, schrieb David Tompsett:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div class="h5">
Dear All,<br>
<br>
I have been attempting to use lapwdm to obtain the
expectation value of<br>
the radial distribution of orbitals i.e. <r>.<br>
<br>
I have attempted the simplest case that I could think
of, the hydrogen<br>
atom in a large cell. I used an RMT of 4.5 au so that a
fraction 0.99 of<br>
the charge lies inside the muffin tin. Then after
converged SCF cycles I<br>
ran "x lapwdm -up" with input:<br>
<br>
-9. Emin cutoff energy<br>
1 number of atoms for which density matrix is
calculated<br>
1 1 0 index of 1st atom, number of L’s, L1<br>
11 1 r-index, (l,s)-index<br>
<br>
From the solution of the Schroedinger equation for a
Hydrogen atom, the<br>
result should be <r> = 1.5 au.<br>
However the result from :XOP in the output of lapwdm is
104670.12. I do<br>
not know what the units are. Also, the value of this
result seems to<br>
vary strongly with the RMT.<br>
<br>
Can anyone help me understand this behaviour or if there
is a problem<br>
with lapwdm? I looked through the source code, but could
not follow it.<br>
<br>
Many thanks,<br>
David.<br>
<br>
<br>
</div>
</div>
_______________________________________________<br>
Wien mailing list<br>
<a moz-do-not-send="true"
href="mailto:Wien@zeus.theochem.tuwien.ac.at"
target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a moz-do-not-send="true"
href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien"
target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><span
class="HOEnZb"><font color="#888888"><br>
</font></span></blockquote>
<span class="HOEnZb"><font color="#888888">
<br>
-- <br>
Peter Blaha<br>
Inst.Materials Chemistry<br>
TU Vienna<br>
Getreidemarkt 9<br>
A-1060 Vienna<br>
Austria<br>
<a moz-do-not-send="true" href="tel:%2B43-1-5880115671"
value="+4315880115671" target="_blank">+43-1-5880115671</a><br>
_______________________________________________<br>
Wien mailing list<br>
<a moz-do-not-send="true"
href="mailto:Wien@zeus.theochem.tuwien.ac.at"
target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a moz-do-not-send="true"
href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien"
target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</font></span></blockquote>
</div>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Wien mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
</pre>
</blockquote>
<br>
</body>
</html>