Optimized ground states parameters for all the three osmium double perovskites.

## Abstract

In this chapter, Osmium-based double perovskites Ba2XOsO6(X = Mg, Zn, Cd) have been investigated for their magnetic structure, electronic, elastic, mechanical and thermodynamic belongings. These materials have been recently reported experimentally for their magnetic structure. Here, we report the first successful ab initio calculations on the physical properties of these materials. The structural optimization for these Ba2XOsO6(X = Mg, Zn, Cd) double perovskite compounds has been finalized within density functional theory via full potential linearized augmented plane wave (FP-LAPW) method. The structural investigation exposes the ferromagnetic phase stability of these compounds. The spin-polarized electronic and magnetic properties were calculated within generalized gradient approximation (GGA), Hubbard approximation (GGA+U) and modified Becke-Johnson approximation (mBJ). The electronic profile establishes the half-metallic nature for all the three compounds. The total spin magnetic moment was found to be an integer value of 2 μb. The elastic constants have been calculated and used to predict mechanical stuffs like Shear modulus (G), Poisson ratio (v) and anisotropic factor. The calculated B/G and Cauchy pressure (C12-C44) both characterize these materials as brittle. The thermodynamic parameters like heat capacity and Debye temperature have been predicted in the temperature range of 0–1000 K.

### Keywords

- Ba2XOsO6 (X = Mg
- Zn
- Cd)
- spintronics
- ferromagnetic
- elastic
- mechanical behavior
- thermodynamics

## 1. Introduction

The need of advanced materials with novel properties for industrial and technological use has strained the material community to have a deep and appropriate understanding of the periodic table elements, along with their combinations. Therefore, materials community consequently observes the vital changes in innovative designing of novel materials. A tremendous increase in simulation power, along with algorithmic improvements in quantum theory allows one to have well-organized and exact quantum mechanical calculations. This has hence stretched the computing power to such extent that those properties of materials which were once observed extremely difficult are now easily being calculated with a great precision [1, 2]. From last few years work on perovskites especially double perovskites has geared up due to their vast technological applications and displaying multifunctional properties. The general formula of perovskites is looked as ABO_{3}, where “A” and “B” are cations and “O” is oxygen anion. The charge of “A” and “B” cations can vary in the original Perovskite. Double perovskites are potential members of this diverse perovskite family having different structures, composition and properties. The double perovskite compounds having a general formula A_{2}BB′O_{6} have benefited the material community because of great technological applications including spintronic materials, multi-ferroic materials, half-metallic materials, ferromagnetic materials, magneto-dielectric materials, magneto-optic materials, insulating ferrimagnetism [3, 4, 5, 6, 7, 8]. The double perovskite family exhibits a wide range of magnetic behaviors, like actually simple antiferromagnets presented by (Ba_{2}LiOsO_{6}), ferromagnets (Ba_{2}MgReO_{6} and Ba_{2}NaOsO_{6}), spin singlet ground states for (Ba_{2}YMoO_{6}), and spin glasses(Sr_{2}CaReO_{6}) [9, 10, 11, 12, 13]. Osmium based double perovskites especially Sr_{2}FeOsO_{6}, SrFeCaOsO_{6}, Ca_{2}FeOsO_{6} have been extensively investigated for puzzling magnetic behavior [14, 15, 16, 17, 18, 19]. Perovskites with ABO_{3} structure like BaPuO_{3}, SrPuO_{3}, BaAmO_{3}, SrAmO_{3}, EuGaO_{3}, EuInO_{3} etc. have been reported for spintronic applications [20, 21, 22, 23, 24]. Further materials like BaMoO_{3}, SrMoO_{3}, XReO_{3} (X = Rb, Cs, Tl), SnTaO_{3}, PbMoO_{3} [25, 26, 27, 28] have been recently reported for fuel cell applications. Numerous investigations have also been reported on halide perovskites for solar cell applications like MAPbBr_{3} and MAPbI_{3} [29, 30, 31, 32, 33, 34, 35, 36].

Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) cubic double perovskites have been recently reported in space group Fm-3 m (225). The complete details of the respective experimental lattice parameters are given in Table 1. The Ba atoms are located at 8c (1/4, 1/4, 1/4) of the unit cell, X atoms are at position 4b (0.5, 0.5, 0.5), Os atoms are positioned at 41a (0, 0, 0) and O atoms at 24e (X, 0, 0) (X = 0.233, 0.238, 0.234) respectively for Mg, Zn, Cd [37]. Further double perovskites of the variant A_{2}BB′O_{6} like Ba_{2}MgReO_{6}, Sr_{2}MnTaO_{6}, Ba_{2}InTaO_{6} and many more have also been reported for electronic, magnetic, mechanical, optical, thermoelectric and thermodynamic investigations [38, 39, 40, 41, 42, 43, 44, 45, 46, 47]. As far Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) compounds are concerned which belongs to the same variant not much attention has been paid towards these perovskites so far, regarding the above mentioned characteristic properties. Hence, in the present work an attempt to predict the properties for these double perovskite has been made and to check out their potential applications. The most successful density functional theory (DFT) has been employed for the investigation of magnetic, electronic, elastic, mechanical and thermo-physical behavior. For the investigation of thermo-physical behavior quasi harmonic Debye model [48, 49] has been used for the prediction of important parameters like specific heat, thermal expansion, Debye temperature, Grüneisen parameter etc.

