Light-induced nonadiabatic photodissociation of the NaH molecule including electron-rotation coupling

Imagine a tiny, dancing pair of atoms: Sodium (Na) and Hydrogen (H). Together, they form a molecule called NaH. In the world of quantum physics, these atoms aren't just sitting still; they are vibrating, spinning, and constantly changing their energy states.

This paper is like a high-speed movie camera recording what happens when we zap this dancing pair with two precise laser pulses. The scientists wanted to see how the molecule breaks apart (dissociates) and, more importantly, how the "spin" of the electrons interacts with the "spin" of the whole molecule.

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

1. The Stage: The Dancing Duo

Think of the NaH molecule as a figure skater.

The Skater's Body: The two atoms (Na and H) are the skater's body. They can stretch and shrink (vibrate) and spin around (rotate).
The Skater's Arms: The electrons are the skater's arms. They can move in different patterns (energy states).
The Rulebook: Usually, physicists use a rule called the Born-Oppenheimer approximation. It's like saying, "The skater's body moves slowly, and the arms move super fast, so let's ignore how the arms affect the body's spin." This works most of the time.

2. The Plot Twist: The Laser "Spotlight"

The researchers hit the molecule with two laser pulses:

The Pump Pulse: This is the "kick." It wakes the molecule up, pushing some of the electrons to a higher energy level. The molecule starts vibrating wildly.
The Probe Pulse: This is the "camera flash" that hits later. It checks on the molecule to see what it's doing.

3. The Magic Trick: Light-Induced Conical Intersections (LICIs)

Here is the coolest part. In the dark (no lasers), two energy levels in a simple molecule like NaH can never touch. They are like two parallel train tracks that never cross.

But when you shine a strong laser on them, the laser acts like a magic wand. It bends the energy tracks so that they actually cross each other!

The Crossing: Where these tracks cross is called a Conical Intersection. It's a "portal" or a "funnel."
The Effect: When the molecule hits this portal, it can instantly jump from one energy track to another. This causes the molecule to break apart (dissociate) very quickly.

4. The Hidden Variable: The "Spin-Spin" Connection

This is what makes this specific paper special.

The Old Way (1D Model): Imagine the molecule is a stick that can only vibrate back and forth. It can't spin. This is a simple model.
The Better Way (2D Model): Now, let the stick spin. The laser pushes the stick to align with the light. This changes how the molecule breaks apart.
The New Discovery (3D Model): The scientists asked, "What if the electrons' spin also talks to the molecule's spin?"
- Imagine the skater's arms (electrons) are spinning in a specific direction, and the skater's body (nuclei) is also spinning. Usually, we think the heavy body doesn't care about the light arms.
- But under intense laser light, the arms start pushing the body. This is the Electron-Rotation Coupling.

5. What Did They Find?

The team ran complex computer simulations to see what happens when they include this "spin-spin" talk.

The Big Picture (Breaking Apart): When they looked at how many molecules broke apart or how fast the pieces flew, the "spin-spin" effect didn't change much. The simple models (1D and 2D) were actually pretty good at predicting this.
The Direction (Where the Pieces Fly): This is where the magic happened. When they looked at which direction the broken pieces flew, the "spin-spin" effect made a huge difference, but only in a very specific spot.
- The Analogy: Imagine throwing a ball. If you ignore the wind (rotation), the ball goes straight. If you account for the wind, it curves.
- The scientists found that when the molecule breaks apart exactly straight along the laser beam (the "Z-axis"), the "spin-spin" coupling acts like a shield. It stops the molecule from breaking apart in that exact direction. Without this coupling, the model predicted it would break apart there.

The Takeaway

This paper is a bit like a detective story. The scientists were looking for a subtle clue (the electron-rotation coupling) in a chaotic crime scene (the breaking molecule).

They found that while this clue doesn't change the number of broken pieces, it does change the direction they fly in, specifically blocking them from flying straight down the laser's path.

In short:

