Math 2.0

Statistical mechanics, structure of matter Example

Hypernetted-chain equation

Example

Consider a classical fluid in ℝ3\R^3 with number density ρ\rho, inverse temperature β\beta and pair potential u(r)u(r). Let gg be the radial distribution function, h=g−1h = g - 1, and cc the direct correlation function of the Ornstein–Zernike equation. The hypernetted-chain closure is g(r)=exp⁡(−βu(r)+h(r)−c(r)),equivalentlyc(r)=h(r)−ln⁡g(r)−βu(r).g(r) = \exp\big(-\beta u(r) + h(r) - c(r)\big), \quad \text{equivalently} \quad c(r) = h(r) - \ln g(r) - \beta u(r). Together with the Ornstein–Zernike equation h(r12)=c(r12)+ρ∫ℝ3c(r13) h(r23) dx3,h(r_{12}) = c(r_{12}) + \rho \int_{\R^3} c(r_{13}) \, h(r_{23}) \, dx_3, where rij=∣xi−xj∣r_{ij} = \abs{x_i - x_j}, it forms a closed nonlinear integral equation for hh: substitute c=h−ln⁡(1+h)−βuc = h - \ln(1 + h) - \beta u.

Check at low density: as ρ→0\rho \to 0 the Ornstein–Zernike equation gives c=hc = h to leading order. Then the closure gives g=e−βug = e^{-\beta u}, the exact low-density limit of the radial distribution function.

Remarks

The closure is usually solved numerically by iteration, with the convolution done by the fast Fourier transform. The Percus–Yevick approximation is usually more accurate for hard spheres, while the hypernetted-chain equation is usually preferred for long-range potentials such as the Coulomb potential.

Further reading

  1. Web https://en.wikipedia.org/wiki/Hypernetted-chain_equation
About this page
TagMRJ2FHX6K1
URIm2://MRJ2FHX6K1
Kindexample
WritingW0 stub
MSC82D15, 82B21, 45G10
Instance ofNonlinear integral equation
Same as

Markdown · JSON · Edit on GitHub

What this page needs for W1

  • an inventory row: the row of the inventory that the node accounts for (0 of 1)
  • a source with a locator (0 of 1)

The writing levels say how complete the text of a page is. They do not measure its truth.

What links here

No page links to this page yet.