Pair distribution function: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
mNo edit summary
mNo edit summary
Line 8: Line 8:
See Eq. 5.10 of Ref. 1:
See Eq. 5.10 of Ref. 1:


:<math>\ln g(r_{12}) + \frac{\Phi(r_{12})}{kT} - E(r_{12}) = n \int \left(g(r_{13}) -1 - \ln g(r_{13}) -  \frac{\Phi(r_{13})}{kT} - E(r_{13})  \right)(g(r_{23}) -1)  ~{\rm d}r_3</math>
:<math>\ln g(r_{12}) + \frac{\Phi(r_{12})}{kT} - E(r_{12}) = n \int \left(g(r_{13}) -1 - \ln g(r_{13}) -  \frac{\Phi(r_{13})}{kT} - E(r_{13})  \right)(g(r_{23}) -1)  ~{\rm d}{\mathbf r}_3</math>


where <math>r_{12} = |{\mathbf r}_2 - {\mathbf r}_1|</math>.
where, ''i.e.'' <math>r_{12} = |{\mathbf r}_2 - {\mathbf r}_1|</math>.
==See also==
==See also==
*[[Radial distribution function]]
*[[Radial distribution function]]

Revision as of 17:09, 10 July 2007

For a fluid of N particles, enclosed in a volume V at a given temperature T (canonical ensemble) interacting via the `central' intermolecular pair potential Φ(r), the two particle distribution function is defined as

gN(2)(r1,r2)=V2∫...∫e−βΦ(r1,...,rN)dr3...drN∫e−βΦ(r1,...,rN)dr1...drN

where β=1/(kBT), where kB is the Boltzmann constant.

Exact convolution equation for g(r)

See Eq. 5.10 of Ref. 1:

lng(r12)+Φ(r12)kT−E(r12)=n∫(g(r13)−1−lng(r13)−Φ(r13)kT−E(r13))(g(r23)−1)dr3

where, i.e. r12=|r2−r1|.

See also

References

  1. J. S. Rowlinson "The equation of state of dense systems", Reports on Progress in Physics 28 pp. 169-199 (1965)