Benchmarking TD-DFT against Vibrationally Resolved Absorption Spectra at Room Temperature: 7-Aminocoumarins as Test Cases

Effect of Pyrrolidinedithiocarbamate Ligand on the Luminescence Properties of Heteroligand Samarium and Europium Complexes: Experimental and Theoretical Study

The photophysical properties of two isostructural heteroligand lanthanide complexes of general formula Ln(pdtc)3(phen) (pdtc = pyrrolidinedithiocarbamate anion, phen = 1,10-phenanthroline), Ln = Sm3+ (1), Eu3+ (2)) were studied in solid state and dichloromethane (DCM) solution. The two lanthanide complexes were investigated by experimental techniques for structural (single-crystal X-ray diffraction analysis of 1, powder XRD, TG-DTA) and spectroscopic [electron paramagnetic resonance (EPR), infrared (IR), ultraviolet-visible (UV-vis), photoluminescence (PL)] characterization. DFT/TDDFT/ωB97xD and multireference SA-CASSCF/NEVPT2 calculations with perturbative spin-orbit coupling corrections were applied to construct the Jablonski energy diagrams and to discuss the excited state energy transfer mechanism with competing excited state processes and possible sensitized mechanism of metal-centered emission. The first excited state (S1) involved in the excited state energy transfer L(antenna)-to-Ln was predicted to have interligand (pdtc-to-phen) charge transfer character in contrast to the previously predicted ligand-to-metal charge transfer character. The theoretical consideration showed similar relaxation paths and luminescence quenching channels and appropriate Donor*(phen)-Acceptor*(Ln3+) energy gap for 1 and 2. The experimental measurements in the solid state, however, showed efficient luminescence and good ability to convert UV to visible light only for the Sm(pdtc)3(phen) complex. The minor emission of 2 was explained by partial reduction of Eu3+, confirmed by EPR and calculated electron density distribution data.

First-Principles Modeling of the Absorption Spectrum of Crystal Violet in Solution: The Importance of Environmentally Driven Symmetry Breaking

Theoretical spectroscopy plays a crucial role in understanding the properties of the materials and molecules. One of the most promising methods for computing optical spectra of chromophores embedded in complex environments from the first principles is the cumulant approach, where both (generally anharmonic) vibrational degrees of freedom and environmental interactions are explicitly accounted for. In this work, we verify the capabilities of the cumulant approach in describing the effect of complex environmental interactions on linear absorption spectra by studying Crystal Violet (CV) in different solvents. The experimental absorption spectrum of CV strongly depends on the nature of the solvent, indicating strong coupling to the condensed-phase environment. We demonstrate that these changes in absorption line shape are driven by an increased splitting between absorption bands of two low-lying excited states that is caused by a breaking of the D3 symmetry of the molecule and that in polar solvents, this symmetry breaking is mainly driven by electrostatic interactions with the condensed-phase environment rather than distortion of the structure of the molecule, in contrast with conclusions reached in a number of previous studies. Our results reveal the importance of explicitly including a counterion in the calculations in nonpolar solvents due to electrostatic interactions between CV and the ion. In polar solvents, these interactions are strongly reduced due to solvent screening effects, thus minimizing the symmetry breaking. Computed spectra in methanol are found to be in reasonable agreement with the experiment, demonstrating the strengths of the outlined approach in modeling strong environmental interactions.

The role of microenvironments on computed vibrationally-resolved emission spectra: The case of oxazines

Oxazine dyes act as reporters of their near environment by the response of their fluorescence spectra. At the same time, their fluorescence spectra exhibit a pronounced vibrational progression. In this work, we computationally investigate the impact of near-environment models consisting of aggregated water as well as betaine molecules on the vibrational profile of fluorescence spectra of different oxazine derivatives. For aggregated betaine and a water molecule located above the plane of the dyes, we observe a distinct modification of the vibrational profile, which is more pronounced than the effect of a continuum description of a solvent environment. Our analysis shows that this effect cannot be explained by a pure change in the electronic structure, but that also vibrational degrees of freedom of the environment can be decisive for the vibrational profile and should, hence, not generally be neglected.

Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels

