Pair Potential

How can we predict the properties of a solid, liquid, or gas, which contain countless interacting atoms? The answer lies in a powerful reductionist approach at the heart of physical science: the pair potential. This concept attempts to explain the behavior of a massive collective by understanding the simplest possible interaction—the force between just two particles. By modeling this fundamental dance of attraction and repulsion, we gain a key to unlock the secrets connecting the microscopic and macroscopic worlds.

This article addresses the challenge of bridging these scales, demonstrating how a simple function for pairwise energy can explain complex, observable phenomena. We will explore how this single idea provides a foundation for much of our understanding of matter. Across the following sections, you will discover the core theory behind pairwise interactions and see its remarkable reach across diverse scientific fields.

The journey begins in "Principles and Mechanisms," where we will dissect the celebrated Lennard-Jones potential and see how its features give rise to the basic properties of matter. We will then expand our view in "Applications and Interdisciplinary Connections" to witness how this concept is applied, adapted, and even transcended in materials science, biology, and beyond.

Imagine trying to understand the grand, chaotic, and beautiful dance of a bustling ballroom by only knowing one simple rule: how any two dancers interact. This is precisely the spirit behind one of the most powerful ideas in physical science: the ​​pair potential​​. We want to explain the behavior of trillions upon trillions of atoms—the properties of a gas, the flow of a liquid, the strength of a solid—by starting with the simplest possible question: how does the energy between just two of them change as we pull them apart?

This idea is a magnificent piece of scientific reductionism. If we can write down a function for the potential energy, $u(r)$, that depends only on the distance $r$ between two particles, we might just have a key to unlock the secrets of matter.

Let’s not be abstract. Let's look at a famous and wonderfully useful model for the interaction between two simple, neutral atoms, like two argon atoms. This is the celebrated ​​Lennard-Jones potential​​:

$u(r) = 4\epsilon \left[ \left(\frac{\sigma}{r}\right)^{12} - \left(\frac{\sigma}{r}\right)^{6} \right]$

At first glance, this might seem a bit intimidating, but it tells a very simple story—a story of a sensitive relationship. It has two parts, each describing a fundamental force.

First, consider the term proportional to $(\sigma/r)^{12}$. The distance $r$ is in the denominator, raised to a very high power. This means if you try to make $r$ very small—if you try to shove the two atoms on top of each other—this term becomes astronomically large and positive. A large positive potential energy is nature's way of saying "Don't do that!" It's an overwhelmingly strong repulsive force. This term models the fundamental principle that two pieces of matter cannot occupy the same space. It's a fantastically steep wall of energy. This is not just a mathematical convenience; it has profound physical consequences. In a liquid, the probability of finding two atoms closer than a certain distance (related to $\sigma$) is practically zero, precisely because the energy cost is too high. This is why the measured ​​radial distribution function​​, $g(r)$—a function that tells us the relative probability of finding a neighbor at distance $r$—plummets to zero for small $r$. This powerful repulsive term carves out a personal space for each atom.

Now, what about the other term, proportional to $-(\sigma/r)^6$? The negative sign means this is an attractive force. It's a gentler, longer-range attraction known as a London dispersion force, arising from fleeting, coordinated fluctuations in the electron clouds of the atoms. This is the "stickiness" that holds condensed matter together. Without this term, there would be no liquids or solids; every substance would just be a gas.

So, the pair potential tells a dynamic story: as two distant atoms approach each other, they feel a subtle pull. As they get closer, the attraction gets stronger. But if they get too close, an incredibly powerful repulsive force suddenly kicks in and pushes them apart.

What happens when these two forces—the short-range repulsion and the longer-range attraction—are perfectly balanced? That's where the potential energy is at its minimum. It's the most comfortable, stable distance for the two particles to be. We can find this "sweet spot" by doing a little calculus, taking the derivative of $u(r)$ and setting it to zero. For the Lennard-Jones potential, this equilibrium distance, let's call it $r_0$, turns out to be $r_0 = 2^{1/6}\sigma$.

