Convergent calculations of double ionization of helium: from ($γ$,2e) to (e,3e) processes

This paper resolves conflicting interpretations of double ionization of helium by demonstrating that the original Kheifets et al. calculations remain valid and can be reproduced either by incorporating the second Born approximation or by using a corrected ground state similar to that proposed by Jones and Madison.

The Gist

Imagine you are trying to predict exactly how a billiard ball will scatter when it hits a cluster of two other balls sitting on a table. In the world of atoms, this is what happens when a fast-moving electron smashes into a helium atom, knocking out its two electrons. This is called an (e,3e) process (one electron in, three electrons out).

For a long time, scientists have been arguing about the best way to calculate the outcome of this cosmic billiard game. This paper by Kheifets and Bray is essentially them saying, "We've re-run the numbers, checked our assumptions, and we're sticking with our original prediction."

Here is the breakdown of the drama, the arguments, and the resolution, using some everyday analogies.

The Setup: The Cosmic Billiard Table

The authors are using a super-precise mathematical tool called the CCC method (Convergent Close-Coupling). Think of this as a high-end video game engine. The more "pixels" (mathematical details) you add to the engine, the more realistic the physics becomes. The beauty of their engine is that if you keep adding pixels, the picture eventually stops changing—it "converges." Once it converges, you know you have the right answer, regardless of what anyone else says.

The Conflict: The "Wrong" Prediction

A few years ago, the authors used their engine to predict what would happen when a fast electron hits helium. They predicted a certain result. However, when experimentalists (the people actually doing the billiard game in a lab) measured it, the real-world results were much bigger (about 3 to 12 times bigger) than the prediction.

This caused a huge debate. Two other groups of scientists stepped in to explain the discrepancy:

1. The "Bad Physics" Argument (Berakdar): This group said, "Your engine is using the wrong rules! You are using the 'First Born Approximation,' which is like assuming the billiard balls are ghosts that don't really interact until they hit. We think you need to use 'Second Born' rules, which account for the balls 'feeling' each other before they touch."
2. The "Bad Starting Position" Argument (Jones and Madison): This group said, "Your rules are fine, but you started with the wrong setup. You used a model of the helium atom that was too simple. The electrons in helium are very close together, and your model didn't handle that 'cusp' (the tight squeeze) correctly. If you use a better starting model, your old rules will work."

The Investigation: Testing the Theories

Kheifets and Bray decided to settle this by running their own "stress test" on their engine. They didn't just guess; they changed the settings to see what happened.

Test 1: The "Second Born" Rules
They tried adding the "Second Born" complexity (the ghost-interaction idea).

• Result: Nothing changed. The prediction stayed exactly the same.
• Analogy: It's like trying to fix a blurry photo by turning up the contrast. The picture was already sharp; the problem wasn't the contrast. This proved that the "First Born" rules were actually fine for this specific speed of electron.

Test 2: The "Better Starting Position"
They tried using the new, fancy starting model (the Pluvinage state) that Jones and Madison recommended.

• Result: Disaster. When they used this new model, their calculations became wildly inaccurate and didn't match the real-world experiments at all.
• Analogy: It's like trying to drive a car with a brand-new, high-tech GPS, but you accidentally plug in a map of Mars instead of Earth. The car is great, but the starting point is wrong.

Test 3: The "Improved" Starting Position
They then tried a modified version of that fancy model (the Le Sech state), which fixed a specific flaw in the original "Mars map."

• Result: Suddenly, the numbers matched their original, simpler model perfectly.
• Analogy: They fixed the GPS coordinates. Now, both the "Simple Map" and the "Fancy Map" lead to the same destination.

The Verdict: Who Was Right?

The authors found that:

1. The "Ghost" argument was wrong. You don't need complex "Second Born" rules for this specific experiment.
2. The "Bad Map" argument was half-right, but the solution was tricky. The fancy map (Pluvinage) was actually worse than the old one. However, if you tweak the fancy map (Le Sech), it works just as well as the old one.
3. The Original Prediction Stands. Their original calculations were correct. The disagreement with the experiment wasn't because their math was wrong; it was likely due to other subtle factors not fully captured in the comparison, or perhaps the experimental scaling needs re-evaluation.

The Big Picture

The authors conclude that their original approach is robust. They argue that the physics of knocking out two electrons with a fast electron (e,3e) is fundamentally the same as knocking them out with light (γ,2e). Since their method works perfectly for light, it should work for electrons too.

In summary:
Imagine a group of architects arguing over why a bridge they built collapsed.

• One says, "The steel was too weak!"
• Another says, "The foundation was laid on sand!"
• The original builders (Kheifets and Bray) say, "We tested the steel (it's fine) and we tested the foundation (it's fine). We think the bridge is actually standing, and the people saying it collapsed are measuring it wrong."

They proved their bridge is solid by showing that changing the "rules" or the "foundation" didn't break it, but rather confirmed that their original design was the most reliable one all along. They are confident in their science and hope this debate encourages more experiments to settle the score once and for all.

Technical Summary

1. Problem Statement

The paper addresses a significant discrepancy in the theoretical understanding of electron-impact double ionization of helium ($(e,3e)$ processes).

• The Conflict: In 1999, Lahmam-Bennani et al. performed the first absolute measurements of the $(e,3e)$ process. Subsequent theoretical attempts to reproduce these results yielded conflicting conclusions:
  • Berakdar (2000): Achieved agreement with experiment using a Faddeev-type approach but concluded that the First Born Approximation (FBA) was invalid for this process, implying higher-order terms were necessary. This indirectly suggested that previous Convergent Close-Coupling (CCC) calculations by Kheifets et al. (which used FBA) were incorrect.
  • Jones and Madison (2003): Argued that the FBA was valid, but the disagreement in previous CCC results was due to the use of an inaccurate initial ground state (specifically, Hylleraas-type wavefunctions failing to satisfy Kato cusp conditions). They claimed that using a Pluvinage ground state (which satisfies cusp conditions) combined with a 3C final state restored agreement with experiment within the FBA framework.
