<div dir="ltr">Dear Wien2k users,<div><br></div><div>Hello, I am currently working on obtaining transition dipole moment values of wurtzite ZnO along a path using optic.</div><div><br></div><div>I have obtained some results but to check whether they are reasonable, I have been trying to obtain the refractive index values using x kram.</div><div><br></div><div>It seems that if I perform x kram after obtaining the momentum matrix elements along a certain path, (specified by xcrysden.klist for bandstructure calculations) I get incorrect refractive index values. That is I performed something like the following commands.</div><div><br></div><div>x lapw1 -band</div><div>x optic</div><div>x joint</div><div>x kram</div><div><br></div><div>and opticplot_lapw to see the refractive index curve. Am I not supposed to do this and perform x optic on a regular k mesh?</div><div><br></div><div>I also tried the same thing without the bandstructure calculation part, i.e.</div><div><br></div><div>x kgen</div><div>x lapw1</div><div>x lapw2 -fermi</div><div>x optic</div><div>x joint</div><div>x kram</div><div><br></div><div>and the results seemed to be more reasonable but still the refractive index values were about 0.2 lower than the known values. Is this within errors of DFT calculation? Or is it too large and I may have done something wrong?</div><div><br></div><div>Also in the end I need to obtain the momentum matrix elements using the first way, and when I do it should I use a finer k mesh before the band calculation? That is,</div><div><br></div><div>x kgen</div><div>x lapw1 -band</div><div>x optic</div><div><br></div><div>Right now I am just using a value of 1000 used in init_lapw. Could there be a reasonable value (or minimum value) of k points that should be used?</div><div><br></div><div>Thank you very much for reading the long post, and any help will be greatly appreciated.</div><div><br></div><div>Yours Sincerely,</div><div>Yong Woo Kim</div><div><br></div></div>