Computation of phase equilibria: Difference between revisions
mNo edit summary |
mNo edit summary |
||
| Line 1: | Line 1: | ||
[CURRENTLY WORKING ON THE PAGE] | |||
The computation of phase equilibria using computer simulation can follow different strategies. | The computation of phase equilibria using computer simulation can follow different strategies. | ||
== Liquid-vapor equilibria of one component systems == | == Liquid-vapor equilibria of one component systems == | ||
The thermodynamic equilibrium implies, for two phases <math> \alpha </math> and <math> \beta </math>: | The thermodynamic equilibrium implies, for two phases <math> \alpha </math> and <math> \beta </math>: | ||
* Equal temperature: <math> T_{\alpha} = T_{\beta} </math> | |||
* Equal pressure: <math> p_{\alpha} = p_{\beta} </math> | |||
* Equal chemical potential: <math> \mu_{\alpha} = \mu_{\beta} </math> | |||
=== Independent simulations for each phase at fixed <math> T </math> in the [[canonical ensemble]] === | |||
The simulations can be carried out either using [[Monte Carlo]] or [[Molecular dynamics]] techhniques. | |||
Let us assume that we have some knowledge on the phase diagram of the system. We could: | |||
- Fix a temperature | |||
- Perform a few simulations in the low density region (where the gas phase density is expected to be) | |||
- Perform a few simulations in the moderate / high density regions (where the liquid phase should appear) | |||
- In these simulations we can compute for each density (at fixed T) the values of the pressure and the | |||
chemical potentials (for instance using the [[Widom test-particle method]]) | |||
Revision as of 17:46, 21 September 2007
[CURRENTLY WORKING ON THE PAGE]
The computation of phase equilibria using computer simulation can follow different strategies.
Liquid-vapor equilibria of one component systems
The thermodynamic equilibrium implies, for two phases Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha } and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta } :
- Equal temperature: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T_{\alpha} = T_{\beta} }
- Equal pressure: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p_{\alpha} = p_{\beta} }
- Equal chemical potential: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu_{\alpha} = \mu_{\beta} }
Independent simulations for each phase at fixed Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T } in the canonical ensemble
The simulations can be carried out either using Monte Carlo or Molecular dynamics techhniques. Let us assume that we have some knowledge on the phase diagram of the system. We could:
- Fix a temperature
- Perform a few simulations in the low density region (where the gas phase density is expected to be)
- Perform a few simulations in the moderate / high density regions (where the liquid phase should appear)
- In these simulations we can compute for each density (at fixed T) the values of the pressure and the chemical potentials (for instance using the Widom test-particle method)