We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.

Octopus-inspired deception and signaling systems from an exceptionally-stable acene variant

Multifunctional platforms that can dynamically modulate their color and appearance have attracted attention for applications as varied as displays, signaling, camouflage, anti-counterfeiting, sensing, biomedical imaging, energy conservation, and robotics. Within this context, the development of camouflage systems with tunable spectroscopic and fluorescent properties that span the ultraviolet, visible, and near-infrared spectral regions has remained exceedingly challenging because of frequently competing materials and device design requirements. Herein, we draw inspiration from the unique blue rings of the Hapalochlaena lunulata octopus for the development of deception and signaling systems that resolve these critical challenges. As the active material, our actuator-type systems incorporate a readily-prepared and easily-processable nonacene-like molecule with an ambient-atmosphere stability that exceeds the state-of-the-art for comparable acenes by orders of magnitude. Devices from this active material feature a powerful and unique combination of advantages, including straightforward benchtop fabrication, competitive baseline performance metrics, robustness during cycling with the capacity for autonomous self-repair, and multiple dynamic multispectral operating modes. When considered together, the described exciting discoveries point to new scientific and technological opportunities in the areas of functional organic materials, reconfigurable soft actuators, and adaptive photonic systems.

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.

Cost-Effective Simulations of Vibrationally-Resolved Absorption Spectra of Fluorophores with Machine-Learning-Based Inhomogeneous Broadening

The results of electronic and vibrational structure simulations are an invaluable support for interpreting experimental absorption/emission spectra, which stimulates the development of reliable and cost-effective computational protocols. In this work, we contribute to these efforts and propose an efficient first-principle protocol for simulating vibrationally-resolved absorption spectra, including nonempirical estimations of the inhomogeneous broadening. To this end, we analyze three key aspects: (i) a metric-based selection of density functional approximation (DFA) so to benefit from the computational efficiency of time-dependent density function theory (TD-DFT) while safeguarding the accuracy of the vibrationally-resolved spectra, (ii) an assessment of two vibrational structure schemes (vertical gradient and adiabatic Hessian) to compute the Franck-Condon factors, and (iii) the use of machine learning to speed up nonempirical estimations of the inhomogeneous broadening. In more detail, we predict the absorption band shapes for a set of 20 medium-sized fluorescent dyes, focusing on the bright ππ^★ S₀ → S₁ transition and using experimental results as references. We demonstrate that, for the studied 20-dye set which includes structures with large structural variability, the preselection of DFAs based on an easily accessible metric ensures accurate band shapes with respect to the reference approach and that range-separated functionals show the best performance when combined with the vertical gradient model. As far as band widths are concerned, we propose a new machine-learning-based approach for determining the inhomogeneous broadening induced by the solvent microenvironment. This approach is shown to be very robust offering inhomogeneous broadenings with errors as small as 2 cm^-1 with respect to genuine electronic-structure calculations, with a total CPU time reduced by 98%.

Modeling Luminescence Spectrum of BaZrO3:Ti Including Vibronic Coupling from First Principles Calculations

Herein, we present a methodology based on constrained density functional theory and vibrational mode computations to simulate and interpret the luminescence spectra of periodic solids. A multi-dimension harmonic model is used to combine electronic and vibrational contributions into an overall vibrationally resolved emission spectrum. We applied it to Ti-doped BaZrO₃ to accurately reproduce its blue luminescence and unambiguously assign the observed luminescence to a Ti³⁺ + O^- → Ti⁴⁺ + O^2- charge transfer.

Estimation of Excited-State Geometries of Benzene and Fluorobenzene through Vibronic Analyses of Absorption Spectra

