Ab Initio Simulation of Femtosecond Time-Resolved Multi-Pulse Spectroscopies applied to the Heptazine$\cdots$H$_2$O Complex

This paper generalizes the quasi-classical doorway-window methodology to simulate multi-pulse spectroscopies, demonstrating through *ab initio* simulations of the heptazine$\cdots$H$_2$O complex that pump-stimulated experiments provide significantly richer insights into ultrafast radiationless relaxation dynamics than conventional pump-probe and 2D techniques.

The Gist

Imagine you are trying to understand how a complex machine works, like a high-performance car engine, but you can't take it apart. All you can do is shine a flashlight on it and watch how it moves.

This paper is about a team of scientists who developed a new, super-powerful way to "flashlight" a tiny molecule called Heptazine (specifically, a Heptazine molecule holding hands with a water molecule). They wanted to see exactly how this molecule reacts when hit by light, which is crucial for understanding how we might use similar materials to split water and create clean hydrogen fuel.

Here is the breakdown of their work using simple analogies:

1. The Problem: The "One-Flash" Camera

For a long time, scientists used a technique called Pump-Probe spectroscopy. Think of this like taking a photo of a speeding bullet.

• The Pump: You hit the molecule with a bright flash of light (the "pump") to wake it up.
• The Probe: A split-second later, you hit it with a second flash (the "probe") to take a picture.
• The Limit: This gives you a snapshot. You see where the molecule is, but you miss the messy, chaotic details of how it got there. It's like trying to understand a dance by only looking at the start and end positions.

2. The Solution: The "Multi-Pulse" Movie

The scientists in this paper invented a way to use multiple pulses of light to create a "movie" instead of just a snapshot. They focused on two new techniques:

• Pump-Push-Probe (PPP): Imagine you hit the molecule with a flash (Pump). Then, just as it starts to settle down, you hit it with a second, different flash (the Push). This "push" kicks the molecule into a higher, more energetic state. Finally, you take a picture (Probe).

  • The Analogy: It's like pushing a child on a swing. The first push gets them moving. The "push" pulse hits them right when they are at the bottom of the swing, sending them even higher. By watching how they react to that second push, you learn things about the swing's mechanics that a single push couldn't reveal.

• Pump-2D Spectroscopy (P-2D): This is even more complex. It uses five pulses of light.

  • The Analogy: Imagine a 3D MRI scan of a moving object. Instead of just seeing the object from one angle, this technique lets you see the molecule from every angle simultaneously, separating the "shape" of the molecule from the "speed" of its movement. It creates a rich, multi-dimensional map of the molecule's energy.

3. The Simulation: The "Virtual Lab"

Since these molecules are too small to see with a microscope and the events happen in femtoseconds (one quadrillionth of a second—faster than a camera can capture), the scientists couldn't just run an experiment in a lab.

Instead, they built a virtual reality simulation on a supercomputer.

• They used "Ab Initio" methods, which means they calculated the physics from the very bottom up (using the laws of quantum mechanics) without guessing.
• They ran 300+ virtual "movies" (trajectories) of the molecule moving.
• They used a clever shortcut called the "Doorway-Window" method.
  • The Analogy: Imagine you want to know what happens in a crowded room. Instead of tracking every single person (which is too hard), you look at the "Doorway" (where people enter) and the "Window" (what they look like when they leave). By mathematically connecting the two, you can figure out exactly what happened inside the room without needing to film every single second. This made the complex math fast enough to run on a computer.

4. What They Found

When they applied this new "multi-flash" technique to the Heptazine-Water complex, they discovered things that the old "single-flash" method missed:

• The Hidden Dance: The molecule has a "bright" state (easy to see) and "dark" states (hard to see). The old method mostly saw the bright state fading away. The new method revealed that the energy was actually hopping through a series of invisible "dark" states before settling down.
• The Hot Potato: When the "Push" pulse hit the molecule, it didn't just sit there; it got incredibly "hot" (vibrationally excited). The new simulations showed exactly how this heat spread out and how the molecule tried to cool down.
• The Fuel Connection: This is important because Heptazine is a building block for materials used to split water into hydrogen fuel. Understanding exactly how these molecules handle energy and protons (hydrogen atoms) helps engineers design better, more efficient solar fuels.

The Bottom Line

This paper is like upgrading from a still camera to a high-speed, 3D, multi-angle movie camera for the microscopic world.

By using advanced computer simulations, the scientists showed that by hitting molecules with a sequence of carefully timed light pulses, we can see the invisible, ultra-fast steps of chemical reactions. This helps us understand how nature handles energy, which is a giant step forward in our quest for clean, renewable energy.

Technical Summary

1. Problem Statement

The research addresses the limitations of conventional time-resolved spectroscopy (such as standard pump-probe and 2D spectroscopy) in capturing the ultrafast, complex nonadiabatic dynamics of molecular chromophores, specifically the heptazine (Hz) molecule in a hydrogen-bonded complex with water (Hz· · · H2O).

• Context: Heptazine is a fundamental building block of graphitic carbon nitride ($g\text{-C}_3N_4$), a material of interest for photocatalysis. However, its precise reaction mechanisms (e.g., proton-coupled electron transfer, PCET) and ultrafast relaxation pathways are difficult to characterize due to the short lifetimes of excited states and the "dark" nature of intermediate states.
• Limitation: Conventional pump-probe (PP) and 2D spectroscopy often fail to resolve the full complexity of relaxation cascades involving multiple excited states, conical intersections, and the re-population of long-lived states from higher-energy manifolds.
• Goal: To develop and apply advanced multi-pulse spectroscopic simulation protocols (Pump-Push-Probe and Pump-2D) to extract richer dynamical information than standard techniques allow.

