Relaxation of a single-particle excitation in a Fermi system within the diffusion approximation of kinetic theory

This paper numerically investigates the time evolution of single-particle excitations in a Fermi system using a nonlinear diffusion equation, proposing a method to distinguish between the relaxation of the excitation and the nuclear core while revealing discrepancies between the resulting relaxation times and previously estimated kinetic coefficients.

The Gist

Here is an explanation of the paper, translated into everyday language using analogies to make the complex physics concepts accessible.

The Big Picture: A Crowd of Particles Trying to Calm Down

Imagine a crowded dance floor inside a nuclear ballroom. The dancers are nucleons (protons and neutrons), and they are part of a Fermi system (like an atomic nucleus).

Usually, these dancers are moving in a very organized, chaotic-but-steady way. They have a "dance floor limit" (the Fermi surface) where they can't really go any faster without bumping into someone.

Now, imagine someone throws a single, energetic dancer onto the floor who is moving much faster than everyone else. This is a single-particle excitation. The paper asks: How long does it take for this fast dancer to slow down and blend in with the rest of the crowd?

The author, Sergiy Lukyanov, uses a mathematical tool called the Diffusion Approximation to answer this. Think of this tool as a way to predict how a drop of ink spreads out in a glass of water, but instead of ink, it's energy and momentum spreading through the nuclear dance floor.

The Main Problem: The "Too Fast" Clock

The author ran a computer simulation to see how long this relaxation (calming down) takes. He found a strange problem:

• What we expect: Based on previous studies, the "relaxation time" (the time it takes to calm down) should be about 10 to 100 attoseconds (a very, very short time, but still measurable).
• What the math said: The simulation showed the dancers calming down in about 1 attosecond.

The simulation was 10 times too fast. The math predicted the system settled down much quicker than reality suggests. The author calls this a "discrepancy."

The New Idea: Separating the "Core" from the "Guest"

In the past, scientists looked at the whole dance floor at once. They couldn't tell if the fast dancer was slowing down because of the music, or if the whole crowd was just getting tired.

Lukyanov's innovation was to separate the two processes:

1. The Core (The Crowd): The background dancers who are just doing their usual routine.
2. The Excitation (The Guest): The single fast dancer who was thrown in.

He created a method to watch the "Guest" slow down independently from the "Crowd."

The Analogy:
Imagine a noisy party.

• Old way: You measure how long it takes for the whole room to get quiet.
• New way: You measure how long it takes for one specific loud person to stop shouting, while ignoring the background chatter.

What He Discovered

By separating the two, he found some interesting patterns:

1. The Guest is faster than the Crowd: The single fast dancer (the excitation) actually settles down faster than the entire room does. This makes sense; it's easier for one person to stop running than to get a whole stadium to stop moving.
2. Energy matters:
   • If the guest is super energetic (very far from the dance floor limit), it takes them longer to calm down.
   • However, for the whole room, adding more energy actually makes the whole system settle down faster.
3. Size matters: In a bigger nucleus (a bigger dance hall), the single dancer calms down slower because there are more people to interact with, but the "Guest" effect is actually weaker in larger groups.

The Big Mystery: Why is the Math Wrong?

The author tried to fix the "too fast" result by tweaking the Kinetic Coefficients.

• Think of these coefficients as the "friction" of the dance floor.
  • High friction = dancers slide less, they stop quickly.
  • Low friction = dancers slide a lot, they take longer to stop.

He realized that to get the math to match reality (to make the relaxation time longer), he would have to assume the "friction" of the nuclear dance floor is 20 to 200 times lower than what other scientists have calculated before.

The Conclusion:
The author concludes that there is a fundamental disagreement in how we calculate this "friction" (the diffusion and drift coefficients).

• Either our current understanding of how particles bump into each other is wrong.
• Or, the "Diffusion Approximation" (the tool used to solve the problem) is missing something important, perhaps some quantum effects that act like invisible springs or magnets between the dancers.

Summary in One Sentence

The author developed a new way to watch a single fast particle calm down inside an atomic nucleus, found that it happens incredibly fast (too fast for our current theories), and realized that we likely need to rethink how we calculate the "friction" between particles in nuclear matter.

Technical Summary

Here is a detailed technical summary of the paper "Relaxation of a single-particle excitation in a Fermi system within the diffusion approximation of kinetic theory" by Sergiy V. Lukyanov.

1. Problem Statement

The study addresses the time evolution and relaxation mechanisms of non-stationary states in Fermi systems, specifically atomic nuclei. While the Landau–Vlasov kinetic equation is the standard tool for describing such systems, the collision integral is analytically complex. The diffusion approximation simplifies this by modeling collisions as effective diffusion and drift processes in momentum space.

However, previous studies (e.g., Refs. [3, 12]) have reported discrepancies in relaxation timescales:

• Theoretical estimates for the two-particle relaxation time ($\tau_{coll}$) in nuclei are typically $10^{-23}$to$10^{-22}$ s.
• Previous calculations using the diffusion approximation yielded effective relaxation times ($\tau_{eq}$) around $10^{-24}$ s, an order of magnitude too short.
• Crucially, previous models treated the system as a whole, failing to distinguish between the relaxation of the single-particle excitation (the "peak") and the relaxation of the nuclear core (the background distribution).

The primary objective of this work is to numerically solve the nonlinear diffusion equation to separate these contributions, analyze their individual relaxation times, and investigate the influence of kinetic coefficients (diffusion $D_p$ and drift $K_p$) on the characteristic timescales.

2. Methodology

