Use mode heat capacity implemented in thermal_property.py
Implementation of thermodynamic Gruneisen parameter
Fundtion of mode heat capacity is taken out from the class