Vibrationally resolved electronic spectra including vibrational pre-excitation: Theory and application to VIPER spectroscopy

Vibrational-Mode-Selective Modulation of Electronic Excitation

Vibrational-mode-selective modulation of electronic excitation is conducted with a synchronized femtosecond (fs) visible (vis) pulse and a picosecond (ps) infrared (IR) pulse. The mechanism of modulation of vibrational and vibronic relaxation behavior of excited state is investigated with ultrafast vis/IR, IR/IR, and vis-IR/IR transient spectroscopy, optical gating experiments and theoretical calculations. An organic molecule, 4'-(N,N-dimethylamino)-3-methoxyflavone (DMA3MHF) is chosen as the model system. Upon 1608 cm-1 excitation, the skeleton stretching vibration of DMA3MHF is energized, which can significantly change the shape of the absorption, facilitate the radiative decay and promote emission from vibrational excited states. As results, a remarkable enhancement and a slight blueshift in fluorescence are observed. The mode-selective modulation of electronic excitation is not limited in luminescence or photophysics. It is expected to be widely applicable in tuning many photochemical processes.

Blue-Shifted and Broadened Fluorescence Enhancement by Visible and Mode-Selective Infrared Double Excitations

Mode-selective vibrational excitations to modify the electronic states of fluorescein dianion in methanol solutions are carried out with a femtosecond visible pulse synchronized with a tunable high-power, narrow-band picosecond infrared (IR) pulse. In this work, simultaneous intensity enhancement, peak blueshift, and line width broadening of fluorescence are observed in the visible/IR double resonance experiments. Comprehensive investigations on the modulation mechanism with scanning the vibrational excitation frequencies, tuning the time delay between the two excitation pulses, theoretical calculations, and nonlinear and linear spectroscopic measurements suggest that the fluorescence intensity enhancement is caused by the increase of the Franck-Condon factor induced by the vibrational excitations at the electronic ground state. Various enhancement effects are observed as vibrations initially excited by the IR photons relax to populate the vibrational modes of lower frequencies. The peak blueshift and line width broadening are caused by both increasing the Franck-Condon factors among different subensembles because of IR pre-excitation and the long-lived processes induced by the initial IR excitation. The results demonstrate that the fluorescence from the visible/IR double resonance experiments is not a simple sum frequency effect, and vibrational relaxations can produce profound effects modifying luminescence.

Molecular factors determining brightness in fluorescence-encoded infrared vibrational spectroscopy

Fluorescence-encoded infrared (FEIR) spectroscopy is a recently developed technique for solution-phase vibrational spectroscopy with detection sensitivity at the single-molecule level. While its spectroscopic information content and important criteria for its practical experimental optimization have been identified, a general understanding of the electronic and nuclear properties required for highly sensitive detection, i.e., what makes a molecule a "good FEIR chromophore," is lacking. This work explores the molecular factors that determine FEIR vibrational activity and assesses computational approaches for its prediction. We employ density functional theory (DFT) and its time-dependent version (TD-DFT) to compute vibrational and electronic transition dipole moments, their relative orientation, and the Franck-Condon factors involved in FEIR activity. We apply these methods to compute the FEIR activities of normal modes of chromophores from the coumarin family and compare these predictions with experimental FEIR cross sections. We discuss the extent to which we can use computational models to predict the FEIR activity of individual vibrations in a candidate molecule. The results discussed in this work provide the groundwork for computational strategies for choosing FEIR vibrational probes or informing the structure of designer chromophores for single-molecule spectroscopic applications.

Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.

Theoretical Insight into the Effect of Phosphorus Oxygenation on Nonradiative Decays: Comparative Analysis of P-Bridged Stilbene Analogs

