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 '''Gibbs-Duhem integration''' refers to a number of methods that couple | The so-called '''Gibbs-Duhem integration''' refers to a number of methods that couple | ||
molecular [[Computer simulation techniques |simulation techniques]] with [[Thermodynamic relations |thermodynamic equations]] in order to draw | molecular [[Computer simulation techniques |simulation techniques]] with [[Thermodynamic relations |thermodynamic equations]] in order to draw | ||
[[Computation of phase equilibria | phase coexistence]] lines. The original method was proposed by David Kofke | [[Computation of phase equilibria | phase coexistence]] lines. The original method was proposed by David Kofke (Refs. 1 and 2). | ||
== Basic Features == | == Basic Features == | ||
Line 28: | Line 27: | ||
where | where | ||
* <math> \beta := 1/k_B T </math>, where <math> k_B </math> is the [[Boltzmann constant]] | * <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 one will have | When a differential change of the conditions is performed one will, 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 73: | Line 72: | ||
The Gibbs-Duhem integration technique, for this example, will be a numerical procedure covering the following tasks: | The Gibbs-Duhem integration technique, for this example, will be a numerical procedure covering the following tasks: | ||
* Computer simulation (for instance using [[Metropolis Monte Carlo]] in the | * Computer simulation (for instance using [[Metropolis Monte Carlo]] in the NpT ensemble) runs to estimate the values of <math> \bar{L}, \bar{V} </math> for both | ||
phases at given values of <math> [\beta, \beta p, \lambda ] </math>. | phases at given values of <math> [\beta, \beta p, \lambda ] </math>. | ||
Line 80: | Line 79: | ||
== Peculiarities of the method (Warnings) == | == Peculiarities of the method (Warnings) == | ||
* A good initial point must be known to start the procedure (See | * A good initial point must be known to start the procedure (See Ref. 3 and the entry: [[computation of phase equilibria]]). | ||
* The ''integrand'' of the differential equation is computed with some numerical uncertainty | * The ''integrand'' of the differential equation is computed with some numerical uncertainty | ||
Line 87: | Line 86: | ||
== 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", Molecular Physics '''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", Journal of Chemical Physics '''98''' pp. 4149-4162 (1993)] | ||
#[http://dx.doi.org/10.1063/1.2137705 A. van 't Hof, S. W. de Leeuw, and C. J. Peters "Computing the starting state for Gibbs-Duhem integration", Journal of Chemical Physics '''124''' 054905 (2006)] | |||
#[http://dx.doi.org/10.1063/1.2137706 A. van 't Hof, C. J. Peters, and S. W. de Leeuw "An advanced Gibbs-Duhem integration method: Theory and applications", Journal of Chemical Physics '''124''' 054906 (2006)] | |||
[[category: computer simulation techniques]] | [[category: computer simulation techniques]] |