And what is the energy at this most stable distance? If we plug $r_0$ back into the potential formula, we find that the minimum energy is exactly $-\epsilon$. This parameter, $\epsilon$, is called the ​​potential well depth​​. It is a measure of how strongly the two particles are bound together. A larger $\epsilon$ means a deeper well and a stronger "bond."

These two parameters, $\sigma$ (the effective size) and $\epsilon$ (the interaction strength), are the soul of the potential. They are microscopic quantities, but as we shall see, they orchestrate the macroscopic world. If we imagine a simple solid as particles sitting in their neighbors' potential wells, the total binding energy is just the sum of all these pairwise energies. For a simple chain of three atoms, with each interacting with its neighbor at the equilibrium distance, the total energy is just $2 \times (-\epsilon) = -2\epsilon$.

This is where the magic happens. How does the tiny, invisible dance of pairs manifest in the world we can see and touch?

Think about boiling water. You put a pot on the stove, add energy in the form of heat, and at 100°C, the liquid turns into steam. What is that energy, the ​​latent heat of vaporization​​, actually doing? It's giving each molecule enough energy to overcome the attractive forces of its neighbors—to climb out of the potential wells created by its fellow molecules. Using a simple model, we can directly relate this macroscopic quantity to our microscopic parameters. The energy required to vaporize one mole of a substance is directly proportional to the well depth, $\epsilon$, and the number of neighbours each molecule has in the liquid. The stickier the molecules (larger $\epsilon$), the more energy it takes to pull them apart into a gas.

Here's another beautiful example: the cooling of a gas. An ideal gas, with no intermolecular forces, wouldn't change its temperature if you let it expand into a vacuum. But a real gas, like argon, does. As the gas expands, the average distance between the argon atoms increases. They have to move "up" the potential energy curve, away from the attractive well. Where does the energy to do this come from? It can't come from outside. It must come from the gas itself. The atoms convert some of their kinetic energy—their motion—into potential energy. Since temperature is a measure of the average kinetic energy, the gas cools down! This phenomenon, the ​​Joule-Thomson effect​​, is used in industrial processes to liquefy gases and is a direct, measurable consequence of the attractive part of the pair potential.

Even the famous ​​Van der Waals equation​​, one of the first and most successful attempts to describe real gases, has its roots in the pair potential. The $a$ parameter in that equation, which corrects the pressure for intermolecular attractions, is no mere fudge factor. It can be derived directly from the integral of the attractive part of a model pair potential, quantifying the total "stickiness" of the gas.

This is all very nice, you might say, but how do we know what the potential $u(r)$ looks like? We can't take two atoms with a pair of tweezers and measure the force between them. The answer is wonderfully indirect: we watch the whole ballroom and infer the rules of the dance.

By scattering X-rays or neutrons off a liquid, physicists can measure the radial distribution function, $g(r)$, which, as we've seen, maps out the average structure of the liquid. Now, there is a profound connection in statistical mechanics between this structure and the underlying potential. In the simple case of a very low-density gas, the relationship is beautifully direct:

$g(r) \approx \exp\left(-\frac{u(r)}{k_B T}\right)$

where $k_B$ is Boltzmann's constant and $T$ is the temperature. This means we can simply rearrange the equation to find the potential: $u(r) \approx -k_B T \ln(g(r))$. By measuring the structure, we can directly map out the potential energy landscape!

For a dense liquid, the situation is more complex, but the principle remains. The total potential energy of a liquid (its ​​excess internal energy​​ compared to an ideal gas) can be calculated if we know both the pair potential $u(r)$ and the structure $g(r)$, through an integral that sums up the energy of all pairs, weighted by the probability of finding them at each distance. The structure and the interactions are two sides of the same coin.

The idea of the pair potential is powerful, but nature, as always, has a few more tricks up her sleeve. The assumption that the total energy is just the sum of all pairs is an approximation—a very good one in many cases, but an approximation nonetheless.

