Femtosecond Nonadiabatic Confinement of Molecular Dication Yield

By combining experimental observations with ab-initio calculations, this study reveals that ultrafast nonadiabatic relaxation competes with strong-field ionization in ethylene, confining the production of molecular dications to a narrow 15-femtosecond temporal window driven by resonant enhancement during bond expansion.

The Gist

Imagine a molecule of ethylene (a simple gas used to ripen fruit) as a tiny, vibrating trampoline made of two carbon atoms and four hydrogen atoms. Scientists wanted to understand what happens when you hit this trampoline with a super-fast, high-energy "punch" and then immediately follow up with a series of rapid "taps."

Here is the story of their discovery, broken down into everyday concepts:

The Setup: The Punch and the Tap

The researchers used two different types of light to play a game of "pump and probe" with the ethylene molecule:

1. The Pump (The Punch): They hit the molecule with an extreme ultraviolet (XUV) pulse. Think of this as a single, incredibly fast, high-energy punch. It knocks one electron out of the molecule, turning it into a positively charged "cation" (a molecule with a missing piece). This punch is so fast it happens in a fraction of a second (attoseconds).
2. The Probe (The Taps): A few femtoseconds later (a femtosecond is a quadrillionth of a second), they hit the now-charged molecule with a near-infrared laser. This isn't one big hit; it's a series of rapid taps. To knock out a second electron and turn the molecule into a "dication" (a molecule with two missing pieces), the molecule has to absorb several of these taps at once.

The Mystery: The 15-Second Sweet Spot

When they varied the time between the punch and the taps, they found something surprising. They didn't get the most dications immediately after the punch, nor did they get them a long time later. Instead, the number of dications created peaked sharply at a delay of about 15 femtoseconds.

It's as if the molecule has a very specific, tiny window of time where it is perfectly "poised" to accept the second hit. Miss that window by a few femtoseconds, and the result is much lower.

The Mechanism: Stretching the Trampoline

Why does this 15-femtosecond window exist? The paper explains it using a race between two competing forces:

1. The Stretch (Nuclear Dynamics): After the first punch, the molecule starts to vibrate and stretch. Specifically, the bond between the two carbon atoms (the C=C double bond) begins to lengthen, like a rubber band being pulled.

   • As this bond stretches, the energy required to knock out the second electron changes.
   • At a specific stretch length (around 1.4 to 1.5 Angstroms), the molecule enters a "resonant" state. This is like finding the perfect rhythm on a swing; the multiple taps from the laser hit the molecule at just the right moment to knock the second electron out very efficiently. This is called Resonance-Enhanced Multi-Photon Ionization (REMPI).

2. The Fade (Non-Adiabatic Relaxation): However, the excited states of the molecule are unstable. They are like a spinning top that is wobbling; they naturally want to settle down or "relax" into a calmer state very quickly. This relaxation happens on the same ultrafast timescale (around 15–20 femtoseconds).

   • If the molecule relaxes too quickly, it loses the specific energy configuration needed to catch the laser taps efficiently.
   • If the bond hasn't stretched enough yet, the taps aren't efficient either.

The Result: The peak at 15 femtoseconds is the "Goldilocks" moment. It is the exact split-second where the bond has stretched enough to make the laser taps super effective, but the molecule hasn't yet relaxed and lost that special configuration.

The Analogy: The Juggling Act

Imagine a juggler (the molecule) trying to catch a ball (the second electron being knocked out).

• The Punch: The juggler is hit, causing them to spin and stretch their arms out.
• The Taps: A machine starts firing balls at them.
• The Window: For the first few seconds, the juggler is spinning too wildly to catch the balls. Then, their arms stretch out to the perfect length, and they are in the perfect rhythm to catch the balls (the 15 fs peak). But immediately after that, they start to calm down and stop spinning, or their arms collapse, and they can no longer catch the balls as well.

The Takeaway

The paper claims that this experiment reveals a general rule for how molecules behave under intense light: Ultrafast relaxation (calming down) competes with strong-field ionization (getting hit).

The researchers used advanced computer simulations to confirm that this "confinement" of the dication yield to a narrow 15-femtosecond window is caused by the tug-of-war between the bond stretching (which helps the ionization) and the electronic states relaxing (which hurts the ionization).

In short, the molecule doesn't just sit there waiting to be hit; it is constantly moving and changing. The laser only works best when it catches the molecule in a fleeting, specific pose that lasts for only a few femtoseconds.

Technical Summary

Technical Summary: Femtosecond Non-Adiabatic Confinement of Molecular Dication Yield

