Dear Pwscf users,<br>Here I am sending you the input file of Li2S (large band gap semi-conductor) for scf and phonon calculations. I have some queries regarding this<br><br>1. For electronic calculation i am using 8x8x8 k-points which results in only 60 points. Is it enough to do correct electronic calculations?<br>
With these 60 points only, I calculated the bulk mobulus and the equilibrium lattice constants, which matches well with the experiment and other theoritical results.<br><br>2. For phonon calculation, I didnt get the LO-TO splitting. Can you please tell me what more parameters i will use to get the splitting. <br>
I set epsil .true., but it gives the error&nbsp; <br><br>&nbsp;&nbsp;&nbsp;&nbsp; from phq_readin : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; no elec. field with metals<br><br>Any suggestion is highly apperciable <br><br><a href="http://li2s.scf.in">li2s.scf.in</a><br>
<br>&amp;CONTROL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = &#39;scf&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = &#39;from_scratch&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = &#39;./OUT&#39;,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;/root/bipul/espresso-4.0.1/pseudo/Li2S&#39;,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = &#39;Li2S&#39;,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0d-5 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; forc_conv_thr = 1.0d-4 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true. ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true. ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dt = 80<br>
&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) =&nbsp; 10.72700,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 3,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 32,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 600,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = &#39;smearing&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = &#39;methfessel-paxton&#39;,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.01<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.5 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-8<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = &#39;david&#39; ,<br>&nbsp;/<br>&nbsp;&amp;IONS<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;S&nbsp;&nbsp; 32.066&nbsp; S.pw91-van_ak.UPF<br>&nbsp;Li&nbsp; 6.914&nbsp;&nbsp; Li.pw91-s-van_ak.UPF<br>ATOMIC_POSITIONS (crystal)<br>S&nbsp; 0.0 0.0 0.0<br>Li 0.25 0.25 0.25<br>
Li 0.75 0.75 0.75<br>K_POINTS automatic<br>&nbsp; 8 8 8 1 1 1<br><br><a href="http://li2s.ph.in">li2s.ph.in</a><br><br>phonons of Li2S<br>&nbsp;&amp;inputph<br>&nbsp;! recover =.true.,<br>&nbsp; tr2_ph=1.0d-12,<br>&nbsp; alpha_mix(1)=0.6,<br>&nbsp; prefix=&#39;Li2S&#39;,<br>
&nbsp;! fildvscf=&#39;li2sdv&#39;,<br>&nbsp; amass(1)=32.066,<br>&nbsp; amass(2)=6.941,<br>&nbsp; outdir=&#39;./OUT&#39;,<br>&nbsp; fildyn=&#39;li2s.dyn&#39;,<br>&nbsp;! elph=.true.,<br>&nbsp;! trans=.true.,<br>&nbsp; ldisp=.true.<br>&nbsp; nq1=4, nq2=4, nq3=4<br>&nbsp;/<br>
<br clear="all"><br>-- <br>Premlata Pandit<br>Ph.d. student,<br>Barkatullah university,<br>Bhopal 462026,<br>MP, India <br>