Accurate Assignments of Excited-State Resonance Raman Spectra: A Benchmark Study Combining Experiment and Theory

The Best Models of Bodipy's Electronic Excited State: Comparing Predictions from Various DFT Functionals with Measurements from Femtosecond Stimulated Raman Spectroscopy

Density functional theory (DFT) and time-dependent DFT (TD-DFT) are pivotal approaches for modeling electronically excited states of molecules. However, choosing a DFT exchange-correlation functional (XCF) among the myriad of alternatives is an overwhelming task that can affect the interpretation of results and lead to erroneous conclusions. The performance of these XCFs to describe the excited-state properties is often addressed by comparing them with high-level wave function methods or experimentally available vertical excitation energies; however, this is a limited analysis that relies on evaluation of a single point in the excited-state potential energy surface (PES). Different strategies have been proposed but are limited by the difficulty of experimentally accessing the electronic excited-state properties. In this work, we have tested the performance of 12 different XCFs and TD-DFT to describe the excited-state potential energy surface of Bodipy (2,6-diethyl-1,3,5,7-tetramethyl-8-phenyldipyrromethene difluoroborate). We compare those results with resonance Raman spectra collected by using femtosecond stimulated Raman spectroscopy (FSRS). By simultaneously fitting the absorption spectrum, fluorescence spectrum, and all of the resonance Raman excitation profiles within the independent mode displaced harmonic oscillator (IMDHO) formalism, we can describe the PES at the Franck-Condon (FC) region and determine the solvent and intramolecular reorganization energy after relaxation. This allows a direct comparison of the TD-DFT output with experimental observables. Our analysis reveals that using vertical absorption energies might not be a good criterion to determine the best XCF for a given molecular system and that FSRS opens up a new way to benchmark the excited-state performance of XCFs of fluorescent dyes.

Spectro-Electrochemical Properties of A New Non-Enzymatic Modified Working Electrode Used for Histamine Assessment in the Diagnosis of Food Poisoning

We successfully prepared a non-enzymatic sensor based on a graphene-thiophene composite for histamine detection. The self-assembling properties of the thiophene onto Au support and the high electrical conductivity of graphene encouraged the choice of this type of composite. The composite was deposited via electrochemical polymerization onto the Au layer of a screen-printed microelectrode. The electropolymerization and electrochemical detection of histamine were both achieved by cyclic voltammetry. Two types of electrolytes were used for the electrochemical detection: (a) phosphate buffer solution (PBS), which showed low-intensity redox peaks for histamine; and (b) trichloroacetic acid (TCA) 0.01 M, which showed improved results over PBS and did not damage the microelectrode. For the concentration range of 100-200 mg/kg, the sensor shows a linear regression pattern for the oxidation peak fitted on the equation Ipa = 123.412 + 0.49933 ×x, with R2 = 0.94178. The lowest limit of detection was calculated to be 13.8 mg/kg and the limit of quantification was calculated at 46 mg/kg. These results are important since by monitoring the amount of histamine in a food product, early onset of spoilage can be easily detected, thus reducing foodborne poisoning and food waste (by recycling products that are still edible).

Solvent Tuning Excited State Structural Dynamics in a Novel Bianthryl

Symmetry breaking charge separation (SBCS) is central to photochemical energy conversion. The widely studied 9,9-bianthryl (9,9'BA) is the prototype, but the role of bianthryl structure is hardly investigated. Here we investigate excited state structural dynamics in a bianthryl of reduced symmetry, 1,9-bianthryl (1,9'BA), through ultrafast electronic and vibrational spectroscopy. Resonance selective Raman in polar solvents reveals a Franck-Condon state mode that disappears concomitant with the rise of ring breathing modes of radical species. Solvent-dependent dynamics show that CS is driven by solvent orientational motion, as in 9,9'BA. In nonpolar solvents the excited state undergoes multistep structural relaxation, including subpicosecond Franck-Condon state decay and biexponential diffusion-controlled structural evolution to a distorted slightly polar state. These data suggest two possible routes to SBCS; the established solvent driven pathway in rapidly relaxing polar solvents and, in slowly relaxing media, initial intramolecular reorganization to a polar structure which drives solvent orientational relaxation.

