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