Incorporation of the phosphorus element into a π-conjugated skeleton offers valuable prospects for adjusting the electronic structure of the resulting functional π-electron systems. Trivalent phosphorus has the potential to decrease the LUMO level through σ*-π* interaction, which is further enhanced by its oxygenation to the pentavalent P center. This study shows that utilizing our computational analysis to examine excited-state dynamics based on radiative/nonradiative rate constants and fluorescence quantum yield (ΦF) is effective for analyzing the photophysical properties of P-containing organic dyes. We theoretically investigate how the trivalent phosphanyl group and pentavalent phosphine oxide moieties affect radiative and nonradiative decay processes. We evaluate four variations of P-bridged stilbene analogs. Our analysis reveals that the primary decay pathway for photoexcited bis-phosphanyl-bridged stilbene is the intersystem crossing (ISC) to the triplet state and nonradiative. The oxidation of the phosphine moiety, however, suppresses the ISC due to the relative destabilization of the triplet states. The calculated rate constants match an increase in experimental ΦF from 0.07 to 0.98, as simulated from 0.23 to 0.94. The reduced HOMO-LUMO gap supports a red shift in the fluorescence spectra relative to the phosphine analog. The thiophene-fused variant with the nonoxidized trivalent P center exhibits intense emission with a high ΦF, 0.95. Our prediction indicates that the ISC transfer is obstructed owing to the relatively destabilized triplet state induced by the thiophene substitution. Conversely, the thiophene-fused analog with the phosphine oxide moieties triggers a high-rate internal conversion mediated by conical intersection, leading to a decreased ΦF.

Multimode vibrational dynamics and orientational effects in fluorescence-encoded infrared spectroscopy. II. Analysis of early-time signals

Developing fluorescence-encoded infrared (FEIR) vibrational spectroscopy for single-molecule applications requires a detailed understanding of how the molecular response and external experimental parameters manifest in the detected signals. In Paper I [L. Whaley-Mayda, A. Guha, and A. Tokmakoff, J. Chem. Phys. 159, 194201 (2023)] we introduced a nonlinear response function theory to describe vibrational dynamics, vibronic coupling, and transition dipole orientation in FEIR experiments with ultrashort pulses. In this second paper, we apply the theory to investigate the role of intermode vibrational coherence, the orientation of vibrational and electronic transition dipoles, and the effects of finite pulse durations in experimental measurements. We focus on measurements at early encoding delays-where signal sizes are largest and therefore of most value for single-molecule experiments, but where many of these phenomena are most pronounced and can complicate the appearance of data. We compare experiments on coumarin dyes with finite-pulse response function simulations to explain the time-dependent behavior of FEIR spectra. The role of the orientational response is explored by analyzing polarization-dependent experiments and their ability to resolve relative dipole angles in the molecular frame. This work serves to demonstrate the molecular information content of FEIR experiments, and develop insight and guidelines for their interpretation.

Multimode vibrational dynamics and orientational effects in fluorescence-encoded infrared spectroscopy. I. Response function theory

Fluorescence-encoded infrared (FEIR) spectroscopy is an emerging technique for performing vibrational spectroscopy in solution with detection sensitivity down to single molecules. FEIR experiments use ultrashort pulses to excite a fluorescent molecule's vibrational and electronic transitions in a sequential, time-resolved manner, and are therefore sensitive to intervening vibrational dynamics on the ground state, vibronic coupling, and the relative orientation of vibrational and electronic transition dipole moments. This series of papers presents a theoretical treatment of FEIR spectroscopy that describes these phenomena and examines their manifestation in experimental data. This first paper develops a nonlinear response function description of Fourier-transform FEIR experiments for a two-level electronic system coupled to multiple vibrations, which is then applied to interpret experimental measurements in the second paper [L. Whaley-Mayda et al., J. Chem. Phys. 159, 194202 (2023)]. Vibrational coherence between pairs of modes produce oscillatory features that interfere with the vibrations' population response in a manner dependent on the relative signs of their respective Franck-Condon wavefunction overlaps, leading to time-dependent distortions in FEIR spectra. The orientational response of population and coherence contributions are analyzed and the ability of polarization-dependent experiments to extract relative transition dipole angles is discussed. Overall, this work presents a framework for understanding the full spectroscopic information content of FEIR measurements to aid data interpretation and inform optimal experimental design.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Non-Phenomenological Description of the Time-Resolved Emission in Solution with Quantum-Classical Vibronic Approaches-Application to Coumarin C153 in Methanol