The parameters used in theoretical modeling of vibrational patterns within Franck-Condon (FC) approximation can be adjusted to match the vibrationally well-resolved experimental absorption spectrum of molecules. These simulation parameters can then be used to reveal the structural changes occurring between the initial and final states assuming the harmonic oscillator approximation holds for both states. Such a theoretical approach has been applied to benzene and fluorobenzene to disclose the first excited-state geometries of both compounds. The carbon-carbon bond length of benzene in the ¹B_2u state has been calculated as 1.430 Å, which is in very good agreement with the experimental bond length of 1.432 Å. The FC spectral fit method has been exploited to reveal the ¹B₂ state of fluorobenzene as well. Commonly employed density functional theory (DFT) and time-dependent DFT methods have been used to calculate the ground- and excited-state geometries of both compounds, respectively. The comparison of geometrical parameters and vibrational frequencies at the relevant states shows that frequently used hybrid functionals perform quite well in the ground state, whereas their performances drop considerably while predicting the excited-state properties. Among the hybrid functionals studied, TD-B3LYP with 6-31+G(d) basis set can be chosen to calculate the excited-state properties of molecules, albeit with much less anticipation of accuracy from the performance that B3LYP usually shows at the ground state.

Synthesis and characterization of polymethine dyes carrying thiobarbituric and carboxylic acid moieties

Polymethine dyes are prepared using a convenient synthesis and characterized by physicochemical and computational methods.

How to Predict Excited State Geometry by Using Empirical Parameters Obtained from Franck-Condon Analysis of Optical Spectrum

Excited state geometries of molecules can be calculated with highly reliable wavefunction schemes. Most of such schemes, however, are applicable to small molecules and can hardly be viewed as error-free for excited state geometries. In this study, a theoretical approach is presented in which the excited state geometries of molecules can be predicted by using vibrationally resolved experimental absorption spectrum in combination with the theoretical modelling of vibrational pattern based on Franck-Condon approximation. Huang-Rhys factors have been empirically determined and used as input for revealing the structural changes occurring between the ground and the excited state geometries upon photoexcitation. Naphthalene molecule has been chosen as a test case to show the robustness of the proposed theoretical approach. Predicted 1B_2u excited state geometry of the naphthalene has similar but slightly different bond length alternation pattern when compared with the geometries calculated with CIS, B3LYP, and CC2 methods. Excited state geometries of perylene and pyrene molecules are also determined with the presented theoretical approach. This powerful method can be applied to other molecules and specifically to relatively large molecules rather easily as long as vibrationally resolved experimental spectra are available to use.

Assessment of long-range corrected density functional theory on the absorption and vibrationally resolved fluorescence spectrum of carbon nanobelts

Time-dependent (TD) density functional theory (DFT) and Franck-Condon Hertzberg-Teller (FCHT) calculations of various DFT functionals [B3LYP, CAM-B3LYP, ωB97XD, and optimally tuned (OT) long-range corrected (LC)-BLYP] were performed to examine how well DFT functionals can predict the experimental absorption and fluorescence spectra of a 12-carbon nanobelt (CNB). OT-LC-BLYP reproduced the experimental absorption spectrum well in terms of the peak position and intensity in the case of using a basis set with a diffuse function, such as 6-31+G(d,p) and 6-311+G(d,p), whereas B3LYP showed a red-shift in the peak positions and CAM-B3LYP and ωB97XD, which have a long-range HF exchange, showed blue shifts. Regarding the fluorescence spectrum calculations with FCHT using 6-311+G(d,p), the OT-LC-BLYP reproduced both the peak intensities and positions closest to the experimental spectrum. B3LYP, however, showed red-shifted peaks, and ωB97XD showed blue-shifted peaks. CAM-B3LYP provided less blue-shifted peaks, but the relative peak intensities mismatched the experimental ones. Furthermore, calculations of the absorption and vibrationally resolved fluorescence spectra of 16-CNB and 24-CNB using OT-LC-BLYP/6-311+G(d,p) showed absorption and fluorescence spectra close to the experimental spectra with high accuracy. Moreover, the application of a polarizable continuum model (dichloromethane) produced a red shift in the peak positions of the absorption spectrum with increasing intensity but an increase in the peak intensities of the fluorescence calculations without shifting the peak position.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

Theoretical study of α, β unsaturated carbonyl thiophene derivatives to investigate optoelectronic properties toward organic photovoltaics

