Quantum chemical comparison of vertical, adiabatic, and 0-0 excitation energies: The PYP and GFP chromophores

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

Computational Study of the Influence of Nitrogen-Containing Unsaturated Heterocyclic Substituents on Electronic and Spectroscopic Properties of Squaraine Derivatives

A comprehensive understanding of the influence of the different structural elements for the molecular properties is crucial for improving the material design procedures in the field of photosensitive dyes. The present study provides a detailed analysis of the influence of the number of heteroatoms (nitrogens) and the lengths of the π-electron skeleton on the one- and two-photon absorption of the symmetric squaraine dyes. Extended computational study covers the conventional vertical excitation calculations within the TD-DFT formalisms as well as several advanced single-reference methods including double excitations such as CIS(D), SAC-CI, and ADC(2). Additionally, the weaknesses of the vertical approach are investigated by including the geometry relaxation upon excitation via adiabatic and 0-0 treatment.

Does the wavelength dependent photoisomerization process of the p‑coumaric acid come out from the electronic state dependent pathways?

Similar to the anion photoactive yellow protein (PYP) chromophore, the neutral form of the PYP chromophore was also found to exhibit a the wavelength-dependent photoisomerization quantum yield. The isomerization quantum yield increases with the increasing excitation energy on the S₁ state, while decreases when being excited to the S₂ state. Does this wavelength dependent product yield come out from the specific reaction pathways of the S₁ and S₂ states? This would mean that, the relaxation pathway of the S₂ state is distinct from that of the S₁ state and does not involve twisting motion. Does it break Kasha's rule by exhibiting a direct transition from the S₂ state to the ground state? The underlying mechanism needs further in. In this article, we employed the on-the-fly dynamics simulations and static electronic structure calculations to reveal the deactivation mechanism of the neutral form of the PYP chromophore. Our results indicated that the CC twisting motion dominates the S₁ state decay process. In contrast, for the decay process of the S₂ state, an ultrafast transition from the S₂ to the S₁ state through a planar conical intersection is observed, and the excess energy activates a new reaction channel to the ground state characterized by a puckering distortion of the ring. This pathway competes with the photoisomerization channel. No direct transition from S₂ to S₀ is observed, hence Kasha's rule is valid for this process. Our calcualtions can provide a reasonable explanation of the wavelength-dependent isomerization quantum yield of neutral PYP chromophore, and we hope it can provide theoretical foundations for comparing the effect of protonation state on the dynamcal behaviors of PYP chromophore.

Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal?

To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S₀ →S₁ and S₀ →S₂ energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.

Why do TD-DFT excitation energies of BODIPY/Aza-BODIPY families largely deviate from experiment? Answers from electron correlated and multireference methods

The vertical excitation energies of 17 boron-dipyrromethene (BODIPY) core structures with a variety of substituents and ring sizes are benchmarked using time-dependent density functional theory (TD-DFT) with nine different functionals combined with the cc-pVTZ basis set. When compared to experimental measurements, all functionals provide mean absolute errors (mean AEs) greater than 0.3 eV, larger than the 0.1-0.3 eV differences typically expected from TD-DFT. Due to the high linear correlation of TD-DFT results with experiment, most functionals can be used to predict excitation energies if corrected empirically. Using the CAM-B3LYP functional, 0-0 transition energies are determined, and while the absolute difference is improved (mean AE = 0.478 eV compared to 0.579 eV), the correlation diminishes substantially (R(2) = 0.961 to 0.862). Two very recently introduced charge transfer (CT) indices, q(CT) and d(CT), and electron density difference (EDD) plots demonstrate that CT does not play a significant role for most of the BODIPYs examined and, thus, cannot be the source of error in TD-DFT. To assess TD-DFT methods, vertical excitation energies are determined utilizing TD-HF, configuration interaction CIS and CIS(D), equation of motion EOM-CCSD, SAC-CI, and Laplace-transform based local coupled-cluster singles and approximate doubles LCC2* methods. Moreover, multireference CASSCF and CASPT2 vertical excitation energies were also obtained for all species (except CASPT2 was not feasible for the four largest systems). The SAC-CI/cc-pVDZ, LCC2*/cc-pVDZ, and CASPT2/cc-pVDZ approaches are shown to have the smallest mean AEs of 0.154, 0.109, and 0.100 eV, respectively; the utility of the LCC2* approach is demonstrated for eight extended BODIPYs and aza-BODIPYs. We found that the problems with TD-DFT arise from difficulties in dealing with the differential electron correlation (as assessed by comparing CCS, CC2, LR-CCSD, CCSDR(T), and CCSDR(3) vertical excitation energies for five compounds) and from contributions of multireference character and double excitations (from analysis of the CASSCF wave functions).

Controlling electron emission from the photoactive yellow protein chromophore by substitution at the coumaric acid group

Understanding how the interactions between a chromophore and its surrounding protein control the function of a photoactive protein remains a challenge. Here, we present the results of photoelectron spectroscopy measurements and quantum chemistry calculations aimed at investigating how substitution at the coumaryl tail of the photoactive yellow protein chromophore controls competing relaxation pathways following photoexcitation of isolated chromophores in the gas phase with ultraviolet light in the range 350-315 nm. The photoelectron spectra are dominated by electrons resulting from direct detachment and fast detachment from the 2(1)ππ* state but also have a low electron kinetic energy component arising from autodetachment from lower lying electronically excited states or thermionic emission from the electronic ground state. We find that substituting the hydrogen atom of the carboxylic acid group with a methyl group lowers the threshold for electron detachment but has very little effect on the competition between the different relaxation pathways, whereas substituting with a thioester group raises the threshold for electron detachment and appears to 'turn off' the competing electron emission processes from lower lying electronically excited states. This has potential implications in terms of tuning the light-induced electron donor properties of photoactive yellow protein.

0-0 Energies Using Hybrid Schemes: Benchmarks of TD-DFT, CIS(D), ADC(2), CC2, and BSE/GW formalisms for 80 Real-Life Compounds

The 0-0 energies of 80 medium and large molecules have been computed with a large panel of theoretical formalisms. We have used an approach computationally tractable for large molecules, that is, the structural and vibrational parameters are obtained with TD-DFT, the solvent effects are accounted for with the PCM model, whereas the total and transition energies have been determined with TD-DFT and with five wave function approaches accounting for contributions from double excitations, namely, CIS(D), ADC(2), CC2, SCS-CC2, and SOS-CC2, as well as Green's function based BSE/GW approach. Atomic basis sets including diffuse functions have been systematically applied, and several variations of the PCM have been evaluated. Using solvent corrections obtained with corrected linear-response approach, we found that three schemes, namely, ADC(2), CC2, and BSE/GW allow one to reach a mean absolute deviation smaller than 0.15 eV compared to the measurements, the two former yielding slightly better correlation with experiments than the latter. CIS(D), SCS-CC2, and SOS-CC2 provide significantly larger deviations, though the latter approach delivers highly consistent transition energies. In addition, we show that (i) ADC(2) and CC2 values are extremely close to each other but for systems absorbing at low energies; (ii) the linear-response PCM scheme tends to overestimate solvation effects; and that (iii) the average impact of nonequilibrium correction on 0-0 energies is negligible.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.