Parameter | Present | Other | Present | Other | Present | Other | |
---|---|---|---|---|---|---|---|

Ba_{2}MgOsO_{6} | Ba_{2}ZnOsO_{6} | Ba_{2}CdOsO_{6} | |||||

Lattice Constant | 8.1548 | 8.07 [37] 8.08 [61] 8.06 [61] | 8.19 | 8.09 [37] 8.09 [61] 8.06 [61] | 8.388 | 8.31 [37] 8.325 [61] | |

Volume | 914.90 | 929.62 | 995.96 | ||||

B | 150.93 | 144.23 | 140.14 | ||||

B^{′} | 4.44 | 4.13 | 4.7 | ||||

Bond length | |||||||

Os-O | 1.95 | 1.95 | 2.00 | ||||

Mg, Zn, Cd-O | 2.11 | 2.15 | 2.18 | ||||

Ba-Ba | 4.07 | 4.10 | 4.19 | ||||

Ba, Mg, Zn, Cd | 3.53 | 3.55 | 3.62 | ||||

Os-Mg, Zn, Cd | 4.07 | 4.10 | 4.19 | ||||

E_{0} | −68425.9 | −71617.3 | −79217.0 |

## 2. Computational details

The computational technique used during the calculations process is based on full-potential linearized augmented plane wave (FP-LAPW) [50] method based upon density functional theory (DFT) [51] as employed in WIEN2K. For structural optimization generalized gradient approximation (GGA) scheme of Perdew, Burke and Ernzerhof (PBE) [52] has been used. For electronic and magnetic calculations in addition to (GGA), Hubbard approximation (GGA + U) [45] and modified Becke-Johnson (mBJ) [53] has been used. For GGA + U approach the incorporation U- term can be done by various methods [54, 55]. In the present work we have used self-interaction correction method (SIC) [56] as implemented in WIEN2K. The value of U_{eff} was varied from 1 to 5 eV and J was set to 0, so as to properly adjust the Os-d in density of states. The final U value used throughout the calculations was set to 2.00 eV [57]. For precise energy convergence the value of R_{MT}K_{max} was taken 7, where R_{MT} is the small atomic radius in unit cell and K_{max} denotes the size of the largest **k** vector in the plane wave expansion. The value of L_{max} was taken as 10, and G_{max} = 12 (a.u.)^{−1}. The energy and charge convergence criterion is considered when the total energy is stable within 0.001 Ry and the charge difference is less than 0.001e/a.u.^{3} per unit cell. A mesh of 1000 K points is considered for Brillouin zone integration via tetrahedral method [58]. The elastic constants were calculated using the scheme developed by Charpin [59] as integrated in WIEN2K package. The thermodynamic parameters have been calculated using quasi-harmonic Debye model [48, 49] for the pressure and temperature dependency of some essential thermodynamic parameters. In this model the Gibbs function takes the form;

where E(V), P(V),

In the above equation

