Woods-Saxon Potential

How do we model the chaotic swarm of protons and neutrons packed into an atomic nucleus? The forces at play are immensely powerful yet incredibly short-ranged, making a particle-by-particle description nearly impossible. Nuclear physics overcomes this complexity through the ​​mean-field approximation​​, which simplifies the countless individual interactions into a single, average potential landscape. While simple models like the square well offer a starting point, they fail to capture the nucleus's true nature: a dense core with a fuzzy, diffuse surface. This is the knowledge gap that the Woods-Saxon potential masterfully fills.

This article delves into this cornerstone of the nuclear shell model. You will learn about the elegant mathematical form of the Woods-Saxon potential and how its parameters directly relate to tangible nuclear properties. The first chapter, ​​Principles and Mechanisms​​, will break down the formula, connect it to physical concepts like nuclear density and the Schrödinger equation, and reveal its crucial partnership with the spin-orbit interaction. Following this, the ​​Applications and Interdisciplinary Connections​​ chapter will demonstrate its predictive power, showing how this single model illuminates a vast range of phenomena, from fusion reactions and alpha decay to the shifting structures of exotic nuclei at the edge of stability.

To understand the heart of an atomic nucleus, we must first understand the world a single proton or neutron—a ​​nucleon​​—inhabits. Imagine being a lone nucleon inside this bustling, dense city of particles. What forces do you feel? What rules govern your motion? The nuclear world is governed by the strong force, an interaction so powerful it easily overcomes the electric repulsion between protons, but also so short-ranged that it vanishes almost completely just outside the nucleus. Our task, as physicists, is to paint a portrait of the average force field, or ​​potential​​, that a nucleon experiences. This is the essence of the ​​mean-field approximation​​: we smooth out the chaotic, individual tugs from every other nucleon into a single, static landscape.

What would be the simplest landscape to imagine? Perhaps a hole in the ground with a flat bottom and perfectly vertical walls—a ​​spherical square well​​. Inside the hole, you feel a constant pull downwards; outside, you feel nothing. This is a fine first guess, a physicist's caricature, but nature is rarely so angular. Experiments that probe the size and shape of nuclei tell a more nuanced story. They reveal that a nucleus isn't a tiny billiard ball with a hard edge; it's more like a droplet of liquid. It has a region of remarkably constant density deep in its interior, but its edge is fuzzy, thinning out over a small but finite distance. This is the concept of a ​​diffuse surface​​.

Our potential landscape must reflect this physical reality. It should be deep and relatively flat in the center, corresponding to the constant pull within the dense nuclear interior. Then, as we approach the edge, it shouldn't stop abruptly but should smoothly slope upwards, becoming zero outside the nucleus where the strong force no longer reaches. This is the challenge: to find a simple, elegant mathematical function that captures this sophisticated shape.

In the 1950s, Roger Woods and David Saxon proposed just such a function, one that has become a cornerstone of nuclear physics. The ​​Woods-Saxon potential​​ is a masterpiece of physical intuition captured in a simple formula:

At first glance, it might look intimidating, but let's break it down piece by piece, for each part tells a story.

• ​​The Depth, $-V_0$​​: This is the potential at the very center of the nucleus. The parameter $V_0$ is a positive number (typically around $50$ MeV, or Mega-electron-Volts), representing the strength of the attraction. The negative sign is the physicist's convention for an attractive force—it's a potential well, not a hill. It's the measure of how strongly the nucleus as a whole pulls on a single nucleon.

• ​​The Radius, $R$​​: This parameter defines the size of the nucleus. It's not the "edge" in a hard sense, but rather a characteristic radius. A wonderful feature of this formula is what happens exactly at $r=R$. The exponential term becomes $\exp(0)=1$, and the potential is precisely $V(R) = -V_0 / (1+1) = -V_0/2$. So, $R$ is the radius at which the potential has faded to exactly half its central depth. This provides a natural and consistent way to define the nuclear surface.

• ​​The Diffuseness, $a$​​: This is the "fuzziness" parameter. It has units of distance and controls how thick the surface is. If $a$ were vanishingly small, the denominator would abruptly jump from $1$ to infinity as $r$ crosses $R$, and we would recover our old, unrealistic square well. But for a finite, physical value of $a$ (typically around $0.65$ fm, or femtometers), the potential makes a smooth, graceful transition from the deep interior to the zero-potential exterior. It's the parameter that gives the potential its realistic, rounded shoulder.

Let’s make this tangible. For a heavy nucleus like Lead-208, the parameters are roughly $V_0 = 50$ MeV, $R \approx 7.4$ fm, and $a = 0.65$ fm. If we calculate the potential a nucleon feels:

