Strong-field ionization of He by elliptically polarized light in attoclock configuration

By solving the 3D time-dependent Schrödinger equation for strong-field ionization of helium, this study supports experimental angular offset values derived from non-adiabatic field intensity calibration, thereby contradicting the adiabatic calibration results and the semiclassical conclusions of Boge et al. regarding tunneling time extraction in attoclock measurements.

The Gist

Imagine you are trying to time a very fast event, like a bullet leaving a gun barrel. But instead of a gun, you have an atom (Helium), and instead of a bullet, you have an electron. You want to know exactly how long it takes for that electron to "tunnel" (a quantum trick where it passes through a barrier it shouldn't be able to cross) and escape the atom.

This is the challenge of Attoclock science. It's like trying to measure the blink of an eye, but the eye is an electron and the blink happens in attoseconds (one quintillionth of a second).

Here is a simple breakdown of what this paper does, using some everyday analogies:

1. The Setup: The Rotating Flashlight

In these experiments, scientists shine a laser on the helium atom. But they don't just use a straight beam; they use elliptically polarized light.

• The Analogy: Imagine a lighthouse beam that doesn't just sweep back and forth, but spins in a circle (or an oval) very quickly.
• The Goal: As the electron escapes the atom, this spinning laser beam acts like a giant, invisible hand that pushes the electron. The direction the electron flies off depends on exactly when it escaped.
  • If it escapes at 12:00, it gets pushed to the right.
  • If it escapes at 12:01, it gets pushed slightly lower.
• By measuring the final angle of the electron, scientists can work backward to figure out the exact moment it left the atom.

2. The Mystery: The "Offset" Angle

According to the simplest theory (called the Strong Field Approximation), if the electron escapes exactly when the laser is strongest, it should fly off in a perfectly predictable direction (let's say, straight up).

However, in real experiments, the electron flies off at a slightly different angle. It's "offset" by a few degrees.

• The Question: Why is there this offset?
  • Theory A: Maybe the electron takes a tiny, finite amount of time to tunnel through the barrier. This "tunneling time" causes the laser to rotate a bit more before the electron is free, changing its angle.
  • Theory B: Maybe the electron is being pulled back slightly by the positive charge of the atom's core (like a magnet) as it leaves, which bends its path.

3. The Controversy: Two Different Maps

A previous study (by Boge et al.) tried to solve this. They measured the angle and tried to calibrate their "map" (the laser intensity) in two different ways:

1. The "Slow" Map (Adiabatic): Assumes the electron moves slowly and the laser field doesn't change much while it tunnels.
2. The "Fast" Map (Non-Adiabatic): Assumes the electron moves so fast that the laser field changes significantly during the tunneling process.

The previous researchers found that their data looked better with the "Slow" Map. They concluded that the tunneling time was effectively zero and the offset was just due to the atom's magnetic pull.

4. The New Study: The Supercomputer Simulation

The authors of this paper (Ivanov and Kheifets) decided to run a massive, ultra-precise simulation on a supercomputer. Instead of using a simplified model or making guesses about how the electron behaves, they solved the fundamental equations of quantum mechanics (the Schrödinger equation) from scratch.

• The Analogy: If the previous study was like a weather forecaster using a simple rule of thumb ("If it's cloudy, it will rain"), this study is like running a billion simulations of air pressure, humidity, and wind speed to predict the storm exactly.

5. The Verdict: The "Fast" Map Wins

When they compared their super-precise simulation results to the experimental data, they found a surprising twist:

• Their simulation did not match the "Slow" Map (Adiabatic) data.
• Their simulation perfectly matched the "Fast" Map (Non-Adiabatic) data.

What does this mean?
It suggests that the previous interpretation might have been wrong. The electron isn't just being pulled by the atom; it seems to be experiencing non-adiabatic effects. In plain English, the electron is moving so fast that the laser field changes while it is tunneling. This changes the rules of the game.

Why This Matters

This paper creates a bit of a headache for the scientific community.

• We have Experiment A (Adiabatic calibration).
• We have Experiment B (Non-adiabatic calibration).
• We have Theory A (TIPIS model).
• We have Theory B (This new supercomputer simulation).

Currently, Theory B matches Experiment B, but Experiment A matches Theory A. Since the new simulation is considered the "gold standard" of accuracy, it implies that the "Slow Map" (Adiabatic) used in previous experiments might be flawed.

The Bottom Line:
The authors are saying, "Our supercomputer says the electron is tunneling in a way that changes the laser field as it goes. If you use the old way of calculating the laser's strength, your results are off. We need to rethink how we measure the time it takes for an electron to escape an atom."

It's a reminder that in the quantum world, things are rarely as simple as they seem, and sometimes, the most accurate way to understand a tiny particle is to let a giant computer do the heavy lifting.

Technical Summary

1. Problem Statement

The paper addresses a significant controversy in attosecond science regarding the measurement of tunneling time in strong-field ionization. Specifically, it investigates the interpretation of recent "attoclock" experiments on Helium (He) conducted by Boge et al. [PRL 111, 103003 (2013)].

• The Core Issue: Attoclock measurements rely on mapping the ionization instant to the final angle of the photoelectron momentum vector. The angular offset ($\theta_m$) from the predicted direction is used to extract the tunneling time.
• The Calibration Dilemma: To interpret the data, the local electric field intensity must be calibrated in situ. This calibration depends on the Keldysh parameter ($\gamma$).
  • Adiabatic Hypothesis: Assumes the electron exits the tunnel with zero momentum. This leads to one set of experimental offset values.
  • Non-Adiabatic Hypothesis: Assumes the electron exits with non-zero momentum (relevant as $\gamma$ increases). This leads to a different set of offset values.
• The Conflict: Boge et al. used a semi-classical model (TIPIS) and found qualitative agreement with the adiabatic calibration, concluding that tunneling time vanishes. However, the authors of this paper argue that this conclusion is contradictory because the TIPIS model predictions for the non-adiabatic case are qualitatively different from the experiment, while the adiabatic TIPIS predictions are numerically closer to the non-adiabatic experimental data. This creates a "triad" of contradictions between the two theories (Adiabatic vs. Non-adiabatic TIPIS) and the experiment.

2. Methodology

The authors perform a rigorous, ab initio numerical simulation to resolve this ambiguity without relying on semi-classical approximations or adjustable parameters.

• Governing Equation: They solve the full 3D Time-Dependent Schrödinger Equation (TDSE) for a Helium atom under the Single Active Electron (SAE) approximation.
  $$i\frac{\partial\Psi}{\partial t} = [\hat{H}_{atom} + \hat{H}_{int}(t)]\Psi$$
• Interaction Hamiltonian: The interaction with the elliptically polarized laser field is treated using both length and velocity gauges to ensure gauge invariance. The velocity gauge was found to converge much faster for high intensities.
• Laser Parameters: The simulation mimics the Boge et al. experiment:
  • Ellipticity ($\epsilon$) = 0.87.
  • Wavelength ($\lambda$) = 735 nm ($\omega$ = 1.69 eV).
  • Pulse duration ($T_1$) = 6 optical cycles (with some checks at 3 cycles).
  • Intensity range: $1.0 \times 10^{14}$ to $2.25 \times 10^{14}$ W/cm$^2$.
• Numerical Technique:
  • The wavefunction is expanded in partial waves ($L_{max}$ up to 70 depending on intensity).
  • The radial part is discretized on a spatial grid within a box.
  • Time propagation uses a matrix iteration method.
  • Ionization amplitudes are obtained by projecting the final wavefunction onto incoming scattering states.
• Validation: The authors performed extensive convergence tests regarding partial wave expansion ($L_{max}$), time step size ($\Delta t$), box size ($R_{max}$), and Carrier Envelope Phase (CEP). The estimated error margin for the angular offset is $\pm 1^\circ$.

3. Key Contributions

• First-Principles Calculation: The study provides the first fully quantum mechanical, parameter-free calculation of the angular offset $\theta_m$ for He in the specific attoclock configuration, avoiding the approximations inherent in the Strong Field Approximation (SFA) or semi-classical models like TIPIS.
• Resolution of the Calibration Controversy: By comparing their TDSE results against the two sets of experimental data (calibrated adiabatically vs. non-adiabatically), the authors provide a definitive theoretical benchmark.
• Analysis of CEP Effects: The paper investigates the impact of the Carrier Envelope Phase (CEP) on the momentum distribution, explaining asymmetries observed in experiments and confirming that the angular offset is robust when corrected for CEP-induced drifts.

4. Results

• Angular Offset ($\theta_m$): The TDSE calculations predict a non-zero angular offset relative to the Strong Field Approximation (SFA) prediction ($\theta_m = 0$).
• Agreement with Non-Adiabatic Data: The calculated offset values show clear quantitative agreement with the set of experimental data calibrated under the non-adiabatic hypothesis.
• Disagreement with Adiabatic Data: The calculations strongly disagree with the experimental data calibrated under the adiabatic hypothesis.
• TIPIS Model Failure:
  • The semi-classical TIPIS model based on the adiabatic hypothesis qualitatively matches the adiabatic experimental trend but is numerically closer to the non-adiabatic data.
  • The TIPIS model based on the non-adiabatic hypothesis predicts an offset that increases with intensity, which is qualitatively opposite to the experimental trend (where offset decreases or stays flat).
• Conclusion on the Triad: The authors conclude that one component of the "triad" (Adiabatic Experiment, Non-Adiabatic Experiment, TIPIS Theory) must be flawed. Since the ab initio TDSE supports the non-adiabatic experimental calibration, the TIPIS model is likely inaccurate for extracting tunneling times in this regime.

5. Significance

• Implications for Tunneling Time: If the agreement between the TDSE and the non-adiabatic calibration is not accidental, it implies that non-adiabatic effects are significant in strong-field ionization of Helium and cannot be ignored.
• Critique of Current Interpretations: The findings challenge the conclusions of Boge et al., suggesting that their claim of "vanishing tunneling time" based on the adiabatic calibration may be an artifact of an incorrect calibration method or an inaccurate theoretical model (TIPIS).
• Methodological Impact: The paper demonstrates that high-precision, fully quantum mechanical TDSE calculations are necessary to interpret attoclock data correctly, as semi-classical models may fail to capture the complex interplay between the ionic potential, non-adiabatic dynamics, and field intensity calibration.
• Future Directions: The results suggest that the interpretation of tunneling time measurements reported in recent literature requires re-evaluation, potentially indicating a finite, non-zero tunneling time influenced by non-adiabatic dynamics.