Changed default values for lattice anisotropy to prevent NaN
Corrected sign error in energy calculations
Implemented coupling of materials to the phonon temperature