Thermal properties & other

Calculation of thermal properties from quasiparticles and MD
Generate a harmonic displacements for testing
Added function to generate MD input for VASP and LAMMPS from POSCAR cell
Fix the error in phonon occupancy calculation
Now fit temperature is calculated from phonon occupancy (using Bose-Einstein distribution instead of Maxwell Boltzmann)
Other minor improvements