I have a question, and I will be really grateful if someone helps me. I have a polynomial equation for internal energy which I calculated by integration an equation of state formula, which is based on density. But, because I calculated this using integration one integration constant which is temperature dependent( based on other articles) that I don't know how can I formulate it to have its magnitude to calculate internal energy at other range of data. My simulation box has contained 200 molecules of ionic liquid with one negative ion( PF6) and a positive one( butyl methyl imidazolium). Because according to internal energy equation at zero density internal energy is equal to the integration constant, we considered it as ionic liquid internal energy at ideal gas state. With all those in mind, how can I use a degree of freedom of rotational, vibrational, and translational to formulate this integration constant dependent of temperature to use it in other range of data? Or, is there any other method to formulate it? Thanks very much in advance