Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck-Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck-Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck-Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.

Calculating Photoabsorption Cross-Sections for Atmospheric Volatile Organic Compounds

Characterizing the photochemical reactivity of transient volatile organic compounds (VOCs) in our atmosphere begins with a proper understanding of their abilities to absorb sunlight. Unfortunately, the photoabsorption cross-sections for a large number of transient VOCs remain unavailable experimentally due to their short lifetime or high reactivity. While structure-activity relationships (SARs) have been successfully employed to estimate the unknown photoabsorption cross-sections of VOCs, computational photochemistry offers another promising strategy to predict not only the vertical electronic transitions of a given molecule but also the width and shape of the bands forming its absorption spectrum. In this work, we focus on the use of the nuclear ensemble approach (NEA) to determine the photoabsorption cross-section of four exemplary VOCs, namely, acrolein, methylhydroperoxide, 2-hydroperoxy-propanal, and (microsolvated) pyruvic acid. More specifically, we analyze the influence that different strategies for sampling the ground-state nuclear density-Wigner sampling and ab initio molecular dynamics with a quantum thermostat-can have on the simulated absorption spectra. We highlight the potential shortcomings of using uncoupled harmonic modes within Wigner sampling of nuclear density to describe flexible or microsolvated VOCs and some limitations of SARs for multichromophoric VOCs. Our results suggest that the NEA could constitute a powerful tool for the atmospheric community to predict the photoabsorption cross-section for transient VOCs.

Systematic improvement of molecular excited state calculations by inclusion of nuclear quantum motion: A mode-resolved picture and the effect of molecular size

The energies of molecular excited states arise as solutions to the electronic Schrödinger equation and are often compared to experiment. At the same time, nuclear quantum motion is known to be important and to induce a redshift of excited state energies. However, it is thus far unclear whether incorporating nuclear quantum motion in molecular excited state calculations leads to a systematic improvement of their predictive accuracy, making further investigation necessary. Here, we present such an investigation by employing two first-principles methods for capturing the effect of quantum fluctuations on excited state energies, which we apply to the Thiel set of organic molecules. We show that accounting for zero-point motion leads to much improved agreement with experiment, compared to "static" calculations that only account for electronic effects, and the magnitude of the redshift can become as large as 1.36 eV. Moreover, we show that the effect of nuclear quantum motion on excited state energies largely depends on the molecular size, with smaller molecules exhibiting larger redshifts. Our methodology also makes it possible to analyze the contribution of individual vibrational normal modes to the redshift of excited state energies, and in several molecules, we identify a limited number of modes dominating this effect. Overall, our study provides a foundation for systematically quantifying the shift of excited state energies due to nuclear quantum motion and for understanding this effect at a microscopic level.

On the origin of the shift between vertical excitation and band maximum in molecular photoabsorption

The analysis of the photoabsorption spectra of molecules shows that the band maximum is usually redshifted in comparison to the vertical excitation. We conducted a throughout analysis of this shift based on low-dimensional analytical and numerical model systems, showing that its origin is rooted in the frequency change between the ground and the excited states in multidimensional systems. Moreover, we deliver a benchmark of ab initio results for the shift based on a comparison of vertical excitations and band maxima calculated with the nuclear ensemble approach for the 28 organic molecules in the Mülheim molecular dataset. The mean value of the shift calculated over 60 transitions is 0.11 ± 0.08 eV. The mean value of the band width is 0.32 ± 0.14 eV.

Graphical abstract

Theoretical insights into UV-Vis absorption spectra of difluoroboron β-diketonates with an extended π system: An analysis based on DFT and TD-DFT calculations

The UV-Vis absorption spectra of difluoroboron β-diketonates with aromatic substituents at the β-carbon are studied thoroughly using DFT and TD-DFT with the CAM-B3LYP functional. The complicated experimental spectra of these dyes can be correctly interpreted by considering their structural features. A closer look at the calculated data shows that the conformational flexibility of these compounds markedly influences their spectral shape. For the complexes with an extended π system, several conformers with significantly different absorption spectra are present in the equilibrium mixture in solution. Introducing a donor group alters the electronic structure of the complexes, so the charge distribution asymmetry in the molecules increases and the nature of the electronic transitions changes. Thus, both types of substituents, aromatic and donor ones, affect the spectral shape. Understanding their roles may help one to explain the absorption spectra of these and similar compounds and predict their response to analytes and other factors.

