Thermodynamic integration: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
mNo edit summary
Line 1: Line 1:
'''Thermodynamic integration''' is used to calculate the difference in the [[Helmholtz energy function]], <math>A</math>, between two states.
'''Thermodynamic integration''' is used to calculate the difference in the [[Helmholtz energy function]], <math>A</math>, between two states.
The path '''must''' be ''continuous'' and ''reversible'' (Ref. 1 Eq. 3.5)
The path '''must''' be ''continuous'' and ''reversible'', i.e., the system must evolve through a succession of equilibrium states (Ref. 1 Eq. 3.5)


:<math>\Delta A = A(\lambda) - A(\lambda_0) = \int_{\lambda_0}^{\lambda}  \left\langle \frac{\partial U(\mathbf{r},\lambda)}{\partial \lambda} \right\rangle_{\lambda} ~\mathrm{d}\lambda</math>
:<math>\Delta A = A(\lambda) - A(\lambda_0) = \int_{\lambda_0}^{\lambda}  \left\langle \frac{\partial U(\mathbf{r},\lambda)}{\partial \lambda} \right\rangle_{\lambda} ~\mathrm{d}\lambda</math>

Revision as of 11:21, 5 July 2011

Thermodynamic integration is used to calculate the difference in the Helmholtz energy function, 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 A} , between two states. The path must be continuous and reversible, i.e., the system must evolve through a succession of equilibrium states (Ref. 1 Eq. 3.5)

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 \Delta A = A(\lambda) - A(\lambda_0) = \int_{\lambda_0}^{\lambda} \left\langle \frac{\partial U(\mathbf{r},\lambda)}{\partial \lambda} \right\rangle_{\lambda} ~\mathrm{d}\lambda}

Isothermal integration

At constant temperature (Ref. 2 Eq. 5):

Isobaric integration

At constant pressure (Ref. 2 Eq. 6):

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 \frac{G(T_2,p)}{Nk_BT_2} = \frac{G(T_1,p)}{Nk_BT_1} - \int_{T_1}^{T_2} \frac{H(T)}{Nk_BT^2} ~\mathrm{d}T }

where 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 G} is the Gibbs energy function 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 H} is the enthalpy.

Isochoric integration

At constant volume (Ref. 2 Eq. 7):

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 \frac{A(T_2,V)}{Nk_BT_2} = \frac{A(T_1,V)}{Nk_BT_1} - \int_{T_1}^{T_2} \frac{U(T)}{Nk_BT^2} ~\mathrm{d}T }

where is the internal energy.

See also

References

  1. J. A. Barker and D. Henderson "What is "liquid"? Understanding the states of matter ", Reviews of Modern Physics 48 pp. 587 - 671 (1976)
  2. C. Vega, E. Sanz, J. L. F. Abascal and E. G. Noya "Determination of phase diagrams via computer simulation: methodology and applications to water, electrolytes and proteins", Journal of Physics: Condensed Matter 20 153101 (2008) (section 4)

Related reading