The non-equilibrium Gibbs function G*(V, P, T) can be minimized with respect to volume V;

Solution of Eq. (5) gives a detailed information about the thermodynamic quantities like thermal expansion _{V}, heat capacity at constant pressure C_{P}, given respectively by;

In Eq. (8)

## 3. Results and discussion

### 3.1 Structural properties

The optimized volume for all the three compounds has been made by fitting the total energy as a function of its cell volume using Birch–Murnaghan’s equation of state [60]. Marjerrison et al. [37] have recently reported all the three compounds in cubic B1-phase space group Fm-3 m (225). The Ba atoms are located at position 8c (0.25, 0.25, 0.25), Mg, Zn, Cd at 4b (0.5, 0.5, 0.5), Os at 4a (0, 0, 0) and O atoms are sited at24e (X, 0, 0) (X = 0.233, 0.238, 0.234) respectively for Mg, Zn, Cd. The geometry and structural optimization has been carried in Non-magnetic (NM), ferromagnetic (FM), and anti-ferromagnetic (AFM) phases. The ground state energy was found lowest for all the three compounds in the ferromagnetic phase as presented in Figure 1(a–c), and thus a stable configuration.

The optimized ground state lattice constants are close to available experimental and theoretical results. The ground state parameters like bulk modulus (B_{0}), lattice constant (

### 3.2 Elastic and mechanical properties

In the present work the Charpin method has been employed for the calculation of elastic constants C_{ij} (C_{11}, C_{12}, C_{44}) values as implemented in WIEN2K. The values of elastic constants were obtained by calculating the total energy as a function of volume-conserving strains. The value of the elastic constants and mechanical properties are summed up in Table 2. The calculated values of elastic constants properly satisfy the criteria for cubic elastic constants and ensures the stability C_{11}–C_{12} > 0, C_{11} > 0, C_{44} > 0, (C_{11} + 2C_{12}) > 0, C_{12} < B < C_{11} [47]. The Poisson’s ratio (*ν*), Young’s modulus (E), and Shear modulus (G) are calculated by using [62, 63, 64] and presented in Table 2. According to Hill [65] average shear modulus, *G* is defined as arithmetic mean of Voigt, *GV* and Reuss, *GR* values. Young’s modulus (E) deals with the stiffness of the material. The obtained value of (E) was found to be 215.75, 190.82, 169.83 GPa respectively for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd). Thus large value of (E) provides a clear indication that these compounds will behave as tough materials. Ba_{2}MgOsO_{6} has the largest value of (E) as compared to other perovskites under consideration in this study. The reason for the decreasing value of (E) is the Bulk modulus which has also a decreasing trend as one goes with X position from Mg to Cd via Zn. The B/G ratio is the measure of ductility and brittleness of a material. According to Pugh [66], a material is brittle if the ratio B/G < 1.75 and is ductile if B/G > 1.75. The B/G ratio for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd), was calculated to be 1.765, 1.95, 2.00 respectively, which is higher than the limit value for all the three compounds, thus all the three compounds will show ductile nature. Cauchy pressure (*C12*–*C44*) also helps to estimate the ductility and brittleness of a material. The positive value of (*C12*–*C44*) portrays a material as ductile and negative value as brittle. The calculated value was also found to be positive for all the three compounds. Hence both B/G value and Cauchy pressure verifies the ductile nature for all the three perovskites Ba_{2}XOsO_{6}(X = Mg, Zn, Cd).

GGA | BaMgOsO_{6} | BaZnOsO_{6} | Ba_{2}CdOsO_{6} |
---|---|---|---|

C_{11} | 262.60 | 240.70 | 232.07 |

C_{12} | 86.90 | 98.85 | 84.32 |

C_{44} | 84.10 | 76.91 | 61.66 |

B | 150 | 145.59 | 132.95 |

G_{V} | 85.61 | 74.51 | 66.54 |

G_{R} | 85.57 | 74.39 | 66.02 |

G | 85.50 | 74.45 | 66.28 |

E | 215.75 | 190.82 | 169.83 |

ν | 0.2617 | 0.2815 | 0.2871 |

B/G | 1.7652 | 1.95 | 2.00 |