In a dense liquid, the force between particle 1 and particle 2 isn't taking place in a vacuum. It is surrounded by a sea of other particles, which jostle, screen, and mediate the interaction. The effective potential energy we would infer from the structure, $g(r)$, is not the "true" two-body potential $u(r)$, but something we call the ​​potential of mean force​​, $w(r)$. It includes not only the direct interaction between 1 and 2 but also all the indirect effects of the surrounding crowd of particles. The two become equal only in the limit of a very dilute gas, where there is no crowd to worry about.

The story gets even deeper. Sometimes, the energy of three particles together is not just the sum of the energies of the three pairs (1-2, 1-3, and 2-3). There can be a genuine ​​three-body force​​ that only exists when all three are present. This is like a conversation among three people that has a different character than three separate two-person conversations. These many-body forces are crucial for accurately describing the properties of dense gases and liquids.

Nowhere is the limitation of the pair potential more stark than in metals. In a metal, the outer electrons are not tied to their parent atoms. They are delocalized, forming a "sea" of electrons in which the positive ions are embedded. The energy of a metal ion doesn't depend on its pairwise interactions with its neighbors. It depends on the local density of the electron sea it finds itself in. This is a profoundly many-body concept. A simple pairwise potential like Lennard-Jones completely fails to describe the elastic properties of metals; for example, it predicts a relationship between elastic constants (the ​​Cauchy relation​​, $C_{12} = C_{44}$) that is simply not observed in most real metals. To model metals correctly, scientists had to invent new approaches, like the ​​Embedded-Atom Method (EAM)​​, which explicitly includes this density-dependent, many-body energy.

And so, our journey from a simple dance of two to the complex ballroom of many reveals a beautiful arc in science. We start with a simple, powerful idea—the pair potential—that explains a vast range of phenomena with stunning elegance. But by pushing its limits, we discover where it falls short, and in doing so, we are forced to uncover a deeper, richer, and more accurate picture of how the world truly works. The pair potential is not the final answer, but it is an essential and beautiful chapter in the story of matter.

In the previous chapter, we dissected the idea of the pair potential, the simple, elegant rule describing how two particles interact. We saw how this concept forms the very foundation of our understanding of matter. But physics is not just about abstract principles; it is about explaining the world we see, touch, and are a part of. Now that we have learned the notes and scales of pairwise interactions, let's step back and listen to the music. In this chapter, we will embark on a journey to see how this one idea—the potential between a pair of particles—blossoms into a rich tapestry of applications, weaving together materials science, chemistry, biology, and even the very fabric of spacetime.

Imagine a crystal, a vast, orderly city of atoms. How do we calculate its total energy? The task seems daunting, like trying to map the social network of the entire city at once. The principle of pairwise interaction gives us a wonderfully simple starting point: the total potential energy is just the grand sum of all the one-on-one "conversations" between every possible pair of atoms in the crystal.

This surprisingly simple accounting has profound consequences. It allows us to calculate real, measurable properties of materials from first principles. Consider the phenomenon of surface tension, the force that makes water bead up and allows insects to walk on its surface. What is a surface? It is a boundary where the perfect, repeating structure of a crystal is broken. An atom at the surface has fewer neighbors than an atom deep inside. Creating a surface means "breaking" the bonds of the atoms that are now exposed. The energy required to break these bonds, per unit area, is precisely the surface tension. Using a model like the Lennard-Jones potential, we can sum the energy of these lost interactions and derive a direct expression for surface tension, $\gamma_{sv}$, in terms of the microscopic potential parameters $\epsilon$ and $\sigma$. Suddenly, a macroscopic property is tied directly to the depth and range of the potential between two atoms.

But is it always so simple? What about metals? In a metal, the outer electrons are not tethered to their parent atoms; they form a delocalized "sea" of charge in which the positive atomic cores are immersed. Here, the interaction between two atomic cores is not a private affair. The presence of other cores and the surrounding electron sea changes their interaction. A simple pairwise sum is no longer enough. Materials scientists have developed more sophisticated models like the Embedded Atom Model (EAM) to capture this. In the EAM, the energy of an atom has two components: the sum of its direct pairwise "chats" with its neighbors, and a second, crucial term called the embedding energy. This embedding energy is like a "community fee" the atom must pay to exist within the local electron sea. It depends on the total electron density contributed by all its neighbors, making it an inherently many-body effect. By extending the simple pairwise concept, the EAM allows for highly accurate calculations of properties like the cohesive energy of metals, demonstrating how our models evolve to embrace the richer physics of real materials.