Lasers can create "portals" (intersections) in molecules that let them break apart.
The spinning of the whole molecule matters a lot.
The spinning of the electrons also matters, but mostly when you look very closely at the direction the pieces fly. It's a tiny effect that only shows up when you look at the "straight-ahead" angle, proving that even tiny electrons can push heavy atoms around if the light is strong enough.

Here is a detailed technical summary of the paper "Light-induced nonadiabatic photodissociation of the NaH molecule including electron-rotation coupling."

1. Problem Statement and Motivation

The paper addresses the complex dynamics of molecular photodissociation under strong laser fields, specifically focusing on the NaH molecule. While the Born-Oppenheimer (BO) approximation typically separates electronic and nuclear motions, this approximation breaks down near conical intersections (CIs), where electronic states become degenerate.

The Challenge: In field-free diatomic molecules, CIs cannot naturally form because there is only one vibrational degree of freedom. However, strong laser fields can induce Light-Induced Conical Intersections (LICIs) by coupling the molecular rotation to the electronic states, effectively providing the necessary second degree of freedom.
The Gap: Previous studies often neglected the coupling between the molecular rotational angular momentum ( $K$ ) and the total electronic angular momentum projected onto the molecular axis ( $\Omega = \Lambda + \Sigma$ ), known as $K$ - $\Omega$ coupling. This effect is usually ignored due to the mass difference between nuclei and electrons but may become significant at high laser intensities or in specific angular distributions.
Objective: The authors aim to investigate the ultrafast photodissociation of NaH using a pump-probe scheme, explicitly comparing three levels of theoretical description:
1. 1D Model: Neglects nuclear rotation (rotation is a parameter).
2. 2D Model: Includes nuclear rotation as a dynamic variable (Light-Induced Avoided Crossing/Conical Intersection).
3. 3D Model: Includes both nuclear rotation and the $K$ - $\Omega$ electron-rotation coupling.

2. Methodology

System and Electronic Structure

Molecule: Sodium Hydride (NaH).
States: The model includes the three lowest-lying singlet electronic states: $X^1\Sigma^+$ , $A^1\Sigma^+$ , and $B^1\Pi$ .
Potential Energy Curves (PECs): Calculated using the MOLPRO package at the MRCI/CAS(2/9)/aug-cc-pV5Z level.
Dipole Moments: Permanent dipole moments (PDMs) and transition dipole moments (TDMs) were computed to define the light-matter interaction.

Theoretical Framework

