I am studying the melting and recrystallization of Ge2Sb2Te5 by running NVT simulations at different temperatures. I would like to calculate the free energy of each output structure and am looking for the best method to do so. I was planning to use the 2pt analysis from the DoSPT code, however, I am studying the hexagonal phase of GST and from my knowledge DoSPT only supports cubic crystals. Are there any other post processing methods that can read my trajectories to get entropy? I could then use F = U - TS to get free energy.