Dear All,<br>I want to calculate the cohesive energy of graphite. I am taking C-pz-vbc.UPF for graphite and C-pz.rrkjus.UPF for isolated carbon atom, is this acceptable?<br>Sincerely,<br>Thanehwor Kaloni<br>SNBNCBS<br><a href="http://www.bose.res.in">www.bose.res.in</a><br>
<br><div class="gmail_quote">On Tue, Jun 15, 2010 at 11:37 PM,  <span dir="ltr">&lt;<a href="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Send Pw_forum mailing list submissions to<br>
        <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:pw_forum-owner@pwscf.org">pw_forum-owner@pwscf.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of Pw_forum digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: Computing Second order response (loc duong ding)<br>
   2. pseudo-potential for virtual atom (German Samolyuk)<br>
   3. Re: total potential plot (Matteo Cococcioni)<br>
   4. Re: Computing Second order response (Stefano Baroni)<br>
   5. Problem in reading PP&#39;s taken from QE site? (sonu kumar)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 15 Jun 2010 05:56:10 -0700 (PDT)<br>
From: loc duong ding &lt;<a href="mailto:mambom1902@yahoo.com">mambom1902@yahoo.com</a>&gt;<br>
Subject: Re: [Pw_forum] Computing Second order response<br>
To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
Message-ID: &lt;<a href="mailto:508640.34950.qm@web38802.mail.mud.yahoo.com">508640.34950.qm@web38802.mail.mud.yahoo.com</a>&gt;<br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
I don&#39;t understand why Raman activity in this case is zero. In practical, metallic system also provide Raman peak (ex. Carbon nannotube).  I try to calculate graphene (no gap) and I get the proper result. Can you give some more explanation or suggest me some reference about this topic?<br>

<br>
 I appreciate to your instructions.<br>
<br>
&gt;So, the the Raman activity should be<br>
&gt; zero ...<br>
&gt;SB<br>
<br>
On Jun 15, 2010, at 3:50 AM, loc duong ding<br>
&gt; wrote:<br>
<br>
&gt; I calculate zigag nano ribbon with hydrogen terminated at the<br>
&gt; edge. I checked DOS that showed no gap.<br>
&gt;<br>
&gt;&gt; yes - SB<br>
&gt;<br>
&gt;<br>
&gt; On Jun 14, 2010, at<br>
&gt;&gt;&gt; 11:03 AM, loc duong ding<br>
&gt; wrote:<br>
&gt;<br>
&gt;&gt;&gt; what&#39;s the gap of the<br>
&gt;&gt; unperturbed<br>
&gt;<br>
&gt;&gt;&gt; system?<br>
&gt;&gt;&gt; SB<br>
&gt;&gt;<br>
&gt;&gt; Actually,<br>
&gt; this<br>
&gt;&gt; is the firs time I calculate the Raman spectrum. I don&#39;t clear<br>
&gt; about the gap of<br>
&gt;&gt; the unperturbed system? It means the electronic<br>
&gt; band gap?<br>
<br>
<br>
 -----------------------------------------------<br>
