Low-Energy Scattering: Scattering Lengths and Effective Ranges

Hello! Welcome to the final lesson in our module on quantum scattering theory.

In our last session, we generalized the S-matrix and T-matrix to handle multi-channel scattering, where particles can change their internal states or even transform into different particles. We also discovered a profound connection: the poles of the S-matrix in the complex energy plane correspond to the bound states of the system.

Today, we will apply this powerful machinery to a very important physical regime: low-energy scattering. We will see that in this limit, the intricate details of the scattering potential become less important, and the interaction can be described by just two universal parameters: the scattering length and the effective range. This analysis is fundamental to diverse fields, from nuclear physics to the study of ultra-cold atomic gases, a vibrant frontier of condensed matter physics.

Our learning outcome for this lesson is to analyze low-energy scattering and learn how to determine scattering lengths and effective ranges.

1. The Low-Energy Limit: s-Wave Dominance

First, what do we mean by "low energy"? In scattering, the relevant energy scale is set by the potential itself. A particle with momentum $p=\hslash k$ has an associated de Broglie wavelength $\lambda =2\pi /k$. If this wavelength is much larger than the range $R$ of the potential, the particle is "too big" to resolve the potential's detailed structure. This is the low-energy limit, mathematically expressed as:

Recall the partial wave expansion. The radial Schrödinger equation for the $l$-th partial wave includes an effective potential:

The second term is the repulsive centrifugal barrier. When $kR\ll 1$, an incoming particle has very little kinetic energy. For any partial wave with $l>0$, this energy is typically insufficient to overcome the centrifugal barrier. The particle is effectively repelled from the scattering center before it can even feel the potential $V(r)$.

Consequently, at low energies, only the $l=0$ partial wave, or s-wave, contributes significantly to the scattering. This has two immediate consequences:

1. The scattering is isotropic (uniform in all directions), since the $l=0$ spherical harmonic $Y_{00}$ is a constant.
2. The problem simplifies dramatically, as we only need to calculate one phase shift, ${\delta }_0(k)$.

2. The Scattering Length

To understand how the s-wave phase shift behaves as $k\to 0$, let's look at the radial equation outside the potential's range ($r>R$), where $V(r)=0$. For the s-wave ($l=0$), this is:

The general solution is $u_0(r)\propto \sin (kr+{\delta }_0)$. In the extreme low-energy limit ($k\to 0$), the $k^2$ term vanishes, and the equation simplifies even further to $u_0^{\prime \prime }(r)=0$. The solution is a straight line:

where $C$ and $a$ are constants. This linear form must match the asymptotic form $\sin (kr+{\delta }_0)$ in the limit of small $k$. By matching their logarithmic derivatives, one arrives at the fundamental definition of the scattering length, $a$:

This single parameter, $a$, which has units of length, characterizes the entire scattering process at zero energy.

To build a deeper intuition for the scattering length and its properties, please study the following resource.

As the reading explains, the scattering length has a simple and powerful geometric interpretation: $a$ is the position where the zero-energy radial wave function outside the potential (or its linear extrapolation) intersects the axis.

This leads to a few key cases for the sign of $a$:

• Repulsive Potential ($V(r)>0$): The wave function is pushed away from the origin, so it crosses the axis at a positive value. Thus, $a>0$.
• Weak Attractive Potential: The wave function is pulled towards the origin but not enough to form a bound state. The slope at $r=0$ is positive, but the function curves over and remains positive. We must extrapolate it backward to find the intercept, which is negative. Thus, $a<0$.
• Strong Attractive Potential: The potential is strong enough to "pull" the wave function over, causing it to have a node. The intercept with the axis is at a positive value. Thus, $a>0$.

Notice the crucial transition: as an attractive potential becomes just strong enough to support a bound state, the scattering length diverges, $|a|\to \infty$. This is precisely the connection to the S-matrix poles we discussed in the last lesson!

3. Effective Range and the Shape-Independent Approximation

The scattering length approximation, k cot(δ₀) = -1/a, is also known as the zero-range approximation. It is exact only at $k=0$. To improve this for small but non-zero energies, we can expand k cot(δ₀) in powers of $k^2$. The next term in this expansion defines the effective range, $r_0$:

