<div dir="ltr"><div><div><div><div><div><div><div>Dear Shiva,<br><br><br><div>>I have 2 questions:<br></div>>1.I want to calculate bulk
 modules for the orthorhombic structure, but ev.x just have sc,bcc,fcc 
and hex, so what shoul I do?<br>
>2.Is the output of ev.x is trustable? I asked this, because it depends on how many numbers are in the input files!!!<br>>please help me in this.<br>>thanks.<br><div><div><div><br></div></div></div><br></div><br>
</div>I think you cannot use ev.x as a black box. The fitted bulk moduli and other parameters depend on the number of points, the accuracy of the calculated energies, and the range of volumes. Fitting is an art. <br><br></div>
Zero temperature calculations of bulk modulus can be compared with values derived from measurements of ultrasound speed, where the volume variations are tiny, much smaller than 1 %. However, for these volue variations the energy changes are also tiny and you need very well converge calculations (very large cutoffs, many k-points). Hence, a compromise must be achieved. <br>
</div><div>On the other hand, many e-v data improve statistics, check the errors given by the fitting program. If ev.x doest not provide errors, then use gnuplot or Origin. <br></div>Use volume variations between 1 and 10 %, and check the dependence of the bulk modulus with the convergence parameters and with the number of k-points and number of data. Check if the e-v fitting function is able to reproduce the e-v data that are outside the fitting range. When the results are carefully converged, everything matches. <br>
<br></div>I use the following methods to evaluate if a fitting is well done. Plot the difference g(V)=E(V)-f(V), where E(V) are the calculated data and f(V) is the fitting function. If g(V) looks random, the fitting proccess is fine, and the statistical errors  given by the fitting program are meaningfull.  If g(V) is not random, then the functional form of the fitting function is not adecuate, or the calculated values are biased (maybe due to poor convergence). <br>
<br></div>Best regards<br></div>eduardo<br><div><div><div><br><div><div><div><div><div><br>Eduardo Menendez Proupin<div> 
        <p>
    <b>Grupo de Cálculos Cuánticos</b><br>
    E.T.S.I.Telecomunicación<br><br>Universidad Politécnica de Madrid<br><br>Telefono: +34 91 549 57 00 Extensión 2310<br><br>URL: <a href="http://www.gdc.tat.upm.es" target="_blank">http://www.gdc.tat.upm.es</a></p></div>
<br clear="all"><div>
<div><br>Direccion permanente:<br>Departamento de Fisica, Facultad de Ciencias, Universidad de Chile <br>URL: <a href="http://www.gnm.cl/emenendez" target="_blank">http://www.gnm.cl/emenendez</a></div></div>
</div></div></div></div></div></div></div></div></div>