<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div>Dear Wien2k Community,<br><br></div>   I have no experience in calculating phonon dispersion, but I am trying to calculate it using Wien2k 17.1 and Phonopy.<br></div>   I am getting "strange" results  for wurtzite ZnO and I would highly appreciate your comments. Basically, I got a nice plot for the 2x2x2 supercell with negligible negative frequencies for the acoustic modes at gamma. Trying to check convergence, I increased the supercell to 3x3x2, but the negative frequencies became not negligible. The next step was to increase the cell to 4x4x3 with more k-points, but the negative frequencies are still there. Moreover, some of the optical modes look like not converging.<br></div>   LO-TO splitting are not taken into account.<br></div>   Well, I will describe my steps and, please, give me some suggestions or indicate the errors.<br><br></div>1) Structural optimization:<br></div>   1.1) I optimized the wurtzite ZnO cell using the "[5]  VARY A and C (2D-case)" option.<br></div>   1.2) 3% of RMT reduction (1.65 for Zn and 1.42 for O)<br></div>   1.3) non-magnetic configuration<br></div>   1.4) LSDA<br></div>   1.5) Core/valence separation energy : -7.0 Ry (valence states for Zn: 3P (local orbital), 3D and 4S // valence states for O: 2S (local orbital) and 2P)<br></div>   1.6) RKmax: 8.0 and 8.5 (convergence tests)<br></div>   1.7) K-grid :  (8 8 4) and (9 9 4) (convergence tests)<br></div>   1.8) Emax: 1.5 Ry (no SOC)<br></div>   1.9) Gmax: 14<br>   1.10) MSR1a <br>   1.11) tolf : 0.05 (case.inM)<br>   1.12) run_lapw -ec 0.0001 -cc 0.0001 -fc 0.01 -i 200 <br><br></div>Results:<br></div>   Sigma ~ 1.7x10^{-3}<br></div>   Differences in optimized "a" and "c" for the two RKmax and K-grid sets ~ 0.004%<br></div>   When optimizing "u" with the optimized "a" and "c", the residual forces where lower than 3.2x10^{-3} mRy/Bohr<br><br></div>2) Supercell:<br>   2.1)0% of RMT reduction (1.97 for Zn and 1.69 for O)<br>   2.2) non-magnetic configuration<br>   2.3) LSDA<br>   2.4) Core/valence separation energy : -7.0 Ry (valence states for Zn: 3P (local orbital), 3D and 4S // valence states for O: 2S (local orbital) and 2P)<br>   2.5) RKmax: 8.0 <br>   2.6) K-grid : (4 4 2) for the 2x2x2 supercell<br>                     (3 3 2) for the 3x3x2 supercell<br>                     (4 4 2) for the 4x4x3 supercell (trying to check if the k-grid plays a role)<br>   2.7) Emax: 1.5 Ry (no SOC)<br>   2.8) Gmax: 14<br>   2.9) run_lapw -p -it -ec 0.0001 -cc 0.0001 -fc 0.01 -i 200<br><br></div>=> All the K-grids had "no shift" !<br><br></div>   Finally, two questions:<br><br></div>1) When I built the P1 1x1x1 supercell, I simply separated the equivalent atoms, preserving positions and the gamma angle equal to 120. Phonopy finds the correct original space group. Is this the correct procedure ?<br><br></div>2) There is a paper where the authors (Prof. Blaha among them) calculate the phonon dispersion for wurtzite ZnO using Wien93 (PRB 53, 11425 (1996)). In this calculation, the 3P states for Zn are core states. Do you believe that this is relevant ?<br><br></div>   I am sorry for the huge email !<br></div><div>   Thank you !<br></div>   All the best,<br></div>                 Luis<br><br></div>