We report a joint experimental and theoretical work on the steady-state spectroscopy and time-resolved emission of the coumarin C153 dye in methanol. The lowest energy excited state of this molecule is characterized by an intramolecular charge transfer thus leading to remarkable shifts of the time-resolved emission spectra, dictated by the methanol reorganization dynamics. We selected this system as a prototypical test case for the first application of a novel computational protocol aimed at the prediction of transient emission spectral shapes, including both vibronic and solvent effects, without applying any phenomenological broadening. It combines a recently developed quantum-classical approach, the adiabatic molecular dynamics generalized vertical Hessian method (Ad-MD|gVH), with nonequilibrium molecular dynamics simulations. For the steady-state spectra we show that the Ad-MD|gVH approach is able to reproduce quite accurately the spectral shapes and the Stokes shift, while a ∼0.15 eV error is found on the prediction of the solvent shift going from gas phase to methanol. The spectral shape of the time-resolved emission signals is, overall, well reproduced, although the simulated spectra are slightly too broad and asymmetric at low energies with respect to experiments. As far as the spectral shift is concerned, the calculated spectra from 4 ps to 100 ps are in excellent agreement with experiments, correctly predicting the end of the solvent reorganization after about 20 ps. On the other hand, before 4 ps solvent dynamics is predicted to be too fast in the simulations and, in the sub-ps timescale, the uncertainty due to the experimental time resolution (300 fs) makes the comparison less straightforward. Finally, analysis of the reorganization of the first solvation shell surrounding the excited solute, based on atomic radial distribution functions and orientational correlations, indicates a fast solvent response (≈100 fs) characterized by the strengthening of the carbonyl-methanol hydrogen bond interactions, followed by the solvent reorientation, occurring on the ps timescale, to maximize local dipolar interactions.

Choose your leaving group: selective photodeprotection in a mixture of pHP-caged compounds by VIPER excitation

Photocages are light-triggerable molecular moieties that can locally release a pre-determined leaving group (LG). Finding a suitable photocage for a particular application may be challenging, as the choice may be limited by for instance the optical or physicochemical properties of the system. Using more than one photocage to release different LGs in a reaction mixture may even be more difficult. In this work an experimental strategy is presented that allows us to hand-pick the release of different LGs, and to do so in any order. This is achieved by using isotopologue photocage-LG mixtures in combination with ultrafast VIbrationally Promoted Electronic Resonance (VIPER) excitation. The latter provides the required molecular selectivity simply by tuning the wavenumber of the used IR pulses to the resonance of a specific photocage isotopologue, as is demonstrated here for the para-hydroxyphenacyl (pHP) photocage. For spectroscopic convenience, we use isotopologues of the infrared (IR) spectroscopic marker -SCN as different LGs. Especially for applications where fast LG release is required, pHP is found to be an excellent candidate, as free LG formation is observed to occur with a 10 ps lifetime. The devised strategy may open up new complex uncaging applications, where multiple LGs can be formed locally on a short time scale and in any sequence.

Vibrationally resolved two-photon electronic spectra including vibrational pre-excitation: Theory and application to VIPER spectroscopy with two-photon excitation

Following up on our previous work on vibrationally resolved electronic absorption spectra including the effect of vibrational pre-excitation [von Cosel et al., J. Chem. Phys. 147, 164116 (2017)], we present a combined theoretical and experimental study of two-photon-induced vibronic transitions in polyatomic molecules that are probed in the VIbrationally Promoted Electronic Resonance experiment using two-photon excitation (2P-VIPER). In order to compute vibronic spectra, we employ time-independent and time-dependent methods based on the evaluation of Franck-Condon overlap integrals and Fourier transformations of time-domain correlation functions, respectively. The time-independent approach uses a generalized version of the FCclasses method, while the time-dependent approach relies on the analytical evaluation of Gaussian moments within the harmonic approximation, including Duschinsky rotation effects. For the Coumarin 6 dye, two-dimensional 2P-VIPER experiments involving excitation to the lowest-lying singlet excited state (S₁) are presented and compared with corresponding one-photon VIPER spectra. In both cases, coumarin ring modes and a CO stretch mode show VIPER activity, albeit with different relative intensities. Selective pre-excitation of these modes leads to a pronounced redshift of the low-frequency edge of the electronic absorption spectrum, which is a prerequisite for the VIPER experiment. Theoretical analysis underscores the role of interference between Franck-Condon and Herzberg-Teller effects in the two-photon experiment, which is at the root of the observed intensity distribution.

FCclasses3: Vibrationally-resolved spectra simulated at the edge of the harmonic approximation

