<div dir="ltr">Dear users,<div>I am trying to calculate critical temperature for Titanium. I used the following</div><div>files:</div><div><a href="http://Ti.scf.in">Ti.scf.in</a></div><div><div> &control</div><div>    calculation='scf'</div><div>    restart_mode='from_scratch',</div><div>    prefix='Ti',</div><div>    pseudo_dir = './',</div><div>    outdir='./'</div><div> /</div><div> &system</div><div>    ibrav =  4,</div><div>    celldm(1) = 5.576203,</div><div>    celldm(3) = 1.587874,</div><div>    nat = 1,</div><div>    ntyp = 1,</div><div>    occupations = 'smearing',</div><div>    degauss = 0.03,</div><div>    smearing = 'cold',</div><div>    ecutwfc = 60.0,</div><div>    la2F = .true.,</div><div> /</div><div> &electrons</div><div>     mixing_beta = 0.7,</div><div>     conv_thr =  1.0d-15,</div><div>     diagonalization = 'cg',</div><div>/</div><div>ATOMIC_SPECIES</div><div> Ti  47.88 Ti.pbe-sp-hgh.UPF</div><div>ATOMIC_POSITIONS ANGSTROM</div><div> Ti 0.33333 0.66667 0.25000</div><div>K_POINTS {automatic}</div><div>8 8 8 0 0 0</div></div><div><br></div><div>Question 1) The calculation is finished without error, but in the output file</div><div>has the message  "No symmetry found</div><div>          (note: 23 additional sym.ops. were found but ignored</div><div>           their fractional translations are incommensurate with FFT grid)".</div><div>This is very interesting because I tested various inputs with different elements with hexagonal structure (Zr, Tc, Re...) the same message appear. So, is there any</div><div>problem or bug with the ibrav = 4 (hexagonal)?</div><div>In addiction, in the ph.x calculation using the input below, I found negative frequencies, these frequencies can be attributed to this message?</div><div><br></div><div>Question 2)<br></div><div>After the scf calculation, I performed ph.x using the file:</div><div><div> &inputph</div><div>  tr2_ph=1.0d-16,</div><div>  prefix='Ti',</div><div>  fildvscf='aldv',</div><div>  amass(1)=47.88,</div><div>  outdir='./',</div><div>  fildyn='Ti.dyn',</div><div>  electron_phonon='interpolated',</div><div>  trans=.true.,</div><div>  ldisp=.true.</div><div>  nq1=4, nq2=4, nq3=4</div><div> /</div></div><div>In this calculation, the file Ti.dy0 has 36 points and it was generate</div><div>36 elph.* files. However, the kpoints.x generate 50 k-points:</div><div><div>bravais lattice  >> 4</div><div>     filout [mesh_k]  >> </div><div>     enter celldm(3)  >> 1.587874 </div><div>     mesh: n1 n2 n3   >> 8 8 8</div><div>     mesh: k1 k2 k3 (0 no shift, 1 shifted)  >> 0 0 0</div><div>     write all k? [f] >> </div><div>     # of k-points   ==    50  of   512<br></div></div><div><br></div><div>If I use mesh: n1 n2 n3   >> 4 4 4, I have 12 k-points</div><div>So, I expected that the number o k-points generated by the kpoints.x is the</div><div>same of the obtained by the ph.x. Is this corrected and how can I improve this?</div><div><br></div><div>Thanks in advance.</div><div>Osmair V. Oliveira</div><div>Brazil</div><div><br></div></div>