The world is full of surfaces and interfaces. How does a single atom "feel" the presence of a vast, solid surface? Again, we can find the answer by summing up pair potentials. Imagine integrating the tiny pull of a Lennard-Jones potential over every atom in a semi-infinite solid. The result is a beautiful and powerful expression for the interaction between an atom and a surface. A fascinating mathematical consequence emerges: the individual $r^{-6}$ and $r^{-12}$ dependencies of the pair potential, when integrated over the three-dimensional half-space of the solid, transform into dependencies of $z^{-3}$ and $z^{-9}$ on the perpendicular distance $z$ from the surface. This is a direct illustration of how geometry shapes the laws of force. This principle is not just a curiosity; it's the foundation of adsorption, catalysis, and the operation of tools like the Atomic Force Microscope, which "feels" surfaces atom by atom.

This idea of summing forces extends to other geometries. If you place a particle inside a hollow sphere, what is the net force from the sphere's atoms? Isaac Newton famously proved for gravity that the net force from a uniform spherical shell on a mass inside is zero. For van der Waals forces, the situation is different. Integrating a pairwise potential like $w(r) = -C/r^6$ over the volume of a hollow shell reveals a net potential that depends on the shell's inner and outer radii. This has tangible implications for designing modern materials, such as drug-delivery systems where a medication-laden nanoparticle is encapsulated within a protective shell.

Perhaps the most captivating example of this collective force comes from the biological world. The ability of a gecko to scamper up a glass wall seems to defy gravity. The secret lies not in suction or glue, but in the collective strength of billions of tiny van der Waals forces. Each of the microscopic, hair-like spatulae on a gecko's foot pad interacts weakly with the surface. Individually, the force is minuscule, but summed together, they are strong enough to support the gecko's entire body weight. We can model this spectacular feat of nature by representing the gecko's foot and the wall as two grids of particles interacting via the Lennard-Jones potential and simply summing all the pairwise forces to find the total adhesion. It's a powerful reminder that the macroscopic world of our experience is built upon the silent, ceaseless summation of microscopic forces.

So far, we have lived in a world where the total interaction is the sum of its parts. But sometimes, three's a crowd. The interaction between particles A and B can be altered by the presence of a third particle, C. This is the realm of many-body forces, where the simple pairwise addition rule breaks down.

A dramatic example comes from the subatomic world, in the theory of Quantum Chromodynamics (QCD). A simple model of a baryon (like a proton) consists of three quarks. The potential energy holding them together is not just the sum of three pairwise interactions. It includes a dominant "confinement" term, envisioned as a network of "flux tubes" or strings of energy connecting the quarks. The energy of this confining potential is proportional to the minimum total length of the strings needed to connect the three quarks. This connection length is not a sum of pairwise distances but a collective property of the three-quark system. This is a fundamental aspect of the strong nuclear force, where pairwise thinking is simply incomplete.

We have already encountered a more subtle example in the Embedded Atom Model for metals. The embedding energy, $F(\rho)$, depends on the electron density $\rho$ at an atom's location. Since $\rho$ is the sum of contributions from all surrounding atoms, the energy of any single atom depends on the collective arrangement of its entire neighborhood. This is a many-body effect, elegantly disguised within a framework that still retains the spirit of pairwise thinking. These examples teach us an important lesson: while the pair potential is an incredibly powerful tool, we must always be prepared for nature to be more communal, where the whole is not merely the sum of its parts.

This brings us to a deep and essential question: Where do pair potentials like Lennard-Jones come from in the first place? Are they fundamental laws of nature? The answer is no. For complex systems, they are effective potentials—brilliant approximations born from a process called "coarse-graining."