We introduce FCclasses3, a code to carry out vibronic simulations of electronic spectra and nonradiative rates, based on the harmonic approximation. Key new features are: implementation of the full family of vertical and adiabatic harmonic models, vibrational analysis in curvilinear coordinates, extension to several electronic spectroscopies and implementation of time-dependent approaches. The use of curvilinear valence internal coordinates allows the adoption of quadratic model potential energy surfaces (PES) of the initial and final states expanded at arbitrary configurations. Moreover, the implementation of suitable projectors provides a robust framework for defining reduced-dimensionality models by sorting flexible coordinates out of the harmonic subset, so that they can then be treated at anharmonic level, or with mixed quantum classical approaches. A set of tools to facilitate input preparation and output analysis is also provided. We show the program at work in the simulation of different spectra (one and two-photon absorption, emission and resonance Raman) and internal conversion rate of a typical rigid molecule, anthracene. Then, we focus on absorption and emission spectra of a series of flexible polyphenyl molecules, highlighting the relevance of some of the newly implemented features. The code is freely available at http://www.iccom.cnr.it/en/fcclasses/.

Resonance conditions, detection quality, and single-molecule sensitivity in fluorescence-encoded infrared vibrational spectroscopy

Fluorescence-encoded Infrared (FEIR) spectroscopy is a vibrational spectroscopy technique that has recently demonstrated the capability of single-molecule sensitivity in solution without near-field enhancement. This work explores the practical experimental factors that are required for successful FEIR measurements in both the single-molecule and bulk regimes. We investigate the role of resonance conditions by performing measurements on a series of coumarin fluorophores of varying electronic transition frequencies. To analyze variations in signal strength and signal to background between molecules, we introduce an FEIR brightness metric that normalizes out measurement-specific parameters. We find that the effect of the resonance condition on FEIR brightness can be reasonably well described by the electronic absorption spectrum. We discuss strategies for optimizing detection quality and sensitivity in bulk and single-molecule experiments.

Symmetry-Breaking Enhanced Herzberg-Teller Effect with Brominated Polyacenes

Molecular symmetry is vital to the selection rule of vibrationally resolved electronic transition, particularly when the nuclear dependence of electronic wave function is explicitly treated by including Franck-Condon (FC) factor, Franck-Condon/Herzberg-Teller (FC/HT) interference, and Herzberg-Teller (HT) coupling. Our present study investigated the light absorption spectra of highly symmetric tetracene, pentacene, and hexacene molecules of point-group D_2h, as well as their monobrominated derivatives with a lower C_s symmetry. It was found that the symmetry-breaking monobromination allows more vibrational normal modes and their pairs to contribute to FC/HT interference and HT coupling, respectively. Through a projection of a molecule's vibrational normal modes to its irreducible representations, a linear relationship between the FC/HT intensity to the polyacene's size was deduced alongside a quadratic dependence of the HT intensity. Both theoretically derived correlations were well justified by our numerical simulations, which also demonstrated an approximately 20% improvement on the agreement with experimental line shape if the HT theory is adopted to replace the FC approximation. Moreover, for these low-symmetry monobrominated polyacenes, the FC intensity was even weaker than its FC/HT and HT counterparts at some excitation energies, making the HT theory imperative to decipher vibronic coupling, a fundamental driving force behind numerous chemical, biological, and photophysical processes.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Quantum dissipative systems beyond the standard harmonic model: Features of linear absorption and dynamics

Current simulations of ultraviolet-visible absorption lineshapes and dynamics of condensed phase systems largely adopt a harmonic description to model vibrations. Often, this involves a model of displaced harmonic oscillators that have the same curvature. Although convenient, for many realistic molecular systems, this approximation no longer suffices. We elucidate nonstandard harmonic and anharmonic effects on linear absorption and dynamics using a stochastic Schrödinger equation approach to account for the environment. First, a harmonic oscillator model with ground and excited potentials that differ in curvature is utilized. Using this model, it is shown that curvature difference gives rise to an additional substructure in the vibronic progression of absorption spectra. This effect is explained and subsequently quantified via a derived expression for the Franck-Condon coefficients. Subsequently, anharmonic features in dissipative systems are studied, using a Morse potential and parameters that correspond to the diatomic molecule H₂ for differing displacements and environment interaction. Finally, using a model potential, the population dynamics and absorption spectra for the stiff-stilbene photoswitch are presented and features are explained by a combination of curvature difference and anharmonicity in the form of potential energy barriers on the excited potential.

Generating Function Approach to Single Vibronic Level Fluorescence Spectra