The role of solvation models on the computed absorption and emission spectra: the case of fireflies oxyluciferin

Absorption and emission energies calculation covering both implicit and explicit solvation models using oxyluciferin as the case of study.

First principles investigation of the spectral properties of neutral, zwitterionic, and bis-cationic azaacenes

An in-depth investigation of the optical properties of recently-synthesized linear azaacene derivatives of various electronic nature (neutral, dicationic, and zwitterionic) is presented.

Ab Initio Prediction of Fluorescence Lifetimes Involving Solvent Environments by Means of COSMO and Vibrational Broadening

The fluorescence lifetime is a key property of fluorophores that can be utilized for microenvironment probing, analyte sensing, and multiplexing as well as barcoding applications. For the rational design of lifetime probes and barcodes, theoretical methods have been developed to enable the ab initio prediction of this parameter, which depends strongly on interactions with solvent molecules and other chemical species in the emitteŕs immediate environment. In this work, we investigate how a conductor-like screening model (COSMO) can account for variations in fluorescence lifetimes that are caused by such fluorophore-solvent interactions. Therefore, we calculate vibrationally broadened fluorescence spectra using the nuclear ensemble method to obtain distorted molecular geometries to sample the electronic transitions with time-dependent density functional theory (TDDFT). The influence of the solvent on fluorescence lifetimes is accounted for with COSMO. For example, for 4-hydroxythiazole fluorophore containing different heteroatoms and acidic and basic moieties in aprotic and protic solvents of varying polarity, this approach was compared to experimentally determined lifetimes in the same solvents. Our results demonstrate a good correlation between theoretically predicted and experimentally measured fluorescence lifetimes except for the polar solvents ethanol and acetonitrile that can specifically interact with the heteroatoms and the carboxylic acid of the thiazole derivative.

Franck-Condon Models for Simulating the Band Shape of Electronic Absorption Spectra

Band shape is an essential ingredient in the simulation of electronic absorption spectra. The excitation of multiple series of vibrational levels during an electronic excitation is a main contributor to band shapes. Here we present two simple models based on the Franck-Condon displaced-harmonic-oscillator model. The models are both derived from the time-dependent formulation of electronic spectroscopy. They assume that the transition dipoles do not depend on geometry and that the potential energy surfaces are locally quadratic; one model is second order in time and is called LQ2, and the other is third order in time and is called LQ3. These models are suitable for simulating the unresolved vibronic band shapes of electronic spectra that involve many vibrational modes. The models are straightforward and can be easily applied to simulate absorption spectra that are composed of many electronic transitions. As compared to carrying out molecular dynamics simulations, they require relatively few electronic structure calculations, and the additional cost for constructing the spectra is negligible. Therefore, the models are suitable for simulating the spectra of complex systems such as transition-metal complexes.

Low-Lying ππ* States of Heteroaromatic Molecules: A Challenge for Excited State Methods

The description of low-lying ππ* states of linear acenes by standard electronic structure methods is known to be challenging. Here, we broaden the framework of this problem by considering a set of fused heteroaromatic rings and demonstrate that standard electronic structure methods do not provide a balanced description of the two (typically) lowest singlet state (L_a and L_b) excitations. While the L_b state is highly sensitive to correlation effects, L_a suffers from the same drawbacks as charge transfer excitations. We show that the comparison between CIS/CIS(D) can serve as a diagnostic for detecting the two problematic excited states. Standard TD-DFT and even its spin-flip variant lead to inaccurate excitation energies and interstate gaps, with only a double hybrid functional performing somewhat better. The complication inherent to a balanced description of these states is so important that even CC2 and ADC(2) do not necessarily match the ADC(3) reference.