Imagine trying to simulate a long, flexible polymer chain in a solvent. Tracking the motion of every single atom in the polymer and the surrounding water molecules would be computationally impossible. The strategy is to "zoom out": we can represent a group of, say, 10 monomers as a single "bead." Now we have a simpler system of beads. But what is the potential between two such beads? It must be an effective potential that implicitly averages over all the intricate, wiggling motions of the underlying atoms.

The science of finding these effective potentials is a cornerstone of modern computational chemistry and physics. A key piece of the puzzle is the radial distribution function, $g(r)$, which tells us the probability of finding a particle at a distance $r$ from another. Henderson's theorem, a profound result from liquid-state theory, provides a theoretical anchor: for a system governed by pair forces, the pair potential $u(r)$ is uniquely determined by the structure $g(r)$ at a given temperature and density. This means that if we can create a model potential that reproduces the correct structure of a fluid, we have found the unique effective pair potential for that system.

In practice, this is a delicate art. An effective potential that correctly reproduces the structure ($g(r)$) might not correctly reproduce thermodynamic properties like the pressure or the free energy, because many-body effects contribute differently to structure and energy. A clever solution is to add a coordinate-independent, but density-dependent, energy term to the model. This term acts like a uniform background energy that adjusts the overall thermodynamics without altering the forces between particles, thereby leaving the structure intact.

So how do we practically find the potential? One powerful method is Iterative Boltzmann Inversion (IBI). The process is a beautiful dialogue between simulation and a target (often from experiment). You start with a guess for the potential—a good first guess is the "potential of mean force," derived directly from the target structure via $u(r) \approx -k_B T \ln g_{\text{target}}(r)$. Then you run a simulation with this potential and calculate the resulting structure, $g_{\text{sim}}(r)$. You compare it to your target. Is your simulated fluid too clumpy at a certain distance? That means your potential is too attractive there; you must adjust it to be more repulsive. Is it too spread out? Your potential is too repulsive; make it more attractive. You apply a correction, $u_{\text{new}}(r) = u_{\text{old}}(r) + k_B T \ln [g_{\text{sim}}(r)/g_{\text{target}}(r)]$, and repeat the process. Iteration by iteration, the simulated structure is driven towards the real one, and the effective potential emerges from this process. This is how we can take experimental X-ray or neutron scattering data, which gives us $g(r)$, and reverse-engineer the effective forces that create the structure we observe.

We end our journey with a final, breathtaking twist that showcases the profound unity of physics. Let's return to something as seemingly mundane as a real gas, described by the van der Waals equation of state. The parameter $a$ in that equation accounts for the attractive forces between molecules. These forces mean that the molecules have a negative potential energy when they are close to each other.

Now, let us invoke the most famous equation in all of science: Albert Einstein's $E = mc^2$. This equation does not just apply to nuclear bombs and particle accelerators. It states that all forms of energy have a mass equivalent. This includes kinetic energy, thermal energy, and, crucially for us, potential energy.

The total energy of a van der Waals gas includes the kinetic energy of its molecules and the negative potential energy from their mutual attractions. Therefore, the total relativistic mass of a container of this gas is the sum of the masses of its molecules plus the mass equivalent of their kinetic energy minus the mass equivalent of their binding energy. In other words, a bottle of real gas weighs slightly less than it would if its molecules felt no attraction to one another. The amount of this "mass defect" is precisely $\Delta M = -a n^2 / (V c^2)$. This is the exact same principle that explains why the mass of a helium nucleus is less than the mass of two separate protons and two separate neutrons. Binding energy reduces mass. To see this principle at work not just in the heart of a star but in a simple container of gas is a stunning testament to the universality of physical law.

From the strength of a crystal to a gecko's grip, from the heart of a proton to the very mass of a gas, the simple idea of a pair potential has proven to be one of the most fertile and far-reaching concepts in science. It shows how the intricate dance of the macroscopic world can often be understood by patiently listening to the quiet, pairwise conversations between its smallest constituents.