An efficient time-dependent generating function method to compute vibronic emission and absorption spectra arising from transitions from a singly excited vibrational initial state is presented. In contrast to existing finite temperature approaches that intrinsically contain these transitions weighted by a Boltzmann factor, the current approach allows one to calculate these transitions individually. Using vibrational frequencies and normal modes computed by the second-order approximate coupled cluster (CC2) method, this formalism is used to compute the single vibronic level (SVL) fluorescence spectra of anthracene. Calculated spectra are in excellent agreement with spectra measured in jet-cooled expansion experiments. Duschinsky mixing is necessary to explain intensities of certain peaks. In a few cases, CC2, however, underestimates Duschinsky mixing, leading to too low peak intensities. An empirical correction of the Duschinsky matrix is presented. The presented method has the potential to facilitate the assignment and interpretation of SVL fluorescence spectra.

Vibrationally Broadened Optical Spectra of Selected Radicals and Cations Derived from Adamantane: A Time-Dependent Correlation Function Approach

Diamondoids are hydrogen-saturated molecular motifs cut out of diamond, forming a class of materials with tunable optoelectronic properties. In this work, we extend previous work on neutral, closed-shell diamondoids by computing with hybrid density functional theory and time-dependent correlation functions vibrationally broadened absorption spectra of cations and radicals derived from the simplest diamondoid, adamantane, namely, the neutral 1- and 2-adamantyl radicals (C₁₀H₁₅), the 1- and 2-adamantyl cations (C₁₀H₁₅⁺), and the adamantane radical cation (C₁₀H₁₆⁺). For selected cases, we also report vibrationally broadened emission, photoelectron, and resonance Raman spectra. Furthermore, the effect of the damping factor on the vibrational fine-structure is studied. The following trends are found: (1) Low-energy absorptions of the adamantyl radicals and cations, and of the adamantane cation, are all strongly red-shifted with respect to adamantane; (2) also, emission spectra are strongly red-shifted, whereas photoelectron spectra are less affected for the cases studied; (3) vibrational fine-structures are reduced compared to those of adamantane; (4) the spectroscopic signals of 1- and 2-adamantyl species are significantly different from each other; and (5) reducing the damping factor has only a limited effect on the vibrational fine-structure in most cases. This suggests that removing hydrogen atoms and/or electrons from adamantane leads to new optoelectronic properties, which should be detectable by vibronic spectroscopy.

Fluorescence-Encoded Infrared Spectroscopy: Ultrafast Vibrational Spectroscopy on Small Ensembles of Molecules in Solution

Fluorescence-encoded infrared (FEIR) spectroscopy is an ultrafast technique that uses a visible pulse to up-convert information about IR-driven vibrations into a fluorescent electronic population. Here we present an updated experimental approach to FEIR that achieves high sensitivity through confocal microscopy, high repetition rate excitation, and single-photon counting. We demonstrate the sensitivity of our experiment by measuring ultrafast vibrational transients and Fourier transform spectra of increasingly dilute solutions of a coumarin dye. We collect high-quality data at 40 μM (∼2 orders of magnitude below the limit for conventional IR) and make measurements down to the 10-100 nM range (∼5 orders of magnitude) before background signals become overwhelming. At 10 nM we measure the average number of molecules in the focal volume to be ∼20 using fluorescence correlation spectroscopy. This level of sensitivity opens up the possibility of performing fluctuation correlation vibrational spectroscopy or-with further improvement-single-molecule measurements.

Infrared pre-excitation grants isotopomer-specific photochemistry

Species-selective photochemistry is often hampered by overlapping UV-Vis spectra. We overcome this long-standing problem by combined vibrational and electronic excitation as demonstrated by isotopomer selection. The influence of various factors on selectivity is discussed.

Controlling Photochemistry via Isotopomers and IR Pre-excitation

It is a photochemist's dream to be able to photoinduce a reaction of a specific molecular species in an ensemble of similar but not identical ones. The problem is that similar molecules often exhibit nearly identical UV-Vis absorption spectra, making them difficult or impossible to distinguish or to select spectroscopically. The ultrafast VIPER (VIbrationally Promoted Electronic Resonance) pulse sequence allows to pick a single species for electronic excitation based on its infrared spectrum. The latter usually shows more features, allowing the discrimination between species than the UV-Vis spectrum. Here, we show that it is possible to induce and monitor species-selective photochemistry even for molecules with virtually identical UV-Vis spectra, which is the case for isotopomers. Next to isotope-selective photochemistry in solution, applications to orthogonal photo-uncaging and species-selective spectroscopy and photochemistry in mixtures are within reach.