Photophysics of the red-form Kaede chromophore

The chromophore responsible for colour switching in the optical highlighting protein Kaede has unexpectedly complicated excited state dynamics, which are measured and analysed here. This will inform the development of new imaging proteins.

Excited-state resonance Raman spectroscopy probes the sequential two-photon excitation mechanism of a photochromic molecular switch

Some diarylethene molecular switches have a low quantum yield for cycloreversion when excited by a single photon, but react more efficiently following sequential two-photon excitation. The increase in reaction efficiency depends on both the relative time delay and the wavelength of the second photon. This paper examines the wavelength-dependent mechanism for sequential excitation using excited-state resonance Raman spectroscopy to probe the ultrafast (sub-30 fs) dynamics on the upper electronic state following secondary excitation. The approach uses femtosecond stimulated Raman scattering (FSRS) to measure the time-gated, excited-state resonance Raman spectrum in resonance with two different excited-state absorption bands. The relative intensities of the Raman bands reveal the initial dynamics in the higher-lying states, S_n, by providing information on the relative gradients of the potential energy surfaces that are accessed via secondary excitation. The excited-state resonance Raman spectra reveal specific modes that become enhanced depending on the Raman excitation wavelength, 750 or 400 nm. Many of the modes that become enhanced in the 750 nm FSRS spectrum are assigned as vibrational motions localized on the central cyclohexadiene ring. Many of the modes that become enhanced in the 400 nm FSRS spectrum are assigned as motions along the conjugated backbone and peripheral phenyl rings. These observations are consistent with earlier measurements that showed higher efficiency following secondary excitation into the lower excited-state absorption band and illustrate a powerful new way to probe the ultrafast dynamics of higher-lying excited states immediately following sequential two-photon excitation.

Ultrafast Dynamics of a Molecular Switch from Resonance Raman Spectroscopy: Comparing Visible and UV Excitation

Resonance Raman spectroscopy probes the ultrafast dynamics of a diarylethene (DAE) molecular switch following excitation into the first two optical absorption bands. Mode-specific resonance enhancements for Raman excitation at visible (750-560 nm) and near-UV (420-390 nm) wavelengths compared with the calculated and experimental off-resonance Raman spectrum at 785 nm reveal different Franck-Condon active vibrations for the two electronically excited states. The resonance enhancements at visible wavelengths are consistent with initial motion on the first excited-state that promotes the cycloreversion reaction, whereas the enhancements for excitation at near-UV wavelengths highlight motions involving conjugated backbone and phenyl ring stretching modes that are orthogonal to the reaction coordinate. The results support a mechanism involving rapid internal conversion from the higher-lying state followed by cycloreversion on the first excited state. These observations provide new information about the reactivity of DAE derivatives following excitation in the visible and near-UV.

Intramolecular charge transfer of a push-pull chromophore with restricted internal rotation of an electron donor

Intramolecular charge transfer (ICT) of 4-(dicyanomethylene)-2-methyl-6-[2-(2,3,6,7-tetrahydro-1H,5H-benzo[ij]quinolizin-9-yl)vinyl]-4H-pyran (LD688) in DMSO solution was investigated by femtosecond stimulated Raman spectroscopy (FSRS) with 403 nm excitation. The molecular structure of LD688 is similar to that of a well-known push-pull chromophore, 4-dicyanomethylene-2-methyl-6-(p-dimethylaminostyryl)-4H-pyran (DCM), except that the internal rotation of the electron-donating dimethylamino group is restricted with the introduction of the julolidine moiety. Upon photo-excitation, LD688 shows an ultrafast (1.0 ps) ICT followed by the vibrational relaxation (3-8 ps) in the charge-transfer (CT) state. Two distinct Raman spectra of LD688 in the locally excited (LE) and CT state of the S₁ state were retrieved from FSRS measurements. Based on the time-dependent density functional theory (TDDFT) simulations, a "twisted" julolidine geometry of LD688 was proposed for the ICT state, which was further confirmed in comparison to the spectral changes of several push-pull chromophores with the π-conjugated backbone of stilbene, biphenyl, styrylpyran, styrylpyridinium, and styrene in terms of the skeletal vibrational modes of ν_19b,py, ν_CC,ph, and ν_CN.