This is the effective range expansion, or the shape-independent approximation. The name highlights its significance: at low energies, the scattering phase shift depends only on the two parameters $a$ and $r_0$, not on the specific shape of the potential (be it a square well, Gaussian, Yukawa, etc.). Two different potentials with the same $a$ and $r_0$ are indistinguishable in a low-energy scattering experiment.

The following reading provides the formal derivation and definition of $r_0$.

The key takeaway is that $r_0$ provides the first-order energy correction to the zero-range theory. While $a$ can be much larger than the potential range $R$, $r_0$ is typically on the order of $R$.

4. Connection to Bound States and Applications

We are now ready to make the connection between these low-energy scattering parameters and the existence of bound states. In the previous lesson, we saw that a bound state corresponds to a pole in the S-matrix on the positive imaginary k-axis, at $k=i\kappa$, where the bound state energy is $E_B=-\frac{{\hslash }^2{\kappa }^2}{2m}$.

The s-wave S-matrix element is $S_0(k)=e^{2i{\delta }_0}$. Using the relation cot(δ₀) = (e^(2iδ₀)+1) / (i(e^(2iδ₀)-1)), we can write $S_0$ in terms of k cot(δ₀):

Now, let's substitute our approximations:

1. Zero-Range Approximation: $k\cot {\delta }_0=-1/a$

   $$S_0(k)=\frac{-1/a-ik}{-1/a+ik}=\frac{k+i/a}{k-i/a}$$

   The S-matrix has a pole when the denominator is zero, i.e., at $k=i/a$. This is our bound state condition, $k=i\kappa$. Therefore, we find $\kappa =1/a$, which gives the bound state energy:

   $$E_B=-\frac{{\hslash }^2}{2m_r{a}^2}$$

   (Here we use reduced mass $m_r$ for the general two-body case). This remarkable formula connects a positive-energy scattering parameter ($a$) to a negative-energy bound state. It only works for a shallow bound state, where $|a|$ is large and positive.

2. Finite-Range Approximation: $k\cot {\delta }_0=-1/a+\frac{1}{2}{r}_0{k}^2$
   Substituting this into the S-matrix denominator and setting $k=i\kappa$ gives the condition for the pole:

   $$-\frac{1}{a}-\frac{1}{2}{r}_0{\kappa }^2+\kappa =0$$

   This is a quadratic equation for $\kappa$ that gives a more accurate estimate for the bound state energy, now including the effective range correction.

These approximations are incredibly powerful in practice. The following reading discusses this connection and shows two famous examples.

The universality is striking: the same formalism accurately describes the Helium dimer (an atomic physics system with energy $\sim 10^{-7}$ eV) and the deuteron (a nuclear system with energy $\sim 10^6$ eV), which differ in energy and length scales by many orders of magnitude.

For Your Interest: Numerical Computation

Given your background in writing simulations, you might find the final sections of the paper by Macêdo-Lima and Madeira particularly interesting. Sections 3 and 4 detail a complete numerical procedure for calculating the scattering length and effective range for any given potential, including the Lennard-Jones potential which is fundamental in condensed matter and molecular physics. This provides a direct path from a theoretical model potential to the experimentally measurable quantities $a$ and $r_0$.

Conclusion

This lesson concludes our module on scattering theory by examining the physically crucial low-energy limit.

Key Takeaways:

• Low-energy scattering ($kR\ll 1$) is dominated by the s-wave ($l=0$) and is therefore isotropic.
• The entire scattering process at near-zero energy is described by a single parameter, the scattering length ($a$). The low-energy cross section is simply ${\sigma }_0=4\pi a^2$.
• The effective range ($r_0$) provides the first energy-dependent correction, leading to the shape-independent approximation: $k\cot {\delta }_0\approx -1/a+\frac{1}{2}{r}_0{k}^2$.
• These two parameters, $a$ and $r_0$, universally describe low-energy interactions, regardless of the potential's specific shape.
• The scattering length is directly related to the energy of any shallow bound states the potential might support, beautifully linking the physics of positive-energy scattering and negative-energy bound states.

Preview of the Next Module:
We have now completed our comprehensive look at symmetries and scattering in quantum mechanics. In the next module, "Open Quantum Systems and Measurement," we will shift our focus. We'll move beyond the idealized case of isolated systems described by wavefunctions and introduce the density matrix. This powerful formalism is essential for describing mixed states, subsystems, entanglement, and the irreversible process of decoherence—topics that are central to modern quantum information science and condensed matter physics.