Herein theoretical study, we designed sixteen conjugated arylated α, β unsaturated carbonyl thiophene based compounds by using density functional theory (DFT) and time-dependent (TD) density functional theory at modified Perdue Wang density functional MPW1PW91 functional with 6-31G (d, p) basis set. Ground and excited state geometries, electronic and photophysical characteristics of designed molecules are evaluated by assuming the electron-donating and electron-withdrawing effects of the substituents that are attached to these newly designed molecules. Furthermore, calculation of vibrational spectra, time-dependent effect, isotopic substitution effect and force constant along with thermodynamic quantities are also carried out by using MOPAC (Molecular Orbital Package) with strong implementation of semi empirical Hamiltonians. The results reveal that our designed molecules can be a good candidates for electroluminescent and optoelectronic devices for further fabrication of solar cell devices.

Energy transfer mechanism in luminescence Eu(III) and Tb(III) complexes of coumarin-3-carboxylic acid: A theoretical study

Excited state energy level diagrams of coumarin-3-carboxylic acid (HCCA) chromophore, Eu(CCA)Cl₂(H₂O)₂ (1), Eu(CCA)₂Cl(H₂O)₂ (2)_, Eu(CCA)₃(H₂O)₃ (3), Tb(CCA)₂Cl(H₂O) (4) and Tb(CCA)₂(NO₃)(H₂O) (5) in gas phase and polar solution have been calculated by means of DFT/TDDFT/ωB97XD methods. Based on these results, the ability of CCA to sensitize Eu(III) and Tb(III) luminescence has been examined. The competitive excited state processes in the complexes - fluorescence, intersystem crossing (ISC) and phosphorescence, were analyzed depending on the environment, number of the ligands, Ln(III) ion type (Eu and Tb) and counteranion (Cl^- and NO₃^-). It has been found that the environment altered the S₁ state energy, oscillator strength, fluorescence lifetime as well as the S₁ character - polar solution stabilized the S₁(ππ*) state, whereas non-polar solution (gas phase, solid state) stabilized the S₁(nπ*) state. The S₁(nπ*) state was decisive for the efficient energy transfer as it suppressed the S₁ emission of CCA and favored ISC or direct transfer to the emitting levels of Eu(III). The HCCA triplet (T₁) state minimum energy (~2.7, ~2.6^ZPE eV) and (ππ*) character were retained in Eu/Tb-CCA complexes regardless of the environment. The energy gap between the higher energy T₁ donor state and the acceptor levels ⁵D₁ of Eu(III) (~0.5 eV) and ⁵D₄ of Tb(III) (~0.1 eV) provided optimal resonance conditions for effective energy transfer for Eu(III), but less probability for Tb(III). The nonradiative energy (CCA → Eu(III)) transfer rates and quantum luminescence yield for 2 and 3 were calculated by a strategy combining DFT geometries, INDO/S excitation energies and calculated Judd-Ofelt parameters. The excitation channel T₁ → ⁵D₀ through an exchange mechanism was predicted as the most probable one to populate the main emissive Eu-centered state in complexes 2 and 3. The more efficient luminescence of 3 than that of 2 was discussed and explained.

Triangulenium dyes: the comprehensive photo-absorption and emission story of a versatile family of chromophores

The triangulenium dyes constitute a family of versatile chromophores whose impressive photo-absorption and emission properties are currently highlighted in numerous novel experimental applications. In this investigation, we provide a comprehensive TDDFT characterization of their spectroscopic properties elucidating the origin of their large and complex absorption and emission vibronic spectra spread over the (whole) visible region. More precisely, by benchmarking the performance of 10 commonly-used exchange-correlation density functionals belonging to different classes of approximation, we develop and validate a computational protocol allowing the accurate modeling of both the position and optical line-shape of their vibrationally-resolved absorption and emission band structures. We find that semilocal approximations provide the best estimate of the structure of the vibronic spectra, however they spuriously and strongly underestimate their position. We finally show that global-hybrid density functionals mixing between 20 and 30% of exact-like exchange are an excellent compromise to get a satisfactory estimate of both of these properties.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules

Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.

Unspecified verticality of Franck–Condon transitions, absorption and emission spectra of cyanine dyes, and a classically inspired approximation

