a function which can collect all the infomation of a subsolution phase has been implemented. the next step is to construct the gibbs energy.