ReactionEquilibrium:Background:Reaction Equilibrium Reaction Equilibrium
From Wiketomica
| Introduction |
Background
|
| Basic Layout |
| Examples |
| Problems |
| Run Simulation |
| Assessment |
| Credits |
| Etomica Modules |
For the elementary reaction
where X and Y may be R or B, chemical equilibrium requires
μXY − μX − μY = 0
If we separate the chemical potential into various contributions, as described in the Standard Properties page, this is
![\mu^0_{XY}(T) + kT \ln(\rho x_{XY}\phi_{XY}) - [\mu^0_X(T) + kT \ln(\rho x_{X}\phi_{X}) + \mu^0_Y(T) + kT \ln(\rho x_{Y}\phi_{Y})] = 0](/w/images/math/3/1/9/319a8aaf2a7f1e5f993581c7e1f3542b.png)
We can collect terms as follows
![\left (\frac {1}{\rho}\right ) \left (\frac{x_{XY}}{x_Xx_Y}\right ) \left (\frac{\phi_{XY}}{\phi_X\phi_Y}\right ) = \exp[-(\mu^0_{XY} - \mu^0_X - \mu^0_Y) / kT]](/w/images/math/b/4/2/b42e2f05a06ddfebc3fd4415988d1339.png)
We define the equilibrium constant K(T) by the right-hand side of this equation; on this Standard properties page the standard-state chemical potentials are given in terms of the potential parameters, so K(T) is known. Then

The fugacity coefficients φi are not known in general. For sufficiently low density they approach unity (1.0). It may also happen that the ratio of the φi is approximately 1.0 even when they are not individually. For our calculations we will make this approximation.
We can write an equilibrium equation for each of the elementary reactions listed on the System page. This gives us three equations to solve for the six unknowns in this system: ρ,xR,xB,xRR,xRB,xBB. Two of the other three equations express the fact that the total number of R atoms and B atoms in the system cannot change
![\rho\left (x_R + 2x_{RR} + x_{RB}\right ) = \left [\rho(x_R + 2x_{RR} + x_{RB})\right ]_{initial}](/w/images/math/3/3/b/33bacb96d2672869792997afde241a17.png)
![\rho\left (x_B + 2x_{BB} + x_{RB}\right ) = \left [\rho(x_B + 2x_{BB} + x_{RB})\right ]_{initial}](/w/images/math/5/d/c/5dc612ed66541b1e938ce060c016ece8.png)
and the other is simply the normalization of the mole fractions
xR + xB + xRR + xRB + xBB = 1
