Phase structure of below-threshold harmonics in aligned molecules: a few-level model system

By means of few-level models, this study demonstrates that sub-threshold harmonics in aligned molecules exhibit pronounced phase alternations and polarization behaviors that depend on their energy relative to transition frequencies, thereby enabling the prediction of a mirrored polarization in higher-order harmonics for systems with orthogonal transition dipole moments.

The Gist

Imagine a molecule as a tiny, complex musical instrument, like a violin with two different strings. When you strike it with a powerful laser (the bow), it produces not just a single tone, but an entire orchestra of new, higher tones known as "overtones."

Normally, scientists focus on the loudest, highest tones. Yet this article is interested in the quieter, lower tones that occur just below a certain "threshold" volume. The researchers wanted to understand the timing (phase) and direction (polarization) of these specific tones when the instrument is perfectly aligned.

Here is the breakdown of their discovery using simple analogies:

1. The Two-Level System: A Single Swing

First, the scientists considered a simplified model: a molecule with only two energy states, like a child on a swing.

• The Setup: They pushed the swing with a laser.
• The Discovery: They found a strange rule regarding the timing of the tones the swing produces.
  • Below the "Sweet Spot": When the tones lie below a certain energy level (the transition energy), the timing of the tones flips back and forth. Imagine a drummer beating a rhythm: Left, Right, Left, Right. The "phase" (the start of the beat) switches by 180 degrees (π) for each new tone.
  • Above the "Sweet Spot": Once the tones become higher than this energy level, the flipping stops. It becomes constant, like a drummer who only strikes Left, Left, Left.

Why does this happen?
The article explains this with a mathematical recipe. It is like a chain reaction. If the "recipe" for generating the next tone contains a negative sign, the tone flips its timing. If the sign is positive, it keeps the same timing. The switch occurs exactly when the energy of the tone exceeds the molecule's natural energy gap.

2. The Four-Level System: The Crossed Strings

Next, they built a more complex model to simulate a real molecule. Imagine a molecule to which two of these "swings" (two-level systems) are attached:

• Swing A is oriented horizontally (like an x-axis).
• Swing B is oriented vertically (like a y-axis).
• They are decoupled, meaning they do not talk to each other, but both are struck by the same laser.

The Magic Trick:
Since the two swings have slightly different natural frequencies, the "Sweet Spot" (where the timing flips) occurs at different tones for each swing.

• Low Tones: For the first few tones, both swings are "below" their Sweet Spots. Both flip their timing synchronously. The resulting light points in the same direction as the laser.
• High Tones: Eventually, the tones become high enough that Swing A is "above" its Sweet Spot (constant timing), while Swing B is still "below" its Sweet Spot (flipping timing).
  • Now one swing says "Left," while the other says "Right" (a phase difference of 180 degrees).
  • When these two opposing signals are combined, the resulting light does not just point in the direction of the laser. It reflects or flips to the opposite side.

3. The Real Implication

The article suggests that real molecules (such as certain organic crystals) possessing these two perpendicular "strings" with different energy gaps should exhibit exactly this behavior.

• When illuminated by a laser, the low-energy overtones point in one direction.
• The high-energy overtones (still below the ionization threshold) suddenly point in a mirrored direction.

Summary

Imagine it as a dance floor with two groups of dancers:

1. Group A and Group B dance to the same music.
2. During the slow songs, both dance in sync.
3. During the fast songs, Group A keeps the rhythm constant, but Group B starts dancing backward.
4. When you look at the whole floor, the combined dance movements suddenly flip direction.

The article claims that by observing how the light (the dance) changes in direction and timing, we can learn something about the hidden energy levels and the structure of the molecule, specifically about how its electrons move between bound states without flying off into space. This offers a new way to "see" the internal structure of molecules with light.

Technical Summary

Technical Summary: Phase Structure of Sub-Threshold Harmonics in Aligned Molecules

Problem Statement
Sub-Threshold Harmonics (BTH), generated via High-Harmonic Generation (HHG) from molecules excited by short, intense laser pulses, exhibit mechanisms fundamentally distinct from those of Above-Threshold Harmonics. While Above-Threshold Harmonics typically originate from a single active electron undergoing tunnel ionization, propagation, and recombination, BTH involves significant contributions from bound electrons propagating under the laser field along various electronic excited states. Consequently, BTH properties are highly target-dependent and rely on molecular energy levels and dipole transition matrix elements. This sensitivity renders BTH a valuable tool for High-Harmonic Spectroscopy and the generation of frequency combs in the Vacuum Ultraviolet (VUV) range. However, analyzing these contributions is complex. Existing semi-classical models, while effective for very strong fields ($I_0 \approx 10^{14}$ W/cm$^2$), fail to accurately describe BTH at moderate intensities and obscure the specific roles of discrete electronic transitions and their transition directions in determining the polarization properties of the emitted harmonics.