Problem Statement
Doubly charged molecular ions (dications) serve as unique probes for electronic correlation and electron-nuclear entanglement within parent cations. While time-resolved studies using XUV pump and near-infrared (NIR) probe pulses have revealed ultrafast dynamics in various systems, the specific influence of the multi-photon nature of the NIR probe on the dication yield has often been overlooked due to simulation challenges. In particular, the interplay between ultrafast non-adiabatic relaxation of the intermediate singly-charged cation and the strong-field (or multi-photon) ionization rate required to produce the dication remains to be fully characterized. This study investigates the formation of ethylene dications ($C_2H_4^{++}$) to understand how these competing dynamics confine the dication yield to a narrow temporal window.

Methodology
The research combines experimental pump-probe spectroscopy with advanced ab-initio theoretical modeling:

• Experimental Setup: The team utilized an XUV pump – NIR probe scheme. The XUV pump (generated via High-Harmonic Generation in Xenon or Krypton, photon energy 17–23.3 eV) photoionizes and excites neutral ethylene ($^{13}C^{12}CH_4$) to various cationic states ($D_0$ to $D_4$). A time-delayed NIR probe pulse (central wavelength ~800 nm, intensity $10^{12}$–$10^{13}$ W/cm$^2$) subsequently drives multi-photon ionization to produce the dication. The experiment operates in the multi-photon regime (Keldysh parameter $\gamma \approx 3.5$–6), well below the tunneling ionization threshold. Time-of-flight mass spectrometry was used to detect the $^{13}C^{12}CH_4^{++}$ dication yield as a function of pump-probe delay.
• Theoretical Framework: To interpret the results, the authors employed a two-step simulation approach:
  1. Non-Adiabatic Dynamics: Trajectory surface-hopping (TSH) simulations were used to model the coupled electron-nuclear dynamics of the ethylene cation following XUV ionization. This tracked the population evolution of cationic states ($D_0$–$D_4$) and the stretching of the C=C bond.
  2. Multi-Photon Ionization: The ASTRA (AttoSecond TRAnsitions) package was used to explicitly calculate the multi-photon ionization probabilities of the intermediate cationic states by the NIR probe. This involved solving the time-dependent Schrödinger equation (TDSE) using a transition-density-matrix time-dependent close-coupling formalism. This step accounted for the non-perturbative nature of the probe and the specific geometry (C=C bond length) of the cation at each time step.

Key Results

• Experimental Observation: The dication yield exhibits a distinct peak at a pump-probe delay of approximately 15 fs ($\Delta t = 14.8 \pm 1.7$ fs). The yield is confined within a narrow temporal window and shows no significant dependence on the NIR probe intensity within the tested range. No dication signal is observed with either the XUV or NIR pulse alone.
• Theoretical Reproduction: The combined theoretical model successfully reproduced the observed delay. The calculated dication yield peaks around 10 fs, showing excellent qualitative agreement with the experimental trend.
• Mechanism of Confinement: The peak arises from a competition between two factors:
  1. Resonant Enhancement: As the C=C bond stretches during non-adiabatic relaxation (reaching lengths of ~1.4–1.5 Å between 10–20 fs), the ionization potential decreases, and resonance-enhanced multi-photon ionization (REMPI) pathways become accessible. This increases the ionization rate ($Y_i(t)$) for specific states.
  2. Population Decay: Simultaneously, the electronic population of the excited cationic states ($D_1$, $D_2$) undergoes rapid non-adiabatic relaxation and decay.
     The dication yield is maximized when the ionization rate is enhanced by nuclear motion (bond stretching) before the electronic population of the relevant states ($D_1$ and $D_2$) decays significantly. States like $D_3$ show favorable bond stretching but contribute little due to low residual population.
• State Selectivity: The analysis reveals that the NIR probe is state-selective. While the XUV pulse populates multiple cationic states, the dication yield is dominated by the ionization of the $D_1$ and $D_2$ states via REMPI processes facilitated by the specific C=C bond lengths achieved at the ~15 fs delay.

Significance and Claims
The paper claims to demonstrate a general mechanism where ultrafast non-adiabatic relaxation of a molecular ion competes with its strong-field (multi-photon) ionization rate. This competition results in the "confinement" of the dication yield to a narrow temporal window of a few femtoseconds.

The authors emphasize that their work highlights the necessity of explicitly including both the non-adiabatic dynamics of the intermediate cation and the non-perturbative, multi-photon nature of the probe pulse to accurately interpret time-resolved dication yields. They suggest that this mechanism is likely relevant for other molecular systems, particularly those studied with high-energy XUV or X-ray pumps that simultaneously ionize and excite the system. The study serves as a benchmark for understanding fragmentation dynamics in polyatomic systems following strong-field ionization and motivates future experiments using intense attosecond sources at even shorter wavelengths.