Quantifying the Performances of DFT for Predicting Vibrationally Resolved Optical Spectra: Asymmetric Fluoroborate Dyes as Working Examples

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.

One- and two-photon absorption spectra of organoboron complexes: vibronic and environmental effects

We synthesized a series of four parent aza-β-ketoiminate organoboron complexes and performed spectroscopic studies using both experimental and computational techniques. We studied how benzannulation influences the vibronic structure of the UV/Vis absorption bands with a focus on the bright lowest-energy π → π* electronic excitation. Theoretical simulations, accounting for inhomogeneous broadening effects using different embedding schemes, allowed gaining in-depth insights into the observed differences in band shapes induced by structural modifications. We observed huge variations in the distributions of vibronic transitions depending on the position of benzannulation. By and large, the harmonic approximation combined with the adiabatic hessian model delivers qualitatively correct band shapes for the one-photon absorption spectra, except in one case. We also assessed the importance of non-Condon effects (accounted for by the linear term in Herzberg-Teller expansion of the dipole moment) for S0 → S1 band shapes. It turned out that non-Condon contributions have no effect on the band shape in one-photon absorption spectra. In contrast, these effects significantly change the Franck-Condon band shapes of the two-photon absorption spectra. For one of the studied organoboron complexes we also performed a preliminary exploration of mechanical anharmonicity, resulting in an increase of the intensity of the 0-0 transition, which improves the agreement with the experimental data compared to the harmonic model.

Pitfall in simulations of vibronic TD-DFT spectra: diagnosis and assessment

In this Communication, we study the effect of spurious oscillations in the profiles of energy derivatives with respect to nuclear coordinates calculated with density functional approximations (DFAs) for formaldehyde, pyridine, and furan in their ground and electronic excited states. These spurious oscillations, which can only be removed using extensive integration grids that increase enormously the CPU cost of DFA calculations, are significant in the case of third- and fourth-order energy derivatives of the ground and excited states computed by M06-2X and ωB97X functionals. The errors in question propagate to anharmonic vibronic spectra computed under the Franck-Condon approximation, i.e., positions and intensities of vibronic transitions are affected to a large extent (shifts as significant as hundreds of cm-1 were observed). On the other hand, the LC-BLYP and CAM-B3LYP functionals show a much less pronounced effect due to spurious oscillations. Based on the results presented herein, we recommend either LC-BLYP or CAM-B3LYP with integration grids (250, 974) (or larger) for numerically stable simulations of vibronic spectra including anharmonic effects.

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%.

Metal complexes of benzimidazole-derived as potential anti-cancer agents: synthesis, characterization, combined experimental and computational studies

In this study, a series of 14 Cu (II), Zn (II), Ni (II) and Ag (I) complexes containing bis-benzimidazole derivatives were successfully designed and synthesized from 2-(1H-benzimidazole-2-yl)-phenol derivatives and corresponding metal salt solutions. The compound structures were identified by FT-IR, ¹H-NMR, powder X-ray diffraction and ESI-MS analyses, and the presence of the metal in the complexes was confirmed by ultraviolet-visible spectroscopy and ICP optical emission spectrometry. Electronic structure calculations were also carried out to describe the detailed structures in addition to the electronic absorption spectra of the ligands. The cytotoxic activity of the complexes was evaluated against three human cancer cell lines: lung (A549), breast (MDA-MB-231) and prostate (PC3) cancer cells. All complexes inhibited anti-proliferative cancer cells better than free ligands, especially Zn (II) and Ag (I) complexes, which are most sensitive to MDA-MB-231 cells. In addition, showing the growth inhibition of three cancer cell lines with IC₅₀ < 10.4 µM, complexes C₁ , C₃ and C₁₄ could be considered potential multi-targeted anti-cancer agents.

Metal complexes of benzimidazole-derived as potential anti-cancer agents: synthesis, characterization, combined experimental and computational studies

