Influence of Vibronic Coupling on Ultrafast Singlet Fission in a Linear Terrylenediimide Dimer

Exploring the Role of Excited States' Degeneracy on Vibronic Coupling with Atomic-Scale Optics

Interactions between molecular electronic and vibrational states manifest themselves in a variety of forms and have a strong impact on molecular physics and chemistry. For example, the efficiency of energy transfer between organic molecules, ubiquitous in biological systems and in organic optoelectronics, is strongly influenced by vibronic coupling. Using an approach based on scanning tunneling microscope-induced luminescence (STML), we reveal vibronic interactions in optical spectra of a series of single phthalocyanine derivative molecules featuring degenerate or near-degenerate excited states. Based on detailed theoretical simulations, we disentangle spectroscopic signatures belonging to Franck-Condon and Herzberg-Teller vibronic progressions in tip-position-resolved STML spectra, and we directly map out the vibronic coupling between the close-lying excited states of the molecules.

Distinct vibrational motions promote disparate excited-state decay pathways in cofacial perylenediimide dimers

A complex interplay of structural, electronic, and vibrational degrees of freedom underpins the fate of molecular excited states. Organic assemblies exhibit a myriad of excited-state decay processes, such as symmetry-breaking charge separation (SB-CS), excimer (EX) formation, singlet fission, and energy transfer. Recent studies of cofacial and slip-stacked perylene-3,4:9,10-bis(dicarboximide) (PDI) multimers demonstrate that slight variations in core substituents and H- or J-type aggregation can determine whether the system follows an SB-CS pathway or an EX one. However, questions regarding the relative importance of structural properties and molecular vibrations in driving the excited-state dynamics remain. Here, we use a combination of two-dimensional electronic spectroscopy, femtosecond stimulated Raman spectroscopy, and quantum chemistry computations to compare the photophysics of two PDI dimers. The dimer with 1,7-bis(pyrrolidin-1'-yl) substituents (5PDI2) undergoes ultrafast SB-CS from a photoexcited mixed state, while the dimer with bis-1,7-(3',5'-di-t-butylphenoxy) substituents (PPDI2) rapidly forms an EX state. Examination of their quantum beating features reveals that SB-CS in 5PDI2 is driven by the collective vibronic coupling of two or more excited-state vibrations. In contrast, we observe signatures of low-frequency vibrational coherence transfer during EX formation by PPDI2, which aligns with several previous studies. We conclude that key electronic and structural differences between 5PDI2 and PPDI2 determine their markedly different photophysics.

Robust vibrational coherence protected by a core-shell structure in silver nanoclusters

Vibrational coherence has attracted considerable research interests because of its potential functions in light harvesting systems. Although positive signs of vibrational coherence in metal nanoclusters have been observed, the underlying mechanism remains to be verified. Here, we demonstrate that robust vibrational coherence with a lifetime of 1 ps can be clearly identified in Ag44(SR)30 core-shell nanoclusters, in which an icosahedral Ag12 core is well protected by a dodecahedral Ag20 cage. Ultrafast spectroscopy reveals that two vibrational modes at around 2.4 THz and 1.6 THz, corresponding to the breathing mode and quadrupolar-like mode of the icosahedral Ag12 core, respectively, are responsible for the generation of vibrational coherence. In addition, the vibrational coherence of Ag44 has an additional high frequency mode (2.4 THz) when compared with that of Ag29, in which there is only one low frequency vibration mode (1.6 THz), and the relatively faster dephasing in two-layer Ag29 relative to that in Ag44 further supports the fact that the robust vibrational coherence in Ag44 is ascribed to its unique matryoshka-like core-shell structure. Our findings not only present unambiguous experimental evidence for a multi-layer core-shell structure protected vibrational coherence under ambient conditions but also offers a practical strategy for the design of highly efficient quantum optoelectronic devices.

Two-Dimensional Electronic Spectroscopy Resolves Relative Excited-State Displacements

Knowledge of relative displacements between potential energy surfaces (PES) is critical in spectroscopy and photochemistry. Information on displacements is encoded in vibrational coherences. Here we apply ultrafast two-dimensional electronic spectroscopy in a pump-probe half-broadband (HB2DES) geometry to probe the ground- and excited-state potential landscapes of cresyl violet. 2D coherence maps reveal that while the coherence amplitude of the dominant 585 cm-1 Raman-active mode is mainly localized in the ground-state bleach and stimulated emission regions, a 338 cm-1 mode is enhanced in excited-state absorption. Modeling these data with a three-level displaced harmonic oscillator model using the hierarchical equation of motion-phase matching approach (HEOM-PMA) shows that the S1 ← S0 PES displacement is greater along the 585 cm-1 coordinate than the 338 cm-1 coordinate, while Sn ← S1 displacements are similar along both coordinates. HB2DES is thus a powerful tool for exploiting nuclear wavepackets to extract quantitative multidimensional, vibrational coordinate information across multiple PESs.

Mastering Femtosecond Stimulated Raman Spectroscopy: A Practical Guide

Femtosecond stimulated Raman spectroscopy (FSRS) is a powerful nonlinear spectroscopic technique that probes changes in molecular and material structure with high temporal and spectral resolution. With proper spectral interpretation, this is equivalent to mapping out reactive pathways on highly anharmonic excited-state potential energy surfaces with femtosecond to picosecond time resolution. FSRS has been used to examine structural dynamics in a wide range of samples, including photoactive proteins, photovoltaic materials, plasmonic nanostructures, polymers, and a range of others, with experiments performed in multiple groups around the world. As the FSRS technique grows in popularity and is increasingly implemented in user facilities, there is a need for a widespread understanding of the methodology and best practices. In this review, we present a practical guide to FSRS, including discussions of instrumentation, as well as data acquisition and analysis. First, we describe common methods of generating the three pulses required for FSRS: the probe, Raman pump, and actinic pump, including a discussion of the parameters to consider when selecting a beam generation method. We then outline approaches for effective and efficient FSRS data acquisition. We discuss common data analysis techniques for FSRS, as well as more advanced analyses aimed at extracting small signals on a large background. We conclude with a discussion of some of the new directions for FSRS research, including spectromicroscopy. Overall, this review provides researchers with a practical handbook for FSRS as a technique with the aim of encouraging many scientists and engineers to use it in their research.

Theory predicts UV/vis-to-IR photonic down conversion mediated by excited state vibrational polaritons

This work proposes a photophysical phenomenon whereby ultraviolet/visible (UV/vis) excitation of a molecule involving a Franck-Condon (FC) active vibration yields infrared (IR) emission by strong coupling to an optical cavity. The resulting UV/vis-to-IR photonic down conversion process is mediated by vibrational polaritons in the electronic excited state potential. It is shown that the formation of excited state vibrational polaritons (ESVP) via UV/vis excitation only involve vibrational modes with both a non-zero FC activity and IR activity in the excited state. Density functional theory calculations are used to identify 1-Pyreneacetic acid as a molecule with this property and the dynamics of ESVP are modeled. Overall, this work introduces an avenue of polariton chemistry where excited state dynamics are influenced by the formation of vibrational polaritons. Along with this, the UV/vis-to-IR photonic down conversion is potentially useful in both sensing excited state vibrations and quantum transduction schemes.

Unveiling the double triplet nature of the 2Ag state in conjugated stilbenoid compounds to achieve efficient singlet fission

In this investigation, the excited-state evolution in a series of all-trans stilbenoid compounds, displaying a low-lying dark singlet state of 2Ag-like symmetry nearly degenerate with the bright 1Bu state, was unveiled by employing advanced ultrafast spectroscopies while probing the effect of solvent polarizability. Together with the dual emission, femtosecond transient absorption and broadband fluorescence up-conversion disclosed the double nature of the 2Ag-like state showing both singlet features, a lifetime typical of a singlet and the ability to emit, and a triplet character, exhibiting a triplet-like absorption spectrum. The ultrafast formation (in hundreds of femtoseconds) from the non-relaxed upper singlet state led to the identification of 2Ag as the correlated triplet pair of singlet fission. The spectral difference obtained by comparison of transient absorption peaks of the 2Ag (1TT) and the triplet states was found to be in remarkable agreement with the observed triplet yield and the 1(TT) separation rate constant. Indeed, this spectral shift provided an experimental method to gain qualitative insight into the ease of separation of the 1(TT) and the relative SF efficiency. The highly conjugated polyene-like structures enable the ultrafast formation of the double triplet, but then the large binding energy prevents the triplet separation and thus the effective completion of singlet fission. Even though thermodynamically feasible for all the investigated stilbenoids according to TD-DFT calculations, singlet fission resulted to occur efficiently in the case of 1-(pyridyl-4-ylethenyl)-4-(p-nitrostyryl)benzene and nitro-styrylfuran with the triplet yield reaching 120% and 140%, respectively, triggered by their greatly enhanced intramolecular charge transfer character relative to the other compounds in the series.

Single-Molecule Ultrafast Fluorescence-Detected Pump-Probe Microscopy

We introduce fluorescence-detected pump-probe microscopy by combining a wavelength-tunable ultrafast laser with a confocal scanning fluorescence microscope, enabling access to the femtosecond time scale on the micrometer spatial scale. In addition, we obtain spectral information from Fourier transformation over excitation pulse-pair time delays. We demonstrate this new approach on a model system of a terrylene bisimide (TBI) dye embedded in a PMMA matrix and acquire the linear excitation spectrum as well as time-dependent pump-probe spectra simultaneously. We then push the technique toward single TBI molecules and analyze the statistical distribution of their excitation spectra. Furthermore, we demonstrate the ultrafast transient evolution of several individual molecules, highlighting their different behavior in contrast to the ensemble due to their individual local environment. By correlating the linear and nonlinear spectra, we assess the effect of the molecular environment on the excited-state energy.

Coherent Vibronic Wavepackets Show Structure-Directed Charge Flow in Host-Guest Donor-Acceptor Complexes

Designing and controlling charge transfer (CT) pathways in organic semiconductors are important for solar energy applications. To be useful, a photogenerated, Coulombically bound CT exciton must further separate into free charge carriers; direct observations of the detailed CT relaxation pathways, however, are lacking. Here, photoinduced CT and relaxation dynamics in three host-guest complexes, where a perylene (Per) electron donor guest is incorporated into two symmetric and one asymmetric extended viologen cyclophane acceptor hosts, are presented. The central ring in the extended viologen is either p-phenylene (ExV²⁺) or electron-rich 2,5-dimethoxy-p-phenylene (ExMeOV²⁺), resulting in two symmetric cyclophanes with unsubstituted or methoxy-substituted central rings, ExBox⁴⁺ and ExMeOBox⁴⁺, respectively, and an asymmetric cyclophane with one of the central viologen rings being methoxylated ExMeOVBox⁴⁺. Upon photoexcitation, the asymmetric host-guest ExMeOVBox⁴⁺ ⊃ Per complex exhibits directional CT toward the energetically unfavorable methoxylated side due to structural restrictions that facilitate strong interactions between the Per donor and the ExMeOV²⁺ side. The CT state relaxation pathways are probed using ultrafast optical spectroscopy by focusing on coherent vibronic wavepackets, which are used to identify CT relaxations along charge localization and vibronic decoherence coordinates. Specific low- and high-frequency nuclear motions are direct indicators of a delocalized CT state and the degree of CT character. Our results show that the CT pathway can be controlled by subtle chemical modifications of the acceptor host in addition to illustrating how coherent vibronic wavepackets can be used to probe the nature and time evolution of the CT states.

Singlet fission as a polarized spin generator for dynamic nuclear polarization

Singlet fission (SF), converting a singlet excited state into a spin-correlated triplet-pair state, is an effective way to generate a spin quintet state in organic materials. Although its application to photovoltaics as an exciton multiplier has been extensively studied, the use of its unique spin degree of freedom has been largely unexplored. Here, we demonstrate that the spin polarization of the quintet multiexcitons generated by SF improves the sensitivity of magnetic resonance of water molecules through dynamic nuclear polarization (DNP). We form supramolecular assemblies of a few pentacene chromophores and use SF-born quintet spins to achieve DNP of water-glycerol, the most basic biological matrix, as evidenced by the dependence of nuclear polarization enhancement on magnetic field and microwave power. Our demonstration opens a use of SF as a polarized spin generator in bio-quantum technology.

Connecting the dots for fundamental understanding of structure-photophysics-property relationships of COFs, MOFs, and perovskites using a Multiparticle Holstein Formalism

Photoactive organic and hybrid organic-inorganic materials such as conjugated polymers, covalent organic frameworks (COFs), metal-organic frameworks (MOFs), and layered perovskites, display intriguing photophysical signatures upon interaction with light. Elucidating structure-photophysics-property relationships across a broad range of functional materials is nontrivial and requires our fundamental understanding of the intricate interplay among excitons (electron-hole pair), polarons (charges), bipolarons, phonons (vibrations), inter-layer stacking interactions, and different forms of structural and conformational defects. In parallel with electronic structure modeling and data-driven science that are actively pursued to successfully accelerate materials discovery, an accurate, computationally inexpensive, and physically-motivated theoretical model, which consistently makes quantitative connections with conceptually complicated experimental observations, is equally important. Within this context, the first part of this perspective highlights a unified theoretical framework in which the electronic coupling as well as the local coupling between the electronic and nuclear degrees of freedom can be efficiently described for a broad range of quasiparticles with similarly structured Holstein-style vibronic Hamiltonians. The second part of this perspective discusses excitonic and polaronic photophysical signatures in polymers, COFs, MOFs, and perovskites, and attempts to bridge the gap between different research fields using a common theoretical construct - the Multiparticle Holstein Formalism. We envision that the synergistic integration of state-of-the-art computational approaches with the Multiparticle Holstein Formalism will help identify and establish new, transformative design strategies that will guide the synthesis and characterization of next-generation energy materials optimized for a broad range of optoelectronic, spintronic, and photonic applications.

Singlet fission dynamics modulated by molecular configuration in covalently linked pyrene dimers, Anti- and Syn-1,2-di(pyrenyl)benzene

Covalently linked dimers (CLDs) and their structural isomers have attracted much attention as potential materials for improving power conversion efficiencies through singlet fission (SF). Here, we designed and synthesized two covalently ortho-linked pyrene (Py) dimers, anti- and syn-1,2-di(pyrenyl)benzene (Anti-DPyB and Syn-DPyB, respectively), and investigated the effect of molecular configuration on SF dynamics using steady-state and time-resolved spectroscopies. Both Anti-DPyB and Syn-DPyB, which have different Py-stacking configurations, form excimers, which then relax to the correlated triplet pair ((T₁T₁)) state, indicating the occurrence of SF. Unlike previous studies where the excimer formation inhibited an SF process, the (T₁T₁)'s of Anti-DPyB and Syn-DPyB are formed through the excimer state. The dissociation of (T₁T₁)'s to 2T₁ in Anti-DPyB is more favorable than in Syn-DPyB. Our results showcase that the molecular configuration of a CLD plays an important role in SF dynamics.

Inducing Singlet Fission in Perylene Thin Films by Molecular Contortion

Singlet fission occurs only in a limited number of molecules, and expanding the molecular toolbox is necessary for progress. Here, we apply the molecular contortion strategy to tune singlet and triplet energies and observe changes in the excited-state dynamics that are consistent with singlet fission playing a role in thin films of contorted perylene. Perylene is a prototypical molecular chromophore, which does not undergo singlet fission in its planar form from its S₁ state. The tuning of the energetics that control singlet fission through molecular contortion can be applied to a large repertoire of established molecular chromophores.

Mapping Simulated Two-Dimensional Spectra to Molecular Models Using Machine Learning

Two-dimensional (2D) spectroscopy encodes molecular properties and dynamics into expansive spectral data sets. Translating these data into meaningful chemical insights is challenging because of the many ways chemical properties can influence the spectra. To address the task of extracting chemical information from 2D spectroscopy, we study the capacity of simple feedforward neural networks (NNs) to map simulated 2D electronic spectra to underlying physical Hamiltonians. We examined hundreds of simulated 2D spectra corresponding to monomers and dimers with varied Franck-Condon active vibrations and monomer-monomer electronic couplings. We find the NNs are able to correctly characterize most Hamiltonian parameters in this study with an accuracy above 90%. Our results demonstrate that NNs can aid in interpreting 2D spectra, leading from spectroscopic features to underlying effective Hamiltonians.

Vibrational response functions for multidimensional electronic spectroscopy in the adiabatic regime: A coherent-state approach

Multi-dimensional spectroscopy represents a particularly insightful tool for investigating the interplay of nuclear and electronic dynamics, which plays an important role in a number of photophysical processes and photochemical reactions. Here, we present a coherent state representation of the vibronic dynamics and of the resulting response functions for the widely used linearly displaced harmonic oscillator model. Analytical expressions are initially derived for the case of third-order response functions in an N-level system, with ground state initialization of the oscillator (zero-temperature limit). The results are then generalized to the case of Mth order response functions, with arbitrary M. The formal derivation is translated into a simple recipe, whereby the explicit analytical expressions of the response functions can be derived directly from the Feynman diagrams. We further generalize to the whole set of initial coherent states, which form an overcomplete basis. This allows one, in principle, to derive the dependence of the response functions on arbitrary initial states of the vibrational modes and is here applied to the case of thermal states. Finally, a non-Hermitian Hamiltonian approach is used to include in the above expressions the effect of vibrational relaxation.

π-Stacking-Dependent Vibronic Couplings Drive Excited-State Dynamics in Perylenediimide Assemblies

Vibronic coupling, the interplay of electronic and nuclear vibrational motion, is considered a critical mechanism in photoinduced reactions such as energy transfer, charge transfer, and singlet fission. However, our understanding of how particular vibronic couplings impact excited-state dynamics is lacking due to the limited number of experimental studies of model molecular systems. Herein, we use two-dimensional electronic spectroscopy (2DES) to launch and interrogate a range of vibronic coherences in two distinct types of perylenediimide slip stacks─along the short and long molecular axes, which form either an excimer or a mixed state between the Frenkel exciton (FE) and charge transfer states. We explore the functionality of these vibronic coherences using quantum beatmaps, which display the Fourier amplitude signal oscillations as a function of pump and probe frequencies, along with knowledge of the characteristic signatures of the FE, ionic, and excimer species. We find that a low-frequency vibrational mode of the short-axis slip stack appears concomitantly with the formation of the excimer state, survives 2-fold longer than in the FE state in the reference monomer, and shows a phase shift compared to other modes. For the long-axis slip stacks, a pair of low-frequency modes coupled to a high-frequency coordinate of the FE state were found to play a critical role in mixed-state generation. Our findings thus experimentally reveal the complex and varying roles of vibronic couplings in tightly packed multimers undergoing a range of photoinduced processes.

Accelerating symmetry-breaking charge separation in a perylenediimide trimer through a vibronically coherent dimer intermediate

Understanding the photophysics and photochemistry of molecular π-stacked chromophores is important for utilizing them as functional photonic materials. However, these investigations have been mostly limited to covalent molecular dimers, which can only approximate the electronic and vibronic interactions present in the higher oligomers typical of functional organic materials. Here we show that a comparison of the excited-state dynamics of a covalent slip-stacked perylenediimide dimer (2) and trimer (3) provides fundamental insights into electronic state mixing and symmetry-breaking charge separation (SB-CS) beyond the dimer limit. We find that coherent vibronic coupling to high-frequency modes facilitates ultrafast state mixing between the Frenkel exciton (FE) and charge-transfer (CT) states. Subsequently, solvent fluctuations and interchromophore low-frequency vibrations promote CT character in the coherent FE/CT mixed state. The coherent FE/CT mixed state persists in 2, but, in 3, low-frequency vibronic coupling collapses the coherence, resulting in ultrafast SB-CS between the distal perylenediimide units.

Exploring a new class of singlet fission fluorene derivatives with high-energy triplets

In this study, we report strong experimental evidence for singlet fission (SF) in a new class of fluorene-based molecules, exhibiting two-branched donor-acceptor structures. The time-resolved spectroscopic results disclose ultrafast formation of a double triplet state (occurring in few picoseconds) and efficient triplet exciton separation (up to 145% triplet yield). The solvent polarity effect and the role of intramolecular charge transfer (ICT) on the SF mechanism have been thoroughly investigated with several advanced spectroscopies. We found that a stronger push-pull character favors SF, as long as the ICT does not act as a trap by opening a competitive pathway. Within the context of other widely-known SF chromophores, the unconventional property of generating high-energy triplet excitons (ca. 2 eV) via SF makes these materials outstanding candidates as photosensitizers for photovoltaic devices.

Computational predictions of adaptive aromaticity for the design of singlet fission materials

The concept of adaptive aromaticity has been demonstrated as an alternative strategy for the design of singlet fission materials.

Towards multistate multimode landscapes in singlet fission of pentacene: the dual role of charge-transfer states

Singlet fission duplicates triplet excitons for improving light harvesting efficiency. The presence of the interaction between electronic and nuclear degrees of freedom complicates the interpretation of correlated triplet pairs. We report a quantum chemistry study on the significance and subtleties of multistate and multimode pathways in forming triplet pair states of the pentacene dimer through a six-state vibronic-coupling Hamiltonian derived from many-electron adiabatic wavefunctions of an ab initio density matrix renormalization group. The resulting spin values of the singlet manifolds on each pentacene center are computed, and the varying spin nature can be distinguished clearly with respect to dimer stacking and vibronic progression. Our monomer spin assignments reveal the coexistence of both lower-lying weak and higher-lying strong charge transfer states which interact vibronically with the triplet pair state, providing important implications for its generation and separation occurring in vibronic regions. This work conveys the importance of the many-electron process requiring close low-lying singlet manifolds to determine the subtle fission details, and represents an important step for understanding vibronically resolved spin states and conversions underlying efficient singlet fission.

Vibronic coupling in energy transfer dynamics and two-dimensional electronic-vibrational spectra

We introduce a heterodimer model in which multiple mechanisms of vibronic coupling and their impact on energy transfer can be explicitly studied. We consider vibronic coupling that arises through either Franck-Condon activity in which each site in the heterodimer has a local electron-phonon coupling or Herzberg-Teller activity in which the transition dipole moment coupling the sites has an explicit vibrational mode-dependence. We have computed two-dimensional electronic-vibrational (2DEV) spectra for this model while varying the magnitude of these two effects and find that 2DEV spectra contain static and dynamic signatures of both types of vibronic coupling. Franck-Condon activity emerges through a change in the observed excitonic structure, while Herzberg-Teller activity is evident in the appearance of significant side-band transitions that mimic the lower-energy excitonic structure. A comparison of quantum beating patterns obtained from analysis of the simulated 2DEV spectra shows that this technique can report on the mechanism of energy transfer, elucidating a means of experimentally determining the role of specific vibronic coupling mechanisms in such processes.

Temperature Tuning of Coherent Mixing between States Driving Singlet Fission in a Spiro-Fused Terrylenediimide Dimer

The excited-state dynamics of a spiro-fused terrylene-3,4:11,12-bis(dicarboximide) (TDI) dimer (sTDI2) in toluene and 2-methyltetrahydrofuran (mTHF) were investigated as a function of temperature using femtosecond- and nanosecond-transient absorption spectroscopy, as well as two-dimensional electronic spectroscopy. The spiro conjugation and the corresponding geometry of this compound guarantee a short intermonomer distance along with a partial orbital overlap between the orthogonal TDI π-electron systems, providing electronic coupling between the TDIs. Photoexcitation of sTDI2 in toluene, a low dielectric solvent, at 295 K, results in the ultrafast formation of a state composed of a coherent mixture of singlet 1(S1S0), multiexciton 1(T1T1), and charge-transfer (CT) electronic characters. This mixed species decays to decorrelated triplet states on the nanosecond timescale, completing the process of intramolecular singlet fission (SF) in sTDI2. Upon decreasing the temperature from 295 to 200 K, the contribution of the 1(T1T1) state to the mixed species decreases concurrently with an increase in the CT state character. We attribute this behavior to the variation in the vibrational energy level alignment between the states comprising the mixture due to changes in the temperature and hence the local dielectric environment. In contrast, photoexcitation of sTDI2 in more polar mTHF at 295 K results in the formation of a mixed singlet and CT state before undergoing symmetry-breaking charge separation, owing to the increased stabilization of the CT state in the medium. However, in glassy mTHF at 85 K, photoexcited sTDI2 exhibits discernible multiexciton character, comparable to that observed in toluene at 200 K, which we rationalize by the similarity of the dielectric constants under these two sets of conditions. These observations of mixed states of varying diabatic contributions over the range of experimental conditions show that the temperature and the static dielectric constant can directly control the composition of the electronically mixed excited state of sTDI2 and thus the fate of the SF process.

Exploitation of Baird Aromaticity and Clar’s Rule for Tuning the Triplet Energies of Polycyclic Aromatic Hydrocarbons

Polycyclic aromatic hydrocarbons (PAH) are a prominent substance class with a variety of applications in molecular materials science. Their electronic properties crucially depend on the bond topology in ways that are often highly non-intuitive. Here, we study, using density functional theory, the triplet states of four biphenylene-derived PAHs finding dramatically different triplet excitation energies for closely related isomeric structures. These differences are rationalised using a qualitative description of Clar sextets and Baird quartets, quantified in terms of nucleus independent chemical shifts, and represented graphically through a recently developed method for visualising chemical shielding tensors (VIST). The results are further interpreted in terms of a 2D rigid rotor model of aromaticity and through an analysis of the natural transition orbitals involved in the triplet excited states showing good consistency between the different viewpoints. We believe that this work constitutes an important step in consolidating these varying viewpoints of electronically excited states.