Composition Derivatives of Activity Coefficients

The derivation below is based on this resource from the people who created the CocoSimulator and this older paper from the same authors.

To write the derivative of the activity coefficient w.r.t to mole fraction, we’ll work from the Gibbs excess energy to find an expression for the activity coefficient in terms of derivatives of Gibbs excess energy.

Therefore, the activity coefficient is:

where $n_{k \neq i}$ indicates the unconstrained derivative where all components are fixed except $i$

• Using chain rule to expand:

• We can expand the second term again using chain rule
$$\frac{\partial Q}{\partial n_i} = \sum_{j=1}^n Q_j \frac{\partial x_j}{\partial n_i} |_{n_{k \neq i}}$$$$Q_j = \frac{\partial Q}{\partial x_j}|_{x_{k\neq j}}$$
1. We can easily find $Q_j$, but how do we get the derivative of the mole fraction w.r.t moles?
• Mole fraction derivative
• We know from the differential for partial moles:
$$dn_i = d(n_tx_i) = n_t dx_i + x_i dn_t$$$$dn_i - x_i \sum_j dn_j = n_t dx_i$$$$dx_i = \sum_j \frac{\delta_{ij}-x_i}{n_t} dn_j$$
• $\delta_{ij}$ is the kronecker delta, so it is 1 when $i=j$ and 0 everywhere else
• We also know from the total differential
$$dx_i = \sum_j \biggl(\frac{\partial x_i}{\partial n_j} \biggr) dn_j$$
• By inspection of the last two equations
$$\biggl (\frac{\partial x_i}{\partial n_j} \biggr )_{n_{k \neq j}} = \frac{1}{n_t}(\delta_{ij} - x_i)$$
• Putting all the expansions together:
$$\ln \gamma_i = Q + n_t \frac{\partial Q}{\partial n_i}|_{T,P,n_{k\neq i}}$$$$\ln \gamma_i = Q + n_t\biggl(\sum_{k=1}^n Q_k \frac{\partial x_k}{\partial n_i} |_{n_{l \neq i}} \biggr)$$$$\ln \gamma_i = Q + n_t\biggl(\sum_{k=1}^n Q_k \frac{1}{n_t}(\delta_{ik} - x_k)\biggr)$$
1. Cancel out $n_t$ and expand the Kronecker delta out of the sum to get our desired expression for the activity coefficient:

• Now we can easily find the derivative of the activity coefficient w.r.t the mole fraction - this is the version that only constrains one mole fraction but doesn’t enforce them all summing to one:
$$\frac{\partial \ln \gamma_i}{\partial x_j}|_{x_{l \neq i}} = Q_{ij} - \sum_{k=1}^n x_kQ_{kj}$$$$Q_{ij} \equiv \frac{\partial Q_i}{\partial x_j} |_{x_k\neq i}$$
• To get the constrained version, we can observe that an increase in the mole fraction of one component will result in an equal and opposite decrease in the mole fraction of another (given all other mole fractions are fixed except one). This then implies the following equations, starting from the total derivative of the activity coefficient:
$$d(\ln \gamma_i(x_1, \dots, x_N)) = \sum_{k=1}^N \frac{\partial \ln \gamma_i}{\partial x_k}|_{x_{j \neq k}}dx_k$$

Only two mole fractions are actually varying so when we take the partial derivative:

$$\frac{\partial}{\partial x_l}(\ln \gamma_i(x_1, \dots, x_N)) = \frac{\partial \ln \gamma_i}{\partial x_l}|_{x_{j \neq l}}+ \frac{\partial \ln \gamma_i}{\partial x_m}|_{x_{j \neq m}} \frac{dx_m}{dx_l}$$

Note that that $l$ can be equal to $i$.

We know by the logic above that $dx_l = -dx_m$, so

$$\frac{\partial \ln \gamma_i}{\partial x_l}|_{x_{j \neq i,l}}= \frac{\partial \ln \gamma_i}{\partial x_l}|_{x_{j \neq l}}- \frac{\partial \ln \gamma_i}{\partial x_m}|_{x_{j \neq m}}$$$$\frac{\partial \ln \gamma_i}{\partial x_l}|_{x_{j \neq i,l}}= Q_{il} - \sum_{k=1}^n x_kQ_{kl} - \biggl(Q_{im} - \sum_{k=1}^n x_kQ_{km} \biggr)$$$$\frac{\partial \ln \gamma_i}{\partial x_l}|_{x_{j \neq i,l}}= Q_{il} - Q_{im} + \sum_{k=1}^n x_k(Q_{kl} -Q_{km})$$

So if we put this in terms of the Gibbs excess energy directly:

From what I understand you can arbitrarily choose $m$. This is the component which will be determined by all the others. In the paper, they choose the $n$th component. From the third page in the PDF: