Thermodynamic properties of rhodium—A first principle study

The high-pressure and high-temperature thermodynamic properties of rhodium (up to 2000 K and 300 GPa) are presented using the first principle approach within the quasi-harmonic approximation. The thermal Helmholtz free energy includes the contribution of both phonon vibrations and electronic excitations. The performance of three popular exchange-correlation functionals—local density approximation [Perdew et al., Phys. Rev. B 23, 5048 (1981)], Perdew–Burke–Ernzerhof generalized gradient approximation (PBE) [Perdew et al., Phys. Rev. Lett. 77, 3865 (1996)], PBE modified for dense solids [Perdew et al., Phys. Rev. Lett. 100, 136406 (2008)] are shown. The simulated thermal expansion coefficient, isobaric heat capacity, mode-Grüneisen parameter, thermodynamic average Grüneisen parameter, and bulk modulus are compared with the available experimental and theoretical reports. The contribution of thermal electronic excitations to the obtained thermodynamic parameters is significant at low pressure and high temperatures, except in bulk modulus, where it is small. The pressure-dependent elastic constant coefficient (Cij) and the Debye temperature are computed at 0 K. The Pugh ratio calculated from Cij indicates that rhodium undergoes brittle to ductile transitions at an average pressure of 7.45 GPa.