Loc Duong Dinh<br>
Ms-Ph.D Student<br>
Sungkyunkwan Advanced Institute of Nanotechnology,<br>
Sungkyunkwan University,<br>
Suwon, 440-746, Korea<br>
Email: <a href="mailto:mambom1902@yahoo.com">mambom1902@yahoo.com</a>, <a href="mailto:ddloc@skku.edu">ddloc@skku.edu</a><br>
<br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 15 Jun 2010 10:27:43 -0400<br>
From: German Samolyuk &lt;<a href="mailto:samolyuk@gmail.com">samolyuk@gmail.com</a>&gt;<br>
Subject: [Pw_forum] pseudo-potential for virtual atom<br>
To: PWSCF Forum &lt;<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:AANLkTilUQNs7ugkqa6DgZzxacqRhOqhcDe0j1Zfi7eNG@mail.gmail.com">AANLkTilUQNs7ugkqa6DgZzxacqRhOqhcDe0j1Zfi7eNG@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Dear all,<br>
<br>
I need to produce pseudo-potential for virtual atom, something like<br>
W_{1-x}Re_x. The pseudo-potentials from the database,<br>
<br>
W.pw91-nsp-van.UPF and Re.pw91-nsp-van.UPF, have different number of valence<br>
electrons states and utility vitual.x does not<br>
work for this two ps because of this reason. However, if I make ps for these<br>
two elements with the same number of wave functions using ld1.x,<br>
virtual.x can not read ps files because UPF format has been changed. I tried<br>
versions 4.1, 4.1.1  and 4.2.<br>
What should I do to solve this problem?<br>
Thank you in advance.<br>
<br>
German Samolyuk,<br>
Oak Ridge National Laboratory, USA<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20100615/05a128e2/attachment-0001.htm" target="_blank">http://www.democritos.it/pipermail/pw_forum/attachments/20100615/05a128e2/attachment-0001.htm</a><br>

