Ultrafast nonadiabatic dynamics of tetraphenylsubstituted nitrogen-based heterocycles

This study employs mixed quantum-classical trajectory simulations to elucidate the distinct ultrafast nonadiabatic deactivation pathways of tetraphenylpyrazine (TPP) and tetraphenylpyrrole (TePP), revealing how differences in intramolecular flexibility govern their contrasting solid-state luminescence enhancement versus dual-state emission behaviors.

The Gist

Imagine you have two identical-looking dancers, TePP and TPP. They are both wearing four heavy, flowing capes (the phenyl rings) attached to a central body. You shine a bright light on them (photoexcitation), and they start to dance. The question scientists wanted to answer is: Why does one dancer glow brightly no matter where they are, while the other only glows when they are crowded in a room full of people?

Here is the story of their dance, broken down simply.

The Two Dancers

1. TePP (The "Dual-State" Dancer):

   • The Vibe: This dancer is naturally stiff and balanced. Whether they are dancing alone on a stage (in a gas/solution) or packed tight in a crowded club (solid state), they keep their rhythm and glow brightly.
   • The Secret: Their central body (a pyrrole ring) is small and tight. This forces their four capes to stay close together, creating a "steric hug." This hug prevents the dancer from twisting into awkward, energy-sapping poses. They stay in a "glowing" pose whether they are alone or crowded.

2. TPP (The "Solid-State" Dancer):

   • The Vibe: This dancer is flexible but tricky. When dancing alone, they immediately twist and contort their body into a weird shape that stops them from glowing. They look dull. But, if you pack them into a crowded room (solid state), the other dancers bump into them, physically stopping them from twisting. Suddenly, they are forced to stay in a glowing pose and shine brightly!
   • The Secret: Their central body (a pyrazine ring) is more open. When alone, they can easily twist their core, which acts like a "leak" that drains their energy before they can glow.

The Experiment: A High-Speed Camera

The scientists didn't just watch the dancers; they used a super-fast, high-tech camera (simulated by a computer) to film the dance in slow motion (femtoseconds, which is a quadrillionth of a second). They wanted to see exactly how the dancers moved right after the light hit them.

They used a method called "Surface Hopping." Imagine the dancers are running on a trampoline. Sometimes, the trampoline has a hole (a "conical intersection"). If the dancer hits the hole, they fall through to the ground and stop glowing (non-radiative decay).

• What happened to TePP?
  The dancer jumped, spun around, and moved their capes in a loose, wavy motion. They never found the hole. They kept bouncing on the trampoline, staying in the "glowing" zone. The computer showed that even when alone, their body was too stiff to twist into the "hole."

• What happened to TPP?
  The dancer jumped, but because their core was flexible, they immediately started twisting their central body. This twisting motion acted like a key, unlocking a door to the "hole." They fell through the trampoline almost instantly, losing their glow. The computer showed that this "twisting" happened so fast that even without a crowd to stop them, they were already dim.

The "Flash Photography" (Observables)

To prove this, the scientists simulated two types of "flash photography":

1. Time-Resolved Fluorescence (The Glow):

   • TePP: The glow stayed bright and steady for the whole video, just getting slightly redder (like a sunset) as the dancer relaxed.
   • TPP: The glow faded away quickly and turned red. This confirmed that the dancer was losing energy by twisting, not by hitting the ground.

2. Ultrafast Electron Diffraction (The Shape):

   • This is like taking a snapshot of the dancer's skeleton.
   • TePP: The snapshot showed the whole dancer moving together, like a wave passing through their capes.
   • TPP: The snapshot showed the center of the dancer twisting violently, while the capes stayed relatively still. This proved the "leak" was in the core, not the edges.

The Big Takeaway

For a long time, scientists thought TPP was dim in solution because the liquid solvent was "messing up" the dance.

This paper proves that's not true.