Theoretical Framework

The study employs the Landau–Vlasov kinetic equation within the diffusion approximation. The collision integral is approximated as:
$$\frac{\partial f}{\partial t} = -\nabla_p \cdot \left[ K_p f(1-f) \frac{p}{m} + f^2 \nabla_p D_p \right] + \nabla_p^2 [f D_p]$$
where $f(\vec{p}, t)$ is the Wigner distribution function.

• System Model: A spherical atomic nucleus modeled as infinite nuclear matter (spatially homogeneous, $f=f(p,t)$).
• Coefficients: Constant kinetic coefficients $D_{p,0}$ and $K_{p,0}$ are used. Their ratio determines the equilibrium temperature: $T_{eq} = -D_{p,0}/K_{p,0}$.
• Parameters: The study uses a renormalized drift coefficient to ensure $T_{eq} = 4$ MeV (consistent with previous phenomenological models), while keeping the diffusion coefficient fixed at values derived from Ref. [14].

Separation of Contributions

A novel methodological contribution is the mathematical separation of the distribution function deviation ($\delta f$) into two components:

1. Core Relaxation ($\delta f_0$): The deviation of the step-like background distribution from its own equilibrium.
2. Single-Particle Relaxation ($\delta f_1$): The deviation of the specific single-particle excitation (a peak above the Fermi surface) from its conditional equilibrium.

The total deviation is expressed as $\delta f_{tot} = \delta f_0 + \delta f_1$.

Numerical Implementation

• Equation: The resulting nonlinear diffusion equation (Eq. 7) is solved numerically.
• Initial Conditions:
  • Step distribution: A sharp Fermi step function.
  • Single-particle excitation: A step function plus a localized peak representing a nucleon with momentum $p > p_F$.
• Relaxation Time Definition: The effective relaxation time $\tau_n$ is calculated using the area under the root-mean-square (RMS) deviation curve:
  $$\tau_n = \frac{1}{\Delta_n(0)} \int_0^{t_{max}} \Delta_n(t) dt$$
  where $\Delta_n(t) = \sqrt{\int [\delta f_n(p,t)]^2 d\vec{p}}$. This definition captures the effective timescale of the entire evolution, not just the asymptotic limit.

3. Key Results

A. Relaxation Timescales

• Total System ($\tau$): For a finite nucleus, the total relaxation time is approximately $\tau \approx 10^{-24}$ s. This is significantly shorter than the expected two-particle collision time ($\tau_{coll} \approx 10^{-23} \text{--} 10^{-22}$ s).
• Single-Particle Excitation ($\tau_1$): The relaxation time of the isolated single-particle peak is shorter than the total system relaxation time ($\tau_1 < \tau$). This confirms that single-particle collisions drive the core relaxation, making the overall system relaxation a composite of faster single-particle processes.
• Core Relaxation ($\tau_0$): The relaxation of the step-like background is slightly longer than the single-particle peak but still in the $10^{-24}$ s range.

B. Dependence on System Parameters

• Excitation Energy ($E_{ex}$):
  • Total System ($\tau$): Decreases non-linearly as $E_{ex}$ increases.
  • Single-Particle ($\tau_1$): Increases as $E_{ex}$ increases. This is physically intuitive: a nucleon further from the Fermi surface takes longer to thermalize.
• Mass Number ($A$):
  • As $A$ increases, the total relaxation time $\tau$ increases, approaching the limit of infinite nuclear matter ($\tau_0 \approx 8.9 \times 10^{-24}$ s).
  • Conversely, the single-particle relaxation time $\tau_1$ decreases with increasing $A$. This is attributed to the width of the single-particle peak being inversely proportional to $A$ (Eq. 35); a narrower peak relaxes faster.

C. Influence of Kinetic Coefficients

The study introduced a scaling parameter $y$ to vary $D_p$ and $K_p$ linearly ($D \to D/y, K \to K/y$) while keeping $T_{eq}$ constant.

• Increasing $y$ (reducing coefficients) increases relaxation times.
• To achieve relaxation times consistent with the expected $\tau_{coll} \approx 10^{-22}$ s, the kinetic coefficients would need to be reduced by a factor of $\sim 20$.
• To reach $\tau_1 \approx 5 \times 10^{-23}$ s, a reduction factor of $\sim 200$ is required.

4. Significance and Conclusions

Main Conclusion:
The paper identifies a significant discrepancy between the relaxation times derived from the standard diffusion approximation (using coefficients from previous literature) and the physically expected two-particle relaxation times in nuclei. The calculated times ($\sim 10^{-24}$ s) are too short by an order of magnitude.

Key Contributions:

1. Decomposition of Relaxation: The work successfully separates the relaxation of the single-particle excitation from the nuclear core, revealing that the single-particle component relaxes faster than the background.
2. Parameter Sensitivity: It demonstrates that the relaxation timescale is strictly governed by the magnitude of the diffusion and drift coefficients.
3. Model Limitations: The inability to reproduce the correct timescale without drastically altering the kinetic coefficients (which contradicts previous microscopic calculations) suggests that the diffusion approximation itself may be insufficient or that the current estimates of $D_p$ and $K_p$ are inaccurate for this specific regime.

Implications:
The results imply that the standard diffusion approximation, while useful for qualitative descriptions, may fail to capture the correct timescales for single-particle relaxation in finite Fermi systems without further refinement. The author suggests that future work must re-examine the assumptions of the diffusion approximation, potentially incorporating quantum effects or refining the interparticle interaction potentials to resolve the discrepancy in kinetic coefficients.