Draw schematically the radial distribution function $g(r)$ for a Lennard-Jones fluid at low and high particle densities. Discuss in both cases the behavior at small $r$. At high densities $g(r)$ has some damped oscillations. Explain their origin and what would happen at long distances.

**What I know:**

The radial distribution function $g(r)$ describes how density varies as a function of distance from a reference particle (*Wikipedia*).

Specifically in Lennard-Jones case, we have (Pasteboard still does not work for me):

Here we can observe that if repulsion outweighs attraction, the curve will die out on the x axis. On the contrary, if it is the other way around, the curve will grow exponentially. For the sake of clarity, I would remark y = 0 as the Ideal Gas behaviour.

I have read that $g(r)$ vanishes at short distances, because the probability of finding two particles close to each other vanishes due to the repulsive part of the potential. At high densities $g(r)$ can show some damped oscillations: https://imgur.com/a/i32D39I.

These oscillations express the preference of the particles to be found at specific distances from a reference particle at the origin. For instance in the LJ case, a first layer of particles will be localized closed to the minimum of LJ's potential, which is the origin of the first peak in $g(r)$.

This layer prevents other particles from getting close to it, which is what causes the first

minimum in $g(r)$.

**What I do not know.**

1)The explanation: 'These oscillations express the preference of the particles to be found at specific distances from a reference particle at the origin'. **How can particles have 'preference' to be at different distances? I think this is not a good physical explanation. I have been thinking about this and I would say that as there is a high density of particles at small distances, the repulsion force is exerted on these, which triggers such oscillations. Finally there is a point where they are so close to each other that the repulsion force provides them with an initial kinetic energy which will be equal to a final electrostatic potential energy of zero (selecting our zero of electrostatic potential energy at y=0) , and that is the moment when the curve dies out on the x axis. Do you agree with this explanation?**

2)Are we dealing with underdamped oscillations at high densities? Could you give an insight into these kind of oscillations (if what I have just said is not enough or incorrect)?

3)When we are dealing with long distances, would $g(r)$ tend to be one? What would we observe, the Ideal Gas behaviour? Why? Is it because there is no interaction between the particles?

As always, I am interested in explaining such phenomena from a physical point of view.

1) The word "preference" suggests conscious choice. Of course the particles do not possess consciousness. What is meant is that the spatial distribution of particles settles into a configuration which minimises potential energy, and this is achieved when there is some structure. For example, **hexagonal close packing**.

Gravitational energy is far too small on the scale of microscopic particles. The balance is between kinetic energy and electrostatic potential energy as defined by the Lennard-Jones potential.

Ideal gas particles are assumed to have no attraction at all, and only "hard sphere" repulsion. The potential function for ideal gas interactions is therefore $U(r)=0$ for $r>2r_0$ and $U(r)=\infty$ for $r\le 2r_0$, where $r_0$ is the radius of particles.

2) I think you have misunderstood the oscillations. These are not **dynamic** oscillations in the motion of particles, like a mass on a spring. So the question about over/under-damping is not relevant - although it is possible some analogy could be made with dynamic oscillations. These are **static** oscillations in the distribution function, as in the ripples of a wave function in quantum mechanics, or the ripples in the distribution of sand on a beach. (For oscillations in radial wave function, see Figure 11 on page 19 here : http://www.umich.edu/~chem461/QMChap7.pdf).

$g(r)$ is a probability distribution function, so the "oscillations" indicate that there are periodic values of particle separation which are more likely than average (peaks) or less likely than average (troughs). These "oscillations" are not necessarily sinusoidal (harmonic).

High density causes the "oscillations" by restricting the space which particles can move about in. When particles are forced close together they slide into relative positions where they can keep as much motion as possible - ie they form structures such as hexagonal close packing. At low densities particles have plenty of freedom to move about and can occupy all relative positions equally. All values of separation $r$ are equally likely - the probability distribution function is uniform.

3) Yes for long distances (= low densities) $g(r)$ becomes a uniform probability density function. That is right, this describes ideal gas particles in which there is no interaction - except at very short distances.

edited Nov 1 by sammy gerbil