Explicit proof of Anderson's orthogonality catastrophe for the one-dimensional Fermi polaron with attractive interaction

This paper provides a fully analytical proof of Anderson's orthogonality catastrophe for the one-dimensional Fermi polaron with attractive interactions by combining Bethe ansatz determinant representations and Cauchy matrix properties to demonstrate that the quasi-particle residue decays algebraically with an exponent determined by the Fermi-edge phase shift.

The Gist

The Big Picture: The "Orthogonality Catastrophe"

Imagine a huge, perfectly organized dance floor filled with $N$ dancers (the Fermi gas). They are all spinning in perfect sync, following a strict rhythm. This is the "ground state" of the system—a calm, orderly party.

Now, imagine a new guest arrives (the impurity or polaron). This guest wants to join the dance, but they are a bit clumsy. They bump into people, change the rhythm, and force everyone to adjust their steps to avoid a collision.

The Question: If you look at the dance floor before the guest arrived and after the guest arrived, do the two scenes look similar?

The Answer (The Catastrophe): In the world of quantum mechanics, the answer is a hard NO. Even if the guest only bumps into a few people, the entire crowd shifts so drastically that the "before" and "after" scenes are completely unrecognizable. In fact, as the dance floor gets infinitely large, the two scenes become mathematically orthogonal—meaning they have zero in common. This is called Anderson's Orthogonality Catastrophe.

The Specific Problem: The "Attractive" Guest

