Grand canonical ensemble

From SklogWiki

Jump to: navigation, search

The grand-canonical ensemble is particularly well suited to simulation studies of adsorption.

Contents

[edit] Ensemble variables

[edit] Grand canonical partition function

The classical grand canonical partition function for a one-component system in a three-dimensional space is given by:

 Q_{\mu VT} = \sum_{N=0}^{\infty} \frac{ \exp \left[ \beta \mu N \right] V^N}{N! \Lambda^{3N} } \int  d (R^*)^{3N} \exp \left[ - \beta U \left( V, (R^*)^{3N} \right) \right]

where:

  • N is the number of particles
  •  \left. \Lambda \right. is the de Broglie thermal wavelength (which depends on the temperature)
  •  \beta = \frac{1}{k_B T} , with kB being the Boltzmann constant
  • U is the potential energy, which depends on the coordinates of the particles (and on the interaction model)
  •  \left( R^*\right)^{3N} represent the 3N position coordinates of the particles (reduced with the system size): i.e.  \int d (R^*)^{3N} = 1

[edit] Helmholtz energy and partition function

The corresponding thermodynamic potential, the grand potential, Ω, for the aforementioned grand canonical partition function is:

 \Omega = \left. A - \mu N \right. ,

where A is the Helmholtz energy function. Using the relation

\left.U\right.=TS -PV + \mu N

one arrives at

 \left. \Omega \right.= -PV

i.e.:

 \left. p V = k_B T \log Q_{\mu V T } \right.

[edit] See also

[edit] References

Personal tools