Mode-Specific Vibrational Analysis of Exciton Delocalization and Structural Dynamics in Conjugated Oligomers

Exciton delocalization in organic semiconducting polymers, affected by structures at a molecular level, plays a crucial role in modulating relaxation pathways, such as charge generation and singlet fission, which can boost device efficiency. However, the structural diversity of polymers and broad signals from typical electronic spectroscopy have their limits when it comes to revealing the interplay between local structures and exciton delocalization. To tackle these problems, we apply femtosecond stimulated Raman spectroscopy in archetypical conjugated oligothiophenes with different chain lengths. We observed Raman frequency dispersions of symmetric bond stretching modes and mode-specific kinetics in the S₁ Raman spectra, which underpins the subtle and complex interplay between exciton delocalization and bond length alternation along the conjugation coordinate. Our results provide a more general picture of exciton delocalization in the context of molecular structures for conjugated materials.

Transient electronic and vibrational signatures during reversible photoswitching of a cyanobacteriochrome photoreceptor

Cyanobacteriochromes (CBCRs) are an emerging class of photoreceptors that are distant relatives of the phytochromes family. Unlike phytochromes, CBCRs have gained popularity in optogenetics due to their highly diverse spectral properties spanning the UV to near-IR region and only needing a single compact binding domain. AnPixJg2 is a CBCR that can reversibly photoswitch between its red-absorbing (^15ZP_r) and green-absorbing (^15EP_g) forms of the phycocyanobilin (PCB) cofactor. To reveal primary events of photoconversion, we implemented femtosecond transient absorption spectroscopy with a homemade LED box and a miniature peristaltic pump flow cell to track transient electronic responses of the photoexcited AnPixJg2 on molecular time scales. The 525 nm laser-induced P_g-to-P_r reverse conversion exhibits a ~3 ps excited-state lifetime before reaching the conical intersection (CI) and undergoing further relaxation on the 30 ps time scale to generate a long-lived Lumi-G ground state intermediate en route to P_r. The 650 nm laser-induced P_r-to-P_g forward conversion is less efficient than reverse conversion, showing a longer-lived excited state which requires two steps with ~13 and 217 ps time constants to enter the CI region. Furthermore, using a tunable ps Raman pump with broadband Raman probe on both the Stokes and anti-Stokes sides, we collected the pre-resonance ground-state femtosecond stimulated Raman spectroscopy (GS-FSRS) data with mode assignments aided by quantum calculations. Key vibrational marker bands at ~850, 1050, 1615, and 1649 cm^-1 of the P_r conformer exhibit a notable blueshift to those of the P_g conformer inside AnPixJg2, reflecting the PCB chromophore terminal D (major) and A (minor) ring twist along the primary photoswitching reaction coordinate. This integrated ultrafast spectroscopy and computational platform has the potential to elucidate photochemistry and photophysics of more CBCRs and photoactive proteins in general, providing the highly desirable mechanistic insights to facilitate the rational design of functional molecular sensors and devices.

Accessing Excited State Molecular Vibrations by Femtosecond Stimulated Raman Spectroscopy