• TePP is a natural superstar. It glows because its own body is built to resist twisting. It doesn't need a crowd to shine.
• TPP is a "crowd-dependent" star. It is dim in solution not because of the liquid, but because it wants to twist itself into a dark shape. It only glows in the solid state because the crowd physically prevents it from twisting.

The Analogy:
Think of TePP as a rigid, well-built car that runs smoothly on both a bumpy dirt road (solution) and a smooth highway (solid).
Think of TPP as a car with a loose steering wheel. On the dirt road, the wheel wobbles so much the car crashes (no light). But on the highway, the guardrails (the other molecules) hold the wheel straight, so the car drives perfectly and shines.

The scientists used this study to show that we don't always need to look at the environment to understand why a molecule is bright or dark; sometimes, the answer is written in the molecule's own DNA (its internal structure).

Technical Summary

1. Problem Statement

The study addresses the fundamental challenge of understanding why structurally similar organic chromophores exhibit drastically different photophysical behaviors in different phases (solution vs. solid state). Specifically, it compares:

• Tetraphenylpyrazine (TPP): A prototypical Solid-State Luminescence Enhancement (SLE) emitter. It shows weak emission in solution (quantum yield $\Phi \approx 0.53\%$) but strong emission in the solid state ($\Phi \approx 8.3\%$).
• 2,3,4,5-Tetraphenyl-1H-pyrrole (TePP): A Dual-State Emission (DSE) emitter. It maintains high quantum yields in both solution ($\Phi \approx 65.6\%$) and the solid state ($\Phi \approx 74.3\%$).

The core scientific question is whether the quenching observed in TPP in solution is driven by intramolecular mechanisms (intrinsic molecular flexibility and relaxation pathways) or extrinsic environmental factors (solvent interactions or aggregation effects). The authors aim to disentangle these factors by analyzing the gas-phase ultrafast dynamics of both molecules, free from environmental constraints.

2. Methodology

The authors employed a mixed quantum-classical trajectory approach to simulate the nonadiabatic dynamics of single molecules in the gas phase.

• Electronic Structure:
  • Method: Time-Dependent Density Functional Theory (TD-DFT) using the TD-B3LYP-D3/def2-SVP functional and basis set.
  • Validation: The method was benchmarked against Equation of Motion Coupled Cluster (EOM-CCSD) and experimental UV-Vis spectra in previous work (Ref. 41). For TePP, high-level EOM-CCSD and DF-LCC2-LR calculations were used to validate oscillator strengths.
  • States Included: 12 singlet excited states were included in the dynamics to capture the dense manifold of states.
• Dynamics Simulation:
  • Algorithm: Tully's Fewest-Switches Surface Hopping (FSSH) implemented in the SHARC code.
  • Initial Conditions: Trajectories were initiated from a Wigner distribution sampled from the ground-state vibrational modes.
  • Ensemble Size: 206 retained trajectories for TPP and 148 for TePP (after energy conservation filtering).
  • Decoherence: Energy-based decoherence corrections (Granucci and Persico) were applied to prevent artificial coherence.
• Observables:
  • Time-Resolved Fluorescence (TR-FL): Calculated from instantaneous emission intensity ($A(t) = \omega(t)^3 f_{osc}(t)$) for trajectories on the $S_1$ state.
  • Gas-Phase Ultrafast Electron Diffraction (GUED): Simulated via the Independent Atom Model to generate time-resolved difference Pair Distribution Functions ($\Delta$PDF) to track structural changes.

3. Key Contributions

• Mechanistic Differentiation: The study provides a rigorous, state-resolved comparison of the intrinsic non-radiative decay pathways for SLE vs. DSE molecules, demonstrating that their contrasting behaviors are rooted in intramolecular electronic structure and nuclear dynamics rather than solely environmental effects.
• Beyond Population Analysis: Instead of relying solely on electronic state populations (which are theoretical constructs), the authors link dynamics to simulated experimental observables (TR-FL and GUED), offering a more physically grounded interpretation of the data.
• Role of Molecular Rigidity: The work elucidates how specific structural motifs (pyrazine vs. pyrrole cores) dictate the accessibility of conical intersections and dark states, validating the Restricted Access to Conical Intersections (RACI) model in a gas-phase context.