C_{12}–C_{44} | 2.8 | 21.94 | 22.66 |

A | 0.957 | 1.08 | 0.8346 |

T_{m} | 2105 ± 300 | 2100 ± 300 | 1925 ± 300 |

Zener anisotropy factor ‘A’ is the property of a material to show altered characteristic in various direction of its structure. As per this a material is isotropic if and only if ‘A’ factor has unit value or otherwise anisotropic. The calculated value of ‘A’ for the compound was found to 0.975, 1.08, 0.83 which is less than unity for Ba_{2}MgOsO_{6} and Ba_{2}CdOsO_{6} and greater than unity for Ba_{2}ZnOsO_{6}, hence in all the three cases deviating from unity, thus the materials will present anisotropic nature. Poisson’s ratio (ν) describes the nature of bonding forces. The upper and lower limits of Poisson’s ratio are 0.25 and 0.50 [20, 21, 22, 23, 24]. The (ν) value varies from material to material. For covalent materials, (ν) has a typical value of 0.1, for ionic materials (ν) = 0.25 and for metallic materials the value (ν) =0.33. The value of Poisson’s ratio for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) was calculated to be 0.261, 0.281 and 0.287respectively, which lies close to 0.25 and hence suggest a higher ionic behavior as inter-atomic bonding for these compounds. The obtained values of elastic constants have also been used to predict, one important thermodynamic parameter known as melting temperature [26, 27, 28]. The calculated value of melting temperature was found 2100 ± 300, 2105 ± 300, 1925 ± 300 K respectively for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd). The calculated values of elastic constants, mechanical properties including melting temperature are grouped in Table 2.

### 3.3 Electronic and magnetic properties

For electronic structure calculations spin resolved band structure and density of states have been plotted using different correlation potentials. GGA calculated lattice parameter has been used to plot band structure and density of states within GGA, GGA + U and mBJ. These band structures and density plots usually deliver a decent understanding of the electronic contour of a material. The combination of different methods for band structure and density of state plots has been done as to understand the variation of results within different correlations. Figures 2(a–c), 3(a–c) and 4(a–c) represent the spin included band structures within GGA, GGA + U and mBJ respectively for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd).

It is clear from these figures that the band profile for all the three compounds at the Fermi level is almost similar for all the approximation, presenting 100% of spin polarization. The Fermi level is set at 0 eV, separating the valance band maximum (VBM) from the conduction band minimum (CBM) in all figures. For spin up states the Fermi level remains fully occupied presenting metallic nature for all the three compounds and for spin down states the Fermi level remains completely vacant falling in a gap and thus generating a gap between(VBM) and (CBM), presenting the semi-conducting nature for the compounds. In case of Ba_{2}MgOsO_{6} within GGA, GGA + U and mBJ respectively, the (VBM) lies on symmetry points ‘Γ’ at −1.30 eV, −1.1 eV and −1.3 eV, and CBM lies on symmetry point ^{‘}X^{’} at 0.001, 0.7 and 1.2 eV respectively within GGA, GGA + U and mBJ. Hence in all the three cases the (VBM) and (CBM) lie on ‘Γ’ and X point making the compound indirect band gap semi-conductor in spin down states. The band gap value changes as we apply U and mBJ, the value of band gap found in GGA, GGA + U and mBJ are 1.3, 1.8 and 2.5 eV respectively for Ba_{2}MgOsO_{6}. For Ba_{2}ZnOsO_{6} the valance band maxima (VBM) lie on symmetry points ‘Γ’ at −1.40, −1.2, −1.4 eV respectively in GGA, GGA + U and mBJ, and the conduction band minimum (CBM) lies on symmetry point ‘X’ at 0.2, 0.7, and 1.0 eV respectively in GGA, GGA + U and mBJ, thus generating an indirect band gap of 1.6, 1.9 and 2.4 eV respectively for GGA, GGA + U and mBJ. Similarly for Ba_{2}CdOsO_{6} the valance band maxima (VBM) lie on symmetry points ‘L’ at −1.00, −0.8, −0.9 eV respectively in GGA, GGA + U and mBJ, the conduction band minimum (CBM) lies on symmetry point ‘X’ at 0.5, 0.8, and 1.4 eV respectively in GGA, GGA + U and mBJ, thus generating an indirect band gap of 1.5, 1.6 and 2.3 eV respectively for GGA, GGA + U and mBJ. Thus from the band structure calculations 100% of spin polarization at Fermi level is observed. The compounds behave as metallic for spin up states and semi-conducting for spin down states. The overall band picture presents half-metallic nature for all the three compounds.

