Lennard-Jones model: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
Line 36: Line 36:
: <math> \frac{r_{min}}{\sigma} = 2^{1/6} \simeq  1.12246 ...  </math>
: <math> \frac{r_{min}}{\sigma} = 2^{1/6} \simeq  1.12246 ...  </math>


== Approximations in simulation: truncation and shifting ==


== Related potential models ==
== Related potential models ==

Revision as of 13:28, 23 March 2007

Lennard-Jones potential

The Lennard-Jones potential is given by

V(r)=4ϵ[(σr)12−(σr)6]

where:

  • V(r) : potential energy of interaction between two particles at a distance r;
  • σ : diameter (length);
  • ϵ : well depth (energy)

Reduced units:

  • Density, ρ*≡ρσ3, where ρ=N/V (number of particles N divided by the volume V.)
  • Temperature; T*≡kBT/ϵ, where T is the absolute temperature and kB is the Boltzmann constant

Argon

The Lennard-Jones parameters for argon are ϵ/kB≈ 119.8 K and σ≈ 0.3405 nm. (Ref. ?)

This figure was produced using gnuplot with the command:

plot (4*120*((0.34/x)**12-(0.34/x)**6))

Features

Special points:

  • V(σ)=0
  • Minimum value of V(r) at r=rmin;
rminσ=21/6≃1.12246...

Approximations in simulation: truncation and shifting

Related potential models

It is relatively common the use of potential functions given by:

V(r)=cm,nϵ[(σr)m−(σr)n].

with m and n being positive integer numbers and m>n, and cm,n is chosen to get the minumum value of V(r) being Vmin=−ϵ

These forms are usually refered to as m-n Lennard-Jones Potential.

The 9-3 Lennard-Jones interaction potential is often use to model the interaction between the atoms/molecules of a fluid and a continuous solid wall. In (9-3 Lennard-Jones potential) a justification of this use is presented.

References

  1. J. E. Lennard-Jones, "Cohesion", Proceedings of the Physical Society, 43 pp. 461-482 (1931)