Electronically excited states in poly(p-phenylenevinylene): vertical excitations and torsional potentials from high-level ab initio calculations

Contrasting the excited state properties of different conformers of trans- and cis-2,2'-bipyridine oligomers in the gas phase

In this article, we present conformation-dependent photophysical and excited state properties of trans- and cis- BPY oligomers. Oligomers up to tetramers for three conformers, namely, o-, m-, and p-, are constructed and optimized at the B3LYP-D3/def2-SVPD level. The photophysical and excited state properties are interpreted in terms of UV and CD spectra at the RI-ADC(2)/def2-TZVPD level. The UV spectra of oligomers of the m-conformer show high-intensity and red-shifted UV bands compared to o- and p-oligomers. The CD spectra of p-oligomers show intense CD bands compared to o- and p-oligomers in the case of trans-structures. In contrast, oligomers of each conformer of cis-structures show high-intensity CD bands. The excited states of (BPY)2 and (BPY)4 are also characterized by analysis of one-electron transition density matrix considering three descriptors: ωCT, dexc, and PRNTO. The ωCT values of dimers are in the range of 0.06-0.32, which indicates the excited states are mainly LE states, whereas, for (BPY)4, the ωCT values range from 0.17 to 0.53, indicating the possibility of partial CT in the excited states. These observations are also explained using the NTOs and e-h correlation plots.

Conformational Effect on the Excitonic States of 2-Phenylpyridine Oligomers: Ab Initio Studies and Analysis

In this work, we report the effect of different conformations of 2-phenylpyridine oligomers ((PhPy)n=1-5) on the excited state properties from the results obtained at the RI-ADC(2)/def2-TZVP level. Three different conformers, namely, A, B, and C, are considered for each oligomer. All the oligomers of conformer A have linear-type structures, whereas conformers B and C form helical structures at n = 5 and n = 3, respectively. The differences in the geometries of the three conformers are reflected in the UV and CD spectra. The UV spectra of conformer A show high-intensity peaks compared to the conformers B and C, for each oligomer. While the helical oligomers of conformers B and C show high-intensity CD bands, the intensities of CD bands for all of the oligomers of conformer A are weaker. Analysis of the properties of the first five excited states in (PhPy)5 is carried out using three descriptors, and the results reveal that these are partially charge transfer states.

Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models

We report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.

Modeling and Characterization of Exciplexes in Photoredox CO2 Reduction: Insights from Quantum Chemistry and Fluorescence Spectroscopy

Interactions between excited-state arenes and amines can lead to the formation of structures with a distinct emission behavior. These excited-state complexes or exciplexes can reduce the ability of the arene to participate in other reactions, such as CO₂ reduction, or increase the likelihood of degradation via Birch reduction. Exciplex geometries are necessary to understand photophysical behavior and probe degradation pathways but are challenging to calculate. We establish a detailed computational protocol for calculation, verification, and characterization of exciplexes. Using fluorescence spectroscopy, we first demonstrate the formation of exciplexes between excited-state oligo-(p-phenylene) (OPP), shown to successfully carry out CO₂ reduction, and triethylamine. Time-dependent density functional theory is employed to optimize the geometries of these exciplexes, which are validated by comparing both emission energies and their solvatochromism with the experiment. Excited-state energy decomposition analysis confirms the predominant role played by charge transfer interactions in the red shift of emissions relative to the isolated excited-state OPP*. We find that although the exciplex emission frequency depends strongly on solvent dielectric, the extent of charge separation in an exciplex does not. Our results also suggest that the formation of solvent-separated ionic radical states upon complete electron transfer competes with exciplex formation in higher-dielectric solvents, thereby leading to reduced exciplex emission intensities in fluorescence experiments.

Effects of Heterocyclic Ring Fusion and Chain Elongation on Chiroptical Properties of Polyaza[9]helicene: A Computational Study

In the present work, the effect of lateral and helical extensions on the physical and chiroptical properties of azahelicenes is reported. Starting with the experimentally reported polyaza[9]helicene (9Ha), three derivatives, two with laterally fused electron-withdrawing rings and the third with larger helical length, are designed. For the excited-state properties such as UV-vis and CD spectra, performances of different DFT functionals are evaluated by comparing the energies and characters of the excited states against the ADC(2) results. CPL properties are calculated at DFT level. Among the three designed systems, pyrazine-based 9HaP shows an improved g_CPL value compared to that for parent 9Ha. However, quinoxaline-based 9HaQ is found to be the worst CPL emitter with the lowest dissymmetry factor. The helically extended derivative, 11Ha, shows good CPL results, but g_CPL remains smaller than that for the parent system. The CPL results are analyzed in terms of electric dipole transition moment (EDTM) and magnetic dipole transition moment (MDTM) vectors, and angles between these two vectors.

