Since there are four states, then the grand partition function of the system is
where α is the chemical potential
Then, the occupancy of the system is
Then using this equation,
and approximating Z_int to be kT/0.00018 eV, the model would look as that attached in the figure. That is the occupancy vs. pressure graph.
There are more occupancies when the oxygen is high (high pressure) especially in the lungs. Heme sites tend to be occupied by oxygen.