Please give a reference for the code if you find it useful.

Thermodynamics properties were calculated using the QHA code written by 
Eyvaz Isaev which is distributed as a part of Quantum Espresso code, see 
http://www.quantum-espresso.org
