T i N converged to 10 -6 eV with the tetrahedron method with Blöchl corrections. The Brillouin zone (BZ) integration was sampled by using a 9×9×9 Monkhorst-Pack k-points grid for the calculations of relaxation and electronic structure. The elastic constants C i j were calculated by an automatic procedure implemented in VASP. The cubic elastic constants for disordered Ti To study the interatomic chemical bonding in multicomponent systems based on TiN, we used crystal orbital Hamilton population (COHP) analysis. This technique adopted for plane-wave electronic structure calculations (projected COHP) and realized in the Local Orbital Basis Suite Towards Electronic-Structure Reconstruction (LOBSTER) code allows deriving an energy-resolved local bonding analysis and partitioning the band-structure energy into bonding and antibonding contributions from the PAW electronic structure calculations. Similar to the density of states (DOS), where the energy integration up to the Fermi level gives the number of electrons, the energy integration of the COHP for a pair of atoms indicates bond strength in terms of their contribution to the band structure energy.