The computed vertical energy, E_v,a/f, from the equilibrium geometry of the initial electronic state is frequently considered as representative of the experimental excitation/emission energy, E_abs/fl = hc/λ_max. Application of the quantum mechanical version of the Franck–Condon principle does not involve precise specification of nuclear positions before, after, or during an electronic transition. Moreover, the duration of an electronic transition is not experimentally accessible in spectra with resolved vibrational structure. It is shown that computed vibronic spectra based on TDDFT methods and application of quantum mechanical FC analysis predict E_abs = hc/λ_max with a 10-fold improvement in accuracy compared to E_v,a for nine cyanine dyes. It is argued that part of the reason for accuracy when this FC analysis is compared to experiment as opposed to E_v,a/f is the unspecified verticality of transitions in the context of the quantum version of the FC principle. Classical FC transitions that preserve nuclear kinetic energy before and after an electronic transition were previously found to occur at a weighted average of final and initial electronic state molecular geometries known as the r-centroid. Inspired by this approach a qualitative method using computed vertical and adiabatic energies and the harmonic approximation is developed and applied yielding a 5-fold improvement in accuracy compared to E_v,a. This improvement results from the dominance of low frequency vibronic transitions in the cyanine dye major band. The model gives insight into the nature of the redshift when qPCR dye EvaGreen is complexed to λDNA and is applicable to the low frequency band of similar non cyanine dyes such as curcumin. It is found that the computed vibronic cyanine dye spectra from time-dependent FC analysis at 0 K and 298 K show decreased intensity at higher temperature suggestive of increased intensity with restricted motion shown when cyanine dyes are used in biomedical imaging. A 2-layer ONIOM model of the DNA minor groove indicates restricted motion of the TC-1 dye excited state in this setting indicative of enhanced fluorescence.

Insight into cyanine dye λ_max from quantum and classical FC principle; high accuracy with classically inspired approach.

Photochemical mechanism of DEACM uncaging: a combined time-resolved spectroscopic and computational study

Combined spectroscopic and computational studies elucidate excited-state photocleavage in DEACM cages, explaining vastly different time scales for different leaving groups.

Simulation of Vacuum Ultraviolet Absorption Spectra: Paraffin, Isoparaffin, Olefin, Naphthene, and Aromatic Hydrocarbon Class Compounds

The advent of a new vacuum ultraviolet (VUV) spectroscopic absorption detector for gas chromatography has enabled applications in many areas. Theoretical simulations of VUV spectra using computational chemistry can aid the new technique in situations where experimental spectra are unavailable. In this study, VUV spectral simulations of paraffin, isoparaffin, olefin, naphthene, and aromatic (PIONA) compounds using time-dependent density functional theory (TDDFT) methods were investigated. Important factors for the simulations, such as functionals/basis sets and formalism of oscillator strength calculations, were examined and parameters for future PIONA compound simulations were obtained by fitting computational results to experimental spectra. The simulations produced satisfactory correlations between experimental observations and theoretical calculations, and enabled potential analysis applications for complex higher distillate fuels, such as diesel fuel. Further improvement of the methods was proposed.

Aggregation Effects on Pigment Coatings: Pigment Red 179 as a Case Study

Here, we have studied, with a combined experimental and computational approach, the effect of the crystal environment and aggregation on the electronic properties of Pigment Red 179, which affect both its color and optical energy gap. Spectra acquired in the near-infrared and visible range of energies suggest that this molecule is indeed a "cool" dye, which can be employed as a red pigment that provides effective color coverage to different substrates without contributing to their heating during light irradiation. Spectra acquired on different polymer mixtures at different pigment concentrations (i.e., 2.5-10 wt %) suggest that absorption features depend on chromophoric arrangements promoted by the strong intermolecular π-π interactions. Calculations, performed at the time-dependent density functional theory level, allowed to both attribute the nature of the electronic transitions causing the observed spectra involved and understand the effect of the environment. Indeed, the visible spectra of the pigment is dominated by two localized transitions, with negligible charge transfer for both a dye monomer and dimer either in vacuum or acetonitrile solution. Instead, models including the crystal environment of the pigment show the presence of a high-wavelength S1 ← S0 charge transfer transition between two adjacent molecules, in quantitative agreement with the experimental absorption energy of the crystal pigment.

