Verlet modified: Difference between revisions
Carl McBride (talk | contribs) mNo edit summary |
No edit summary |
||
| Line 2: | Line 2: | ||
in terms of the [[cavity correlation function]], is (Eq. 3) | in terms of the [[cavity correlation function]], is (Eq. 3) | ||
:<math>y(r) = \gamma (r) - A \frac{1}{2} \gamma^2(r) \left[ \frac{1}{1+ B \gamma(r) /2} \right]</math> | :<math>\ln y(r) = \gamma (r) - A \frac{1}{2} \gamma^2(r) \left[ \frac{1}{1+ B \gamma(r) /2} \right]</math> | ||
where several sets of values are tried for ''A'' and ''B'' (Note, when ''A=0'' the [[HNC| hyper-netted chain]] is recovered). | where several sets of values are tried for ''A'' and ''B'' (Note, when ''A=0'' the [[HNC| hyper-netted chain]] is recovered). | ||
| Line 8: | Line 8: | ||
by minimising the difference between the pressures obtained via the [[Pressure equation | virial]] and [[Compressibility equation | compressibility]] routes: | by minimising the difference between the pressures obtained via the [[Pressure equation | virial]] and [[Compressibility equation | compressibility]] routes: | ||
:<math>y(r) = \gamma (r) - A \frac{1}{2} \gamma^2(r) \left[ \frac{1+ \lambda \gamma(r)}{1+ \mu \gamma(r)} \right]</math> | :<math>\ln y(r) = \gamma (r) - A \frac{1}{2} \gamma^2(r) \left[ \frac{1+ \lambda \gamma(r)}{1+ \mu \gamma(r)} \right]</math> | ||
with <math>A= 0.80</math>, <math>\lambda= 0.03496</math> and <math>\mu = 0.6586</math>. | with <math>A= 0.80</math>, <math>\lambda= 0.03496</math> and <math>\mu = 0.6586</math>. | ||
Revision as of 18:00, 1 September 2015
The Verlet modified (1980) (Ref. 1) closure relation for hard sphere fluids, in terms of the cavity correlation function, is (Eq. 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 \ln y(r) = \gamma (r) - A \frac{1}{2} \gamma^2(r) \left[ \frac{1}{1+ B \gamma(r) /2} \right]}
where several sets of values are tried for A and B (Note, when A=0 the hyper-netted chain is recovered). Later (Ref. 2) Verlet used a Padé (2/1) approximant (Eq. 6) fitted to obtain the best hard sphere results by minimising the difference between the pressures obtained via the virial and compressibility routes:
- 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 \ln y(r) = \gamma (r) - A \frac{1}{2} \gamma^2(r) \left[ \frac{1+ \lambda \gamma(r)}{1+ \mu \gamma(r)} \right]}
with 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= 0.80} , 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 \lambda= 0.03496} 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 \mu = 0.6586} .