Ornstein-Zernike relation: Difference between revisions
Carl McBride (talk | contribs) mNo edit summary |
Carl McBride (talk | contribs) mNo edit summary |
||
| Line 8: | Line 8: | ||
*<math>B(r)</math> is the [[Closures | bridge]] function. | *<math>B(r)</math> is the [[Closures | bridge]] function. | ||
*<math>\omega(r)</math> is the [[thermal potential]], <math>\omega(r) \equiv \gamma(r) + B(r)</math>. | *<math>\omega(r)</math> is the [[thermal potential]], <math>\omega(r) \equiv \gamma(r) + B(r)</math>. | ||
*<math>f(r)</math> is the [[Mayer <math>f</math>-function | *<math>f(r)</math> is the [[Mayer]] <math>f</math>-function, defined as <math>f(r) \equiv e^{-\beta \Phi(r)} -1</math>. | ||
| Line 19: | Line 19: | ||
Because of this <math>h</math> must be determined [[self-consistently]]. | Because of this <math>h</math> must be determined [[self-consistently]]. | ||
This need for self-consistency is characteristic of all many-body problems. | This need for self-consistency is characteristic of all many-body problems. | ||
(Hansen | (Hansen and McDonald, section 5.2 p. 106) For a system in an external field, the OZ has the form (5.2.7) | ||
:<math>h(1,2) = c(1,2) + \int \rho^{(1)}(3) c(1,3)h(3,2) d3</math> | :<math>h(1,2) = c(1,2) + \int \rho^{(1)}(3) c(1,3)h(3,2) d3</math> | ||
If the system is both homogeneous and isotropic, the OZ relation becomes (Ref. 1Eq. 6) | If the system is both homogeneous and isotropic, the OZ relation becomes (Ref. 1Eq. 6) | ||
<math>\gamma (r) \equiv h(r) - c(r) = \rho \int h(r')~c(|r - r'|) dr'</math> | <math>\gamma (r) \equiv h(r) - c(r) = \rho \int h(r')~c(|r - r'|) dr'</math> | ||
In words, this equation (Hansen | |||
In words, this equation (Hansen and McDonald, section 5.2 p. 107) | |||
``...describes the fact that the ''total'' correlation between particles 1 and 2, represented by <math>h(1,2)</math>, | ``...describes the fact that the ''total'' correlation between particles 1 and 2, represented by <math>h(1,2)</math>, | ||
is due in part to the ''direct'' correlation between 1 and 2, represented by <math>c(1,2)</math>, but also to the ''indirect'' correlation, | is due in part to the ''direct'' correlation between 1 and 2, represented by <math>c(1,2)</math>, but also to the ''indirect'' correlation, | ||
| Line 73: | Line 74: | ||
==References== | ==References== | ||
#[ | #[L. S. Ornstein and F. Zernike "Accidental deviations of density and opalescence at the critical point of a single substance", Koninklijke Nederlandse Akademie van Wetenschappen Amsterdam Proc. Sec. Sci. '''17''' pp. 793- (1914)] | ||
#[ | #[http://dx.doi.org/10.1103/PhysRevA.45.816 James A. Given "Liquid-state methods for random media: Random sequential adsorption", Physical Review A '''45''' pp. 816 - 824 (1992)] | ||
#[ | #[http://dx.doi.org/10.1063/1.470724 Der-Ming Duh and A. D. J. Haymet "Integral equation theory for uncharged liquids: The Lennard-Jones fluid and the bridge function", Journal of Chemical Physics '''103''' pp. 2625-2633 (1995)] | ||
#[Hansen and MacDonald "Theory of Simple Liquids"] | |||
Revision as of 15:59, 20 February 2007
Notation:
- 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(r)} is the pair distribution 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 \Phi(r)} is the pair potential acting between pairs.
- is the total correlation 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 h(1,2) \equiv g(r) -1} .
- 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 c(1,2)} is the direct correlation 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 \gamma (r)} is the indirect (or series or chain) correlation 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 \gamma (r) \equiv h(r) - c(r)} .
- 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 y(r_{12})} is the cavity correlation functionFailed 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 y(r) \equiv g(r) /e^{-\beta \Phi(r)}}
- 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 B(r)} is the bridge 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 \omega(r)} is the thermal 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 \omega(r) \equiv \gamma(r) + B(r)} .
- 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 f(r)} is the Mayer Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle f} -function, defined as 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 f(r) \equiv e^{-\beta \Phi(r)} -1} .
The Ornstein-Zernike relation (OZ) integral equation is
- 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=h\left[c\right]}
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 h[c]} denotes a functional of 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 c} . This relation is exact. This is complemented by the closure relation
- 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 c=c\left[h\right]}
Note that 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} depends on 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 c} , 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 c} depends on 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} . Because of this 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} must be determined self-consistently. This need for self-consistency is characteristic of all many-body problems. (Hansen and McDonald, section 5.2 p. 106) For a system in an external field, the OZ has the form (5.2.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 h(1,2) = c(1,2) + \int \rho^{(1)}(3) c(1,3)h(3,2) d3}
If the system is both homogeneous and isotropic, the OZ relation becomes (Ref. 1Eq. 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 \gamma (r) \equiv h(r) - c(r) = \rho \int h(r')~c(|r - r'|) dr'}
In words, this equation (Hansen and McDonald, section 5.2 p. 107)
``...describes the fact that the total correlation between particles 1 and 2, represented by 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(1,2)}
,
is due in part to the direct correlation between 1 and 2, represented by 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 c(1,2)}
, but also to the indirect correlation,
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 \gamma (r)}
, propagated via increasingly large numbers of intermediate particles."
Notice that this equation is basically a convolution, i.e.
- 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 \equiv c + \rho h\otimes c }
(Note: the convolution operation written here as 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 \otimes} is more frequently written as 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 *} ) This can be seen by expanding the integral in terms of 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(r)} (here truncated at the fourth iteration):
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(r) = c(r) + \rho \int c(|r - r'|) c(r') dr' + \rho^2 \int \int c(|r - r'|) c(|r' - r''|) c(r'') dr''dr' + \rho^3 \int\int\int c(|r - r'|) c(|r' - r''|) c(|r'' - r'''|) c(r''') dr'''dr''dr' + \rho^4 \int \int\int\int c(|r - r'|) c(|r' - r''|) c(|r'' - r'''|) c(|r''' - r''''|) h(r'''') dr'''' dr'''dr''dr'}
etc. Diagrammatically this expression can be written as (Ref. 2):
where the bold lines connecting root points denote 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 c} functions, the blobs denote 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} functions. An arrow pointing from left to right indicates an uphill path from one root point to another. An `uphill path' is a sequence of Mayer bonds passing through increasing particle labels. The OZ relation can be derived by performing a functional differentiation of the grand canonical distribution function (HM check this).
OZ equation in Fourier space
The OZ equation may be written in Fourier space as (Eq. 5 in Ref. 3):
- 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 \hat{\gamma} = (I - \rho \hat{c})^{-1} \hat{c} \rho \hat{c}}
The carets denote the three-dimensional Fourier transformed quantities which reduce explicitly to:
- 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 \hat{\gamma} (k) = \frac{4 \pi}{k} \int_0^\infty r~\sin (kr) \gamma(r) dr}
- 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 \gamma (r) = \frac{1}{2 \pi^2 r} \int_0^\infty k~\sin (kr) \hat{\gamma}(r) dk}
Note:
- 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 \hat{h}(0) = \int h(r) dr}
- 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 \hat{c}(0) = \int c(r) dr}
References
- [L. S. Ornstein and F. Zernike "Accidental deviations of density and opalescence at the critical point of a single substance", Koninklijke Nederlandse Akademie van Wetenschappen Amsterdam Proc. Sec. Sci. 17 pp. 793- (1914)]
- James A. Given "Liquid-state methods for random media: Random sequential adsorption", Physical Review A 45 pp. 816 - 824 (1992)
- Der-Ming Duh and A. D. J. Haymet "Integral equation theory for uncharged liquids: The Lennard-Jones fluid and the bridge function", Journal of Chemical Physics 103 pp. 2625-2633 (1995)
- [Hansen and MacDonald "Theory of Simple Liquids"]
