Phase transformation for si gas hydrates


I am simulating the thermal dissociation of si methane hydrates using lammps. To know if the results of the simulation agrees the physical behavior of the system, I need to know if the melting of si hydrate is a first-order or a second-order phase transformation. Does anyone have any information on that?