starting work on phase equilibrium calculations