9-3 Lennard-Jones potential

[EN CONSTRUCCION]

Functional form

The 9-3 Lennard-Jones potential is related to the standard Lennard-Jones potential.

It takes the form:

$V(r)=\frac{3\sqrt{3}}{2}ϵ[{(\frac{\sigma }{r})}^9-{(\frac{\sigma }{r})}^3].$

The minimum value of $V(r)$ is obtained at $r=r_{min}$, with

• $V(r_{min})=-ϵ$,

• $\frac{r_{min}}{\sigma }=3^{1/6}$

Applications

It is commonly used to model the interaction between the particles of a fluid with a flat structureless solid wall.

Interaction between a solid and a fluid molecule

Let us consider the space divided in two regions:

• $x<0$: this region is occupied by a diffuse solid with density ${\rho }_s$ composed of 12-6 Lennard-Jones atoms

with paremeters ${\sigma }_s$ and ${ϵ}_a$

Our aim is to compute the total interaction between this solid and a molecule located at a position $x_f>0$. Such an interaction can be computed using cylindrical coordinates ( I GUESS SO, at least).

The interaction will be:

$V_W(x)=4{ϵ}_{sf}{\rho }_s\underset{0}{\overset{2\pi }{\int }}d\varphi \underset{-\mathrm{\infty }}{\overset{-x}{\int }}dz\underset{0}{\overset{\mathrm{\infty }}{\int }}\text{dr}[{\sigma }^{12}\frac{r}{(r^2+z^2{)}^6}-{\sigma }^6\frac{r}{(r^2+z^2{)}^3}].$

$V_W(x)=8\pi {ϵ}_{sf}{\rho }_s\underset{-\mathrm{\infty }}{\overset{-x}{\int }}\text{d}z{[\frac{{\sigma }^{12}}{10(r^2+z^2{)}^5}-\frac{{\sigma }^6}{4(r^2+z^2{)}^2}]}_{r=\mathrm{\infty }}^{r=0}.$

$V_W(x)=8\pi {ϵ}_{sf}{\rho }_s\underset{-\mathrm{\infty }}{\overset{-x}{\int }}\text{d}z[\frac{{\sigma }^{12}}{10z^{10}}-\frac{{\sigma }^6}{4z^4}];$

$V_W(x)=8\pi {ϵ}_{sf}{\rho }_s{[-\frac{{\sigma }^{12}}{90z^9}+\frac{{\sigma }^6}{12z^3}]}_{z=-\mathrm{\infty }}^{z=-x};$

$V_W(x)=\frac{4\pi {ϵ}_{sf}{\rho }_s{\sigma }^3}{3}[\frac{{\sigma }^9}{15x^9}-\frac{{\sigma }^3}{2x^3}]$

[TO BE CONTINUED]