Elastic constants of cubic crystals

In this paper we present details of our developed open source software, cubic-elastic, for the calculation of the elastic constants (ECs) of cubic crystals. The comparison of the calculated ECs for various types of cubic systems by this software with those from the other available softwares as well as experimentally measured results confirms that our code can predict reliable results. The success of our code originates from its use of single deformation. The other codes usually use rhombohedral strain (RS). RS leads to 3B(0) + 4C(44) expression. Hence, RS systematically adds error to the C-44 through the bulk modulus calculations, and thereby may not be mathematically an appropriate approach. The total energy is accurately calculated by the WIEN2k within the highly accurate full-potential (linearized) augmented plane-waves plus local orbitals method. The ECs are calculated by the second-order derivatives of the fitted polynomials to the calculated total energies with respect to the elements of strain tensors at zero strains. We have presented the theoretical background and methodology of the cubic-elastic. We have validated the software by taking a variety of cubic samples into consideration and calculated their ECs. The zero bulk error calculations show that the results obtained from the cubic-elastic are in good agreement with the available experimental data and the previous theoretical results and predicts the sign of elastic constants correctly. The calculated Cauchy's pressure (C '') and Poisson's ratio (v) of LaS predict that it is an ionic compound. This prediction is in agreement (disagreement) with the previous ionic (covalent) bonds prediction deduced from previous v (C ''). (C) 2014 Elsevier B.V. All rights reserved.