2. Methodology

The authors employed a quasi-classical Doorway-Window (DW) approximation combined with ab initio on-the-fly nonadiabatic molecular dynamics.

• Electronic Structure Calculations:
  • Ground State: Calculated using the second-order Møller-Plesset (MP2) method.
  • Excited States: Calculated using the algebraic-diagrammatic construction of second order (ADC(2)).
  • Basis Set: cc-pVDZ augmented with aug-cc-pVDZ functions on water and relevant nitrogen atoms.
  • Software: TURBOMOLE for electronic structure; ZagHop for trajectory propagation.
• Trajectory Simulations:
  • Initial Conditions: Sampled from a zero-temperature harmonic Wigner distribution based on the MP2 Hessian.
  • Dynamics: 300 trajectories propagated using the Velocity Verlet algorithm with a 0.5 fs time step.
  • Nonadiabatic Transitions: Handled via the Landau-Zener (LZ) surface-hopping algorithm (using Belyaev-Lebedev formula) within a manifold of 8 low-lying excited states.
• Spectroscopic Protocols:
  • Pump-Push-Probe (PPP): A three-pulse sequence.
    • Pump: Excites ground state to $S_1$ ($\hbar\omega_{pu} = 4.29$ eV).
    • Push: Delayed by $T_1 = 100$ fs, re-excites the $S_1$ population to higher excited states (Manifold II, 22 states) with $\hbar\omega_{push} = 2.8$ eV.
    • Probe: Delayed by $T_2$ relative to the push, detecting the subsequent relaxation.
  • Pump-2D (P-2D): A five-pulse sequence (Pump + FWM).
    • Uses phase-coherent pulse pairs for excitation and detection to generate a fifth-order signal.
    • Fourier transformed with respect to coherence times ($\tau_{23}$ and $\tau_{45}$) to spread the signal over two frequency axes.
• Signal Calculation: The DW approximation separates signals into Ground-State Bleach (GSB), Stimulated Emission (SE), and Excited-State Absorption (ESA) components, calculated as averages of doorway distributions and window functions.

3. Key Contributions

1. Generalization of DW Methodology: The authors successfully generalized the quasi-classical DW approximation, previously used for pump-probe, to multi-pulse spectroscopies (specifically 3-pulse PPP and 5-pulse P-2D).
2. First Ab Initio Simulations: This work presents the first ab initio simulations of PPP and P-2D spectra for a molecular system, demonstrating the feasibility of these complex protocols for analyzing ultrafast dynamics.
3. Disentanglement of Signal Components: The simulation allows for the separate analysis of GSB, SE, and ESA contributions, which are experimentally inseparable in the total signal.
4. Insight into "Hot" Ensembles: The study demonstrates how multi-pulse techniques can probe the dynamics of "hot" excited-state ensembles (re-excited from $S_1$) rather than just the ground-state ensemble, relaxing standard Franck-Condon and symmetry selection rules.

4. Key Results

The simulations were performed on the Hz· · · H2O complex with a fixed pump-push delay ($T_1 = 100$ fs) and a variable push-probe delay ($T_2$).

• Pump-Push-Probe (PPP) Results:

  • Stimulated Emission (SE): The PPP SE signal reveals the ultrafast (sub-20 fs) decay of higher excited states (populated by the push pulse) back to the long-lived $S_1$ state. A distinct red-shift and oscillation are observed, reflecting the vibrational cooling and relaxation cascade.
  • Excited-State Absorption (ESA): Unlike standard PP where the ESA signal decays rapidly (~30 fs), the PPP ESA signal remains quasi-stationary. This indicates a continuous repopulation of the bright $1\pi\pi^*$ states ($S_5/S_6$) via a cascade of slower relaxation from higher states populated by the push pulse.
  • Comparison: PPP provides a "movie" of the relaxation cascade that standard PP misses, showing how the system returns to the $S_1$ state from higher manifolds.

• Pump-2D (P-2D) Results:

  • Rich Structural Information: The P-2D SE signals exhibit significantly richer peak structures compared to standard 2D signals.
  • Broadening: The P-2D signals show a much wider distribution in both excitation and emission frequencies. This reflects the broad distribution of nuclear geometries and electronic states in the "hot" $S_1$ ensemble at $T_1 = 100$ fs, whereas standard 2D signals are narrow because they originate from the compact ground-state distribution.
  • Dynamics Visualization: The evolution of the P-2D signals (e.g., peak shifting and shape changes from $T_2 = 20$ to $80$ fs) captures the ultrafast radiationless decay through conical intersections "in action."
  • ESA in P-2D: The P-2D ESA signal is less structured but broader in excitation frequency, reflecting transitions from a diverse set of electronic states and geometries.

5. Significance

• Superior Information Content: The study proves that pump-stimulated experiments (PPP and P-2D) provide much richer information on ultrafast radiationless relaxation dynamics than conventional PP and 2D experiments. They can resolve the complex interplay between bright and dark states and the continuous repopulation of intermediate states.
• Methodological Validation: It validates the quasi-classical DW approximation as a computationally efficient and conceptually simple tool for simulating high-order nonlinear spectroscopies, bypassing the extreme computational cost of calculating high-order polarizations directly.
• Guidance for Experiments: The simulations offer a roadmap for experimentalists. By predicting how pulse parameters (delays, frequencies) affect the signal, these ab initio simulations can assist in the planning and optimization of complex multi-pulse experiments for studying photocatalysts and other molecular systems.
• Mechanistic Insight: The results clarify the dynamics of the Hz· · · H2O complex, showing that while the hydrogen bond has weak effects on the internal conversion dynamics, the multi-pulse techniques are essential for visualizing the full relaxation cascade and the role of higher excited states in the system's photophysics.