MixPred

MixPred is the module for mixture isotherm prediction for the given pressure and composition. MixPred GUI is shown below:

MixPred

The isotherm fitting results from IsoFit and HeatFit can be directly loaded into MixPred module, and the user can specify the desired pressure range and mixture composition to calculate the mixture loadings.

Watch how to use MixPred here

MixPred has two models for mixture adsorption loading prediction which are as follows:

  • Extended Dual-site Langmuir Model (EDSL)

  • Ideal Adsorbed Solution Theory (IAST)

Extended Dual-site Langmuir Model (EDSL)

The extended dual-site Langmuir (EDSL) model is a generalization of the Langmuir model for multi-component mixture adsorption. The model is expressed as,

\[q_{mix,i}^{*} = \frac{q_{sat,1,i}b_{1,i}P_{i}}{1+\sum_{1}^{N}b_{1,j}P_{j}} + \frac{q_{sat,2,i}b_{2,i}P_{i}}{1+\sum_{1}^{N}b_{2,j}P_{j}}\]

where \(N\) is the number of components involved; \(q_{sat,1,i}\), \(q_{sat,2,i}\), \(b_{1,i}\), \(b_{2,i}\), and \(P_{i}\) are the single or pure component Langmuir isotherm parameters and partial pressure of component \(i\), respectively. The pure component Langmuir isotherm parameters for individual components can be obtained by fitting Langmuir isotherm model to the pure component isotherm data for the given component using either IsoFit or HeatFit. Note that the EDSL model is thermodynamically consistent only when the saturation capacities for each component are equal.

\[\begin{split}q_{sat, 1, 1}=q_{sat, 1, 2}=q_{sat, 1, 3}= ... q_{sat, 1, N} \\ q_{sat, 2, 1}=q_{sat, 2, 2}=q_{sat, 2, 3}= ... q_{sat, 2, N}\end{split}\]

Note

The EDSL model in MixPred is applicable only when the isotherm of each individual component is described using the Langmuir isotherm model.

Ideal Adsorbed Solution Theory (IAST)

IAST is a thermodynamic framework to calculate the mixture isotherms using pure component isotherm. The IAST is based on three fundamental assumptions:

  • The surface area of the adsorbent is equally accessible to all adsorbates.

  • The adsorbed phase is an ideal solution.

  • The adsorbent is homogeneous.

The solution of the IAST involves solving non-linear equations consisting of a rediced grand potential. For the component \(i\), the reduced grand potential \(\psi_{i}^{*}\) is given as,

\[\psi_{i}^{*} = \int_{0}^{P_{i}^{*}} \frac{q_{i}^{*}(P)}{P} \,dP\]

where \(P_{i}^{*}\) is the fictitious pressure and \(q_{i}^{*}\) is the equilibrium loading of component \(i\) as obtained from pure component isotherm model. The fictitious pressure \(P_{i}^{*}\) is the pressure for component \(i\) at which it exerts the same reduced grand potential as the other components in the mixture. \(P_{i}^{*}\) is related to partial pressure \(P_{i}\) of component \(i\) given as,

\[P_{i} = x_{i}\,P_{i}^{*} \qquad \mathrm{for}\, i=1,2,3,\ldots, N.\]

where \(N\) is the number of adsorbing components in the mixture and \(x_{i}\) is the adsorbed mole fraction of component \(i\). Note that,

\[\sum_{i}\,x_{i} = 1\]

IAST states that adsorption equilibrium is achieved when the reduced grand potential of each component is same,

\[\psi_{1}^{*} = \psi_{2}^{*} = \psi_{3}^{*} = ... \psi_{N}^{*}\]

In MixPred the IAST equations above are solved for \(2N\) unknows, \(P_{i}^{*}\) and \(x_{i}\). The mixture adsorption loading for each component are then calculated using

\[\begin{split}q_{tot} = \frac{1}{\sum_{i}^{N}\,\left(\frac{x_{i}}{q_{i}^{*}(P_{i}^{*})}\right)} \\ q_{mix, i}^{*} = x_{i}\,q_{tot}\end{split}\]

Note

The IAST equations are highly non-linear and often extremely sensitive to the initial guess used. The solution can often be computationally expensive and time-consuming.