On the simulation of vibrationally resolved electronic spectra of medium-size molecules: the case of styryl substituted BODIPYs

BODIPY dyes are used in a variety of applications because of their peculiar spectroscopic and photo-physical properties that vary depending on the stereochemistry of the functional groups attached to the boron-dipyrromethene core structure. In this work, we have applied several computational methods, adapted for semi-rigid molecules based on the Franck-Condon principle, for the study of the optical properties of BODIPY systems and for the understanding of the influence of functional groups on their spectroscopic features. We have analyzed the electronic spectra of two styryl substituted BODIPY molecules of technological interest, properly taking into account the vibronic contribution. For comparison with recently recorded experimental data in methanol, the vibrationally resolved electronic spectra of these systems were computed using both Time-Independent (TI) and Time-Dependent (TD) formalisms. The first step toward the analysis of optical properties of the styryl modified BODIPYs was a benchmark of several density functionals, to select the most appropriate one. We have found that all benchmarked functionals provide good results in terms of band shape but some of them show strong discrepancies in terms of band position. Beyond the issue of the electronic structure calculation method, different levels of sophistication can be adopted for the calculation of vibronic transitions. In this study, the effect of mode couplings and the influence of the Herzberg-Teller terms on the theoretical spectra has been investigated. It has been found that all levels of theory considered give reproducible results for the investigated systems: band positions and shapes are similar at all levels and little improvements have been found in terms of band shape with the inclusion of Herzberg-Teller effect. Inclusion of temperature effects proved to be challenging due to the important impact of large amplitude motions. Better agreement can be achieved by adopting a suitable set of coordinates coupled with a reduced-dimensionality scheme.

First-principles investigation of the double ESIPT process in a thiophene-based dye

Explanation of the experimental triple emission with theoretical tools requires advanced solvent models.

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.

Analytical Evaluation for Calculation of Two-Center Franck–Condon Factor and Matrix Elements

The Franck–Condon (FC) factor is defined as squares of the Franck–Condon (FC) overlap integral and represents one of the principle fundamental factors of molecular physics. The FC factor is used to determine the transition probabilities in different vibrational levels of the two electronic states and the spectral line intensities of diatomic and polyatomic molecules. In this study, new analytical formulas were derived to calculate Franck–Condon integral (FCI) of harmonic oscillators and matrix elements (xη, e−2cx, and e−cx2) including simple finite summations of binomial coefficients. These formulas are valid for arbitrary values. The results of formulas are in agreement with the results in the literature.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Using non-empirically tuned range-separated functionals with simulated emission bands to model fluorescence lifetimes

Fluorescence lifetimes were evaluated using TD-DFT under different approximations for the emitting molecule and various exchange-correlation functionals, such as B3LYP, BMK, CAM-B3LYP, LC-BLYP, M06, M06-2X, M11, PBE0, ωB97, ωB97X, LC-BLYP*, and ωB97X* where the range-separation parameters in the last two functionals were tuned in a non-empirical fashion. Changes in the optimised molecular geometries between the ground and electronically excited states were found to affect the quality of the calculated lifetimes significantly, while the inclusion of vibronic features led to further improvements over the assumption of a vertical electronic transition. The LC-BLYP* functional was found to return the most accurate fluorescence lifetimes with unsigned errors that are mostly within 1.5 ns of experimental values.

Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules

We present an investigation of the excited-state structural parameters determined for a large set of small compounds with the dual goals of defining reference values for further works and assessing the quality of the geometries obtained with relatively cheap computational approaches. In the first stage, we compare the excited-state geometries obtained with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large atomic basis sets. It is found that CASPT2 and CC3 results are generally in very good agreement with one another (typical differences of ca. 3 × 10^–3 Å) when all electrons are correlated and when the aug-cc-pVTZ atomic basis set is employed with both methods. In a second stage, a statistical analysis reveals that, on the one hand, the excited-state (ES) bond lengths are much more sensitive to the selected level of theory than their ground-state (GS) counterparts and, on the other hand, that CCSDR(3) is probably the most cost-effective method delivering accurate structures. Indeed, CCSD tends to provide too compact multiple bond lengths on an almost systematic basis, whereas both CC2 and ADC(2) tend to exaggerate these bond distances, with more erratic error patterns, especially for the latter method. The deviations are particularly marked for the polarized CO and CN bonds, as well as for the puckering angle in formaldehyde homologues. In the last part of this contribution, we provide a series of CCSDR(3) GS and ES geometries of medium-sized molecules to be used as references in further investigations.