4. Key Results

A. Electronic Structure at the Franck-Condon Region

• TePP: Excited states are dominated by Intramolecular Charge Transfer (CT) character. The donor is the pyrrole $\pi$-orbital, and acceptors are delocalized over phenyl rings. The states form a dense, strongly coupled manifold.
• TPP: Excited states are primarily Locally Excited (LE) $\pi\pi^*$ on the pyrazine core, with higher-lying $n\pi^*$ states. The states are more symmetry-stabilized and less coupled initially compared to TePP.

B. Population Dynamics (0–200 fs)

• TePP: Exhibits ultrafast, diffusive mixing within the CT manifold. Population spreads rapidly across multiple states (half-life of initial character $\approx 11$ fs) but remains in emissive configurations. No significant population reaches the ground state ($S_0$) within 200 fs.
• TPP: Shows slower, sequential relaxation. Population remains in the LE $\pi\pi^*$ manifold longer (half-life $\approx 33$ fs) before progressively transferring to dark $n\pi^*$ states.
• Crucial Finding: Neither molecule accesses the ground state ($S_0$) within the 200 fs simulation window. This implies that the quenching of TPP in solution is not due to ultrafast internal conversion to $S_0$, but rather trapping in dark excited states ($n\pi^*$).

C. Structural Dynamics (Dihedrals and Normal Modes)

• TePP: Dynamics are driven by low-frequency, peripheral motions (phenyl wagging, N-H distortions). These motions are loosely correlated and allow efficient redistribution without driving the system into quenched geometries.
• TPP: Dynamics are driven by mid-frequency, core-localized distortions (pyrazine ring deformations and symmetry-breaking). These specific collective motions facilitate access to the dark $n\pi^*$ states.

D. Simulated Observables

• Time-Resolved Fluorescence (TR-FL):
  • TePP: Maintains high emission intensity with a modest red-shift ($\sim 0.5$ eV) over 200 fs, consistent with DSE behavior.
  • TPP: Shows a rapid loss of intensity and a large red-shift ($\sim 1$ eV), indicating population trapping in non-emissive states. This explains the low solution quantum yield.
• Ultrafast Electron Diffraction (GUED):
  • TePP: $\Delta$PDF signals extend to large interatomic distances ($r > 8-9$ Å), indicating correlated long-range motion transmitting structural changes from the core to distant phenyls.
  • TPP: $\Delta$PDF signals are localized at intermediate distances ($r \approx 4-6$ Å), reflecting localized core distortions with minimal long-range correlation.

5. Significance and Conclusions

The paper concludes that the contrasting photoluminescence of TPP and TePP is an intrinsic property of their molecular frameworks:

1. TePP (DSE): Its pyrrole core and specific phenyl arrangement create steric constraints that stabilize emissive CT states and hinder access to dark states, regardless of the environment (gas, solution, or solid).
2. TPP (SLE): Its pyrazine core possesses an intrinsic relaxation pathway involving core distortions that funnel population into dark $n\pi^*$ states. In solution, this pathway is active, causing quenching. In the solid state, molecular packing (steric hindrance) restricts these specific core distortions, effectively blocking the non-radiative funnel and restoring emission (SLE).

Implication: The study demonstrates that gas-phase dynamics can accurately predict solution-phase behavior for these systems, as the dominant quenching mechanism is intramolecular. It reinforces the RACI model, showing that solid-state enhancement arises not from creating new emissive pathways, but from blocking the intrinsic non-radiative decay coordinates available in the isolated molecule. This provides a robust theoretical framework for designing future organic emitters for OLEDs and bioimaging.