thermal_expansion = phonopy_qha.get_thermal_expansion()
volume_temperature = phonopy_qha.get_volume_temperature()
heat_capacity_P_numerical = phonopy_qha.get_heat_capacity_P_numerical()
volume_expansion = phonopy_qha.get_volume_expansion()
gibbs_temperature = phonopy_qha.get_gibbs_temperature()
#phonopy_qha.plot_bulk_modulus()
#plt.show()
phonopy_qha.plot_qha()
plt.show()
phonopy_qha.plot_gruneisen_temperature()
plt.show()
phonopy_qha.plot_gibbs_temperature()
plt.show()
phonopy_qha.plot_heat_capacity_P_numerical()
plt.show()
phonopy_qha.write_gibbs_temperature()
phonopy_qha.write_heat_capacity_P_numerical()
qha_output = ArrayData()
qha_output.set_array('temperature', np.array(qha_temperatures))
qha_output.set_array('helmholtz_volume', np.array(helmholtz_volume))
qha_output.set_array('thermal_expansion', np.array(thermal_expansion))
qha_output.set_array('volume_temperature', np.array(volume_temperature))
qha_output.set_array('heat_capacity_P_numerical',
np.array(heat_capacity_P_numerical))
qha_output.set_array('volume_expansion', np.array(volume_expansion))
qha_output.set_array('gibbs_temperature', np.array(gibbs_temperature))
#qha_output.store()