<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 15 Jun 2010 10:34:34 -0500<br>
From: Matteo Cococcioni &lt;<a href="mailto:matteo@umn.edu">matteo@umn.edu</a>&gt;<br>
Subject: Re: [Pw_forum] total potential plot<br>
To: PWSCF Forum &lt;<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:4C179D8A.6050503@umn.edu">4C179D8A.6050503@umn.edu</a>&gt;<br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
<br>
<br>
Thanks Paolo.<br>
<br>
<br>
Paolo Giannozzi wrote:<br>
&gt; On Jun 14, 2010, at 22:50 , Matteo Cococcioni wrote:<br>
&gt;<br>
&gt;<br>
&gt;&gt; (V_ext + V_H + V_xc)<br>
&gt;&gt;<br>
&gt;<br>
&gt; V_ext = nuclear (pseudo-)potential acting on electrons?<br>
&gt;<br>
&gt;<br>
<br>
yes, that&#39;s what I meant. In fact I wasn&#39;t sure this corresponds to<br>
&quot;V_bare&quot; or to &quot;local ionic potential&quot;.<br>
I will use plot_num = 1.<br>
<br>
regards,<br>
<br>
Matteo<br>
<br>
&gt;&gt;                       1  = total potential V_bare+V_H + V_xc<br>
&gt;&gt;<br>
&gt;<br>
&gt; then this is what you need.<br>
&gt;<br>
&gt;<br>
&gt;&gt;                       2  = local ionic potential<br>
&gt;&gt;<br>
&gt;<br>
&gt; The latest QE version is more explicit:<br>
&gt;<br>
&gt;<br>
&gt;&gt;                       2  = local ionic potential V_bare<br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;&gt; 2) is there any known reason why plotrho fails when trying to make<br>
&gt;&gt; a 2D<br>
&gt;&gt; plot of the potential?<br>
&gt;&gt;<br>
&gt;<br>
&gt; it doesn&#39;t fail. Run it interactively.<br>
&gt;<br>
&gt; P.<br>
&gt; ---<br>
&gt; Paolo Giannozzi, Dept of Physics, University of Udine<br>
&gt; via delle Scienze 208, 33100 Udine, Italy<br>
&gt; Phone +39-0432-558216, fax +39-0432-558222<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Tue, 15 Jun 2010 18:01:36 +0200<br>
From: Stefano Baroni &lt;<a href="mailto:baroni@sissa.it">baroni@sissa.it</a>&gt;<br>
Subject: Re: [Pw_forum] Computing Second order response<br>
To: PWSCF Forum &lt;<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:0953D55F-8EEF-47DA-8FE4-BB8FFF770F1F@sissa.it">0953D55F-8EEF-47DA-8FE4-BB8FFF770F1F@sissa.it</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
well, you probably know more than I do, but the cross section for (off-resonant) Raman scattering is basically the derivative of the dielectric constant with respect to the amplitude of the vibrational distortion. I do not know what exactly is measured (you probably do), but this is what is calculated, and the dielectric tensor of a metal is infinite ... in any case, it is always a good idea to know what is measured (which I don&#39;t, in this case), and what is calculated (which you do not seem to do, in this case ...).<br>

<br>
graphene may be another story, because the density of states at the fermi energy is zero ...<br>
<br>
Stefano B<br>
<br>
On Jun 15, 2010, at 2:56 PM, loc duong ding wrote:<br>
<br>
&gt; I don&#39;t understand why Raman activity in this case is zero. In practical, metallic system also provide Raman peak (ex. Carbon nannotube).  I try to calculate graphene (no gap) and I get the proper result. Can you give some more explanation or suggest me some reference about this topic?<br>

&gt;<br>
&gt; I appreciate to your instructions.<br>
&gt;<br>
&gt;&gt; So, the the Raman activity should be<br>
&gt;&gt; zero ...<br>
&gt;&gt; SB<br>
&gt;<br>
&gt; On Jun 15, 2010, at 3:50 AM, loc duong ding<br>
&gt;&gt; wrote:<br>
&gt;<br>
&gt;&gt; I calculate zigag nano ribbon with hydrogen terminated at the<br>
&gt;&gt; edge. I checked DOS that showed no gap.<br>
&gt;&gt;<br>
&gt;&gt;&gt; yes - SB<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; On Jun 14, 2010, at<br>
&gt;&gt;&gt;&gt; 11:03 AM, loc duong ding<br>
&gt;&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt;&gt;&gt; what&#39;s the gap of the<br>
&gt;&gt;&gt; unperturbed<br>
&gt;&gt;<br>
&gt;&gt;&gt;&gt; system?<br>
&gt;&gt;&gt;&gt; SB<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Actually,<br>
&gt;&gt; this<br>
&gt;&gt;&gt; is the firs time I calculate the Raman spectrum. I don&#39;t clear<br>
&gt;&gt; about the gap of<br>
&gt;&gt;&gt; the unperturbed system? It means the electronic<br>
&gt;&gt; band gap?<br>
&gt;<br>
&gt;<br>
&gt; -----------------------------------------------<br>
&gt; Loc Duong Dinh<br>
&gt; Ms-Ph.D Student<br>
&gt; Sungkyunkwan Advanced Institute of Nanotechnology,<br>
&gt; Sungkyunkwan University,<br>
&gt; Suwon, 440-746, Korea<br>
&gt; Email: <a href="mailto:mambom1902@yahoo.com">mambom1902@yahoo.com</a>, <a href="mailto:ddloc@skku.edu">ddloc@skku.edu</a><br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br>
---<br>
Stefano Baroni - SISSA  &amp;  DEMOCRITOS National Simulation Center - Trieste<br>
<a href="http://stefano.baroni.me" target="_blank">http://stefano.baroni.me</a> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)<br>
<br>
La morale est une logique de l&#39;action comme la logique est une morale de la pens?e - Jean Piaget<br>
<br>
Please, if possible, don&#39;t  send me MS Word or PowerPoint attachments<br>
Why? See:  <a href="http://www.gnu.org/philosophy/no-word-attachments.html" target="_blank">http://www.gnu.org/philosophy/no-word-attachments.html</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20100615/c8219b04/attachment-0001.htm" target="_blank">http://www.democritos.it/pipermail/pw_forum/attachments/20100615/c8219b04/attachment-0001.htm</a><br>

<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Tue, 15 Jun 2010 23:21:56 +0530<br>
From: sonu kumar &lt;<a href="mailto:1009ukumar@gmail.com">1009ukumar@gmail.com</a>&gt;<br>
Subject: [Pw_forum] Problem in reading PP&#39;s taken from QE site?<br>
To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
Message-ID:<br>
        &lt;<a href="mailto:AANLkTilFnFCb6WJWDS6RX3EH8Ul00C5DtqvhC-qJR14U@mail.gmail.com">AANLkTilFnFCb6WJWDS6RX3EH8Ul00C5DtqvhC-qJR14U@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Dear all QE users,<br>
<br>
I m having problem with reading of pp&#39;s taken from QE site<br>
<a href="http://www.quantum-espresso.org/pseudo.php" target="_blank">http://www.quantum-espresso.org/pseudo.php</a><br>
<br>
I did many scf runs e.g. on Al, AlAs, PbTiO3 etc., and found the problem in<br>
reading the pp file.<br>
<br>
Error of not finding the  pp file comes( like file ./Al.pz-vbc_s1.UPF not<br>
found), while with exactly same pp present in &quot;pseudo&quot;, no error comes. e.g<br>
with   Al, As  pp&#39;s : Al.pz-vbc.UPF, As.pz-bhs.UPF etc.<br>
<br>
I, further, checked  by reading the first two lines of pp&#39;s. PP&#39;s present in<br>
&quot;pseudo&quot; are read  but  that taken from site are not read.<br>
<br>
small trivial program used to read is:<br>
  program red<br>
      implicit none<br>
      character :: char1*10, char2*50<br>
      open(unit=1, file=&quot;As.pz-bhs.UPF&quot;, status=&quot;old&quot;)<br>
<br>
      read (1,&quot;(A10)&quot;) char1<br>
      read (1,&quot;(A50)&quot;) char2<br>
<br>
      write(*,&quot;(A10)&quot;) char1<br>
      write(*,&quot;(A50)&quot;) char2<br>
  end program red<br>
<br>
Output of above program:<br>
<br>
1)For pp As.pz-bhs.UPF in &quot; pseudo &quot; is :<br>
<br>
&lt;PP_INFO&gt;<br>
Generated using unknown code<br>
<br>
2) For EXACTLY SAME pp As.pz-bhs.UPF from site is:<br>
<br>
forrtl: No such file or directory<br>
forrtl: severe (29): file not found, unit 1, file<br>
/home/hcg/Desktop/DESKTOP_sonu/CODES/pwscf/sonu-exercise_QE/AlAs/As.pz-bhs_s1.UPF<br>
Image              PC                Routine            Line<br>
Source<br>
a.out              000000000046E6BD  Unknown               Unknown  Unknown<br>
a.out              000000000046D1C5  Unknown               Unknown  Unknown<br>
a.out              0000000000447EE9  Unknown               Unknown  Unknown<br>
a.out              000000000043596F  Unknown               Unknown  Unknown<br>
a.out              00000000004351A2  Unknown               Unknown  Unknown<br>
a.out              0000000000408FDD  Unknown               Unknown  Unknown<br>
a.out              0000000000402CC4  Unknown               Unknown  Unknown<br>
a.out              0000000000402C2C  Unknown               Unknown  Unknown<br>
libc.so.6          00007FDAD9BB6586  Unknown               Unknown  Unknown<br>
a.out              0000000000402B29  Unknown               Unknown  Unknown<br>
<br>
<br>
<br>
I don&#39;t know why the pp&#39;s are  not read by the program?<br>
<br>
I really need your help !!!<br>
Please help me to find the solution.<br>
<br>
With Kind Regards,<br>
Sonu Kumar<br>
<br>
Phd Student<br>
Physics Department<br>
Indian Institute of Technology<br>
Delhi-110016, India<br>
web:-<a href="http://www.iitd.ac.in/" target="_blank">http://www.iitd.ac.in/</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20100615/fd26a9a9/attachment.htm" target="_blank">http://www.democritos.it/pipermail/pw_forum/attachments/20100615/fd26a9a9/attachment.htm</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br>
<br>
End of Pw_forum Digest, Vol 36, Issue 45<br>
****************************************<br>
</blockquote></div><br><br clear="all"><br>-- <br>###################################<br>******************************************<br>Thaneshwor Prashad Kaloni<br>TWAS-BOSE FELLOW /JRF<br>S N Bose national Centre for Basic Sciences<br>
Kolkata 700098,<br>INDIA<br>