• At the very center ($r=0$), the exponential term is minuscule, and the potential is almost exactly $-50$ MeV.
• At the half-potential radius ($r=R=7.4$ fm), the potential is exactly $-25$ MeV.
• Just a short distance further out, at $r=R+2a \approx 8.7$ fm, the potential has already faded to about $-6$ MeV. You can see how the force field is strong inside but dies away rapidly and smoothly at the surface.

There's another beautiful piece of physics hidden in the radius $R$. It is not a fixed number but scales with the total number of nucleons, $A$, according to the rule $R = r_0 A^{1/3}$, where $r_0$ is a constant. This is the same scaling you'd find for the radius of a water drop as you add more molecules! It implies that the volume of a nucleus is directly proportional to the number of particles in it. In other words, ​​nuclear matter has a constant density​​, a property known as saturation. This simple scaling law connects the abstract potential to the tangible, liquid-drop nature of the nucleus.

A potential is more than a static landscape; it dictates motion by creating forces. The force is related to the slope of the potential. Deep inside the nucleus, where the potential is nearly flat like a vast plateau, the force is very weak. Far outside, where the potential is also flat (at zero), the force is zero. So, where is the pull on the nucleon strongest? It is strongest where the slope is steepest. For the Woods-Saxon potential, a little bit of calculus shows that the force is maximal exactly at $r=R$, the half-potential radius. This confirms our intuition: the "action" is at the surface.

Now for the most important step: we must remember that a nucleon is not a classical particle. It is a quantum object, a wave governed by the ​​Schrödinger equation​​. To find the allowed energy states for a nucleon, we must solve this equation using the Woods-Saxon potential. For a particle with orbital angular momentum $l$, the radial part of the Schrödinger equation takes on a beautifully simple, one-dimensional form:

Here, $u_l(r)$ is the (reduced) radial wave function, $E$ is the energy, and $V_{WS}(r)$ is our Woods-Saxon potential. The equation represents a delicate balance. The first term is related to the kinetic energy, the "wiggliness" of the wave. The second term is the potential energy, consisting of two parts: our attractive Woods-Saxon well pulling the nucleon inward, and a repulsive ​​centrifugal barrier​​ that pushes the nucleon outward, the quantum "price" for carrying angular momentum.

For a nucleon to be truly ​​bound​​ to the nucleus, its wave function must be contained. This imposes strict boundary conditions: the wave function $u_l(r)$ must be zero at the center ($r=0$) and must decay to zero at infinity. It is trapped. Only certain discrete, quantized energies $E$ will allow for solutions that satisfy these conditions.

This realism is a key advantage of the Woods-Saxon potential over simpler models like the harmonic oscillator. A harmonic oscillator potential, $V(r) \propto r^2$, grows infinitely large with distance. It is a "cosmic prison" from which no particle, no matter its energy, can ever escape. It can only host bound states, and its wave functions decay in a Gaussian manner ($\sim \exp(-\alpha r^2)$), which is too fast to accurately describe the tenuous wave function tails of real, weakly bound nucleons. The Woods-Saxon potential, by correctly going to zero at large distances, not only yields the correct exponential decay for its bound states ($\sim \exp(-\kappa r)$) but can also describe unbound ​​continuum states​​ ($E>0$). These are essential for modeling nuclear reactions, where a particle can come in from the outside, interact with the nucleus, and fly away again.

One might still ask: is this elegant formula just a clever fit to the data, or is there a deeper reason for its success? The answer is a resounding "yes". The mean-field potential that a nucleon feels is, after all, the collective result of its interactions with all the other nucleons. It has been shown that if you start with the experimentally observed fact that the density of nucleons follows a shape very similar to the Woods-Saxon form (called a Fermi distribution), and you then calculate the average potential generated by this density, the result is, to a very good approximation, the Woods-Saxon potential itself!. The shape of the potential is a direct and beautiful reflection of the distribution of the matter that creates it.

This picture of a simple central potential is remarkably successful, but it is incomplete. It correctly predicts the first few "magic numbers" of nuclear physics ($2, 8, 20$)—the nucleon counts at which nuclei are exceptionally stable—but it fails for all the heavier ones. The model is missing a crucial ingredient.

This ingredient is the ​​spin-orbit interaction​​. Each nucleon has an intrinsic spin, like a tiny spinning top. A profound relativistic effect, related to an effect called ​​Thomas precession​​, causes the nucleon to feel an additional force that depends on the orientation of its spin relative to its orbital motion. The story of this interaction contains a wonderful twist. A naive application of the theory gives a force with the wrong sign—one that would make the nuclear shells even less realistic! The correct, and much stronger, interaction only emerges from a more complete relativistic treatment. The outcome is that states where the nucleon's spin and orbit are aligned ($j = l + 1/2$) have their energy significantly lowered, while states where they are anti-aligned ($j=l-1/2$) are raised. This effect is most dramatic for states with high orbital angular momentum, $l$.

