&gt;Dear all,<br>&gt;I have a problem in calculation force constant with q2r.x.<br>&gt;Here are my input files of pw.x ph.x q2r.x. I am so sad. Please tell me the reasons &gt;Thank you so much.<br>+++++++++++++++++++++++++++++++scf.in<br>&amp;CONTROL<br>&nbsp;&nbsp; calculation =&nbsp;&nbsp; 'scf' ,<br>&nbsp;&nbsp; restart_mode = 'from_scratch' ,<br>&nbsp;&nbsp; outdir='/home/user/ecnuphyk/SUN/tmp/tmp-p1/' ,<br>&nbsp;&nbsp; pseudo_dir = '/home/user/ecnuphyk/SUN/pseudo/' ,<br>&nbsp;&nbsp; prefix='alh3' ,<br>&nbsp;&nbsp; tstress = .true. ,<br>&nbsp;&nbsp; tprnfor = .true. ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0E-5&nbsp; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; forc_conv_thr = 1.0D-4 ,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp; ibrav = 0 ,<br>&nbsp;&nbsp; celldm(1) = 1.889725988578923203102406997093&nbsp;&nbsp;
 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat&nbsp; =&nbsp; 16 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp&nbsp; =&nbsp; 2 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin&nbsp; =&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc&nbsp; =&nbsp; 44.0982 ,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr&nbsp; =&nbsp; 1.0d-7&nbsp; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.7
 ,<br>&nbsp;/<br>CELL_PARAMETERS (alat)<br>...<br>ATOMIC_SPECIES<br>...<br>ATOMIC_POSITIONS {crystal}<br>...<br>K_POINTS automatic<br>&nbsp; 4 4 3 0 0&nbsp; 0<br>++++++++++++++++++++++++++++ph.in<br>&amp;inputph<br>&nbsp; tr2_ph=1.0d-12,<br>&nbsp; prefix='alh3',<br>&nbsp; ldisp=.true.,<br>&nbsp; nq1=4, nq2=4, nq3=3<br>&nbsp; epsil=.true. ,<br>&nbsp; amass(1)=26.98,<br>&nbsp; amass(2)=1.0,<br>&nbsp; outdir= '/home/user/ecnuphyk/SUN/tmp/tmp-p1/' ,<br>&nbsp; fildyn='alh3.dyn',<br>&nbsp;/<br>+++++++++++++++++++++++++++++q2r.in<br>&nbsp;&amp;input<br>&nbsp;&nbsp; fildyn='alh3.dyn', zasr='simple', flfrc='alh3444.fc'<br>&nbsp;/<br style="background-color: rgb(255, 0, 0);"><span style="background-color: rgb(255, 0, 0);">++++++++++++++++++++++++++++q2r.in.out&nbsp;&nbsp; (The error )<br>&nbsp; reading grid info from file alh3.dyn0<br>&nbsp; reading force constants from file alh3.dyn1<br>&nbsp;macroscopic fields =&nbsp; T<br>&nbsp;&nbsp; 4.51585&nbsp;&nbsp; 0.00000&nbsp;
 -0.08800<br>&nbsp;&nbsp; 0.00000&nbsp;&nbsp; 4.37104&nbsp;&nbsp; 0.00000<br>&nbsp; -0.08800&nbsp;&nbsp; 0.00000&nbsp;&nbsp; 4.43104<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp;&nbsp; 2.74722&nbsp; -0.27565&nbsp; -0.08418<br>&nbsp;&nbsp; 0.24590&nbsp;&nbsp; 2.66281&nbsp;&nbsp; 0.41340<br>&nbsp; -0.11507&nbsp; -0.40383&nbsp;&nbsp; 2.71996<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<br>&nbsp;&nbsp; 2.74722&nbsp;&nbsp; 0.27565&nbsp; -0.08418<br>&nbsp; -0.32097&nbsp;&nbsp; 2.64842&nbsp; -0.40555<br>&nbsp; -0.11506&nbsp;&nbsp; 0.40383&nbsp;&nbsp; 2.71996<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3<br>&nbsp;&nbsp; 2.74722&nbsp; -0.27565&nbsp; -0.08418<br>&nbsp;&nbsp; 0.24590&nbsp;&nbsp; 2.66281&nbsp;&nbsp; 0.41340<br>&nbsp; -0.11507&nbsp; -0.40383&nbsp;&nbsp; 2.71996<br>&nbsp;
 na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4<br>&nbsp;&nbsp; 2.74722&nbsp;&nbsp; 0.27565&nbsp; -0.08418<br>&nbsp; -0.32097&nbsp;&nbsp; 2.64842&nbsp; -0.40555<br>&nbsp; -0.11507&nbsp;&nbsp; 0.40383&nbsp;&nbsp; 2.71996<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5<br>&nbsp; -0.62424&nbsp;&nbsp; 0.00001&nbsp;&nbsp; 0.05784<br>&nbsp; -0.00008&nbsp; -0.63566&nbsp; -0.00071<br>&nbsp;&nbsp; 0.06224&nbsp; -0.00088&nbsp; -1.42838<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6<br>&nbsp; -1.06619&nbsp; -0.39916&nbsp;&nbsp; 0.01343<br>&nbsp; -0.40231&nbsp; -0.99561&nbsp;&nbsp; 0.02342<br>&nbsp;&nbsp; 0.01086&nbsp;&nbsp; 0.02141&nbsp; -0.62960<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7<br>&nbsp; -1.06619&nbsp;&nbsp; 0.39915&nbsp;&nbsp; 0.01342<br>&nbsp;&nbsp; 0.40230&nbsp; -0.99561&nbsp; -0.02341<br>&nbsp;&nbsp; 0.01086&nbsp;
 -0.02142&nbsp; -0.62960<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8<br>&nbsp; -0.62424&nbsp; -0.00002&nbsp;&nbsp; 0.05784<br>&nbsp;&nbsp; 0.00007&nbsp; -0.63566&nbsp;&nbsp; 0.00072<br>&nbsp;&nbsp; 0.06224&nbsp;&nbsp; 0.00091&nbsp; -1.42838<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 9<br>&nbsp; -1.06613&nbsp;&nbsp; 0.39958&nbsp;&nbsp; 0.01388<br>&nbsp;&nbsp; 0.40245&nbsp; -0.99702&nbsp; -0.02390<br>&nbsp;&nbsp; 0.01136&nbsp; -0.02197&nbsp; -0.62981<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10<br>&nbsp; -1.06613&nbsp; -0.39958&nbsp;&nbsp; 0.01389<br>&nbsp; -0.40245&nbsp; -0.99702&nbsp;&nbsp; 0.02391<br>&nbsp;&nbsp; 0.01136&nbsp;&nbsp; 0.02197&nbsp; -0.62981<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 11<br>&nbsp; -0.62425&nbsp;&nbsp; 0.00002&nbsp;&nbsp; 0.05784<br>&nbsp; -0.00007&nbsp; -0.63567&nbsp;
 -0.00072<br>&nbsp;&nbsp; 0.06224&nbsp; -0.00091&nbsp; -1.42827<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12<br>&nbsp; -1.06613&nbsp;&nbsp; 0.39959&nbsp;&nbsp; 0.01389<br>&nbsp;&nbsp; 0.40248&nbsp; -0.99704&nbsp; -0.02391<br>&nbsp;&nbsp; 0.01135&nbsp; -0.02197&nbsp; -0.62985<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 13<br>&nbsp; -1.06613&nbsp; -0.39959&nbsp;&nbsp; 0.01388<br>&nbsp; -0.40247&nbsp; -0.99704&nbsp;&nbsp; 0.02390<br>&nbsp;&nbsp; 0.01136&nbsp;&nbsp; 0.02198&nbsp; -0.62985<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14<br>&nbsp; -1.06619&nbsp; -0.39916&nbsp;&nbsp; 0.01342<br>&nbsp; -0.40233&nbsp; -0.99564&nbsp;&nbsp; 0.02341<br>&nbsp;&nbsp; 0.01085&nbsp;&nbsp; 0.02143&nbsp; -0.62963<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 15<br>&nbsp; -1.06619&nbsp;&nbsp; 0.39916&nbsp;&nbsp; 0.01343<br>&nbsp;&nbsp; 0.40233&nbsp;
 -0.99563&nbsp; -0.02342<br>&nbsp;&nbsp; 0.01085&nbsp; -0.02142&nbsp; -0.62963<br>&nbsp; na=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 16<br>&nbsp; -0.62425&nbsp; -0.00001&nbsp;&nbsp; 0.05784<br>&nbsp;&nbsp; 0.00008&nbsp; -0.63567&nbsp;&nbsp; 0.00071<br>&nbsp;&nbsp; 0.06223&nbsp;&nbsp; 0.00088&nbsp; -1.42827<br>&nbsp; nqs=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp;q=&nbsp;&nbsp; 0.00000000&nbsp; 0.00000000&nbsp; 0.00000000<br>&nbsp; reading force constants from file alh3.dyn2<br>&nbsp; nqs=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<br>&nbsp;q=&nbsp;&nbsp; 0.00095372&nbsp; 0.00000000&nbsp; 0.05148507<br>&nbsp;&nbsp; 1&nbsp;&nbsp; 1&nbsp;&nbsp; 1<br><br><span style="background-color: rgb(185, 185, 185);">&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span><br style="background-color: rgb(185, 185, 185);"><span style="background-color:
 rgb(185, 185, 185);">&nbsp;&nbsp;&nbsp;&nbsp; from init : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1</span><br style="background-color: rgb(185, 185, 185);"><span style="background-color: rgb(185, 185, 185);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nc already filled: wrong q grid or wrong nr</span><br style="background-color: rgb(185, 185, 185);"><span style="background-color: rgb(185, 185, 185);">&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span><br style="background-color: rgb(185, 185, 185);"><br style="background-color: rgb(185, 185, 185);"><span style="background-color: rgb(185, 185, 185);">&nbsp;&nbsp;&nbsp;&nbsp; stopping ...</span><br><br style="background-color: rgb(255, 0, 0);"></span>&gt;************************************************************************<br><BR><BR><p>&#32;



      <hr size=1><a href="https://member.cn.yahoo.com/cnreg/reginfo.html?id=89034" target=blank>ÑÅ»¢ÓÊÏ䣬ÖÕÉú»ï°é£¡</a>