Methodology
To isolate and analyze the contribution of bound-bound electronic transitions to BTH, the authors employ multi-level model systems where specific transitions, dipole orientations, and resonances can be precisely controlled. The study proceeds in two stages:

1. Analysis of the Two-Level System (TLS): The authors model a single TLS with a field-free Hamiltonian and a transition dipole moment oriented along the x-axis. They solve the time-dependent Schrödinger equation (TDSE) numerically using an explicit Runge-Kutta method of order 5.(4) to obtain time-dependent populations and the resulting dipole acceleration. The harmonic spectrum is derived via Fourier transformation. The study focuses on moderate field strengths and long wavelengths to avoid continuum transitions and ensure that the dynamics are dominated by bound-bound interactions.
2. Construction of a Four-Level Model: Building on the TLS results, the authors construct a four-level model consisting of two decoupled TLSs aligned along orthogonal directions (x and y). These subsystems possess different transition frequencies ($\omega_{21x} = 0.1$ and $\omega_{30y} = 0.126$) but identical dipole strengths. The system is excited by a linearly polarized laser field with a variable polarization angle $\alpha$. The authors solve the TDSE for this 2D system numerically to investigate how the phase properties of individual TLSs manifest in the collective harmonic response of a molecule with two different axes.

Additionally, the authors provide an analytical derivation (in the Supplementary Material) based on the adiabatic approximation and Floquet theory. This derivation yields a continued fraction expression for the Fourier-transformed dipole moment, enabling the extraction of harmonic phases under periodic excitation.

Main Results

• Phase Alternation in TLS: In a single TLS, the phase of the emitted harmonics shows a pronounced dependence on the driving frequency relative to the transition energy. For harmonics with photon energies below the transition energy between the dominant field-dressed states, the phase alternates by $\pi$ between consecutive odd harmonic orders. In contrast, the phase remains constant for harmonics above the transition energy. This behavior is attributed to the sign change in the prefactor of the recursive relationship governing harmonic generation, specifically the term $(\omega_{10}^2 - n^2\omega_L^2)$.
• Polarization Mirroring in Orthogonal Systems: In the four-level model with orthogonal transition dipoles, the phase difference between the two subsystems leads to distinct polarization patterns. For harmonic orders below the resonance of the lower-frequency subsystem (orders 1–7), the phases of both subsystems are identical, resulting in harmonics that follow the polarization of the driving field. However, for higher harmonic orders (starting at order 9), the subsystem with the higher transition frequency operates in the "constant phase" regime, while the lower-frequency subsystem remains in the "alternating phase" regime. This generates a phase difference of $\pi$ between the x- and y-components of the emitted radiation.
• Counter-Rotation: This $\pi$ phase shift causes the polarization vector of higher harmonics to rotate in the opposite direction to the polarization of the driving field as the driving angle $\alpha$ increases. This "counter-rotation" effectively mirrors the polarization vector across the x-axis.
• Distinction from Interference: The authors note that while the resulting amplitude minima and phase jumps resemble interference effects from two centers, as observed in diatomic molecules (e.g., $H_2^+$), the microscopic origin here is different. The observed signatures arise from the interplay of two decoupled TLSs with different transition frequencies and orthogonal dipoles, not from spatial interference between emission centers.

Significance and Claims
The work claims that these findings demonstrate how laser-frequency-dependent phase jumps in simple multi-level systems translate directly into distinct polarization patterns in the harmonic response of aligned molecules with orthogonal transition dipoles. The authors suggest that aligned systems with orthogonal transition dipoles (such as certain organic molecules like acenes or perylene-bismides) could exhibit analogous phase and polarization features in the sub-threshold regime.

The significance of this work lies in providing a simplified framework for investigating the mechanisms responsible for the complex polarizations of BTH. The authors argue that an improved understanding of these anisotropic phase features could advance existing HHG-based orbital tomography in aligned molecular targets. Specifically, focusing on BTH could enable more detailed reconstructions that are sensitive to excited-state dynamics. The work acknowledges that in real physical systems, continuum contributions might mask these features; however, the model offers a clear theoretical prediction for systems where bound-bound transitions dominate.