In this study, a series of 14 Cu (II), Zn (II), Ni (II) and Ag (I) complexes containing bis-benzimidazole derivatives were successfully designed and synthesized from 2-(1H-benzimidazole-2-yl)-phenol derivatives and corresponding metal salt solutions. The compound structures were identified by FT-IR, ¹H-NMR, powder X-ray diffraction and ESI-MS analyses, and the presence of the metal in the complexes was confirmed by ultraviolet-visible spectroscopy and ICP optical emission spectrometry. Electronic structure calculations were also carried out to describe the detailed structures in addition to the electronic absorption spectra of the ligands. The cytotoxic activity of the complexes was evaluated against three human cancer cell lines: lung (A549), breast (MDA-MB-231) and prostate (PC3) cancer cells. All complexes inhibited anti-proliferative cancer cells better than free ligands, especially Zn (II) and Ag (I) complexes, which are most sensitive to MDA-MB-231 cells. In addition, showing the growth inhibition of three cancer cell lines with IC₅₀ < 10.4 µM, complexes C₁ , C₃ and C₁₄ could be considered potential multi-targeted anti-cancer agents.

FB-REDA: fragment-based decomposition analysis of the reorganization energy for organic semiconductors

We present a fragment-based decomposition analysis tool (FB-REDA) for the reorganisation energy (λ). This tool delivers insights on how to rationally design low-λ organic semiconductors. The contribution of the fragment vibrational modes to the reorganization energy is exploited to identity the individual contributions of the molecular building blocks. The usefulness of the approach is demonstrated by offering three strategies to reduce the reorganization energy of a promising dopant-free hole transport material (TPA1PM, λ = 213 meV). A reduction of nearly 50% (TPD3PM, λ = 108 meV) is achieved. The proposed design principles are likely transferable to other organic semiconductors exploiting common molecular building blocks.

Intramolecular Hydrogen Bonds in Selected Aromatic Compounds: Recent Developments

A review of intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl Mannich bases, dipyrrins, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, and 4-Bora-3a,4a-diaza-s-indacene (BODIPY) dyes with tautomeric sensors as substituents is presented in this paper. Ortho-hydroxy Schiff and Mannich base derivatives are known as model molecules for analysing the properties of intramolecular hydrogen bonding. The compounds under discussion possess physicochemical features modulated by the presence of strong intramolecular hydrogen bonds. The equilibrium between intra- and inter-molecular hydrogen bonds in BODIPY is discussed. Therefore, the summary can serve as a knowledge compendium of the influence of the hydrogen bond on the molecular properties of aromatic compounds.

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.

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.

Combining the ensemble and Franck-Condon approaches for calculating spectral shapes of molecules in solution

The correct treatment of vibronic effects is vital for the modeling of absorption spectra of many solvated dyes. Vibronic spectra for small dyes in solution can be easily computed within the Franck-Condon approximation using an implicit solvent model. However, implicit solvent models neglect specific solute-solvent interactions on the electronic excited state. On the other hand, a straightforward way to account for solute-solvent interactions and temperature-dependent broadening is by computing vertical excitation energies obtained from an ensemble of solute-solvent conformations. Ensemble approaches usually do not account for vibronic transitions and thus often produce spectral shapes in poor agreement with experiment. We address these shortcomings by combining zero-temperature vibronic fine structure with vertical excitations computed for a room-temperature ensemble of solute-solvent configurations. In this combined approach, all temperature-dependent broadening is treated classically through the sampling of configurations and quantum mechanical vibronic contributions are included as a zero-temperature correction to each vertical transition. In our calculation of the vertical excitations, significant regions of the solvent environment are treated fully quantum mechanically to account for solute-solvent polarization and charge-transfer. For the Franck-Condon calculations, a small amount of frozen explicit solvent is considered in order to capture solvent effects on the vibronic shape function. We test the proposed method by comparing calculated and experimental absorption spectra of Nile red and the green fluorescent protein chromophore in polar and non-polar solvents. For systems with strong solute-solvent interactions, the combined approach yields significant improvements over the ensemble approach. For systems with weak to moderate solute-solvent interactions, both the high-energy vibronic tail and the width of the spectra are in excellent agreement with experiments.

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.