Functional group introduction and aromatic unit variation in a set of π-conjugated macrocycles: revealing the central role of local and global aromaticity

π-Conjugated macrocycles are molecules with unique properties that are increasingly exploited for applications and the question of whether they can sustain global aromatic or antiaromatic ring currents is particularly intriguing. However, there are only a small number of experimental studies that investigate how the properties of π-conjugated macrocycles evolve with systematic structural changes. Here, we present such a systematic experimental study of a set of [2.2.2.2]cyclophanetetraenes, all with formally Hückel antiaromatic ground states, and combine it with an in-depth computational analysis. The study reveals the central role of local and global aromaticity for rationalizing the observed optoelectronic properties, ranging from extremely large Stokes shifts of up to 1.6 eV to reversible fourfold reduction, a highly useful feature for charge storage/accumulation applications. A recently developed method for the visualization of chemical shielding tensors (VIST) is applied to provide unique insight into local and global ring currents occurring in different planes along the macrocycle. Conformational changes as a result of the structural variations can further explain some of the observations. The study contributes to the development of structure-property relationships and molecular design guidelines and will help to understand, rationalize, and predict the properties of other π-conjugated macrocycles. It will also assist in the design of macrocycle-based supramolecular elements with defined properties.

Predicting phosphorescence energies and inferring wavefunction localization with machine learning

Phosphorescence is commonly utilized for applications including light-emitting diodes and photovoltaics. Machine learning (ML) approaches trained on ab initio datasets of singlet-triplet energy gaps may expedite the discovery of phosphorescent compounds with the desired emission energies. However, we show that standard ML approaches for modeling potential energy surfaces inaccurately predict singlet-triplet energy gaps due to the failure to account for spatial localities of spin transitions. To solve this, we introduce localization layers in a neural network model that weight atomic contributions to the energy gap, thereby allowing the model to isolate the most determinative chemical environments. Trained on the singlet-triplet energy gaps of organic molecules, we apply our method to an out-of-sample test set of large phosphorescent compounds and demonstrate the substantial improvement that localization layers have on predicting their phosphorescence energies. Remarkably, the inferred localization weights have a strong relationship with the ab initio spin density of the singlet-triplet transition, and thus infer localities of the molecule that determine the spin transition, despite the fact that no direct electronic information was provided during training. The use of localization layers is expected to improve the modeling of many localized, non-extensive phenomena and could be implemented in any atom-centered neural network model.

Assessing the Performance of DFT Functionals for Excited-State Properties of Pyridine-Thiophene Oligomers

In this article, we have examined the accuracy of various density functional theory (DFT) functionals to reproduce the absorption and CD spectra of pyridine-thiophene oligomers. The performance of different levels of approximations in DFT functionals is discussed with reference to the ADC(2) results. Starting from a linear system, like monomer, calculations are carried out at ADC(2) and DFT levels till a helical system, like pentamer, is formed. For vertical excitation energies, results obtained with functionals, like CAM-B3LYP, ωB97XD, and M06-2X, are closer to the ADC(2) results. However, analysis of excited-state properties shows that the state ordering patterns or results regarding natural transition orbitals from these DFT functionals sometimes differ from the ADC(2) results. Global hybrid functionals like B3LYP and PBE0 produce excitation energies which are far away from the ADC(2) benchmark results. Similarly, pure functionals and their long-range corrected versions produce either redshifted or blueshifted energies. For the CD spectra, the above three mentioned functionals, CAM-B3LYP, ωB97XD, and M06-2X, again produce spectra closer to the benchmark spectra.

Assessing the effects of increasing conjugation length on exciton diffusion: from small molecules to the polymeric limit

