<br><br>
<div class="gmail_quote">On Mon, Apr 4, 2011 at 9:29 AM, <span dir="ltr">&lt;<a href="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</a>&gt;</span> wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">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. Van der Waals (Adriano Mosca Conte)<br>  2. Shift of Energy (pari shok)<br>
  3. Re: Shift of Energy (Duy Le)<br>  4. Re: 4.3 compile error. (?S)<br><br><br>----------------------------------------------------------------------<br><br>Message: 1<br>Date: Mon, 4 Apr 2011 15:27:08 +0200 (CEST)<br>From: &quot;Adriano Mosca Conte&quot; &lt;<a href="mailto:Adriano.Mosca.Conte@roma2.infn.it">Adriano.Mosca.Conte@roma2.infn.it</a>&gt;<br>
Subject: [Pw_forum] Van der Waals<br>To: <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>Message-ID: &lt;<a href="mailto:59184.141.108.248.78.1301923628.squirrel@141.108.248.78">59184.141.108.248.78.1301923628.squirrel@141.108.248.78</a>&gt;<br>
Content-Type: text/plain;charset=iso-8859-1<br><br>I have two questions about the calculations with VdW forces<br>(input_dft       = &#39;vdW-DF&#39;):<br>1- Is it possible to perform VdW calculations with any kind of functional?<br>
Or is it mandatory to use revPBE?<br>2- I generated a VdW table on two different machines (sp6 of CINECA,<br>matrix of  CASPUR). I got different numbers in the files vdW_kernel_table<br>(sign<br>and order of magnitude are different in many cases). Could it depend on the<br>
modules I loaded on matrix? Which are the libraries used to generate<br>the table?<br>Adriano<br><br><br>--<br>Universit? degli Studi Tor Vergata<br>via della Ricerca Scientifica 1<br>00133 Roma, Italia<br>Tel. +39 06 7259 4741<br>
Fax  +39 06 2023507<br><br><br>------------------------------<br><br>Message: 2<br>Date: Mon, 4 Apr 2011 10:32:53 -0400<br>From: pari shok &lt;<a href="mailto:parishok@gmail.com">parishok@gmail.com</a>&gt;<br>Subject: [Pw_forum] Shift of Energy<br>
To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>Message-ID: &lt;BANLkTi=<a href="mailto:EzJF-KgYL0Cf2Yxe45cWJte4y9g@mail.gmail.com">EzJF-KgYL0Cf2Yxe45cWJte4y9g@mail.gmail.com</a>&gt;<br>Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>Dear All,<br>Thank you very much for your help.<br>Please let me ask my final question w.r.t DOS shift of energy of SiC<br>supercell.<br>The total energy of the 8-atom cell of SiC is: -76.276 Ry, while that of<br>72-atom SiC is -694 Ry. Neglecting almost -7 Ry, the total energy is scaling<br>
accordingly.<br>However, I still see the shift of energy in DOS diagram of two slabs (the<br>max and min of bandgap shifts, although the gap is almost 2.1 eV for both<br>4H-SiC slabs). I was wondering whether this shift of energy is inevitable or<br>
it is as the result of my input file.<br>I kept k points (w.r.t. crystal), and cut off energy the same for both<br>structures. Should I scale the cut off energy as well. If positive, is there<br>any relation between scaling the cell and cut off energy.<br>
As you know, I need to find a reference for my further calculations.<br>I really really appreciate your help.<br>Yours<br>P Shok<br><br>Dear  All,<br>    let me clarify a little point.<br>    The potential (hence the eigenvalues) of  any calculation in<br>
periodic boundary conditions, not just a plane-wave one , is defined up<br>to  an additive constant... In a finite system a reeference can be fixed<br>by requiring that the potential goes to zero at infinite distance but in<br>
a periodic system this cannot be defined as there is nothing like a<br>point at infinite distance ...<br>   In most codes (and in pw.x) this arbitrary constant is fixed by<br>setting the average electrostatic potential to zero. If the supercell is<br>
properly constructed, with the properly scaled dimensions, the same<br>cutoff, equivalent k-points etc.. the two calculations should be<br>equivalent and the eigenvalues should be exactly mapped, the average<br>potential should be the same, the total energy should be an exact<br>
multiple of the one of the original cell.<br><br>If this does not happen one is NOT doing the supercell calculation properly.<br><br>So if when doing a supercell you do not get the scaled result check your<br>system definition.<br>
<br>stefano<br><br><br>On 04/04/2011 11:46 AM, Giuseppe Mattioli wrote:<br>&gt;* Dear all<br>*&gt;* I suspect that a a little misunderstanding has been going on here...<br>*&gt;* I try to reformulate the question: in the case of a 8-atoms cell the valence<br>
*&gt;* band maximum (VBM) has been found around 10 eV and the conduction band<br>*&gt;* minimum (CBM) around 12 eV with a 2 eV band gap. So there is a shift in the<br>*&gt;* position of VBM and CBM when calculated in a 72-atoms supercell,<br>
but the band<br>*&gt;* gap value is the same. Then, the answer to P Shok could be: plane wave<br>*&gt;* pseudopotential calculations do not ensure an universal alignment of<br>*&gt;* eigenvalues. You should use some &quot;internal&quot; reference, like the 1s<br>
eigenvalue<br>*&gt;* of an He atom, in your cell; or trust that the two VBM values are<br>aligned &quot;de<br>*&gt;* facto&quot;.<br>*&gt;* HTH<br>*&gt;*<br>*&gt;* Giuseppe<br>*&gt;*<br>*&gt;* On Monday 04 April 2011 09:39:18 Gabriele Sclauzero wrote:<br>
*&gt;&gt;* Yes,<br>*&gt;&gt;*<br>*&gt;&gt;*     I perfectly agree with Eyvaz! Moreover the experimental bandgap is<br>*&gt;&gt;* underestimated within LDA or GGA, and for 4H-SiC you would get about<br>*&gt;&gt;* 2.2-2.3 eV with LDA.<br>
*&gt;&gt;*<br>*&gt;&gt;*     How do you measure the bandgap? Beware that in the fundamental cell the<br>*&gt;&gt;* bandgap is not a direct one, but it might become such when you use a<br>*&gt;&gt;* supercell (because of the refolding of the Brillouin zones).<br>
*&gt;&gt;*<br>*&gt;&gt;*<br>*&gt;&gt;* HTH<br>*&gt;&gt;*<br>*&gt;&gt;* GS<br>*&gt;&gt;*<br>*&gt;&gt;* Il giorno 03/apr/2011, alle ore 10.59, Eyvaz Isaev ha scritto:<br>*&gt;&gt;&gt;* Hi,<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;&gt;* The bandgap of 8-atoms is around 10-12 eV,<br>
*&gt;&gt;&gt;&gt;* but that of 72-atom is around 5-7 eV.<br>*&gt;&gt;&gt;* This is unbelievable, as SiC is an indirect band gap semiconductor and<br>*&gt;&gt;&gt;* the bandgap is around (2.5 - 3)eV depending on crystal modification<br>
*&gt;&gt;&gt;* (3C-SiC, 2H-SiC; 4H-SiC; 6H-SiC, etc.).<br>*&gt;&gt;&gt;* See <a href="http://www.matprop.ru/SiC_bandstr" target="_blank">http://www.matprop.ru/SiC_bandstr</a><br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* So, check carefully your input/output files.<br>
*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* Bests,<br>*&gt;&gt;&gt;* Eyvaz.<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* -------------------------------------------------------------------<br>*&gt;&gt;&gt;* Prof. Eyvaz Isaev,<br>*&gt;&gt;&gt;* Department of Physics, Chemistry, and Biology (IFM), Linkoping<br>
*&gt;&gt;&gt;* University, Sweden Theoretical Physics Department, Moscow State Institute<br>*&gt;&gt;&gt;* of Steel&amp;  Alloys, Russia, isaev at <a href="http://ifm.liu.se/" target="_blank">ifm.liu.se</a><br>&lt;<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>&gt;, eyvaz_isaev at<br>
<a href="http://yahoo.com/" target="_blank">yahoo.com</a> &lt;<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>&gt;<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;*<br>
*&gt;&gt;&gt;* From: pari shok&lt;parishok at <a href="http://gmail.com/" target="_blank">gmail.com</a><br>&lt;<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>&gt;&gt;<br>
*&gt;&gt;&gt;* To: pw_forum at <a href="http://pwscf.org/" target="_blank">pwscf.org</a><br>&lt;<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>&gt;<br>
*&gt;&gt;&gt;* Sent: Sun, April 3, 2011 1:34:51 AM<br>*&gt;&gt;&gt;* Subject: [Pw_forum] shift of energy<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* Dear Paolo,<br>*&gt;&gt;&gt;* DOS of 72-atom SiC shows a shift of energy with respect to 8-atom SiC.The<br>
*&gt;&gt;&gt;* bandgap of 8-atoms is around 10-12 eV, but that of 72-atom is around 5-7<br>*&gt;&gt;&gt;* eV. Would you please help me to understand this shift.<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* Thanks again.<br>*&gt;&gt;&gt;* P Shok<br>
*&gt;&gt;&gt;* UMD<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* On Apr 1, 2011, at 21:36 , pari shok wrote:<br>*&gt;&gt;&gt;&gt;* The DOS of 72-atom SiC (supercell) shows a shift of energy.<br>*&gt;&gt;&gt;* a shift with respect to what?<br>
*&gt;&gt;&gt;* ---<br>*&gt;&gt;&gt;* Paolo Giannozzi, Dept of Chemistry&amp;Physics&amp;Environment,<br>*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>*&gt;&gt;&gt;* Phone +39-0432-558216, fax +39-0432-558222<br>
*&gt;&gt;&gt;*<br>*&gt;&gt;&gt;* _______________________________________________<br>*&gt;&gt;&gt;* Pw_forum mailing list<br>*&gt;&gt;&gt;* Pw_forum at <a href="http://pwscf.org/" target="_blank">pwscf.org</a> &lt;<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>&gt;<br>
*&gt;&gt;&gt;* <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>*&gt;&gt;* ? Gabriele Sclauzero, EPFL SB ITP CSEA<br>*&gt;&gt;*     PH H2 462, Station 3, CH-1015 Lausanne<br>
*&gt;*<br>*&gt;*<br>*<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20110404/6c0c645c/attachment-0001.htm" target="_blank">http://www.democritos.it/pipermail/pw_forum/attachments/20110404/6c0c645c/attachment-0001.htm</a><br>
<br>------------------------------<br><br>Message: 3<br>Date: Mon, 4 Apr 2011 11:04:51 -0400<br>From: Duy Le &lt;<a href="mailto:ttduyle@gmail.com">ttduyle@gmail.com</a>&gt;<br>Subject: Re: [Pw_forum] Shift of Energy<br>To: PWSCF Forum &lt;<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>&gt;<br>
Cc: pari shok &lt;<a href="mailto:parishok@gmail.com">parishok@gmail.com</a>&gt;<br>Message-ID: &lt;BANLkTi=<a href="mailto:keZQ5tpRzvHs7t_d4x7V7Jrw7Qw@mail.gmail.com">keZQ5tpRzvHs7t_d4x7V7Jrw7Qw@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=ISO-8859-1<br><br>On Mon, Apr 4, 2011 at 10:32 AM, pari shok &lt;<a href="mailto:parishok@gmail.com">parishok@gmail.com</a>&gt; wrote:<br>&gt; Dear All,<br>&gt; Thank you very much for your help.<br>
&gt; Please let me ask my final question w.r.t DOS shift of energy of SiC<br>&gt; supercell.<br>&gt; The total energy of the 8-atom cell of SiC is: -76.276 Ry, while that of<br>&gt; 72-atom SiC is -694 Ry. Neglecting almost -7 Ry, the total energy is scaling<br>
&gt; accordingly.<br>If you are doing correctly, the results must be identical. So make<br>sure to check the structures, and the convergence of total energy vs.<br>cutoff energy and k-point sampling.<br>&gt; However, I still see the shift of energy in DOS diagram of two slabs (the<br>
&gt; max and min of bandgap shifts, although the gap is almost 2.1 eV for both<br>&gt; 4H-SiC slabs). I was wondering whether this shift of energy is inevitable or<br>&gt; it is as the result of my input file.<br>Slab? I though you are performing bulk calculation. Anyway, the shift<br>
of energies&#39; levels is normal. You should set the reference as center<br>of the gap (set it to be zero) then compare the DOS or band<br>structures.<br>By the way, if your bigger slab experiences the reconstruction, the<br>
two slabs (8 atoms and 72 atoms) won&#39;t be identical =&gt; you should not<br>get the same results in this case. But 7Ry seems to be huge.<br>&gt; I kept k points (w.r.t. crystal), and cut off energy the same for both<br>
&gt; structures. Should I scale the cut off energy as well. If positive, is there<br>Okie, but you should scale the k-point sampling accordingly because<br>the BZ of 72 atoms cell is smaller. But this won&#39;t harm the results if<br>
total energy is already converged vs. k-point sampling.<br>&gt; any relation between scaling the cell and cut off energy.<br>No. Cut off energy depends on only pseudo potential.<br>&gt; As you know, I need to find a reference for my further calculations.<br>
&gt; I really really appreciate your help.<br>&gt; Yours<br>&gt; P Shok<br>&gt;<br><br>--------------------------------------------------<br>Duy Le<br>PhD Student<br>Department of Physics<br>University of Central Florida.<br>
<br>&quot;Men don&#39;t need hand to do things&quot;<br>&gt; Dear  All,<br>&gt;      let me clarify a little point.<br>&gt;      The potential (hence the eigenvalues) of  any calculation in<br>&gt; periodic boundary conditions, not just a plane-wave one , is defined up<br>
&gt;<br>&gt; to  an additive constant... In a finite system a reeference can be fixed<br>&gt; by requiring that the potential goes to zero at infinite distance but in<br>&gt; a periodic system this cannot be defined as there is nothing like a<br>
&gt;<br>&gt; point at infinite distance ...<br>&gt;     In most codes (and in pw.x) this arbitrary constant is fixed by<br>&gt; setting the average electrostatic potential to zero. If the supercell is<br>&gt; properly constructed, with the properly scaled dimensions, the same<br>
&gt;<br>&gt; cutoff, equivalent k-points etc.. the two calculations should be<br>&gt; equivalent and the eigenvalues should be exactly mapped, the average<br>&gt; potential should be the same, the total energy should be an exact<br>
&gt; multiple of the one of the original cell.<br>&gt;<br>&gt;<br>&gt; If this does not happen one is NOT doing the supercell calculation properly.<br>&gt;<br>&gt; So if when doing a supercell you do not get the scaled result check your<br>
&gt; system definition.<br>&gt;<br>&gt; stefano<br>&gt;<br>&gt;<br>&gt; On 04/04/2011 11:46 AM, Giuseppe Mattioli wrote:<br>&gt;<br>&gt;&gt; Dear all<br>&gt;&gt; I suspect that a a little misunderstanding has been going on here...<br>
&gt;&gt; I try to reformulate the question: in the case of a 8-atoms cell the<br>&gt;&gt; valence<br>&gt;&gt; band maximum (VBM) has been found around 10 eV and the conduction band<br>&gt;<br>&gt;&gt; minimum (CBM) around 12 eV with a 2 eV band gap. So there is a shift in<br>
&gt;&gt; the<br>&gt;&gt; position of VBM and CBM when calculated in a 72-atoms supercell, but the<br>&gt;&gt; band<br>&gt;&gt; gap value is the same. Then, the answer to P Shok could be: plane wave<br>&gt;<br>&gt;&gt; pseudopotential calculations do not ensure an universal alignment of<br>
&gt;&gt; eigenvalues. You should use some &quot;internal&quot; reference, like the 1s<br>&gt;&gt; eigenvalue<br>&gt;&gt; of an He atom, in your cell; or trust that the two VBM values are aligned<br>&gt;&gt; &quot;de<br>&gt;<br>
&gt;&gt; facto&quot;.<br>&gt;&gt; HTH<br>&gt;&gt;<br>&gt;&gt; Giuseppe<br>&gt;&gt;<br>&gt;&gt; On Monday 04 April 2011 09:39:18 Gabriele Sclauzero wrote:<br>&gt;&gt;&gt; Yes,<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;     I perfectly agree with Eyvaz! Moreover the experimental bandgap is<br>
&gt;<br>&gt;&gt;&gt; underestimated within LDA or GGA, and for 4H-SiC you would get about<br>&gt;&gt;&gt; 2.2-2.3 eV with LDA.<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;     How do you measure the bandgap? Beware that in the fundamental cell<br>
&gt;&gt;&gt; the<br>&gt;<br>&gt;&gt;&gt; bandgap is not a direct one, but it might become such when you use a<br>&gt;&gt;&gt; supercell (because of the refolding of the Brillouin zones).<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>
&gt;&gt;&gt; HTH<br>&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; GS<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Il giorno 03/apr/2011, alle ore 10.59, Eyvaz Isaev ha scritto:<br>&gt;&gt;&gt;&gt; Hi,<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt;&gt; The bandgap of 8-atoms is around 10-12 eV,<br>
&gt;<br>&gt;&gt;&gt;&gt;&gt; but that of 72-atom is around 5-7 eV.<br>&gt;&gt;&gt;&gt; This is unbelievable, as SiC is an indirect band gap semiconductor and<br>&gt;&gt;&gt;&gt; the bandgap is around (2.5 - 3)eV depending on crystal modification<br>
&gt;<br>&gt;&gt;&gt;&gt; (3C-SiC, 2H-SiC; 4H-SiC; 6H-SiC, etc.).<br>&gt;&gt;&gt;&gt; See <a href="http://www.matprop.ru/SiC_bandstr" target="_blank">http://www.matprop.ru/SiC_bandstr</a><br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; So, check carefully your input/output files.<br>
&gt;<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; Bests,<br>&gt;&gt;&gt;&gt; Eyvaz.<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; -------------------------------------------------------------------<br>&gt;&gt;&gt;&gt; Prof. Eyvaz Isaev,<br>
&gt;<br>&gt;&gt;&gt;&gt; Department of Physics, Chemistry, and Biology (IFM), Linkoping<br>&gt;&gt;&gt;&gt; University, Sweden Theoretical Physics Department, Moscow State<br>&gt;&gt;&gt;&gt; Institute<br>&gt;&gt;&gt;&gt; of Steel&amp;  Alloys, Russia, isaev at <a href="http://ifm.liu.se/" target="_blank">ifm.liu.se</a>, eyvaz_isaev at <a href="http://yahoo.com/" target="_blank">yahoo.com</a><br>
&gt;<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; From: pari shok&lt;parishok at <a href="http://gmail.com/" target="_blank">gmail.com</a>&gt;<br>&gt;&gt;&gt;&gt; To: pw_forum at <a href="http://pwscf.org/" target="_blank">pwscf.org</a><br>
&gt;<br>&gt;&gt;&gt;&gt; Sent: Sun, April 3, 2011 1:34:51 AM<br>&gt;&gt;&gt;&gt; Subject: [Pw_forum] shift of energy<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; Dear Paolo,<br>&gt;&gt;&gt;&gt; DOS of 72-atom SiC shows a shift of energy with respect to 8-atom<br>
&gt;&gt;&gt;&gt; SiC.The<br>&gt;<br>&gt;&gt;&gt;&gt; bandgap of 8-atoms is around 10-12 eV, but that of 72-atom is around 5-7<br>&gt;&gt;&gt;&gt; eV. Would you please help me to understand this shift.<br>&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; Thanks again.<br>&gt;<br>&gt;&gt;&gt;&gt; P Shok<br>&gt;&gt;&gt;&gt; UMD<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; On Apr 1, 2011, at 21:36 , pari shok wrote:<br>&gt;&gt;&gt;&gt;&gt; The DOS of 72-atom SiC (supercell) shows a shift of energy.<br>
&gt;<br>&gt;&gt;&gt;&gt; a shift with respect to what?<br>&gt;&gt;&gt;&gt; ---<br>&gt;&gt;&gt;&gt; Paolo Giannozzi, Dept of Chemistry&amp;Physics&amp;Environment,<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
&gt;<br>&gt;&gt;&gt;&gt; Phone +39-0432-558216, fax +39-0432-558222<br>&gt;&gt;&gt;&gt;<br>&gt;&gt;&gt;&gt; _______________________________________________<br>&gt;&gt;&gt;&gt; Pw_forum mailing list<br>&gt;&gt;&gt;&gt; Pw_forum at <a href="http://pwscf.org/" target="_blank">pwscf.org</a><br>
&gt;<br>&gt;&gt;&gt;&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>&gt;&gt;&gt; ? Gabriele Sclauzero, EPFL SB ITP CSEA<br>&gt;&gt;&gt;     PH H2 462, Station 3, CH-1015 Lausanne<br>
&gt;<br>&gt;&gt;<br>&gt;&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>&gt;<br><br><br>------------------------------<br><br>Message: 4<br>Date: Tue, 5 Apr 2011 00:14:21 +0800 (CST)<br>From: ?S &lt;<a href="mailto:flux_ray12@163.com">flux_ray12@163.com</a>&gt;<br>Subject: Re: [Pw_forum] 4.3 compile error.<br>
To: &quot;PWSCF Forum&quot; &lt;<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>&gt;<br>Message-ID: &lt;<a href="mailto:174f5d1.7e22.12f214b1d61.Coremail.flux_ray12@163.com">174f5d1.7e22.12f214b1d61.Coremail.flux_ray12@163.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;gbk&quot;<br><br>Hi, Vit.<br>I got the same problem like you.<br><br>You can open the file TDDFPT/src/lr_main.f90 and find the line 252. It should be written as:<br>IF ( lr_io_level &gt; 0 .and. (mod(LR_iteration,restart_step)==0 .OR. LR_iteration==itermax .OR. LR_iteration==1) ) CALL lr_write_restart()<br>
Then, you can add a &#39;&amp;&#39; after &#39;==1) )&#39; and put &#39;CALL lr_write_restart()&#39; to the next line.<br>Do not forget add enough space before &#39;CALL lr_write_restart()&#39;, otherwise it would become annotate instead of code.<br>
The modified one could like:<br>IF ( lr_io_level &gt; 0 .and. (mod(LR_iteration,restart_step)==0 .OR. LR_iteration==itermax .OR. LR_iteration==1) )  &amp;<br>          CALL lr_write_restart()<br><br><br>At 2011-04-03 19:41:28?Vit &lt;<a href="mailto:vitruss@gmail.com">vitruss@gmail.com</a>&gt; wrote:<br>
<br>&gt;Dear QE users!<br>&gt;I&#39;m having troubles compiling new release. I&#39;m trying to use ACML where<br>&gt;apllicable, but I&#39;m getting the same error even with distro (debian 6.0.1)<br>&gt;blas and lapack and internal fftw.<br>
&gt;Could you please help me resolving this issue?<br>&gt;<br>&gt;mpif90 -g -o turbo_lanczos.x \<br>&gt;                lr_variables.o lr_charg_resp.o bcast_lr_input.o lr_readin.o<br>&gt;lr_alloc_init.o lr_calc_dens.o lr_dot.o lr_dealloc.o lr_ortho.o lr_read_wf.o<br>
&gt;lr_normalise.o lr_lanczos.o lr_apply_liouvillian.o lr_main.o lr_dv_setup.o<br>&gt;lr_setup_dgc.o lr_solve_e.o lr_dvpsi_e.o lr_ch_psi_all.o lr_cgsolve_all.o<br>&gt;lr_h_psiq.o lr_sm1_psi.o stop_lr.o lr_read_d0psi.o lr_restart.o<br>
&gt;lr_write_restart.o print_clock_lr.o sd0psi.o lr_set_boxes_density.o<br>&gt;lr_init_nfo.o  ../../PH/libph.a ../../PW/libpw.a ../../Modules/libqemod.a<br>&gt;../../flib/ptools.a ../../flib/flib.a ../../clib/clib.a<br>&gt;../../iotk/src/libiotk.a   -L/opt/acml4.4.0/gfortran64/lib -lacml -<br>
&gt;L/opt/acml4.4.0/gfortran64/lib -lacml -L/opt/acml4.4.0/gfortran64/lib -lacml<br>&gt;lr_main.o: In function `lr_main&#39;:<br>&gt;/opt/qe/4.3/TDDFPT/src/lr_main.f90:252: undefined reference to `lr_wri_&#39;<br>&gt;collect2: ld returned 1 exit status<br>
&gt;make[2]: *** [turbo_lanczos.x] Error 1<br>&gt;make[2]: Leaving directory `/opt/qe/4.3/TDDFPT/src&#39;<br>&gt;make[2]: Entering directory `/opt/qe/4.3/TDDFPT/tools&#39;<br>&gt;test -n &quot;&quot; &amp;&amp; ( cd ../.. ; make -w  || exit 1) || :<br>
&gt;mpif90 -O3 -g -x f95-cpp-input -D__GFORTRAN -D__STD_F95 -D__ACML -D__MPI -<br>&gt;D__PARA -I../../include -I../../iotk/src -I../../Modules -I. -c<br>&gt;tddfpt_calculate_spectrum.f90<br>&gt;mpif90 -g -o turbo_spectrum.x \<br>
&gt;        tddfpt_calculate_spectrum.o ../../Modules/libqemod.a ../../PW/libpw.a<br>&gt;../../flib/ptools.a ../../flib/flib.a ../../clib/clib.a<br>&gt;../../iotk/src/libiotk.a   -L/opt/acml4.4.0/gfortran64/lib -lacml -<br>
&gt;L/opt/acml4.4.0/gfortran64/lib -lacml -L/opt/acml4.4.0/gfortran64/lib -lacml<br>&gt;( cd ../../bin ; ln -fs ../TDDFPT/tools/turbo_spectrum.x . )<br>&gt;if [ -d ../bin ] ; then  ( cd ../bin ; ln -fs ../tools/turbo_spectrum.x . );<br>
&gt;fi<br>&gt;make[2]: Leaving directory `/opt/qe/4.3/TDDFPT/tools&#39;<br>&gt;make[1]: Leaving directory `/opt/qe/4.3/TDDFPT&#39;<br>&gt;<br>&gt;With best regards,<br>&gt;Koroteev Victor.<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>
-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20110405/a5f0eb48/attachment.htm" target="_blank">http://www.democritos.it/pipermail/pw_forum/attachments/20110405/a5f0eb48/attachment.htm</a><br>
<br></blockquote>
<div>  Hi dear users.</div>
<div>  is there anybody who&#39;s know about Teter pseudopotential.I used it for Sr.Sr has 38 electrons and valence electrons should be 4s2 4p6 5s2 4d0 in ground state but in this pseudopotential zion=10 and the valence electrons are 4s2 4p6 5s0 4d0!whate&#39;s the maning of it?can i use it for ground state alculation?</div>

<div>Ok! this pseudopotential is used in abinit but i can&#39;t ask my question of them.</div>
<div>thank&#39;s alot.</div>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">------------------------------<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 46, Issue 10<br>
****************************************<br></blockquote></div><br>