This is not a minor tweak; it is a seismic shift in the nuclear energy landscape. Consider the $1f_{7/2}$ state, which has $l=3$. The flat-bottomed shape of the Woods-Saxon potential already lowers high-$l$ states relative to a harmonic oscillator. Then, the powerful spin-orbit force grabs this state and, because its spin and orbit are aligned ($j = 3 + 1/2 = 7/2$), yanks its energy dramatically downwards. It is pulled down so far that it leaves its original shell and becomes an "intruder" in the shell below, creating the new magic number at $28$. This mechanism, the interplay between the realistic shape of the Woods-Saxon potential and the strong spin-orbit force, is the key to understanding the full sequence of magic numbers and the structure of almost all atomic nuclei. The journey that began with a simple, intuitive picture of a fuzzy-edged potential well has led us to the very heart of the modern nuclear shell model.

Now that we have acquainted ourselves with the form and function of the Woods-Saxon potential, we can begin to appreciate its true power. It is far more than a tidy mathematical formula; it is a key that unlocks a vast array of phenomena in the nuclear world. Like a well-crafted map, it may not capture every pebble and blade of grass of the underlying territory, but it reliably guides us through the mountains, valleys, and rivers of nuclear structure and dynamics. Its beauty lies not in being an exact replica of reality, but in being an exceptionally effective and insightful approximation.

The remarkable "transferability" of this potential is where our journey begins. Why should one simple formula work for the helium nucleus, with its four nucleons, all the way up to lead, with over two hundred? The secret lies in one of the most profound properties of nuclear matter: saturation. Nuclei are not like drops of water that can be compressed further and further; they maintain a nearly constant interior density. This means that a nucleus's volume is directly proportional to the number of particles, $A$, inside it. A simple geometric argument then tells us that its radius, $R$, must scale with the cube root of the mass number, $A^{1/3}$. Furthermore, the "fuzziness" of the nuclear surface—its diffuseness, $a$—isn't determined by the overall size of the nucleus, but by the short range of the nuclear force itself and the quantum jitters of nucleons at the edge. Thus, the radius $R$ grows predictably with $A$, while the diffuseness $a$ remains nearly constant for most nuclei. This elegant scaling is the foundation of global models that allow us to use a single, unified framework to describe reactions across the entire chart of nuclides.

Let us first watch the grand drama of two nuclei interacting. Imagine firing a projectile nucleus, like oxygen, at a heavy target, like lead. Two great forces are at play. From afar, the two positively charged nuclei feel a powerful electrostatic repulsion—the Coulomb force—that tries to push them apart. This force grows stronger and stronger as they approach. But if they get close enough, the intensely attractive, short-range strong nuclear force, described beautifully by our Woods-Saxon potential, springs into action and tries to pull them together.

The result is a dramatic competition: a potential energy landscape with a huge hill—the Coulomb barrier—that the projectile must climb. The peak of this hill, the fusion barrier, is found precisely where the outward push of the Coulomb force is momentarily balanced by the inward pull of the nuclear force's attraction, a pull whose strength is governed by the gradient of the Woods-Saxon potential. Only by surmounting this barrier can the two nuclei merge in a fusion reaction, releasing enormous amounts of energy. The height of this barrier, a critical parameter for understanding stellar nucleosynthesis and terrestrial fusion energy, is born directly from the interplay of the well-understood Coulomb potential and the phenomenological genius of the Woods-Saxon form.

A similar, though simpler, drama unfolds when a single neutron approaches a nucleus. Having no charge, the neutron feels no Coulomb repulsion. It sees only the inviting well of the Woods-Saxon potential. Does it get captured? Here, classical mechanics gives us a wonderful intuition. The outcome depends on the neutron's energy and its "impact parameter"—how far off-center its trajectory is aimed. The combination of the attractive nuclear potential and the "centrifugal barrier" (a fictitious repulsive force that accounts for angular momentum) creates an effective potential landscape. If the incoming neutron has just the right energy and impact parameter to skim the top of this effective barrier, it will be captured. By analyzing the shape of this barrier, which is dominated by the Woods-Saxon potential, we can calculate the maximum impact parameter for capture, which in turn tells us the capture cross-section—the effective "size" of the target nucleus for that reaction.

The nucleus is a quantum realm, and the Woods-Saxon potential truly shines when we consider its quantum mechanical consequences. Consider alpha decay, the process by which a heavy nucleus like uranium spits out a helium nucleus (an alpha particle). The puzzle that baffled early physicists was how the alpha particle, with its relatively low energy, could escape the immensely deep attractive potential well holding the nucleus together.