Excited state vibrations are crucial for determining the photophysical and photochemical properties of molecular compounds. Stimulated Raman scattering can coherently stimulate and probe molecular vibrations with optical pulses, but it is generally restricted to ground state properties. Working under resonance conditions enables cross-section enhancement and selective excitation to a targeted electronic level but is hampered by an increased signal complexity due to the presence of overlapping spectral contributions. Here, we show how detailed information about ground and excited state vibrations can be disentangled by exploiting the relative time delay between Raman and probe pulses to control the excited state population, combined with a diagrammatic formalism to dissect the pathways concurring with the signal generation. The proposed method is then exploited to elucidate the vibrational properties of the ground and excited electronic states in the paradigmatic case of cresyl violet. We anticipate that the presented approach holds the potential for selective mapping of the reaction coordinates pertaining to transient electronic stages implied in photoactive compounds.

Ab Initio Approach to Femtosecond Stimulated Raman Spectroscopy: Investigating Vibrational Modes Probed in Excited-State Relaxation of Quaterthiophenes

Femtosecond stimulated Raman spectroscopy (FSRS) is an ultrafast pump-probe technique designed to elucidate excited-state molecular dynamics by means of vibrational spectroscopy. We present a first-principles protocol for the simulation of FSRS that integrates ab initio molecular dynamics with computational resonance Raman spectroscopy. Theoretical calculations can monitor the time-dependent evolution of specific vibrational modes and thus provide insight into the nature of the motion responsible for the experimental FSRS signal, and we apply this technique to study quaterthiophene derivatives. The S₁ state of two different quaterthiophene derivatives relaxes via in-phase and out-of-phase stretching modes whose frequencies are coupled to the dihedral backbone angle, such that the spectral evolution reflects the excited-state relaxation toward a planar conformation. The simulated spectra aid in confirming the experimental assignment of the vibrational modes that are probed in the existing FSRS experiments on quaterthiophenes.

Simulation of time-resolved x-ray absorption spectroscopy of ultrafast dynamics in particle-hole-excited 4-(2-thienyl)-2,1,3-benzothiadiazole

To date, alternating co-polymers based on electron-rich and electron-poor units are the most attractive materials to control functionality of organic semiconductor layers in which ultrafast excited-state processes play a key role. We present a computational study of the photoinduced excited-state dynamics of the 4-(2-thienyl)-2,1,3-benzothiadiazole (BT-1T) molecule, which is a common building block in the backbone of π-conjugated polymers used for organic electronics. In contrast to homo-polymer materials, such as oligothiophene, BT-1T has two non-identical units, namely, thiophene and benzothiadiazole, making it attractive for intramolecular charge transfer studies. To gain a thorough understanding of the coupling of excited-state dynamics with nuclear motion, we consider a scenario based on femtosecond time-resolved x-ray absorption spectroscopy using an x-ray free-electron laser in combination with a synchronized ultraviolet femtosecond laser. Using Tully's fewest switches surface hopping approach in combination with excited-state calculations at the level of configuration interaction singles, we calculate the gas-phase x-ray absorption spectrum at the carbon and nitrogen K edges as a function of time after excitation to the lowest electronically excited state. The results of our time-resolved calculations exhibit the charge transfer driven by non-Born-Oppenheimer physics from the benzothiadiazole to thiophene units during relaxation to the ground state. Furthermore, our ab initio molecular dynamics simulations indicate that the excited-state relaxation processes involve bond elongation in the benzothiadiazole unit as well as thiophene ring puckering at a time scale of 100 fs. We show that these dynamical trends can be identified from the time-dependent x-ray absorption spectrum.

On the Discrepancy between Experimental and Calculated Raman Intensities for Conjugated Phenyl and Thiophene Derivatives