Hamiltonian: The time-dependent nuclear Hamiltonian was formulated in the Wigner $D$ -function basis $|JM\Omega\rangle$ $∣ J M Ω ⟩$ to properly handle the $K$ $K$ - $\Omega$ $Ω$ coupling (Hund's case (a)).
- The interaction term includes both PDMs (diagonal) and TDMs (off-diagonal), with angular dependence determined by Euler angles ( $\theta, \phi, \chi$ ).
Laser Setup: A pump-probe configuration was used.
- Pump Pulse: $\hbar\omega_{pm} = 3.39$ eV (3rd harmonic), $\tau_{pm} = 22$ fs, $I_{pm} = 1 \times 10^{12}$ W/cm². Excites population from $X$ to $A$ .
- Probe Pulse: $\hbar\omega_{pr} = 1.13$ eV, $\tau_{pr} = 30$ fs. Two intensities were tested: $1 \times 10^{11} $W/cm² (linear regime) and$ 1 \times 10^{12}$ W/cm² (multiphoton regime).
- LICI Formation: The probe frequency was tuned to create two LICIs between $X/A$ and two between $A/B$ via Floquet theory (dressed states).
Propagation: The Time-Dependent Schrödinger Equation (TDSE) was solved using the Multiconfigurational Time-Dependent Hartree (MCTDH) method (Heidelberg package).
- Grid: Vibrational coordinate ( $R$ ) discretized with 512 points; Rotational coordinates ( $\theta, \phi, \chi$ ) expanded in Wigner $D$ -functions.
- Absorbing Potential: A Complex Absorbing Potential (CAP) was applied to calculate flux and dissociation probabilities without spurious reflections.

Observables Calculated

Dissociation Probabilities ( $P_{diss}$ ): Total probability of dissociation into channels $X$ , $A$ , and $B$ .
Kinetic Energy Release (KER) Spectra: Energy distribution of the photofragments.
Photofragment Angular Distributions (PAD): The angular probability of fragments relative to the laser polarization axis.

3. Key Contributions

Explicit Inclusion of $K$ - $\Omega$ Coupling: The study provides a rigorous 3D simulation framework that incorporates the coupling between nuclear rotation and electronic angular momentum, a factor often omitted in standard 2D models.
Systematic Comparison of Dimensions: The paper offers a direct comparison of 1D (parameterized rotation), 2D (dynamic rotation), and 3D (rotation + $K$ - $\Omega$ coupling) models to isolate the specific impact of electron-rotation coupling.
Characterization of LICI Dynamics: It elucidates how the position and nature of light-induced degeneracies (LICIs vs. LIACs) depend on the level of theoretical description and the delay time between pump and probe pulses.

4. Key Results

A. Dissociation Probabilities

1D vs. 2D/3D: Significant differences were observed between the 1D model and the 2D/3D models, particularly at longer delay times. The 1D model fails to capture the alignment of molecules induced by the pump pulse, leading to incorrect coupling strengths.
2D vs. 3D: For total dissociation probabilities, the 2D and 3D results were indistinguishable. The inclusion of $K$ - $\Omega$ coupling did not significantly alter the total yield or the timing of dissociation peaks.
Intensity Effects: At higher probe intensities ($10^{12} $W/cm²), saturation effects distorted the peaks, but the linear scaling of the$ B$-state dissociation confirmed the system remained largely within the perturbative regime.

B. Kinetic Energy Release (KER) Spectra

1D vs. 2D: Differences emerged at long delay times due to the blurring of the wavepacket and rotational alignment effects.
2D vs. 3D: The relative difference between 2D and 3D KER spectra was minimal (< 2-3%), even at high intensities. The $K$ - $\Omega$ coupling does not measurably alter the energy distribution of the fragments.

C. Photofragment Angular Distributions (PAD)

Critical Finding: This is the only observable where the 3D model showed a distinct difference from the 2D model.
- Small Angles ( $\theta \approx 0^\circ$ ): In the 2D model, the PAD for the $B$ state showed non-zero dissociation at $\theta = 0^\circ$ due to the alignment of the wavepacket and the location of the LICI ( $\theta=0$ ).
- 3D Suppression: The inclusion of $K$ - $\Omega$ coupling in the 3D model entirely suppressed dissociation at $\theta = 0^\circ$ . The coupling modifies the effective potential surfaces, preventing the wavepacket from accessing the dissociation channel along the polarization axis.
Integration Effect: When the angular distribution is integrated over the azimuthal angle $\phi$ (to match experimental total angular distributions), the differences between 2D and 3D models vanish. This is because the integration weights (proportional to $\sin\theta$ ) diminish the impact of the small-angle region where the discrepancy exists.

5. Significance and Conclusion

Validation of Approximations: The study confirms that for total dissociation yields and energy spectra, the computationally cheaper 2D model (which includes nuclear rotation but neglects $K$ - $\Omega$ coupling) is sufficient. The 1D model is insufficient for accurate dynamics.
Angular Sensitivity: The $K$ - $\Omega$ coupling is crucial for accurately predicting angular distributions, specifically in the forward/backward scattering directions ( $\theta \approx 0^\circ$ ). Neglecting this coupling leads to qualitative errors in the angular profile of the photofragments.
Physical Insight: The results suggest that while electron-rotation coupling is a subtle effect, it acts as a "filter" for specific angular channels. It does not change the overall energy or probability of dissociation but fundamentally alters where the fragments are ejected.
Future Directions: The authors note that the mechanism behind the suppression of the 3D model at small angles requires further investigation, particularly as laser intensities increase and the coupling becomes steeper.

In summary, this paper establishes that while electron-rotation coupling is negligible for global observables like total yield and KER in NaH, it is a critical factor for high-resolution angular distribution measurements, necessitating the use of 3D models for precise spectroscopic predictions in strong-field molecular dynamics.