For a long time, physicists knew this happened if the guest was a "static" obstacle (like a heavy boulder that couldn't move). But what if the guest is a Fermi Polaron? This is a particle that has a finite mass (it can move) and interacts with the crowd.

There are two types of interactions:

1. Repulsive: The guest pushes the dancers away.
2. Attractive: The guest pulls the dancers close.

The "Repulsive" case was already solved. The "Attractive" case is much trickier. Why? Because when the guest is attractive, they don't just push the crowd; they actually grab two dancers and form a tight little trio (a bound state). This changes the math completely, making it hard to prove that the "Catastrophe" still happens.

What This Paper Does: The "Math Detective" Work

The author, Giuliano Orso, provides a rigorous, step-by-step proof that the Catastrophe does happen even in this tricky "attractive" scenario.

Here is how he did it, using simple analogies:

1. The Two Lists (The Determinants)

To prove the dance floors are different, the author had to compare two massive lists of numbers (matrices):

• List A (The Norm): How "big" the new dance floor is.
• List B (The Overlap): How much the new dance floor looks like the old one.

The "Quasi-particle residue" ($Z$) is the ratio of these two lists. If $Z$ is 1, they look the same. If $Z$ is 0, they are totally different.

2. The Magic Tool: Cauchy Matrices

Calculating these lists for millions of dancers is impossible with normal math. The author used a special mathematical tool called Cauchy Matrices.

• Analogy: Imagine trying to count every possible handshake in a stadium. It's impossible. But if you realize the seating chart follows a specific pattern (like a grid), you can use a shortcut formula to get the answer instantly.
• The author realized the "dance floor" math followed this specific pattern, allowing him to use a shortcut to calculate the size of the lists.

3. The Result: The Algebraic Decay

The author calculated what happens when the number of dancers ($N$) goes to infinity.

• He found that the "similarity" ($Z$) doesn't just drop a little; it drops to zero following a specific power law: $Z \sim 1/N^\theta$.
• The exponent $\theta$ (the Greek letter theta) tells us how fast the similarity vanishes.
• The Surprise: The author proved that $\theta$ depends only on how much the guest changes the rhythm of the dancers at the edge of the crowd (the "Fermi edge"). It doesn't matter that the guest formed a tight trio with two dancers; the "Catastrophe" still happens exactly as predicted.

Why This Matters

1. It's a "Proof," not just a Guess: Before this, scientists had to use supercomputers to simulate the dance floor and guess that the Catastrophe happened. This paper provides a pure mathematical proof that it must happen, no matter how you look at it.
2. It Unifies the Rules: It shows that whether the guest pushes (repulsive) or pulls (attractive), the fundamental rule of quantum mechanics remains the same: Adding a single impurity to a sea of fermions destroys the memory of the original state.
3. Real World Applications: This helps us understand:
   • X-ray absorption: Why metals absorb light in specific ways.
   • Quantum Computers: How errors (impurities) spread through a system.
   • Ultracold Gases: Scientists use lasers to create these "dance floors" in labs. This paper helps them predict exactly what they will see when they introduce a single atom into the mix.

The "Takeaway" Metaphor

Think of the Fermi gas as a perfectly synchronized flash mob.

• Before: Everyone is doing the exact same move.
• After: One person (the impurity) joins and starts a new move.
• The Result: Because everyone is connected, that one new move ripples through the entire crowd. By the time the song ends, the original flash mob is gone. The new formation is so different that if you tried to overlay the "Before" video on the "After" video, they would cancel each other out completely.

This paper proves mathematically that no matter how you try to hide the new move (even if the new guy grabs a partner), the flash mob will always be completely transformed.

Technical Summary

1. Problem Statement

The paper addresses the Anderson Orthogonality Catastrophe (AOC) in the context of the one-dimensional (1D) Fermi polaron with attractive contact interactions.

• System: A system of $N$ spin-up fermions interacting with a single spin-down fermion (impurity) via an attractive delta-function potential. The model is integrable (Yang-Gaudin model) and solvable via the Bethe Ansatz.
• The Phenomenon: AOC describes the vanishing overlap between the ground state of a non-interacting Fermi gas ($|\psi_{NI}\rangle$) and the ground state of the same system with an impurity ($|\psi\rangle$) in the thermodynamic limit ($N, L \to \infty$).
• The Challenge: While AOC is well-understood for static impurities and repulsive 1D Fermi polarons, the attractive case is theoretically more complex. The attraction leads to the formation of a two-body bound state, causing two Bethe-ansatz quasi-momenta to become complex. This invalidates previous "shortcuts" (like the Edwards representation) used for repulsive cases, necessitating the use of the more complex Takahashi representation for the many-body wavefunction.
• Goal: To provide a fully analytical proof that the quasi-particle residue $Z$ (the squared overlap) decays algebraically ($Z \sim N^{-\theta}$) and to derive the exact expression for the Anderson exponent $\theta$ and the prefactor $W$.

2. Methodology

The author employs a rigorous analytical approach combining Bethe Ansatz solutions with linear algebra properties of specific matrix structures.

• Quasi-Particle Residue Definition:
  The residue is defined as:
  $$Z = \frac{|\langle \psi_{NI} | \psi \rangle|^2}{\langle \psi_{NI} | \psi_{NI} \rangle \langle \psi | \psi \rangle} = \frac{|\det(V)|^2}{\det(S)}$$
  where $S$ is the norm matrix and $V$ is the overlap matrix, both of size $N \times N$.

• Takahashi Representation:
  The wavefunctions are expressed using the Takahashi representation, which accounts for the complex quasi-momenta associated with the bound state. This allows $Z$ to be written entirely in terms of determinants of matrices $S$ and $V$.

• Cauchy Matrix Properties:
  A central mathematical tool is the identification of the overlap matrix $V$ (and its core component $C$) as a Cauchy matrix (entries of the form $1/(x_i + y_j)$). The author utilizes known determinant formulas and inverse matrix properties of Cauchy matrices to decompose the problem.

• Asymptotic Analysis:
  The core of the proof involves deriving the leading asymptotics of $\ln(\det S)$ and $\ln(\det V)$ in the thermodynamic limit. The author expands these determinants into sums, converts them into integrals where appropriate, and carefully handles subleading logarithmic corrections using the Euler-Maclaurin formula and properties of the Digamma function ($\psi$) and Hurwitz Zeta function.

3. Key Contributions and Derivations

A. Norm Matrix ($S$) Asymptotics

The determinant of the norm matrix $S$ is derived explicitly. The author shows that in the large-$N$ limit:
$$\ln \det(S) \simeq \ln(k_F^{2N}) + N a - \ln N$$
where $a$ is a coefficient depending on the interaction parameter $\alpha = g'/k_F$, derived analytically as:
$$a(\alpha) = 2\alpha \arctan(1/\alpha) - 2 + \ln(1+\alpha^2)$$

B. Overlap Matrix ($V$) Asymptotics

This is the most complex part of the work. The determinant is factored as $\det(V) = (g'/\pi)^N (1 + \mathbf{v} \cdot \boldsymbol{\beta}) \det(C)$.

1. Vector Scaling: The scalar product term $(1 + \mathbf{v} \cdot \boldsymbol{\beta})$ scales as $N^{1+\xi_N}$, where $\xi_N$ is the reduced phase shift at the Fermi edge.
2. Cauchy Determinant: The determinant of the Cauchy matrix $C$ is decomposed into four terms ($F_1, F_2, F_3, F_4$).
   • $F_1$ converges to a constant.
   • $F_2, F_3, F_4$ are analyzed to extract linear ($N$) and logarithmic ($\ln N$) terms.
   • Crucial cancellations occur between the linear terms of the norm and overlap matrices.

The final asymptotic form for the overlap determinant is:
$$\ln \det(V) \simeq N \left( \ln(g'/\pi) - \sum a_r \right) - \left( \xi_N^2 + \frac{1}{2} \right) \ln N$$

C. Cancellation of Exponential Terms

A major theoretical achievement is demonstrating that the exponential terms (proportional to $N$) in the numerator and denominator of $Z$ cancel out exactly.

• The sum of coefficients $\sum a_r$ from the overlap matrix exactly matches the coefficient $a$ from the norm matrix (adjusted for interaction constants).
• This cancellation leaves only the algebraic dependence on $N$.

4. Results

The paper derives the final expression for the quasi-particle residue:
$$Z = W N^{-\theta}$$

1. The Anderson Exponent ($\theta$):
The exponent is derived analytically as:
$$\theta = 2\xi_N^2 = \frac{2}{\pi^2} \left[ \arctan\left( \frac{-gm}{2k_F} \right) \right]^2$$
where $\xi_N$ is the scattering phase shift at the Fermi edge.

• Significance: This confirms that $\theta$ depends solely on the scattering phase shift at the Fermi surface, regardless of the presence of the two-body bound state (complex quasi-momenta). The result holds for both attractive and repulsive interactions, with the sign of the interaction strength not affecting the form of the exponent.

2. The Prefactor ($W$):
While the analytical derivation confirms the power-law decay, the exact prefactor $W$ (which depends on the interaction strength) is obtained numerically.

• $W$ decreases monotonically as the attraction increases.
• In the strong attraction limit ($|\alpha| \gg 1$), $W \sim \alpha^{-1}$.

5. Significance and Conclusion

• Rigorous Proof: This work provides the first fully analytical proof of the AOC for the attractive 1D Fermi polaron, resolving previous reliance on numerical fitting or variational ansatzes (which failed to capture the AOC).
• Universality: It confirms the universality of the Anderson exponent in 1D integrable systems, showing it is determined entirely by the Fermi-edge phase shift, even in the presence of bound states.
• Methodological Advancement: The paper successfully adapts the Takahashi representation and Cauchy matrix techniques to handle complex quasi-momenta, offering a blueprint for analyzing other 1D integrable models with bound states or unequal masses.
• Future Directions: The authors suggest extending this analytical framework to calculate the prefactor $W$ analytically (by including $1/N$ corrections to phase shifts) and applying the method to systems with unequal masses.

In summary, the paper bridges the gap between numerical observations and rigorous theory, confirming that the orthogonality catastrophe is an intrinsic feature of the 1D Fermi polaron, driven by the accumulation of phase shifts at the Fermi surface, even when the impurity forms a bound state with a fermion.