For the further explanation of the band picture, total density of states (TDOS) and partial density of states (PDOS) have been plotted. The spin included combined TDOS shown in Figure 5(a–c) disclose the same results as presented by band structure plots presenting metallic nature for spin up states and semi-conducting for spin down states for all the three approximations and hence overall half-metallic nature. The DOS peaks are found to increase in case of GGA + U and mBJ. The contribution to the TDOS picture has been represented by the partial contribution to the DOS diagram as depicted in Figure 6(a–c) for both spin up and down states within GGA + U. The (PDOS) has been plotted for (Ba-‘s’, ‘p’, ‘d’, ‘f’), (Mg- ‘s’, ‘p’, Zn-‘s’, ‘p’, ‘d’), (Cd-‘s’, ‘p’, ‘d’), (Os- ‘s’, ‘p’, ‘d’, ‘f’) and O-‘s’, ‘p’ states. From these figures. it is clear that the metallic nature in spin up states for the compound in all the approximations is due to the Os-‘d’ states which are present at Fermi level with a small contribution of O-‘p’ states hybridized with one another and in case of spin down these ‘d’- states of Os and ‘p’-states of O are pulled inside the conduction band, thereby generating a gap in spin down states. Thus the spin included band profile, TDOS and PDOS results display that Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) all present half-metallic nature.

In order to check the magnetic nature of the compounds the total and partial magnetic moments have been calculated with GGA, GGA + U and mBJ. The total magnetic contribution is obtained as the summation of the partial moments of individual elements and the interstitial moments. The total magnetic moment obtained in all approximations is nearly same for all the compounds Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) equal to an integer value 2 _{B} shown in Table 3.The main contribution to the total magnetic moment is mostly found from Osmium atoms. The partial moment of Os element shows a great variation on the application of Hubbard U and mBJ potentials. Hence it is clear that the ferromagnetic nature and large value of total magnetic moment for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) is mainly due to Os atoms. The values of interstitial, partial and total magnetic moments are present in Table 3. Thus the large and integer value of magnetic moment of 2 _{B} further verifies the half-metallic and ferromagnetic nature for Ba_{2}MgOsO_{6}. The integer value of magnetic moment is one of the criteria for the half-metallic nature of a compound [22, 23].

Compound | Method | M_{int} | M_{Ba} | M/_{Mg,Zn,Cd} | M_{Os} | M_{O} | M_{Tot} |
---|---|---|---|---|---|---|---|

Ba_{2}MgOsO_{6} | GGA GGA + U mBJ | 0.38 0.35 0.22 | 0.01 0.01 0.00 | 0.00 0.00 0.00 | 1.14 1.25 1.41 | 0.07 0.06 0.05 | 2.0 2.0 2.0 |

Ba_{2}ZnOsO_{6} | GGA GGA + U mBJ | 0.37 0.34 0.21 | 0.01 0.01 0.00 | 0.00 0.00 0.00 | 1.15 1.26 1.42 | 0.07 0.06 0.05 | 2.0 2.0 2.0 |

Ba_{2}CdOsO_{6} | GGA GGA + U mBJ | 0.37 0.33 0.20 | 0.01 0.00 0.00 | 0.00 0.00 0.00 | 1.18 1.31 1.49 | 0.06 0.05 0.04 | 2.0 2.0 2.0 |

### 3.4 Thermodynamic properties

In order to check the thermodynamic behavior quasi-harmonic Debye approximation [26, 27, 28] has been employed to check the temperature and pressure variation of some noteworthy thermodynamic quantities like heat at constant volume (C_{v}), thermal expansion (α), Grüneisen parameter (γ), Debye temperature (θ_{D}) and also the bulk modulus variation for these double perovskites. The variation of these parameters has been investigated under pressure and temperature. The temperature has been varied from 0 to 1000 K and pressure ranges from 0 to 15 GPa, with the step size of 5 GPa pressure. In this range of temperature quasi harmonic Debye model remains unconditionally valid.

