Editing Gibbs-Duhem integration
Jump to navigation
Jump to search
The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.
Latest revision | Your text | ||
Line 1: | Line 1: | ||
The so-called | CURRENTLY THIS ARTICLE IS UNDER CONSTRUCTION | ||
molecular | == History == | ||
The so-called Gibbs-Duhem Integration referes to a number of methods that couple | |||
molecular simulation techniques with thermodynamic equations in order to draw | |||
phase coexistence lines. | |||
The method was proposed by Kofke (Ref 1-2). | |||
== Basic Features == | == Basic Features == | ||
Consider two thermodynamic phases: <math> a </math> and <math> b </math>, at thermodynamic equilibrium at certain conditions. | Consider two thermodynamic phases: <math> a </math> and <math> b </math>, at thermodynamic equilibrium at certain conditions. | ||
The thermodynamic equilibrium implies: | |||
* Equal | * Equal temperature in both phases: <math> T = T_{a} = T_{b} </math>, i.e. thermal equilbirum. | ||
* Equal | * Equal pressure in both phases <math> p = p_{a} = p_{b} </math>, i.e. mechanical equilbrium. | ||
* Equal | * Equal chemical potentials for the components <math> \mu_i = \mu_{ia} = \mu_{ib} </math>, i.e. ''material'' equilibrium. | ||
In addition | In addition if we are dealing with a statistical mechanics model, with certain parameters that we can represent as <math> \lambda </math> , the | ||
model should be the same in both phases. | model should be the same in both phases. | ||
== Example: phase equilibria of one- | == Example: phase equilibria of one-compoment system == | ||
Notice: The derivation that follows is just a particular route to perform the integration | Notice: The derivation that follows is just a particular route to perform the integration | ||
Line 27: | Line 31: | ||
where | where | ||
* <math> \beta | * <math> \beta = 1/k_B T </math>, where <math> k_B </math> is the [[Boltzmann constant]] | ||
When a differential change of the conditions is performed | When a differential change of the conditions is performed we wil have for any phase: | ||
: <math> d \left( \beta\mu \right) = \left[ \frac{ \partial (\beta \mu) }{\partial \beta} \right]_{\beta p,\lambda} d \beta + | : <math> d \left( \beta\mu \right) = \left[ \frac{ \partial (\beta \mu) }{\partial \beta} \right]_{\beta p,\lambda} d \beta + | ||
Line 35: | Line 39: | ||
</math> | </math> | ||
Taking into account that <math> \mu </math> is the [[Gibbs energy function]] per particle | Taking into account that <math> \mu </math> is the [[Gibbs energy function|Gibbs free energy]] per particle | ||
: <math> d \left( \beta\mu \right) = \frac{E}{N} d \beta + \frac{ V }{N } d (\beta p) + | : <math> d \left( \beta\mu \right) = \frac{E}{N} d \beta + \frac{ V }{N } d (\beta p) + | ||
\left[ \frac{ \partial (\beta \mu) }{\partial \lambda} \right]_{\beta,\beta p} d \lambda. | \left[ \frac{ \partial (\beta \mu) }{\partial \lambda} \right]_{\beta,\beta p} d \lambda. | ||
</math> | </math> | ||
Let us use a bar to design quantities divided by the number of particles: e.g. <math> \bar{E} = E/N; \bar{V} = V/N </math>; | Let us use a bar to design quantities divided by the number of particles: e.g. <math> \bar{E} = E/N; \bar{V} = V/N </math>; | ||
and taking into account the definition: | and taking into account the definition: | ||
: <math> \bar{L} \equiv \left[ \frac {\partial (\beta \mu )}{\partial \lambda }\right]_{\beta,\beta p} </math> | : <math> \bar{L} \equiv \frac{1}{N} \left[ \frac {\partial (\beta \mu )}{\partial \lambda }\right]_{\beta,\beta p} </math> | ||
Again, let us suppose that we have a phase coexistence at a point given by <math>\left[ \beta_0, (\beta p)_0, \lambda_0 \right]</math> and that | Again, let us suppose that we have a phase coexistence at a point given by <math>\left[ \beta_0, (\beta p)_0, \lambda_0 \right]</math> and that | ||
Line 63: | Line 57: | ||
constrained to fulfill: | constrained to fulfill: | ||
<math> \left( \Delta \bar{E} \right) d \beta + \left( \Delta \bar{V} \right) d (\beta p) + \left(\Delta \bar{L} \right) d \lambda = 0 </math> | |||
== References == | == References == | ||
#[http://dx.doi.org/10.1080/00268979300100881 David A. Kofke, Gibbs-Duhem integration: a new method for direct evaluation of phase coexistence by molecular simulation, Mol. Phys. '''78''' , pp 1331 - 1336 (1993)] | |||
#[http://dx.doi.org/10.1063/1.465023 David A. Kofke, Direct evaluation of phase coexistence by molecular simulation via integration along the saturation line, J. Chem. Phys. '''98''' ,pp. 4149-4162 (1993) ] | |||