Dear All:<div><br></div><div>I used ph.x to calculate the phonon frequencies at GAMMA point of Ta4Pd3Te16 (its structure can be found at http://www.springermaterials.com).  I had relaxed the structure so that all forces was little than 1mRy/bohr. </div><div><br></div><div>The first time I used parameters :</div><div><span class="Apple-tab-span" style="white-space:pre">       </span>...</div><div> <span class="Apple-tab-span" style="white-space:pre">    </span>conv_thr=1.0d-7 </div><div><span class="Apple-tab-span" style="white-space:pre">        </span>ecutwfc=40</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>ecutrho=200</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>...</div><div>in self consistent calculation and parameters:</div><div><span class="Apple-tab-span" style="white-space:pre">     </span>...</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>tr2_ph=1.0d-12,  <br><span class="Apple-tab-span" style="white-space:pre">    </span>ldiag = .true.</div><div><span class="Apple-tab-span" style="white-space:pre">       </span>start_irr = 1</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>last_irr = 33</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>...</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>/</div><div><span class="Apple-tab-span" style="white-space:pre">    </span>0 0 0</div><div>in phonon calculations. Note that irrep1-irrep22 are Ag and irrep23-irrep33 are Bg (which are raman active), the rest irreducible representations are not raman active.  I got a reasonable results which can be comparable with experiment but with an exception that the lowest frequency was imaginary:</div><div><div>     freq (  1 -  1) =        -28.6  [cm-1]   --> A_g             R  </div><div>     freq ( 38 - 38) =         30.0  [cm-1]   --> A_g  
           R  </div><div>     freq ( 39 - 39) =         36.3  [cm-1]   --> B_g             R  </div></div><div>     ...</div><div><br></div><div>Therefore I enhanced my precision as:</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>...</div><div> <span class="Apple-tab-span" style="white-space: pre;">  </span>conv_thr=1.0d-8 </div><div><span class="Apple-tab-span" style="white-space: pre;">      </span>ecutwfc=40</div><div><span class="Apple-tab-span" style="white-space: pre;"> </span>ecutrho=300</div><div><span class="Apple-tab-span" style="white-space: pre;">        </span>...</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>------------------</div><div><div><span class="Apple-tab-span" style="white-space: pre;">      </span>...</div><div><span class="Apple-tab-span" style="white-space: pre;">        </span>tr2_ph=1.0d-14,  <br><span class="Apple-tab-span" style="white-space: pre;">  </span>ldiag = .true.</div><div><span class="Apple-tab-span" style="white-space: pre;">     </span>start_irr = 1</div><div><span class="Apple-tab-span" style="white-space: pre;">      </span>last_irr = 33</div><div><span class="Apple-tab-span" style="white-space: pre;">      </span>...</div><div><span class="Apple-tab-span" style="white-space: pre;">        </span>/</div><div><span class="Apple-tab-span" style="white-space: pre;">  </span>0 0 0</div><span>I expected to get a better result, but actually I got a worse one:</span></div><div><span><div>     freq (  1 -  1) =       -504.0  [cm-1]   --> A_g             R  </div><div>     freq (  2 -  2) =       -333.9  [cm-1]   --> A_g             R  </div><div>     freq (  3 -  3) =       -251.0  [cm-1]   --> A_g   &nbs
p;         R  </div><div>     freq (  4 -  4) =       -140.0  [cm-1]   --> A_g             R  </div><div>     freq ( 41 - 41) =         44.2  [cm-1]   --> A_g             R  </div><div><span class="Apple-tab-span" style="white-space:pre">        </span>......</div><div><br></div><div>I completely can't understand why. Can anybody help me ?</div><div><br></div><div>Thanks</div><div><br></div>--<br><span style="font-size: medium;">Zhida Song</span><br style="font-size: medium;"><span style="font-size: medium;">Institute of Physics</span><br style="font-size: medium;"><span style="font-size: medium;">Chinese Academy of Sciences</span><br style="font-size: medium;"><span style="font-size: medium;">P.O. Box 603, Beijing 100190, China</span></span></div><br><br><br>