Figure 7(a–c) presents the variation of bulk modulus (B) with temperature at different pressure points respectively at 0, 5, 10, and 15 GPa. Our results present a clear decrease in bulk modulus with temperature and an increase is observed with pressure at different temperature values. The reason for this decrease of bulk modulus with temperature is that temperature reduces the hardness of a material, while pressure tends to increase the same.

Figure 8(a–c) depicts the variation of specific heat at constant volume (C_{V}) with temperature and pressure. One can have a clear understanding from the Figure 8(a–c) that the escalation of C_{V} is rapid under the lower temperature values of 0–300 K, but above 300 K a lethargic increase in C_{V} can be seen, which further becomes constant at high temperature at about 800 K beyond which, it follows the famous Dulong-Petit limit [67]. This variation of C_{V} for solids is a common observation. The calculated value of C_{V} at 300 K and 0 GPa of pressure for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) was found to be 223.89, 223.16, 225.04 J mol^{−1} K respectively.

Figure 9(a–c) shows the pressure and temperature dependence of thermal expansion coefficient, ‘α’ respectively for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd). It is clear from the figure that the value of ‘α’ increase with increasing temperature, the increase in ‘α’ is found to be rapid under low temperatures values and under higher temperatures values a sluggish increase in ‘α’ is observed. The main reason for the sluggish increase of ‘α’ under high temperature values may be the saturation of ‘α’ beyond 300 K. Pressure has a reverse effect on ‘α’, increasing pressure decreases the ‘α’. Under high pressure values ‘α’ falls rapidly, same as the increase is observed under low temperatures.

Grüneisen parameter (γ) describes the variation in vibrational frequency of a lattice under the influence of temperature and pressure [68]. Pressure and temperature variation of (γ) for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd) is plotted in Figure 10(a–c). The value of (γ) increases with increasing temperature and under pressure a reverse is observed, pressure decreases the value of (γ) and has a lowest value at 15 GPa of pressure. The predicted value of (γ) at 300 K and 0 GPa of pressure is 2.088, 1.965, 2.099 respectively for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd).

Debye temperature (*θD)* one of the most important thermodynamic parameter helps to exposes accurate presentation of material properties like specific heat capacity and thermal expansion and also provides the decent understanding about the features of a material under the influence of temperature and pressure. The Debye temperature variation with respect to temperature/pressure is presented in Figure 11(a–c) . From these figures it is clear that Debye temperature shows a decreasing trend with increasing temperature and an increasing trend with increasing pressure. The calculated value of Debye temperature at 300 K and 0 GPa is 446.08, 452.97, 435.12 GPa respectively for Ba_{2}XOsO_{6}(X = Mg, Zn, Cd).Some part of this work has been recently reported [53].

## 4. Conclusions

Ab initio calculations on electronic structure, magnetic, elastic, mechanical and thermodynamic properties of cubic double perovskite oxides Ba_{2}XOsO_{6} (X = Mg, Zn, Cd) Ba_{2}ZnOsO_{6} have been reported within density functional theory via full potential linearized augmented plane wave (FP-LAPW) method. The structural investigation reveals the ferromagnetic phase stability for these compounds. The spin polarized electronic and magnetic properties were calculated within generalized gradient approximation (GGA), Hubbard approximation (GGA + U) and mBJ (modified Becke-Johnson approximation). The electronic profile establishes half-metallic nature for these compounds and hence can strength the modern technological domain in terms of spintronic materials. The calculated total spin magnetic moment was found equal to 2 μ_{B} for all the three compounds. Thus these materials are also looked for magnetic materials. The elastic constants have been calculated and used to predict mechanical stuffs like Shear modulus (G) Poisson ratio (ν) and anisotropic factor (A). The calculated B/G and Cauchy pressure (C_{12}–C_{44}) both characterize the material as brittle. The thermodynamic parameters like heat capacity and Debye temperature have also been predicted in the temperature range of 0–1000 K.