Organic solar cells (OSC) generally contain long-chain π-conjugated polymers as donor materials, but, more recently, small-molecule donors have also attracted considerable attention. The nature of these compounds is of crucial importance concerning the various processes that determine device performance, among which singlet exciton diffusion is one of the most relevant. The efficiency of the diffusion mechanism depends on several aspects, from system morphology to electronic structure properties, which vary importantly with molecular size. In this work, we investigated the effects of conjugation length on the exciton diffusion length through electronic structure calculations and an exciton diffusion model. By applying extrapolation procedures to thiophene and phenylene vinylene oligomer series, we investigate their electronic and optical properties from the small-molecule point of view to the polymeric limit. Several properties are calculated as a function of oligomer size, including transition energies, absorption and emission spectra, reorganization energies, exciton coupling and Förster radii. Finally, an exciton diffusion model is used to estimate diffusion lengths as a function of oligomer size and for the polymeric limit showing agreement with experimental data. Results also show that longer conjugation lengths correlate with longer exciton diffusion lengths in spite of also being associated with shorter exciton lifetimes.

Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature

We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.

First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature

We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σ_s ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10^-2 cm² s^-1 to D ∼ 1 × 10^-3 cm² s^-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.

Computational Analysis of Electron Transfer Kinetics for CO2 Reduction with Organic Photoredox Catalysts

We present a fundamental description of the electron transfer (ET) step from substituted oligo(p-phenylene) (OPP) radical anions to CO₂, with the larger goal of assessing the viability of underexplored, organic photoredox routes for utilization of anthropogenic CO₂. This work varies the electrophilicity of para-substituents to OPP and probes the dependence of rate coefficients and interfragment interactions on the substituent Hammett parameter, σ_p, using constrained density functional theory (CDFT) and energy decomposition analysis (EDA). Large electronic couplings across substituents indicates an adiabatic electron transfer process for reactants at contact. As one might intuitively expect, free energy changes dominate trends in ET rate coefficients in most cases, and rates increase with substituent electron-donating ability. However, we observe an unexpected dip in rate coefficients for the most electron-donating groups, due to the combined impact of flattening free energies and a steep increase in reorganization energies. Our analysis shows that, with decreasing σ_p, flattening OPP LUMO levels lower the marginal increase in free energy. EDA reveals trends in electrostatics and charge transfer interactions between the catalyst and substrate fragments that influence free energy changes across substituents. Reorganization energies do not exhibit a direct dependence on σ_p and are largely similar across systems, with the exception of substituents containing lone pairs of electrons that exhibit significant deformation upon electron transfer. Our study therefore suggests that while a wide range of ET rates are observed, there is an upper limit to rate enhancements achievable by only tuning the substituent electrophilicity.

First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel-Holstein Hamiltonian

A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].

TheoDORE: A toolbox for a detailed and automated analysis of electronic excited state computations

The advent of ever more powerful excited-state electronic structure methods has led to a tremendous increase in the predictive power of computation, but it has also rendered the analysis of these computations much more challenging and time-consuming. TheoDORE tackles this problem through providing tools for post-processing excited-state computations, which automate repetitive tasks and provide rigorous and reproducible descriptors. Interfaces are available for ten different quantum chemistry codes and a range of excited-state methods implemented therein. This article provides an overview of three popular functionalities within TheoDORE, a fragment-based analysis for assigning state character, the computation of exciton sizes for measuring charge transfer, and the natural transition orbitals used not only for visualization but also for quantifying multiconfigurational character. Using the examples of an organic push-pull chromophore and a transition metal complex, it is shown how these tools can be used for a rigorous and automated assignment of excited-state character. In the case of a conjugated polymer, we venture beyond the limits of the traditional molecular orbital picture to uncover spatial correlation effects using electron-hole correlation plots and conditional densities.

First-principles quantum simulations of exciton diffusion on a minimal oligothiophene chain at finite temperature

High-dimensional multiconfigurational quantum dynamics simulations are carried out at finite temperature to simulate exciton diffusion on an oligothiophene chain, representative of a segment of the poly(3-hexylthiophene) (P3HT) polymer. The ab initio parametrized site-based Hamiltonian of Binder et al. [Phys. Rev. Lett., 2018, 120, 227401] is employed to model a 20-site system, including intra-ring and inter-ring high-frequency modes as well as torsional modes which undergo thermal fluctuations induced by an explicit harmonic oscillator bath. The system-bath dynamics is treated within the setting of a stochastic mean-field Schrödinger equation. For the 20-site excitonic system, a total of 20 Frenkel states and 248 modes are propagated using the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method. The resulting dynamics can be interpreted in terms of the coherent motion of an exciton-polaron quasi-particle stochastically driven by torsional fluctuations. This dynamics yields a near-linear mean squared displacement (MSD) as a function of time, from which a diffusion coefficient can be deduced which increases with temperature, up to 5.7 × 10-3 cm2 s-1 at T = 300 K.

