<div dir="ltr"><span style="font-size:13px">Hi Jiqiang,</span><div style="font-size:13px"><br></div><div style="font-size:13px">Thank you very much for your help! </div><div style="font-size:13px"><br></div><div style="font-size:13px">I tried using only gamma point, but I got an error message when I calculated phonon saying "can't start from gamma tricks...". Besides, I need to calculate the raman frequencies after phonon calculation, so I think it's better if I use more k points.</div><div style="font-size:13px"><br></div><div style="font-size:13px">Have you tried Raman calculation? Could you give me more advice?</div><div style="font-size:13px"><br></div><div style="font-size:13px">Best,</div><div style="font-size:13px"><br></div><div style="font-size:13px">Xiaolin</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 28, 2015 at 10:21 PM, Xiaolin Xu <span dir="ltr"><<a href="mailto:xxu2@caltech.edu" target="_blank">xxu2@caltech.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Jiqiang,<div><br></div><div>Thank you very much for your help! </div><div><br></div><div>I tried using only gamma point, but I got an error message when I calculated phonon saying "can't start from gamma tricks...". Besides, I need to calculate the raman frequencies after phonon calculation, so I think it's better if I use more k points.</div><div><br></div><div>Have you tried Raman calculation? Could you give me more advice?</div><div><br></div><div>Best,</div><div><br></div><div>Xiaolin</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 28, 2015 at 6:20 PM, jqli14 <span dir="ltr"><<a href="mailto:jqli14@fudan.edu.cn" target="_blank">jqli14@fudan.edu.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">(1)employ k_points gamma instead of k_pionts automatic 1 1 1 0 0 0, although both give the same result, but the later does not call the gamma point specific algorithm to speed the calculation.<br>
(2)if it can not converge within handreds of electronic steps, reduce the value of mixing_beta or change the mixing_mode.<br>
<br>
Jiqiang Li<br>
<div><div><br>
<br>
在 2015-10-29 07:57:58，"Xiaolin Xu" <<a href="mailto:xiaolin0314@gmail.com" target="_blank">xiaolin0314@gmail.com</a>> 写道：<br>
<br>
>Dear all,<br>
><br>
><br>
>I was doing phonon dispersion calculations for lithium iron phosphate. However, the calculation speed is really slow and it does not converge. Is this because my system is too big? There are 28 atoms in the unit cell. My input file is like this:<br>
>(1)<br>
><br>
><br>
><br>
>&CONTROL<br>
><br>
>calculation = "scf",<br>
><br>
>pseudo_dir = "/home/xiaolin/pseudo",<br>
><br>
>outdir = "/home/xiaolin/test2",<br>
><br>
>tstress=.true.,<br>
><br>
>tprnfor=.true.,<br>
><br>
>nstep=100,<br>
><br>
>/<br>
><br>
>&SYSTEM<br>
><br>
>ibrav = 8,<br>
><br>
>celldm(1) = 19.749526<br>
><br>
>celldm(2) = 0.582432,<br>
><br>
>celldm(3) = 0.454693,<br>
><br>
>nat = 28,<br>
><br>
>ntyp = 4,<br>
><br>
>ecutwfc = 80.D0,<br>
><br>
>/<br>
><br>
>&ELECTRONS<br>
><br>
>mixing_beta = 0.7,<br>
><br>
>/<br>
><br>
>&IONS<br>
><br>
>/<br>
><br>
>ATOMIC_SPECIES<br>
><br>
>Li 6.941 Li.pbe-s-hgh.UPF<br>
><br>
>Fe 55.845 Fe.pbe-sp-hgh.UPF<br>
><br>
>P 30.97376 P.pbe-hgh.UPF<br>
><br>
>O 15.9994 O.pbe-hgh.UPF<br>
><br>
>ATOMIC_POSITIONS (alat)<br>
><br>
>Li 0 0 0<br>
><br>
>...<br>
><br>
>K_POINTS (automatic)<br>
><br>
>1 1 1 0 0 0<br>
><br>
><br>
><br>
><br>
>(2)<br>
><br>
><br>
><br>
><br>
>Normal modes for lifepo4<br>
><br>
>&inputph<br>
><br>
>tr2_ph=1.5d-5,<br>
><br>
>amass(1)=6.941,<br>
><br>
>amass(2)=55.845,<br>
><br>
>amass(3)=30.97376,<br>
><br>
>amass(4)=15.9994,<br>
><br>
>outdir = '/home/xiaolin/test2',<br>
><br>
>trans=.true.,<br>
><br>
>asr=.true.,<br>
><br>
>nogg=.true.,<br>
><br>
>fildyn='dmat.lifepo4'<br>
><br>
>/<br>
><br>
>0.0 0.0 0.0<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
>Thank you in advance!<br>
><br>
><br>
><br>
><br>
>Lin<br>
<br>
<br>
<br>
<br>
<br>
<br>
</div></div>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></blockquote></div><br></div>
</div></div></blockquote></div><br></div>