• The Goal: The authors investigate these claims to determine if the original CCC results were flawed due to the approximation order or the initial state, and to establish a predictive theory for complex collision processes.

2. Methodology

The authors employ the Convergent Close-Coupling (CCC) method, a formalism designed to test convergence by varying input parameters (basis size, state descriptions) until results stabilize.

• Formalism:
  • The total wavefunction is expanded using a complete Laguerre basis.
  • The final state wavefunction ($\Psi_{k_1k_2}$) is treated as a product of a true continuum state and positive energy pseudostates, solving the Lippmann-Schwinger equation to fully account for electron-electron interactions.
  • The projectile-target interaction is treated in both the First Born and Second Born approximations to test the necessity of higher-order terms.
• Initial State Variations:
  • Hylleraas: A standard multi-configuration wavefunction used in earlier CCC work.
  • Pluvinage: A wavefunction satisfying Kato cusp conditions but lacking screening effects.
  • Le Sech: An improved version of the Pluvinage state that includes screening of electron-electron interaction by the nucleus and shielding of the electron-nucleus interaction by the other electron.
• Validation Strategy:
  • Gauge Invariance: The authors calculate Double Photoionization (DPI or $(\gamma,2e)$) cross-sections using three gauges of the electromagnetic operator (Length, Velocity, Acceleration). If the wavefunctions are exact, results in all three gauges should be identical. Discrepancies indicate inaccuracies in the ground state description.
  • Scaling: They compare the $(\gamma,2e)$ and $(e,3e)$ results, noting that in the optical limit (low momentum transfer), the physics should be closely related.

3. Key Contributions

1. Systematic Testing of Approximations: The authors explicitly test the Second Born term for the specific kinematics of the experiment and find it indistinguishable from the First Born term, challenging Berakdar's claim that the FBA is insufficient.
2. Ground State Sensitivity Analysis: They rigorously test the Pluvinage, Le Sech, and Hylleraas ground states. They demonstrate that the Pluvinage state, while satisfying cusp conditions, fails to provide gauge-invariant results for DPI, particularly in the Length gauge.
3. Re-evaluation of the "3C" Final State: The authors challenge the notion that the 3C final state (used by Jones and Madison) is superior to the CCC final state, arguing that the CCC final state is more rigorous for the kinematics considered.
4. Convergence Verification: By using a significantly larger basis set ($N_l = 60-l$) compared to previous work ($N_l = 17-l$), they ensure the results are numerically converged and not artifacts of interpolation or basis truncation.

4. Results

• Gauge Dependence in DPI:
  • Pluvinage State: Shows severe gauge dependence (Length, Velocity, and Acceleration gauges yield vastly different results). The Length gauge, in particular, fails to match experiment.
  • Le Sech State: Significantly reduces gauge variation, showing good agreement with experiment at 20 eV excess energy.
  • Hylleraas State: Yields excellent agreement between all three gauges and experimental data.
• $(e,3e)$ Calculations:
  • When using the Pluvinage ground state with the CCC final state, the results are qualitatively different from the experiment and require a scaling factor of ~6 to match magnitudes.
  • When using the Le Sech or Hylleraas ground states with the CCC final state, the results match the experimental shape and magnitude (after applying scaling factors of 1.8 and 2.2, respectively).
  • Crucial Finding: The First Born and Second Born approximations yield indistinguishable results for this kinematics, refuting the claim that the FBA is the source of the error.
• Comparison with Other Theories:
  • The 3C-based calculations of Jones and Madison (using Pluvinage) and Berakdar (using 4-body Green's functions) agree with each other and the experiment.
  • However, the authors argue this agreement is coincidental or fragile; changing the projectile charge in the 3C model destroys the agreement, whereas the CCC method remains robust.
  • The authors conclude that the success of Jones and Madison's model is likely due to error cancellation (an inaccurate ground state compensating for an approximate final state) rather than physical correctness.

5. Significance

• Validation of the CCC Method: The paper reaffirms the predictive power of the CCC method. The authors maintain that their original calculations (using Hylleraas ground states and FBA) were fundamentally correct, and the discrepancies with early experiments were due to experimental uncertainties or specific kinematic factors, not theoretical flaws.
• Refutation of the "Cusp Condition" Necessity: While satisfying Kato cusp conditions is theoretically desirable, the paper demonstrates that a wavefunction satisfying these conditions (Pluvinage) can still yield poor physical results if it lacks proper screening effects. The "improved" Le Sech state works, but so does the standard Hylleraas state, suggesting the Pluvinage state is not the "magic bullet" proposed by Jones and Madison.
• First Born Validity: The study confirms that for fast projectile electron impact double ionization of helium, the First Born approximation is sufficient, provided the initial and final states are treated with high accuracy (convergent close-coupling).
• Scientific Merit over Agreement: The authors emphasize a philosophy of "predictive theory," arguing that results should be trusted based on convergence and internal consistency (gauge invariance) even if they initially disagree with experiment or other theories.

In conclusion, Kheifets and Bray demonstrate that the original CCC approach to $(e,3e)$ is robust. The apparent success of alternative theories relying on the Pluvinage ground state is attributed to the specific combination of an approximate final state (3C) and a specific ground state, rather than a fundamental failure of the First Born approximation or the CCC method.