Synthesis and Density Functional Theory Studies of Azirinyl and Oxiranyl Functionalized Isoindigo and (3Z,3'Z)-3,3'-(ethane-1,2-diylidene)bis(indolin-2-one) Derivatives

The design and synthesis of functionalized isoindigo compounds by reaction of isoindigo with (S)-glycidyl tosylate, epibromohydrin, 2-(bromomethyl)-1-(arylsulfonyl)aziridine, and 2-(bromomethyl)-1-(alkylsulfonyl)aziridine in the presence of MeONa proceed under mild conditions in moderate yields. (3Z,3’Z)-3,3’-(Ethane-1,2-diylidene)bis(1-(oxiran-2-ylmethyl)indolin-2-one), with an extended central olefin π-conjugated moiety was also reacted with methyl-oxiranes to give the corresponding N,N’-disubstituted derivative. Calculations with DFT and TD-DFT of hypothetical isoindigo-thiophene DA molecules with various electron withdrawing substituents, including aziridine, oxirane, nitrile, carbonyl, and sulfonate, indicated that the proximity and strength of the functional group have a significant effect on the HOMO, LUMO, vertical excitation energy, and oscillator strength of the π–π* transitions.

Ground-State Geometry and Vibrations of Polyphenylenevinylene Oligomers

Conformational space of polyphenylenevinylene oligomers is systematically investigated computationally at energies relevant for room temperature dynamics in a solvent and in a solid state. Our calculations show that optimal oligomer structures are essentially planar. However, lack of a deep minimum at the planar geometry allows for large molecular deformations even at very low temperatures. At larger angles, rotational motion of dihedrals intermix with two orthogonal bending motions of the entire molecule. In a crystalline environment these degrees of freedom intermix with translational and rotational motions, whereas purely intramolecular modes are well separated. The reliability of our calculations is confirmed by an excellent match of the theoretical and experimental Raman spectra of crystalline stilbene in the entire spectral range including the low-frequency part. Obtained results provide important insights into nature of low-frequency vibrations, which play a key role in charge transport in organic semiconductors.

Ultrafast excitation energy transfer in a benzimidazole-naphthopyran donor-acceptor dyad

The excited-state dynamics of a donor-acceptor dyad composed of 1-propyl-2-pyridinyl-benzimidazole (PPBI) as donor and the photochromic molecular switch diphenylnaphthopyran (DPNP) as acceptor linked via an ester bridge has been investigated by a combination of static and time-resolved spectroscopies and quantum chemical calculations. The UV absorption spectrum of the dyad is virtually identical to the sum of the spectra of its individual constituents, indicating only weak electronic coupling between the donor and acceptor in the electronic ground state. After selective photoexcitation of the PPBI chromophore in the dyad at λpump = 310 nm, however, a fast electronic energy transfer (EET) from the donor to the acceptor is observed, by which the lifetime of the normally long-lived excited state of PPBI is reduced to a few ps. Enabled by the EET, the acceptor switches from its ring-closed naphtopyran form to its ring-opened merocyanine form. The singular value decomposition-based global analyses of the measured femtosecond time-resolved transient absorption spectra of the dyad and its two building blocks as reference compounds allowed us to determine a value for the EET time constant in the dyad of τ = 2.90 ± 0.60 ps. For comparison, Förster theory predicts characteristic FRET times between 1.2 ps ≤ τ ≤ 4.2 ps, in good agreement with the experimental result.

Conformational Dynamics Guides Coherent Exciton Migration in Conjugated Polymer Materials: First-Principles Quantum Dynamical Study

We report on high-dimensional quantum dynamical simulations of photoinduced exciton migration in a single-chain oligothiophene segment, in view of elucidating the controversial nature of the elementary exciton transport steps in semiconducting polymers. A novel first-principles parametrized Frenkel J aggregate Hamiltonian is employed that goes significantly beyond the standard Frenkel-Holstein Hamiltonian. Departing from a nonequilibrium state created by photoexcitation, these simulations provide evidence of an ultrafast two-timescale process at low temperatures, involving exciton-polaron formation within tens of femtoseconds (fs), followed by torsional relaxation on an ∼400  fs timescale. The second step is the driving force for exciton migration, as initial conjugation breaks are removed by dynamical planarization. The quantum coherent nature of the elementary exciton migration step is consistent with experimental observations highlighting the correlated and vibrationally coherent nature of the dynamics on ultrafast timescales.

The correspondence between the conformational and chromophoric properties of amorphous conjugated polymers in mesoscale condensed systems

For π-conjugated polymers, the notion of spectroscopic units or "chromophores" provides illuminating insights into the experimentally observed absorption/emission spectra and the mechanisms of energy/charge transfer. To date, however, no statistical analysis has revealed a direct correspondence between chromophoric and conformational properties-with the latter being fundamental to polymer semiconductors. Herein, we propose a "persistence length" calculation to re-evaluate chain conformation over a full conjugation length. The mesoscale condensed systems of MEH-PPV and MEH-PPV/C₆₀ hybrid (system size ∼10 × 10 × 10 nm³) are utilized as two prototypical model systems, along with a full range of segmental lengths (2-20-mer) and five lowest singlet excited states to hint at the generality of the features presented. We demonstrate, for the first time, that two properly re-defined conformational factors that characterize chain folding and planarity, respectively, capture excellently the population distribution of chromophores in both systems investigated. In contrast, the conventional strategy of utilizing two adjacent monomer units to characterize (local) chain conformation results in only an inconspicuous correlation between the two, as previously reported. It is further shown that chain folding-and not chain planarity-is more relevant in capturing the associated oscillator strength for the first excited state, where the transient dipole moments are known to align with the chain conformation, although the corresponding excitation energy and exciton size seem relatively unaffected. The observed effects of C₆₀ on the MEH-PPV adsorption spectra also agree with recent experimental trends. Overall, the present findings are expected to aid future multiscale computer simulations and spectroscopy-data interpretations for polymer semiconductors and their hybrid systems.

Charge-transfer states in triazole linked donor-acceptor materials: strong effects of chemical modification and solvation

A series of 1,2,3-triazole linked donor-acceptor chromophores are prepared by Click Chemistry from ene-yne starting materials. The effects of three distinct chemical variations are investigated: enhancing the acceptor strength through oxidation of the sulphur atom, alteration of the double bond configuration, and variation of the triazole substitution pattern. A detailed photophysical characterization shows that these alterations have a negligible effect on the absorption while dramatically altering the emission wavelengths. In addition, strong solvatochromism is found leading to significant red shifts in the case of polar solvents. The experimental findings are rationalized and related to the electronic structure properties of the chromophores by time-dependent density functional theory as well as the ab initio algebraic diagrammatic construction method for the polarization propagator in connection with a new formalism allowing to model the influence of solvation onto long-lived excited states and their emission energies. These calculations highlight the varying degree of intramolecular charge transfer character present for the different molecules and show that the amount of charge transfer is strongly modulated by the conducted chemical modifications, by the solvation of the chromophores, and by the structural relaxation in the excited state. It is, furthermore, shown that enhanced charge separation, as induced by chemical modification or solvation, reduces the singlet-triplet gaps and that two of the investigated molecules possess sufficiently low gaps to be considered as candidates for thermally activated delayed fluorescence.

Investigating the optical properties of BOIMPY dyes using ab initio tools

Using a computational approach combining Time-Dependent Density Functional Theory (TD-DFT) and second-order Coupled Cluster (CC2) approaches, we investigate the spectral properties of a large panel of BOIMPY dyes.

Excitons in poly(para phenylene vinylene): a quantum-chemical perspective based on high-level ab initio calculations

Exciton analyses of high-level quantum-chemical computations for poly(paraphenylene vinylene) reveal the nature of the excitonic bands in PPV oligomers.

Excited-state deactivation in 8-oxo-deoxyguanosine: comparison between anionic and neutral forms

Ab initio explorations of excited-state potential-energy surfaces show that a radiationless deactivation mechanism via intramolecular excited-state proton transfer is available in neutral 8-oxo-deoxyguanosine, whereas it is not available in the anionic form.

Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.

Statistical analysis of electronic excitation processes: Spatial location, compactness, charge transfer, and electron-hole correlation

We report the development of a set of excited-state analysis tools that are based on the construction of an effective exciton wavefunction and its statistical analysis in terms of spatial multipole moments. This construction does not only enable the quantification of the spatial location and compactness of the individual hole and electron densities but also correlation phenomena can be analyzed, which makes this procedure particularly useful when excitonic or charge-resonance effects are of interest. The methods are first applied to bianthryl with a focus on elucidating charge-resonance interactions. It is shown how these derive from anticorrelations between the electron and hole quasiparticles, and it is discussed how the resulting variations in state characters affect the excited-state absorption spectrum. As a second example, cytosine is chosen. It is illustrated how the various descriptors vary for valence, Rydberg, and core-excited states, and the possibility of using this information for an automatic characterization of state characters is discussed.

New tools for the systematic analysis and visualization of electronic excitations. II. Applications

The excited states of a diverse set of molecules are examined using a collection of newly implemented analysis methods. These examples expose the particular power of three of these tools: (i) natural difference orbitals (the eigenvectors of the difference density matrix) for the description of orbital relaxation effects, (ii) analysis of the one-electron transition density matrix in terms of an electron-hole picture to identify charge resonance and excitonic correlation effects, and (iii) state-averaged natural transition orbitals for a compact simultaneous representation of several states. Furthermore, the utility of a wide array of additional analysis methods is highlighted. Five molecules with diverse excited state characteristics are chosen for these tasks: pyridine as a prototypical small heteroaromatic molecule, a model system of six neon atoms to study charge resonance effects, hexatriene in its neutral and radical cation forms to exemplify the cases of double excitations and spin-polarization, respectively, and a model iridium complex as a representative metal organic compound. Using these examples a number of phenomena, which are at first sight unexpected, are highlighted and their physical significance is discussed. Moreover, the generality of the conclusions of this paper is verified by a comparison of single- and multireference ab initio methods.

New tools for the systematic analysis and visualization of electronic excitations. I. Formalism

A variety of density matrix based methods for the analysis and visualization of electronic excitations are discussed and their implementation within the framework of the algebraic diagrammatic construction of the polarization propagator is reported. Their mathematical expressions are given and an extensive phenomenological discussion is provided to aid the interpretation of the results. Starting from several standard procedures, e.g., population analysis, natural orbital decomposition, and density plotting, we proceed to more advanced concepts of natural transition orbitals and attachment/detachment densities. In addition, special focus is laid on information coded in the transition density matrix and its phenomenological analysis in terms of an electron-hole picture. Taking advantage of both the orbital and real space representations of the density matrices, the physical information in these analysis methods is outlined, and similarities and differences between the approaches are highlighted. Moreover, new analysis tools for excited states are introduced including state averaged natural transition orbitals, which give a compact description of a number of states simultaneously, and natural difference orbitals (defined as the eigenvectors of the difference density matrix), which reveal details about orbital relaxation effects.

Absorption and fluorescence spectra of poly(p-phenylenevinylene) (PPV) oligomers: an ab initio simulation

The absorption and fluorescence spectra of poly(p-phenylenevinylene) (PPV) oligomers with up to seven repeat units were theoretically investigated using the algebraic diagrammatic construction method to second order, ADC(2), combined with the resolution-of-the-identity (RI) approach. The ground and first excited state geometries of the oligomers were fully optimized. Vertical excitation energies and oscillator strengths of the first four transitions were computed. The vibrational broadening of the absorption and fluorescence spectra was studied using a semiclassical nuclear ensemble method. After correcting for basis set and solvent effects, we achieved a balanced description of the absorption and fluorescence spectra by means of the ADC(2) approach. This fact is documented by the computed Stokes shift along the PPV series, which is in good agreement with the experimental values. The experimentally observed band width of the UV absorption and fluorescence spectra is well reproduced by the present simulations showing that the nuclear ensemble generated should be well suitable for consecutive surface hopping dynamics simulations.

Electronic excitation and structural relaxation of the adenine dinucleotide in gas phase and solution

The excited states and potential surfaces of the adenine dinucleotide are analyzed in gas phase and in solution using a correlated ab initio methodology in a QM/MM framework. In agreement with previous studies, a rather flat S1 surface with a number of minima of different character is found. Specifically, our results suggest that exciplexes with remarkably short intermolecular separation down to ~2.0 Å are formed. A detailed analysis shows that due to strong orbital interactions their character differs significantly from any states present in the Franck-Condon region. The lowest S1 energy minimum is a ππ* exciplex with only a small amount of charge transfer. It possesses appreciable oscillator strength with a polarization almost perpendicular to the planes of the two adenine molecules.