<div dir="ltr"><div>Dear Sir</div><div><br></div><div>I tried to reproduce results for elastic properties of CsCl with Wien2k_16.1 as mentioned in original paper  of Prof. R. Golesorkhtarbar et al. 2013</div><div>As mentioned in the paper the authors suggested to use RMT=2.0. k-mesh= 15 x 15 x 15, RKmax=9 with PBESol.</div><div>They suggested an optimal value of Langragian deformation between [0.05-0.08] (for higher order any value of neta is possible).</div><div><br></div><div><br></div><div>1. The original structure is: <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CsCl.struct_orig">CsCl.struct_orig</a>  optimised with run_lapw -ec 0.00001   -p -cc 0.0001  -red 5%. kmesh 15 15 15. Forces are not present in struct file so I left -fc switch.<div><br></div><div>2. The well optimised structure is: <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CsCl.struct_orig">CsCl.struct</a></div></div><div><br></div><div>3. The entire optimised directory: <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CsCl.tar.gz">CsCl.tar.gz</a></div><div><br></div><div>I tested two  case with the optimised struct mentioned in link 2 to calculate the elastic properties:</div><div><br></div><div>A. RMT=2.5, k-mesh= 15 x 15 x 15, RKmax=9 with PBESo, -ec 0.0001 eta=0.05, deformation =7 with energy approach.</div><div><br></div><div>I got the result: <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/ElaStic_2nd.out">ElaStic_2nd.out</a> with strain and order <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/ElaStic_2nd.in">ElaStic_2nd.in</a>. these obtained results are very far from the results reported in paper: Second last colum of Table10.</div><div><br></div><div>The entire directory is: <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CsCl_elast.tar.gz">CsCl_elast.tar.gz</a></div><div><br></div><div>B. For second approach I used RMT=2.0 at suggested in the paper and eta=0.08, deformation =7 with  -ec 0.00001 -cc 0.0001</div><div>The summery of results are here: <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/Energy-vs-Strain.tar.gz">Energy-vs-Strain.tar.gz</a> with input and output parameters.</div><div><br></div><div>In this case results are somewhat improved but still far from the reported data.</div><div><br></div><div><br></div><div>Could you please help me why I am not getting results in the vicinity?</div><div>It would be much helpful if Prof. Peter pass the information to the developers of Elastic code.</div><div><br></div><div>Let me know if I can provide any additional information.</div><div><br></div><div>Warm regards</div><div><br></div><div>Fatima <br></div><div><br></div><div><br></div></div>