UV characterization and photodegradation mechanism of the fungicide chlorothalonil in the presence and absence of oxygen

An experimental and theoretical study of the UV spectrum of chlorothalonil (CT) was carried out and the vibrationally resolved HOMO→LUMO transition is presented for the first time. The fluorescence spectrum has also been recorded. Furthermore, preparative photolysis allowed a detailed study of the photoproducts formed with recognition of different isomers. In the presence of oxygen only the first reductive dechlorination-decyanation occurred, while in its absence a successive dechlorination-decyanation takes place.

Nonempirical Simulations of Inhomogeneous Broadening of Electronic Transitions in Solution: Predicting Band Shapes in One- and Two-Photon Absorption Spectra of Chalcones

We have examined several approaches relying on the Polarizable Embedding (PE) scheme to predict optical band shapes for two chalcone molecules in methanol solution. The PE-TDDFT and PERI-CC2 methods were combined with molecular dynamics simulations, where the solute geometry was kept either as rigid, flexible or partly-flexible (restrained) body. The first approach, termed RBMD-PE-TDDFT, was employed to estimate the inhomogeneous broadening for subsequent convolution with the vibrationally-resolved spectra of the molecule in solution determined quantum-mechanically (QM). As demonstrated, the RBMD-PE-TDDFT/QM-PCM approach delivers accurate band widths, also reproducing their correct asymmetric shapes. Further refinement can be obtained by the estimation of the inhomogeneous broadening using the RBMD-PERI-CC2 method. On the other hand, the remaining two approaches (FBMD-PE-TDDFT and ResBMD-PE-TDDFT), which lack quantum-mechanical treatment of molecular vibrations, lead to underestimated band widths. In this study, we also proposed a simple strategy regarding the rapid selection of the exchange-correlation functional for the simulations of vibrationally-resolved one- and two-photon absorption spectra based on two easy-to-compute metrics.

Theoretical simulations for vibrationally-resolved absorption spectra of naphthalenediimide cyclophane derivatives

In the present work we systematically investigated the vibrationally-resolved absorption spectra of three core substituted naphthalenediimide cyclophane derivatives. It has been performed by time-dependent density functional theory calculations using three different exchange-correlation functionals, including the conventional B3LYP functional and two long-range corrected functionals: CAM-B3LYP and ωB97XD. The solvent effects were also considered with the polarizable continuum model. Calculation results showed that long range corrections are needed to correctly describe the optical properties of the three molecules because of the strong charge transfer characteristic of the excited states. The core substitution induced red shift to the first absorption band is nicely explained by the theoretical calculations. It is found that this band mainly involves the transitions within the core substituted naphthalenediimide chromophore. The high energy absorption band, on the other hand, is generated mainly from the un-substituted chromophore. These characters result in different substitution dependence for those two main absorption bands. Furthermore, the simulated vibrational profiles of the first two absorption bands nicely reproduce the observed vibrational features in the measured spectra. The accuracy of the calculated results from different functionals and basis sets has been discussed.

Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach

The electronic absorption spectra of pyridine and nicotine in aqueous solution have been computed using a multistep approach. The computational protocol consists in studying the solute solvation with accurate molecular dynamics simulations, characterizing the hydrogen bond interactions, and calculating electronic transitions for a series of configurations extracted from the molecular dynamics trajectories with a polarizable QM/MM scheme based on the fluctuating charge model. Molecular dynamics simulations and electronic transition calculations have been performed on both pyridine and nicotine. Furthermore, the contributions of solute vibrational effect on electronic absorption spectra have been taken into account in the so called vertical gradient approximation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.