George Gamow solved the mystery with quantum tunneling: the alpha particle doesn't climb over the barrier; it "tunnels" through it. His initial model used a simplified square-well potential, but we know the edge of a nucleus is not a sharp cliff. It is a gentle slope, a "diffuse" surface described by the Woods-Saxon potential. Replacing the crude square well with the more realistic Woods-Saxon form provides a crucial correction to the tunneling probability. The "fuzzier" the edge of the potential (the larger the diffuseness parameter $a$), the thinner the barrier is at its base, making it slightly easier for the alpha particle to leak out. This seemingly small change in the potential's shape has an exponential effect on the decay rate, altering the calculated half-life of the nucleus by orders of magnitude. The realistic shape of the potential is not just an aesthetic improvement; it is essential for quantitative prediction.

We have a tendency to picture nuclei as perfect little spheres, but nature is more creative than that. Many nuclei, particularly those far from the "magic number" closed shells, are deformed, taking on shapes more like an American football or a doorknob. The Woods-Saxon model accommodates this reality with remarkable elegance. Instead of a fixed radius $R$, we can let the radius depend on the angle relative to the nucleus's symmetry axis.

This "deformed Woods-Saxon potential" is the heart of models that describe the collective behavior of nuclei, such as rotation and vibration. When a projectile, like a proton, scatters off a deformed nucleus, it can transfer energy and set the nucleus spinning, exciting it to a higher rotational energy level. The interaction responsible for this kick is not the spherical part of the potential, but the small, angle-dependent deformation. By expanding the potential in terms of the deformation parameter $\beta_2$, we can isolate a "transition potential" whose strength depends directly on the derivative of the original Woods-Saxon potential. It is the physics of the nuclear surface, the region of greatest change, that governs these collective excitations.

Perhaps the most profound application of the Woods-Saxon potential is its role in the nuclear shell model, our most successful picture of nuclear structure. To explain the observed nuclear "magic numbers"—the nucleon counts at which nuclei are exceptionally stable—one must include a strong spin-orbit interaction. This is a purely quantum mechanical force that depends on the coupling between a nucleon's orbital motion ($L$) and its intrinsic spin ($S$). Its strength is not constant but depends on location, being proportional to the gradient of the central potential: $\frac{1}{r}\frac{dV}{dr}$.

This means the spin-orbit force is a surface phenomenon; it is strongest where the potential is steepest. Now, let us venture to the frontiers of nuclear existence, to the "drip lines" where nuclei are so saturated with neutrons that they can barely hold onto the last one. These exotic nuclei develop a "neutron skin," or even a "halo," meaning their surface is far more diffuse and spread out than that of a stable nucleus. In the language of our model, they have a larger diffuseness parameter $a$.

The consequence is breathtaking. A larger $a$ means a more gradual slope for the potential. A more gradual slope means a smaller gradient $\frac{dV}{dr}$. And a smaller gradient means a weaker, or "quenched," spin-orbit force. This quenching can be so significant that the traditional magic numbers that dictate stability in familiar nuclei actually vanish, and new magic numbers appear. The simple "fuzziness" parameter in our potential is directly linked to the changing rules of nuclear structure at the limits of stability.

Throughout this exploration, we have treated the Woods-Saxon potential as a given. We fit its parameters—$V_0$, $R$, $a$—to experiment and use it to predict other phenomena. But what is it, really? Is it a fundamental law? The answer is no, and the truth is even more beautiful.

In a full, microscopic many-body theory, a nucleon moving through a nucleus interacts with all other nucleons in a dizzyingly complex dance. The optical model potential is the result of averaging over all these complex interactions to produce a simple, effective one-body potential that describes the elastic scattering channel. This true microscopic potential is a formidable object: it is nonlocal (the force at point $\mathbf{r}$ depends on the wavefunction at point $\mathbf{r}'$), it is complex (its imaginary part accounts for reactions that remove the nucleon from the elastic channel), and it is fundamentally energy-dependent.

The local, real Woods-Saxon potential is a brilliant phenomenological approximation to this much more complicated reality. Its success tells us that, for many purposes, the nonlocal effects can be mimicked by a local potential with an energy-dependent depth. It is not the final theory, but an indispensable tool. Physicists must still navigate practical challenges, such as "ambiguities" where different sets of Woods-Saxon parameters can describe the same low-energy scattering data, requiring careful analysis to disentangle them.

The story of the Woods-Saxon potential is a perfect example of the physicist's art. It is an admission that we cannot solve the full nuclear many-body problem from first principles in every case. But it is also a declaration that by capturing the essential physics—a strong attraction, a finite size, and a diffuse surface—in a simple, elegant form, we can build a model of astonishing predictive power. It is a testament to how a thoughtful approximation can illuminate the profound and beautiful unity of the physical world.