Compared with experimental spectra, calculations for conjugated phenyl and thiophene oligomers tend to overestimate the ground state Raman intensities of higher-frequency vibrations (1200-1800 cm^-1) relative to the intensities at lower frequencies (<1200 cm^-1). The discrepancy was observed in previous benchmarking work that examined the method dependence of the calculated Raman spectra for a series of aromatic molecules. This paper further investigates the nature of the discrepancy by examining the role of anharmonic corrections and the dependence of the calculated Raman spectra on the inter-ring torsion angle for the representative molecules biphenyl (BP), 2-phenylthiophene (PT), and 2,2'-bithiophene (BT). Perturbative anharmonic corrections to the spectra calculated using density functional theory (DFT) provide only slightly better agreement with experiment. On the other hand, calculations at larger torsion angles give up to 30% improvement in the relative Raman intensities compared with the spectra calculated at the optimized geometries. The torsion-angle dependence of the Raman intensities is most pronounced for delocalized C-C and C-S stretching modes, and less pronounced for bending and ring distortion modes that do not involve inter-ring stretching. Higher-level calculations using the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] method indicate that DFT underestimates the energy barrier for torsion isomerization at small angles, and it overestimates the barriers at large angles, thus predicting minimum geometries at torsion angles that are too small. Therefore, the results suggest that the discrepancy in relative Raman intensities may be related to an overestimation of inter-ring conjugation by DFT, which also tends to favor geometries that are too planar.

Femtosecond-to-nanosecond dynamics of flavin mononucleotide monitored by stimulated Raman spectroscopy and simulations

Flavin mononucleotide (FMN) belongs to the large family of flavins, ubiquitous yellow-coloured biological chromophores that contain an isoalloxazine ring system. As a cofactor in flavoproteins, it is found in various enzymes and photosensory receptors, like those featuring the light-oxygen-voltage (LOV) domain. The photocycle of FMN is triggered by blue light and proceeds via a cascade of intermediate states. In this work, we have studied isolated FMN in an aqueous solution in order to elucidate the intrinsic electronic and vibrational changes of the chromophore upon excitation. The ultrafast transitions of excited FMN were monitored through the joint use of femtosecond stimulated Raman spectroscopy (FSRS) and transient absorption spectroscopy encompassing a time window between 0 ps and 6 ns with 50 fs time resolution. Global analysis of the obtained transient visible absorption and transient Raman spectra in combination with extensive quantum chemistry calculations identified unambiguously the singlet and triplet FMN populations and addressed solvent dynamics effects. The good agreement between the experimental and theoretical spectra facilitated the assignment of electronic transitions and vibrations. Our results represent the first steps towards more complex experiments aimed at tracking structural changes of FMN embedded in light-inducible proteins upon photoexcitation.

Probing Dynamics in Higher-Lying Electronic States with Resonance-Enhanced Femtosecond Stimulated Raman Spectroscopy

Femtosecond stimulated Raman scattering (FSRS) measurements typically probe the structural dynamics of a molecule in the first electronically excited state, S₁. While these measurements often rely on an electronic resonance condition to increase signal strength or enhance species selectivity, the effects of the resonance condition are usually neglected. However, mode-specific enhancements of the vibrational transitions in an FSRS spectrum contain detailed information about the resonant (upper) electronic state. Analogous to ground-state resonance Raman spectroscopy, the relative intensities of the Raman bands reveal displacements of the upper potential energy surface due to changes in the bonding pattern upon S _n ← S₁ electronic excitation, and therefore provide a sensitive probe of the ultrafast dynamics in the higher-lying state, S _n. Raman gain profiles from the wavelength-dependent FSRS spectrum of the model compound 2,5-diphenylthiophene (DPT) reveal several modes with large displacement in the upper potential energy surface, including strong enhancement of a delocalized C-S-C stretching and ring deformation mode. The experimental results provide a benchmark for comparison with calculated spectra using time-dependent density functional theory (TD-DFT) and equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD), where the calculations are based on the time-dependent formalism for resonance Raman spectroscopy. The simulated spectra are obtained from S₁-S _n transition strengths and the energy gradients of the upper (S _n) potential energy surfaces along the S₁ normal mode coordinates. The experimental results provide a stringent test of the computational approach, and indicate important limitations based on the level of theory and basis set. This work provides a foundation for making more accurate assignments of resonance-enhanced excited-state Raman spectra, as well as extracting novel information about higher-lying excited states in the transient absorption spectrum of a molecule.