<div> </div>  <div>Dear Stefano,</div>  <div>Thank you, I will study the relevant part for pert. potential., do the suggested tests and modifications,</div>  <div>and send the result back to close the case.</div>  <div>But a fast point that I wanted to note that is,</div>  <div>InN turns out to be a metal according to my calculations,</div>  <div>without several corrections used, like GW or SIC.</div>  <div>That was the reason of the k-mesh.</div>  <div>Thank you for your help.</div>  <div><BR>Emine Kucukbenli</div>  <div><BR><B><I>Stefano de Gironcoli <degironc@sissa.it></I></B> wrote:</div>  <BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid"><BR>If I remember correctly there used to be a term in the phonon <BR>calculation in which the perturbing potential was written as DeltaV(r) = <BR>exp(+iqr)*\DeltaVperiodic(r) and \DeltaVperiodic was expanded up to<BR>the charge density cutoff ... this term was
 therefore slightly different<BR>for two q wavevectors that would otherwise be equivalent. This never<BR>gave any noticeable difference in real calculations and I would be <BR>surprised if this would be the cause of your results but a way to check<BR>this is to increase the charge density cutoff and see what happens as<BR>in this case the problem should disappear when the potential description<BR>become complete.<BR><BR>From your input it looks to me that you are using the default value for <BR>the selfconsistent threshold (conv_thr=1.d-6 in the &electrons namelist).<BR>This is usually fine for total energy which is variational, but already <BR>for structural relaxation it is often reduced so as to get more accurate <BR>forces.<BR><BR>Phonon calculations are even more sensitive to the quality of <BR>self-consistency of your unperturbed system. Try to reduce conv_thr to <BR>1.d-8 in the scf (and nscf) steps and see if this improves things.<BR><BR>InN is a (wide band gap)
 semiconductor: there is no real need to use a <BR>metallic setup which, by the way, is incompatible with the calculation<BR>of dielectric constant and effective charges, that you will probably need<BR>to get correct phonons frequencies at gamma.<BR>I would expect that 2 or 10 Chadi-Cohen points should be enough for this <BR>system (automatic grids 2 2 2 1 1 1 and 4 4 4 1 1 1, respectively).<BR><BR>hope this helps,<BR><BR>stefano<BR><BR>On Sun, 10 Dec 2006, e kb wrote:<BR><BR>> Dear Stefano , QE users,<BR>> Thank you for your interest,<BR>><BR>> Here are the scf, nscf and ph.x input files respectively, for X(0,0,1).<BR>> For X(0,1,1) the scf file is the same, in nscf input I just changed the phonon namelist to (0.00 1.00 1.00) so not to make the mail longer and unreadable, I am skipping it and just adding ph.X input for X(0,1,1).<BR>><BR>> I should note that this was a part of test calculations for the new pseudo for Indium that I recently generated
 using the 'atomic' code. The cut off energy that I used for these calculations (34 Ryd) is chosen to be enough to represent correct lattice constants, but still the pseudo and the value of cut off energy needs further examination. That is why I would not be surprised to see any discrepiencies in the frequencies. But again, I could not understand the differences that occured for the same point.<BR>> With a verbose output I checked the phonon polarization vectors,<BR>> As far as I understand, the calculation is done for the same modes for the two points.<BR>><BR>> I do not want to bother you with the whole output files,<BR>> if this is a very simple problem that an experienced/careful user can notice from the inputs,<BR>> otherwise I can attach the output files.<BR>><BR>> Buona domenica,<BR>> Emine Kucukbenli<BR>> ------------<BR>> &control<BR>> calculation = 'scf'<BR>> restart_mode='from_scratch',<BR>> prefix='34_vib',<BR>>
 tprnfor = .true.<BR>> pseudo_dir = '$PSEUDO_DIR' ,<BR>> etot_conv_thr = 1.0D-5<BR>> /<BR>> &system<BR>> ibrav= 2,<BR>> celldm(1)=9.593140000<BR>> nat= 2,<BR>> ntyp= 2,<BR>> ecutwfc =34.0,<BR>> ecutrho =272.0,<BR>> occupations = 'smearing'<BR>> smearing = 'gaussian'<BR>> degauss = 0.004,<BR>> /<BR>> &electrons<BR>> mixing_mode = 'plain'<BR>> mixing_beta = 0.7<BR>> /<BR>> ATOMIC_SPECIES<BR>> In 114.818 In_us_2.upf<BR>> N 14.0067 N.pbe-van_ak.UPF<BR>> ATOMIC_POSITIONS { crystal }<BR>> In 0.00000000000 0.00000000000 0.0000000000000<BR>> N 0.25000000000 0.25000000000 0.2500000000000<BR>> K_POINTS { automatic }<BR>> 8 8 8 0 0 0<BR>> ------------<BR>> &control<BR>> calculation = 'phonon'<BR>> restart_mode='from_scratch',<BR>> prefix='34_vib',<BR>> tprnfor = .true.<BR>> pseudo_dir = '$PSEUDO_DIR' ,<BR>> etot_conv_thr = 1.0D-5<BR>> /<BR>> &system<BR>>
 ibrav= 2,<BR>> celldm(1)=9.593140000<BR>> nat= 2,<BR>> ntyp= 2,<BR>> ecutwfc =34.0,<BR>> ecutrho =272.0,<BR>> occupations = 'smearing'<BR>> smearing = 'gaussian'<BR>> degauss = 0.004,<BR>> /<BR>> &electrons<BR>> mixing_mode = 'plain'<BR>> mixing_beta = 0.7<BR>> /<BR>> &phonon<BR>> xqq(1) = .00, xqq(2) = .00, xqq(3) = 1.00<BR>> /<BR>> ATOMIC_SPECIES<BR>> In 114.818 In_us_2.upf<BR>> N 14.0067 N.pbe-van_ak.UPF<BR>> ATOMIC_POSITIONS { crystal }<BR>> In 0.00000000000 0.00000000000 0.0000000000000<BR>> N 0.25000000000 0.25000000000 0.2500000000000<BR>> K_POINTS { automatic }<BR>> 8 8 8 0 0 0<BR>> -------------------------<BR>><BR>> &inputph<BR>> tr2_ph=1.0d-14,<BR>> amass(1)=114.818,<BR>> amass(2)=14.0067,<BR>> prefix='34_vib',<BR>> fildyn='34_vib.dyn',<BR>> /<BR>> 0.00 0.00 1.00<BR>> --------------------------------<BR>><BR>><BR>> &inputph<BR>>
 tr2_ph=1.0d-14,<BR>> amass(1)=114.818,<BR>> amass(2)=14.0067,<BR>> prefix='34_vib',<BR>> fildyn='34_vib.dyn',<BR>> /<BR>> 0.00 1.00 1.00<BR>> ----------------------------------<BR>> X(0,0,1)<BR>> omega( 1) = 3.163928 [THz] = 105.537987 [cm-1]<BR>> omega( 2) = 3.163928 [THz] = 105.537987 [cm-1]<BR>> omega( 3) = 6.030058 [THz] = 201.142433 [cm-1]<BR>> omega( 4) = 14.775784 [THz] = 492.870383 [cm-1]<BR>> omega( 5) = 14.775784 [THz] = 492.870383 [cm-1]<BR>> omega( 6) = 16.372030 [THz] = 546.115769 [cm-1]<BR>> X(0,1,1)<BR>> omega( 1) = 2.153190 [THz] = 71.823161 [cm-1]<BR>> omega( 2) = 2.153190 [THz] = 71.823161 [cm-1]<BR>> omega( 3) = 6.239891 [THz] = 208.141736 [cm-1]<BR>> omega( 4) = 14.879198 [THz] = 496.319913 [cm-1]<BR>> omega( 5) = 14.879198 [THz] = 496.319913 [cm-1]<BR>> omega( 6) = 16.354833 [THz] = 545.542147 [cm-1]<BR>> ----------------------------------<BR>><BR>> Stefano Baroni
 <BARONI@SISSA.IT>wrote:<BR>> Dear Emin, dear all: all my apologies!<BR>><BR>> Nicola was of course right, and I was WRONG! (011) IS indeed equivalent to (001), because the latter is equivalent to (-1,0,0) and the difference between (011) and (-100) is (111) which IS an FCC reciprocal-lattice vector. The coordinates of the K point are (3/4,3/4,0). With reference to the figure mentioned in the link below, K is the intersection among the planes x+y+z=3/2, z=0, and x=y, hence the above coordinates.<BR>><BR>> To come to your original question. Are the symmetries/degeneracies of the results what you expect? Can you make the input/output of your calculation accessible to us?<BR>><BR>><BR>> Stefano<BR>><BR>> On Dec 10, 2006, at 2:46 PM, e kb wrote:<BR>><BR>> Thank you for your answer Stefano,<BR>> Of course this explains the different frequencies,<BR>> and that also means that I misunderstand the kpoints given in example 5,<BR>> about
 band structure for an fcc,<BR>> and Dr. Marzari's explanation to another QE user,<BR>> who was confused like me.<BR>> http://www.democritos.it/pipermail/pw_forum/2006-June/004365.html<BR>><BR>> Thank you again,<BR>> Emine Kucukbenli<BR>><BR>> Stefano Baroni <BARONI@SISSA.IT>wrote:<BR>> Emine: (001) and (011) are *NOT* two "X" points in the FCC Brillouin zone. The difference between the second and the first (010) is *NOT* a reciprocal-space vector. If it was, X would also be, which you will agree it is not. As a matter of fact, (011) has its own nickname (it is called "K"). No surprise, therefore, that frequencies @K are different from frequencies @X. Are these arguments clear and convincing enough? Stefano<BR>><BR>> On Dec 9, 2006, at 11:10 PM, e kb wrote:<BR>><BR>> Dear QE users,<BR>> I am performing phonon calculations for a zinc blende structure,<BR>> I calculated the frequencies at high symmetry points, evaluating a single
 point at each time.<BR>> So I followed the procedure of scf, nscf and ph.x calculation one by one for L(0.5 0.5 0.5) - Gamma - X(0,0,1) - X(0,1,1) points.<BR>> In these calculations every parameter (else than q-vector in &phonon list in nscf calculation and phonon card for ph.x input of course) is kept fixed.<BR>><BR>> Although I expect the frequencies at two X points - (0,0,1) and (0,1,1) - to be the same,<BR>> the difference in TA modes' frequencies are around 30 cm-1,<BR>> in the other ones, around 5 cm-1.<BR>><BR>> I would like to attribute this to unconverged lattice parameters, cut off energies or thresholds etc. but what I do not understand is why the program treats these two points differently, although they should be equivalent. So I expect to see the same value for frequency at these points, even if it is wrong, highly negative, or sth else that it should not be.<BR>><BR>> I will appreciate any answers.<BR>> Thank
 you.<BR>><BR>> Emine Kucukbenli<BR>><BR>><BR>><BR>><BR>><BR>> ---------------------------------<BR>> Everyone is raving about the all-new Yahoo! Mail beta.<BR>><BR>> ---<BR>> Stefano Baroni - SISSA & DEMOCRITOS National Simulation Center - Trieste<BR>> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)<BR>><BR>><BR>> Please, if possible, don't send me MS Word or PowerPoint attachments<BR>> Why? See: http://www.gnu.org/philosophy/no-word-attachments.html<BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>> ---------------------------------<BR>> Check out the all-new Yahoo! Mail beta - Fire up a more powerful email and get things done faster.<BR>><BR>> ---<BR>> Stefano Baroni - SISSA & DEMOCRITOS National Simulation Center - Trieste<BR>> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)<BR>><BR>><BR>> Please, if possible, don't send
 me MS Word or PowerPoint attachments<BR>> Why? See: http://www.gnu.org/philosophy/no-word-attachments.html<BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>><BR>> ---------------------------------<BR>> Have a burning question? Go to Yahoo! Answers and get answers from real people who know.<BR>_______________________________________________<BR>Pw_forum mailing list<BR>Pw_forum@pwscf.org<BR>http://www.democritos.it/mailman/listinfo/pw_forum<BR></BLOCKQUOTE><BR><p>
<hr size=1>Cheap Talk? <a href="http://us.rd.yahoo.com/mail_us/taglines/postman8/*http://us.rd.yahoo.com/evt=39663/*http://voice.yahoo.com">Check out</a> Yahoo! Messenger's low PC-to-Phone call rates.